Direct method of solution for general boundary value problem of the Laplace equation

A general boundary value problem for two-dimensional Laplace equation in the domain enclosed by a piecewise smooth curve is considered. The Dirichlet and the Neumann data are prescribed on respective parts of the boundary, while there is the second part of the boundary on which no boundary data are given. There is the third part of the boundary on which the Robin condition is prescribed. This problem of finding unknown values along the whole boundary is ill posed. In this sense we call our problem an inverse boundary value problem. In order for a solution to be identified the inverse problem is reformulated in terms of a variational problem, which is then recast into primary and adjoint boundary value problems of the Laplace equation in its conventional form. A direct method for numerical solution of the inverse boundary value problem using the boundary element method is presented. This method proposes a non-iterative and unified treatment of conventional boundary value problem, the Cauchy problem, and under- or over-determined problems.     

Regularized MFS-Based Boundary Identification in Two-Dimensional Helmholtz-Type Equations

We study the stable numerical identification of an unknown portion of the boundary on which a given boundary condition is provided and additional Cauchy data are given on the remaining known portion of the boundary of a two-dimensional domain for problems governed by either the Helmholtz or the modified Helmholtz equation. This inverse geometric problem is solved using the method of fundamental solutions (MFS) in conjunction with the Tikhonov regularization method. The optimal value for the regularization parameter is chosen according to Hansen's L-curve criterion. The stability, convergence, accuracy and efficiency of the proposed method are investigated by considering several examples.     

The MFS for numerical boundary identification in two-dimensional harmonic problems

In this study, we briefly review the applications of the method of fundamental solutions to inverse problems over the last decade. Subsequently, we consider the inverse geometric problem of identifying an unknown part of the boundary of a domain in which the Laplace equation is satisfied. Additional Cauchy data are provided on the known part of the boundary. The method of fundamental solutions is employed in conjunction with regularization in order to obtain a stable solution. Numerical results are presented and discussed.     

The method of fundamental solutions for a biharmonic inverse boundary determination problem

In this paper, a nonlinear inverse boundary value problem associated to the biharmonic equation is investigated. This problem consists of determining an unknown boundary portion of a solution domain by using additional data on the remaining known part of the boundary. The method of fundamental solutions (MFS), in combination with the Tikhonov zeroth order regularization technique, are employed. It is shown that the MFS regularization numerical technique produces a stable and accurate numerical solution for an optimal choice of the regularization parameter. A. Zeb on study leave visiting the University of Leeds.     

Numerical modeling of micro- and macro-behavior of viscoelastic porous materials

This paper presents a numerical method for solving the two-dimensional problem of a polygonal linear viscoelastic domain containing an arbitrary number of non-overlapping circular holes of arbitrary sizes. The solution of the problem is based on the use of the correspondence principle. The governing equation for the problem in the Laplace domain is a complex hypersingular boundary integral equation written in terms of the unknown transformed displacements on the boundaries of the holes and the exterior boundaries of the finite body. No specific physical model is involved in the governing equation, which means that the method is capable of handling a variety of viscoelastic models. A truncated complex Fourier series with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the boundaries of the holes. A truncated complex series of Chebyshev polynomials with coefficients dependent on the transform parameter is used to approximate the unknown transformed displacements on the straight boundaries of the finite body. A system of linear algebraic equations is formed using the overspecification method. The viscoelastic stresses and displacements are calculated through the viscoelastic analogs of the Kolosov–Muskhelishvili potentials, and an analytical inverse Laplace transform is used to provide the time domain solution. Using the concept of representative volume, the effective viscoelastic properties of an equivalent homogeneous material are then found directly from the corresponding constitutive equations for the average field values. Several examples are given to demonstrate the accuracy of the method. The results for the stresses and displacements are compared with the numerical solutions obtained by commercial finite element software (ANSYS). The results for the effective properties are compared with those obtained with the self-consistent and Mori–Tanaka schemes.     

Numerical boundary identification for Helmholtz-type equations

We study the identification of an unknown portion of the boundary of a two-dimensional domain occupied by a material satisfying Helmholtz-type equations from additional Cauchy data on the remaining known portion of the boundary. This inverse geometric problem is approached using the boundary element method (BEM) in conjunction with the Tikhonov first-order regularization procedure, whilst the choice of the regularization parameter is based on the L-curve criterion. The numerical results obtained show that the proposed method produces a convergent and stable solution     

Indirect Trefftz method for boundary value problem of Poisson equation

Trefftz method is the boundary-type solution procedure using regular T-complete functions satisfying the governing equation. Until now, it has been mainly applied to numerical analyses of the problems governed with the homogeneous differential equations such as the two- and three-dimensional Laplace problems and the two-dimensional elastic problem without body forces. On the other hand, this paper describes the application of the indirect Trefftz method to the solution of the boundary value problems of the two-dimensional Poisson equation. Since the Poisson equation has an inhomogeneous term, it is generally difficult to determine the T-complete function satisfying the governing equation. In this paper, the inhomogeneous term containing an unknown function is approximated by a polynomial in the Cartesian coordinates to determine the particular solutions related to the inhomogeneous term. Then, the boundary value problem of the Poisson equation is transformed to that of the Laplace equation by using the particular solution. Once the boundary value problem of the Poisson equation is solved according to the ordinary Trefftz formulation, the solution of the boundary value problem of the Poisson equation is estimated from the solution of the Laplace equation and the particular solution. The unknown parameters included in the particular solution are determined by the iterative process. The present scheme is applied to some examples in order to examine the numerical properties.     

Numerical solutions of boundary detection problems using modified collocation Trefftz method and exponentially convergent scalar homotopy algorithm

In this paper, the boundary detection problem, which is governed by the Laplace equation, is analyzed by the modified collocation Trefftz method (MCTM) and the exponentially convergent scalar homotopy algorithm (ECSHA). In the boundary detection problem, the Cauchy data is given on part of the boundary and the Dirichlet boundary condition on the other part of the boundary, whose spatial position is unknown a priori. By adopting the MCTM, which is meshless and integral-free, the numerical solution is expressed by a linear combination of the T-complete functions of the Laplace equation. The use of a characteristic length in MCTM can stabilize the numerical procedure and ensure highly accurate solutions. Since the coefficients of MCTM and the position of part of the boundary are unknown, to collocate the boundary conditions will yield a system of nonlinear algebraic equations; the ECSHA, which is exponentially convergent, is adopted to solve the system of nonlinear algebraic equations. Several numerical examples are provided to demonstrate the ability and accuracy of the proposed meshless scheme. In addition, the consistency of the proposed scheme is validated by adding noise into the boundary conditions.     

A comparison of different methods to solve inverse biharmonic boundary value problems

The Boundary Element Method (BEM) is applied to solve numerically some inverse boundary value problems associated to the biharmonic equation which involve over‐ and under‐specified boundary portions of the solution domain. The resulting ill‐conditioned system of linear equations is solved using the regularization and the minimal energy methods, followed by a further application of the Singular Value Decomposition Method (SVD). The regularization method incorporates a smoothing effect into the least squares functional, whilst the minimal energy method is based on minimizing the energy functional for the Laplace equation subject to the linear constraints generated by the BEM discretization of the biharmonic equation. The numerical results are compared with known analytical solutions and the stability of the numerical solution is investigated by introducing noise into the input data. Copyright © 1999 John Wiley & Sons, Ltd.     

The Spring-Damping Regularization Method and the Lie-Group Shooting Method for Inverse Cauchy Problems

The inverse Cauchy problems for elliptic equations, such as the Laplace equation, the Poisson equation, the Helmholtz equation and the modified Helmholtz equation, defined in annular domains are investigated. The outer boundary of the annulus is imposed by overspecified boundary data, and we seek unknown data on the inner boundary through the numerical solution by a spring-damping regularization method and its Lie-group shooting method (LGSM). Several numerical examples are examined to show that the LGSM can overcome the ill-posed behavior of inverse Cauchy problem against the disturbance from random noise, and the computational cost is very cheap.     

Sensitivity analysis scheme of boundary value problem of 2D Poisson equation by using Trefftz method

This paper describes the sensitivity analysis of the boundary value problem of two-dimensional Poisson equation by using Trefftz method. A non-homogeneous term of two-dimensional Poisson equation is approximated with a polynomial function in the Cartesian coordinates to derive a particular solution. The unknown function of the boundary value problem is approximate with the superposition of the T-complete functions of Laplace equation and the derived particular solution with unknown parameters. The parameters are determined so that the approximate solution satisfies boundary conditions. Since the T-complete functions and the particular solution are regular, direct differentiation of the expression results leads to the sensitivity expressions. The boundary-specified value and the shape parameter are taken as the variables of the sensitivity analysis to formulate the sensitivity analysis methods. The present scheme is applied to some numerical examples in order to confirm the validity of the present algorithm.     

Boundary integral equation solution of moving boundary phase change problems

Boundary integral equation methods are presented for the solution of some two-dimensional phase change problems. Convection may enter through boundary conditions, but cannot be considered within phase boundaries. A general formulation based on space-time Green's functions is developed using the complete heat equation, followed by a simpler formulation using the Laplace equation. The latter is pursued and applied in detail. An elementary, noniterative system is constructed, featuring linear interpolation over elements on a polygonal boundary. Nodal values of the temperature gradient normal to a phase change boundary are produced directly in the numerical solution. The system performs well against basic analytical solutions, using these values in the interphase jump condition, with the simplest formulation of the surface normal at boundary vertices. Because the discretized surface changes automatically to fit the scale of the problem, the method appears to offer many of the advantages of moving mesh finite element methods. However, it only requires the manipulation of a surface mesh and solution for surface variables. In some applications, coarse meshes and very large time steps may be used, relative to those which would be required by fixed grid domain methods. Computations are also compared to original lab data, describing two-dimensional soil freezing with a time-dependent boundary condition. Agreement between simulated and measured histories is good.     

Potential discontinuity simulations with the numerical Green’s function in steady and transient problems

This paper applies the numerical Green’s function (NGF) boundary element formulation (BEM) first in standard form to solve the Laplace equation and then, coupled to the operational quadrature method (OQM), to solve time domain problems (TD-BEM). Both involve the analysis of potential discontinuities in the respective scalar model simulation. The implementation of the associated Green’s function acting as the fundamental solution is advantageous since element discretization of actual discontinuity surfaces are no longer required. In the OQM the convolution integral is substituted by a quadrature formula, whose weights are computed using the fundamental solution in the Laplace domain, producing the direct solution to the problem in the time domain. Applications of the NGF to problems involving the Laplace equation and its transient counterpart are presented for two-dimensional potential flow examples, confirming that the formulation is stable and accurate.     

Identification of a Corroded Boundary and Its Robin Coefficient

An inverse geometric problem for two-dimensional Helmholtz-type equations arising in corrosion detection is considered. This problem involves determining an unknown corroded portion of the boundary of a two-dimensional domain and possibly its surface heat transfer (impedance) Robin coefficient from one or two pairs of boundary Cauchy data (boundary temperature and heat flux), and is solved numerically using the meshless method of fundamental solutions. A nonlinear unconstrained minimisation of the objective function is regularised when noise is added into the input boundary data. The stability of the numerical results is investigated for several test examples, with respect to noise in the input data and various values of the regularisation parameters.     

Application of the boundary element method to three-dimensional potential problems in heterogeneous media

The present work discusses a solution procedure for heterogeneous media three-dimensional potential problems, involving nonlinear boundary conditions. The problem is represented mathematically by the Laplace equation and the adopted numerical technique is the boundary element method (BEM), here using velocity correcting fields to simulate the conductivity variation of the domain. The integral equation is discretized using surface elements for the boundary integrals and cells, for the domain integrals. The adopted strategy subdivides the discretized equations in two systems: the principal one involves the calculation of the potential in all boundary nodes and the secondary which determines the correcting field of the directional derivatives of the potential in all points. Comparisons with other numerical and analytical solutions are presented for some examples.     

Identification of injection and extraction wells from overspecified boundary data

In this paper, we focus on the identification of wells’ positions and fluxes/flows from the knowledge of overspecified data: hydraulic head and flux, on a part of the domain boundary. The used method is based on minimizing a constitutive law gap functional. We consider two inverse problems: in the first one overspecified conditions are available throughout the entire domain boundary; in the second inverse problem, in addition to the wells, boundary condition are also unknown on an inaccessible part of the domain boundary.     

Numerical operational methods for time-dependent linear problems

A general and systematic discussion on the use of the operational method of Laplace transform for numerically solving complex time-dependent linear problems is presented. Application of Laplace transform with respect to time on the governing differential equations as well as the boundary and initial conditions of the problem reduces it to one independent of time, which is solved in the transform domain by any convenient numerical technique, such as the finite element method, the finite difference method or the boundary integral equation method. Finally, the time domain solution is obtained by a numerical inversion of the transformed solution. Eight existing methods of numerical inversion of the Laplace transform are systematically discussed with respect to their use, range of applicability, accuracy and computational efficiency on the basis of some framework vibration problems. Other applications of the Laplace transform method in conjunction with the finite element method or the boundary integral equation method in the areas of earthquake dynamic response of frameworks, thermaliy induced beam vibrations, forced vibrations of cylindrical shells, dynamic stress concentrations around holes in plates and viscoelastic stress analysis are also briefly described to demonstrate the generality and advantages of the method against other known methods.     

An alternative method for degenerate scale problems in boundary element methods for the two-dimensional Laplace equation

The boundary integral equation approach has been shown to suffer a nonunique solution when the geometry is equal to a degenerate scale for a potential problem. In this paper, the degenerate scale problem in boundary element method for the two-dimensional Laplace equation is analytically studied in the continuous system by using degenerate kernels and Fourier series instead of using discrete system using circulants [Engng Anal. Bound. Elem. 25 (2001) 819]. For circular and multiply-connected domain problems, the rank-deficiency problem of the degenerate scale is solved by using the combined Helmholtz exterior integral equation formulation (CHEEF) concept. An additional constraint by collocating a point outside the domain is added to promote the rank of influence matrix. Two examples are shown to demonstrate the numerical instability using the singular integral equation for circular and annular domain problems. The CHEEF concept is successfully applied to overcome the degenerate scale and the error is suppressed in the numerical experiment.     

A formulation and solution for boundary element analysis of inhomogeneous-nonlinear problem

A new formulation is presented in this paper for the boundary element analysis of a nonlinear potential-type problem wherein the linear term is governed by the Laplace operator, and the nonlinear term is a function of the spatial coordinates as well as the unknown solution function. The formulation aims to transform the domain integral relevant to the inhomogeneous-nonlinear term to a corresponding boundary integral. The proposed approach is different from the more popular schemes for the purpose, such as the Dual Reciprocity and Multiple Reciprocity Methods. The inhomogeneous-nonlinear term is first approximated by a polynomial in terms of the space coordinates with unknown coefficients. Integral equations on the selected points (referred to “computing points”) on the boundary as well as inside domain are employed to determine the above-mentioned unknown coefficients using the least square method. The number of computing points affects the accuracy of the result, which is discussed through some numerical examples in two-dimensional space.     

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.