3. 湍流模型基础

外流场仿真中，绝大多数工程问题都处于湍流状态。机翼绕流、汽车尾流、建筑风载，其流动都表现出强烈的非定常、多尺度与随机性。湍流模型的选择直接决定计算成本与结果可信度，因此理解各类模型的物理假设与适用边界，是CFD工程师的核心能力。

湍流的物理特征

湍流最本质的特征是涡的级串过程：大尺度涡从平均流中抽取动能，通过惯性作用逐级破碎为小尺度涡，最终在Kolmogorov尺度上被黏性耗散为热能。这一过程可用能量级串图描述，大涡含能、小涡耗散，中间存在惯性子区。湍流还具有扩散性、耗散性和连续性，脉动速度使动量、热量与组分在流场中快速混合。

描述湍流需要区分三个速度尺度：平均速度、脉动速度和瞬时速度。雷诺分解将瞬时量写为时均量与脉动量之和，例如 u = ū + u'。这一分解是后续所有统计湍流模型的基础。湍流强度 I = u'/ū 是外流场入口边界的关键参数，自由来流中通常取0.1%～1%，风洞试验段约0.5%，大气边界层可达10%以上。

雷诺平均Navier-Stokes方程

将雷诺分解代入不可压缩N-S方程并取时间平均，得到RANS方程。动量方程中多出一项 -ρ·(u'ᵢu'ⱼ)̄，即雷诺应力张量。它由9个分量组成，因对称性实际有6个独立分量。RANS方程本身不封闭：未知量多于方程数，必须引入湍流模型来建立雷诺应力与平均速度场的关系，这就是湍流封闭问题。

所有工程湍流模型都在回答同一个问题：如何用已知的平均量表达未知的雷诺应力。不同模型对这一问题给出不同层次的近似，由此形成从涡粘模型到雷诺应力模型再到尺度解析模拟的谱系。

涡粘模型：Boussinesq假设下的工程主力

涡粘模型的核心是Boussinesq假设：将雷诺应力类比为黏性应力，用涡粘系数 μₜ 与平均应变率张量关联。这样封闭问题就转化为如何计算 μₜ。根据求解附加输运方程的数量，涡粘模型分为零方程、一方程和二方程模型。

Spalart-Allmaras（SA）是一方程模型的代表，直接求解修正涡粘量 ν̃ 的输运方程。它计算量小、鲁棒性好，对附着边界层和中等分离流精度可靠，是航空外流场（如机翼、翼身组合体）的常用选择。但SA对自由剪切流和强分离的预测能力有限。

k-ε模型求解湍动能 k 与耗散率 ε 两个方程，是工业界应用最广的模型。标准k-ε在自由剪切流中表现良好，但近壁区需要壁面函数处理，且对逆压梯度分离流预测偏晚。RNG k-ε和Realizable k-ε通过修正系数改善了旋流和强应变适应性。

k-ω SST模型由Menter提出，巧妙地在近壁区使用k-ω模型（无需壁面函数，能捕捉边界层细节），在远场切换为k-ε模型（避免k-ω对自由来流敏感的问题），并通过剪切应力输运限制器防止涡粘过度预测。对于存在逆压梯度、流动分离的外流场问题，如汽车尾部、翼型大攻角，SST通常是首选。

雷诺应力模型（RSM）

涡粘模型的根本局限在于Boussinesq假设的各向同性涡粘假设，无法描述雷诺应力的各向异性。RSM直接求解6个雷诺应力分量和耗散率的输运方程，不依赖涡粘假设。它能更准确地捕捉强旋流、浮力流和曲率效应，但方程数量多、计算成本高、收敛难度大，工程外流场中多用于特定验证场景。

大涡模拟（LES）简介

LES采用空间滤波将湍流场分为大尺度（直接求解）和小尺度（模化）两部分。大涡携带主要能量且依赖几何边界，小涡近似各向同性，更易模化。亚格子应力模型（如Smagorinsky、WALE）负责封闭滤波后的方程。LES能捕捉非定常涡结构和瞬态分离，但要求近壁 y⁺≈1、网格量巨大、时间步受CFL限制，工程外流场中常采用DES/IDDES等混合方法，在近壁用RANS、分离区用LES，兼顾精度与成本。

工程实践中的关键参数

以某翼型绕流为例，来流速度30 m/s，弦长1 m，空气运动黏度1.5×10⁻⁵ m²/s，雷诺数约2×10⁶。若采用SST模型，入口湍流强度取0.5%，湍流长度尺度取弦长的0.07倍，即0.07 m。由 k = 1.5(I·U)² 得 k≈0.034 m²/s²，由 ω = k^0.5 / (Cμ^0.25 · L) 得 ω≈180 1/s。近壁首层网格高度按 y⁺≈1 估算，约 2×10⁻⁵ m。这些参数直接决定边界层解析质量与分离预测精度。

判断结果是否可信，应检查：残差是否降至10⁻⁴以下并平稳、关键监测量（升力系数、阻力系数）是否收敛、y⁺分布是否满足所选模型要求、网格无关性是否通过。若分离区预测与试验偏差大，可尝试SST替代SA，或加密尾缘网格，而非盲目更换模型。