Multigrid methods are quite efficient for solving the pressure Poisson equation in simulations of incompressible flow. However, for viscous liquids, geometric multigrid turned out to be less efficient for solving the variational viscosity equation. In this contribution, we present an Unsmoothed Aggregation Algebraic MultiGrid (UAAMG) method with a multi-color Gauss-Seidel smoother, which consistently solves the variational viscosity equation in a few iterations for various material parameters. Moreover, we augment the OpenVDB data structure with Intel SIMD intrinsic functions to perform sparse matrix-vector multiplications efficiently on all multigrid levels. Our framework is 2.0 to 14.6 times faster compared to the state-of-the-art adaptive octree solver in commercial software for the large-scale simulation of both non-viscous and viscous flow.
Bibliographical noteKAUST Repository Item: Exported on 2022-09-14
Acknowledgements: This work has been supported and funded by KAUST through individual baseline funding. The authors are grateful to Xinxin Zhang for providing the code of the first version of the unsmoothed aggregation-based multigrid Poisson solver. The anonymous reviewers’ valuable comments that improved the manuscript are gratefully acknowledged.
ASJC Scopus subject areas
- Computer Graphics and Computer-Aided Design