The analytic continuation of solutions to elliptic boundary value problems in two independent variables

An explicit representation is derived for the continuation across an analytic boundary of the solution to a boundary value problem for an analytic elliptic equation of second order in two independent variables. The representation is in terms of Cauchy data on the boundary and the complex Riemann function. This is equivalent to a representation for the solution to Cauchy's problem given by Henrici in 1957. It is confirmed that the method of complex characteristics is satisfactory for locating real singularities in the solution provided that the Riemann function is entire in its four arguments. Applications to Laplace's and Helmholtz's equations are discussed. By inserting known, simple solutions to the latter equation into the representation formula, several nontrivial integral relations involving the Bessel function J₀, and a possibly new series expansion for J_μ(x), are found.     

A subregion DRBEM formulation for the dynamic analysis of two-dimensional cracks

The dual reciprocity boundary element method employing the step by step time integration technique is developed to analyse two-dimensional dynamic crack problems. In this method the equation of motion is expressed in boundary integral form using elastostatic fundamental solutions. In order to transform the domain integral into an equivalent boundary integral, a general radial basis function is used for the derivation of the particular solutions. The dual reciprocity boundary element method is combined with an efficient subregion boundary element method to overcome the difficulty of a singular system of algebraic equations in crack problems. Dynamic stress intensity factors are calculated using the discontinuous quarter-point elements. Several examples are presented to show the formulation details and to demonstrate the computational efficiency of the method.     

Application of Wiener-Hopf technique to linear nonhomogeneous integral equations for a new representation of Chandrasekhar's H-function in radiative transfer, its existence and uniqueness

In this paper, the linear nonhomogeneous integral equation of H-functions is considered to find a new form of H-function as its solution. The Wiener-Hopf technique is used to express a known function into two functions with different zones of analyticity. The linear nonhomogeneous integral equation is thereafter expressed into two different sets of functions having the different zones of regularity. The modified form of Liouville's theorem is thereafter used, Cauchy's integral formulae are used to determine functional representation over the cut region in a complex plane. The new form of H-function is derived both for conservative and nonconservative cases. The existence of solution of linear nonhomogeneous integral equations and its uniqueness are shown. For numerical calculation of this new H-function, a set of useful formulae are derived both for conservative and nonconservative cases.     

A Symmetric Galerkin Boundary Element Method for Linear Poroelasticity

In linear poroelasticity so far only collocation boundary element methods have been available. However, in some applications, e.g., when coupling with finite elements is desired, a symmetric formulation is preferable. Choosing a Galerkin approach which involves the second boundary integral equation, such a formulation is possible. Here, a previously presented integration by part technique for the regularization of the first boundary integral equation is extended to the second boundary integral equation as well. While the weakly singular representation of the double layer operator has been presented before, the emphasis lies here on the so called hyper-singular boundary integral operator. Due to the regularization, this operator can be evaluated numerically and, hence, be used within a numerical scheme for the first time. Different numerical studies will be presented to show the behavior of the established symmetric Galerkin boundary element method, also comparing it with collocation boundary element methods. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hammerstein Integral Equations with Nontrivial Solutions

By means of critical point theory, existence theorems for nontrivial solutions to the Hammerstein equation x = KFx are given, where K is a compact linear integral operator and F is a nonlinear superposition operator. To this end, appropriate conditions on the spectrum of the linear parte are combined with growth and representation conditions on the nonlinear part to ensure the applicability of the mountain — pass lemma. The abstract existence theorems are applied to nonlinear elliptic equations and systems subject to Dirichlet boundary conditions.     

Generalized Dirichlet‐to‐Neumann Map in Time‐Dependent Domains

We study the heat, linear Schrödinger (LS), and linear KdV equations in the domain l(t) < x < ∞ , 0 < t < T , with prescribed initial and boundary conditions and with l(t) a given differentiable function. For the first two equations, we show that the unknown Neumann or Dirichlet boundary value can be computed as the solution of a linear Volterra integral equation with an explicit weakly singular kernel. This integral equation can be derived from the formal Fourier integral representation of the solution. For the linear KdV equation we show that the two unknown boundary values can be computed as the solution of a system of linear Volterra integral equations with explicit weakly singular kernels. The derivation in this case makes crucial use of analyticity and certain invariance properties in the complex spectral plane. The above Volterra equations are shown to admit a unique solution.     

A fast boundary element algorithm for time-dependant potential problems

The numerical solution of time-dependant potential problems via the boundary element has been crippled by the high computational cost due to the inherent time history constraint in the integral representation. Using a boundary-only formulation, the time integrations, at any instant in time, have to be evaluated starting from the initial time. This time-history dependence becomes impractical and inadequate for problems where computations are to be performed for extended times. This also made the boundary element uncompetitive compared to the domain-mesh based methods, such as finite difference and finite element methods, for the solution of transient potential problems. Generally, the evaluation of the potential at N domain points using M boundary points at the Kth time step requires an amount of computer operations of the order O(KM²+KNM). This paper presents an algorithm which requires a computational cost of the order of only O(M²+NM), where the dependence from the past K-steps is removed. The algorithm combines the boundary element method and a scheme, which uses virtual collocation points and radial basis functions to approximate the domain integral.     

The Poisson and the Biharmonic Equations in the Interior of a Convex Polygon

A new method for analyzing linear elliptic partial differential equations in the interior of a convex polygon was developed in the late 1990s. This method does not rely on the classical approach of separation of variables and on the use of classical integral transforms and therefore is well suited for the investigation of the biharmonic equation. Here, we present a novel integral representation of the solution of the biharmonic equation in the interior of a convex polygon. This representation contains certain free parameters and therefore is more general than the one presented in [1]. For a given boundary value problem, by choosing these free parameters appropriately, one can obtain the simplest possible representation for the solution. This representation still involves certain unknown boundary values, thus for this formula to become effective it is necessary to characterize the unknown boundary values in terms of the given boundary conditions. This requires the investigation of certain relations refereed to as the global relations. A general approach for analyzing these relations is illustrated by solving several problems formulated in the interior of a semistrip. In addition, for completeness, similar results are presented for the Poisson equation by employing an integral representation for the Laplace equation which is more general than the one derived in the late 1990s.     

Boundary element method and wave equation

In this paper the boundary integral expression for a one-dimensional wave equation with homogeneous boundary conditions is developed. This is done using the time dependent fundamental solution of the corresponding hyperbolic partial differential equation. The integral expression developed is a generalized function with the same form as the well-known D'Alembert formula. The derivatives of the solution and some useful invariants on the characteristics of the partial differential equation are also calculated.The boundary element method is applied to find the numerical solution. The results show excellent agreement with analytical solutions.A multi-step procedure for large time steps which can be used in the boundary element method is also described.In addition, the way in which boundary conditions are introduced during the time dependent process is explained in detail. In the Appendix the main properties of Dirac's delta function and the Heaviside unit step function are described.     

Laplace方程边值问题的边界积分方程法

§ 1.　Introduction　　Inengineeringandtechnology ,theproblemofstaticelectricfieldscanbeattributedtotheboundaryproblemofLaplaceequationofstaticeletricpotentialfunction .Themethodsofclassi calmathematicalphysicscanbeonlyusedtosolveboundaryproblemofverysimpledomainandspecialboundarycondition .Althoughthemethodsoflimitedelementscanbeusedtosolvetheproblemsonarbitrarydomain ,butitneedstopartitionthewholedomainandtocalculateverycomplex .Theapproachofboundaryintegralequationistosolverelatedproblemsb…     

Boundary element method for spatially-periodic potential problems

In this paper, we propose a boundary element method for two-dimensional potential problems with one-dimensional spatial periodicity, which have been difficult to be solved by the ordinary boundary element method. In the presented method, we reduce the potential problems with Dirichlet and Neumann boundary conditions to integral equation problems with the periodic fundamental solution of the Laplace operator and, then, obtain approximate solutions by solving linear systems given by discretizing the integral equations. Numerical examples are also included. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Linear integral equations for conformal mapping of bounded multiply connected regions onto a disk with circular slits

In this paper we present a boundary integral equation method for the numerical conformal mapping of bounded multiply connected region Ω onto a disk with circular slits. The method is based on some uniquely solvable boundary integral equations with classical adjoint and generalized Neumann kernels. These boundary integral equations are constructed from a boundary relationship satisfied by a function analytic on a multiply connected region. Some numerical examples are presented to illustrate the efficiency of the presented method.     

Local mass conservative Hermite interpolation for the solution of flow problems by a multi-domain boundary element approach

This paper presents a local Hermite radial basis function interpolation scheme for the velocity and pressure fields. The interpolation for velocity satisfies the continuity equation (mass conservative interpolation) while the pressure interpolation obeys the pressure equation. Additionally, the Dual Reciprocity Boundary Element method (DRBEM) is applied to obtain an integral representation of the Navier-Stokes equations. Then, the proposed local interpolation is used to obtain the values of the field variables and their partial derivatives at the boundary of the sub-domains. This interpolation allows one to obtain the boundary values needed for the integral formulas for velocity and pressure at some nodes within the sub-domains. In the proposed approach the boundary elements are merely used to parameterize the geometry, but not for the evaluation of the integrals as it is usually done. The presented multi-domain approach is different from the traditional ones in boundary elements because the resulting integral equations are non singular and the boundary data needed for the boundary integrals are approximated using a local interpolation. Some accurate results for simple Stokes problems and for the Navier-Stokes equations at low Reynolds numbers up to Re = 400 were obtained.     

Dual bivariational principles for linear problems

This paper presents dual bivariational principles which yield upper and lower bounds for 〈g, φ〉, where g is an arbitrary function and φ is the solution of the linear equation Aφ = f with general mixed boundary conditions. Variational principles associated with 〈f, φ〉 are taken as the starting-point, and the results generalize those of recent authors for linear integral equations.     

Fast integral equation methods for the Laplace-Beltrami equation on the sphere

Integral equation methods for solving the Laplace-Beltrami equation on the unit sphere in the presence of multiple “islands” are presented. The surface of the sphere is first mapped to a multiply-connected region in the complex plane via a stereographic projection. After discretizing the integral equation, the resulting dense linear system is solved iteratively using the fast multipole method for the 2D Coulomb potential in order to calculate the matrix-vector products. This numerical scheme requires only O(N) operations, where N is the number of nodes in the discretization of the boundary. The performance of the method is demonstrated on several examples.     

A Class of Quadratic Integral-Differential Equations

A nonlinear integro-ordinary differential equation built up by a linear ordinary differential operator of n th order with constant coefficients and a quadratic integral term is dealt with. The integral term represents the so-called autocorrelation of the unknown function. Applying the Fourier cosine transformation, the integral-differential equation is reduced to a quadratic boundary value problem for the complex Fourier transform of the solution in the upper half-plane. This problem in turn is reduced to a linear boundary value problem which can be solved in closed form. There are infinitely many solutions of the integral-differential equation depending on the prescribed zeros of a function related to the complex Fourier transform.     

Numerical evaluation of various levels of singular integrals, arising in BEM and its application in hydrofoil analysis

The sound implementation of the boundary element method (BEM) is highly dependent on an accurate numerical integration of singular integrals. In this paper, a set of various types of singular domain integrals with three-dimensional boundary element discretization is evaluated based on a transformation integration technique. In the BEM, the integration domain (body surface) needs to be discretized into small elements. For each element, the integral I(x^p, x) is calculated on the domain dS. Several types of integrals IB_α and IC_α are numerically and analytically computed and compared with the relative error. The method is extended to evaluate singular integrals which arise in the solution of the three-dimensional Laplace’s equation. An example of the elliptic hydrofoil is performed to study the physical accuracy. The results obtained using both numerical and analytical methods are shown in good agreement with the experimental data.     

On the coupling of BEM and FEM for exterior problems for the Helmholtz equation

This paper deals with the coupled procedure of the boundary element method (BEM) and the finite element method (FEM) for the exterior boundary value problems for the Helmholtz equation. A circle is selected as the common boundary on which the integral equation is set up with Fourier expansion. As a result, the exterior problems are transformed into nonlocal boundary value problems in a bounded domain which is treated with FEM, and the normal derivative of the unknown function at the common boundary does not appear. The solvability of the variational equation and the error estimate are also discussed.

     

A BEM for transient thermoelastic fracture analysis of FGMs

A boundary element method for the transient thermoelastic fracture analysis in isotropic, continuously non-homogeneous and linear elastic functionally graded materials subjected to a thermal shock is presented. The material parameters are assumed to be continuous functions of the Cartesian coordinates. Laplace-domain fundamental solutions of linear coupled thermoelasticity for infinite, isotropic, homogeneous and linear elastic solids are applied to derive the boundary-domain integral equation formulation. The numerical implementation is performed by using a collocation method for the spatial discretization. Numerical results for the dynamic stress intensity factors are presented and discussed. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)