A flux-split algorithm applied to conservative models for multicomponent compressible flows
A flux-split algorithm applied to conservative models for multicomponent compressible flows
复制标题
DOI：
10.1016/s0021-9991(02)00050-5
复制
发表时间：
2003-02
影响因子：
4.1
通讯作者：
A. Marquina;P. Mulet
中科院分区：
文献类型：
--
作者：
A. Marquina;P. Mulet
本文考虑了气体动力学欧拉方程的保守推广，以在笛卡尔坐标系下描述双组分可压缩流动。众所周知，应用于保守模型的经典激波捕获方案在两种气体之间的界面附近是振荡的。一些作者已经解决了这个问题，提出了一种原始的一致性算法[J]。第一版。拉格朗日成分（鬼流体法，Fedkiw等）[J]。第一版。科学通报，2004(2):1 - 4。第一版。物理学报，169(2001):594]。由于第一作者（参见[J]）的缘故，我们用一种通量分裂算法直接求解这个保守模型。第一版。物理学报，125(1996)42])，以及高阶（WENO5）通量重建[J]。第一版。物理学报115 (1994)200；[32]。该算法似乎以一种不影响实验物理的方式减少了界面附近的振荡。在Haas和Sturtevant研究的相同条件下，我们用1.22马赫激波撞击空气中的氦气泡的相互作用的数值模拟来验证我们的算法。流体力学，181(1987)41]，由Quirk和Karni成功模拟[J]。流体力学。318(1996)129]。
In this paper we consider a conservative extension of the Euler equations for gas dynamics to describe a two-component compressible flow in Cartesian coordinates. It is well known that classical shock-capturing schemes applied to conservative models are oscillatory near the interface between the two gases. Several authors have addressed this problem proposing either a primitive consistent algorithm [J. Comput. Phys. 112 (1994) 31] or Lagrangian ingredients (Ghost Fluid Method by Fedkiw et al. [J. Comput. Phys. 152 (1999) 452] and [J. Comput. Phys. 169 (2001) 594]). We solve directly this conservative model by a flux-split algorithm, due to the first author (see [J. Comput. Phys. 125 (1996) 42]), together with a high-order (WENO5) flux reconstruction [J. Comput. Phys. 115 (1994) 200; 83 (1989) 32]. This algorithm seems to reduce the oscillations near the interfaces in a way that does not affect the physics of the experiments. We validate our algorithm with the numerical simulation of the interaction of a Mach 1.22 shock wave impinging a helium bubble in air, under the same conditions studied by Haas and Sturtevant [J. Fluid Mech. 181 (1987) 41] and successfully simulated by Quirk and Karni [J. Fluid Mech. 318 (1996) 129].