针栓式喷注单元膜束撞击雾化混合过程数值模拟

为了全面认识针栓式喷注器喷雾场结构，基于自适应网格加密技术和分三相计算的PLIC VOF（Piecewise Linear Interface Calculation Volume of Fluid）方法对针栓式喷注单元膜束撞击雾化混合过程进行了仿真分析，通过对两路推进剂分别进行界面追踪，获得了膜束撞击雾化混合过程的详细结构特征，与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。以此为基础对膜束撞击的喷雾场结构、撞击变形过程、流场涡结构、雾化破碎典型特征及破碎后的雾化混合分布特征进行了识别分析，结果表明：膜束撞击形成了液束未穿透液膜和液束穿透液膜2种不同的喷雾扇结构。膜束撞击形成的喷雾扇呈"Ω"形，膜束同时发生弯曲变形和横截面变形。另外，膜束撞击同时受到正压和剪切应力作用，导致了一系列复杂涡流现象，使得相互作用增强，雾化混合均增强，这也是膜束撞击喷注构型优于膜膜撞击的本质原因。最后，还发现膜束撞击喷雾场液滴分布呈现分区结构特征，分别是液束控制主导的上雾化区、液膜控制主导的下雾化区及夹在中间的混合区，实际中应兼顾雾化特性和混合特性，选取中等动量比膜束撞击，这可为针栓式喷注器的理论研究和工程设计提供重要参考。     

Transient simulation of regression rate on thrust regulation process in hybrid rocket motor

The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB(hydroxyl-terminated polybutadiene) hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.     

基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

动基座近舰面流场数值模拟

航母尾迹流场对舰载机的着舰有较大影响,所以需要对其流场特点进行研究,分析不同状态下舰载机气动特性的变化。采用嵌套网格技术对航母处于垂荡状态下无人机的着舰进行了模拟。首先,利用SFS2舰船进行数值计算,验证了舰船流场的数值模拟方法。然后,对比了单独无人机定常与非定常计算结果,表明所建立的嵌套网格适用于无人机流场的模拟。接着,对航母单相流和两相流的流场结果进行了分析,结果显示甲板下方的流动对甲板上方流场没有大的影响。因此,忽略了水的影响只对航母在空气流场中的特性进行研究,结果表明航母尾迹非定常特性明显,静止航母下滑轨迹上的速度均处于周期性波动状态,且波动幅值随着远离航母而逐渐衰减;而在垂荡情况下航母尾迹变得更加紊乱,水平方向速度波动的周期性减弱,但垂向速度的波动幅值进一步增大。对于静止航母,无人机在不同时刻着舰气动特性的变化也存在差异;当航母处于垂荡状态时,无人机的升力和俯仰力矩在短时间内会有更大的波动。     

涡轮平面叶栅非轴对称端壁优化设计

开发了一套造型灵活直观、网格生成速度快的涡轮平面叶栅非轴对称端壁优化设计工具,该工具的核心技术是非均匀有理B样条(NURBS)曲面造型和网格变形.在此基础上以商业软件Isight为优化驱动器,以CFX为求解器,搭建了非轴对称端壁优化设计流程.以Pack B涡轮平面叶栅为例,对其进行了非轴对称端壁优化设计.优化后涡轮平面叶栅总压损失系数减小了12.96%.结果表明:涡轮平面叶栅端部的静压分布改变削弱了涡轮平面叶栅通道中马蹄涡、通道涡的强度,提高了涡轮平面叶栅的气动性能.      

基于非结构网格的燃面推进算法

根据惠更斯波传播原理,构建了一种在通用CFD软件的非结构网格系统上直接计算燃面推进的数值方法。该方法将燃面对应的计算网格面作为一系列平面波源,利用网格几何关系对边界节点奇异性进行判定,通过所有平面波源影响区域的叠加求解新的燃面位置。数值实现过程中,根据节点奇异性和波源影响区域以三维矢量运算得到节点的推进矢量。结合网格的重分,在通用CFD软件环境中实现了燃面的推进,所用数值算法具有良好的精度。另外,对三维动态内弹道模拟中燃面计算的稳定性和精度等问题进行了研究。     

流固耦合数据交换的插值精度影响因素研究

提出了一种投影-形函数插值的新方法,用发展了三维线性插值方法。又包括使用优化插值和常体积转化等6种插值方法,从表征多种不同网格质量和待插值函数值情况的角度,进行了插值精度验证的应用研究,探索了网格质量和待插值函数值对插值精度的影响程度。研究结果表明,投影-形函数插值方法的效果很好;流/固的网格质量及其匹配状况在很多情况下对插值精度的影响要远大于插值方法本身。     

星载大型反射面天线的刚-柔-姿控一体化在轨振动分析方法

随着大型可展开环形天线在航天器上的应用,口径越来越大,指标更加严格,卫星姿控、轨控、太阳翼驱动等导致的机械运动必然会引起大型反射面天线的振动,从而造成电性能降低,影响任务完成质量。提供了一种获取大型环形天线在轨振动影响的刚–柔–姿控一体化分析方法,建立了集扰动源、整星刚柔耦合动力学模型、姿态控制系统、天线振动影响分析的一体化仿真分析模型,实现了在典型扰动模式下的环形天线的振动响应计算、环形天线整体指向和变形计算。分析结果为天线在轨振动影响分析、性能指标预示、振动传递机理及抑制措施提供支持。