Adaptive Fluid Cohomology on Surfaces

Abstract

Simulating inviscid, incompressible fluids on non-simply-connected curved surfaces requires careful treatment of the flow's local and global behavior. While recent theoretical advancements have established the critical dynamics of the harmonic component in such flows, practical applications remain computationally restricted by a lack of spatial and temporal adaptivity. Furthermore, simulations on poor-quality meshes often lead to numerical instability and a failure to preserve the flow's underlying harmonic component when using naive interpolation methods. In this paper, we introduce Adaptive Fluid Cohomology, a framework that integrates dynamic spatial and temporal refinement into the simulation of the Euler equations. We leverage a posteriori error estimation to adjust spatial resolution on the fly, alongside a standard Dormand-Prince 5(4) time-stepping scheme for temporal accuracy. To ensure stability during mesh mutations, we develop a novel method that robustly transfers the harmonic basis during remeshing. While our experimental evaluation focuses on 2D surface flows, the underlying theoretical formulation is presented to capture the 3D setting as well. Our evaluation demonstrates that this adaptive approach accurately recreates the dynamics of high-resolution simulations while reducing the memory footprint by up to 86% and maintaining numerical stability even on poor-quality triangulations where static methods fail.