技术分享︱多重参考系模型在风扇通风仿真中的自动化实现：精度与效率的工程平衡

引言在厂房通风、数据中心散热等封闭空间的流体力学仿真中旋转风扇对全流场流速分布的驱动作用是计算的核心。传统的仿真工作流通常依赖人工进行几何处理、网格划分与求解器设置难以满足海量工况的快速评估需求。为解决这一痛点神工坊®技术团队开发了一套基于 HPC 集群的SimForge™仿真应用平台。该平台将底层复杂的网格拓扑与物理模型封装实现了从参数输入到结果输出的全自动化。在构建该平台的核心求解逻辑时如何在“计算精度”与“求解效率”之间取得最佳平衡是算法选型的关键。01 MRF模型选型论证多重参考系模型应用示例在处理包含旋转机械的计算流体力学问题时静止域与旋转域动静干涉边界的处理是求解的核心难点。当前工业界针对该类问题的主流处理模型主要分为两类瞬态滑移网格模型Sliding Mesh Model, SMM与稳态多重参考系模型Multiple Reference Frame, MRF。瞬态滑移网格模型SMM基于网格的真实物理运动进行非定常求解能够高保真地捕捉转子与定子交互时的强瞬态气动效应如尾迹脱落与流场脉动。然而SMM 求解严格受限于 Courant 稳定性条件要求极小的物理时间步长。这导致计算资源与时间成本呈指数级增长难以匹配 SaaS 平台对高并发与计算结果高效流转的工程需求。多重参考系模型MRF则采用稳态近似策略将全局计算域在空间上严格解耦为静止域如厂房宏观环境与旋转域如风扇叶轮近场区域。在数学处理上静止域内的纳维-斯托克斯Navier-Stokes方程在绝对参考系下封闭求解而旋转域内的控制方程则转换至跟随转子同步旋转的相对参考系下求解最终通过网格交界面实现动量与质量通量的守恒传递。单旋转叶轮几何模型基于宏观封闭空间通风的工程场景特性全流场的时间平均分布特征如宏观风流分布、有效风速覆盖率远比局部的强瞬态动静干涉现象更为关键。因此本平台底层求解器选型确立为 MRF 模型并由此引出“冰冻转子法”的工程论证。在物理等效性上MRF 模型不考虑旋转域与相邻静止域之间的真实相对位移相当于将转动部件在其某一瞬态几何相位上“冻结”进而求解该位置下的稳态流场。基于此物理近似特征该理论方法在工程计算中被称为“冰冻转子法”Frozen Rotor Approach。该方法在确保宏观流场计算精度的前提下大幅降低了计算开销实现了求解效率的最大化是支撑大规模云端自动化仿真落地的理想核心算法。02 速度变换与算子映射基于前述的MRF模型冰冻转子法旋转域与静止域采用不同的参考系描述。为实现静止网格上的旋转效应模拟其核心在于速度场的参考系变换以及动量方程的重构。与瞬态滑移网格方法不同MRF 求解器并非物理移动网格而是通过算子对控制方程进行运动学与动力学变换将旋转坐标系下的惯性力如科里奥利力、离心力以源项形式引入。2.1相对速度公式化与附加源项算子静止与运动参考系在旋转非惯性系中流体质点的绝对速度v 与相对速度r 之间由牵连速度算子联系其中ω为旋转角速度矢量r为位置矢量。若以相对速度v作为动量方程的求解变量由于参考系的旋转必须在动量方程中显式引入两个附加的体积力源项算子以维持力学平衡。其一为科里奥利力算子Coriolis Force Operator表达式为该算子表征流体在旋转参考系中相对运动时因参考系旋转而受到的切向与法向偏转效应。其二为离心力算子Centrifugal Force Operator表达式为该算子表征流体随参考系转动时受到的背离旋转轴向外的径向惯性力。2.2绝对速度公式与算子简化在风扇全流场模拟中主流底层求解器通常默认采用绝对速度公式化。该方法以绝对速度v作为求解变量但在计算跨越网格面的对流质量通量时依然使用相对速度v。通过严格的代数推导上述复杂的科里奥利力算子与离心力算子可以被合并与化简最终在动量方程中表现为单一的简化旋转源项算子化简后的稳态动量方程形式为该公式化的工程优势在于算子的简化不仅降低了计算复杂度更重要的是在具有强烈旋流的区域如风扇尾流区绝对速度公式化能够有效避免离散化过程中的数值耗散从而显著提高全场求解的鲁棒性和精度。2.3MRF交界面公式化与通量传递由于计算域被划分为静止域如厂房和旋转域如风扇叶轮两者的控制方程参考系不同。实现“冰冻转子法”的关键在于交界面Interface处的数据传递。首先是局部参考系转换算子的应用。当底层求解器在旋转域边界计算通量并向静止域传递时系统会自动应用速度矢量的空间转换算子。其次是通量守恒的严格保证。无论是采用相对速度还是绝对速度公式化流经交界面的质量通量Mass flux、动量通量Momentum flux和其他标量通量都会通过关系严格转换为相邻域所在参考系下的表达形式确保物理量在非一致性Non-conformal网格交界面上的绝对守恒。03 自动化实现与调度在 SaaS 平台的后台架构中上述多重参考系理论通过 Python 脚本与底层求解器的高效交互得以实现。该自动化工作流精确映射了MRF 模型的关键边界设置与求解需求其核心落地流程涵盖以下三个技术环节。3.1 区域智能识别与交界面构建算法自动解析几何空间的拓扑特征将全局计算域自适应划分为静止域与多个相互独立的旋转域。在此基础上后台脚本自动生成底层交互指令创建非一致性网格交界面。该机制确保了动量通量与质量通量在不同参考系边界处的数据插值与严格守恒传递。3.2 运动参考系的动态映射针对多风扇阵列的复杂工况自动化脚本动态读取前端输入的风扇空间坐标与转速参数精准定位每个转子的旋转轴与原点。通过这一映射过程旋转角速度 ω 的物理定义被精确转化为底层求解器所识别的运动参考系边界条件实现了非惯性系算子的自动施加。3.3 HPC 并行调度与高效求解针对多风扇流场引入的千万级网格规模平台架构在底层深度集成了 Slurm 等高性能计算集群调度系统。通过脚本动态配置节点资源与通讯文件系统调用 OpenMPI 驱动底层 CFD 求解器执行多节点跨核并行计算。该架构在保证流场宏观计算精度的前提下实现了仿真任务的高并发与高效流转。04 结语本文从工程痛点出发深度解构了 SaaS 化通风仿真平台底层的核心计算理论。通过引入稳态多重参考系MRF模型并利用非惯性系下速度场与动量方程的算子变换神工坊®技术团队成功在计算精度与求解效率之间找到了最优解。结合 HPC 集群与 Python 自动化脚本复杂的物理算子与底层交互被完全黑盒化。这不仅大幅降低了多风扇联合仿真的门槛也为同类流体计算的 SaaS 化落地提供了可靠的技术参考。