高超声速流动因其独特的物理化学特征,难以采用传统实验方法或连续介质模拟方法全面研究。本论文通过改进硬球-拟颗粒模拟(HS-PPM),显著降低了离散模拟的计算量、提高了计算效率,从而改进了该方法模拟高超声速流动的适用性。HS-PPM原本是针对气固多相流动提出的一种离散模拟方法,已成功应用于纳微尺度气体流动-扩散耦合过程的模拟。本论文为探索其在高超声速流动领域的适用性,首先分析了高超声速飞行器再入过程的流动特征,并调研了国内外在该领域的研究进展。在调研基础上,改进了HS-PPM中的壁面描述、气体表面作用模型、边界条件设置等,使其更适合描述高速流动。在此基础上论文对比了硬球模型(HS)、HS-PPM和直接模拟蒙特卡洛(DSMC)方法模拟同一高超声速圆球绕流问题的结果,表明HS-PPM的结果较DSMC更接近HS。对雷诺数为100、马赫数从4到19的三维圆球绕流,在HS-PPM中分别采用完全热边界和完全滑移边界得到了阻力系数的上下限,并与相应的HS结果吻合。上述研究表明了HS-PPM处理高超声速稀薄气体流动的可行性。通过外推填充率趋于0的圆球绕流结果、模拟马赫数为24而努森数从0.11到4.55的零攻角三维尖锥绕流及马赫数为20而努森数为0.03的长尾钝锥绕流等实例,分析其误差并指出了HS-PPM改进的方向。作为初步的误差修正,对HS-PPM引入了Larsen-Borgnakke能量模型。模拟的粒子具有平动能和包括转动、振动的内能,并根据平衡分布赋予粒子能量初值。在每次碰撞中通过随机取样更新粒子能量使碰撞后粒子不同能量模式的宏观统计值符合各自的平衡分布。对长尾钝锥的对比模拟证明能量模型的引入基本修正了无反应HS-PPM结果中激波温度偏高的问题。但对Apollo返回舱模型以马赫数33再入努森数约为0.081的大气层的模拟表明,当真实气体效应显著时仅考虑能量转化不足以修正温度误差。本论文继而提出了简单的氮氧离解反应模型,将其中的能量引入Larsen-Borgnakke 能量模型并初步实现了粒子分裂模型及其算法,从而显著改进了Apollo返回舱模型算例中HS-PPM的温度场分布。本论文的研究表明并改进了 HS-PPM对高超声速稀薄流动的适用性。但本论文的反应模型只考虑了比较简单的空气组分和反应网络,未来的研究还应从复杂反应模型、可变碰撞截面等多方面继续改进HS-PPM。