流场可压缩性对涡相互作用影响的数值研究

涡相互作用作为冲压发动机喷注装置的典型抽象流动现象,研究其可压缩性影响对于认识包含化学反应的真实燃料喷射场具有一定基础理论价值。基于经过数值验证的可压缩Navier-Stokes算法与压力泊松方程初始条件设置方法研究流场可压缩性对涡相互作用演化过程的影响。结果表明,以较高涡旋马赫数表征的可压缩性在改变涡对形态的同时,具有延缓涡旋相互靠近,迟滞融合进程的作用。关于涡对系统开始融合的临界条件,可压缩相互作用开始的临界展弦比与无量纲时间相比于低速涡对明显提高。为在无量纲时间意义下统一不同马赫数涡对相互作用的进程,从涡心密度随时间变化规律出发在不可压涡对特征时间的基础上,初步构建了考虑可压缩性的时间尺度修正关系。     

采用基于火焰面的燃烧模型研究部分预混燃烧

考虑到火焰面模型的优点,采用基于RANS的稳态火焰面模型,稳态火焰面/反应进度变量模型和非稳态火焰面/反应进度变量模型对部分预混燃烧室进行了数值计算,并将这三种燃烧模型的计算结果和实验值进行对比研究。发现稳态和非稳态火焰面/反应进度变量模型均成功地预测到了部分预混燃烧中的三岔火焰结构和火焰抬举现象,分析了三岔火焰的形成机理及抬举高度。文中稳态火焰面模型计算部分预混燃烧完全失效,无法捕捉到火焰抬举现象,稳态火焰面/反应进度变量方法计算的火焰抬举高度仅仅为10,基于反应进度变量为水、二氧化碳、一氧化碳和氢气的质量分数之和的非稳态火焰面/反应进度变量方法计算的火焰抬举高度仅为20,和实验值之间的误差为42.8%,而基于反应进度变量为一氧化碳和二氧化碳质量分数之和的非稳态火焰面/反应进度变量方法计算的火焰抬举高度大致32,和实验测量值35非常接近,误差仅为8%。对燃烧热力学标量温度和组分的计算,可以发现非稳态火焰面/反应进度变量方法和实验结果吻合最好,其次是稳态火焰面/反应进度变量模型,最差的是稳态火焰面模型。     

高速复杂流动PIV技术研究实践与挑战

粒子图像测速技术目前已经发展成为实验流体力学领域应用最广泛的非接触激光测试方法之一,为认知复杂流动机理提供直观的流场信息.本文基于超声速流场PIV技术研究实践,针对示踪粒子布撒器设计、粒子松弛特性模型构建、激波流场测试分析、超声速平板湍流边界层结构分析等方面具体问题的研究和认识,从理论、定量化的角度深入分析了应用于超声速流场PIV技术现阶段依然存在的问题.从应用于超声速流场PIV技术的原理出发,针对高速复杂流场的PIV测试现状,总结了应用于超声速流场PIV技术发展过程中的光学部件、示踪粒子及布撒系统所遇到的一系列挑战,以及国内外利用PIV技术在高速复杂流场研究中所取得的成就,针对PIV技术能否适用于高超声速流场的测量做了系统化地探索.并根据实践经验提出了应用于超声速流场PIV技术未来的发展方向:通用的精确的PIV方法不存在,必须从具体研究的流动机理角度改造相应的PIV测试手段.     

亚声速横向气流中液体射流破碎过程的直接模拟

为研究亚声速横向气流中液体射流柱的变形、弯曲以及其破碎过程,利用LES结合VOF的方法,对射流破碎过程进行了直接模拟。通过计算观察得到射流柱进入到横流气体中后由于RayleighTaylor(RT)不稳定性和Kelvin-Helmholtz(K-H)不稳定性的共同作用迅速发生变形,并形成表面波,同时伴随着细小液滴的脱落,气流在射流柱后方形成涡状结构。在射流柱的迎风面上同时存在两种表面波,即RT表面波和KH表面波。模拟所得液柱纵向初始表面波波长为0.22mm,与K-H表面波波长理论公式计算所得相接近。通过对比射流柱的背风面和迎风面两侧,发现在横向来流中背风面上的大量细小液滴主要是由横向来流对液柱的不稳定性扰动造成迎风面上的液体向背风面方向挤压并剪切,最终脱落造成。     

姿控发动机喷管内流场及推力计算研究

运用时间相关法数值求解定常Ｎ－Ｓ方程和采用有限体积方法离散控制方程,对姿控发动机喷管内流场及其产生的推力进行了计算,得到了较为精细的流场图象。计算结果表明,不同壁面温度条件可以显著改变喷管壁面附近的流场结构及其产生的推力。     

基于IDDES框架的γ-Reθ转捩模型

针对高超声速流动的γ-Re_θ转捩模型在模拟强制转捩时存在捕获边界层内扰动不足的缺点，将RANS（Reynolds Averaged Navier-Stokes）框架改为IDDES（Improved Delayed Detached Eddy Simulation）框架，既能像基于RANS框架的转捩模型一样模拟复杂构型的自然转捩，又能发挥IDDES能够捕获更多脉动信息的优点，较为准确地模拟粗糙颗粒诱导强制转捩。通过对一系列简单构型的自然转捩及来流马赫数为6条件下平板上单个粗糙颗粒诱导强制转捩的模拟表明，模型既能体现γ-Re_θ转捩模型的优点，在自然转捩模拟中具有较强的鲁棒性，能够反映雷诺数等因素对转捩位置的影响规律；也能体现IDDES方法的优点，能够捕捉粗糙颗粒诱导的扰动及涡结构，从而较为准确地刻画出强制转捩的整个流程。     

直流喷射首次破碎的形变过程研究

为深入了解直流喷射过程中射流柱的细部结构以及其脱落过程,利用LES结合VOF的方法,对静止大气中垂直射流现象进行模拟。模拟得到的液相喷雾结构和试验结果能够很好地吻合。通过计算观察得到射流柱进入到大气中后由于Rayleigh-Taylor(RT)不稳定性迅速形成伞状的头部,头部的边缘在不断变薄失稳的过程中脱落形成液带,液带脱落产生的液滴具有等距性,间距为0.106mm,并在伞状边缘的下方形成一个气涡,气涡与上游的射流柱相互作用,促进射流柱表面的脱落。射流柱表面由于Kelvin-Helmholtz(K-H)表面波的作用呈现鱼鳞状的结构,并导致射流柱整体断裂、破碎,其表面波波长由初始的0.26mm迅速增长到0.78mm。     

Scramjet燃烧室流场的三维并行数值模拟及试验比较

作者采用时间相关法，通过有限体积离散，运用带化学反应的全N－S方程，在神州巨型机上，针对试验模型，对油气比Φ＝0．0和0．35的喷氢Scramjet燃烧室流场进行了三维并行数值模拟，得到了流场的精细结构。并行模拟所得壁面压力分布与试验所测得的壁面压力分布吻合较好。     

合并式电离—复合反应的DSMC模拟

采用ＤＳＭＣ方法研究分子碰撞反应机理,在分析水平上模拟绝热封闭体系事并式电离－复合反应的驰豫过程。     

Scramjet燃烧室流场的三维并行数值模拟及试验比较

作者采用时间相关法,通过有限体积离散,运用带化学反应的全N S方程,在神州巨型机上,针对试验模型,对油气比Φ=0.0和0.35的喷氢Scramjet燃烧室流场进行了三维并行数值模拟,得到了流场的精细结构。并行模拟所得壁面压力分布与试验所测得的壁面压力分布吻合较好。