ACM Transactions on Graphics (SIGGRAPH) 2025

Authors
Bernhard Braun, Technical University of Munich
Jan Bender, RWTH Aachen University
Nils Thuerey, Technical University of Munich

Abstract
Capturing the visually compelling features of large-scale water phenomena,such as the spray clouds of crashing waves, stormy seas, or waterfalls, involves simulating not only the water but also the motion of the air interacting with it. However, current solutions in the visual effects industry still largely rely on single-phase solvers and non-physical “white-water” heuristics. To address these limitations, we present Phase-Field-FLIP (PF-FLIP), a hybrid Eulerian/Lagrangian method for the fully physics-based simulation of very large-scale, highly turbulent multiphase flows at high Reynolds numbers and high fluid density contrasts. PF-FLIP transports mass and momentum in a consistent, non-dissipative manner and, unlike most existing multiphase approaches, does not require a surface reconstruction step. Furthermore, we employ spatial adaptivity across all critical components of the simulation algorithm, including the pressure Poisson solver. We augment PF-FLIP with a dual multiresolution scheme that couples an efficient treeless adaptive grid with adaptive particles, along with a fast adaptive Poisson solver tailored for high-density-contrast multiphase flows. Our method enables the simulation of two-phase flow scenarios with a level of physical realism and detail previously unattainable in graphics, supporting billions of particles and adaptive 3D resolutions with thousands of grid cells per dimension on a single workstation.

Links
Preprint
Code for the MSBG framework
Video

Our simulation of a large dam discharge showing violently turbulent water-air interaction induced by a powerful water jet at high speed (≈ 40𝑚/𝑠). A real photograph of a dam discharge event at similar scale (Three Gorges Dam, China, © Picture Alliance ) is included in the lower right corner.

Necessity of Two-Phase Simulation.
Water simulations in Graphics often neglect the air phase, assuming constant zero air pressure. In these free surface methods, the governing Navier-Stokes equations are thus solved only for the water phase. This simplification, however, becomes insufficient for large-scale phenomena where the violent motion of large masses of water induces strong movements of the surrounding air and aerodynamic forces, which scale with the square of the relative air-water velocity, become dominant. The resulting high Weber numbers inhibit surface tension from stabilizing the liquid, leading to fragmentation into spray droplets. These droplets do not simply follow ballistic trajectories; instead, they form intricate, dynamic patterns driven by the complex interplay of gravity, inertial forces and—crucially—the turbulent air velocity field. Accurately capturing these phenomena requires both modeling the coupled liquid-air velocity and a numerical method powerful enough to resolve scales ranging from tiny spray droplets to large gravity-driven waves.

Fully simulating not only the water but also the surrounding air velocity field is crucial for visually convincing large-scale water simulations: An ocean wave breaking on the beach was simulated with our proposed method. Shown is the two-phase velocity field (vorticity). Upper left corner: The same simulation rendered in 3D. Effective resolution: 4096 × 1024 × 512.
Large waves hitting the coast at high speed. Comparison between our simulation and a real photograph for qualitative reference. Photograph © AP
Photo / Ben Birchall.
Large waves hitting the coast at high speed. Comparison between our simulation and a real photograph for qualitative reference. Photograph © AP Photo / Ben Birchall.