A two-dimensional geometric multigrid model for Poisson equation with interface on structured adaptive mesh refinement grid. (June 2021)

• Record Type: Journal Article
• Authors: Zhang, Yunxing; Ma, Shan; Liao, Kangping; Duan, Wenyang
• Is Part Of: Applied ocean research. Volume 111(2021)
• Publication Date: 2021-06
• Subjects: Poisson equation -- Interface flows -- Structured adaptive mesh refinement -- Geometric multigrid method -- Galerkin coarse grid approximation; Ocean engineering -- Periodicals; 620.416205
• DOI: 10.1016/j.apor.2021.102655
• Languages: English
• ISSNs: 0141-1187

Abstract

Highlights: A geometric multigrid model for Poisson equation with interface on structured adaptive mesh refinement grid is developed. Galerkin coarse grid approximation method is adopted in generating coefficients on coarse grid to enhance the robustness and efficiency of the model, and a simple form of the method is given which is easy to program. The model is demonstrated to be 2nd-order accuracy with acceptable efficiency for even the density ratio reaches 10 4 . 4 The model is implemented in a developed Navier-Stokes solver to validate its performance in simulating physical problems. Abstract: In this study, a two-dimensional Geometric Multigrid (GMG) model for Poisson equation with interface on Structured Adaptive Mesh Refinement (SAMR) grid is developed. The model is designed to be combined with Navier-Stokes solvers with staggered arrangements of variables, although solvers with collocated arrangement are also applicable. Based on a general GMG method, special attention for flux conservation is taken on the coarse-fine interface. As large density ratios commonly exist in interface flows such as free surface flows, Galerkin Coarse grid Approximation (GCA) method is adopted to generate coefficients on coarse grids to enhance the robustness and efficiency of the model. Benchmark case is carried out for the validation of the model, and it is demonstrated to obtain 2nd-order accuracy and acceptable efficiency for even the density ratio reaches 10 4 . Point iteration and line relaxation iteration of Gauss-Seidel methods are compared to obtain a better performance. Furthermore, the model is implemented in a developed Navier-Stokes solver to validate its performance in simulating interface problems. The numerical results are compared with theoretical solution or experimental data from reliable sources.