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

将基于风洞网格测力试验数据建立的气动力模型与刚体运动方程进行耦合求解得到多体分离轨迹-时间特性，建立了一种多体分离的离线轨迹预测方法。为了提高气动力模型的预测精度，针对移动最小二乘法（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）框架下，模拟了静态液滴在圆柱上的铺展、动态液滴撞击圆柱、圆柱在水池中沉降三个问题。结果表明该方法能够有效地模拟曲面上的两相流动，具有良好的精度、鲁棒性和质量守恒性。     

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

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

高焓流动中的可压缩颗粒求解器（第1部分）：考虑多物理效应的点力颗粒两相流理论方程

从含烧蚀颗粒的高超声速边界层到航空发动机燃烧室或火箭发动机燃烧室里的气液、气固两相燃烧等问题，都存在弥散相颗粒调制携带流体的两相耦合物理过程。尽管不可压颗粒两相流的理论、数值工具和实验都比较成熟，但可压缩颗粒两相流与不可压缩的情况是截然不同的。本文首先对可压缩颗粒两相流动力学和热力学进行理论研究，理论分析了经典的可压缩流体动力学和热力学方程，把颗粒动力学对不可压流动的调制理论修正拓展到了可压缩情况，发现了可压缩颗粒两相流特有的动量和能量调制无量纲参数；推导了考虑多物理效应的颗粒动力学方程，并将其与携带流体相的控制方程耦合。研究发现，尽管多物理效应下的颗粒动力学是复杂的，但其影响携带流体相的路径是唯一确定的，即：多物理效应的颗粒动力学仅能通过相间阻力调制携带流体相的动量，通过相间阻力做功和热对流调制携带流体相的能量。