Journal
IMA JOURNAL OF NUMERICAL ANALYSIS
Volume 27, Issue 2, Pages 262-292Publisher
OXFORD UNIV PRESS
DOI: 10.1093/imanum/drl023
Keywords
finite elements; evolving surfaces; conservation; diffusion; existence; error estimates; computations
Categories
Funding
- Engineering and Physical Sciences Research Council [GR/S70883/01] Funding Source: researchfish
Ask authors/readers for more resources
In this article, we define a new evolving surface finite-element method for numerically approximating partial differential equations on hypersurfaces Gamma(t) in < Ropf >(n+1) which evolve with time. The key idea is based on approximating Gamma(t) by an evolving interpolated polyhedral (polygonal if n = 1) surface Gamma(h)(t) consisting of a union of simplices (triangles for n = 2) whose vertices lie on Gamma(t). A finite-element space of functions is then defined by taking the set of all continuous functions on Gamma(h)(t) which are linear affine on each simplex. The finite-element nodal basis functions enjoy a transport property which simplifies the computation. We formulate a conservation law for a scalar quantity on Gamma(t) and, in the case of a diffusive flux, derive a transport and diffusion equation which takes into account the tangential velocity and the local stretching of the surface. Using surface gradients to define weak forms of elliptic operators naturally generates weak formulations of elliptic and parabolic equations on Gamma(t). Our finite-element method is applied to the weak form of the conservation equation. The computations of the mass and element stiffness matrices are simple and straightforward. Error bounds are derived in the case of semi-discretization in space. Numerical experiments are described which indicate the order of convergence and also the power of the method. We describe how this framework may be employed in applications.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available