TY - GEN
T1 - A fourth-order accurate compact difference scheme for solving the three-dimensional poisson equation with arbitrary boundaries
AU - Hosseinverdi, Shirzad
AU - Fasel, Hermann F.
N1 - Funding Information:
This work was supported by the Air Force Office of Scientific Research (AFOSR) under grant numbers FA9550-14-1-0184 & FA9550-19-1-0174, with Drs. D. Smith & G. Abate serving as the program managers, respectively.
Publisher Copyright:
© 2020, American Institute of Aeronautics and Astronautics Inc, AIAA. All rights reserved.
PY - 2020
Y1 - 2020
N2 - This paper presents an efficient high-order sharp-interface method for solving the three-dimensional (3D) Poisson equation with Dirichlet boundary conditions on a nonuniform Carte-sian grid with irregular domain boundaries. The new approach is based on the combina-tion of the fourth-order compact finite difference scheme and the preconditioned stabilized biconjugate-gradient (BiCGSTAB) method. Contrary to the original immersed interface method by LeVeque & Li [1], the new method does not require jump corrections, instead, the (regular) compact finite difference stencil is adjusted at the irregular grid points (in the vicinity of the interfaces of the immersed bodies) to obtain a solution that is sharp across the interface while keeping the fourth-order global accuracy. The contribution of the present work is the design of a fourth-order-accurate 3D Poisson solver whose accuracy and efficiency does not deteriorate in the presence of an immersed boundary. This is attributed to (i) the modification of the discrete operators near immersed boundaries does not lead to a wide grid stencil due to the compact nature of the discretization and (ii) a preconditioning technique whose efficiency and cost are independent of the complexity of the geometry and the presence or not of an immersed boundary. The accuracy and computational efficiency of the proposed algorithm is demonstrated and validated over a range of problems including smooth and ir-regular boundaries. The test cases show that the new method is fourth-order accurate in the maximum norm whether an immersed boundary is present or not, on uniform or nonuniform grids. Furthermore, the efficiency of the preconditioned BiCGSTAB is demonstrated with re-gard to convergence rate and “extra” floating-point operation (FLOPextra ) which is due to the presence of immersed boundaries. It is shown that the solution method is equally efficient for domains with and without irregular boundaries, with a negligible FLOPextra in the presence of immersed boundaries.
AB - This paper presents an efficient high-order sharp-interface method for solving the three-dimensional (3D) Poisson equation with Dirichlet boundary conditions on a nonuniform Carte-sian grid with irregular domain boundaries. The new approach is based on the combina-tion of the fourth-order compact finite difference scheme and the preconditioned stabilized biconjugate-gradient (BiCGSTAB) method. Contrary to the original immersed interface method by LeVeque & Li [1], the new method does not require jump corrections, instead, the (regular) compact finite difference stencil is adjusted at the irregular grid points (in the vicinity of the interfaces of the immersed bodies) to obtain a solution that is sharp across the interface while keeping the fourth-order global accuracy. The contribution of the present work is the design of a fourth-order-accurate 3D Poisson solver whose accuracy and efficiency does not deteriorate in the presence of an immersed boundary. This is attributed to (i) the modification of the discrete operators near immersed boundaries does not lead to a wide grid stencil due to the compact nature of the discretization and (ii) a preconditioning technique whose efficiency and cost are independent of the complexity of the geometry and the presence or not of an immersed boundary. The accuracy and computational efficiency of the proposed algorithm is demonstrated and validated over a range of problems including smooth and ir-regular boundaries. The test cases show that the new method is fourth-order accurate in the maximum norm whether an immersed boundary is present or not, on uniform or nonuniform grids. Furthermore, the efficiency of the preconditioned BiCGSTAB is demonstrated with re-gard to convergence rate and “extra” floating-point operation (FLOPextra ) which is due to the presence of immersed boundaries. It is shown that the solution method is equally efficient for domains with and without irregular boundaries, with a negligible FLOPextra in the presence of immersed boundaries.
UR - http://www.scopus.com/inward/record.url?scp=85091928965&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85091928965&partnerID=8YFLogxK
U2 - 10.2514/6.2020-0805
DO - 10.2514/6.2020-0805
M3 - Conference contribution
AN - SCOPUS:85091928965
SN - 9781624105951
T3 - AIAA Scitech 2020 Forum
BT - AIAA Scitech 2020 Forum
PB - American Institute of Aeronautics and Astronautics Inc, AIAA
T2 - AIAA Scitech Forum, 2020
Y2 - 6 January 2020 through 10 January 2020
ER -