学术研究报告：空间站舱内对称通风方式的数值模拟研究
一、 作者、机构与发表信息
本研究由哈尔滨工业大学市政环境工程学院的郑忠海、张吉礼、梁珍三位研究者共同完成。相关研究成果以论文《空间站舱内侧上送侧下回对称通风方式数值模拟》的形式，发表于《哈尔滨工业大学学报》（Journal of Harbin Institute of Technology）第39卷第2期，出版时间为2007年2月。
二、 研究背景与目的
本研究属于航天器环境控制与生命保障（ECLSS）工程与计算流体力学（CFD）交叉的学术领域。随着载人航天事业的发展，为航天员创造一个安全、舒适且高效的舱内工作与生活环境至关重要。舱内环境需要维持特定的温度、湿度、气流速度与气体成分。在微重力条件下，由密度差引起的自然对流几乎消失，使得舱内空气的热量传递、污染物扩散以及湿度分布主要依赖于强迫通风系统。因此，设计合理的通风方式，以实现舱内空气的有效组织流动、均匀的温度分布和高效的换热，是航天器环境控制的关键技术问题。
本研究的核心目的是通过数值模拟方法，深入探究一种特定的通风方式——“侧上送侧下回对称通风”——在空间站舱内微重力环境下的流动与换热特性。具体研究目标包括：1）验证所采用的CFD数值方法在模拟微重力对流换热问题上的准确性；2）分析该对称通风方式在不同工况（如送风速度、角度）下的流场与温度场特征；3）特别关注一种可能的非线性流动现象——分叉解（Bifurcation）的出现条件与特征；4）量化评估微重力水平对舱内对流换热的影响程度，为实际工程设计与参数选择提供理论依据和具体建议。
三、 详细研究流程与方法
本研究的工作流程主要包含两个相互关联的数值模拟部分：基础验证与目标工况模拟。
第一部分：基础验证——封闭方腔自然对流模拟 为了确保后续对复杂舱内通风模拟的可靠性，研究首先选择了一个经典的基准问题：二维封闭方腔内的自然对流换热。该方腔上下壁绝热，左右壁分别为恒定高温和低温壁面。研究者建立了该问题的无量纲控制方程（连续性方程、动量方程、能量方程），其中引入了瑞利数（Ra）和努塞尔数（Nu）等关键无量纲参数来表征流动与换热强度。 * 研究对象与工况设置： 研究对象为一个二维方形计算域。设置了从Ra=1×10³到1×10⁶的不同工况，以覆盖从弱到强的自然对流范围。同时，为了研究微重力的影响，设置了从1g₀（地面重力）到10⁻⁶g₀的不同重力水平进行模拟。 * 数值方法： 使用商业CFD软件FLUENT进行求解。计算域采用80×80的结构化网格。对于层流工况（Ra < 10^8），直接采用层流模型；对于更高Ra数可能出现的过渡或湍流，采用标准k-ε湍流模型，并配合壁面函数法处理近壁区域。求解器设置方面，压力-速度耦合采用SIMPLE算法，动量与能量方程的离散采用二阶迎风格式，压力插值采用PRESTO格式。密度处理采用Boussinesq假设，即仅在动量方程的浮力项中考虑密度随温度的变化。 * 数据处理： 主要计算结果为沿热壁高度方向的平均努塞尔数Nu_ave，用于衡量换热强度。将模拟得到的Nu_ave与已有文献（如Davis的基准解）进行对比，以验证数值方法的准确性。
第二部分：目标研究——舱内侧上送侧下回通风模拟 在验证了数值方法的基础上，研究转向核心目标：对空间站舱内的强迫对流换热进行模拟。 * 物理模型： 将舱内截面简化为一个2m×2m的二维空间。通风方式定义为“侧上送侧下回对称”，即两个送风口对称布置于顶部侧边，两个回风口对称布置于底部侧边。送风口尺寸为0.1m。定义了关键参数：送风速度v_in、送风角度α（射流方向与顶壁的夹角）、雷诺数Re和格拉晓夫数Gr。 * 控制方程与假设： 研究分别考虑了层流和湍流两种流态。控制方程在基础验证部分方程的基础上，增加了来流惯性项。引入了无量纲数Gr/Re²来衡量自然对流与强迫对流的相对重要性。当Gr/Re² << 1时，自然对流影响可忽略；当Gr/Re² >> 1时，自然对流主导；两者量级相当时为混合对流。分析指出，在典型的微重力（10⁻³g₀）、低速（0.2 m/s）舱内环境中，Gr/Re²值极小（约0.002），表明自然对流作用微弱。 * 数值模拟方案： 针对不同的送风速度（v_in = 0.7 m/s, 1.0 m/s）和送风角度（α = 30°, 45°, 60°）组合成多种工况进行模拟。同样使用FLUENT软件，根据计算雷诺数判断流态并选择相应模型（层流或标准k-ε模型）。边界条件设定为左右壁面等温，上下壁面绝热。研究特别关注了在特定工况下流场可能出现的“分叉解”现象，即虽然几何和边界条件完全对称，但最终的稳定流场却可能呈现两种不对称的形态。为了捕捉这两种可能的稳定解，计算时需采用不同的初始扰动条件（如分别以左侧或右侧送风口为主导的初始流场）进行全场计算。 * 数据分析： 对每个工况，分析流场形态、温度场分布、全场平均流速以及壁面平均努塞尔数等结果。重点比较了不同重力水平（1g₀ vs. 10⁻³g₀）下结果的差异，以及分叉解出现时的流场特征。
四、 主要研究结果与分析
1. 基础验证结果： 数值模拟得到的封闭方腔在不同Ra数下的平均努塞尔数Nu_ave与经典的Davis基准解吻合度极高，最大误差不超过3%。这充分证明了研究所采用的CFD模型、网格划分策略和求解设置对于模拟自然对流换热问题是准确可靠的。此外，对不同微重力水平下的模拟数据进行拟合，得到了Nu_ave与Ra数的关系式（Nu_ave=0.142Ra^0.298），与Fusegi等人提出的关系式趋势一致。模拟结果明确显示，当微重力水平低于10⁻³g₀时，Nu数接近于1（纯导热的理论值），表明自然对流效应已非常微弱，可以忽略不计；当重力水平在10⁻³g₀时，其换热能力仅为地面重力（1g₀）下的约15%。
2. 舱内通风模拟核心发现： * 分叉解现象及其机理： 研究在特定工况（如Re_in=5000, α=45°，微重力10⁻³g₀）下，成功模拟并观察到了“分叉解”现象。计算结果得到了两种稳定的、镜像对称的不对称流场。例如，“分叉解1”表现为左侧送风射流贴附顶壁，舱内形成的大涡旋中心偏右；而“分叉解2”则恰好相反。研究者用“柯恩达效应”（Coanda Effect）来解释这一现象：射流与壁面之间形成的负压区会使气流产生贴附壁面的趋势。一个微小的初始扰动（在数值计算中体现为不同的初始条件）导致其中一股射流率先偏离轴线贴附壁面，就会破坏原有的对称性，并诱导整个流场形成一种稳定的不对称模式。这一发现对舱内气流组织设计具有重要意义，提示在对称通风设计中仍需考虑不对称流场出现的可能性。 * 通风参数的影响与优化建议： 通过对不同送风速度和角度的组合工况进行模拟，研究得出了以下定量结论： * 重力影响判据： 自然对流的影响程度可用无量纲数Gr/Re²来表征。模拟结果表明，当Gr/Re² < 5时，忽略重力影响（即认为微重力与常重力结果一致）所带来的误差不超过5%。这为工程应用提供了一个简便的判据。 * 微重力等效性： 当微重力水平低于10⁻³g₀时，其产生的流动与换热效应已非常微弱，在工程上可近似视为零重力条件。 * 送风角度优化： 比较不同送风角度下的流场平均流速发现，当送风角度α为45°时，舱内获得的平均流速最大。这表明45°角可能在促进空气混合和循环方面效率较高。 * 参数建议： 综合流场均匀性和换热效果，研究建议送风角度在30°至60°之间选择，其中45°为佳；同时，送风速度不宜过高，建议小于1.0 m/s，以利于形成稳定、均匀的流场。
五、 研究结论与价值
本研究的核心结论是：采用计算流体力学（CFD）方法对空间站舱内环境进行数值模拟，是研究微重力下通风换热问题的有效手段。针对“侧上送侧下回对称通风”方式，研究发现其在特定条件下会出现分叉解，导致不对称的流场。研究明确了自然对流影响的可忽略条件（Gr/Re² < 5），并给出了微重力（<10⁻³g₀）可视为零重力的工程判断。最终，为实际工程应用提出了关键的通风参数选择建议：送风角度宜采用30°~60°，最佳为45°；送风速度应控制在1.0 m/s以下。
本研究具有重要的科学价值与应用价值。在科学层面，它深入揭示了微重力环境下封闭空间内强迫对流与残余自然对流耦合的复杂流动现象，特别是对“分叉解”这一非线性动力学行为在工程通风系统中的出现条件与特征进行了实证分析，丰富了微重力流体物理的知识体系。在应用层面，该研究为航天器舱内通风系统的设计提供了直接的理论依据和具体的参数优化指导。所建立的数值模拟方法与分析判据，可用于预测和评估不同通风方案下的舱内环境品质，减少对昂贵且复杂的地面微重力模拟实验或空间实验的依赖，从而降低设计成本，提高设计可靠性。
六、 研究亮点
七、 其他有价值的发现
研究还指出，分叉解的出现需要特定条件，如送风口为对称长缝风口且空间内部较为开阔无障碍。如果送风口是分散布置（三维结构）或者空间中存在较大障碍物，则可能抑制分叉解的出现。这一补充说明有助于更全面地理解该现象的应用边界，为设计避免或利用分叉解提供了更细致的思路。