受平行壁限制的燃油喷雾二相流特征

通过试验诊断了燃油喷雾在模拟环形燃烧室火焰筒内、外壁面间的二相流型。其结果显示 ,与不受壁面限制的结果相比 ,其基本流型在下游有明显改变     

预混可燃流体在大长径比圆管内流动传热点燃研究

通过对预混可燃流体在大长径比圆管内的流动和传热的分析,首次提出一种用于脉冲火箭发动机多次点火的方法。新型发动机中,预混可燃流体在流动过程中被加热,并自动点火,按需要通过迅速散热保证熄火,从而实现了任意脉冲的自动点火和熄火,并可实现推力调节。对这种点火方法进行了理论分析和实验研究,成功地得到了5个工作脉冲。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

典型流阻元件的流量和熵产模型研究

为揭示流阻元件的流动损失机理,本文从热力学理论出发,建立了流阻元件的质量流量模型和熵产模型,分析了质量流量与熵产之间的关系,研究了熵产随系统压比和进口总压的变化情况,并通过预旋喷嘴和篦齿封严的实验结果对两个模型进行了分析评估。结果表明,质量流量模型与实验结果的最大偏差不超过2.8%;熵产模型与实验结果的最大偏差不超过1.9%。系统熵产随系统压比的增加而增大,随进口总压的增加而减小。当系统的进口总压,进口总温和出口静压不变时,模型中的参数a是衡量不同元件熵产大小的唯一量度。     

煤油燃料超燃冲压发动机燃烧特性实验研究

为研究空气节流时序对超燃冲压发动机点火和火焰稳定的影响,本文通过实验方法研究了13个状态的煤油燃料超燃冲压发动机的燃烧特性,煤油燃烧通过先锋氢气和节流空气增强稳定性。通过两个固定位置的压力传感器来监测火焰稳定状态,采用纹影和OH-PLIF相结合的测量手段,获得了流场结构和火焰发展信息。发动机入口来流条件为Ma = 2.0,总温950 K,总压0.82 MPa。在空气节流的作用下,煤油被先锋火焰引燃;在先锋氢撤除后,煤油仍然可以稳定燃烧。在扩张段中,空气节流和燃烧共同作用产生的激波串移动速度约为52 m/s,但在凹槽内其速度仅为3.7 m/s。通过监测点压力变化情况可以区分所研究状态的火焰稳定与否,通过对13个研究状态的考察,获得了火焰稳定临界曲线。当所研究状态点在临界曲线右上方区域时,火焰状态稳定;当所研究状态点在临界曲线左下方区域时,火焰将被吹熄;当所研究状态点在临界曲线上时,火焰不稳定,在空气节流撤除之前将被吹熄。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

不同舵翼操纵面形式的水动力及艇后流场特性数值模拟研究

为研究共翼型舵和非共翼型舵两种操纵面形式的艇后流场特性及螺旋桨推进特性,基于SST(Menter)湍流模型建立了SUBOFF标准潜艇模型尾流场数值预报模型。经试验结果验证,所建立数值模型进行潜艇尾流场及螺旋桨推进性能预报具有较高精度。将SUBOFF潜艇模型的水平舵改进为共翼型舵及非共翼型舵,对艇后流场及水动力性能进行了预报。数值结果表明：在舵角小于10°时,共翼型舵使潜艇的俯仰力矩和垂向力相对非共翼型舵提升20%以上;在舵角超过10°时,共翼型舵的水动力优势随着舵角增大而减小。尾部流场预报显示：共翼型舵在小舵角时可以有效的消除舵翼结合处的涡流,同时共翼型舵可以有效的降低桨盘面伴流的不均匀性,对尾流品质的改善效果优于非共翼型舵。螺旋桨计算结果显示：共翼型舵螺旋桨推进系数在大部分舵角下都小于非共翼型舵,在舵角为20°时,共翼型舵相对非共翼型舵推力系数下降3.5%,扭矩系数下降2.4%;同时共翼型舵的桨盘面流场均匀度要优于非共翼型舵,舵角为5°时,共翼型舵桨盘面处流场不均匀度相对于非共翼型舵要降低7.1%,舵角为25°时则降低25.1%。     

稀薄流航天器鼻锥迎风凹腔气动力和气动热性能研究

本文采用直接模拟蒙特卡罗(DSMC)方法,对高超声速稀薄流中航天器鼻锥迎风凹腔气动力与气动热性能进行了数值研究。得到了鼻锥外壁面、凹腔侧壁面以及凹腔底面的热流密度分布,分析了不同凹腔深宽比对鼻锥冷却效率以及凹腔腔体内气体参数的影响;以深宽比为1的凹腔为基准,研究了凹腔唇口钝化半径对航天器气动热与气动力的影响。数值结果表明,高超声速稀薄流中迎风凹腔能够降低鼻锥外壁面的热流密度;当凹腔深宽比达到1之后,凹腔侧壁面热流变化趋于一致,热流密度最低点的轴向位置不随深宽比改变,且凹腔底部热流很小;凹腔近底部气体均由稀薄流转化为连续流,腔内气体压力不断振荡;唇口钝化没有明显优势,虽然可以降低鼻锥峰值热流,但是会带来严重的气动力性能下降。     

基于有序BDL树集的船舶电网潮流算法及执行时间分析

针对嵌入式船舶电力模拟训练系统在频繁模拟操作指令下对潮流计算的快速性需求以及嵌入式模拟训练系统极其有限的运算资源之间的矛盾,提出一种基于有序BDL树集的船舶电网潮流算法。该算法在常规回推前推法基础上对节点编号进行优化,将多电站船舶电力系统在实际开环运行模式下等效为由多个辐射状电网构成的有序BDL树集,树集中的每个有序BDL树包含单个电站中的主配电板层、配电层和负载层,并把发电机潮流整合入主配电板层,避免了常规潮流算法中节点编号突变问题。对该算法的时间复杂度分析结果及实际运行结果进行比较分析,证明该算法比常规船舶电网潮流算法更优。     

A self-similar solution of a curved shock wave and its time-dependent force variation for a starting flat plate airfoil in supersonic flow

The problem of aeroelasticity and maneuvering of command surface and gust wing interaction involves a starting flow period which can be seen as the flow of an airfoil attaining suddenly an angle of attack. In the linear or nonlinear case, compressive Mach or shock waves are generated on the windward side and expansive Mach or rarefaction waves are generated on the leeward side. On each side, these waves are composed of an oblique steady state wave, a vertically-moving one-dimensional unsteady wave, and a secondary wave resulting from the interaction between the steady and unsteady ones. An analytical solution in the secondary wave has been obtained by Heaslet and Lomax in the linear case, and this linear solution has been borrowed to give an approximate solution by Bai and Wu for the nonlinear case. The structure of the secondary shock wave and the appearance of various force stages are two issues not yet considered in previous studies and has been studied in the present paper. A self-similar solution is obtained for the secondary shock wave, and the reason to have an initial force plateau as observed numerically is identified. Moreover, six theoretical characteristic time scales for pressure load variation are determined which explain the slope changes of the time-dependent force curve.