2025/12/06 by Blomquist, Matthew, Gaudreault, Stéphane, Theillard, Maxime
Engineering · Physics and Astronomy · #Advanced Numerical Methods in Computational Mathematics #Computational Fluid Dynamics and Aerodynamics #FOS: Mathematics #FOS: Physical sciences #Fluid Dynamics (physics.flu-dyn) #Model Reduction and Neural Networks #Numerical Analysis (math.NA)
paper · doi:10.48550/arxiv.2512.06455
openalex publication_date 2025/12/06 · openalex created_date 2025/12/10 · openalex updated_date 2026/07/28
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.