关于三种CFD软件在非布辛涅斯克条件下方形空腔自然对流模拟的基准测试研究报告
一、 主要作者、机构与发表信息 本项研究由Sneha R. Sondur, Brek Meuris, 和 Ann M. Mescher共同完成。三位作者均来自Department of Mechanical Engineering, University of Washington, Seattle, USA。该研究成果以学术论文形式发表，最终在线发表于2022年7月18日，刊载于学术期刊 Numerical Heat Transfer, Part A: Applications 2023年第84卷第3期（页码252–268）。
二、 研究的学术背景 本研究属于计算流体力学（Computational Fluid Dynamics， CFD）领域，具体聚焦于自然对流（Natural Convection）的数值模拟验证与软件性能评估。自然对流广泛存在于诸多工程应用中，例如建筑通风、电子设备冷却、太阳能接收器传热以及高温制造工艺（如玻璃熔化和聚合物纤维拉制）等。随着高性能计算集群的普及，CFD已成为研究和设计相关系统的重要工具，其预测精度直接关系到工程设计的可靠性。
然而，绝大多数现有的自然对流CFD基准测试和文献参考解（Reference Solutions）均基于布辛涅斯克近似（Boussinesq Approximation）。该近似假设流体物性（除密度在浮力项外）为常数，仅适用于驱动温差（ΔT/T）趋近于零（即dt/T -> 0）的工况。对于许多实际高温差应用（例如温差超过28°C的空气），流体密度、粘度、导热系数等物性随温度变化显著，布辛涅斯克近似不再适用。此时，必须采用完全可压缩形式的Navier-Stokes方程，并考虑物性的温度依赖性，才能准确预测流场和传热特性。尽管已有Martineau等人（2010年）的研究表明，对高压差条件下的自然对流空气流，使用不可压缩方程求解会导致显著误差，但针对非布辛涅斯克条件下低马赫数可压缩流的可靠CFD基准测试研究仍然相对匮乏。
为解决这一问题，INRIA与MAB机构于2004年组织了一次数值研讨会，旨在为高温差驱动的方形空腔内自然对流空气流建立可靠的参考解。这些参考解（即文中的T1, T2, T3三个测试案例）为后续CFD算法的验证提供了宝贵基础。
基于此背景，本研究设定了三个具体目标：1）通过将当前模拟结果与2004年研讨会建立的基准解对比，评估七种不同CFD算法对于低马赫数可压缩空气流的相对精度。2）从中央处理器（CPU）耗时和精度（包括空间平均量与局部量）两方面，对比不同CFD算法的性能。3）提供具有变物性的方形空腔空气流（T2和T3案例）的稳态速度分布、局部努塞尔数、压力和温度等详细数据，以作为未来CFD代码基准测试的额外参考标准。
三、 研究详细工作流程 本研究是一项系统的数值基准测试研究，其核心工作流程为：定义标准问题 → 设置求解器 → 进行网格收敛性测试与误差估计 → 执行计算并获取结果 → 与权威参考解进行多维度对比。
1. 基准问题定义： 研究模拟了一个方形空腔（边长L）内的稳态自然对流空气流。边界条件为：两个垂直壁面分别保持高温（Th）和低温（Tc），两个水平壁面为绝热。重力加速度g沿y轴负方向。初始条件为：腔内空气静止，初始温度To = (Th+Tc)/2，初始压力Po = 101,325 Pa，所有物性按初始温度To计算。研究中定义了三个基准测试案例：T1（瑞利数Ra=10⁶，无量纲温差ε=0.6，变密度但粘度和导热系数为常数）；T2（Ra=10⁶， ε=0.6，密度、粘度、导热系数均为温度的函数）；T3（Ra=10⁷， ε=0.6，所有物性均为温度的函数）。所有模拟均假设空气为理想气体，粘度由萨瑟兰定律（Sutherland’s Law）计算，导热系数通过普朗特数Pr（恒定=0.71）关联得到，比热容及比热比为常数。
2. 求解器与算法选择： 本研究共测试了来自三个软件包的七种不同算法，涵盖了密度基（Density-Based）和压力基（Pressure-Based）两大类求解器，以及稳态（Steady）和瞬态（Transient）两种求解方式。具体包括： * Simcenter STAR-CCM+ 求解器 1：密度基、稳态求解器（采用伪瞬态项）。 * Simcenter STAR-CCM+ 求解器 2：密度基、瞬态隐式求解器（采用双时间步长法）。 * ANSYS-Fluent 求解器 3：密度基、稳态求解器（采用伪时间步长和牛顿型线性化）。 * ANSYS-Fluent 求解器 4：密度基、瞬态求解器（采用二阶后向差分和双时间步长法）。 * ANSYS-Fluent 求解器 5：压力基、稳态求解器（无时间项，使用欠松弛）。 * ANSYS-Fluent 求解器 6：压力基、瞬态求解器（二阶隐式/迭代时间格式，动量方程与压力连续性方程耦合求解）。 * OpenFOAM 求解器 7：压力基、稳态求解器（buoyantSimpleFoam，采用SIMPLE算法）。
所有模拟均采用二阶迎风离散格式处理对流项，二阶中心差分格式处理扩散项。
3. 网格与收敛性研究： 为评估数值误差并确定算法的渐近收敛阶，研究采用了逐级加密的均匀正交网格进行网格收敛性研究。使用了三种网格分辨率：160x160（Grid 3）、320x320（Grid 2）和640x640（Grid 1）。网格细化比为2。对于瞬态求解器，时间步长统一设置为0.001秒。收敛判据为：空腔左右壁面的换热量差异小于10⁻⁷瓦特。同时监测了连续性方程、动量方程和能量方程的残差。对于瞬态求解器，还监控了每个时间步的局部速度和温度以确保收敛。
4. 验证与误差估计方法： 研究从多个维度验证数值解的可靠性和准确性： * 质量守恒：计算初始总质量与最终总质量的差异。 * 能量守恒：验证左右垂直壁面（绝热）的平均努塞尔数是否相等，并计算热通量不平衡度。 * 阶次收敛：利用三种网格的解，根据Roache提出的方法计算算法的实际观察收敛阶（Observed Order of Convergence），并与理论收敛阶（二阶）进行比对。 * 基准对比：将模拟得到的关键结果（包括空间平均的努塞尔数Nu、热力学压力比Pt/Po、以及沿加热壁和冷却壁特定位置的局部努塞尔数、空腔内最大压力Pmax/Po等）与INRIA和MAB会议（特别是Vierendeels等人使用2048x2048网格得出的最精确解）公布的权威参考解进行详细对比，并计算百分比误差。
所有计算均在美国华盛顿大学的Hyak高性能计算集群上完成，每个模拟使用16个核心和120GB内存。
四、 主要研究结果 1. 稳态与瞬态求解器的性能对比： 这是本研究的核心发现之一。使用320x320网格时，两个商业软件的瞬态密度基求解器（STAR-CCM+ Solver 2和Fluent Solver 4）在预测T1, T2, T3三个案例的平均努塞尔数Nu和热力学压力比Pt/Po时，误差均小于0.7%（使用640x640网格时误差小于0.2%）。然而，商业软件的稳态密度基求解器（STAR-CCM+ Solver 1和Fluent Solver 3）则产生了显著误差：预测Nu的误差在T1案例中高达9%，T2和T3案例中约为5%；预测Pt/Po的误差在T2和T3案例中约为8%，T1案例中高达17%。误差根源在于质量不守恒：稳态求解器导致了非物理的质量积累（例如T1案例中质量增加了117%），进而引发了压力的错误预测。这验证了对于高温差非布辛涅斯克流动，稳态密度基求解器的不适用性。
2. 高精度解与收敛性分析： 采用STAR-CCM+瞬态隐式求解器和640x640网格，可以得到与基准解高度吻合的结果。对于T3案例，平均Nu和Pt/Po的最大误差仅分别为0.16%和0.07%。局部努塞尔数的最大误差出现在冷却壁面（1.4%），这归因于冷却壁附近的边界层更薄，需要更细的网格来精确捕捉。收敛性分析表明：对于常物性案例T1，观察到的收敛阶与理论值（2阶）相符；而对于变物性案例T2和更高瑞利数的T3案例，观察收敛阶分别降至1.9和1.8。这表明对于物性变化显著或高瑞利数的复杂流动，需要更精细的网格才能达到网格无关解。
3. Fluent内部压力基与密度基瞬态求解器对比： 使用640x640网格时，Fluent的瞬态压力基求解器（Solver 6）和瞬态密度基求解器（Solver 4）在精度上表现相当，预测平均量的误差均小于0.15%，预测局部Nu的误差也处于相近水平（如T3案例中冷却壁最大Nu误差分别为1.4%和1.3%）。两者预测的沿壁面局部Nu分布曲线几乎完全重合。
4. OpenFOAM稳态求解器的表现： OpenFOAM的稳态压力基求解器（buoyantSimpleFoam）在320x320网格上表现尚可。对于T3案例，平均Nu误差约为2.5%，Pt/Po误差约为0.27%。虽然精度略低于商业瞬态求解器，但该算法的关键优势在于其严格的质量守恒特性，且计算效率最高。
5. 计算效率（CPU时间）对比： 计算成本是工程应用中的重要考量因素。研究详细对比了不同算法的CPU时间： * 商业稳态求解器虽然计算时间最短（例如Fluent稳态密度基求解器仅需数分钟到数十分钟），但这是以牺牲精度和违反质量守恒为代价的，其结果不可信。 * 商业瞬态求解器精度高，但耗时最长。其中，Fluent的瞬态压力基求解器比其瞬态密度基求解器平均节省约35%的CPU时间，在保持精度的同时更具效率优势。 * OpenFOAM稳态求解器在320x320网格上的计算时间，比Fluent瞬态压力基求解器平均节省约55%，在精度可接受（尤其对于初步分析）的情况下，提供了显著的效率优势。 * 研究还提出并验证了一种实用的加速策略：首先运行一个（即便不准确的）稳态求解器以获取一个接近最终稳态的温度场，然后将其作为瞬态求解器的初始条件。 对于T3案例，这种“两步法”使STAR-CCM+瞬态求解器的收敛时间缩短了约3.4倍。
6. 详细的流场与温度场数据： 研究提供了由Fluent瞬态求解器预测的T3案例（Ra=10⁷, ε=0.6）在640x640网格上的详细结果，包括速度矢量图、无量纲温度等值线图、沿水平和垂直中心线的温度分布、以及沿不同水平线的速度和垂直速度分量分布。这些数据清晰展示了非布辛涅斯克流动的非对称特性（例如上下角涡结构不对称），并揭示了加热壁和冷却壁附近热边界层的差异，为未来的CFD验证提供了宝贵的局部参考数据。
五、 结论 本研究对开源代码OpenFOAM和两个商业有限体积软件包（STAR-CCM+和ANSYS-Fluent）在模拟高温差驱动方形空腔自然对流空气流方面的准确性和收敛速度进行了系统的基准测试。结论如下：对于所研究的非布辛涅斯克、低马赫数可压缩流动，商业软件中的瞬态密度基和压力基求解器预测结果最为准确，但计算成本最高。其中，瞬态压力基求解器在精度相当的情况下，比瞬态密度基求解器更具计算效率优势。开源软件OpenFOAM的稳态压力基求解器（buoyantSimpleFoam）虽然精度略低，但严格满足质量守恒，且计算速度最快，是进行初步、快速分析的可行工具。而商业软件的稳态密度基求解器由于存在严重的质量不守恒和压力预测错误，不适用于此类高温差自然对流问题。此外，研究表明，采用从（不准确的）稳态解获取的温度场作为瞬态求解的初始条件，可以大幅减少获得高精度解所需的CPU时间，为实际工程应用提供了有效的加速策略。
六、 研究亮点 1. 明确的针对性：首次系统地对多种主流CFD软件（包括开源和商业）在非布辛涅斯克近似条件下的自然对流模拟能力进行了基准测试，填补了该特定领域系统性软件评估的空白。 2. 全面的算法覆盖：研究同时涵盖了密度基和压力基、稳态和瞬态共七种不同的求解算法，提供了全方位的性能对比图景，结论具有广泛的参考价值。 3. 深刻的机理揭示：不仅报告了性能差异，更深入分析了差异的根源。明确指出商业稳态密度基求解器产生大误差的根本原因在于算法导致了非物理的质量积累，从而引发连锁错误。 4. 多维度的评估指标：评估不仅限于传统的精度（平均Nu、压力）对比，还扩展至CPU时间成本、局部量（局部Nu、速度/温度分布）精度、实际收敛阶、以及质量/能量守恒性等多个维度，评估体系全面。 5. 实用的工程建议：研究不仅止于“测评选优”，还提出了“利用稳态解温度场初始化瞬态求解” 这一具体且高效的计算策略，并量化验证了其加速效果（~3.4倍），对工程实践具有直接指导意义。 6. 丰富的参考数据：研究不仅验证了已有基准解，还生成了高分辨率网格下的详细流场、温度场和局部Nu分布数据，为后续研究者和工程师提供了新的、更丰富的验证参考标准。
七、 其他有价值内容 研究引用了INRIA和MAB 2004年研讨会建立的基准解作为权威验证标准，并详细说明了该基准问题的定义、边界条件和物性模型，确保了本研究工作流程的严谨性和可复现性。文中对密度基和压力基求解器在处理低马赫数可压缩流时的不同策略（如预处理、渐近近似、压力修正等）进行了简要的背景介绍，有助于读者理解不同算法差异的理论基础。最后，所有计算均基于华盛顿大学Hyak超级计算集群完成，并获得了美国国家科学基金会（NSF）的资助（Grant CMMI-1301491），体现了研究的学术规范性和资源支持情况。