An adaptive ghost fluid finite volume method for compressible gas–water simulations

• Published in: Journal of computational physics Vol. 227; no. 12; pp. 6385 - 6409
• Main Authors: Wang, Chunwu, Tang, Huazhong, Liu, Tiegang
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier Inc 01.06.2008; Elsevier
• Subjects: Approximate Riemann solver; Computational techniques; Exact sciences and technology; Finite volume method; Gas–liquid Riemann problem; Ghost fluid method; Level-set equation; Mathematical methods in physics; Moving mesh method; Physics; Ghost fluid method; Level-set equation; Finite volume method; Moving mesh method; Gas–liquid Riemann problem; Approximate Riemann solver; Discontinuity; Hamilton-Jacobi equations; Riemann problem; Strong interactions; Digital simulation; Flow pattern; Two dimensional flow; Adaptive method; Finite volume methods; Gas- liquid Riemann problem; Pressure gradients; Calculation methods; Conservation laws; Shock waves; Algorithms; Gas flow; Density gradient; Calculation; Mesh method; Compressible flow; Finite volume method: Ghost fluid method
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2008.03.005

An adaptive ghost fluid finite volume method is developed for one- and two-dimensional compressible multi-medium flows in this work. It couples the real ghost fluid method (GFM) [C.W. Wang, T.G. Liu, B.C. Khoo, A real-ghost fluid method for the simulation of multi-medium compressible flow, SIAM J. Sci. Comput. 28 (2006) 278–302] and the adaptive moving mesh method [H.Z. Tang, T. Tang. Moving mesh methods for one- and two-dimensional hyperbolic conservation laws, SIAM J. Numer. Anal. 41 (2003) 487–515; H.Z. Tang, T. Tang, P.W. Zhang, An adaptive mesh redistribution method for non-linear Hamilton–Jacobi equations in two- and three-dimensions, J. Comput. Phys. 188 (2003) 543–572], and thus combines their advantages. This work shows that the local mesh clustering in the vicinity of the material interface can effectively reduce both numerical and conservative errors caused by the GFM around the material interface and other discontinuities. Besides the improvement of flow field resolution, the adaptive GFM also largely increases the computational efficiency. Several numerical experiments are conducted to demonstrate robustness and efficiency of the current method. They include several 1D and 2D gas–water flow problems, involving a large density gradient at the material interface and strong shock-interface interactions. The results show that our algorithm can capture the shock waves and the material interface accurately, and is stable and robust even for solutions with large density and pressure gradients.