Numerical investigation on the influence of gas-particle two-way coupling to the shock fluid in the two-dimensional Lagrangian framework

Chen Da-Wei^{1,2}, Sun Hai-Quan^{1}, Wang Pei^{1}, Yu Xi-Jun^{1}, Ma Dong-Jun^{1}

1. National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100094, China;
2. Graduate School of China Academy of Engineering Physics, Beijing 100088, China

When an extreme shock wave releases from the free surface of the material, some high speed particulate matters will be ejected from the material body and enter into the background gas. This induced multiphase mixing phenomenon is known as the ejecta mixing. Ejecta mixing is one of the most important problems in the scope of inner explosive compression engineering, and it is also a frontier research subject of the impact dynamics, multiphase fluid dynamics, computational mathematics, etc. The properties of ejecta mixing have been investigated experimently and analytically for many years. However, the results of numerical simuliation are very rare. At present, the ejecta mixing study mainly focuses on the gas particle one-way coupling, that is, the interests of existing works are in the characteristics of the ejected particulate matter transport in the gas. In fact, after a large number of particles entering into the gas, the gas and the particles will interact with each other. So it is necessary to consider the feedback of particles to the gas. In this paper, the theoretical modeling of gas particle two-way coupling, the discrete algorithm of the mathematical model and the particle phase feedback effects on the gas shock wave propagation are investigated in the framework of Lagrangian coordinates. In order to obtain the details of ejecta movement, the particle trajectory model is chosen as the basic model, and then the governing equations including interactions between gas phase and particle phase are derived. For giving the specific calculation formula, the physical meanings of the coupled interaction source terms in the Lagrangian framework are analyzed and a stable numerical scheme is given based on the staggered strategy. We also devise two different computing models of ejecta mixing, the planar and the column configurations, and then the numerical simulations are carried out. The phenomenon of gas shock speed acceleration caused by particle feedback is found, and the distributions of the physical quantities, such as density, velocity, specific internal energy, pressure, in the gas area are changed. Especially for the convergent configuration, the feedback effects will be amplified further by the geometrical shrinkage, which may have a significant influence on the performance of the inner explosion compression, owing to the obvious uniformity variation of the gas flow field and the gas shock rebound in advance. The mathematical model, the numerical method and the new physical findings in this paper, will provide an important theoretical support for the in-depth understanding of the ejecta mixing and also for the solving of the corresponding engineering application problems.

(Mechanical effects of light on material media, microstructures and particles)

Fund:

Project supported by the National Natural Science Foundation of China (Grant Nos. U1530261, 11571002), the Science Foundation of China Academy of Engineering Physics, China (Grant Nos. 2015B0101021, 2015B0201043, 2013A0202011), the National Defense Basic Scientific Research Program of China (Grant No. B1520133015), and the Foundation of State Key Laboratory of Computational Physics, China (Grant No. 42601-03-02).

Chen Da-Wei,Sun Hai-Quan,Wang Pei et al.. Numerical investigation on the influence of gas-particle two-way coupling to the shock fluid in the two-dimensional Lagrangian framework[J]. Acta Physica Sinica, 2016, 65(8):
.
doi:10.7498/aps.65.084703.

Wang P, Sun H Q, Shao J L, Qin C S, Li X Z 2012 Acta Phys. Sin. 61 234703 (in Chinese) [王裴, 孙海权, 邵建立, 秦承森, 李欣竹 2012 物理学 报 61 234703]

[3]

Zellner M B, McNeil W V, Hammerberg J E, Hixson R S, Obst A W, Olson R T, Payton J R, Rigg P A, Routley N, Stevens G D, Turley W D, Veeser L, Buttler W T 2008 J. Appl. Phys. 103 123502

[4]

Buttler W T, Oró D M, Preston D L, Mikaelian K O, Cherne F J, Hixson R S, Mariam F G, Morris C, Stone J B, Terrones G, Tupa D 2012 J. Fluid Mech. 703 60

[5]

Fung J, Harrison A K, Chitanvis S, Margulies J 2013 Comput. Fluids 83 177

[6]

Sun H Q, Wang P, Chen D W, Qin C S 2014 Explosion and Shock Waves 34 392 (in Chinese) [孙海权, 王裴, 陈大伟, 秦承森 2014 爆炸与冲击 34 392]

[7]

Zhao X W, Li X Z, Wang X J, Song P, Zhang H Z, Wu Q 2015 Acta Phys. Sin. 64 124701 (in Chinese) [赵信文, 李欣竹, 王学军, 宋萍, 张汉钊, 吴 强 2015 物理学报 64 124701]

[8]

Ye Y, Li J, Zhu P F, Qian W X, Liu Z Q, Wang X, Li X Z, Li Z Y, Li Z R, Zhong J 2013 Chin. J. High Pressure Phys. 27 398 (in Chinese) [叶雁, 李军, 朱鹏飞, 钱伟新, 刘振清, 王晓, 李欣竹, 李作友, 李泽仁, 钟杰 2013 高压物理学报 27 398]

[9]

Chen Y T, Hu H B, Tang T G, Ren G W, Li Q Z, Wang R B, Buttler W T 2012 J. Appl. Phys. 111 053509

[10]

Sorenson D S, Minich R W, Romero J L, Tunnell T W, Malone R M 2002 J. Appl. Phys. 92 5830

[11]

Elias P, Chapron P, Mondot M 1989 In Shock Compression of Condensed Matter (Elsevier Science Publishers) p783

[12]

Ogorodnikov V A, Ivanov A G, Mikhailov A L 1998 Combust. Explo. Shock 34 696

[13]

Oró D M, Hammerberg J E, Buttler W T, Mariam F G, Morris C, Rousculp C, Stone J B 2012 AIP Conf. Proc. 1426 1351

[14]

Furnish M D 2012 Sandia Laboratories SAND2012-1552 C

[15]

Rousculp C L, Oro D M, Morris C, Saunders A, Reass W, Griego J R, Turchi P J, Reinovsky R E 2015 Damaged Surface Hydrodynamics (DSH) Flash Report LA-UR-15-22889

[16]

Hu M B, Dang S C, Ma Q, Xia W D 2015 Chin. Phys. B 24 074502

[17]

Gu X 2004 M. S. Dissertation (Harbin: Harbin Engineering University) (in Chinese) [顾旋 2004 硕士学位论文(哈尔滨: 哈尔滨工程大学)]

[18]

Xiong Q G, Aramideh S, Passalacqua A, Kong S C 2015 J. Heat. Trans. 137 061008

[19]

Yu H 2003 Explosion and Shock Waves 23 493 (in Chinese) [喻虹 2003 爆炸与冲击 23 493]

[20]

Moshfegh A, Shams M, Ebrahimi R, Farnia M A 2009 Int. J. Heat Fluid Fl. 30 1142