再冲击载荷作用下流动混合的数值模拟 numerical simulation of flow mixing impacted by reshock.pdfVIP

第29卷第3期 爆炸与冲击 V01．29。No．3 2009年5月 EXPLOSl0NANDSHOCKWAVES May．2009 文章编号：1001—1455(2009)03·0243—06 再冲击载荷作用下流动混合的数值模拟。 王 涛1，柏劲松1’2，李平1’2，钟 敏1 (1．中国工程物理研究院流体物理研究所，四川绵阳621900； 2．中国工程物理研究院流体物理研究所冲击波物理与爆轰物理国防科技重点实验室．四川绵阳621900) of 摘要：基于多介质的流体体积分数VOF(volume method)方法， fluid)方法和PPM(piecewiseparabolic piecewiseparabolicmethod)。为了检验和验证此计算代码，对某激波管实验——再冲击载荷作用下的流体动 力学不稳定性及其导致的混合过程进行了数值模拟。计算结果与实验结果一致。同时还研究了激波反射冲击 作用下流体混合区的演化情况，在反射激波和混合区相互作用的瞬间。混合区的宽度明显减小。之后又迅速增 大；另外。混合增长率与初始扰动的频谱有很大关联。通过对有粘性(分子动力学粘性)和元粘性结果的对比， 发现粘性对混合区的影响很小。 关键词：流体力学；界面不稳定性；MVPPMf流动混合；分子动力学粘性 中图分类号：0357 国标学科代码：130·2531文献标志码：A 1 引 言 分层物质的界面将发生不稳定，他用线性化的Euler方程对交界面的运动进行了数值模拟，同时还给出 征：气泡和尖钉的构型及增长。所谓气泡就是侵入到重物质中的轻物质部分，尖钉就是侵入到轻物质中 的重物质部分，两者向相反的方向增长。同时，R-M不稳定性和R—T不稳定性存在着本质的区别：当加 速度的方向从轻流体指向重流体时，R—T界面是不稳定的；而激波无论是从重流体至轻流体方向还是从 轻流体至重流体方向冲击物质界面，R—M界面总是不稳定的。在冲击波的作用下，界面上的初始扰动 在初始阶段以线性方式增长，后期又以非线性的方式增长，最后导致混合而进入湍流混合阶段。R—M 不稳定性在自然界和工程中经常遇到，由于它在科学研究和工程实践上的重要意义，自被发现以来一直 都是冲击波动力学中一个异常活跃的研究领域。这种冲击波和交界面相互作用的基本过程对研究超音 confinementfusion) 速冲压发动机(seramjet)中的混合和燃烧、激光驱动的惯性约束聚变ICF(inertial 等很有意义，它还可以用来解释天体物理学中的某些现象如超新星爆发，它后期的非线性过程对于揭示 可压缩湍流的本质十分重要。 目前，对激波单次冲击作用下的湍流混合演化已经有了比较全面的认识；但是在多次冲击作用下， 对其发展特征的认识及相关信息的获得还不全面。针对这个问题，本文中在已有的欧拉计算程序MF— PPM(multi-fluid pieeewiseparabolicmethod)t¨]基础上，添加粘性部分的计算(此处主要考虑分子动力 -收稿日期：2007-12—14；修回日期：2008—05-17 中国国防科技工业技术基础质量和可靠性项目 作者简介：王涛(1979一)，男．硕士，助理研究员． 万方数据 244 爆 炸 与 冲 击 第29卷 模拟。对这种情况下R—M不稳定性及其导致的混合过程的发展做了比较细致的研究，计算结果和实验 结果一致性较好，流体混合区的增长对初始条件(初始扰动波长和振幅等)的依赖性很强，但是分子动力 学粘性的影响作用很有限。 2控制方程和计算方法 在考虑热传导和粘性的情况下，以张量形式描述的多介质流场控制方程为 f害+罄一o，警+警一罄+罄，