A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.

RESEARCH ON THE COMPANION SOLUTION FOR A THIN PLATE IN THE MESHLESS LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and element free Galerkin method(EFGM), and is a truly meshless method possessing wide prospects in engineering applications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method （LBIEM）, singularities in the local boundary integrals need to be treated specially. In the current paper, local integral equations are adopted for the nodes inside the domain and moving least square approximation （MLSA） for the nodes on the global boundary, thus singularities will not occur in the new al- gorithm. At the same time, approximation errors of boundary integrals are reduced significantly. As applications and numerical tests, Laplace equation and Helmholtz equation problems are considered and excellent numerical results are obtained. Furthermore, when solving the Helmholtz problems, the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions. Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.

EQUIVALENT BOUNDARY INTEGRAL EQUATIONS WITH INDIRECT UNKNOWNS FOR THIN ELASTIC PLATE BENDING THEORY

Equivalent Boundary Integral Equations (EBIE) with indirect unknowns for thin elastic plate bending theory, which is equivalent to the original boundary value problem, is established rigorously by mathematical technique of non-analytic continuation and is fully proved by means of the variational principle. The previous three kinds of boundary integral equations with indirect unknowns are discussed thoroughly and it is shown that all previous results are not EBIE.

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.

Analytical Solution of Boundary Integral Equations for 2-D Steady Linear Wave Problems

Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed,free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigensolution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore,an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.

BOUNDARY INTEGRAL EQUATIONS FOR THE BENDING PROBLEM OF PLATES ON TWO-PARAMETER FOUNDATION

By means of Fourier integral transformation of generalized function, the fundamental solution for the bending problem of plates on two-parameter foundation is derived in this paper, and the fundamental solution is expanded into a uniformly convergent series. On the basis of the above work, two boundary integral equations which are suitable to arbitrary shapes and arbitrary boundary conditions are established by means of the Rayleigh-Green identity. The content of the paper provides the powerful theories for the application of BEM in this problem.

Roll Flattening Analytical Model in Flat Rolling by Boundary Integral Equation Method

In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening calculation based on semi infinite body model, especially near the two roll barrel edges, a new and more accurate roll flattening model is proposed. Based on boundary integral equation method, an analytical model for solving a finite length semi infinite body is established. The lateral surface displacement field of the finite length semi-infinite body is simulated by finite element method （FEM） and lateral surface displacement decay functions are established. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distribu ted force is obtained and an accurate roll flattening model is established. Different from the traditional semi-infinite body model, the matrix form of the new roll flattening model is established through the mathematical derivation. The result from the new model is more consistent with that by FEM especially near the edges.

SPLITTING EXTRAPOLATIONS FOR SOLVING BOUNDARY INTEGRAL EQUATIONS OF LINEAR ELASTICITY DIRICHLET PROBLEMS ON POLYGONS BY MECHANICAL QUADRATURE METHODS

Taking hm as the mesh width of a curved edge Гm （m = 1, ..., d ） of polygons and using quadrature rules for weakly singular integrals, this paper presents mechanical quadrature methods for solving BIES of the first kind of plane elasticity Dirichlet problems on curved polygons, which possess high accuracy O（h0^3） and low computing complexities. Since multivariate asymptotic expansions of approximate errors with power hi^3 （i = 1, 2, ..., d） are shown, by means of the splitting extrapolations high precision approximations and a posteriori estimate are obtained.

ERROR ANALYSIS FOR A FAST NUMERICAL METHOD TO A BOUNDARY INTEGRAL EQUATION OF THE FIRST KIND

For two-dimensional boundary integral equations of the first kind with logarithmic kernels, the use of the conventional boundary element methods gives linear systems with dense matrix. In a recent work [J. Comput. Math., 22 （2004）, pp. 287-298], it is demonstrated that the dense matrix can be replaced by a sparse one if appropriate graded meshes are used in the quadrature rules. The numerical experiments also indicate that the proposed numerical methods require less computational time than the conventional ones while the formal rate of convergence can be preserved. The purpose of this work is to establish a stability and convergence theory for this fast numerical method. The stability analysis depends on a decomposition of the coefficient matrix for the collocation equation. The formal orders of convergence observed in the numerical experiments are proved rigorously.

Fast Multipole Accelerated Boundary Integral Equation Method for Evaluating the Stress Field Associated with Dislocations in a Finite Medium

In this paper,we develop an efficient numericalmethod based on the boundary integral equation formulation and new version of fast multipole method to solve the boundary value problem for the stress field associated with dislocations in a finite medium.Numerical examples are presented to examine the influence from material boundaries on dislocations.

The Collocation Method and the Splitting Extrapolation for the First Kind of Boundary Integral Equations on Polygonal Regions

In this paper,the collocation methods are used to solve the boundary integral equations of the first kind on the polygon.By means of Sidi’s periodic transformation and domain decomposition,the errors are proved to possess the multi-parameter asymptotic expansion at the interior point with the powers h^(3)/_(i)(i=1,...,d),which means that the approximations of higher accuracy and a posteriori estimation of the errors can be obtained by splitting extrapolations.Numerical experiments are carried out to show that the methods are very efficient.

ASYMPTOTIC ERROR EXPANSION FOR THE NYSTROM METHOD OF NONLINEAR VOLTERRA INTEGRAL EQUATION OF THE SECOND KIND

While the numerical solution of one-dimensional Volterra integral equations of the second kind with regular kernels is well understood, there exist no systematic studies of asymptotic error expansion for the approximate solution. In this paper,we analyse the Nystrom solution of one-dimensional nonlinear Volterra integral equation of the second kind and show that approkimate solution admits an asymptotic error expansion in even powers of the step-size h, beginning with a term in h2. So that the Richardson's extrapolation can be done. This will increase the accuracy of numerical solution greatly.

SINGULAR INTEGRAL EQUATIONS AND BOUNDARY ELEMENT METHOD OF CRACKS IN THERMALLY STRESSED PLANAR SOLIDS

Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, and then,its numerical method is proposed by the numerical method of the singular integral equations combined with boundary element method. Moreover, the singular nature of temperature gradient field near the crack front is proved by the main-part analysis method of the singular integral equation, and the singular temperature gradients are exactly obtained. Finally, several typical examples calculated.

Iterative convergence of boundary-volume integral equation method

The boundary-volume integral equation numerical technique can be a powerful tool for piecewise heterogeneous media, but it is limited to small problems or low frequencies because of great computational cost. Therefore, a restarted GMRES method is applied to solve large-scale boundary-volume scattering problems in this paper to overcome the computational barrier. The iterative method is firstly applied to responses of dimensionless frequency to a semicircular alluvial valley filled with sediments, compared with the standard Gaussian elimination method. Then the method is tested by a heterogeneous multilayered model to show its applicability. Numerical experiments indicate that the preconditioned GMRES method can significantly improve computational efficiency especially for large Earth models and high frequencies, but with a faster convergence for the left diagonal preconditioning.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

A moving Kriging interpolation-based boundary node method for two-dimensional potential problems

In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.

AN IMPROVED HYBRID BOUNDARY NODE METHOD IN TWO-DIMENSIONAL SOLIDS

The hybrid boundary node method （HBNM） is a promising method for solving boundary value problems with the hybrid displacement variational formulation and shape functions from the moving least squares（MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the latter. Following its application in solving potential problems, it is further developed and numerically implemented for 2D solids in this paper. The rigid movement method is employed to solve the hyper-singular integrations. Numerical examples for some 2D solids have been given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method are studied through numerical examples.

WAVE SUPERPOSITION METHOD BASED ON VIRTUAL SOURCE BOUNDARY WITH COMPLEX RADIUS VECTOR FOR SOLVING ACOUSTIC RADIATION PROBLEMS

By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at the eigenfrequencies of the corresponding interior problems is proposed. By adding appropriate damp to the virtual source system of the wave superposition method, the unique solutions for all wave numbers can be ensured. Based on this thought, a novel method-wave superposition method with complex radius vector is constructed. Not only is the computational time of this method approximately equal to that of the standard wave superposition method, but also the accuracy is much higher compared with other correlative methods. Finally, by taking the pulsating sphere and oscillating sphere as examples, the results of calculation show that the present method can effectively overcome the non-uniqueness problem.

Interaction of Linear Waves with Infinitely Long Horizontal Cylinders Studied by Boundary Element Method

The two-dimensional problems concerning the interaction of linear water waves with cylinders of arbitrary shape in two-layer deep water are investigated by use of the Boundary Integral Equation Method （BIEM）. Simpler new expressions for the Green functions are derived, and verified by comparison of results obtained by BIEM with those by an analytical method. Examined are the radiation and scattering of linear waves by two typical configurations of cylinders in two- layer deep water. Hydrodynamic behaviors including hydrodynamic coefficients, wave forces, reflection and transmission coef- ficients and energies are analyzed in detail, and some imeresting physical phenomena are observed.