Closest point method

From HandWiki

The closest point method (CPM) is an embedding method for solving partial differential equations on surfaces. The closest point method uses standard numerical approaches such as finite differences, finite element or spectral methods in order to solve the embedding partial differential equation (PDE) which is equal to the original PDE on the surface. The solution is computed in a band surrounding the surface in order to be computationally efficient. In order to extend the data off the surface, the closest point method uses a closest point representation. This representation extends function values to be constant along directions normal to the surface.

Definitions

Closest Point function: Given a surface [math]\displaystyle{ \mathcal S, cp(\mathbf x) }[/math] refers to a (possibly non-unique) point belonging to [math]\displaystyle{ \mathcal S }[/math], which is closest to [math]\displaystyle{ \mathbf x }[/math] [SE].

Closest point extension: Let [math]\displaystyle{ \mathcal S }[/math], be a smooth surface in [math]\displaystyle{ \mathbb R^d }[/math]. The closest point extension of a function [math]\displaystyle{ u : \mathcal S \rightarrow \mathbb R }[/math], to a neighborhood [math]\displaystyle{ \Omega }[/math] of [math]\displaystyle{ \mathcal S }[/math], is the function [math]\displaystyle{ v: \Omega \rightarrow \mathbb R }[/math], defined by [math]\displaystyle{ v(\mathbf x) = u(cp(\mathbf x)) }[/math].

Closest point method

Initialization consists of these steps [EW]:

  • If it is not already given, a closest point representation of the surface is constructed.
  • A computational domain is chosen. Typically this is a band around the surface.
  • Replace surface gradients by standard gradients in [math]\displaystyle{ \mathbb R^3 }[/math].
  • Solution is initialized by extending the initial surface data on to the computational domain using the closest point function.

After initialization, alternate between the following two steps:

  • Using the closest point function, extend the solution off the surface to the computational domain.
  • Compute the solution to the embedding PDE on a Cartesian mesh in the computational domain for one time step.

Banding

The surface PDE is extended into [math]\displaystyle{ \mathbb R^3 }[/math] however it is only necessary to solve this new PDE near the surface. Hence, we solve the PDE in a band surrounding the surface for efficient computational purposes. [math]\displaystyle{ \Omega_c {x : \| x - cp(x) \|_2 \leq \lambda} }[/math] where [math]\displaystyle{ \lambda }[/math] is the bandwidth.

Example: Heat equation on a circle

Using initial profile [math]\displaystyle{ u_S (\theta , t) = \sin (\theta) }[/math] leads to the solution [math]\displaystyle{ u_S (\theta, t) = \exp (-t)\sin (\theta) }[/math] for the heat equation. Forward Euler time-stepping is used with relation [math]\displaystyle{ \Delta t = 0.1 \Delta x^2 }[/math] and degree-four interpolation polynomials for the interpolations. Second-order centered differences are used for the spatial discretization. The CPM results in the expected second order error in the solution [math]\displaystyle{ u }[/math].

Applications

The closest point method can be applied to various PDEs on surfaces. Reaction–diffusion problems on point clouds [RD], eigenvalue problems [EV], and level set equations [LS] are a few examples.

See also

References

  • [EM] Ruuth, S. J., & Merriman, B. (2008). A simple embedding method for solving partial differential equations on surfaces.Journal of Computational Physics,227(3), 1943–1961 here
  • [RD] Macdonald, C. B., Merriman, B., & Ruuth, S. J. (2013). Simple computation of reaction-diffusion processes on point clouds. Proceedings of the National Academy of Sciences, 110(23), 9209–9214 here
  • [EV] Macdonald, C. B., Brandman, J., & Ruuth, S. J. (2011). Solving eigenvalue problems on curved surfaces using the Closest Point Method. Journal of Computational Physics, 230(22), 7944–7956. here
  • [LS] Macdonald, C. B., & Ruuth, S. J. (2008). Level set equations on surfaces via the Closest Point Method. Journal of Scientific Computing, 35(2–3), 219–240. here