倾斜射流撞壁铺展的数值仿真

为了加深对射流撞壁铺展形成液膜的认识，开展倾斜射流撞壁数值仿真研究。采用网格自适应加密技术对射流撞壁后的液膜铺展过程开展两相数值仿真研究，获得并分析了典型工况下液膜铺展的过程、流场结构以及射流撞壁区局部流动特征。从数值仿真结果中能清晰地分辨出液膜的关键特征，与试验结果的对比也表明了数值仿真方法的可行性与准确性。通过数值仿真发现，射流撞壁后，流动以滞止点为中心，呈辐射状向四周铺展，汇入液膜边缘处水跃区后流动方向偏转，并继续向下游流动，这是射流撞壁铺展形成液膜的基本过程；液膜的惯性力驱动着液膜呈辐射状向外铺展，而在液膜边缘位置处的表面张力和壁面接触角的影响下，液膜边缘形成高压区推动液膜收缩，液膜惯性力在壁面的剪切作用下逐渐减小，直至减小到与液膜边缘处表面张力等其他作用力相平衡，从而确定液膜的边界；数值仿真结果也验证了撞壁区流动的滞止点处于射流与壁面的椭圆形接触面的一个焦点附近。     

气液比对气液同轴离心式喷嘴缩进室内流动过程的影响

为了研究气液同轴离心式喷嘴缩进室内部非定常流动过程,采用Level Set和VOF相耦合的方法结合网格自适应技术对缩进长度为8 mm的液体中心式气液同轴离心式喷嘴流动过程进行了数值仿真研究,计算得到了较为精细的液膜一次破碎过程、流场结构和压力振荡特性。结果表明：液膜的破碎模式受气液比的影响较大,随气液比的增加,液膜破碎模式由曲张表面波主导的破碎变为穿孔破碎。此外,清晰获得了自激振荡过程,分析了缩进室内部压力场及速度场分布特征,研究发现随着气体压降增加,气体环缝出口会出现膨胀波和激波,形成一个“扇环形”的超声速流场区域；激波后气流分离,出现旋涡,形成局部高压区,旋涡中心随激波面周期性地上下移动,致使局部压力出现周期性振荡。     

气液两相流混合模型代数重构方法

为提高国家数值风洞工程气液两相流软件混合模型求解器的自由界面模拟分辨率，增强多介质流体界面问题的模拟能力，基于代数重构方法，应用双曲正切函数界面捕捉格式对体积分数进行了插值重构，实现两相流混合模型与自由界面代数重构方法的耦合，抑制原混合模型下自由界面的非物理耗散现象，保持界面清晰尖锐；引入预处理技术，解决低速流动中由于特征值量级差距导致的计算刚性问题，提高收敛速度；应用双时间进行非定常时间推进。最后，采用方块对流、Zalesak盘刚性旋转和溃坝算例进行考核验证，数值结果显示，耦合了双曲正切函数界面捕捉格式的两相流混合模型在保持质量守恒的同时，获得了更高的自由界面分辨率，以及更丰富的多相流动细节。     

充液挠性航天器姿态机动终端滑模控制

对存在未知外界干扰、参数不确定问题的刚–液–柔多体耦合航天器姿态控制进行了研究。将液体燃料的晃动等效为球摆模型，挠性附件假设为欧拉–伯努利梁，运用拉格朗日方法建立航天器的动力学方程。将外界干扰、航天器转动惯量的参数不确定性以及液体晃动和挠性附件振动带来的耦合干扰归结为集总干扰，设计干扰观测器对其进行补偿；在干扰观测器的基础上，设计一种模糊滑模控制律。在原有的终端滑模控制基础上采用模糊控制对切换增益进行改进，达到抑制系统抖动的目的。数值仿真结果表明：所设计的模糊终端滑模控制律不仅能够实现充液挠性航天器的姿态机动，而且能够有效抑制液体晃动和挠性附件的振动，具有更好的控制性能。