PDE-BASED MULTIDIMENSIONAL EXTRAPOLATION OF SCALAR FIELDS OVER INTERFACES WITH KINKS AND HIGH CURVATURES

摘要

We present a PDE-based approach for the multidimensional extrapolation of smooth scalar quantities across interfaces with kinks and regions of high curvature. Unlike the commonly used method of [T. Aslam, T. Comput. Phys., 193 (2004), pp. 349-355], in which normal derivatives are extrapolated, the proposed approach is based on the extrapolation and weighting of Cartesian derivatives. As a result, second- and third-order accurate extensions in the L-infinity norm are obtained with linear and quadratic extrapolations, respectively, even in the presence of sharp geometric features. The accuracy of the method is demonstrated on a number of examples in two and three spatial dimensions and compared to the approach of [T. Aslam, T. Comput. Phys., 193 (2004), pp. 349-355]. The importance of accurate extrapolation near sharp geometric features is highlighted on an example of solving the diffusion equation on evolving domains.