不同初始磁场对激波冲击R22重气柱过程影响的数值模拟

基于磁流体动力(magneto-hydrodynamic,MHD)方程,采用CTU+CT方法,对在不同初始磁场作用下的平面入射激波与磁化R22重质气柱作用过程进行了数值研究。数值结果清晰地描述了不同初始磁场条件时激波诱导R22气柱界面不稳定性的过程,揭示了磁场控制界面不稳定性的机理。另外,还分析了磁感应强度对界面不稳定性的影响,发现在磁场较小时,涡层附着于界面,但随着磁感应强度的增大,平均涡量随之增大,涡层与界面逐渐分离,最终更好地抑制了界面不稳定性。同时,还发现平均涡度拟能随着磁感应强度的增大而减小,而垂直磁场比平行磁场更能降低平均涡度拟能,因而平均涡度拟能可较好地反映磁场对不稳定性的影响效果。     

高超声速自适应激波针数值研究

针对传统的与钝体轴线共线安装的固定式激波针方法在有攻角状态所存在的问题, 在前人工作基础上得到一种新型高超声速飞行器减阻/降热方法------自适应激波针方法. 将该方法应用于三维高超声速轴对称钝锥外形以及扁平楔外形, 并采用数值模拟的方法对其进行了概念验证. 在0○~120○攻角范围内, 对不同L/D参数的激波针外形流场以及前缘壁面的压力、热流分布等进行了对比分析. 结果表明, 这种新型自适应激波针方法无论在无攻角还是有攻角状态, 均可有效降低高超声速飞行器头部壁面的压力和热流, 可以有效解决传统激波针方法在较大攻角情况状态下失效的问题.      

高超声速流动中液滴蒸发对流动参数的影响研究

采用带有非平衡相变模型的数值方法对高超声速流动中液滴蒸发的影响进行了研究,重点探讨了液滴蒸发对斜激波下游区域流场和劈面参数的影响。研究结果表明:液滴在激波后区域的蒸发存在弛豫现象,蒸发弛豫过程与激波强度相关,并影响模型表面气动参数的均匀性;当偏转角从3°增大至16°时,研究发现特别是偏转角为10°时,激波强度刚好能使液滴在整个激波下游区域完全蒸发,此时蒸发引起的流场参数不均性影响最大,与不考虑蒸发时相比,液滴质量含量、压强、静温减小幅度分别达到27.2g/kg、401Pa、54K,该三项参数变化量的比分别达到96.8%、4.8%、14.2%。     

一种改进混合迎风格式在翼型流场激波捕捉中的应用

波阻是飞行器超音速飞行的关键设计因素,精确捕捉激波在流场中的位置,是数值模拟含激波流场和精确计算波阻的一个重要研究内容.本文基于网格节点有限体积空间离散方法,采用AUSM格式与FVS格式的混合格式(MAUSM方法)计算对流通量,从而抑制在数值模拟流场出现的激波处振荡和过冲现象,确保AUSM准确捕获接触间断的特性和FVS格式捕捉激波的能力.本文使用MAUSM方法分别计算了在跨声速和超声速条件下的NACA0012翼型流场,并与中心差分格式的计算结果进行比较.结果表明,对于存在激波的翼型流场,MAUSM方法是有效的.     

高速弹丸诱导斜爆轰激波结构实验研究

斜爆轰推进系统在高超声速推进领域具有广阔的应用前景,其释热迅速、比冲高、燃烧室结构简单的优点吸引研究人员的持续关注.然而,斜爆轰的地面试验同时涉及到高速试验环境模拟、燃料与氧化剂混合、高温燃烧流场结构测量等技术难点,当前国内外系统的试验研究仍然十分有限,难以支撑斜爆轰发动机的研制.为了研究自持传播的斜爆轰激波结构与波面流动特性,基于爆轰驱动二级轻气炮开展了高速弹丸诱导斜爆轰实验研究,使用直径30 mm球头圆柱形弹丸发射进入充满氢/氧可燃混合气体的实验舱中以起爆斜爆轰波,并采用两种阴影技术对实验流动结构进行测量.实验中在不同速度、不同充气压力下观察到三种弹丸诱导激波结构,即激波诱导燃烧、弹丸起爆爆轰波和相对弹丸驻定的斜爆轰波,实验舱充气压力下降则会造成爆轰横波尺度增加与波面流动失稳.实验中,斜爆轰激波角与理论分析结果吻合较好,弹丸气动不稳定带来较大的弹丸攻角会对激波角测量带来一定偏差.通过对斜爆轰波波面法向传播速度的测量发现,随着远离弹丸,斜爆轰传播速度由弹丸飞行速度衰减至接近实验气体CJ速度,弹丸速度的降低会加速斜爆轰波传播速度的衰减.     

环形激波绕射、反射和聚焦流场的间断有限元模拟

采用间断有限元方法对环形激波在圆柱形激波管内绕射、反射和聚焦流场进行了数值模拟。将二维守恒方程的间断有限元方法发展到轴对称Euler方程,并对环形激波绕后台阶流动进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉运动激波在圆柱形激波管内传播的复杂流场结构;在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力,在聚焦点附近不会产生振荡或抹平间断现象。     

磁场对激波冲击R22重气柱作用过程影响的数值模拟

为研究平面入射激波与磁化R22重质圆形气柱的作用过程,首先通过数值方法得到了不同初始条件下激波诱导R22气柱的Kelvin-Helmholtz (KH)及Richtmyer-Meshkov (RM)不稳定性导致的重气柱变形过程,并详细讨论了不同情况下透射激波在气柱内聚焦诱导射流的过程;然后在加入磁场的情况下,采用CTU+CT算法进行数值模拟,以保证数值结果满足任意时刻磁场的散度为零。计算结果表明:磁场对激波诱导R22气柱不稳定性具有抑制作用;法向磁场和流向磁场都可以很好地抑制RM不稳定性;对于KH不稳定性,法向磁场的控制效果更好,不仅可以抑制界面上涡串的卷起,还可以阻止主涡的发展,而流向磁场做不到后者;磁场对射流影响不大,射流处的磁能量可以一定程度上抑制射流的衰减,同时法向磁场可以减小聚焦时压力及速度峰值。     

圆球诱发斜爆轰波的数值研究

斜爆轰发动机是飞行器在高马赫数飞行条件下的一种新型发动机,具有结构简单、成本低和比冲高等优点.但是斜爆轰发动机的来流马赫数范围广,来流条件复杂,为实现斜爆轰波的迅速、可靠引发,采用钝头体来诱发.利用Euler方程和氢氧基元反应模型,对超声速氢气/空气混合气体中圆球诱导的斜爆轰流场进行了数值研究.不同于楔面诱发的斜爆轰波,球体首先会在驻点附近诱发正激波/爆轰波,然后在稀疏波作用下发展为斜激波/爆轰波.模拟结果显示,经过钝头体压缩的预混气体达到自燃温度后,会出现两种流场:当马赫数较低时,由于稀疏波的影响,燃烧熄灭,钝头体下游不会出现燃烧情况;而当马赫数较高时,燃烧阵面能传到下游.分析表明,当钝头体的尺度较小时,驻点附近的能量不足以诱发爆轰波,只会形成明显的燃烧带与激波非耦合结构;当钝头体的尺度较大时,流场中不会出现燃烧带与激波的非耦合现象,且这一特征与马赫数无关.通过调整球体直径,获得了激波和燃烧带部分耦合的燃烧流场结构,这一流场结构在楔面诱发的斜爆轰波中并不存在,说明稀疏波与爆轰波面的相互作用是决定圆球诱发斜爆轰波的关键.     

基于ALE方程的动网格膛口流场数值研究

基于ALE(Arbitrary Lagrangian-Eulerian Equation)方程的有限体积法,采用高精度Roe方法及结构化动网格,利用嵌入网格技术及动边界条件,对弹丸由膛内高压气体推动射出到完全飞高初始流场的整个过程进行了数值模拟.根据数值结果绘制了膛口流场密度和压力分布的时序图,形象地再现了膛口流场中初始激波、弓形激波及膛口冲击波,以及接触闻断、漩涡和剪切层等的动力学发展过程.结果基本反映了膛口流场的典型变化特征,为进一步研究真实膛口流场(如带化学反应、湍流等)奠定了基础.     

MESHLESS METHOD FOR NUMERICAL SIMULATION OF SHOCK-INDUCED COMBUSTION

发展了基于无网格方法的激波诱导燃烧流场数值模拟算法. 该算法采用二维多组分Euler方程,在点云离散的基础上采用曲面逼近计算空间导数,引入多组分HLLC （Harten-Lax-van Leer-contact） 格式计算无黏通量,运用四阶Runge-Kutta 法进行时间显式推进,化学动力学采用有限速率反应模型. 对不同预混气体中的激波诱导燃烧流场进行了数值模拟,结果同相关文献吻合较好,验证了算法的正确性.