Lorentz force velocimetry

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 contrast to the BEM, the MFS avoids the numerical integration of singular fundamental solution and is an inherent meshfree method. The method, however, is compromised by requiring a controversial fictitious boundary outside the physical domain to circumvent the singularity of fundamental solution, which has seriously restricted its applicability to real-world problems. But nevertheless the MFS has been found very competitive to some application areas such as infinite domain problems.

The MFS is also known by quite a few different names in the literature. Among these are the charge simulation method, the superposition method, the desingularized method, the indirect boundary element method, and the virtual boundary element method, just to name a few.

MFS formulation

Consider a partial differential equation governing certain type of problems

$Lu=f\left(x,y\right),\ \ \left(x,y\right)\in \Omega ,$

$u=g\left(x,y\right),\ \ \left(x,y\right)\in \partial \Omega _{D},$

${\frac {\partial u}{\partial n}}=h\left(x,y\right),\ \ h\left(x,y\right)\in \partial \Omega _{N},$

where $L$ is the differential partial operator, $\Omega$ represents the computational domain, $\partial \Omega _{D}$ and $\partial \Omega _{N}$ denote the Dirichlet and Neumann boundary, respectively, $\partial \Omega _{D}\cup \partial \Omega _{N}=\partial \Omega$ and $\partial \Omega _{D}\cap \partial \Omega _{N}=\varnothing$.

The MFS employs the fundamental solution of the operator as its basis function to represent the approximation of unknown function u as follows

${{u}^{*}}\left(x,y\right)=\sum \limits _{i=1}^{N}\alpha _{i}\phi \left(r_{i}\right)$

where $r_{i}=\left\|\left(x,y\right)-\left(sx_{i},sy_{i}\right)\right\|$ denotes the Euclidean distance between collocation points $\left(x,y\right)$ and source points $\left(sx_{i},sy_{i}\right)$, $\phi \left(\cdot \right)$ is the fundamental solution which satisfies

$L\phi =\delta \,$

where $\delta$ denotes Dirac delta function, and ${{\alpha }_{i}}$ are the unknown coefficients.

With the source points located outside the physical domain, the MFS avoid the fundamental solution singularity. Substituting the approximation into boundary condition yields the following matrix equation

$\left[{\begin{matrix}\phi \left(\left.r_{j}\right|_{x_{i},y_{i}}\right)\\{\frac {\partial \phi \left(\left.r_{j}\right|_{x_{k},y_{k}}\right)}{\partial n}}\\\end{matrix}}\right]\ \cdot \ \alpha =\left({\begin{matrix}g\left(x_{i},y_{i}\right)\\h\left(x_{k},y_{k}\right)\\\end{matrix}}\right),$

where $\left(x_{i},y_{i}\right)$ and $\left(x_{k},y_{k}\right)$ denote the collocation points, respectively, on Dirichlet and Neumann boundaries. The unknown coefficients $\alpha _{i}$ can uniquely be determined by the above algebraic equation. And then we can evaluate numerical solution at any location in physical domain.

History and recent developments

The ideas behind the MFS have been around for a few decades and were developed primarily by V. D. Kupradze and M. A. Alexidze in the late 1950s and early 1960s. However, the method was proposed as a computational technique much later by R. Mathon and R. L. Johnston in the late 1970s, followed by a number of papers by Mathon, Johnston and Graeme Fairweather with applications. Slowly but surely the MFS becomes a useful tool for the solution of a large variety of physical and engineering problems.

A major obstacle was overcome when, in the 1990s, M. A. Golberg and C. S. Chen extended the MFS to deal with inhomogeneous equations and time-dependent problems. Recent developments indicate that the MFS can be used to solve partial differential equations with variable coefficients. The MFS has proved particularly effective for certain classes of problems such as inverse, unbounded domain, and free-boundary problems.

Some new techniques have recently been developed to cure the fictitious boundary problem in the MFS, such as the boundary knot method, singular boundary method, and regularized meshless method.

See also

• Radial basis function
• Boundary element method
• Boundary knot method
• Boundary particle method
• Singular boundary method

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

External links

• International Center for Numerical Simulation Software in Engineering & Sciences