大涡模拟 (LES) — CAE术语解说
LES是什么 — 空间滤波的思路
好问题。RANS是对湍流的「全部尺度」进行时间平均后模型化。也就是说把所有涡的信息都压入模型中。而LES是大涡（携带能量的尺度）用Navier-Stokes方程直接计算，仅对网格无法解析的小涡进行SGS模型近似。所以LES是「大涡直接计算，小涡仅模型化」的混合思路。
怎样区分「大涡」和「小涡」呢？
这就是空间滤波。对Navier-Stokes方程用滤波函数 $G$ 进行卷积积分，只提取大于网格宽度 $\Delta$ 的分量。数学上写成这样：
滤波后的NS方程中出现SGS应力张量 $\tau_{ij}^{\mathrm{sgs}} = \overline{u_i u_j} - \bar{u}_i \bar{u}_j$。这变成了未知量，需要用某种模型来闭合。实际上网格本身通常承担滤波的作用，这叫「隐式滤波（implicit filtering）」。
SGS模型（Smagorinsky·WALE·动态模型）
SGS应力的模型化方法有哪些种类？
代表的有3种，给你介绍一下。
1. Smagorinsky模型（1963年）
最经典的SGS模型。定义SGS涡粘度为：
$C_s$ 是Smagorinsky常数，对于各向同性湍流约0.17，通道流约0.1。简单稳定，但壁面附近涡粘度过大。通常用van Driest阻尼函数 $f_d = 1 - \exp(-y^+/A^+)$ 进行修正。
壁面附近涡粘度太大是指湍流被实际过度扩散了吧？
完全正确。壁面附近湍流涡小且弱，但Smagorinsky模型因为 $|\bar{S}|$ 大（速度梯度急陡）就增大涡粘度，这是问题所在。下面的模型改进了这个。
2. WALE模型（Wall-Adapting Local Eddy-viscosity, Nicoud & Ducros 1999）
用速度梯度张量的无迹对称部分：
这里 $\mathcal{S}_{ij}^d$ 是速度梯度张量平方的无迹对称部分。$C_w \approx 0.325$ 是标准值。WALE的优点是壁面附近自然满足 $\nu_{\mathrm{sgs}} \to 0$，不需要van Driest阻尼函数。在层流区也能正确降为零，很方便。
第三种动态模型是什么？
3. 动态Smagorinsky模型（Germano et al. 1991, Lilly 1992）
不把 $C_s$ 当作常数，而在计算中自动调整。用网格滤波 $\Delta$ 和测试滤波 $\hat{\Delta}$（通常 $\hat{\Delta} = 2\Delta$）进行两层滤波，从Germano恒等式
局部求出 $C_s^2$。模型常数随流场变化，在分离流、尾流、壁面附近都能取合适值。但计算容易不稳定，需要空间平均或 $C_s^2$ 不能为负的裁剪处理。OpenFOAM中有 dynamicKEqn 或 dynamicLagrangian 实现。
LES的网格要求和计算成本
做LES时网格要多细？ 和RANS完全不一样吧？
完全不同。LES需要直接解析涡，所以壁面附近的网格要求严格得多。大致标准是这样：
- 壁面垂直方向：$y^+ \approx 1$（第一层单元）
- 流向：$\Delta x^+ \approx 50 \sim 100$
- 展向：$\Delta z^+ \approx 15 \sim 40$
比如汽车外形空气动力学中 $\text{Re}_L = 10^6$ 级别，RANS用几百万网格就够，但壁面解析LES（WRLES）要几亿网格。而且是三维非稳态计算，时间步长也要保持 $\text{CFL} \le 1$，所以计算成本是RANS的100～1000倍。
1000倍？还有值得做的价值吗？
有的。RANS只能看到「时间平均」的流场，LES能看到涡的动力学。比如汽车后方湍流区的非稳态压力波动、建筑物周边风荷载的峰值、燃烧室内火焰与旋流的相互作用——这些「非稳态性本质重要的现象」RANS原理上无法捕捉。喷气发动机燃烧CFD已经几乎都用LES了，风工程也在快速推广。
LES vs RANS vs DNS — 使用判断标准
还有DNS（直接数值模拟）吧。LES和DNS怎么区别？
DNS是对湍流的全部尺度——从最大涡到Kolmogorov尺度 $\eta = (\nu^3/\varepsilon)^{1/4}$ ——都用网格解析。完全不用模型。所以最精确，但网格数按 $N \propto \text{Re}^{9/4}$ 增长，工业高Reynolds数流实际上用不了。
整理如下：
| 方法 | 模型化范围 | 网格数尺度律 | 实用Re上限 |
|---|---|---|---|
| DNS | 无（全尺度解析） | $\sim \text{Re}^{9/4}$ | $\sim 10^4$ |
| LES | 仅SGS尺度 | $\sim \text{Re}^{13/7}$（壁面解析） | $\sim 10^6$ |
| RANS | 全湍流尺度 | Re依赖小 | 无限制 |
实务中怎么区分选用？ 什么时候选LES？
判断标准简化后是这样：
- RANS就够：定常流、附着流的阻力/升力预测、管道内压降计算等。成本最便宜。
- 需要LES：大规模分离、尾流非稳态变动、混合·燃烧、空气动力学噪声（Aeroacoustics）、风工程突风响应等。
- 需要DNS：湍流基础研究、模型验证数据生成、过渡现象详细解明。主要学术目的。
实务上「先用RANS试试，精度不够再切到LES」这个阶梯式做法更现实。一开始就上LES的话计算成本读不准，会很困扰。
DES/DDES — LES和RANS的混合方法
壁面网格那么麻烦的话，壁面附近用RANS，远处用LES的方法行不行？
正是这个想法产生了DES（分离涡模拟）。Spalart（1997年）提出的方法，壁面附近表现为RANS模型（SA、SST等），离开壁面的区域则LES式解析涡。切换由网格尺寸和壁面距离自动判定。
但初期DES有问题。叫「灰色区问题」，RANS和LES区域的边界上，网格中等细密时，RANS涡粘度下降不足，LES又解析不了涡，变成「中不溜秋」的区域。
有改进的办法吗？
改进版是DDES（延迟DES）和IDDES（改进DDES）。DDES是Spalart等（2006年）提出的，边界层内部引入屏蔽函数 $f_d$ 强制RANS模式：
这样能防止附着流边界层意外切到LES模式的问题（MSD：建模应力耗尽）。IDDES还集成了壁面模型LES（WMLES）功能，壁面网格足够细的话就用WMLES运行。
OpenFOAM中有 SpalartAllmarasDDES 或 kOmegaSSTDDES 可直接用。FLUENT也支持DDES。实务上「先用RANS → 精度不足再用DDES → 还是不行就用壁面解析LES」这个逐级上升是王道。
实务中的LES运用建议
实务中用LES要注意什么重点？
经常被问到，重要的整理如下。
- 初值和助走时间：LES是非稳态计算，需要「助走时间（wash-out time）」让统计达到定常。从流通过计算域数次的时间（流通时间）的5～10倍开始统计。
- 入口边界条件：湍流入口条件极其重要。均匀流进入的话等边界层发育要很长的助走区。用Synthetic Eddy Method（SEM）或Recycling/Rescaling人工给湍流脉动效率更高。
- 数值格式：中心差分系为基础。风上差分的数值扩散大，会压平涡。但中心差分容易不稳定，通常混合少量（5～10%）风上分量或用低扩散格式（LUST、Gamma等）。
- 时间步长：Courant数保持在0.5～1.0以下。即使隐格法也是CFL > 5的话时间精度下降，LES就没意义了。
- 统计处理：计算平均值·RMS·功率谱密度（PSD），确认结果不依赖时间窗口。统计不足的话「恰好的快照」就变成了结论，风险很大。
非常清楚了。总结一下，LES就是「成本高但能看到涡的动力学，非稳态现象必须用」对吧。
完美总结。再补充一点，GPU计算的进步让LES成本年年下降。NVIDIA的AmgX、PETSc的GPU对应、基于Lattice Boltzmann的GPU求解器（PowerFLOW、XFlow、ProLB等）也出来了。汽车OEM的空力开发已经几乎都用LES系手法了。再过5年说不定LES比RANS用得更广的领域会增加。
CAE术语的准确理解是团队沟通的基础。— Project NovaSolver也关注实务者的学习支持。
分享您在LES实务中的课题
Project NovaSolver致力于解决CAE工程师日常面对的难题——配置的复杂性、计算成本、结果解读。您的实务经验将成为开发更好工具的驱动力。
咨询（准备中）相关主题
帮助
入讲解
错误