Trefftz boundary elements—multi‐region formulations

This paper is concerned with an effective numerical implementation of the Trefftz boundary element method, for the analysis of two‐dimensional potential problems, defined in arbitrarily shaped domains. The domain is first discretized into multiple subdomains or regions. Each region is treated as a single domain, either finite or infinite, for which a complete set of solutions of the problem is known in the form of an expansion with unknown coefficients. Through the use of weighted residuals, this solution expansion is then forced to satisfy the boundary conditions of the actual domain of the problem, leading thus to a system of equations, from which the unknowns can be readily determined. When this basic procedure is adopted, in the analysis of multiple‐region problems, proper boundary integral equations must be used, along common region interfaces, in order to couple to each other the unknowns of the solution expansions relative to the neighbouring regions. These boundary integrals are obtained from weighted residuals of the coupling conditions which allow the implementation of any order of continuity of the potential field, across the interface boundary, between neighbouring regions. The technique used in the formulation of the region‐coupling conditions drives the performance of the Trefftz boundary element method. While both of the collocation and Galerkin techniques do not generate new unknowns in the problem, the technique of Galerkin presents an additional and unique feature: the size of the matrix of the final algebraic system of equations which is always square and symmetric, does not depend on the number of boundary elements used in the discretization of both the actual and region‐interface boundaries. This feature which is not shared by other numerical methods, allows the Galerkin technique of the Trefftz boundary element method to be effectively applied to problems with multiple regions, as a simple, economic and accurate solution technique. A very difficult example is analysed with this procedure. The accuracy and efficiency of the implementations described herein make the Trefftz boundary element method ideal for the study of potential problems in general arbitrarily‐shaped two‐dimensional domains. Copyright © 1999 John Wiley & Sons, Ltd.     

The boundary element solution of the Laplace and biharmonic equations subjected to noisy boundary data

This study investigates the numerical solution of the Laplace and biharmonic equations subjected to noisy boundary data. Since both equations are linear, they are numerically discretized using the Boundary Element Method (BEM), which does not use any solution domain discretization, to reduce the problem to solving a system of linear algebraic equations for the unspecified boundary values. It is shown that when noisy, lower-order derivatives are prescribed on the boundary, then a direct approach, e.g. Gaussian elimination, for solving the resulting discretized system of linear equations produces an unstable, i.e. unbounded and highly oscillatory, numerical solution for the unspecified higher-order boundary derivatives data. In order to overcome this difficulty, and produce a stable solution of the resulting system of linear equations, the singular value decomposition approach (SVD), truncated at an optimal level given by the L-curve method, is employed. © 1998 John Wiley & Sons, Ltd.     

Iterative solution of Hermite boundary integral equations

An efficient iterative method for solution of the linear equations arising from a Hermite boundary integral approximation has been developed. Along with equations for the boundary unknowns, the Hermite system incorporates equations for the first‐order surface derivatives (gradient) of the potential, and is therefore substantially larger than the matrix for a corresponding linear approximation. However, by exploiting the structure of the Hermite matrix, a two‐level iterative algorithm has been shown to provide a very efficient solution algorithm. In this approach, the boundary function unknowns are treated separately from the gradient, taking advantage of the sparsity and near‐positive definiteness of the gradient equations. In test problems, the new algorithm significantly reduced computation time compared with iterative solution applied to the full matrix. This approach should prove to be even more effective for the larger systems encountered in three‐dimensional analysis, and increased efficiency should come from pre‐conditioning of the non‐sparse matrix component. Published in 2007 by John Wiley & Sons, Ltd.     

Local discretization error bounds using interval boundary element method

In this paper, a method to account for the point‐wise discretization error in the solution for boundary element method is developed. Interval methods are used to enclose the boundary integral equation and a sharp parametric solver for the interval linear system of equations is presented. The developed method does not assume any special properties besides the Laplace equation being a linear elliptic partial differential equation whose Green's function for an isotropic media is known. Numerical results are presented showing the guarantee of the bounds on the solution as well as the convergence of the discretization error. Copyright © 2008 John Wiley & Sons, Ltd.     

Accurate radiation boundary conditions for the time‐dependent wave equation on unbounded domains

Asymptotic and exact local radiation boundary conditions (RBC) for the scalar time‐dependent wave equation, first derived by Hagstrom and Hariharan, are reformulated as an auxiliary Cauchy problem for each radial harmonic on a spherical boundary. The reformulation is based on the hierarchy of local boundary operators used by Bayliss and Turkel which satisfy truncations of an asymptotic expansion for each radial harmonic. The residuals of the local operators are determined from the solution of parallel systems of linear first‐order temporal equations. A decomposition into orthogonal transverse modes on the spherical boundary is used so that the residual functions may be computed efficiently and concurrently without altering the local character of the finite element equations. Since the auxiliary functions are based on residuals of an asymptotic expansion, the proposed method has the ability to vary separately the radial and transverse modal orders of the RBC. With the number of equations in the auxiliary Cauchy problem equal to the transverse mode number, this reformulation is exact. In this form, the equivalence with the closely related non‐reflecting boundary condition of Grote and Keller is shown. If fewer equations are used, then the boundary conditions form high‐order accurate asymptotic approximations to the exact condition, with corresponding reduction in work and memory. Numerical studies are performed to assess the accuracy and convergence properties of the exact and asymptotic versions of the RBC. The results demonstrate that the asymptotic formulation has dramatically improved accuracy for time domain simulations compared to standard boundary treatments and improved efficiency over the exact condition. Copyright © 2000 John Wiley & Sons, Ltd.     

Two-dimensional Legendre wavelets for solving fractional Poisson equation with Dirichlet boundary conditions

In this paper, the two-dimensional Legendre wavelets are applied for numerical solution of the fractional Poisson equation with Dirichlet boundary conditions. In this way, a new operational matrix of fractional derivative for the Legendre wavelets is derived and then this operational matrix has been employed to obtain the numerical solution of the above-mentioned problem. The main characteristic behind this approach is that it reduces such problems to those of solving a system of algebraic equations which greatly simplifies the problem. The convergence of the two-dimensional Legendre wavelets expansion is investigated. Also the power of this manageable method is illustrated.     

Solving initial and boundary value problems using learning automata particle swarm optimization

In this article, the particle swarm optimization (PSO) algorithm is modified to use the learning automata (LA) technique for solving initial and boundary value problems. A constrained problem is converted into an unconstrained problem using a penalty method to define an appropriate fitness function, which is optimized using the LA-PSO method. This method analyses a large number of candidate solutions of the unconstrained problem with the LA-PSO algorithm to minimize an error measure, which quantifies how well a candidate solution satisfies the governing ordinary differential equations (ODEs) or partial differential equations (PDEs) and the boundary conditions. This approach is very capable of solving linear and nonlinear ODEs, systems of ordinary differential equations, and linear and nonlinear PDEs. The computational efficiency and accuracy of the PSO algorithm combined with the LA technique for solving initial and boundary value problems were improved. Numerical results demonstrate the high accuracy and efficiency of the proposed method.     

A fast boundary cloud method for 3D exterior electrostatic analysis

An accelerated boundary cloud method (BCM) for boundary‐only analysis of 3D electrostatic problems is presented here. BCM uses scattered points unlike the classical boundary element method (BEM) which uses boundary elements to discretize the surface of the conductors. BCM combines the weighted least‐squares approach for the construction of approximation functions with a boundary integral formulation for the governing equations. A linear base interpolating polynomial that can vary from cloud to cloud is employed. The boundary integrals are computed by using a cell structure and different schemes have been used to evaluate the weakly singular and non‐singular integrals. A singular value decomposition (SVD) based acceleration technique is employed to solve the dense linear system of equations arising in BCM. The performance of BCM is compared with BEM for several 3D examples. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical instability in linearized planing problems

The hydrodynamics of planing ships are studied using a finite pressure element method. In this method, a boundary value problem (BVP) is formulated based on linear planing theory; the planing ship is represented by the pressure distribution acting on the wetted bottom of the ship, and the magnitude of this pressure distribution is evaluated using a boundary element method. The pressure is discretized using overlapping pressure pyramids, known as tent functions, so that the resulting distribution is piece‐wise continuous in both longitudinal and transverse directions. A set of linear algebraic equations for the determination of the pressure is then established using a collocation technique. It is found that the matrix of the linear equations is ill conditioned; this leads to oscillatory behaviour of the predicted pressure distribution if the direct solution method of LU decomposition or Gaussian elimination is used to solve the system of linear equations. In the current study, this numerical instability is analysed in detail. It is found that the problem can be addressed by adopting singular value decomposition (SVD) technique for the solution of the linear equations. Using this method, the hydrodynamic results for flat‐bottomed and prismatic planing ships are calculated and a good agreement is demonstrated with Savitsky's empirical relations. Copyright © 2006 John Wiley & Sons, Ltd.     

An extended boundary element method formulation for the direct calculation of the stress intensity factors in fully anisotropic materials

We propose a formulation for linear elastic fracture mechanics in which the stress intensity factors are found directly from the solution vector of an extended boundary element method formulation. The enrichment is embedded in the boundary element method formulation, rather than adding new degrees of freedom for each enriched node. Therefore, a very limited number of new degrees of freedom is added to the problem, which contributes to preserving the conditioning of the linear system of equations. The Stroh formalism is used to provide boundary element method fundamental solutions for any degree of anisotropy, and these are used for both conventional and enriched degrees of freedom. Several numerical examples are shown with benchmark solutions to validate the proposed method. Copyright © 2016 John Wiley & Sons, Ltd.     

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 boundary layer phenomena in two‐dimensional transversely isotropic piezoelectric media by exact symplectic expansion

A separable variable method is introduced to find the exact homogeneous solutions of a two‐dimensional transversely isotropic piezoelectric media to handle general boundary conditions. The usual method of separable variables for partial differentiation equations cannot be readily applicable due to the tangling of the unknowns and their derivatives. Introducing dual variables of stresses, we obtain a set of first‐order Hamiltonian equations whose eigensolutions are symplectic spanning over the solution space to cover all possible boundary conditions. The solutions consist of two parts. The first part is the derogative zero‐eigenvalue solutions of the Saint Venant type together with all their Jordan chains. The second part is the decaying non‐zero‐eigenvalue solutions describing the boundary layer effects. The classical solutions are actually the zero‐eigenvalue solutions representing the simple extension, bending, equipotential field, and the uniform electric displacement. On the other hand, the non‐zero‐eigenvalue solutions represent the localized solutions, which are sensitive to the boundary conditions and are decaying rapidly with respect to the distance from the boundaries. Some rate‐of‐decay curves of the newly found non‐zero‐eigenvalue solutions are shown by numerical examples. Finally, the complete boundary layer effects are quantified for the first time. Copyright © 2006 John Wiley & Sons, Ltd.     

Explicit expressions for 3D boundary integrals in potential theory

On employing isoparametric, piecewise linear shape functions over a flat triangular domain, exact expressions are derived for all surface potentials involved in the numerical solution of three‐dimensional singular and hyper‐singular boundary integral equations of potential theory. These formulae, which are valid for an arbitrary source point in space, are represented as analytic expressions over the edges of the integration triangle. They can be used to solve integral equations defined on polygonal boundaries via the collocation method or may be utilized as analytic expressions for the inner integrals in the Galerkin technique. In addition, the constant element approximation can be directly obtained with no extra effort. Sample problems solved by the collocation boundary element method for the Laplace equation are included to validate the proposed formulae. Published in 2008 by John Wiley & Sons, Ltd.     

Microscale heat conduction in homogeneous anisotropic media: a dual-reciprocity boundary element method and polynomial time interpolation approach

In this paper, a dual-reciprocity boundary element method based on some polynomial interpolations to the time-dependent variables is presented for the numerical solution of a two-dimensional heat conduction problem governed by a third order partial differential equation (PDE) over a homogeneous anisotropic medium. The PDE is derived using a non-Fourier heat flux model which may account for thermal waves and/or microscopic effects. In the analysis, discontinuous linear elements are used to model the boundary and the variables along the boundary. The systems of algebraic equations are set up to solve all the unknowns. For the purpose of evaluating the proposed method, some numerical examples with known exact solutions are solved. The numerical results obtained agree well with the exact solutions.     

A boundary element method solution for anisotropic nonhomogeneous elasticity

A new BEM approach is presented for the plane elastostatic problem for nonhomogeneous anisotropic bodies. In this case the response of the body is described by two coupled linear second order partial differential equations in terms of displacement with variable coefficient. The incapability of establishing the fundamental solution of the governing equations is overcome by uncoupling them using the concept of analog equation, which converts them to two Poisson’s equations, whose fundamental solution is known and the necessary boundary integral equations are readily obtained. This formulation introduces two additional unknown field functions, which physically represent the two components of a fictitious source. Subsequently, they are determined by approximating them globally with radial basis functions series. The displacements and the stresses are evaluated from the integral representation of the solution of the substitutes equations. The presented method maintains the pure boundary character of the BEM. The obtained numerical results demonstrate the effectiveness and accuracy of the method.     

An indirect elastodynamic boundary element method with analytic bases

An indirect time‐domain boundary element method (BEM) is presented here for the treatment of 2D elastodynamic problems. The approximated solution in this method is formulated as a linear combination of a set of particular solutions, which are called bases. The displacement and stress fields of a basis are analytically derived by means of solving Lame's displacement potentials. A semi‐collocation method is proposed to be the time‐stepping algorithm. This method is equivalent to a displacement discontinuity method with piecewise linear discontinuities in both space and time. The resulting time‐stepping scheme is explicit. The BEM is implemented to solve three numerical examples, Lamb's problem, half‐plane with a buried crack and Selberg's problem. Though Lamb's problem is considered a difficult problem for numerical methods, the current numerical results for the surface displacements show accurately the characteristics of the Rayleigh wave. This method is efficient and accurate. Copyright © 2003 John Wiley & Sons, Ltd.     

Finite element analysis of time‐dependent semi‐infinite wave‐guides with high‐order boundary treatment

A new finite element (FE) scheme is proposed for the solution of time‐dependent semi‐infinite wave‐guide problems, in dispersive or non‐dispersive media. The semi‐infinite domain is truncated via an artificial boundary ??, and a high‐order non‐reflecting boundary condition (NRBC), based on the Higdon non‐reflecting operators, is developed and applied on ??. The new NRBC does not involve any high derivatives beyond second order, but its order of accuracy is as high as one desires. It involves some parameters which are chosen automatically as a pre‐process. A C⁰ semi‐discrete FE formulation incorporating this NRBC is constructed for the problem in the finite domain bounded by ??. Augmented and split versions of this FE formulation are proposed. The semi‐discrete system of equations is solved by the Newmark time‐integration scheme. Numerical examples concerning dispersive waves in a semi‐infinite wave guide are used to demonstrate the performance of the new method. Copyright © 2003 John Wiley & Sons, Ltd.     

A meshless solution for two dimensional density-driven groundwater flow

Density-driven groundwater flow is a complicated nonlinear problem in groundwater hydraulics. The local boundary integral method is a promising meshless scheme that is used for solving several difficult problems in different areas. This method applies the boundary integral equations to the local domain around every node. The nodes can be randomly distributed in the domain and on the global boundary. Therefore, this method is characterised as meshless. The unknown potentials and concentrations in all of the nodes are approximated by interpolation to obtain a system of linear equations. Solving this system of equations leads to the numerical solution for the main problem. In this paper, a combination of the radial basis function interpolation and the local boundary element method is used to solve groundwater flow problem combined with the transport of pollution, which also influences the density of groundwater.     

A fictitious domain decomposition method for the solution of partially axisymmetric acoustic scattering problems. Part 2: Neumann boundary conditions

We present a fictitious domain decomposition method for the fast solution of acoustic scattering problems characterized by a partially axisymmetric sound‐hard scatterer. We apply this method to the solution of a mock‐up submarine problem, and highlight its computational advantages and intrinsic parallelism. A key component of our method is an original idea for addressing a Neumann boundary condition in the general framework of a fictitious domain method. This idea is applicable to many other linear partial differential equations besides the Helmholtz equation. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical transfer-method with boundary conditions for arbitrary curved beam elements

A novel numerical transfer-method is presented to solve a system of linear ordinary differential equations with boundary conditions. It is applied to determine the structural behaviour of the classical problem of an arbitrary curved beam element. The approach of this boundary value problem yields a unique system of differential equations. A Runge–Kutta scheme is chosen to obtain the incremental transfer expression. The use of a recurrence strategy in this equation permits to relate both ends in the domain where boundary conditions are defined. Semicircular arch, semicircular balcony and elliptic–helical beam examples are provided for validation.