界面捕捉中耦合Level-set与VOF算法

针对传统的level-set方法在处理界面问题时守恒性差的缺点,提出了一种耦合level-set方法和VOF方法的算法(CLSVOF)。不仅保持了level-set界面光滑性好的特点,还结合了VOF高守恒性的优点。Level-set方程的求解采用五阶WENO格式,而VOF的体积分数输运方程采用Gueyffier的界面重构方法求解。三维数值计算结果表明这种耦合算法对level-set的守恒性有明显的改善。     

超声速燃烧凹腔质量交换特性的混合RANS/LES模拟

对超声速冷流及燃烧流条件下超燃冲压发动机火焰稳定凹腔的质量交换特性进行了研究.采用混合RANS/LES(Reynolds-averaged Navier-Stokes/large-eddy simulation)方法对非定常流场进行数值模拟,系统研究了凹腔结构参数和横向喷流燃烧放热对凹腔质量交换特性的影响.首次采用"原子追踪法"对化学反应流条件下的凹腔质量交换过程进行了刻画.结果表明:在冷流条件和燃烧流条件下,凹腔驻留时间均随长深比的增大而增加;凹腔后缘倾角对驻留时间的影响在冷流条件下比较显著,而在燃烧流条件下明显减弱;相同结构的凹腔,在燃烧流条件下的驻留时间相对于冷流条件明显减小.      

凹腔稳定的超声速火焰结构的实验研究

凹腔是高湍流度超声速燃烧的火焰稳定装置。为研究其尺寸结构及喷注位置对不同燃料燃烧性能的影响,设计了不同凹腔、不同位置横向喷注的氢及乙烯的超燃实验,拍摄了凹腔内OH.基与CH.基自发辐射火焰图像。并利用非线性分形的研究方法,分析了火焰分形维数与燃烧程度及火焰稳定性的关系。用大涡模拟方法对二维凹腔氢超燃流场进行了数值仿真,分析了凹腔流场结构,并对凹腔后缘压力时间序列进行了相空间重构,计算了其分形维数。实验及仿真结果均表明:燃烧的分形维数高时,燃料与空气接触充分,燃烧剧烈,易于火焰稳定。     

超声速来流稳焰凹腔上游气体燃料横向喷注的流动混合特征

对Ma=1.7的空气超声速来流稳焰凹腔上游气体燃料横向喷注的流动混合特征进行了实验和数值模拟研究。在喷注的模拟燃料氦气中添加丙酮,采用平面激光诱导荧光技术(PLIF)观测了设置开式凹腔时燃料在展向和流向的空间分布。针对下游布置不同凹腔时燃料分布相似的实验结果,选取L/D=7的凹腔对其燃料混合流场进行了大涡模拟。实验和仿真结果表明:实验来流条件下喷注燃料大部分直接由喷流引起的上升反转旋涡对带走,只有小部分燃料经由喷流与凹腔流动涡结构的相互作用进入凹腔剪切层,并随剪切层运动对流进入凹腔内部。     

二维凹腔超声速流动的混合RANS/LES模拟

基于Menter的k-ω SST湍流模型构建了一种混合RANS/LES模拟方法,通过采用一个与到壁面距离相关的衔接函数将处理近壁区的SST模型过渡到处理主流区的Yoshizawa一方程亚格子模型.利用此方法对用于超燃冲压发动机的凹腔的二维超声速流动进行模拟,模拟结果再现了二维凹腔剪切层的拟序结构,计算得到的凹腔自激振荡频率、时均统计的压力分布以及压力剖面与实验结果吻合较好.     

平面激光诱导荧光技术在超声速燃烧火焰结构可视化中的应用

平面激光诱导荧光（PLIF）技术能够高时空分辨成像火焰结构并用于研究超声速燃烧机理。利用OH-PLIF与CH-PLIF技术研究了超声速燃烧的火焰结构。其中，利用OH-PLIF技术对燃烧室中3个展向截面与2个流向截面的凹腔稳定火焰反应区结构进行成像，利用CH-PLIF技术观测凹腔火焰放热区结构。实验结果表明:全局当量比较低时燃烧主要发生在凹腔中，OH沿中轴线对称分布；高当量比时火焰位置更高，OH主要沿燃烧室两侧壁面分布；CH所存在的超声速燃烧放热区呈现高度褶皱和破碎结构，放热区分布在比反应区更窄的区域。     

复杂湍流流动的混合RANS/LES方法研究

混合RANS/LES方法是近年来复杂湍流模拟的重要方法。简要回顾了混合RANS/LES模拟方法的发展历程,着重总结和分析了分离涡混合模型、类Menter SST混合模型以及湍流能量谱一致混合方法三类混合模型的构造方法和发展历程,对这三种方法的应用以及优缺点进行了简要评述,最后指出了混合RANS/LES方法发展应该兼顾的问题,为下一步的混合RANS/LES模拟提供了参考。     

超声速气流中凹腔主动喷注的强迫点火过程实验研究

为了研究在总温846K,总压0.7MPa,入口来流马赫数2.1的超声速来流条件下凹腔主动喷注对点火过程的影响,利用高速摄影相机进行观测,对比了凹腔前壁喷注、后壁面喷注等组合喷注方式下开展的发动机乙烯点火试验。基于对高速摄影图像处理的统计分析,研究发现在点火工况条件下,凹腔后壁面喷注相比于前壁喷注更容易使整个凹腔内形成稳定的火焰,当凹腔主动喷注当量比达到0.03时就能使凹腔内形成稳定的火焰;在凹腔后壁面喷注的条件下,当凹腔主动喷注当量比达到0.06,全局当量比达到0.17时,火焰就会穿过剪切层引燃凹腔下游横向射流,在整个发动机中稳定燃烧;凹腔前壁喷注和后壁面喷注相结合的喷注方式能进一步促进初始火核的形成与传播,当凹腔主动喷注当量比达到0.05,全局当量比达到0.16时,火焰就能穿过剪切层引燃凹腔下游横向射流,在整个发动机中稳定燃烧。     

脉冲等离子体推力器推进剂烧蚀传热计算

运用磁流体动力学的方法,建立脉冲等离子体推力器工作过程模型,展开三维仿真研究,对推进剂的烧蚀传热过程进行分析。通过该模型计算得到的推进剂烧蚀质量和元冲量能够反映出实验结果的大体趋势,同时针对不同电容以及两种不同推力室构型的情况进行了评估。结果显示,推进剂的温度变化波形反映着放电波形的变化波形,影响着推力室等离子体的流动,良好的推进剂表面温度分布更有助于推进剂的喷射过程。     

旋转爆震波发动机二维数值模拟

采用一种改进的化学非平衡流解耦方法对Euler反应流方程解耦处理,对流项采用五阶WENO格式离散,化学反应源项的刚性采用简化的隐式方法处理,时间步进采用二阶精度的Runge-Kutta方法,对H2/Air预混气旋转爆震发动机内流场进行了二维数值模拟。模拟结果给出了不同发动机尺寸条件下的详细起爆过程,结果表明当发动机尺寸小于临界直径时无法成功起爆;详细分析了流场结构和爆震波形状,旋转爆震波的传播速度与理论预测值吻合;性能分析结果表明在喷注总压低于燃烧室平均压力时仍可实现推进剂喷注,没装尾喷管的情况下发动机比冲达到176.5s。