Sawtooth cycle revisited

Solving the pressure Poisson equation dominates a large portion of computational time for incompressible fluid flow simulations. To solve the pressure Poisson equation efficiently, geometric multigrid methods are used directly or used as the preconditioner for the conjugate gradient (CG) method. Conventionally, the V‐cycle multigrid method is widely employed, and little attention has been paid to other cycles. In this paper, we introduce the sawtooth cycle multigrid method and its simple extension called N‐cycle, which provides better convergence in equal time comparison with the V‐cycle as a direct solver of the pressure Poisson equation. We also apply the N‐cycle to the preconditioner of the CG method and show that the N‐cycle multigrid CG method can provide better convergence than the V‐cycle multigrid CG method.

[1]  Matthias Müller,et al.  Real-time Eulerian water simulation using a restricted tall cell grid , 2011, SIGGRAPH 2011.

[2]  Greg Humphreys,et al.  A multigrid solver for boundary value problems using programmable graphics hardware , 2003, HWWS '03.

[3]  Eftychios Sifakis,et al.  Power diagrams and sparse paged grids for high resolution adaptive liquids , 2017, ACM Trans. Graph..

[4]  SzeliskiRichard,et al.  Multigrid and multilevel preconditioners for computational photography , 2011 .

[5]  GrinspunEitan,et al.  Sparse matrix solvers on the GPU , 2003 .

[6]  Satteluri R. K. Lyengar,et al.  Comparison of S and V cycles in multigrid method for linear elliptic equations with variable coefficients , 1992 .

[7]  Eitan Grinspun,et al.  Sparse matrix solvers on the GPU: conjugate gradients and multigrid , 2003, SIGGRAPH Courses.

[8]  Ming C. Lin,et al.  A Multilevel SPH Solver with Unified Solid Boundary Handling , 2016, Comput. Graph. Forum.

[9]  Richard Szeliski,et al.  Multigrid and multilevel preconditioners for computational photography , 2011, ACM Trans. Graph..

[10]  Dieter W. Fellner,et al.  A Cut‐Cell Geometric Multigrid Poisson Solver for Fluid Simulation , 2015, Comput. Graph. Forum.

[11]  Matthias Müller,et al.  A Multigrid Fluid Pressure Solver Handling Separating Solid Boundary Conditions , 2011, IEEE Transactions on Visualization and Computer Graphics.

[12]  Wei-Wen Feng,et al.  A fast multigrid algorithm for mesh deformation , 2006, ACM Trans. Graph..

[13]  Michael M. Kazhdan,et al.  Streaming multigrid for gradient-domain operations on large images , 2008, ACM Trans. Graph..

[14]  Stephen F. McCormick,et al.  Smoothed aggregation multigrid for cloth simulation , 2015, ACM Trans. Graph..

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

[16]  Kun Zhou,et al.  Out-of-core multigrid solver for streaming meshes , 2009, ACM Trans. Graph..

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

[18]  Eftychios Sifakis,et al.  SPGrid: a sparse paged grid structure applied to adaptive smoke simulation , 2014, ACM Trans. Graph..

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

[20]  Jun-yong Noh,et al.  A heterogeneous CPU–GPU parallel approach to a multigrid Poisson solver for incompressible fluid simulation , 2013, Comput. Animat. Virtual Worlds.