Abstract
We present an embedded ghost-fluid method for numerical solutions of the compressible Navier Stokes (CNS) equations in arbitrary complex domains. A PDE multidimensional extrapolation approach is used to reconstruct the solution in the ghost-fluid regions and imposing boundary conditions on the fluid-solid interface, coupled with a multi-dimensional algebraic interpolation for freshly cleared cells. The CNS equations are numerically solved by the second order multidimensional upwind method. Block-structured adaptive mesh refinement, implemented with the Chombo framework, is utilized to reduce the computational cost while keeping high resolution mesh around the embedded boundary and regions of high gradient solutions. The versatility of the method is demonstrated via several numerical examples, in both static and moving geometry, ranging from low Mach number nearly incompressible flows to supersonic flows. Our simulation results are extensively verified against other numerical results and validated against available experimental results where applicable. The significance and advantages of our implementation, which revolve around balancing between the solution accuracy and implementation difficulties, are briefly discussed as well.
Original language | English (US) |
---|---|
Pages (from-to) | 339-378 |
Number of pages | 40 |
Journal | Journal of Computational Physics |
Volume | 337 |
DOIs | |
State | Published - Feb 25 2017 |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Acknowledged KAUST grant number(s): URF/1/1394-01
Acknowledgements: This research was supported by the KAUST Office of Competitive Research Funds under Award No. URF/1/1394-01. The Cray XC40, Shaheen II, at KAUST was utilized for some of the simulations. We thank Dr. Narsimha Rapaka for his assistance in computing the three-dimensional cases on Shaheen II.