Curl-Flow: Boundary-Respecting Pointwise Incompressible Velocity Interpolation for Grid-Based Fluids

Abstract

We propose to augment standard grid-based fluid solvers with pointwise divergence-free velocity interpolation, thereby ensuring exact incompressibility down to the sub-cell level. Our method takes as input a discretely divergence-free velocity field generated by a staggered grid pressure projection, and first recovers a corresponding discrete vector potential. Instead of solving a costly vector Poisson problem for the potential, we develop a fast parallel sweeping strategy to find a candidate potential and apply a gauge transformation to enforce the Coulomb gauge condition and thereby make it numerically smooth. Interpolating this discrete potential generates a pointwise vector potential whose analytical curl is a pointwise incompressible velocity field. Our method further supports irregular solid geometry through the use of level set-based cut-cells and a novel Curl-Noise-inspired potential ramping procedure that simultaneously offers strictly non-penetrating velocities and incompressibility. Experimental comparisons demonstrate that the vector potential reconstruction procedure at the heart of our approach is consistently faster than prior such reconstruction schemes, especially those that solve vector Poisson problems. Moreover, in exchange for its modest extra cost, our overall Curl-Flow framework produces significantly improved particle trajectories that closely respect irregular obstacles, do not suffer from spurious sources or sinks, and yield superior particle distributions over time.

References in corpus (8)

• Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
• Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
• High order exactly divergence-free Hybrid Discontinuous Galerkin Methods for unsteady incompressible flows
• Divergence-Free Reconstruction of Magnetic Fields and WENO Schemes for Magnetohydrodynamics
• A hybridizable discontinuous Galerkin method for the Navier--Stokes equations with pointwise divergence-free velocity field
• An Immersed Boundary Method with Divergence-Free Velocity Interpolation and Force Spreading
• Divergence-free Interpolation of Vector Fields From Point Values - Exact divB=0 in Numerical Simulations
• Numerical generation of vector potentials from specified magnetic fields

Cited by in corpus (2)

• Fluid Implicit Particles on Coadjoint Orbits
• Cell-Constrained Particles for Incompressible Fluids