Characteristic Bending in Incompressible Flows

TL;DR

The paper introduces the Characteristic Bending (CB) method, a novel approach that enhances the accuracy and robustness of advecting quantities in incompressible flows by systematically removing compressible errors through a volume-preserving projection.

Contribution

The CB method provides a new framework for advecting quantities in incompressible flows by interpreting characteristic reconstruction as a reference map and applying a volume-preserving projection to improve accuracy.

Findings

Improves mass and geometric feature preservation in simulations.

Serves as a drop-in replacement for traditional semi-Lagrangian schemes.

Demonstrates robustness in 2D and 3D incompressible flow benchmarks.

Abstract

We present the Characteristic Bending (CB) method, a general framework for advecting quantities under incompressible velocity fields. The method builds on standard semi-Lagrangian advection by interpreting the backward-in-time characteristic reconstruction as the construction of a reference map, a diffeomorphism between the current and initial geometries of the advected space. From this viewpoint, the CB method applies a volume-preserving projection to the map, systematically removing spurious compressible errors arising from time integration, interpolation, or from velocity fields that are only approximately divergence-free. This projection bends the characteristics toward the divergence-free space, preserving mass and geometric features of the advected fields, even in the presence of significant error. We demonstrate the method in both two and three dimensions using benchmark problems and for multiphase flows governed by the incompressible Navier-Stokes equations. The results show that the CB method serves as a drop-in replacement for traditional semi-Lagrangian schemes and as an augmentation of reference map formulations, offering improved robustness and accuracy in incompressible flow simulations.