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

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

发射扰动与初始弹道耦合模型研究

针对火箭、导弹发射扰动与初始弹道互相耦合引起的弹道散布问题,提出发射扰动与弹道解算相耦合的计算分析模型。该模型以多体系统动力学为基础,建立能够模拟弹架相互作用和弹体初始扰动的发射动力学模型,并将弹体受到的气动载荷转化到弹体坐标系下进行刚体动力学计算以获得弹道参数。通过滚转弹应用实例分析表明,采用此模型能够有效模拟发射扰动与初始弹道相互耦合状态;弹架间隙扰动与气动载荷作用都会对弹体在飞行时的姿态角及飞行位置产生较大影响。当存在1 mm的弹架间隙且有气动载荷作用的影响下,与无弹架间隙和气动载荷的作用影响的结果对比发现,存在弹架间隙扰动的影响会使得弹体在飞行过程中的俯仰角和弹道倾角的幅值范围减小4°左右,也使得弹体在飞行过程中的Y向位移量在1.5 s时刻减小6 m左右;存在气动载荷作用的影响,会使得弹体在0.5 s撤去推力后的姿态角成波动式变化,滚转弹稳定飞行,也会使得弹体Y向位移量在撤去推力后持续的平稳增加。     

前体尾流对降落伞工作性能的影响

为研究前体尾流对降落伞工作性能的非定常影响，基于Realizablek-ε湍流模型采用PISO算法开展了物伞系统的非定常绕流数值计算，获得了精细的流场旋涡结构。在此基础上，研究了不同拖曳比下物伞系统的尾涡演变规律、流场分布规律以及伞衣气动特性变化。结果表明:前体尾涡导致伞衣入口处的涡量大小和方向时刻变化，随拖曳比增加，涡量黏性耗散增强，进入伞衣的旋涡强度逐渐减弱，伞衣入口形成稳定的负涡量区，伞衣尾涡脱离周期随之延长；拖曳比对尾涡区后端（伞衣入口处）流场压力的影响远大于前端，随拖曳比增加，流动形式逐渐由闭式转变为开式，流场的速度分布和压力分布更为对称，伞衣入口形成稳定的正压区，内外压差增加；当拖曳比大于9时，前体尾流对降落伞阻力系数和表面压强系数的影响减小。      

大型飞机平尾翼尖涡对后体涡系影响的实验研究

大型飞机在飞行过程中机身后体会产生一对反向旋转的脱体涡（后体主涡），该涡与平尾翼尖涡共同构成飞机后体的涡系结构。在风洞中，利用激光粒子测速（PIV）方法，对单独后体和加装不同展长平尾的后体，分别研究涡系结构的动力学特征。结果表明:后体主涡的涡核中心沿流向明显向上移动；加装平尾后，涡系呈现典型的四涡结构，平尾翼尖涡对后体主涡影响显著，加大了后者向上移动的趋势，同时使其沿展向外移，并显著削弱其涡旋强度；平尾展长增加后，后体主涡受到的影响有所减弱。在低速环境下，来流速度对后体涡系结构的无量纲动力学参数影响较小。     

双模块下颌式内转进气道/圆锥前体一体化布局研究

为实现三维内转进气道的内收缩流场与圆锥前体的外压缩流场的良好匹配，提出了一种双模块下颌式内转进气道/圆锥前体(Double-modules Chin Inward-turning Inlet and Conical Forebody，DCII/CF)一体化设计方法，获得一种新颖的双发并置、侧向安装的DCII/CF一体化布局。针对该布局形式，开展了DCII/CF一体化构型与传统的单模块内转进气道/圆锥前体(Single-module Inward-turning Inlet and Conical Forebody，SII/CF)一体化构型的数值对比研究。结果表明：DCII/CF一体化布局不仅为内转进气道提供了优秀的前体附面层排移效果，还有效避免了传统SII/CF布局中前体附面层与进气道内部流场之间的相互干扰。在Ma∞=6.0设计状态，DCII/CF一体化布局的进气道总压恢复系数相较传统的SII/CF布局有了显著提高，从0.403提高至0.482；但由于前体附面层的排移，该布局的捕获流量略有降低， SII/CF的流量系数为0.956，该布局则为0.917。而在非设计状态，该布局形式同样具备较好的总压恢复性能，在Ma∞=5.0与Ma∞=4.0的总压恢复系数分别达到了0.586和0.682，明显高于SII/CF的总压恢复系数0.507和0.619。     

细长旋成体大迎角绕流的背涡结构与压力脉动特性

在亚临界流动范围内,通过尖拱细长旋成体在无侧滑状态下的脉动压力测量,对大迎角下细长旋成体非对称背涡结构沿轴向的压力脉动特征进行了研究分析,揭示了压力的脉动幅度与频率沿旋成体轴向的变化规律及其与背涡结构之间的对应关系,并对压力脉动特征沿轴变化的原因进行了初步探讨.     

用亚稳态正电子素研究络合的瞬态反应（Ⅲ）

本文用自制的装有 BaF_2探头的高探测效率的正电子湮没寿命谱仪测量了在室温下正电子素 o-Ps 与 p-RPhG 及 m-RPhG(R=NO_2,Cl,H,CH_3,cH_3O)在20%(V/V)二氧六环水溶液中的络合形成速度常数。实验结果表明,这些o-Ps-M 络合形成速度常数 K 值均与作为配位体的 p-RPhG、m-RPhG 的碱性解离常数 K 的负对数 pK_2 呈线性函数。这表明上述瞬态络合反应满足线性自由能关系,由此可推断出上述络合反应的一些重要性质。     

民机后体试验技术研究

本文介绍了在低速风洞中民机后体试验的方法和试验技术。试验结果表明本项技术可以满足民机型号研制需求。     

求解含多个右端向量非对称线性方程组的自适应块拟最小残量IOM(q)算法

许多实际应用问题需要求解含多个右端向量的大型非对称线性方程组 ,通常是把原来方程组分成单独几个含一个右端向量的方程组 ,再用某种迭代法分别单个求解 ,而更加经济有效的方法是应用能同时产生几个迭代向量的块迭代法来直接求解。本文在 IOM(q)算法的基础上 ,提出一种求解此类方程组的块拟最小残量 IOM(q)算法 ,讨论了如何收缩掉已收敛的部分方程组以及如何从产生的块 Krylov序列中删除线性相关或几乎线性相关向量的自适应技术。数值试验表明 ,此新的自适应块算法比块 GMRES算法及其他相关算法具有更好的收敛行为、更少的计算量和 CPU计算时间 ,是求解此类方程组的一种更加经济有效的算法。