Solving Partial Differential Equations on Closed Surfaces with Planar Cartesian Grids
share
We present a general purpose method for solving partial differential equations on a closed surface, using a technique for discretizing the surface introduced by Wenjun Ying and Wei-Cheng Wang (J. Comput. Phys., 2013), which uses projections on coordinate planes. The surface is represented by a set of points at which it intersects the intervals between grid points in a three-dimensional grid. They are designated as primary or secondary. Discrete functions on the surface have independent values at primary points, with values at secondary points determined by an equilibration process. Each primary point and its neighbors have projections to regular grid points in a coordinate plane where the equilibration is done and finite differences are computed. The solution of a p.d.e. can be reduced to standard methods on Cartesian grids in the coordinate planes, with the equilibration allowing seamless transition from one system to another. We obtain second order accuracy for solutions of a variety of equations, including surface diffusion determined by the Laplace-Beltrami operator and the shallow water equations on a sphere.
READ FULL TEXT