Necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions

A necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions is established in the present paper. A super-determined problem of harmonic functions is proposed for the first time. Then a necessary and sufficient condition for the existence of solution of the super-determined problem is proved. At the same time, it is a necessary and sufficient condition for the correct formulation of boundary integral equations with direct unknown quantities. A relation between boundary integral equations and variational principles is discovered for the first time. And a one-to-one correspondence between boundary integral equations with direct and indirect unknown quantities is indicated. The concept and route of the present paper can be applied to other boundary value problems possessing variational principles.     

Variational boundary integral approach for asymmetric impinging jets of arbitrary two‐dimensional nozzle

We consider asymmetric impinging jets issuing from an arbitrary nozzle. The flow is assumed to be two‐dimensional, inviscid, incompressible, and irrotational. The impinging jet from an arbitrary nozzle has a couple of separated infinite free boundaries, which makes the problem hard to solve. We formulate this problem using the stream function represented with a specific single layer potential. This potential can be extended to the surrounding region of the jet flow, and this extension can be proved to be a bounded function. Using this fact, the formulation yields the boundary integral equations on the entire nozzle and free boundary. In addition, a boundary perturbation produces an extraordinary boundary integral equation for the boundary variation. Based on these variational boundary integral equations, we can provide an efficient algorithm that can treat with the asymmetric impinging jets having arbitrarily shaped nozzles. Particularly, the proposed algorithm uses the infinite computational domain instead of a truncated one. To show the convergence and accuracy of the numerical solution, we compare our solutions with the exact solutions of free jets. Numerical results on diverse impinging jets with nozzles of various shapes are also presented to demonstrate the applicability and reliability of the algorithm.     

Volume preserving immersed boundary methods for two‐phase fluid flows

In this article, we propose a simple area‐preserving correction scheme for two‐phase immiscible incompressible flows with an immersed boundary method (IBM). The IBM was originally developed to model blood flow in the heart and has been widely applied to biofluid dynamics problems with complex geometries and immersed elastic membranes. The main idea of the IBM is to use a regular Eulerian computational grid for the fluid mechanics along with a Lagrangian representation of the immersed boundary. Using the discrete Dirac delta function and the indicator function, we can include the surface tension force, variable viscosity and mass density, and gravitational force effects. The principal advantage of the IBM for two‐phase fluid flows is its inherent accuracy due in part to its ability to use a large number of interfacial marker points on the interface. However, because the interface between two fluids is moved in a discrete manner, this can result in a lack of volume conservation. The idea of an area preserving correction scheme is to correct the interface location normally to the interface so that the area remains constant. Various numerical experiments are presented to illustrate the efficiency and accuracy of the proposed conservative IBM for two‐phase fluid flows. Copyright © 2011 John Wiley & Sons, Ltd.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

A new boundary integral equation for plane harmonic functions

In this paper, we propose a new boundary integral equation for plane harmonic functions. As a new approach, the equation is derived from the conservation integrals. Every variable in the integral equation has direct engineering interest. When this integral equation is applied to the Dirichlet problem, one will get an integral equation of the second kind, so that the algebraic equation system in the boundary element method has diagonal dominance. Finally, this equation is applied to elastic torsion problems of cylinders of different sections, and satisfactary numerical results are obtained.     

用面力边界积分方程求解断裂力学J积分

证明面力边界积分方程被积函数的散度等于零，应用Stokes公式，对平面线弹性问题，将面力边界积分的求解转化为边界点的位移势函数的点值计算。应用边界积分方程的求解结果，推导出J积分亦可表示为边界点的积分势函数的点值计算，无需进行数值积分，实例计算说明该方法的有效性。     

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.     

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.     

Two-dimensional polynomial eigenstrain formulation of boundary integral equation with numerical verification

The low-order polynomial-distributed eigenstrain formulation of the boundary integral equation (BIE) and the corresponding definition of the Eshelby tensors are proposed for the elliptical inhomogeneities in two-dimensional elastic media. Taking the results of the traditional subdomain boundary element method (BEM) as the control, the effectiveness of the present algorithm is verified for the elastic media with a single elliptical inhomogeneity. With the present computational model and algorithm, significant improvements are achieved in terms of the efficiency as compared with the traditional BEM and the accuracy as compared with the constant eigenstrain formulation of the BIE.     

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.     

弹性力学平面问题的无奇异边界积分方程

将弹性力学平面问题归化成无奇异边界积分方程,避免了传统的边界元法中的柯西主值（CPV）积分和Hadamard-Finite-Parts(HFP)积分的计算,建立完整的数值求解体系。     

热弹性平面问题的规则化边界积分方程

对于热弹性平面问题,过去广泛集中在直接变量边界元法研究,本文研究间接变量规则化边界元法,建立了间接变量规则化边界积分方程。和直接边界元法相比,间接法具有降低密度函数的连续性要求、位移梯度方程中的热载荷体积分具有较弱奇异性等优点。数值实施中,用精确单元描述边界几何,不连续插值函数逼近边界量。算例表明,本文方法效率高,所得数值结果与精确解相当吻合。     

Axisymmetric Powell-Eyring fluid flow with convective boundary condition: optimal analysis

The effects of axisymmetric flow of a Powell-Eyring fluid over an impermeable radially stretching surface are presented. Characteristics of the heat transfer process are analyzed with a more realistic condition named the convective boundary condition. Governing equations for the flow problem are derived by the boundary layer approximations. The modeled highly coupled partial differential system is converted into a system of ordinary differential equations with acceptable similarity transformations. The convergent series solutions for the resulting system are constructed and analyzed. Optimal values are obtained and presented in a numerical form using an optimal homotopy analysis method (OHAM). The rheological characteristics of different parameters of the velocity and temperature profiles are presented graphically. Tabular variations of the skin friction coefficient and the Nusselt number are also calculated. It is observed that the temperature distribution shows opposite behavior for Prandtl and Biot numbers. Furthermore, the rate of heating/cooling is higher for both the Prandtl and Biot numbers.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

A new velocity–vorticity boundary integral formulation for Navier–Stokes equations

In the present work, an indirect boundary integral method for the numerical solution of Navier–Stokes equations formulated in velocity–vorticity dependent variables is proposed. This wholly integral approach, based on Helmholtz's decomposition, deals directly with the vorticity field and gives emphasis to the establishment of appropriate boundary conditions for the vorticity transport equation. The coupling between the vorticity and the vortical velocity fields is expressed by an iterative procedure. The present analysis shows the usefulness of an integral formulation not only in providing a potentially more efficient computational tool, but also in giving a better understanding to the physics of the phenomenon. Copyright © 2000 John Wiley & Sons, Ltd.     

On the uniqueness of boundary integral equation for the exterior Helmholtz problem

From the point of view of energy analysis,the cause that the uniqueness of theboundary integral equation induced from the exterior Helmholtz problem does not hold isinvestigated in this paper.It is proved that the Sommerfeld’s condition at the infinity ischanged so that it is suitable not only for the radiative wave but also for the absorptive wavewhen we use the boundary integral equation to describe the exterior Helmholtz problem.Therefore,the total energy of the system is conservative.The mathematical dealings toguarantee the uniqueness are discussed based upon this explanation     

各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程。数值计算结果表明在退化尺度附近，充要的随机边界积分方程较习用的随机边界积分方程有较大的优越性。     

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.     

Development of temporal and spatial high‐order schemes for two‐fluid seven‐equation two‐pressure model and its applications in two‐phase flow benchmark problems

Current existing main nuclear thermal‐hydraulics (T‐H) system analysis codes, such as RALAP5, TRACE, and CATHARE, play a crucial role in the nuclear engineering field for the design and safety analysis of nuclear reactor systems. However, two‐fluid model used in these T‐H system analysis codes is ill posed, easily leading to numerical oscillations, and the classical first‐order methods for temporal and special discretization are widely employed for numerical simulations, yielding excessive numerical diffusion. Two‐fluid seven‐equation two‐pressure model is of particular interest due to the inherent well‐posed advantage. Moreover, high‐order accuracy schemes have also attracted great attention to overcome the challenge of serious numerical diffusion induced by low‐order time and space schemes for accurately simulating nuclear T‐H problems. In this paper, the semi‐implicit solution algorithm with high‐order accuracy in space and time is developed for this well‐posed two‐fluid model and the robustness and accuracy are verified and assessed against several important two‐phase flow benchmark tests in the nuclear engineering T‐H field, which include two linear advection problems, the oscillation problem of the liquid column, the Ransom water faucet problem, the reversed water faucet problem, and the two‐phase shock tube problem. The following conclusions are achieved. (1) The proposed semi‐implicit solution algorithm is robust in solving two‐phase flows, even for fast transients and discontinuous solutions. (2) High‐order schemes in both time and space could prevent excessive numerical diffusion effectively and the numerical simulation results are more accurate than those of first‐order time and space schemes, which demonstrates the advantage of using high‐order schemes.     

Pseudo‐characteristic formulation and dynamic boundary conditions for computational aeroacoustics

The present paper deals with the use of the pseudo‐characteristic formulation of the Navier–Stokes and Euler equations recently introduced by Sesterhenn (Comput. Fluid. 2001; 30 :37–67) for the simulation of acoustic wave propagation. The emphasis is put on the formulation of an efficient method on structured curvilinear grids, along with the definition and implementation of efficient boundary conditions. The cases of inflow, outflow, rigid/compliant walls and walls with prescribed impedance are addressed. The proposed boundary conditions are assessed on generic cases. The pseudo‐characteristic formulation enables a straightforward and optimal use of high‐order upwind dispersion‐relation‐preserving schemes, yielding an efficient method. Copyright © 2006 John Wiley & Sons, Ltd.