A fast and accurate semi-Lagrangian particle level set method

Abstract

In this paper, we present an efficient semi-Lagrangian based particle level set method for the accurate capturing of interfaces. This method retains the robust topological properties of the level set method with- out the adverse effects of numerical dissipation. Both the level set method and the particle level set method typically use high order accurate numerical discretizations in time and space, e.g. TVD Runge–Kutta and HJ-WENO schemes. We demonstrate that these computationally expensive schemes are not required. Instead, fast, low order accurate numerical schemes suffice. That is, the addition of particles to the level set method not only removes the difficulties associated with numerical diffusion, but also alleviates the need for computationally expensive high order accurate schemes. We use an efficient, first order accurate semi-Lagrangian advection scheme coupled with a first order accurate fast marching method to evolve the level set function. To accurately track the underlying flow characteristics, the particles are evolved with a second order accurate method. Since we avoid complex high order accurate numerical methods, extending the algorithm to arbitrary data structures becomes more feasible, and we show preliminary results obtained with an octree-based adaptive mesh.