1. 湍流模拟的工程挑战与Fluent求解器定位

湍流是工程流体中最常见、也最昂贵的物理现象。它的本质是多尺度涡的随机叠加:大尺度涡从平均流中抽取能量,通过级串过程逐级破碎,最终在Kolmogorov尺度上被黏性耗散为热。这一过程使湍流同时具备两个让数值模拟头疼的特征——尺度跨度极大与强非线性耦合。以管道内空气流动为例,若管径0.1 m、平均流速10 m/s,按Re≈6.7×10⁴估算,最大涡尺度约为管径量级(10⁻¹ m),而最小耗散尺度η≈D·Re^(-3/4)约为10⁻⁵ m,两者相差四个数量级。若要在网格上直接解析全部尺度,三维网格量级将达到10¹²以上,远超当前工程算力。

因此,工程湍流模拟的核心矛盾不是"能不能算",而是"解析到什么尺度、把哪些尺度交给模型"。围绕这一取舍,形成了三条主要技术路线:RANS、LES与DES。RANS对时间做雷诺平均,把所有尺度的脉动都交给湍流模型(如k-ε、k-ω SST、Spalart-Allmaras),只求解平均场,计算量最小,适合稳态、高雷诺数、以平均量(压降、升阻力、换热系数)为目标的场景。LES对空间做滤波,只对大尺度涡直接求解,小尺度用亚格子模型(如Smagorinsky、WALE)封闭,能捕捉瞬态涡结构与宽频脉动,但近壁区网格要求极高(y⁺≈1),代价随Re迅速上升。DES则是二者的混合:在附面层内用RANS以节省网格,在分离区与尾流中切换为LES以解析非定常涡,是外流分离、气动噪声等问题的折中方案。

湍流模拟方法谱系与Fluent求解器定位 解析尺度越细 → 精度越高、代价越大;模型承担尺度越多 → 代价越低、适用面越窄 大尺度涡(含能区) 小尺度涡(耗散区) RANS 雷诺平均 · 全尺度建模 模型:k-ε / k-ω SST / SA 网格:y⁺ 30~300(壁面函数) 目标:稳态平均量、压降、升阻力 代价:★☆☆☆☆ LES 空间滤波 · 小尺度建模 模型:Smagorinsky / WALE 网格:y⁺≈1,Δ≤惯性子区 目标:瞬态涡结构、宽频脉动 代价:★★★★☆ DES / SBES RANS-LES 混合 · 分区切换 附面层:RANS(省网格) 分离/尾流:LES(解析涡) 目标:分离流、气动噪声 代价:★★★☆☆ Fluent 湍流建模整体工作流 ① 几何与网格 y⁺ 预估、边界层 ② 物理与模型 RANS/LES/DES 选型 ③ 边界条件 湍流强度/尺度 ④ 求解与收敛 残差、监测量 ⑤ 验证 网格/实验 验证不通过 → 返回网格或模型重新评估

在Fluent中,上述路线并非抽象概念,而是通过具体设置落地的。选择RANS时,通常在Viscous面板中启用k-ω SST,并配合壁面函数或增强壁面处理;选择LES时,需将瞬态求解器打开、时间步长取到涡翻转时间的一小部分(例如特征长度/特征速度的1/20~1/10),并保证近壁y⁺≈1;选择DES时,则需在相同面板中启用DES或SBES,并注意RANS-LES交界面附近的网格尺度不能过粗,否则会触发"网格诱导分离"(GIS)误差。

选型决策可以按以下顺序自问:第一,关心的是时均量还是瞬态脉动?若只关心压降、升阻力、换热系数等时均量,且流动无大范围分离,RANS是性价比最高的选择。第二,是否存在显著的非定常分离、涡脱落或宽频噪声?若是,且计算资源允许,优先考虑DES或SBES;若分离区尺度相对整体较小,DES的RANS附面层能显著节省网格。第三,是否必须解析小尺度湍流结构本身?如燃烧、气动噪声频谱、湍流混合细节,才值得付出LES的代价。第四,壁面处理方式是否与模型匹配?k-ε通常配标准壁面函数,k-ω SST可解析到壁面,LES必须解析到壁面,混用会导致系统性偏差。

工程经验提示:在正式计算前,先用一个粗网格RANS算例快速获得流场拓扑(分离区位置、回流区大小、驻点位置),再据此决定是否升级到DES/LES,以及加密区域应放在哪里。这一"先RANS后高保真"的策略,通常能避免把LES网格浪费在并不关键的区域。

需要强调的是,湍流模型不是越复杂越好。RANS的k-ε在自由剪切流中表现稳健,但在逆压梯度与强曲率流中容易高估湍动能;k-ω SST通过混合函数在近壁用k-ω、远场用k-ε,改善了分离预测,但对转捩与强非平衡效应仍有限。LES的亚格子模型在近壁区若网格不足,会低估壁面切应力。DES的RANS-LES切换依赖网格尺度与湍流尺度的比较,若交界面落在附面层内,会人为制造"对数层失配"。因此,模型选型的本质是让模型的假设与你的流动特征、网格能力、精度目标三者对齐,而不是追求方法谱系上的"高级"。

本章后续将围绕这一决策逻辑展开:先建立湍流尺度与建模谱系的定量认识,再逐一拆解Fluent中各模型的适用边界、关键参数与校准方法,最终形成一套可复用的选型与验证流程。