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 trod 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 equa- tion problems are considered and excellent numerical results are obtained.Furthermore, when solving the Hehnholtz 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.     

调和函数边界积分方程的充要条件

本文以调和函数的边值问题为例,探讨了边界积分方程的充要条件.文中首次提出了超定问题的概念,并建立了超定问题有解的一个充要条件,它也就是直接变量边界积分方程的一个充要条件.文中首次阐明了边界积分方程与变分原理的内在的联系,还指出了间接变量与直接变量两类边界积分方程之间存在着一一对应的关系.文中的慨念、思路和论点不难用于其它有变分原理的问题的边界积分方程.     

Steady Free-surface Flow Past an Uneven Channel Bottom

New free-surface flows past a semi-infinite ‘step’ in the bottom of a channel are considered. Surface tension is neglected but gravity is included in the dynamic boundary condition. Fully nonlinear solutions are computed by boundary integral equation methods. Additional weakly nonlinear solutions are derived analytically. A thorough analysis of the weakly nonlinear problem provides a systematic approach to identify all the possible types of solutions and the number of independent parameters.     

A complete boundary integral formulation for steady compressible inviscid flows governed by non-linear equations

A complete boundary integral formulation for steady compressible inviscid flows governed by non-linear equations is established by using the specific mass flux as a dependent variable. Thus, the dimensionality of the problem to be solved is reduced by one and the computational mesh to be generated is needed only on the boundary of the domain. It is shown that the boundary integral formulation developed in this paper is equivalent of the results of distributions of the fundamental solutions of the Laplacian operator equation with a different order along the boundaries of the domain. Hence, we have succeeded in establishing the fundamental-solution method for compressible inviscid flows governed by non-linear equations.     

Boundary integral equations of unique solutions in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＢＥＭ）ｐｒｏｖｉｄｅｓａｎａｔｔｒａｃｔｉｖｅａｌｔｅｒｎａｔｉｖｅｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｉｔｓｍａｉｎａｄｖａｎｔａｇｅｓａｒｅｅｃｏｎｏｍｉｃａｌａｎｄｐａｒｔｉｃｕｌａｒｌｙｃｏｎｖｅｎｉｅｎｔｆｏｒｕｎｂｏｕｎｄｅｄｄｏｍａｉｎａｎｄｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｐｒｏｂｌｅｍｓ．Ｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎ（ＢＩＥ）ｉｓｔｈｅ…     

Boundary integral/spectral element approaches to the Navier-Stokes equations

Numerical algorithms are presented which combine spectral expansions on elemental subdomains with boundary integral formulations for solving viscous flow problems. Three distinct algorithms are described. The first demonstrates the use of spectral elements for the classic boundary integral method for steady Stokes flow. The second extends this algorithm to include domain integrals for solution of the unsteady Navier-Stokes equations. The third algorithm explores the use of boundary integrals as a means of consolidating uncoupled elemental solutions in a domain decomposition approach. Numerical results demonstrating high-order convergence are presented in each case.     

Wave kinematic fields from the boundary integral method

The predictive potential of interior domain solutions from the boundary integral method for 2D extreme wave kinematics is explored. Comparisons with analytical solutions for near‐limit waves confirms the susceptibility of the boundary integral method to poor precision at near‐boundary locations. Additionally, these comparisons identify a domain‐wide precision challenge that is associated with the relatively rapid changes in water surface geometry and kinematics that are typical of extreme waves. A numerical evaluation of Green's integral around the boundary addresses these precision issues through formulation of the integration as a simultaneous system of ordinary differential equations at a cubic level of approximation. Careful attention is given to consistent interpolation of all contributions to the Green's integral. Copyright © 2005 John Wiley & Sons, Ltd.     

Boundary integral equations in quasisteady problems of capillary fluid mechanics Part 2: Application of the stress-stream function

In this paper planar viscous flows with a free boundary are further studied using the quasisteady approximation [1]. The introduction of the bianalytical stress-stream function provides an opportunity to adopt the theory of analytical functions. The mode of construction of the Fredholm boundary integral equations is here proposed through the explicit solutions of two Hilbert problems for holomorphic functions with the application of the conformal mappings. The stabiligy of the equilibrium of the annulus liquid layer is investigated by way of example.

---

Sommario Si prosegue lo studio di flussi piani viscosi con frontiera libera applicando l'approssimazione quasistazionaria [1]. L'introduzione della funzione stress-stream bianalitica consente l'uso della teoria delle funzioni olomorfe. La costruzione delle equazioni integrali di Fredholm al contorno proposta qui si basa sulla risoluzione esplicita di due problemi di Hilbert per funzioni analitiche mediante applicazione della tecnica delle trasformazioni conformi. Come esempio si studia la stabilità dell'equilibrio di uno strato liquido anulare.

     

Bézier curves in the modification of boundary integral equations (BIE) for potential boundary-values problems

The paper presents a modification of the classical boundary integral equation method (BIEM) for two-dimensional potential boundary values problems. The proposed modification consists in describing the boundary geometry by means of Bézier curves. As a result of this analytical modification of the BIEM, a new parametric integral equation system (PIES) was obtained. The kernels of these equations include the geometry of the boundary. This new PIES is no longer defined on the boundary, as in the case of the BIEM, but on the straight line for any given domain. The solution of the new PIES does not require a boundary discretization since it can be reduced merely to an approximation of boundary functions. To solve this PIES a pseudospectral method has been proposed and the results obtained were compared with exact solutions.     

Boundary element calculation of shallow water waves

A boundary element method is proposed for studying periodic shallow water problems. The numerical model is based on the shallow water equation. The key feature of this method is that the boundary integral equations are derived using the weighted residual method and the fundamental solutions for shallow water wave problems are obtained by solving the simultaneous singular equations. The accuracy of this method is studied for the wave reflection problem in a rectangular tank. As a result of this test, it has been shown that the number of element divisions and the distribution of nodes are significant to the accuracy. For numerical examples of external problems, the wave diffraction problems due to single cylindrical, double cylindrical and plate obstructions are analysed and compared with the exact and other numerical solutions. Relatively accurate solutions are obtained.     

边界积分方程法在动态断裂力学中的数值计算研究

本文将D.Nardini和C.A.Brebbia所提出的一种动态边界积分方程新解法应用于动态断裂力学数值计算，对数值实现问题，尤其对数值稳定性及精度问题进行了详细研究，得到了能保证数值稳定性的数值解法，给出了动态断裂力学计算实例，同己有的数值结果比较，表明本文的计算是成功的。     

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。     

A numerical method based on boundary integral equations and radial basis functions for plane anisotropic thermoelastostatic equations with general variable coefficients

A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coefficients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coefficients as long as all the requirements of the laws of physics are satisfied. To check the validity and accuracy of the proposed numerical method, some specific test problems with known solutions are solved.     

Variational principles for hydrodynamic impact problems

We first establish the rigorous field equations of the two continuous stages before and after entering water. Then correspondently, we obtain the specific variational principles, bounded theorems, and boundary integral equations of the second stage problems. The existence of solutions are proved and the scheme of solving the solutions are provided. Finally, as a numerical example, the ship's wave resistence problem is used to demonstrate the specific application of the second stage problems and its accuracy. Then we provide a rigorous and sound theoretical basis of variational finite element method and boundary element method for calculating the accurately fundamental equations.     

Nonsingular kernel boundary element method for thin-plate bending problems

In this paper,the nonsingular fundamental solutions are obtained from Fourierseries under some given conditions.These solutions can be taken as the kernels ofintegral equation.So a new boundary element method is presented,with which allkinds of thin-plate bending problems can’be solved,even with complicated loadings andsinuous boundaries.The calculation is much simpler and more accurate.     

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 a 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.     

Three-dimensional desingularized boundary integral methods for potential problems

The concept of desingularization in three-dimensional boundary integral computations is re-examined. The boundary integral equation is desingularized by moving the singular points away from the boundary and outside the problem domain. We show that the desingularization gives better solutions to several problems. As a result of desingularization, the surface integrals can be evaluated by simpler techniques, speeding up the computation. The effects of the desingularization distance on the solution and the condition of the resulting system of algebraic equations are studied for both direct and indirect versions of the boundary integral method. Computations show that a broad range of desingularization distances gives accurate solutions with significant savings in the computation time. The desingularization distance must be carefully linked to the mesh size to avoid problems with uniqueness and ill-conditioning. As an example, the desingularized indirect approach is tested on unsteady non-linear three-dimensional gravity waves generated by a moving submerged disturbance; minimal computational difficulties are encountered at the truncated boundary.     

A boundary element formulation using the simple method

A new boundary element method is described for calculation of the steady incompressible laminar flows. The method is based on the well-known SIMPLE algorithm. The new boundary element method allows one to find the fields of the pressure and velocity corrections without inner iterations, thus reducing the computational time drastically. This makes it different from the method developed by Patankar and Spalding.³² However, the new method demands a much larger computer strorage. The boundary integral equations are discretized with the help of constant boundary elements and constant cells. The values of the integrals along the boundary elements and the cells for the two-dimensional domain are found analytically. To preserve the stability in the iteration process, under-relaxation for the convection terms is used. This paper gives the results of calculations of the flows between two plane parallel plates at Re = 20 and Re = 200, the flows in a square cavity with a moving upper lid at Re = 1 and Re = 100 and the flow in a plane channel with sudden symmetric expansion at Re =46·6.     

Boundary integral equations for 2D elasticity and its application in discrete dislocation simulation in finite body: 2. Numerical implementation

In the previous paper by Yu and Diab (2013), several sets of boundary integral equations are derived for general anisotropic materials and corresponding equations for materials with different classes of symmetry are deduced. The work presented herein implements two sets of boundary element schemes to numerically solve the stress field. The integration on the element that has the singular point of the kernel is bounded and can be evaluated analytically. Four benchmark elastic problems are solved numerically to show the advantage of the two schemes over the conventional boundary element formulation in eliminating the boundary layer effect. The one with the weaker singularity has better convergence and gives more accurate results. The presented formulation also provides a direct approach to solve for stress field in a finite solid body in the presence of dislocations. Combined with discrete dislocations dynamics, boundary value problems with dislocations in finite bodies can be solved. Two examples, bending of a single crystal beam and pure shearing of a polycrystalline solid, are simulated by discrete dislocation dynamics using the scheme that has the weaker singularity. The comparisons with the published results using the well-established superposition technique validate the proposed formulation and show its quick convergence.