An adaptive semi-Lagrangian advection scheme and its parallelization

Abstract

A semi-Lagrangian advection scheme has been implemented for the linear passive advection equation. The advected scalar field is represented in terms of finite-element bases on a self-adaptive grid. Results concerning the accuracy of the numerical scheme compare well to those with finite-difference schemes on regular meshes found in the literature. Quasi-monotone and quasi-conservative interpolation schemes have been implemented. A parallelized version of the code shows good performance on a KSR-1 virtual shared memory computer.