二维拉格朗日坐标系下气粒混合双向耦合对激波流场影响的计算  

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

作　　者：陈大伟[1,2] 孙海权[1] 王裴[1] 蔚喜军[1] 马东军[1] 

机构地区：[1]北京应用物理与计算数学研究所,计算物理国家重点实验室,北京100094 [2]中国工程物理研究院研究生院,北京100088

出　　处：《物理学报》2016年第8期256-264,共9页Acta Physica Sinica

基　　金：国家自然科学基金(批准号:U1530261,11571002);中国工程物理研究院科学技术发展基金(批准号:2015B0101021,2015B0201043,2013A0202011);国防基础科研计划(批准号:B1520133015);计算物理重点实验室基金(批准号:42601-03-02)资助的课题~~

摘　　要：喷射颗粒与气体混合是内爆压缩领域的热点和难点.针对喷射混合中的气粒双向耦合问题,开展了理论建模、离散算法以及颗粒反馈对激波流场的影响研究.建立了拉格朗日计算框架下的数学模型;给出了耦合源项的离散算法;开展了平面及汇聚构型条件下,气粒双向耦合的数值模拟研究;发现了颗粒反馈导致气体激波提速现象以及气区流场物理量分布形态的改变,初步获得了量化分析结果.本文建立的数学模型、计算方法和获得的新的物理认识,为深入理解喷射混合现象、解决相关工程应用问题提供了重要理论支撑.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, su

关 键 词：内爆混合 拉格朗日坐标系 颗粒轨道模型 双向耦合 

分 类 号：O354.5[理学—流体力学]