激波与火焰相互作用发展过程的数值模拟

摘要

本文主要研究激波与火焰相互作用过程中爆轰的产生问题,课题源于国家自然科学基金青年科学基金项目(10802077/A020403)——激波与火焰相互作用与爆燃转爆轰强化机制的数值模拟。本文分别对此进行了理论研究,模型建立,程序编制和数值模拟。
　　 理论研究部分介绍了对激波与火焰相互作用起重要作用的三种不稳定性现象和成因,Richtmyer-Meshkov不稳定性、Kelvin-Helmholtz不稳定性和Rayleigh-Taylor不稳定性,尤其是R-M不稳定性。分析了斜压效应的产生以及它对涡量和R-M不稳定性的发展影响,推导说明了K-H和R-T联合不稳定性对流体界面失稳的影响,还介绍了界面不稳定性的五个发展阶段,着重讨论了界面不稳定的非线性发展阶段。
　　 计算模型是常温常压下,在一个无限长的激波管中,充满一定化学当量比的氢气和空气可燃混合气,在管道中轴线上有一个正在燃烧的球形火焰,管道左侧有一入射激波向右传播,当到达火焰面时,激波开始与火焰相互作用。文章还采用了在流场中加入小火焰和障碍物的计算模型。
　　 本文采用二维非定常带化学反应的Euler方程作为控制方程,详细氢氧9组分20步基元反应为化学反应模型,时空守恒方法为计算方法。利用FORTRAN语言编制计算程序,并在曙光PHPC100型高性能计算机上运行程序。
　　 通过对结果的研究,找出了爆轰形成关键参数:激波马赫数和火焰半径。马赫数为1.4不会发生爆轰现象;马赫数为2.0时,半径大的火焰可以产生爆轰。马赫数为2.3时,所有半径的火焰都可以产生爆轰。火焰半径决定了反向透射激波的形成和传播方向。能量聚集区和激波多重诱导区最容易产生爆轰。流场加入小火焰和障碍物增大了两种区域的产生机会,缩短爆轰形成的时间和空间跨度。