激波-界面作用产生射流的计算模拟研究.doc

第八届全国流体力学学术会议 2014年9月18～21日 甘肃兰州 文章编号： CSTAM2014-B01-0066 标题：激波-界面作用产生射流的计算模拟Copyright ? 2014 版权所有 中国力学学会地址: 北京市北四环西路15号　邮政编码:100190 　Address: No.15 Beisihuanxi Road, Beijing 10激波-界面作用产生射流的计算模拟在很多科学和工程问题中，可以观察到两种密度比很高的流体（如气体-液体、气体-金属）的自由界面上产生射流。对射流研究的实验包括表面张力驱动的气泡收缩tler et al. 2012)或者气体-液体界面(Tagawa et al. 2012)，可以产生较强的高速射流。 激波-界面作用产生射流的发生发展整个过程和Richtmyer-Meshkov（R-M）、Kelvin-Helmholtz(K-Hateau-Rayleigh甚至-T)等当激波加载到自由界面上，界面上的扰动会增长，不稳定性。R-M不稳定性和R-T不稳定性在很多方面相似，其小扰动线性解也和R-T不稳定性的处理一样。R-M不稳定性和R-T不稳定性的根本差别在于激波的存在。激波会有穿透、反射等活动，使得对R-M不稳定性的理论分析变得尤其困难uttler等人通过实验研究气体-金属界面的射流过程，并试图用R-M不稳定性解释，发现在空泡速度上不符。 本文将用基于切割网格和vel Set方法的陡峭界面可压缩多相流计算方法（ang, Deng Theofanous 2013）对激波-界面作用产生射流的问题进行模拟研究。该方法完全满足守恒性和相容性；采用曲线切割网格线和适当的在界面上完全满足压力、表面张力、粘性力等各种力在gawa实验中，界面是接触角决定的圆弧形，。首先和普遍测试过的气体-气体实验结果进行对比，验证算法的有效性；再，获取从低速到高速的射流速度曲线，研究激波强度、流体密度、可压缩性、界面曲率、气体压力等不同参数对射流速度和形态的影响。Richtmyer-Meshkov不稳定性射流问题在自然界和工程问题中广泛存在，其对应的物理机制也多种多样。对于激波冲击界面产生射流的这类问题，Richtmyer-Meshkov（R-M）的情况（发生相逆转）。 在R-M不稳定性研究中，若足够强的激波由高密度物质传向低密度物质的时候，在界面上会形成射流现象。Tagawa et al.(2012)和(2013)分别从实验和计算方面研究了脉冲作用下弯曲液体界面的微射流。Buttler et al. (2012) 给出了一些由R-M不稳定性主导的金属界面的微射流的实验数据，Dimonte et al.(2013) 详细比对了实验和数值结果，这方面的理论模型由Guy Dimonte 和 P.Ramaprabhu (2010)给出。 本文采用基于切割网格和Level Set 方法的陡峭界面可压缩多相流计算流体力学方法（Chang，Deng Theofanous 2013），通过数值模拟对R-M不稳定性和激波加载到自由界面上产生射流的现象进行模拟研究。该方法完全满足守恒性和相容性；采用曲线切割网格线和适当的在界面上完全满足压力、表面张力、粘性力等各种力ptive Mesh Refinement, AMR）技术提高效率；在最密的一层网格上，应用切割网格方法产生非结构网格，自由界面由一系列的曲线网格线描述。对界面的处理，我们使用Level set 函数帮助演化界面： (1) 其中，是Level set值，表示带符号的点到界面的距离；是界面速度。界面离散格式是高阶迎风中心差分法HOUC3，可以达到3阶精度，下标i表示离散点的位置： (2) 对流场的处理我们使用有限体积法处理完整的 Navier-Stokes (N-S) 方程： (3) 其中： ， (4) 是流场速度，是控制体边界的单位外法向量，是边界速度，是流体密度，总能量 ，粘性应力张量是 ，是传热量，，是粘性和传热系数，。我们用的状态方程是刚性气体状态方程（stiffened-gas model），如下 (5) 对时间项的离散，我们使用的是4阶Runge-Kutta方法。在数值计算到边界及其附近时，我们针对不同问题分别使用了边界对称条件，无粘边界条件。在自由界面垂直于边界时，如air-SF6(Zhang Sohn 1996,1997; Yang et al.1994; Ho