Lagrangian Covector Fluid with Free Surface
Zhiqi Li, Barnabás Börcsök, Duowen Chen, Yutong Sun, Bo Zhu, Greg Turk
TL;DR
The paper presents a novel Lagrangian covector fluid framework that uses flow maps and path integrals to decouple the advection and projection steps for incompressible free surface flows. By introducing long range and short range mapping strategies and a unifying LMSP scheme, the authors transform challenging integral boundary Poisson problems into standard Dirichlet/Neumann solves, improving robustness and vorticity preservation. A Voronoi based meshfree implementation enables accurate gradient/divergence operators on moving particles, and results across 2D and 3D benchmarks demonstrate superior performance over purely particle based methods like the Power Particle Method. The approach opens new avenues for handling free surfaces in covector flow maps and offers a practical route toward robust, vortex preserving fluid simulations with complex boundaries.
Abstract
This paper introduces a novel Lagrangian fluid solver based on covector flow maps. We aim to address the challenges of establishing a robust flow-map solver for incompressible fluids under complex boundary conditions. Our key idea is to use particle trajectories to establish precise flow maps and tailor path integrals of physical quantities along these trajectories to reformulate the Poisson problem during the projection step. We devise a decoupling mechanism based on path-integral identities from flow-map theory. This mechanism integrates long-range flow maps for the main fluid body into a short-range projection framework, ensuring a robust treatment of free boundaries. We show that our method can effectively transform a long-range projection problem with integral boundaries into a Poisson problem with standard boundary conditions -- specifically, zero Dirichlet on the free surface and zero Neumann on solid boundaries. This transformation significantly enhances robustness and accuracy, extending the applicability of flow-map methods to complex free-surface problems.
