fluent大涡模拟

一、项目背景与科学问题

槽道湍流是流体力学中最基本的壁面湍流模型，广泛存在于管道输运、热交换器与航空发动机冷却通道中。湍流的多尺度特性使得直接数值模拟（DNS）代价极高，而雷诺平均（RANS）方法虽计算量小但难以捕捉大尺度涡结构的瞬态演化。大涡模拟（LES）通过在物理空间过滤湍流场，对大尺度涡直接求解、对小尺度涡用亚格子模型模拟，在计算代价与精度间取得平衡，已成为工业湍流研究的主流方法。

本项目采用 Fluent 对 Re_τ = 180 的充分发展槽道湍流进行大涡模拟，分析平均速度剖面、脉动速度场与近壁湍流结构，并与 Moser-Kim-Mansour 的 DNS 基准数据对比，评估 LES 的预测能力。

二、计算方法与理论背景

2.1 大涡模拟方程

LES 对不可压 N-S 方程施加空间滤波，得到滤波后的动量方程，其中亚格子应力 τ_sgs = uū − ūū 需建模。本项目采用 Smagorinsky-Lilly 模型（Cs = 0.1），并结合壁面阻尼函数（Van Driest）修正近壁行为。

2.2 槽道湍流与摩擦雷诺数

槽道流动由上下平行壁面驱动，充分发展后形成统计定常的湍流。摩擦雷诺数 Re_τ = u_τ·δ/ν，其中 u_τ 为摩擦速度，δ 为半槽宽。Re_τ = 180 是 LES/DNS 验证的经典工况。平均速度剖面用壁面单位表示为 u⁺ = f(y⁺)，其近壁区满足黏性底层 u⁺ = y⁺，对数区满足 u⁺ = (1/κ)ln y⁺ + B（κ = 0.41，B ≈ 5.0）。

2.3 数值设置

• 网格：展向与流向均匀网格，近壁第一层网格 y⁺ < 1。
• 时间推进：二阶隐式；压力-速度耦合：SIMPLEC。
• 边界条件：流向/展向周期性，壁面无滑移；恒定质量流率驱动。

三、计算设置与模型构建

• 软件：ANSYS Fluent；亚格子模型：Smagorinsky-Lilly。
• 计算域：2πδ × 2δ × πδ（流向×法向×展向），Re_τ = 180。
• 网格量：约 128 × 97 × 128。
• 时间步长：Δt⁺ = u_τ²Δt/ν ≈ 0.4；统计取样时长 ≥ 1000 个涡翻转时间。

四、结果分析与案例图解读

图 11 的左面板给出 LES 计算的平均速度剖面 u⁺(y⁺)。蓝色实线为 LES 结果，灰色虚线与点线分别对应黏性底层 u⁺ = y⁺ 与对数律 u⁺ = (1/0.41)ln y⁺ + 5.0。LES 曲线在 y⁺ < 5 区域与 u⁺ = y⁺ 完美重合，在 30 < y⁺ < 300 区域与对数律吻合，说明近壁网格与亚格子模型设置正确，湍流边界层结构被准确再现。这一剖面是论文中”LES 验证”章节的标准基准图，通常需与 DNS 数据叠加对比以证明可靠性。

图 11 的右面板给出瞬时流向脉动速度 u’ 的等值线云图（流向-法向平面）。红色与蓝色区域分别对应高速与低速流体条带，可以看到明显的近壁条带结构与由壁面向核心区延伸的大尺度相干涡。这些瞬态结构是 LES 相比 RANS 的核心优势——能够解析涡的时空演化，为研究湍流猝发、传热强化机制提供直接图像。

结合两面板，平均剖面证明统计结果的正确性，瞬态云图展示 LES 解析湍流结构的能力，二者共同验证了本项目 LES 设置对 Re_τ = 180 槽道湍流的适用性。

五、结论与展望

本项目在 Fluent 中完成了 Re_τ = 180 槽道湍流的大涡模拟，平均速度剖面与黏性底层、对数律定量吻合，瞬态场清晰展示了近壁条带与相干涡结构。案例图以统计剖面与瞬态云图双面板呈现，兼顾定量验证与机理展示。

后续可拓展：(1) 提高 Re_τ（如 395、590）并对比 DNS 数据，评估模型的雷诺数敏感性；(2) 加入温度方程，研究湍流传热与热条带结构；(3) 引入动态亚格子模型（Dynamic Smagorinsky），提升分离流工况的精度。