A Fourth-Order Kernel-Free Boundary Integral Method for Interface Problems

This paper presents a fourth-order Cartesian grid based boundary integral method(BIM)for heterogeneous interface problems in two and three dimensional space,where the problem interfaces are irregular and can be explicitly given by parametric curves or implicitly defined by level set functions.The method reformulates the governing equation with interface conditions into boundary integral equations(BIEs)and reinterprets the involved integrals as solutions to some simple interface problems in an extended regular region.Solution of the simple equivalent interface problems for integral evaluation relies on a fourth-order finite difference method with an FFT-based fast elliptic solver.The structure of the coefficient matrix is preserved even with the existence of the interface.In the whole calculation process,analytical expressions of Green’s functions are never determined,formulated or computed.This is the novelty of the proposed kernel-free boundary integral(KFBI)method.Numerical experiments in both two and three dimensions are shown to demonstrate the algorithm efficiency and solution accuracy even for problems with a large diffusion coefficient ratio.

Accuracy Analysis of the Boundary Integral Method for Steady Navier-Stokes Equations around a Rotating Obstacle

This paper deals with the boundary integral method to study the Navier-Stokes equations around a rotating obstacle. The detail of this method is that the exterior domain is truncated into a bounded domain and a new exterior domain by introducing some open ball BR, and the nonlinear problem in the bounded domain and the linearized problem in the new exterior domain are considered and the approximation coupled problem is obtained. We show that the error between the solution u of Navier-Stokes equations around a rotating obstacle and the solution ue of the approximation coupled problem is O（R-1/4） in the Hl-seminorm when Iwl does not exceed some constant.

A Spectrally Accurate Boundary Integral Method for Interfacial Velocities in Two-Dimensional Stokes Flow

We present a new numerical method for solving two-dimensional Stokes flow with deformable interfaces such as dynamics of suspended drops or bubbles.The method is based on a boundary integral formulation for the interfacial velocity and is spectrally accurate in space.We analyze the singular behavior of the integrals(single-layer and double-layer integrals)appearing in the equations.The interfaces are formulated in the tangent angle and arc-length coordinates and,to reduce the stiffness of the evolution equation,the marker points are evenly distributed in arc-length by choosing a proper tangential velocity along the interfaces.Examples of Stokes flow with bubbles are provided to demonstrate the accuracy and effectiveness of the numerical method.

A Kernel-Free Boundary Integral Method for Variable Coefficients Elliptic PDEs

This work proposes a generalized boundary integral method for variable coefficients elliptic partial differential equations(PDEs),including both boundary value and interface problems.The method is kernel-free in the sense that there is no need to know analytical expressions for kernels of the boundary and volume integrals in the solution of boundary integral equations.Evaluation of a boundary or volume integral is replaced with interpolation of a Cartesian grid based solution,which satisfies an equivalent discrete interface problem,while the interface problem is solved by a fast solver in the Cartesian grid.The computational work involved with the generalized boundary integral method is essentially linearly proportional to the number of grid nodes in the domain.This paper gives implementation details for a secondorder version of the kernel-free boundary integral method in two space dimensions and presents numerical experiments to demonstrate the efficiency and accuracy of the method for both boundary value and interface problems.The interface problems demonstrated include those with piecewise constant and large-ratio coefficients and the heterogeneous interface problem,where the elliptic PDEs on two sides of the interface are of different types.

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.

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.

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.

Turbulent flow in converging nozzles,part one:boundary layer solution

The boundary layer integral method is used to investigate the development of the turbulent swirling flow at the entrance region of a conical nozzle. The governing equations in the spherical coordinate system are simplified with the boundary layer as- sumptions and integrated through the boundary layer. The resulting sets of differential equations are then solved by the fourth-order Adams predictor-corrector method. The free vortex and uniform velocity profiles are applied for the tangential and axial velocities at the inlet region, respectively. Due to the lack of experimental data for swirling flows in converging nozzles, the developed model is validated against the numerical simulations. The results of numerical simulations demonstrate the capability of the analytical model in predicting boundary layer parameters such as the boundary layer growth, the shear rate, the boundary layer thickness, and the swirl intensity decay rate for different cone angles. The proposed method introduces a simple and robust procedure to investigate the boundary layer parameters inside the converging geometries.

A MESHLESS METHOD OF A THIN PLATE

A meshless approach to analysis of arbitrary Kirechhoff plates bythe local boundary integral equation (LBIF) method is presented. Themethod combines the advantageous features of all the three meth- ods:the Galerkin finite element method (GFEM), the boundary elementmethod (BEM) and the element- free Galerkin method (EFGM). It is atruly meshless method, which means that the discretization is inde-pendent of geometric subdivision into elements or cells, but is onlyboundary integration, however over a local boundary cen- tered) overa domain in question.

A Non-Singular Boundary Element Method for Interactions between Acoustical Field Sources and Structures

Localized point sources(monopoles)in an acoustical domain are implemented to a three dimensional non-singular Helmholtz boundary element method in the frequency domain.It allows for the straightforward use of higher order surface elements on the boundaries of the problem.It will been shown that the effect of the monopole sources ends up on the right hand side of the resulting matrix system.Some carefully selected examples are studied,such as point sources near and within a concentric spherical core-shell scatterer(with theoretical verification),near a curved focusing surface and near a multi-scale and multi-domain acoustic lens.

A Fast Accurate Boundary IntegralMethod for Potentials on Closely Packed Cells

Boundary integral methods are naturally suited for the computation of harmonic functions on a region having inclusions or cells with different material properties.However,accuracy deteriorates when the cell boundaries are close to each other.We present a boundary integralmethod in two dimensions which is specially designed tomaintain second order accuracy even if boundaries are arbitrarily close.Themethod uses a regularization of the integral kernel which admits analytically determined corrections to maintain accuracy.For boundaries with many components we use the fast multipolemethod for efficient summation.We compute electric potentials on a domain with cells whose conductivity differs from that of the surrounding medium.We first solve an integral equation for a source term on the cell interfaces and then find values of the potential near the interfaces via integrals.Finally we use a Poisson solver to extend the potential to a regular grid covering the entire region.A number of examples are presented.We demonstrate that increased refinement is not needed to maintain accuracy as interfaces become very close.

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.

Numerical quadrature for singular and near-singular integrals of boundary element method and its applications in large-scale acoustic problems

The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed, by which the conventional and fast multipole BEMs （boundary element methods） for 3D acoustic problems based on constant elements are improved. To solve the problem of singular integrals, a Hadamard finite-part integral method is presented, which is a simplified combination of the methods proposed by Kirkup and Wolf. The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART （Projection and Angular ＆ Radial Transformation）. The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab. In addition, the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution. The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations. A large-scale acoustic scattering problem, whose degree of freedoms is about 340,000, is implemented successfully. The results show that, the near singularity is primarily introduced by the hyper-singular kernel, and has great influences on the precision of the solution. The precision of fast multipole BEM is the same as conventional BEM, but the computational complexities are much lower.

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.

3-D NUMERICAL SIMULATION OF UNDERWATER EXPLOSION BUBBLE

The problem is solved by use of the boundary integral equation method and the high order finite element discretization on the assumption that the flow of water is inviscid, incompressible and irrotational, and the bubble gas obeys the isoentropic rule. The evolution of bubble is simulated by means of mixed Euler-Lagrangian method. The comparison of the numerical results with the 'exact solution' of a spheric bub ble and experimental data shows that the mathematical model is reasonable and the 3-D numerical method has high accuracy. Important behaviour of bubble, such as immigration and jetting, and dynamic load on an ellipsoid near explosion are derived.

Effect of a rigid structure on the dynamics of a bubble beneath the free surface

The coupling effect between a pulsating bubble and a free surface near a rigid structure is a complicated physical process.In this study,the evolution of an underwater explosion bubble and the free surface near a rigid structure is modeled by the boundary integral method.An approach of“double-node method”is used to maintain the stability of fluid-structure junction in the simulations,and meshes on the free surface and the structure are transformed to an open domain to ensure the calculation accuracy and efficiency.Validations are conducted against an underwater explosion experiment near a rigid structure.As a result,the simulations trace the jetting behavior of the bubble and the rise of the free surface.Finally,the bubble migration and the height of the free surface for different structure draughts are analyzed.

Bidirectional reflectance distribution function modeling of one-dimensional rough surface in the microwave band

An efficiently iterative analytical-numerical method is proposed for two-dimensional （2D） electromagnetic scattering from a perfectly electric conducting （PEC） target buried under a dielectric rough surface. The basic idea is to employ the Kirchhoff approximation （KA） to accelerate the boundary integral method （BIM）. Below the rough surface, an iterative system is designed between the rough surface and the target. The KA is used to simulate the initial field on the rough surface based on the Fresnel theory, while the target is analyzed by the boundary integral method to obtain a precise result. The fields between the rough surface and the target can be linked by the boundary integral equations below the rough surface. The technique presented here is highly efficient in terms of computational memory, time, and versatility. Numerical simulations of two typical models are carried out to validate the method.

MOTION AND DEFORMATION OF VISCOUS DROP IN STOKES FLOW NEAR RIGID WALL

A boundary integral method was developed for simulating the motion and deformation of a viscous drop in an axisymmetric ambient Stokes flow near a rigid wall and for direct calculating the stress on the wall. Numerical experiments by the method were performed for different initial stand-off distances of the drop to the wall, viscosity ratios, combined surface tension and buoyancy parameters and ambient flow parameters. Numerical results show that due to the action of ambient flow and buoyancy the drop is compressed and stretched respectively in axial and radial directions when time goes. When the ambient flow action is weaker than that of the buoyancy the drop raises and bends upward and the stress on the wall induced by drop motion decreases when time advances. When the ambient flow action is stronger than that of the buoyancy the drop descends and becomes flatter and flatter as time goes. In this case when the initial stand-off distance is large the stress on the wall increases as the drop evolutes but when the stand-off distance is small the stress on the wall decreases as a result of combined effects of ambient flow, buoyancy and the stronger wall action to the flow. The action of the stress on the wall induced by drop motion is restricted in an area near the symmetric axis, which increases when the initial stand-off distance increases. When the initial stand-off distance increases the stress induced by drop motion decreases substantially. The surface tension effects resist the deformation and smooth the profile of the drop surfaces. The drop viscosity will reduce the deformation and migration of the drop.