A schur complement preconditioner for scalable parallel fluid simulation

We present an algorithmically efficient and parallelized domain decomposition based approach to solving Poisson’s equation on irregular domains. Our technique employs the Schur complement method, which permits a high degree of parallel efficiency on multicore systems. We create a novel Schur complement preconditioner which achieves faster convergence, and requires less computation time and memory. This domain decomposition method allows us to apply different linear solvers for different regions of the flow. Subdomains with regular boundaries can be solved with an FFT-based Fast Poisson Solver. We can solve systems with 1,0243 degrees of freedom, and demonstrate its use for the pressure projection step of incompressible liquid and gas simulations. The results demonstrate considerable speedup over preconditioned conjugate gradient methods commonly employed to solve such problems, including a multigrid preconditioned conjugate gradient method.

[1]  C.-C. Jay Kuo,et al.  A Domain Decomposition Preconditioner Based on a Change to a Multilevel Nodal Basis , 1991, SIAM J. Sci. Comput..

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

[3]  James F. O'Brien,et al.  Liquid simulation on lattice-based tetrahedral meshes , 2007 .

[4]  Yongning Zhu,et al.  Animating sand as a fluid , 2005, SIGGRAPH 2005.

[5]  R.D. Falgout,et al.  An Introduction to Algebraic Multigrid Computing , 2006, Computing in Science & Engineering.

[6]  Rüdiger Westermann,et al.  Narrow Band FLIP for Liquid Simulations , 2016, Comput. Graph. Forum.

[7]  Yousef Saad,et al.  ILUT: A dual threshold incomplete LU factorization , 1994, Numer. Linear Algebra Appl..

[8]  Ronald D. Henderson.  Scalable fluid simulation in linear time on shared memory multiprocessors , 2012, DigiPro '12.

[9]  Rüdiger Westermann,et al.  Solving the Fluid Pressure Poisson Equation Using Multigrid—Evaluation and Improvements , 2016, IEEE Transactions on Visualization and Computer Graphics.

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

[11]  Eftychios Sifakis,et al.  Physical simulation for animation and visual effects: parallelization and characterization for chip multiprocessors , 2007, ISCA '07.

[12]  Christopher Batty,et al.  Tetrahedral Embedded Boundary Methods for Accurate and Flexible Adaptive Fluids , 2010, Comput. Graph. Forum.

[13]  Pradeep Dubey,et al.  Large-scale fluid simulation using velocity-vorticity domain decomposition , 2012, ACM Trans. Graph..

[14]  Christopher Wojtan,et al.  Highly adaptive liquid simulations on tetrahedral meshes , 2013, ACM Trans. Graph..

[15]  Rüdiger Westermann,et al.  Large-Scale Liquid Simulation on Adaptive Hexahedral Grids , 2014, IEEE Transactions on Visualization and Computer Graphics.

[16]  Tony F. Chan,et al.  The Interface Probing Technique in Domain Decomposition , 1992, SIAM J. Matrix Anal. Appl..

[17]  Robert Bridson,et al.  A Multipreconditioned Conjugate Gradient Algorithm , 2005, SIAM J. Matrix Anal. Appl..

[18]  Tae-Yong Kim,et al.  Coupling 3D Eulerian, Heightfield and Particle Methods for Interactive Simulation of Large Scale Liquid Phenomena , 2014, Symposium on Computer Animation.

[19]  R. Henderson.  Nonlinear dynamics and pattern formation in turbulent wake transition , 1997, Journal of Fluid Mechanics.

[20]  Gene H. Golub,et al.  A Fast Poisson Solver for the Finite Difference Solution of the Incompressible Navier-Stokes Equations , 1998, SIAM J. Sci. Comput..

[21]  J. Pasciak,et al.  The Construction of Preconditioners for Elliptic Problems by Substructuring. , 2010 .

[22]  Ken Museth,et al.  VDB: High-resolution sparse volumes with dynamic topology , 2013, TOGS.

[23]  James F. O'Brien,et al.  Fluid animation with dynamic meshes , 2006, SIGGRAPH 2006.

[24]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

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

[26]  James Demmel,et al.  Applied Numerical Linear Algebra , 1997 .

[27]  Robert D. Falgout,et al.  Multigrid preconditioned conjugate gradients for the numerical simulation of groundwater flow on the Cray T3D , 1994 .

[28]  Nicole Spillane,et al.  An Adaptive MultiPreconditioned Conjugate Gradient Algorithm , 2016, SIAM J. Sci. Comput..

[29]  Eftychios Sifakis,et al.  A scalable schur-complement fluids solver for heterogeneous compute platforms , 2016, ACM Trans. Graph..

[30]  Ronald Fedkiw,et al.  Chimera grids for water simulation , 2013, SCA '13.

[31]  Nuttapong Chentanez,et al.  A Multigrid Fluid Pressure Solver Handling Separating Solid Boundary Conditions , 2012, IEEE Trans. Vis. Comput. Graph..

[32]  Christopher Wojtan,et al.  Liquid surface tracking with error compensation , 2013, ACM Trans. Graph..

[33]  O. Widlund,et al.  Iterative Methods for the Solution of Elliptic Problems on Regions, Partitioned Into Substructures , 2015 .

[34]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[35]  X. Cai.  Overlapping Domain Decomposition Methods , 2003 .

[36]  Frank Losasso,et al.  Simulating water and smoke with an octree data structure , 2004, SIGGRAPH 2004.

[37]  Robert Bridson,et al.  Fluid Simulation for Computer Graphics , 2008 .

[38]  Jos Stam,et al.  Stable fluids , 1999, SIGGRAPH.

[39]  Yang Yang,et al.  A Fast Iterated Orthogonal Projection Framework for Smoke Simulation , 2016, IEEE Transactions on Visualization and Computer Graphics.

[40]  Christopher Wojtan,et al.  A Practical Method for High‐Resolution Embedded Liquid Surfaces , 2016, Comput. Graph. Forum.

[41]  J. Mandel.  Balancing domain decomposition , 1993 .

[42]  William F. Mitchell,et al.  A Collection of 2D Elliptic Problems for Testing Adaptive Algorithms , 2010 .

[43]  Hyeong-Seok Ko,et al.  Stretching and wiggling liquids , 2009, SIGGRAPH 2009.

[44]  Osamu Tatebe,et al.  The multigrid preconditioned conjugate gradient method , 1993 .

[45]  Robert Bridson,et al.  A fast variational framework for accurate solid-fluid coupling , 2007, SIGGRAPH 2007.

[46]  Tarek P. Mathew,et al.  Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations , 2008, Lecture Notes in Computational Science and Engineering.

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