📚 课程目录第3章 Spalart-Allmaras模型

3. Spalart-Allmaras模型

Spalart-Allmaras（简称S-A）模型是工程外流计算中使用频率极高的一方程湍流模型。它不求解湍动能k和耗散率ε，而是直接求解一个与涡黏性直接相关的输运变量，因此在壁面附近和逆压梯度区域表现出较好的鲁棒性。理解S-A模型的关键，不是记住公式，而是理解它“只输运一个量、靠经验关系构造涡黏性”的设计哲学，以及这种设计带来的适用边界。

一、理论基础：单变量输运与涡黏性构造

S-A模型求解的输运变量通常记为 \(\tilde{\nu}\)，量纲为运动黏度（m²/s）。它与湍流运动黏度 \(\nu_t\) 之间通过一个代数关系连接：

ν_t = ν̃ · f_v1 ,  f_v1 = χ³ / (χ³ + C_v1³) ,  χ = ν̃ / ν

其中 ν 为分子运动黏度，C_v1 为常数（Fluent默认约7.1）。当 χ 较大（湍流充分发展）时 f_v1→1，ν_t≈ν̃；在壁面附近 χ→0，f_v1→0，ν_t被自然抑制。这一构造使S-A模型无需像k-ε那样额外引入壁面衰减函数，壁面处理更简洁。

输运方程本身包含对流项、扩散项、生成项和耗散项。生成项由涡量大小与平均应变共同决定，耗散项则包含壁面阻塞效应。S-A模型的一个显著特点是：它不依赖剪切应力输运的显式建模，而是通过涡量生成项在边界层内产生湍流黏度，因此对自由剪切流（如远场尾迹）的预测能力弱于对附着边界层的预测能力。

S-A模型属于“涡黏性模型”范畴，假设雷诺应力与平均应变率呈线性关系。它不适用于强各向异性湍流、强旋流或分离区内部精细结构预测，但在附着流和轻度分离的外流中性价比很高。

二、壁面处理：低雷诺数直接积分

S-A模型的设计初衷是“从壁面一直积分到边界层外缘”，因此它属于低雷诺数模型，不需要壁面函数。在Fluent中，使用S-A模型时通常要求壁面第一层网格的 y+ 控制在 1 量级（建议 y+ ≤ 1～2）。如果 y+ 过大，壁面附近的 ν̃ 梯度无法被正确解析，涡黏性会被高估或低估，导致摩阻和分离位置偏差。

与k-ε模型不同，S-A模型在壁面处 ν̃=0，且方程本身包含壁面阻塞项，因此不需要额外的壁面函数修正。但在Fluent中，如果网格 y+ 无法做到1，可以启用“壁面函数”选项吗？答案是否定的——S-A模型没有标准壁面函数。若必须使用粗网格，可考虑改用k-ε或SST模型，而不是强行给S-A加壁面函数。

实际工程中，若翼型表面 y+ 控制在0.5～1.5，S-A模型对升力系数和阻力系数的预测通常与实验偏差在3%～8%以内（示例范围，具体取决于网格和来流条件）。若 y+ 超过5，阻力预测可能偏大10%以上，分离位置也会提前。

三、适用场景：外流、翼型、涡轮机械

S-A模型最擅长的领域是空气动力学外流：

翼型与机翼：附着流、轻度逆压梯度、激波-边界层干扰（需配合转捩模型或固定转捩位置）。
涡轮机械：压气机叶片、涡轮叶片外部流道，尤其是叶片表面边界层和尾迹的初步评估。但叶顶间隙泄漏、强旋流区域需谨慎。
导弹、无人机、汽车外流：以附着流为主、分离区不大的场景。
不适用场景：强分离、钝体尾迹、射流、旋流燃烧器、浮力驱动流、多相流中的连续相湍流。

在涡轮机械中，S-A模型常被用于叶栅通道的初步设计迭代，因为其计算量小、收敛快。但若需要精确预测叶顶泄漏涡或端壁二次流，应改用SST k-ω或雷诺应力模型。

四、Fluent设置要点与操作示例

在Fluent中启用S-A模型：Models → Viscous → Spalart-Allmaras (1 eqn)。关键设置包括：

Production Limiter：建议开启，防止驻点区域湍流黏性过度生成。
Curvature Correction：若存在明显流线弯曲（如涡轮叶片），可开启曲率修正。
Options：默认使用“Vorticity-Based”生成项；若流动以应变为主，可切换为“Strain-Based”。

边界条件设置：

Inlet:  ν̃ = 3～5 × ν_∞  (来流湍流黏度比)
Outlet:  ν̃ = 0 或按自由流外推
Wall:    ν̃ = 0  (无滑移壁面)
Symmetry: 零法向梯度

初始化时，建议先用混合初始化，再对 ν̃ 场进行patch：Solution → Initialize → Patch → ν̃ = 3ν_∞。若直接使用默认零场，收敛可能变慢。

求解设置：压力-速度耦合推荐Coupled，压力离散用Second Order，动量用Second Order Upwind，ν̃ 方程用First Order Upwind起步，收敛后再换Second Order。松弛因子：ν̃ 方程默认0.7～0.8，若发散可降至0.5。

S-A模型对网格质量敏感。若壁面第一层网格长宽比过大（>1000）或存在负体积，容易导致 ν̃ 场出现非物理负值。Fluent会限制 ν̃ ≥ 0，但负值截断会破坏守恒性，建议检查网格正交质量 > 0.2。

五、结果判读与校准示例

以NACA 0012翼型为例（来流马赫数0.15，雷诺数6×10⁶，攻角4°），使用S-A模型计算。若 y+ 控制在1左右，升力系数预测约0.44，阻力系数约0.008（示例值，非实测）。若将 y+ 放大到8，阻力系数可能升至0.0095，升力系数变化不大但分离点前移。判断依据：监测翼型表面摩阻系数分布，若前缘附近摩阻异常偏高，说明 y+ 过大或生成项过强。

校准建议：若已知实验升阻比，可微调来流 ν̃/ν 比值（通常3～10），但该参数对附着流影响较小，对分离流影响较大。更有效的校准手段是转捩位置设定——S-A模型本身不预测转捩，需配合γ-Reθ转捩模型或手动固定转捩点。

总结来说，S-A模型是外流工程中“快而稳”的选择，但它的能力边界清晰：附着流和轻度分离是它的主场，强分离和旋流是它的短板。使用时的核心动作是控制 y+、开启生成项限制、合理设置来流 ν̃/ν，并在结果判读时重点关注摩阻分布和分离位置，而不是只看升力系数。