Boundary particle method

Last updated

In applied mathematics, the boundary particle method (BPM) is a boundary-only meshless (meshfree) collocation technique, in the sense that none of inner nodes are required in the numerical solution of nonhomogeneous partial differential equations. Numerical experiments show that the BPM has spectral convergence. Its interpolation matrix can be symmetric.

Contents

History and recent developments

In recent decades, the dual reciprocity method (DRM) [1] and multiple reciprocity method (MRM) [2] have been emerging as promising techniques to evaluate the particular solution of nonhomogeneous partial differential equations in conjunction with the boundary discretization techniques, such as boundary element method (BEM). For instance, the so-called DR-BEM and MR-BEM are popular BEM techniques in the numerical solution of nonhomogeneous problems.

The DRM has become a common method to evaluate the particular solution. However, the DRM requires inner nodes to guarantee the convergence and stability. The MRM has an advantage over the DRM in that it does not require using inner nodes for nonhomogeneous problems.[ citation needed ] Compared with the DRM, the MRM is computationally more expensive in the construction of the interpolation matrices and has limited applicability to general nonhomogeneous problems due to its conventional use of high-order Laplacian operators in the annihilation process.

The recursive composite multiple reciprocity method (RC-MRM), [3] [4] was proposed to overcome the above-mentioned problems. The key idea of the RC-MRM is to employ high-order composite differential operators instead of high-order Laplacian operators to eliminate a number of nonhomogeneous terms in the governing equation. The RC-MRM uses the recursive structures of the MRM interpolation matrix to reduce computational costs.

The boundary particle method (BPM) is a boundary-only discretization of an inhomogeneous partial differential equation by combining the RC-MRM with strong-form meshless boundary collocation discretization schemes, such as the method of fundamental solution (MFS), boundary knot method (BKM), regularized meshless method (RMM), singular boundary method (SBM), and Trefftz method (TM). The BPM has been applied to problems such as nonhomogeneous Helmholtz equation and convection–diffusion equation. The BPM interpolation representation is of a wavelet series.

For the application of the BPM to Helmholtz, [3] Poisson [4] and plate bending problems, [5] the high-order fundamental solution or general solution, harmonic function [6] or Trefftz function (T-complete functions) [7] are often used, for instance, those of Berger, Winkler, and vibrational thin plate equations. [8] The method has been applied to inverse Cauchy problem associated with Poisson [9] and nonhomogeneous Helmholtz equations. [10]

Further comments

The BPM may encounter difficulty in the solution of problems having complex source functions, such as non-smooth, large-gradient functions, or a set of discrete measured data. The solution of such problems involves:[ citation needed ]

(1) The complex functions or a set of discrete measured data can be interpolated by a sum of polynomial or trigonometric function series. Then, the RC-MRM can reduce the nonhomogeneous equation to a high-order homogeneous equation, and the BPM can be implemented to solve these problems with boundary-only discretization.

(2) The domain decomposition may be used to in the BPM boundary-only solution of large-gradient source functions problems.

See also

Related Research Articles

Partial differential equation differential equation that contains unknown multivariable functions and their partial derivatives

In mathematics, a partial differential equation (PDE) is a differential equation that contains unknown multivariable functions and their partial derivatives. PDEs are used to formulate problems involving functions of several variables, and are either solved by hand, or used to create a computer model. A special case is ordinary differential equations (ODEs), which deal with functions of a single variable and their derivatives.

Computational fluid dynamics branch of fluid mechanics that uses numerical analysis and data structures to solve and analyze problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

Numerical methods for partial differential equations is the branch of numerical analysis that studies the numerical solution of partial differential equations (PDEs).

The boundary element method (BEM) is a numerical computational method of solving linear partial differential equations which have been formulated as integral equations. including fluid mechanics, acoustics, electromagnetics, fracture mechanics, and contact mechanics.

Computational electromagnetics (CEM), computational electrodynamics or electromagnetic modeling is the process of modeling the interaction of electromagnetic fields with physical objects and the environment.

Meshfree methods Methods in numerical analysis not requiring knowledge of neighboring piints

In the field of numerical analysis, meshfree methods are those that do not require connection between nodes of the simulation domain, i.e. a mesh, but are rather based on interaction of each node with all its neighbors. As a consequence, original extensive properties such as mass or kinetic energy are no longer assigned to mesh elements but rather to the single nodes. Meshfree methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort. The absence of a mesh allows Lagrangian simulations, in which the nodes can move according to the velocity field.

The beam propagation method (BPM) is an approximation technique for simulating the propagation of light in slowly varying optical waveguides. It is essentially the same as the so-called parabolic equation (PE) method in underwater acoustics. Both BPM and the PE were first introduced in the 1970s. When a wave propagates along a waveguide for a large distance, rigorous numerical simulation is difficult. The BPM relies on approximate differential equations which are also called the one-way models. These one-way models involve only a first order derivative in the variable z and they can be solved as "initial" value problem. The "initial" value problem does not involve time, rather it is for the spatial variable z.

In mathematics, the Trefftz method is a method for the numerical solution of partial differential equations named after the German mathematician Erich Trefftz(de) (1888–1937). It falls within the class of finite element methods.

Light scattering by particles is the process by which small particles cause optical phenomena, such as rainbows, the blue color of the sky, and halos.

The material point method (MPM) is a numerical technique used to simulate the behavior of solids, liquids, gases, and any other continuum material. Especially, it is a robust spatial discretization method for simulating multi-phase (solid-fluid-gas) interactions. In the MPM, a continuum body is described by a number of small Lagrangian elements referred to as 'material points'. These material points are surrounded by a background mesh/grid that is used only to calculate gradient terms such as the deformation gradient. Unlike other mesh-based methods like the finite element method, finite volume method or finite difference method, the MPM is not a mesh based method and is instead categorized as a meshless/meshfree or continuum-based particle method, examples of which are smoothed particle hydrodynamics and peridynamics. Despite the presence of a background mesh, the MPM does not encounter the drawbacks of mesh-based methods which makes it a promising and powerful tool in computational mechanics.

Singular boundary method

In numerical analysis, the singular boundary method (SBM) belongs to a family of meshless boundary collocation techniques which include the method of fundamental solutions (MFS), boundary knot method (BKM), regularized meshless method (RMM), boundary particle method (BPM), modified MFS, and so on. This family of strong-form collocation methods is designed to avoid singular numerical integration and mesh generation in the traditional boundary element method (BEM) in the numerical solution of boundary value problems with boundary nodes, in which a fundamental solution of the governing equation is explicitly known.

In scientific computation and simulation, the method of fundamental solutions (MFS) is getting a growing attention. The method is essentially based on the fundamental solution of a partial differential equation of interest as its basis function. The MFS was developed to overcome the major drawbacks in the boundary element method (BEM) which also uses the fundamental solution to satisfy the governing equation. Consequently, both the MFS and the BEM are of a boundary discretization numerical technique and reduce the computational complexity by one dimensionality and have particular edge over the domain-type numerical techniques such as the finite element and finite volume methods on the solution of infinite domain, thin-walled structures, and inverse problems.

In numerical mathematics, the boundary knot method (BKM) is proposed as an alternative boundary-type meshfree distance function collocation scheme.

In numerical mathematics, the regularized meshless method (RMM), also known as the singular meshless method or desingularized meshless method, is a meshless boundary collocation method designed to solve certain partial differential equations whose fundamental solution is explicitly known. The RMM is a strong-form collocation method with merits being meshless, integration-free, easy-to-implement, and high stability. Until now this method has been successfully applied to some typical problems, such as potential, acoustics, water wave, and inverse problems of bounded and unbounded domains.

The Kansa method is a computer method used to solve partial differential equations. Partial differential equations are mathematical models of things like stresses in a car's body, air flow around a wing, the shock wave in front of a supersonic airplane, quantum mechanical model of an atom, ocean waves, socio-economic models, digital image processing etc. The computer takes the known quantities such as pressure, temperature, air velocity, stress, and then uses the laws of physics to figure out what the rest of the quantities should be like a puzzle being fit together. Then, for example, the stresses in various parts of a car can be determined when that car hits a bump at 70 miles per hour.

The discrete least squares meshless (DLSM) method is a meshless method based on the least squares concept. The method is based on the minimization of a least squares functional, defined as the weighted summation of the squared residual of the governing differential equation and its boundary conditions at nodal points used to discretize the domain and its boundaries. While most of the existing meshless methods need background cells for numerical integration, DLSM did not require a numerical integration procedure due to the use of the discrete least squares method to discretize the governing differential equation. A Moving least squares (MLS) approximation method is used to construct the shape function, making the approach a fully least squares-based approach.

The finite point method (FPM) is a meshfree method for solving partial differential equations (PDEs) on scattered distributions of points. The FPM was proposed in the mid-nineties in, and with the purpose to facilitate the solution of problems involving complex geometries, free surfaces, moving boundaries and adaptive refinement. Since then, the FPM has evolved considerably, showing satisfactory accuracy and capabilities to deal with different fluid and solid mechanics problems.

The generalized-strain mesh-free (GSMF) formulation is a local meshfree method in the field of numerical analysis, completely integration free, working as a weighted-residual weak-form collocation. This method was first presented by Oliveira and Portela (2016), in order to further improve the computational efficiency of meshfree methods in numerical analysis. Local meshfree methods are derived through a weighted-residual formulation which leads to a local weak form that is the well known work theorem of the theory of structures. In an arbitrary local region, the work theorem establishes an energy relationship between a statically-admissible stress field and an independent kinematically-admissible strain field. Based on the independence of these two fields, this formulation results in a local form of the work theorem that is reduced to regular boundary terms only, integration-free and free of volumetric locking.

In the field of numerical analysis, particle methods discretize fluid into particles. Particle methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort.

References

  1. Partridge PW, Brebbia CA, Wrobel LC, The dual reciprocity boundary element method. Computational Mechanics Publications, 1992
  2. Nowak AJ, Neves AC, The multiple reciprocity boundary element method. Computational Mechanics Publication, 1994
  3. 1 2 Chen W, "Meshfree boundary particle method applied to Helmholtz problems". Engineering Analysis with Boundary Elements 2002,26(7): 577–581
  4. 1 2 Chen W, Fu ZJ, Jin BT, "A truly boundary-only meshfree method for inhomogeneous problems based on recursive composite multiple reciprocity technique". Engineering Analysis with Boundary Elements 2010,34(3): 196–205
  5. Fu ZJ, Chen W, Yang W, Winkler plate bending problems by a truly boundary-only boundary particle method. Computational Mechanics 2009,44(6): 757–563
  6. Hon YC, Wu ZM, "A numerical computation for inverse boundary determination problem" Engineering Analysis with Boundary Elements 2000,24(7–8): 599–606
  7. Chen W, Fu ZJ, Qin QH, "Boundary particle method with high-order Trefftz functions". CMC: Computers, Materials & Continua 2010,13(3): 201–217
  8. Chen W, Shen ZJ, Shen LJ, Yuan GW, "General solutions and fundamental solutions of varied orders to the vibrational thin, the Berger, and the Winkler plates" Engineering Analysis with Boundary Elements 2005,29(7): 699–702
  9. Fu ZJ, Chen W, Zhang CZ, "Boundary particle method for Cauchy inhomogeneous potential problems". Inverse Problems in Science and Engineering 2012,20(2): 189–207
  10. Chen W, Fu ZJ, "Boundary particle method for inverse Cauchy problems of inhomogeneous Helmholtz equations". Journal of Marine Science and Technology–Taiwan 2009,17(3): 157–163

Free software and Matlab codes