A least-squares implicit RBF-FD closest point method and applications to PDEs on moving surfaces

The closest point method (Ruuth and Merriman, J. Comput. Phys. 227(3):1943-1961, [2008]) is an embedding method developed to solve a variety of partial differential equations (PDEs) on smooth surfaces, using a closest point representation of the surface and standard Cartesian grid methods in the embedding space. Recently, a closest point method with explicit time-stepping was proposed that uses finite differences derived from radial basis functions (RBF-FD). Here, we propose a least-squares implicit formulation of the closest point method to impose the constant-along-normal extension of the solution on the surface into the embedding space. Our proposed method is particularly flexible with respect to the choice of the computational grid in the embedding space. In particular, we may compute over a computational tube that contains problematic nodes. This fact enables us to combine the proposed method with the grid based particle method (Leung and Zhao, J. Comput. Phys. 228(8):2993-3024, [2009]) to obtain a numerical method for approximating PDEs on moving surfaces. We present a number of examples to illustrate the numerical convergence properties of our proposed method. Experiments for advection-diffusion equations and Cahn-Hilliard equations that are strongly coupled to the velocity of the surface are also presented

Elastic flow interacting with a lateral diffusion process : the one-dimensional graph case

A finite element approach to the elastic flow of a curve coupled with a diffusion equation on the curve is analysed. Considering the graph case, the problem is weakly formulated and approximated with continuous linear finite elements, which is enabled thanks to second-order operator splitting. The error analysis builds up on previous results for the elastic flow. To obtain an error estimate for the quantity on the curve a better control of the velocity is required. For this purpose, a penalty approach is employed and then combined with a generalized Gronwall lemma. Numerical simulations support the theoretical convergence results. Further numerical experiments indicate stability beyond the parameter regime with respect to the penalty term that is covered by the theory

Trace Finite Element Methods for PDEs on Surfaces

In this paper we consider a class of unfitted finite element methods for discretization of partial differential equations on surfaces. In this class of methods known as the Trace Finite Element Method (TraceFEM), restrictions or traces of background surface-independent finite element functions are used to approximate the solution of a PDE on a surface. We treat equations on steady and time-dependent (evolving) surfaces. Higher order TraceFEM is explained in detail. We review the error analysis and algebraic properties of the method. The paper navigates through the known variants of the TraceFEM and the literature on the subject