Richtmyer-Meshkov不稳定性引起湍流混合的模型研究
作者单位：北京应用物理与计算数学研究所
会议名称：《第十四届全国激波与激波管学术会议》
会议日期：2010年
学科分类：080704[工学-流体机械及工程] 080103[工学-流体力学] 08[工学] 0807[工学-动力工程及工程热物理] 0801[工学-力学（可授工学、理学学位）]
基 金：国家自然科学基金委员会-中国工程物理研究院联合基金“界面运动不稳定性及湍流混合研究”(NSAF10676005) 中国工程物理研究院科学技术发展基金“界面不稳定性诱导湍流混合的多相湍流模型研究”(课题批准号:2008B0101006)的资助下完成的
关 键 词：Richtmyer-Meshkov不稳定性 湍流混合 湍流模型 数值模拟 湍流产生项
摘 要：将传统的k-ε模型应用到存在激波的Richtmyer-Meshkov不稳定性(RMI)引起混合的研究中,采用代数关系式补充模型的封闭假设,给出了合理描述RMI特征的湍动能产生项的表达式。对封闭假设和初始化过程进行了详细研究。采用具有总变差递减保持性质的三阶Runge-Kutta方法进行时间积分,对流项采用高阶WENO方法重构网格边界的通量,扩散项采用交替方向隐式迭代(ADI)法求解。在此基础上,对法国原子能委员会(CEA)***等的激波管实验进行了模拟。该实验首次采用激光多普勒测速仪(LDA)测量了激波管中产生于冲击波的气相混合物的瞬态速度,给出了混合区与反射激波相互作用前后混合区速度脉动的发展。它第一次给出了湍流的直接测量,为校核湍流混合模型提供了重要数据。研究发现:计算结果与实验吻合得较好,且较好地再现了第一次反射激波与湍流混合区相互作用之前湍流的衰减和它被反射激波增强的特征。这些结果表明本研究所采用的模型封闭、模型常数、数值算法和程序实现是合适的。