DDFKs: Fluid Simulation with Dynamic Divergence-Free Kernels

TL;DR

This paper introduces DDFKs, a grid-free, divergence-free kernel-based solver for incompressible fluid simulation that maintains incompressibility, reduces numerical dissipation, and efficiently models complex vortex-rich flows.

Contribution

The paper proposes a novel grid-free solver using divergence-free kernels that effectively enforce incompressibility and preserve flow details in fluid simulations.

Findings

Achieves comparable accuracy and robustness to state-of-the-art methods

Maintains incompressibility and preserves vortices effectively

Exhibits lower numerical dissipation and high efficiency

Abstract

Fluid simulations based on memory-efficient spatial representations like implicit neural spatial representations (INSRs) and Gaussian spatial representation (GSR), where the velocity fields are parameterized by neural networks or weighted Gaussian functions, has been an emerging research area. Though advantages over traditional discretizations like spatial adaptivity and continuous differentiability of these spatial representations are leveraged by fluid solvers, solving the time-dependent PDEs that governs the fluid dynamics remain challenging, especially in incompressible fluids where the divergence-free constraint is enforced. In this paper, we propose a grid-free solver Dynamic Divergence-Free Kernels (DDFKs) for incompressible flows based on divergence-free kernels (DFKs). Each DFK is incorporated with a matrix-valued radial basis function and a vector-valued weight, yielding a divergence-free vector field. We model the continuous flow velocity as the sum of multiple DFKs, thus enforcing incompressibility while being able to preserve different level of details. Quantitative and qualitative results show that our method achieves comparable accuracy, robustness, ability to preserve vortices, time and memory efficiency and generality across diverse phenomena to state-of-the-art methods using memory-efficient spatial representations, while excels at maintaining incompressibility. Though our first-order solver are slower than fluid solvers with traditional discretizations, our approach exhibits significantly lower numerical dissipation due to reduced discretization error. We demonstrate our method on diverse incompressible flow examples with rich vortices and various solid boundary conditions.