8.2 湍流模型从RANS到LES

8.2 湍流模型从RANS到LES

本节摘要：第 6 章看到湍流把能量从大涡传给小涡，直至粘性耗散——这意味着直接求解要解析所有尺度，代价随雷诺数爆炸。工程的三条路按算力定价：RANS 时均建模（便宜、工程日用）、LES 大涡模拟（解析大涡、建模小涡）、DNS 直接求解（学术基准、天价）。本节讲清雷诺应力封闭问题、k-epsilon 与 k-omega-SST 两方程模型的分工、壁面函数与 Y+ 的配合，最后算一笔成本账帮你选型。

飞机巡航雷诺数上亿，含能涡的尺度比最小耗散尺度大五六个数量级——DNS 要解析全部尺度，网格量近似按 Re 的三次方增长，工程时限内毫无指望。出路只有一条：不追每个涡的细节，只统计或只解析大涡。湍流模型的层级就是"你愿意花多少算力、换多少细节"的价目表。

一、封闭问题：时均的代价

把速度做雷诺分解 u = ū + u'（第 6.2 节的时均加脉动），代入 N-S 方程再平均，多出一项雷诺应力 −ρ⟨u'ᵢu'ⱼ⟩——脉动之间的动量输运，形状与粘性应力一样但来源完全不同。方程未知数多于方程数，必须补充模型让账本闭合，这就是"封闭问题"，湍流建模一百年来的主战场。

涡粘假设是应用最广的 Closure：雷诺应力类比特粘性，写成 μ_t 乘以应变率，湍流被简化成"粘性变大了的层流"。μ_t 不是物性，随流动变化，于是有了求它的两方程模型。

二、两方程模型与壁面处理

k-epsilon 家族解两个输运方程：湍动能 k 与耗散率 ε，μ_t = ρC_μk²/ε。它远场表现好、近壁失灵；k-omega 家族把第二个变量换成比耗散率 ω，近壁稳健但远场敏感。k-omega-SST 把两者缝合：近壁用 k-omega、远场切 k-epsilon，并限制涡粘防止分离被高估——航空外流、汽车风阻的默认选择。

壁面附近的配合规则浓缩成一个数：Y+（第一层网格中心的无量纲距离）。壁面函数路线要求 Y+ 在 30 到 300 之间，把粘性底层用对数律"代答"；低雷诺数直接解析路线要求 Y+ 接近 1，网格量陡增。画网格前先估 Y+：

import numpy as np

# 平板估计第一层网格高度 目标 Y+ 依模型路线而定
rho, U, L, mu = 1.225, 50.0, 3.0, 1.8e-5      # 来流条件 小型无人机尺度
Re_L = rho*U*L/mu
Cf = 0.026 / Re_L**0.139                      # 经验摩擦系数
tau_w = Cf * 0.5 * rho * U**2
u_tau = np.sqrt(tau_w / rho)                  # 摩擦速度
y_target = 50.0                               # 壁面函数路线 Y+ 目标 50
y1 = y_target * mu / (rho * u_tau)
print(f"Re_L = {Re_L:.2e}  u_tau = {u_tau:.2f} m/s")
print(f"第一层网格高度约 {y1*1000:.3f} mm (壁面函数 Y+ 50)")
print(f"若改低雷诺数解析 Y+ 1 则需 {y1/50*1000:.5f} mm 网格量成倍上涨")
​

三、层级价目表：RANS、LES、DNS

四、成本账演练

用估算式感受三级价差（网格量按解析尺度的幂次缩放）：

Re = 1e6                       # 低雷诺数算例 以 RANS 网格 200 万为基准
n_rans = 2e6
n_les  = n_rans * 12           # 大涡解析要求 粗略一个量级上浮
n_dns  = n_rans * Re**0.9 / 1e4 # Kolmogorov 尺度随 Re 增长 经验缩放演示
for name, n in [("RANS", n_rans), ("LES", n_les), ("DNS", n_dns)]:
    print(f"{name:5s} 网格量约 {n:.1e}")
​

同一段流动，从两百万网格到上亿甚至数十亿——这不是精度升级，是物理内涵的升级：RANS 只有时均场，LES 有真实的旋涡脱落与非定常力谱，DNS 连最小耗散涡的生死都看得见。预算决定你能回答哪一层问题。

⚠️ 报告合规提醒：写明模型名称、网格量、Y+ 范围与验证依据，"用了 CFD"四个字不构成任何证据。

变式思考

混合 RANS-LES（如 DES 系列）在壁面用 RANS 省网格、远离壁面切 LES 捕捉分离，正在接管大分离外流工程；而数据驱动的机器学习封闭模型试图从 DNS 数据库学出更准的应力-应变关系。方向在变，底层逻辑不变：封闭问题的答案永远在"你要的统计量"与"你付的算力"之间找平衡。

追问：模型误差到底有多大、怎么报

涡粘模型在三类流动里系统性失准，报告里要主动交代而非回避：

• 强分离与回流：时均涡粘假设把雷诺应力方向"锁死"在应变方向上，分离点普遍预报偏晚、回流区偏短；旋涡脱落频率可差两成。对策是换雷诺应力模型或上 LES，至少在报告中注明不确定度区间。
• 曲率与旋转：管道二次流、弯曲通道里标准 k-epsilon 误差显著，SST 与雷诺应力模型明显更好；旋转机械还要在源项里补旋转修正。
• 转捩位置：全湍流模型默认处处湍流，层流段摩擦被高估、升力线斜率略偏；低雷诺数机翼、涡轮叶片要开转捩模型（如 gamma-Re theta）。

# 报告不确定度的朴素做法 敏感性矩阵
cases = {
    "SST 全湍流": 0.315,
    "SST 加转捩模型": 0.298,
    "雷诺应力模型": 0.307,
}
base = cases["SST 加转捩模型"]
spread = (max(cases.values()) - min(cases.values())) / 2
print(f"设计阻力系数取 {base:.3f} 模型级不确定度正负 {spread:.3f}")
print("网格级不确定度另算 两者平方和开根 才是总不确定度")
​

把模型散布当成误差条写进报告，比给出一个假装精确的四位小数字专业得多；评审看的就是你对自己工具边界的自知之明。

本节要点回顾

• 封闭问题：雷诺分解多出雷诺应力，必须补模型闭合方程组；
• 涡粘假设加两方程模型是工程主力，SST 缝合了近壁与远场；
• Y+ 是壁面网格与模型的对接协议，30 到 300 或接近 1，画网格前先估；
• RANS/LES/DNS 是价目表不是排行榜，按需选型；
• 结果可信度三件套：模型名、网格量、验证依据。

模型选好，下一章末节把方法带进三个工程现场验收。