【摘 要】
：
冲击波加速不同流体间带有预制初始扰动的密度间断界面会导致Richtmyer-Meshkov (RM)不稳定性发生，后期又诱发湍流混合。该问题在诸如惯性约束聚变、超音速燃烧、天体物理等领域有广泛应用，所以具有重要的研究价值。本文利用可压缩多介质粘性流动和湍流大涡模拟程序MVFT(multi-viscousflow and turbulence)，数值模拟研究了带有反射冲击的多模态RM不稳定性及其诱发
【机 构】
：
中国工程物理研究院流体物理研究所,四川绵阳621900;南京理工大学能源与动力工程学院,江苏南京210014
论文部分内容阅读
冲击波加速不同流体间带有预制初始扰动的密度间断界面会导致Richtmyer-Meshkov (RM)不稳定性发生，后期又诱发湍流混合。该问题在诸如惯性约束聚变、超音速燃烧、天体物理等领域有广泛应用，所以具有重要的研究价值。本文利用可压缩多介质粘性流动和湍流大涡模拟程序MVFT(multi-viscousflow and turbulence)，数值模拟研究了带有反射冲击的多模态RM不稳定性及其诱发的湍流混合演化的动力学特性。模拟结果显示湍流混合区宽度在RM不稳定性发展的不同阶段具有不同的增长方式，混合区湍动能、耗散率、拟涡能以和混合区宽度增长规律相同的方式衰减，说明湍流混合区发展具有统计相似性。湍流混合区脉动量的高阶矩的演化表明RM不稳定性在前期发展中具有强烈的各向异性，随着混合程度的加深，湍流混合区发展逐渐表现出各向同性特征，并出现波数k-5/3标度关系。虽然湍流混合区中心区域脉动强度最大，但是混合程度却最低。
其他文献
本报告基于曲面形态连续介质一般运动的涡量动力学理论，数值研究固定曲面上二维不可压缩与可压缩流动，运动曲面上的二维流动。首先，数值研究固定曲面上的不可压缩流动。具体事例包括含局部高斯凸起曲面的自由流，以及曲面上圆柱绕流。数值结果发现，曲率可以产生涡量，并显著改变圆柱尾迹形态。本文对雷诺数、曲面凸起形态、位置、面摩擦系数等进行了系统的参数演化研究，并对流场空间动力学行为进行了分析，包括涡街形态、最大拉
大气边界层处于对流层底部，厚度一般为数百米到一两公里.由于受到地面的摩擦和热力作用，大气边界层内大气的流动形式为湍流，湍流造成了垂直方向上的热量、动量、物质的交换，是地面和高层大气联系的纽带.早期，人们多认为湍流是不规则的随机过程，上世纪60年代在壁面湍流边界层发现大尺度拟序结构，改变了人们对湍流的认识，而就在稍早的50年代末在大气边界层中发现了对流云街，一种轴向沿近顺风方向的滚轴涡，是最早在大气
为探讨尾迹与平板层流边界层干涉而导致的流动转捩现象,将方柱垂直来流平行放置在平板前缘上方(平板长度58D,方柱中心距平板5D,D为柱体边长),采用独立发展的大涡模拟程序对来流Ma=0.3,雷诺数ReD=3900条件下的流动进行了细致模拟.对结果的初步分析表明,计算实现了层流边界层在方柱尾迹非定常扰动下旁路转捩为湍流全过程的仿真,从20D开始,50D处结束,转捩区长度约30D.尾迹作用下,层流边界层
钝体与端壁连接位置所形成的角区流动(Junction flow)在流体力学领域广泛存在,其时空演化过程具有强三维性和强非定常性.为揭示柱体根部复杂角区流动的物理图像,运用独立发展的大涡模拟程序对一侧存在端壁边界层、高宽比h/d=8的方柱绕流进行了细致模拟.计算的来流马赫数Ma=0.5,雷诺数Red=12000,采用Rescaling-Reintroducing (R-R)方法生成端壁的进口湍流边界
前台阶或后台阶流动都是经典的分离流动，也是工程中的常见流动形态.为考察前/后台阶对平板湍流边界层发展的影响，本文采用大涡模拟方法对来流马赫数0.3条件下的平板、前台阶和后台阶流动进行了对比研究.三个算例的计算域尺寸一致，进口湍流边界层由Rescaling-Reintroducing (R-R)方法生成，保证具有相同的来流条件，前/后台阶的流向位置均设置在平板算例Reθ=1200处，台阶高度h=5θ
飞机座舱内空气循环系统是保障乘客舒适性、降低污染物扩散的重要环节，是民航客机重要考量设计之一。双向送风系统是一种典型的座舱内送风结构，会在座舱内的流动中出成了对扰动及其敏感的异宿轨道。针对座舱内流场进行的实验测量和数值模拟都表明，在出现双向送风速度或舱内几何物形不对称等扰动情况下，异宿轨道将被撕裂，拓扑结构将发生多模态的变化，甚至在座舱不同横截面间出现不同的拓扑模态，以及会出现不同模态间的切换，造
横流失稳是后掠流动中的典型失稳机制，存在定常涡与行进波两种失稳模态。在高空飞行环境中，自由流湍流度较低，定常涡模态主导转捩过程。边界层内的定常涡扰动先后经历线性及非线性增长阶段，非线性作用使扰动幅值达到饱和，但不能促发边界层的转捩。实验及理论研究发现，二次失稳过程是横流失稳导致边界层转捩的关键阶段。二维全局(bi-global)稳定性分析是研究横流二次失稳的主要理论方法。传统的理论分析框架建立在涡
工程实际遇到的三维边界层中，存在着一类边界层，其展向变化远大于流向变化，但又小于法向变化.如小攻角圆锥的绕流问题，见下图1.在背风面的附近存在流向涡，用LST方法预测扰动的增长率，与数值模拟比较，误差较大.对于小攻角组合体飞行器的扰动，在迎风面对称轴附近，用流动稳定性(LST)进行分析，可以发现，存在展向较窄的增长区，如图2.在迎风面对称轴附近用LST方法预测扰动的增长率与数值模拟比较的误差也较大
高超声速飞行器前体压缩面加装强制转捩装置是保证流动进入进气道前转捩为湍流的有效手段，不同的粗糙单元外形促进转捩的效果各有差异，其带来的额外阻力和热流增量也不同，因此在确定外形时需要权衡多方面因素.Berry[1]等对不同外形的边界层强制转捩装置(菱形、斜坡型和三角锥型等)进行了实验研究，最终选出了一种后掠斜坡型粗糙单元作为试飞采用的转捩装置.张明华等[2]的研究结果表明，对于斜坡型粗糙单元，几何参
本文利用可压缩多介质粘性流动和湍流大涡模拟程序MVFT-2D(Multi-viscous Flow and Turbulence)，有效地研究了初始非均匀流场密度为高斯分布，激波Ma=1.27情况下，高、低密度区中多组不同初始振幅组合扰动对双模态Richtmyer-Meshkov不稳定性发展的影响.数值模拟结果表明，处于非均匀流场低密度区的初始大振幅界面的扰动增长最快，高密度区的初始小振幅界面扰动