通道流DNS — 故障排除指南
本文已并入整合版
更充实的内容请参阅 channel-flow-dns.html。
更充实的内容请参阅 channel-flow-dns.html。
通道流DNS — 故障排除指南
故障排除
通道流DNS(或LES验证)中的常见问题是什么?
让我们来看看典型的问题。
统计量与DNS参考值不匹配
原因1: 分辨率不足
- $\Delta x^+$, $\Delta z^+$ 过大,无法解析条纹结构
- 确认: 计算能量谱,检查高波数侧是否存在混叠能量堆积
原因2: 采样时间不足
- 湍流统计量(特别是 $\overline{v'^2}$, $\overline{u'v'}$)的收敛不充分
- 对策: 将统计采样时间延长到 $Tu_\tau/\delta > 30$
原因3: 计算域大小不足
- 流向·跨向过短,最大涡结构被切割
- 确认: 二点相关函数是否在计算域内一半处衰减至零
可以通过二点相关确认计算域的充分性啊。
$Re_\tau$ 与设计值不同
在固定体积力驱动情况下,$Re_\tau$ 是输入参数,但在固定流量驱动情况下,实现的 $Re_\tau$ 事先无法知晓。
- 对策: 助跑计算中监测 $\tau_w$,稳定后确认 $Re_\tau$。如需要可调整流量
LES中壁面附近的精度
用LES计算通道流时,壁面附近的速度分布出现偏差。
SGS模型可能在壁面附近给出过度的耗散。对策:
- WALE或动态Smagorinsky: 使用在壁面附近能自动减少SGS粘性的模型
- 网格分辨率: 确保 $\Delta y^+_{wall} < 2$
- van Driest阻尼: 如果使用静态Smagorinsky,应用van Driest阻尼函数
检查清单
| 检查项目 | 基准 |
|---|---|
| 壁面第一层单元 $\Delta y^+$ | < 1 (DNS), < 2 (wall-resolved LES) |
| 流向分辨率 $\Delta x^+$ | < 10 (DNS), < 40 (LES) |
| 跨向分辨率 $\Delta z^+$ | < 5 (DNS), < 20 (LES) |
| 计算域 $L_x/\delta$ | > $4\pi$ (低Re), > $8\pi$ (高Re) |
| 计算域 $L_z/\delta$ | > $2\pi$ |
| 采样时间 $Tu_\tau/\delta$ | > 20 |
有数值基准,自检就很方便了。
正因为通道流的边界条件简单,分辨率和计算域大小的问题就被纯粹地暴露出来了。在这里掌握基础知识,今后应用到复杂形状的LES也就不难了。
咖啡小憩 闲聊天地
"湍流消失变成层流了"——统计不收敛的陷阱
通道DNS中常见的故障是"计算过程中湍流结构消失,退化成层流"的现象。大多数原因是"初始扰动不足"或"计算域太小,无法容纳维持湍流所需的涡尺度"。湍流通过"自我维持循环"(Hamilton等人的理论)来持续,如果周期长度太短,这个循环就无法建立。通常需要流向长度为8π×壁间距,缩短这个长度就会出现"看似收敛,其实是假的层流解"。如果用"速度场停止波动"来判断收敛,就容易漏掉,所以要用湍流统计量(如u'v'相关)的时间序列来确认。
相关主题
本文评价
感谢您的反馈!
有参考
价值
价值
需要
更详细
更详细
报告
错误
错误