可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

SiC/CVD SiC复合涂层的抗氧化及抗热震性能研究

为了提高C/SiC复合材料耐高温性能,采用泥浆浸渍裂解与真空化学气相沉积(CVD)在材料表面制备了SiC/CVD SiC复合涂层,通过XRD、SEM分析了涂层组成与结构;研究了复合涂层的高温抗氧化(700～1 500℃)和抗热震性能。结果表明,泥浆浸渍法制备的SiC涂层具有一定的封孔效果,可使材料开孔率下降,但高温抗氧化效果并不佳,1 200℃氧化10 min后材料弯曲强度保留率下降明显仅有86%。CVD SiC涂层结构致密,与SiC封孔涂层结合较好,在700～1 500℃具有较好的抗氧化效果,随着氧化温度的升高,氧化后涂层完好,表面O元素逐渐增加,材料失重率缓慢增加但不大于0.5%,且材料性能并未下降。涂层材料在1 200℃-10 min短时热震5次后材料弯曲强度保留率仍有95%以上,且未出现开裂、剥落等热震损伤。在1 200℃-30 min长时热震10次后,涂层材料基本被完全氧化,材料失去保护作用,弯曲强度下降至90%左右。     

基于测试数据的卡箍非线性等效建模方法

为了快速准确地计算卡箍管路系统的振动响应，提出了基于测试数据的卡箍非线性等效建模方法。对卡箍直管系统进行低幅值激励的模态测试，利用试验数据对有限元模型进行修正，识别卡箍线性刚度、线性阻尼参数。对卡箍直管系统开展不同激励水平下的恒力测试，通过样条插值多项式法将恒力频响函数转换为恒位移、恒速度频响函数，并进一步识别卡箍非线性等效刚度、非线性等效阻尼，建立了考虑卡箍非线性的等效动力学模型。研究表明:该卡箍等效动力学模型的预测频响函数与试验频响函数的最大误差仅为-6.9%，可以有效地反映卡箍实际的刚度、阻尼特性，对解决卡箍管路系统的振动响应分析具有重要的指导意义。      

短螺旋型燃烧室旋流结构及气动边界发展规律数值模拟研究

短螺旋型燃烧室的头部轴线沿发动机周向与发动机轴线偏转一定角度,可以在保证燃烧效率不变的条件下有效缩短燃烧室轴向长度。短螺旋型燃烧室流场的最大特点是旋流流动单侧受限。为了研究头部安装角α变化对燃烧室旋流流动特性的影响,基于数值方法对短螺旋型燃烧室进行了计算分析。结果表明:随着α变化,旋流器下游旋涡依次出现对称环状、马蹄状、环状结构;随着α增加,气动边界逐渐出现并抑制旋流的周向扩张,导致流场出现不同的旋涡结构;不同α下的切向角动量随轴向距离增加而衰减,但α为35°和45°时,气动边界在非受限侧出现并对旋流产生约束,角动量衰减变慢;当α为0°,15°,35°,45°时,燃油液滴依次集中分布在旋流器下游两旋涡边缘、侧壁面和头部端壁、非受限旋涡边缘。本文研究了不同α下的旋涡结构及气动边界沿轴向的演化过程,为短螺旋型燃烧室进一步的设计与优化提供基础。     

Unsteady wakes-secondary flow interactions in a high-lift low-pressure turbine cascade

Detailed experimental measurements were conducted to study the interactions between incoming wakes and endwall secondary flow in a high-lift Low-Pressure Turbine (LPT) cascade. All of the measurements were conducted in both the presence and absence of incoming wakes, and numerical analysis was performed to elucidate the flow mechanism. With increasing Reynolds number, the influence of the incoming wakes on suppressing the secondary flow gradually increased owing to the greater influence of incoming wakes on reducing the negative incidence angle at higher Reynolds numbers, leading to a lower blade loading near the leading edge and suppression of the Pressure Side (PS) leg of the horseshoe vortex. However, the effect of unsteady wakes on suppressing the profile losses gradually became weaker owing to the reduced size of the Suction Side (SS) separation bubble and increased mixing loss in the free-flow region at high Reynolds numbers. Incoming wakes clearly improved the aerodynamic performance of the low-pressure turbine cascade at low Reynolds numbers of 25,000 and 50,000. In contrast, at the high Reynolds number of 100,000, the profile loss at the midspan and mass-averaged total losses downstream of the cascade were higher in the presence of wakes than in the absence of wakes, and the unsteady wakes exerted a negative influence on the aerodynamic performance of the LPT cascade.     

基于变马赫数的五孔探针三维插值方法

探讨了五孔探针气动数据的插值方法,为提高插值精度,开发了基于传统线性插值法的三维线性插值法。该方法把同一探针在不同马赫数下的校准图形成三维数据库,将实验数据通过三维图进行插值。并使用改进前后的两种插值方法分别对校准风洞测得的数据进行整理,对比结果证明:在实验工况连续变化的情况下,三维线性插值法在插值精度上要优于传统线性插值法。特别是在来流马赫数变化较大时,该方法可以改善单一校准文件处理造成的数据精度问题,可用于自动化流场采集系统,为流场高速高精度采集奠定基础。     

基于双环滑模控制的直/气复合控制系统设计

直/气复合控制导弹具有强干扰、强非线性以及强耦合等特点，传统的姿态控制器难以适用于该种复杂干扰并存的情况，文章提出了基于双环滑模控制的直/气复合控制器。首先采用有限时间收敛趋近律分别构造内外环滑模控制器，并将角速度回路的滑模变量量化为直接力指令，以解决空气舵与姿控发动机之间的耦合问题。接着使用非线性扩张状态观测器估计综合干扰，从而设计控制器补偿侧向喷流干扰及模型不确定性等综合干扰的影响。然后基于李雅普诺夫方法证明了控制系统闭环稳定，分析了干扰补偿对控制器收敛域的影响。最后仿真结果表明，该方法跟踪速度快，动态过程平稳，具有较强的干扰抑制能力，具有很强的鲁棒性。     

Which way is the bookstore? A closer look at the judgments of relative directions task

We present a detailed analysis of a widely used assay in human spatial cognition, the judgments of relative direction (JRD) task. We conducted three experiments involving virtual navigation interspersed with the JRD task, and included confidence judgments and map drawing as additional metrics. We also present a technique for assessing the similarity of the cognitive representations underlying performance on the JRD and map-drawing tasks. Our results support the construct validity of the JRD task and its connection to allocentric representation. Additionally, we found that chance performance on the JRD task depends on the distribution of the angles of participants’ responses, rather than being constant and 90 degrees. Accordingly, we present a method for better determining chance performance.     

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

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