A projection iterative algorithm boundary element method for the Signorini problem

We propose a projection iterative algorithm based on a fixed point equation for solving a certain class of Signorini problem. The satisfaction of the Signorini boundary conditions is verified in a projection iterative manner, and at each iterative step, an elliptic mixed boundary value problem is solved by a boundary element method which is suitable for any domain. We prove the convergence of the algorithm by the property of projection. The advantage of this algorithm is that it is easy to be implemented and converge quickly. Some numerical results show the accuracy and effectiveness of the algorithm.     

A coupled complex boundary method for the Cauchy problem

Considered in this paper is a Cauchy problem governed by an elliptic partial differential equation. In the Cauchy problem, one wants to recover the unknown Neumann and Dirichlet data on a part of the boundary from the measured Neumann and Dirichlet data, usually contaminated with noise, on the remaining part of the boundary. The Cauchy problem is an inverse problem with severe ill-posedness. In this paper, a coupled complex boundary method (CCBM), originally proposed in [Cheng XL, Gong RF, Han W, et al. A novel coupled complex boundary method for solving inverse source problems. Inverse Prob. 2014;30:055002], is applied to solve the Cauchy problem stably. With the CCBM, all the data, including the known and unknown ones on the boundary are used in a complex Robin boundary on the whole boundary. As a result, the Cauchy problem is transferred into a complex Robin boundary problem of finding the unknown data such that the imaginary part of the solution equals zero in the domain. Then the Tikhonov regularization is applied to the resulting new formulation. Some theoretical analysis is performed on the CCBM-based Tikhonov regularization framework. Moreover, through the adjoint technique, a simple solver is proposed to compute the regularized solution. The finite-element method is used for the discretization. Numerical results are given to show the feasibility and effectiveness of the proposed method.     

The boundary element method applied to a moving free boundary problem

In this paper a boundary problem is considered for which the boundary is to be determined as part of the solution. A time‐dependent problem involving linear diffusion in two spatial dimensions which results in a moving free boundary is posed. The fundamental solution is introduced and Green’s Theorem is used to yield a non‐linear system of integral equations for the unknown solution and the location of the boundary. The boundary element method is used to obtain a numerical solution to this system of integral equations which in turn is used to obtain the solution of the original problem. Graphical results for a two‐dimensional problem are presented. Published in 1999 by John Wiley & Sons, Ltd.     

The singular function boundary integral method for a two-dimensional fracture problem

The singular function boundary integral method (SFBIM) originally developed for Laplacian problems with boundary singularities is extended for solving two-dimensional fracture problems formulated in terms of the Airy stress function. Our goal is the accurate, direct computation of the associated stress intensity factors, which appear as coefficients in the asymptotic expansion of the solution near the crack tip. In the SFBIM, the leading terms of the asymptotic solution are used to approximate the solution and to weight the governing biharmonic equation in the Galerkin sense. The discretized equations are reduced to boundary integrals by means of Green's theorem and the Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers. The numerical results on a model problem show that the method converges extremely fast and yields accurate estimates of the leading stress intensity factors.     

Modified nonlocal boundary value problem method for an ill-posed problem for the biharmonic equation

In this paper, we propose a modified nonlocal boundary value problem method for an homogeneous biharmonic equation in a rectangular domain. We show that the considered problem is ill-posed in the sense of Hadamard, i.e. the solution does not depend continuously on the given data. Convergence estimates for the regularized solution are obtained under a priori bound assumptions for the exact solution. Some numerical results are given to show the effectiveness of the proposed method.     

A boundary integral equation method for a Neumann boundary problem for force-free fields

Summary A Neumann boundary value problem for the equation rot –=0 is considered in 29-1 and 29-2. The approach is by transforming the boundary value problem into an equivalent boundary integral equation deduced from a representation formula for solutions of rot –=0 based on the fundamental solution of the Helmholtz equation. In particular, for the two-dimensional case a detailed discussion of the integral equation is carried out including the approximate solution by numerical integration.     

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.     

A computational method for inverse free boundary determination problem

Based on the method of fundamental solutions and discrepancy principle for the choice of location for source points, we extend in this paper the application of the computational method to determine an unknown free boundary of a Cauchy problem of parabolic‐type equation from measured Dirichlet and Neumann data with noises. The standard Tikhonov regularization technique with the L‐curve method for an optimal regularized parameter is adopted for solving the resultant highly ill‐conditioned system of linear equations. Both one‐dimensional and two‐dimensional numerical examples are given to verify the efficiency and accuracy of the proposed computational method. Copyright © 2007 John Wiley & Sons, Ltd.     

A fast collocation method for an inverse boundary value problem

In this paper, we present an implementation of a fast multiscale collocation method for boundary integral equations of the second kind, and its application to solving an inverse boundary value problem of recovering a coefficient function from a boundary measurement. We illustrate by numerical examples the insensitive nature of the map from the coefficient to measurement, and design and test a Gauss–Newton iteration algorithm for obtaining the best estimate of the unknown coefficient from the given measurement based on a least‐squares formulation. Copyright © 2004 John Wiley & Sons, Ltd.     

A multipole expansion-based boundary element method for axisymmetric potential problem

The multipole expansion is an approximation technique used to evaluate the potential field due to sources located in the far field. Based on the multipole expansion, we describe a new technique to calculate the far potential field due to ring sources which are encountered in the boundary element method (BEM) formulation of axisymmetric problems. As the sources in the near field are processed by the slower conventional BEM, it is important to maximize the amount of multipole calculations taking advantage of both interior and exterior multipole expansions. Numerical results are presented for an axisymmetric potential test problem with Neumann and Dirichlet boundary conditions. The complexity of the proposed method remains O(N²), which is equal to that of the conventional BEM. However, the proposed technique coupled with an iterative solver speeds up the solution procedure. The technique is significantly advantageous when medium and large numbers of elements are present in the domain.     

A local boundary integral-based meshless method for Biot's consolidation problem

Traditional numerical techniques such as FEM and BEM have been successfully applied to the solutions of Biot's consolidation problems. However, these techniques confront some difficulties in dealing with moving boundaries. In addition, pre-designing node connectivity or element is not an easy task. Recently, developed meshless methods may overcome these difficulties. In this paper, a meshless model, based on the local Petrov–Galerkin approach with Heaviside step function as well as radial basis functions, is developed and implemented for the numerical solution of plane strain poroelastic problems. Although the proposed method is based on local boundary integral equation, it does not require any fundamental solution, thus avoiding the singularity integral. It also has no domain integral over local domain, thus largely reducing the computational cost in formulation of system stiffness. This is a truly meshless method. The solution accuracy and the code performance are evaluated through one-dimensional and two-dimensional consolidation problems. Numerical examples indicate that this meshless method is suitable for either regular or irregular node distributions with little loss of accuracy, thus being a promising numerical technique for poroelastic problems.     

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.     

Two-model iteration method for the solution of an inverse boundary heat-exchange problem

A method is proposed for restoration of the boundary condition, with an iterative correction of the initial data used in this method, involving the utilization of both exact and approximate heat-transfer models.Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 56, No. 2, pp. 313–319, February, 1989.     

A contact algorithm for the Signorini problem using space–time finite elements

The most common approach in the finite‐element modelling of continuum systems over space and time is to employ the finite‐element discretization over the spatial domain to reduce the problem to a system of ordinary differential equations in time. The desired time integration scheme can then be used to step across the so‐called time slabs, mesh configurations in which every element shares the same degree of time refinement. These techniques may become inefficient when the nature of the initial boundary value problem is such that a high degree of time refinement is required only in specific spatial regions of the mesh. Ideally one would be able to increase the time refinement only in those necessary regions. We achieve this flexibility by employing space–time elements with independent interpolation functions in both space and time. Our method is used to examine the classic contact problem of Signorini and allows us to increase the time refinement only in the spatial region adjacent to the contact interface. We also develop an interface‐tracking algorithm that tracks the contact boundary through the space–time mesh and compare our results with those of Hertz contact theory. Copyright 2004 John Wiley & Sons, Ltd.     

Signorini’s method for live loads and second-order effects

In this paper we deal with Signorini’s method for live loads to derive approximate solutions of some pure traction boundary value problems in finite elastostatics. Then, the obtained solutions are used to determine the second-order constitutive constants of a homogeneous and isotropic elastic material.     

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.     

A fast multipole boundary element method for solving the thin plate bending problem

A fast multipole boundary element method (BEM) for solving large-scale thin plate bending problems is presented in this paper. The method is based on the Kirchhoff thin plate bending theory and the biharmonic equation governing the deflection of the plate. First, the direct boundary integral equations and the conventional BEM for thin plate bending problems are reviewed. Second, the complex notation of the kernel functions, expansions and translations in the fast multipole BEM are presented. Finally, a few numerical examples are presented to show the accuracy and efficiency of the fast multipole BEM in solving thin plate bending problems. The bending rigidity of a perforated plate is evaluated using the developed code. It is shown that the fast multipole BEM can be applied to solve plate bending problems with good accuracy. Possible improvements in the efficiency of the method are discussed.     

Solution of the Helmholtz eigenvalue problem via the boundary element method

The numerical solution of the Helmholtz eigenvalue problem is considered. The application of the boundary element method reduces it to that of a non-linear eigenvalue problem. Through a polynomial approximation with respect to the wavenumber, the non-linear eigenvalue problem is reduced to a standard generalized eigenvalue problem. The method is applied to the test problems of a three-dimensional sphere with an axisymmetric boundary condition and a two-dimensional square.