基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

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

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

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

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

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

具有复杂围护结构光学窗口的传热特性分析

针对具有复杂围护结构光学窗口传热特性,建立了高热流作用下的光学窗口及其复杂围护结构的辐射-导热传热模型,分别讨论了窗口的辐射半透明、光谱选择性、外部加热热流以及窗口尺寸对其瞬态温度响应的影响。结果表明,忽略窗口半透明特征会产生最大31.4%误差,忽略光谱选择性会产生最大40.4%误差;窗口外部气动加热热流影响显著,温度响应随热流的增大而急剧增大;窗口厚度对其传热特性的影响较大,随着厚度的增大窗口温度响应减小;而窗口半径产生的影响可忽略不计。     

翼尖铰接复合飞行器动力学特性研究

航空飞行器通过翼尖铰接机构复合飞行时的气动耦合效应，会造成飞机产生不同于其单独飞行时的动力学特性，出现复合飞行安全问题。为研究翼尖铰接复合飞行器的动力学特性，使用Newton-Euler方法和Robberson-Wittenburg方法建立了双机组成的翼尖铰接复合飞行器多刚体系统整体和内部的7自由度非线性动力学和运动学方程组。在气动准定常假设下建立双机复合系统非耦合气动力表达式，基于CFD方法开展复合飞行器系统的三维实体建模和非结构网格划分工作，获取复合飞行器系统的气动力数据。搭建动力学仿真平台，开展准配平方案和全配平方案下的动力学仿真。仿真结果表明:准配平方案下飞行器无法持续稳定飞行，而全配平方案下复合飞行器系统各运动参数在仿真时间内始终趋于稳定。在全配平方案下，使用小扰动假设的非解耦线性化方法重新整理7自由度动力学方程组，研究复合飞行器系统运动模态的特征值中出现的2个发散新模态，根据对应的特征向量发现2个发散模态分别由相对滚转角度和角速度主导，同时也比较分析了其他运动模态相比单机飞行时的特性变化规律。      

空间载荷机柜热控设计的仿真研究

空间站主动热控系统的作用是提供和保持航天员、仪器设备和结构部件所要求的温度、湿度和风速等热环境。本文以空间载荷机柜为热控系统的基本单元，建立了以单相流体循环回路为基础、采用液冷结合风冷的冷却方式的载荷机柜数学模型；介绍了模糊增量控制算法的工作原理；提出了采用模糊增量控制算法，以液冷支路进出口温差和风冷支路冷却量占总冷却量的比值为控制目标的反馈控制策略。同时文章采用数值模拟方法研究了控制系统的动态响应效果。仿真结果显示，模糊增量控制算法具有响应迅速、超调量小、无静态误差的特点，控制策略满足有效性、稳定性和高精确度的要求。     

R141b在圆形微通道内的沸腾换热实验研究

为提高换热强度、解决设备内部高热流密度散热问题,采用实验方法研究R141b在不同直径(D=0.5mm和1.0mm)水平圆形微通道内的沸腾换热特性,分析了热流密度(q=2.0kW/m~2～47.6kW/m~2)、质量干度(x=0～0.6)、质量流速(G=111.11kg/(m~2·s)～333.33kg/(m~2·s))的变化对平均传热系数h的影响,探究不同情况下影响沸腾换热的主导因素。实验研究表明:平均传热系数h随热流密度q的增加而减小,在不同范围内减小速率有明显差异;热流密度q=2kW/m~2～5kW/m~2时质量流速G对平均传热系数h影响较明显,热流密度较高时质量流速G对换热影响很小;在质量流速G=111.11kg/(m~2·s)～333.33kg/(m~2·s),质量干度x0.3时,平均传热系数h随质量干度x增加而明显下降,在设计微通道换热器时应尽量使R141b处于初始沸腾阶段以获得更好换热效果,并采取一定措施预防干度过高引起的换热恶化。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

气膜孔形状对高超声速飞行器对撞流气膜冷却效果的影响

文章通过仿真分析手段研究飞行高度50 km、飞行马赫数15的飞行条件下,不同孔型对对撞流的影响,得到不同孔型对气膜冷却效果的影响规律。采用计算流体动力学(CFD)方法,对在入口压力0.5 MPa、质量流量22.5 g/s的稳定短模态(SPM)工作模态下,气膜孔为圆柱直孔、收缩孔、连续扩张孔、分段扩张孔等工况开展对比研究,结果显示,扩张孔气膜冷却的壁面热流最大,圆柱孔的次之,收缩孔的最小。这表明,通过改变对撞流气膜孔的形状可以改变气流流动特性,进而产生不同的气膜冷却效果,在SPM工作模态下收缩孔的气膜冷却效果最好。