《Aerospace》：DSMC Analysis of DiskSat in Very Low Earth Orbits
Máximo Castillo Rivas,
Diego Vera Sepúlveda and
Rodrigo Cassineli Palharini
编辑推荐：
为探究极低地球轨道(VLEO)中新兴平板卫星DiskSat的气动行为，研究人员采用直接模拟蒙特卡洛(DSMC)方法，分析了其在100-200 km高度下的宏观流场结构和气动特性。结果表明，DiskSat在不同高度下呈现从压缩激波层到高度扩散流场的演变，并存在热冻结近真空尾迹区，这为未来VLEO任务设计及轨道寿命分析提供了关键的气动热力学基准。
在太空探索的新浪潮中，卫星正变得越来越“接地气”——这里说的“接地气”是指它们正飞向越来越低的轨道。极低地球轨道（Very Low Earth Orbit, VLEO），通常指高度在150-450公里之间的轨道，因其能带来更高分辨率的地球观测、更低的通信延迟以及通过大气阻力自然清除轨道碎片等战略优势，正成为航天工业的新宠。然而，天下没有免费的午餐。VLEO环境中显著增高的高层大气密度带来了严峻挑战：严重的气动阻力会缩短卫星寿命，而高活性的原子氧（O）则可能导致表面材料被侵蚀。更棘手的是，地面上的真空设备很难完全复现这种极端环境下的气体-表面相互作用，使得预测卫星在轨性能变得异常困难。
与此同时，卫星的设计理念也在发生变革。以CubeSat（立方星）为代表的标准化平台曾大大降低了进入太空的门槛，但随着任务需求向高功率载荷（如合成孔径雷达SAR）演进，立方体构型在体积和表面积上已显捉襟见肘。行业开始将目光投向新型的平面化架构，其中一种极具代表性的设计便是DiskSat。想象一个直径1米、厚度仅2.5厘米的“飞盘”，这种极致的扁平化设计最大化其面积质量比，特别适合搭载需要大孔径的先进载荷，并易于在运载火箭整流罩内堆叠发射，代表了卫星“容器化”的发展趋势。尽管NASA和航空航天公司等机构的近期演示已证实了DiskSat机械结构的飞行准备状态，但一个核心科学问题依然悬而未决：这样一个巨大的平板结构，在VLEO高度下稀薄的大气中飞行时，其周围的气流究竟是何模样？那里的气体如此稀薄，以至于传统的流体连续性假设已然失效，必须从分子层面去理解其行为。这个问题直接关系到卫星的轨道寿命、姿态稳定性和热控设计，是DiskSat迈向VLEO应用必须闯过的“气动关”。
为了回答这个关键问题，由Máximo Castillo Rivas, Diego Vera Sepúlveda 和 Rodrigo Cassineli Palharini组成的研究团队在《Aerospace》期刊上发表了一项研究，他们采用了一种名为直接模拟蒙特卡洛（Direct Simulation Monte Carlo, DSMC）的先进数值模拟方法，对DiskSat在100公里、150公里和200公里三个VLEO高度下的气动热力学行为进行了首次系统的刻画，旨在为其建立一个基础的气动数据库。
本研究的核心是应用直接模拟蒙特卡洛（DSMC）这一粒子数值方法。DSMC适用于模拟稀薄气体流动（其稀薄程度由Knudsen数Kn量化），它通过跟踪大量模拟粒子的运动与碰撞来统计气体的宏观特性。研究使用了开源求解器dsmcFoam+，并采用了变径硬球（Variable Hard Sphere, VHS）碰撞模型、No Time Counter (NTC)碰撞对抽样算法、Larsen–Borgnakke内能交换模型，并在卫星表面应用了完全热适应的麦克斯韦漫反射模型，以建立气动力学的保守估计（即最坏情况）。研究首先对一个高超音速平板绕流的标准算例进行了验证，确保所用方法与代码的可靠性，随后对时间步长、网格、粒子数等参数进行了灵敏度分析，验证了数值结果的稳健性。最终，针对一个部分展开太阳能帆板、质量11公斤的DiskSat三维半模，在0度攻角下，分别对三个目标高度进行了DSMC模拟，系统分析了其流场结构。
4. Simulation Description
研究人员构建了一个包含DiskSat三维半模的计算域，其边界包括自由来流入口、真空出口、对称平面以及设定为固定温度的卫星壁面。计算网格在卫星表面附近进行了局部加密，以捕捉可能的高梯度流动区域，例如在100公里高度下更为致密的压缩区域。
3. Validation and Verification
为确保结果可靠，研究团队首先对一个高超音速平板绕流的标准基准案例进行了模拟，并将获得的流场速度与旋转温度云图与已有文献结果进行对比，结果显示二者高度吻合，验证了所采用DSMC方法和代码设置的正确性。进一步的参数灵敏度分析表明，所采用的时间步长、网格分辨率、单位网格内粒子数（PPC）以及时间平均采样区间均能获得网格无关且统计收敛的解，确保了后续DiskSat模拟结果的数值可信度。
宏观流场结构与气动特性
通过对三个高度（100公里、150公里、200公里）的模拟，研究清晰地揭示了DiskSat气动行为的演变规律：
- •
100公里高度（Kn≈0.14，过渡流区）：在此相对稠密的高度，DiskSat前方观察到一个压缩的激波层，其温度显著升高，表明气体分子在撞击卫星前表面时发生了强烈的压缩和能量交换。
- •
200公里高度（Kn≈240，自由分子流区）：由于气体分子平均自由程极大，气体与卫星的相互作用变得高度扩散，其影响范围向上游显著延伸，不再有明显的集中激波结构。
- •
普遍存在的尾迹区：在所有研究的高度下，DiskSat的尾迹区域都呈现出热力学状态近乎冻结、接近真空的状态。这是因为流向卫星后方的气体分子极度稀少，它们之间几乎不发生碰撞，无法通过碰撞达到热平衡，从而形成了一个低温、低压的尾流区。
这项研究系统地刻画了DiskSat平台在VLEO环境下的宏观气动热力学行为，建立了一个关键的基础数据集。主要结论是，DiskSat的流场结构强烈依赖于轨道高度（即Knudsen数Kn）。在较低高度（如100公里），其前方会形成压缩的激波高温区；而在较高高度（如200公里），气体相互作用变得高度扩散。一个共性的重要发现是，在所有高度下，DiskSat后方均存在一个热力学冻结、近真空的尾迹区，这是其平面几何构型在稀薄流中的典型特征。
这些发现具有重要的科学意义与工程价值。首先，它们首次详细描绘了这种大尺寸平板卫星在VLEO稀薄流中的独特气动“足迹”，弥补了现有文献的空白，现有基于较小几何体（如10厘米CubeSat）的标度律难以准确预测DiskSat的流场。其次，本研究获得的流场结构和气动数据，为后续精确计算DiskSat在VLEO中的轨道寿命、分析其姿态稳定性（特别是在过渡流和自由分子流区）提供了不可或缺的输入。研究者指出，尽管为了建立保守基线，本研究暂未考虑化学反应和高温气体效应（如振动激发），也未采用更复杂的气体-表面相互作用模型，但当前工作为未来纳入这些复杂物理过程、开展更精细化研究奠定了坚实的基础。总之，这项研究为DiskSat及其同类平面构型卫星在VLEO这一战略新兴轨道的任务设计与性能评估，迈出了至关重要的第一步。
