喷水射流预冷对加力燃烧室特性影响分析

为了探究采用射流预冷技术之后加力燃烧室性能，开展了不同喷嘴布置方案、喷水量和来流温度对预冷效果的影响研究。对射流预冷发动机工作过程进行了简化，建立了加力燃烧室进口前段射流预冷喷水特性计算的数学模型。同时搭建了小型试验台，通过与试验结果的比对验证了该模型的准确性，并利用该模型对射流预冷效果进行了仿真预测。结果表明：提高喷嘴数量与布置均匀性能够小幅度改善预冷效果；当来流温度不变时，射流预冷喷射腔室出口处的液态水蒸发量随着喷水量的增加而提高，但蒸发率却处于下降的趋势；当喷水量达到2%时，加力燃烧室燃烧效率对比不喷水工况会有一定的提升；喷水量达到4%以后，加力燃烧室出口温度及燃烧效率随着喷水量的提高而降低；喷水量大于8%以后，恶化了加力燃烧室（V型火焰稳定器）贫油熄火极限与燃烧效率；喷水量达到最大10%时，油气比需从原来设计工况的0.052上升到0.064才能保持稳定点火且对比不喷水时工况，加力燃烧室出口温度由1860K下降到1373K，燃烧效率由80.2%下降到69.2%。     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

电磁力作用下槽道流速度响应的放大机理

在电解质溶液中,电磁场产生的电磁力可以控制流体的运动,从而达到很好的减阻、增升、减振等效果。但由于所施加的电磁力较大,导致控制效率很低,因此以较小的电磁力诱导出大的流动响应成为提高流动控制效率的关键。以层流槽道流动作为研究对象,在槽道的下壁面施加沿展向余弦分布的展向电磁力,推导了线性条件下流向响应速度的解析解,并通过直接数值模拟对非线性条件下的响应进行了计算。结合解析解和数值解,揭示了流场中速度响应的放大机制,讨论了电磁力和流场参数对响应放大效果的影响。结果表明:当振幅较小时,速度响应处于线性范围内,其放大倍数与Re~2成正比,随着渗透深度的增大,先迅速增大后缓慢减小;随着展向波数Kz的增大单调减小。随着振幅的增大,放大倍数进入非线性范围,其值逐渐减小,但速度响应值先增大后减小。在振幅处于10~(-3)～10~(-2)量级时,速度响应可达到的最大值超过0.2,此时的放大倍数在102量级。因此,利用流场的放大效应,是实现高效流动控制的重要环节。     

SAS与PANS模型在圆柱绕流中的应用比较

分别运用SST(Shear Stress Transport)、SST-SAS(Scale-Adaptive Simulation)、两种变fk(模化湍动能的比例)函数的SST-PANS(Partially Averaged Navier-Stokes)湍流模型对Re=3900的圆柱绕流进行了数值研究,重点从湍流结构捕捉、气动力计算、涡黏性控制等方面,比较了SAS与PANS两类RANS/LES混合模型的计算能力,并通过不同网格计算分析了模型的网格敏感性。数值结果表明:SAS及两种变fk方法的PANS模型均具有求解小尺度涡运动的能力,并能较好地反映出绕流尾迹的三维非定常特性,同时PANS模型能捕捉到更多的非定常结构;SAS模型中自适应尺度Lvk立足于当地流动,对网格依赖较小,计算的湍动粘度分布更合理,能够更好地计算剪切层及回流区;两种PANS模型网格独立性较差,出现了雷诺应力不足的现象;类DES可变fk函数构造相对简单,所得fk分布更准确,使用tanh函数计算的尾迹区fk值偏低,对流场调控能力稍差。     

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符；反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。      

圆弧激波缓聚点火及后续燃烧传播

采用激波管实验和准一维数值模拟的方法，对预混可燃气体中圆弧汇聚激波的自点火现象及后续燃烧波的传播特性进行研究。其中圆弧汇聚激波由平面运动激波通过精确设计的弧形过渡管段转变得到。研究表明:收缩段中圆弧汇聚激波波后的非均匀梯度环境由激波在平直段、弧形过渡段和扇形收缩段中传播所分别诱导的3个梯度区共同构成。随着圆弧汇聚激波的不断增强，圆弧激波后某处首先形成一个无激波的温和反应区。该反应区逆流锋面的初期运动速度远超Chapman-Jouguet（CJ）爆轰波速，而反应产物区流动则呈现出一定弱爆轰波特征。进一步分析发现，该反应锋面本质上是一种"自发反应波"（spontaneous reaction wave），而非常规意义上的动力学波，其速度与汇聚激波波后气流点火时间梯度的倒数吻合。而后，反应区的扩张速度很快降至CJ爆轰波速以下，伴随反应锋面附近激波的产生以及激波-火焰复合结构的形成。激波-火焰结构最终加速演变为反向传播的爆轰波。在一定的条件下，由于入射激波转变过程和汇聚所构造的特定点火环境，自发反应波可再次赶超爆轰波，成为新的燃烧波前；而当自发反应波速度再次低于CJ爆轰波速时，它将再次转变为爆轰波；在此过程中，原先的爆轰波阵面蜕变为反应产物中传播的激波。     

高超声速风洞测力数据的关联研究

通过比较分析法研究AGARD-B标模在BIA的FD-07风洞中气动力测量数据与其它高超声速风洞设备测量结果的相互关联。其目的，既可以综合鉴定FD-07风洞的流场性能和测试水平，也可以从另外意义上为数值方法的程序校验提供必要的技术支持。结果表明，在Ma=4.94-7.96范围内，所测得AGARD-B标模的力、力矩和压心位置，不仅随Ma和α的变化趋势与其它风洞一致，而且其变化量比较接近于各风洞测力数据拟合的逼近曲线。     

ADI方法求解完全跨声速非定常小扰动方程

本文修改了Whitlow的ADI格式,以便求解另一变形的完全跨声速小扰动方程。本文给出的NACA64A006和NACA0012翼型的非定常跨声速流动的计算结果与其它数值结果和实验结果的比较表明本文格式是合理的且便于实际应用的。