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

• Published in: Journal of scientific computing Vol. 93; no. 2; p. 37
• Main Authors: Sun, Zhengjie, Gao, Wenwu, Yang, Ran
• Format: Journal Article
• Language: English
• Published: New York Springer US 01.11.2022; Springer Nature B.V
• Subjects: Algorithms; Approximation; Computational Mathematics and Numerical Analysis; Convergence; Derivatives; Discrete functions; Exact solutions; Fourier transforms; Interpolation; Mathematical analysis; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematics; Mathematics and Statistics; Method of lines; Methods; Numerical analysis; Partial differential equations; Reaction-diffusion equations; Sampling; Singularity (mathematics); Theoretical; Time dependence; Time integration; Toruses; Quasi-interpolation; Derivative approximation; 41A25; 41A30; 41A63; Torus; Multiquadric trigonometric kernel; Surface PDE; 41A05
• ISSN: 0885-7474; 1573-7691
• DOI: 10.1007/s10915-022-01998-2

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.