A Convergent Iterated Quasi-interpolation for Periodic Domain and Its Applications to Surface PDEs

The paper first provides an iterated quasi-interpolation scheme for function approximation over periodic domain and then attempts its applications to solve time-dependent surface partial differential equations (PDEs). The key feature of our scheme is that it gives an approximation directly just by taking a weighted average of the available discrete function values evaluated at sampling centers in the periodic domain. As such, it is simple, easy to compute, and the implementation process is stable. Moreover, if the sampling centers distribute uniformly over the periodic domain, it even preserves the same convergence order to all the derivatives. To employ the iterated quasi-interpolation scheme for solving time-dependent surface PDEs, we follow the framework of the method-of-lines. More precisely, we first reformulate the PDEs in terms of the parametric form of the surface. Then we use our quasi-interpolation scheme to approximate both the analytic solution and its spatial derivatives in the reformulated form (of PDEs) to get a semi-discrete ordinary differential equation (ODE) system. Finally, we adopt an appropriate time-integration technique to obtain the full-discrete scheme. As concrete examples, we take the torus for illustration and solve some benchmark reaction-diffusion equations imposed on the torus at the end of the paper. However, the proposed method is general and works on time-dependent PDEs defined on any smooth closed parameterized surfaces without coordinate singularity.