火星是距离地球最近的行星，目前各国的深空探测计划都将火星进入、火星取样返回等任务纳入规划。火星进入探测器面临严峻的气动热环境，须对其设计有效的热防护系统。与地球大气不同，火星大气密度较低且含有高浓度的CO₂，其进入过程存在更为显著的热化学非平衡效应；此外，火星会存在不同程度的沙尘暴，持续时间长且覆盖了探测器的下降高度，对探测器造成额外的侵蚀损伤，给探测器的热防护系统设计带来巨大挑战。为了对飞行器进入过程中含尘流场有更加深入的认识，本文采用数值模拟法方法，针对飞行器含颗粒高速气体来流开展气固耦合问题研究。本文主要研究内容概括如下：

1. 发展了二维轴对称五组分CO₂热化学非平衡计算程序，并针对不同运行工况下的飞行器外流场进行数值模拟计算，与该领域高知名度求解器LUARA(NASA)、LeMANS(UMich)进行了结果的对比与分析。使用力学所JFX激波风洞进行了自由流均匀性以及进入舱模型试验，通过对自由流皮托压力、激波脱体距离进行试验与计算结果对比分析。通过试验验证了低焓(1.3 ~ 1.4 MJ/kg)情况下热非平衡计算的准确性，并分析了喷管出口的流动状态；进一步通过与文献数据对比验证了中高焓(5.6 ~ 12.3 MJ/kg)热化学非平衡计算的准确性，为后续章节发展流固耦合计算程序提供支撑。

2. 基于Euler-Lagrangian框架，建立了非平衡流场与颗粒的单向耦合计算方法，针对颗粒计算的阻力模型以及颗粒—壁面碰撞模型进行相关研究及讨论，选取Henderson阻力模型作为颗粒的阻力计算模型；此外，采用模态半径为0.35 μm的火星大气颗粒分布模型，研究了不同尺寸颗粒在流场中的运动轨迹，并探究了高温相变模型对颗粒运动的影响规律以及不同粒径颗粒的撞击能量分布。

3. 在单向耦合的基础上，建立了双向耦合计算模型，详细介绍了双向耦合假设下的控制方程、颗粒源项具体计算方法、MPI并行条件下颗粒信息储存及传递方式，通过数值模拟方法研究了来流含颗粒前后的激波脱体距离、壁面热流变化结果，针对不同质量浓度的颗粒流条件，对流场结构的影响结果进行讨论与分析。

综上所述，本文的主要意义在于建立了针对火星进入环境下含尘气体与外流场流动的流固耦合计算方法，发展了E-L框架下的单向、双向耦合模型，并完成了数值模拟对含颗粒高速气体与飞行器相互作用的初步分析，为后续深入研究奠定了基础。研究结果有助于理解颗粒—流场—飞行器的多物理场耦合，推动高超声速火星进入飞行器相关研究的发展。