Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate f...Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.展开更多
In this paper,we define arbitrarily high-order energy-conserving methods for Hamilto-nian systems with quadratic holonomic constraints.The derivation of the methods is made within the so-called line integral framework...In this paper,we define arbitrarily high-order energy-conserving methods for Hamilto-nian systems with quadratic holonomic constraints.The derivation of the methods is made within the so-called line integral framework.Numerical tests to illustrate the theoretical findings are presented.展开更多
In this paper, we consider the problem of optimization of adaptive direct methods of operator equations. Adaptivity of a direct method is understood in the sense that the subspace on the basis of which it is construct...In this paper, we consider the problem of optimization of adaptive direct methods of operator equations. Adaptivity of a direct method is understood in the sense that the subspace on the basis of which it is constructed is chosen depending on the operator of the concrete equation (otherwise, nonadaptive direct method is then concerned), which would essentially let us increase the precision. For some classes of the second kind of Fredhlom integral equations with anisotropic smooth kernels we determine the exact order of the error of adaptive direct methods, and we also give an optimal algorithm.展开更多
This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic So...This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.展开更多
Consider the scattering of a time-harmonic electromagnetic plane wave by an arbitrarily shaped and filled cavity embedded in a perfect electrically conducting infinite ground plane.A method of symmetric coupling of fi...Consider the scattering of a time-harmonic electromagnetic plane wave by an arbitrarily shaped and filled cavity embedded in a perfect electrically conducting infinite ground plane.A method of symmetric coupling of finite element and boundary integral equations is presented for the solutions of electromagnetic scattering in both transverse electric and magnetic polarization cases.Given the incident field,the direct problem is to determine the field distribution from the known shape of the cavity;while the inverse problem is to determine the shape of the cavity from the measurement of the field on an artificial boundary enclosing the cavity.In this paper,both the direct and inverse scattering problems are discussed based on a symmetric coupling method.Variational formulations for the direct scattering problem are presented,existence and uniqueness of weak solutions are studied,and the domain derivatives of the field with respect to the cavity shape are derived.Uniqueness and local stability results are established in terms of the inverse problem.展开更多
A detailed fracture mechanics analysis of bridge-toughening in a fiber reinforced composite is presented in this paper. The integral equation governing bridge-toughening as well as crack opening displacement (COD) for...A detailed fracture mechanics analysis of bridge-toughening in a fiber reinforced composite is presented in this paper. The integral equation governing bridge-toughening as well as crack opening displacement (COD) for the composite with interfacial layer is derived from the Castigliano's theorem and interface shear-lag model. A numerical result of the COD equation is obtained using the iteration solution of the second Fredholm integral equation. In order to investigate the effect of various parameters on the toughening, an approximate analytical solution of the equation is present and its error analysis is performed, which demonstrates the approximate solution to be appropriate. A parametric study of the influence of the crack length, interfacial shear modules, thickness of the interphase, fiber radius, fiber volume fraction and properties of materials on composite toughening is therefore carried out. The results are useful for experimental demonstration and toughening design including the fabrication process of the composite.展开更多
In the current work, transient heat conduction in a semi-infinite medium is considered for its many applications in various heat fields. Here, the homotopy analysis method （HAM） is applied to solve this problem and ...In the current work, transient heat conduction in a semi-infinite medium is considered for its many applications in various heat fields. Here, the homotopy analysis method （HAM） is applied to solve this problem and analytical results are compared with those of the exact and integral methods results. The results show that the HAM can give much better approximations than the other approximate methods： Changes in heat fluxes and profiles of temperature are obtained at different times and positions for copper, iron and aluminum.展开更多
This work is a continuation of the earlier article [1]. We establish new numerical methods for solving systems of Volterra integral equations with cardinal splines. The unknown functions are expressed as a linear comb...This work is a continuation of the earlier article [1]. We establish new numerical methods for solving systems of Volterra integral equations with cardinal splines. The unknown functions are expressed as a linear combination of horizontal translations of certain cardinal spline functions with small compact supports. Then a simple system of equations on the coefficients is acquired for the system of integral equations. It is relatively straight forward to solve the system of unknowns and an approximation of the original solution with high accuracy is achieved. Several cardinal splines are applied in the paper to enhance the accuracy. The sufficient condition for the existence of the inverse matrix is examined and the convergence rate is investigated. We demonstrated the value of the methods using several examples.展开更多
From the dislocation type solution of the torsion of single crack, by using the concept of finite part integrals. The torsion problem of cylinder with a single crack was reduced into an integral equation with strong s...From the dislocation type solution of the torsion of single crack, by using the concept of finite part integrals. The torsion problem of cylinder with a single crack was reduced into an integral equation with strong singularity. The numerical method was also obtained and several numerical examples were calculated successfully at the end of this paper.,展开更多
In this paper, the exact Bayesian limits, taking conjugate and noninformative prior distribution, and the exact fiducial limits for the mean of the lognormal distribution are presented. They can be found iteratively b...In this paper, the exact Bayesian limits, taking conjugate and noninformative prior distribution, and the exact fiducial limits for the mean of the lognormal distribution are presented. They can be found iteratively by one-dimension integral on a finite interval. The new algorithm is very convenient and with high accuracy. It can meet the practical engineering need excellently. However, the primitive algorithm is rather cumbersome. By the way, the very close approximate limits with a simple algorithm are derived. They can be applied immediately to engineering. Otherwise, they can also be used as a search interval to find the root of equation for the exact limits.展开更多
Abstract The reconstruction of cylindrically layered media is investigated in this article. The inverse problem is modeled using a source-type integral equation with a series of cylindrical waves as incidences, and a ...Abstract The reconstruction of cylindrically layered media is investigated in this article. The inverse problem is modeled using a source-type integral equation with a series of cylindrical waves as incidences, and a conventional Born iterative procedure is modified for solving the integral equation. In the modified iterative procedure, a conventional single-point approximation for the calculation of the field inside media is replaced by a multi-points approximation to improve the numerical stability of its solution. Numerical simulations for different permittivity distributions are demonstrated in terms of artificial scattering data with the procedure. The result shows that the procedure enjoys both accuracy and stability in the numerical computation.展开更多
基金Aeronautical Science Foundation of China (99A52007)
文摘Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.
文摘In this paper,we define arbitrarily high-order energy-conserving methods for Hamilto-nian systems with quadratic holonomic constraints.The derivation of the methods is made within the so-called line integral framework.Numerical tests to illustrate the theoretical findings are presented.
基金This work is supported by the Special Funds for Major State Basic Research Projects (Grant No. G19990328)the Zhejiang Provincial Natural Science Foundation (Grant No. 100002).
文摘In this paper, we consider the problem of optimization of adaptive direct methods of operator equations. Adaptivity of a direct method is understood in the sense that the subspace on the basis of which it is constructed is chosen depending on the operator of the concrete equation (otherwise, nonadaptive direct method is then concerned), which would essentially let us increase the precision. For some classes of the second kind of Fredhlom integral equations with anisotropic smooth kernels we determine the exact order of the error of adaptive direct methods, and we also give an optimal algorithm.
基金Project supported by the Natural Science Foundation of China(10371009)Research Fund for the Doctoral Program Higher Education
文摘This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.
基金the NSF grants DMS-0908325,CCF-0830161,EAR-0724527,DMS-0968360the ONR grant N00014-09-1-0384 and a special research grant from Zhejiang University.The research of PL was supported in part by the NSF grants EAR-0724656,DMS-0914595,and DMS-1042958.
文摘Consider the scattering of a time-harmonic electromagnetic plane wave by an arbitrarily shaped and filled cavity embedded in a perfect electrically conducting infinite ground plane.A method of symmetric coupling of finite element and boundary integral equations is presented for the solutions of electromagnetic scattering in both transverse electric and magnetic polarization cases.Given the incident field,the direct problem is to determine the field distribution from the known shape of the cavity;while the inverse problem is to determine the shape of the cavity from the measurement of the field on an artificial boundary enclosing the cavity.In this paper,both the direct and inverse scattering problems are discussed based on a symmetric coupling method.Variational formulations for the direct scattering problem are presented,existence and uniqueness of weak solutions are studied,and the domain derivatives of the field with respect to the cavity shape are derived.Uniqueness and local stability results are established in terms of the inverse problem.
基金National Natural Science Foundatjon and China Postdoctoral Scjence Fbundation
文摘A detailed fracture mechanics analysis of bridge-toughening in a fiber reinforced composite is presented in this paper. The integral equation governing bridge-toughening as well as crack opening displacement (COD) for the composite with interfacial layer is derived from the Castigliano's theorem and interface shear-lag model. A numerical result of the COD equation is obtained using the iteration solution of the second Fredholm integral equation. In order to investigate the effect of various parameters on the toughening, an approximate analytical solution of the equation is present and its error analysis is performed, which demonstrates the approximate solution to be appropriate. A parametric study of the influence of the crack length, interfacial shear modules, thickness of the interphase, fiber radius, fiber volume fraction and properties of materials on composite toughening is therefore carried out. The results are useful for experimental demonstration and toughening design including the fabrication process of the composite.
文摘In the current work, transient heat conduction in a semi-infinite medium is considered for its many applications in various heat fields. Here, the homotopy analysis method （HAM） is applied to solve this problem and analytical results are compared with those of the exact and integral methods results. The results show that the HAM can give much better approximations than the other approximate methods： Changes in heat fluxes and profiles of temperature are obtained at different times and positions for copper, iron and aluminum.
文摘This work is a continuation of the earlier article [1]. We establish new numerical methods for solving systems of Volterra integral equations with cardinal splines. The unknown functions are expressed as a linear combination of horizontal translations of certain cardinal spline functions with small compact supports. Then a simple system of equations on the coefficients is acquired for the system of integral equations. It is relatively straight forward to solve the system of unknowns and an approximation of the original solution with high accuracy is achieved. Several cardinal splines are applied in the paper to enhance the accuracy. The sufficient condition for the existence of the inverse matrix is examined and the convergence rate is investigated. We demonstrated the value of the methods using several examples.
基金Project supported by P.H.D.Foundation of the State Education Commission of China
文摘From the dislocation type solution of the torsion of single crack, by using the concept of finite part integrals. The torsion problem of cylinder with a single crack was reduced into an integral equation with strong singularity. The numerical method was also obtained and several numerical examples were calculated successfully at the end of this paper.,
文摘In this paper, the exact Bayesian limits, taking conjugate and noninformative prior distribution, and the exact fiducial limits for the mean of the lognormal distribution are presented. They can be found iteratively by one-dimension integral on a finite interval. The new algorithm is very convenient and with high accuracy. It can meet the practical engineering need excellently. However, the primitive algorithm is rather cumbersome. By the way, the very close approximate limits with a simple algorithm are derived. They can be applied immediately to engineering. Otherwise, they can also be used as a search interval to find the root of equation for the exact limits.
基金the National Natural Science Foundation of China（60671065）.
文摘Abstract The reconstruction of cylindrically layered media is investigated in this article. The inverse problem is modeled using a source-type integral equation with a series of cylindrical waves as incidences, and a conventional Born iterative procedure is modified for solving the integral equation. In the modified iterative procedure, a conventional single-point approximation for the calculation of the field inside media is replaced by a multi-points approximation to improve the numerical stability of its solution. Numerical simulations for different permittivity distributions are demonstrated in terms of artificial scattering data with the procedure. The result shows that the procedure enjoys both accuracy and stability in the numerical computation.
