A Multigrid Fluid Pressure Solver Handling Separating Solid Boundary Conditions

We present a multigrid method for solving the linear complementarity problem (LCP) resulting from discretizing the Poisson equation subject to separating solid boundary conditions in an Eulerian liquid simulation's pressure projection step. The method requires only a few small changes to a multigrid solver for linear systems. Our generalized solver is fast enough to handle 3D liquid simulations with separating boundary conditions in practical domain sizes. Previous methods could only handle relatively small 2D domains in reasonable time, because they used expensive quadratic programming (QP) solvers. We demonstrate our technique in several practical scenarios, including nonaxis-aligned containers and moving solids in which the omission of separating boundary conditions results in disturbing artifacts of liquid sticking to solids. Our measurements show, that the convergence rate of our LCP solver is close to that of a standard multigrid solver.

[1]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[2]  J. A. Viecelli,et al.  A method for including arbitrary external boundaries in the MAC incompressible fluid computing technique , 1969 .

[3]  John E. Howland,et al.  Computer graphics , 1990, IEEE Potentials.

[4]  Dimitris N. Metaxas,et al.  Realistic Animation of Liquids , 1996, Graphics Interface.

[5]  M. Ferris,et al.  Complementarity problems in GAMS and the PATH solver 1 This material is based on research supported , 2000 .

[6]  Ronald Fedkiw,et al.  Practical animation of liquids , 2001, SIGGRAPH.

[7]  R. Fedkiw,et al.  Robust Treatment of Interfaces for Fluid Flows and Computer Graphics , 2003 .

[8]  Ben Houston,et al.  A unified approach for modeling complex occlusions in fluid simulations , 2003, SIGGRAPH '03.

[9]  Duc Quang Nguyen,et al.  Directable photorealistic liquids , 2004, SCA '04.

[10]  Joachim Schöberl,et al.  Minimizing Quadratic Functions Subject to Bound Constraints with the Rate of Convergence and Finite Termination , 2005, Comput. Optim. Appl..

[11]  Robert Bridson,et al.  Animating sand as a fluid , 2005, ACM Trans. Graph..

[12]  Erik Hagersten,et al.  Multigrid and Gauss-Seidel smoothers revisited: parallelization on chip multiprocessors , 2006, ICS '06.

[13]  Lin Shi,et al.  A fast multigrid algorithm for mesh deformation , 2006, ACM Trans. Graph..

[14]  James F. O'Brien,et al.  Eurographics/acm Siggraph Symposium on Computer Animation (2007) Liquid Simulation on Lattice-based Tetrahedral Meshes , 2022 .

[15]  Robert Bridson,et al.  A fast variational framework for accurate solid-fluid coupling , 2007, ACM Trans. Graph..

[16]  Ronald Fedkiw,et al.  Two-way coupling of fluids to rigid and deformable solids and shells , 2008, ACM Trans. Graph..

[17]  Robert Bridson,et al.  Accurate viscous free surfaces for buckling, coiling, and rotating liquids , 2008, SCA '08.

[18]  Jonathan M. Cohen,et al.  Low viscosity flow simulations for animation , 2008, SCA '08.

[19]  Matthias Müller,et al.  Hierarchical Position Based Dynamics , 2008, VRIPHYS.

[20]  R. Fedkiw,et al.  A novel algorithm for incompressible flow using only a coarse grid projection , 2010, ACM Trans. Graph..

[21]  Ming C. Lin,et al.  Free-flowing granular materials with two-way solid coupling , 2010, ACM Trans. Graph..

[22]  Eftychios Sifakis,et al.  A parallel multigrid Poisson solver for fluids simulation on large grids , 2010, SCA '10.

[23]  Eftychios Sifakis,et al.  An efficient multigrid method for the simulation of high-resolution elastic solids , 2010, TOGS.

[24]  Matthias Müller,et al.  Real-time Eulerian water simulation using a restricted tall cell grid , 2011, ACM Trans. Graph..