基于煤油的旋转爆轰发动机流场数值模拟

为了深入研究气液两相旋转爆轰发动机的流场结构，建立了非定常两相爆轰的Eulerian-Lagrangian模型，使用SST（shear-stress transport） k-ω模型，采用一步反应机理的化学反应模型，进行了煤油/空气非预混的二维数值模拟。结果表明，采用30 μm粒径的液滴颗粒，在来流总温1 000 K的空气流中，液滴经历雾化破碎、蒸发、混合过程，在当量比0.70~1.15范围，形成稳定单个旋转爆轰波；煤油液滴被爆轰波扫过后未完全燃烧，部分煤油组分混杂在高温产物中沿下游排出；在燃烧室入口处，爆轰波前形成的空气三角形区域面积大于液滴颗粒三角形区域。     

一种用于TSTO级间分离CFD计算的网格动态优化技术

两级入轨空天运输系统具有难度较低和经济性好的优点，级间分离是该运输系统设计的一项关键技术。相比前后串联级间分离，两级入轨上下两级间存在严重的激波相互干扰，直接影响分离物体的气动力矩。为了更好地模拟分离过程中两级之间的运动激波多次反射问题，本文结合非结构混合网格分布优化技术和压力比值激波识别技术，建立了用于级间分离运动过程的并行网格动态优化技术。基于计算流体力学软件NNW-FlowStar，完成了两级入轨运输系统并联分离过程的数值模拟。结果表明，采用网格动态优化技术，可实现级间复杂激波的精细捕捉，更精准地获得了强干扰下飞行器的气动特性，显著提高了并联分离过程中运动轨迹和姿态的预测精度。     

超声速气流中液体喷雾流动数值模拟研究

超声速气流中液体喷雾流动特性对超声速燃烧基础研究及其工程应用具有重要意义。为了定量探索超声速气流中液体横向射流雾化特性，本文对超声速气流中液体喷雾流动进行了数值模拟研究。数值模拟方法基于Eulerian-Lagrangian两相流计算架构，考虑气液双向耦合，采用KH/RT液滴二次破碎模型计算液滴雾化过程，采用大涡模拟计算气相流动。结果表明，该数值模拟方法所获得的液雾场突起结构、穿透高度、液滴平均速度分布等液雾特性均能与试验结果较符合；初始液滴直径分布对破碎后液滴平均速度影响较小而对破碎后液滴平均直径及液相平均体积分数影响较大，初始液滴直径分布需在后续的建模与模拟中进行更多研究。     

基于网格测力数据的多体分离轨迹预测方法研究

将基于风洞网格测力试验数据建立的气动力模型与刚体运动方程进行耦合求解得到多体分离轨迹-时间特性，建立了一种多体分离的离线轨迹预测方法。为了提高气动力模型的预测精度，针对移动最小二乘法（MLS）模型提出一种新的权函数形式，针对Kriging气动力数学模型通过加入样点预处理提出了Kriging-Pre数学模型。研究方法应用于来流马赫数6条件下，某并联两级入轨飞行器标模的分离特性研究。研究表明采用改进的两种气动力数学模型均可有效提高分离轨迹预测精度，得到与CFD以及风洞试验定性一致的结论。验证了本文提出的离线轨迹预测方法可以满足当前多体分离特性定性分析需求，具有较高时效性。     

激波风洞两级入轨飞行器纵向级间分离试验技术

高超声速多体分离问题是航天多体飞行器研发中的关键技术问题，基于分离过程中高速流动的复杂性，对高速多体分离的风洞试验研究极具挑战性，特别是激波风洞分离试验。激波风洞具有高速、高焓试验气流特点，更准确评估高超声速分离气动力/热特性，但是其有效试验时间短（ms量级），进行主动式动态级间分离试验极其困难。提出一种应用于激波风洞主动式多体分离试验的高速气动发射系统（HPELS），使得模型在短试验时间内完成主动分离测试，详细介绍了HPELS延迟时间、模型分离时间等精确的时间标定及时序控制方法。针对分离过程中模型的运动轨迹及气动力参数的高性能评估，发展了基于纹影图像的非接触式分离运动轨迹捕获及气动力参数测量技术。两级入轨（TSTO）飞行器的安全级间分离是典型的高速多体分离问题，设计了并联式TSTO飞行器并针对作者提出的纵向分离方案，在JF-12复现飞行条件激波风洞验证了高速动态多体分离试验技术应用的有效性，同时首次在激波风洞对TSTO纵向分离方案进行了原理性验证。初步对比结果显示，试验结果与数值计算结果具有良好的一致性。     

一种并联两级入轨飞行器纵向分离方案的数值研究

两级入轨（TSTO）飞行器或将成为下一代天地运输往返系统，其具有低成本、高效率和多用途等优点，但是两级分离成功与否将直接决定入轨任务的成败。目前的并联式TSTO飞行器多采用横向级间分离，该方法会在两级间产生复杂强气动干扰而直接增加了分离风险，所以探索一种可以避免或减弱两级强气动干扰的新分离方式是十分必要的。提出并着重分析了一种并联式TSTO纵向级间分离（LSS）方案，即轨道级在助推级背面沿着飞行方向分离，对其进行了动态分离过程的数值研究。针对新分离方案，设计了一种由宽速域乘波体和可重复使用空天飞机分别作为助推级和轨道级的TSTO组合飞行器，在高超声速条件下，采用重叠动网格技术分析了不同来流攻角（AOA）下的纵向分离流动机理、非定常壁面压力分布及气动特性变化规律。结果表明：TSTO纵向分离过程中仅存在VI型激波干扰和激波汇聚等简单的弱干扰类型，两级间无明显的激波反射或激波边界层干扰；非定常压力分布特性表明助推级前缘激波是轨道级受力变化的主要影响因素；纵向分离过程中，助推级受到的气动干扰力载荷小于轨道级。此外，不同来流攻角下，两级气动干扰流场结构具有相似性，并给出了实现安全纵向分离的攻角条...     

冰晶粒子运动过程中的相变特性

针对冰晶在温暖环境下运动时固液汽耦合的相变现象，基于欧拉法建立了冰晶粒子的运动相变模型和计算方法。计算了冰晶粒子在强迫对流环境下的融化相变过程，与实验结果对比验证了运动相变模型和计算方法的准确性。针对NACA0012翼型计算了冰晶绕流运动时的热力学特性，得到了冰晶粒子到达撞击表面时的融化状态与收集系数。研究了冰晶粒径大小、初始球形度、气流相对湿度和温度对运动相变的影响。结果表明：冰晶粒子运动相变模型可以有效地评估冰晶结冰风险，冰晶粒子的融化速率主要取决于粒子直径、球形度、气流温度、湿度等因素，环境温度为288 K时冰晶粒子的融化时间为27.5 s，而相同条件下环境温度为302 K时的融化时间仅有5.2 s。     

曲壁面两相流的简单扩散界面浸没边界法

针对在欧拉-拉格朗日网格框架下曲边界上的两相流润湿边界条件实施问题，本文提出了一种简单扩散界面浸没边界法。在此方法的框架下，静止和移动曲边界的湿润边界条件都可以采用平板边界方式处理。在该方法中，我们在实际边界内外一个拉格朗日网格间距处添加了两层虚拟贴体拉格朗日网格，并通过插值得到三层拉格朗日点上的序参数，从而得到实际的曲边界润湿条件。该方案实施思路简单，不仅适用于静止曲边界，还适用于复杂的运动曲边界。为了验证该方法的有效性，本文在广义守恒相场格子玻尔兹曼方法（LBM）框架下，模拟了静态液滴在圆柱上的铺展、动态液滴撞击圆柱、圆柱在水池中沉降三个问题。结果表明该方法能够有效地模拟曲面上的两相流动，具有良好的精度、鲁棒性和质量守恒性。