3、CFD基础：计算流体力学基本方程、湍流模型选择、边界条件设置
各位工程师朋友，咱们今天聊聊CFD的基础。说实话，很多人一上来就怼软件，点鼠标点得飞起，结果算出来的东西跟实测差十万八千里。为什么？因为底层的物理没搞明白。你想想看，你连方程长什么样都不知道，怎么判断收敛曲线对不对？
我做了十几年热仿真，踩过的坑比你们见过的网格都多。今天我就把CFD最核心的三个东西——基本方程、湍流模型、边界条件——掰开了揉碎了讲清楚。别怕，咱们不搞数学推导，只讲工程直觉。
核心观点：CFD不是黑箱，是物理。你理解了多少物理，仿真就有多准。
3.1 三大守恒方程：流体的“宪法”
CFD求解器到底在算什么？说白了，就是在解三个方程：质量守恒、动量守恒、能量守恒。这三个方程是流体力学的“宪法”，所有流动现象都得遵守。
质量守恒方程（连续性方程）
这个方程讲的是：流入控制体的质量，等于流出控制体的质量，加上内部质量的变化。用大白话说就是——流体不会凭空消失，也不会凭空产生。
∂ρ/∂t + ∇·(ρU) = 0对于不可压缩流动（比如水，或者低速空气），密度ρ是常数，方程简化为：
∇·U = 0我在做储能风道设计时，经常用这个方程来检查入口和出口的流量是否匹配。如果入口流量和出口流量对不上，那肯定是边界条件设错了，或者网格有问题。
动量守恒方程（Navier-Stokes方程）
这个方程是CFD的“灵魂”。它描述的是：流体微团的动量变化率，等于作用在它上面的各种力之和——压力、粘性力、体积力（比如重力）。
∂(ρU)/∂t + ∇·(ρUU) = -∇p + ∇·τ + ρg嗯，这里要注意。左边第二项∇·(ρUU)是对流项，它是方程非线性的根源。为什么CFD这么难收敛？就是因为它。我刚开始做仿真时，经常遇到残差震荡，后来发现是对流项离散格式选得太激进。
能量守恒方程
对于储能热管理，这个方程至关重要。它描述的是：流体携带的热量，加上热传导、热对流、热辐射的贡献。
∂(ρh)/∂t + ∇·(ρUh) = ∇·(k∇T) + S_h其中S_h是源项，比如电池的发热量。我在做电池包热仿真时，最头疼的就是这个源项怎么给。给大了算出来温度过高，给小了又太乐观。后来我学乖了，直接用实测的发热曲线拟合进去。
我的经验：对于储能系统，能量方程里的粘性耗散项通常可以忽略。但如果你在做高速流动（比如风道入口风速超过20m/s），那就要小心了，粘性耗散产生的热量不可忽视。
3.2 湍流模型选择：别迷信“高级”
说到湍流，很多新手就头大。其实你只要记住一句话：没有完美的湍流模型，只有合适的湍流模型。
为什么会这样？因为NS方程本身是精确的，但直接求解（DNS）需要超级计算机，工程上根本用不起。所以我们要用雷诺平均（RANS）或者大涡模拟（LES）来近似。不同的近似方法，就产生了不同的湍流模型。
常用湍流模型对比
| 模型 | 适用场景 | 计算量 | 精度 |
|---|---|---|---|
| 标准k-ε | 高雷诺数、远离壁面 | 低 | 中等 |
| RNG k-ε | 中等旋流、分离流 | 低 | 中等偏高 |
| Realizable k-ε | 射流、旋转流 | 低 | 高 |
| SST k-ω | 壁面边界层、分离流 | 中等 | 高 |
| LES | 大尺度涡、瞬态细节 | 极高 | 很高 |
我个人习惯，做储能风道优化时首选SST k-ω模型。为什么？因为风道里有大量的壁面边界层和局部分离区，SST k-ω对壁面附近的处理比k-ε好得多。我曾经用标准k-ε算一个风道弯头，结果压损比实测低了30%，换成SST k-ω后误差缩小到5%以内。
避坑指南：我曾经犯过一个错误——在低雷诺数区域用了标准k-ε模型。结果算出来的湍流粘度异常大，流动几乎变成了层流。后来才发现，标准k-ε要求y+在30-300之间，而我当时的网格y+只有5。记住：模型选对了，网格也得跟上。
怎么选？给你三个原则：
- 看雷诺数：储能风道一般Re在10^4-10^5，属于中等湍流。SST k-ω或Realizable k-ε都行。
- 看几何复杂度：如果风道里有大量弯头、变截面、障碍物，用SST k-ω。如果只是直管道，标准k-ε就够了。
- 看计算资源：LES虽然准，但一个储能集装箱的LES仿真，用32核工作站也得跑一周。你老板等得起吗？
3.3 边界条件设置：仿真准不准，全看这一步
我经常跟团队说：边界条件设对了，仿真就成功了一半。另一半是网格质量。边界条件设错了，求解器再牛也白搭。
入口边界条件
最常见的是速度入口和质量流量入口。对于储能风道，我建议用质量流量入口。为什么？因为风机的特性曲线通常给出的是流量-压头关系，而不是速度。用质量流量入口，可以更准确地匹配风机性能。
# 质量流量入口设置示例（Fluent TUI）
/define/boundary-conditions/inlet
mass-flow-inlet
mass-flow-rate 2.5 ; 单位 kg/s
direction normal-to-boundary
turbulence-intensity 5%
turbulence-viscosity-ratio 10湍流强度怎么设？对于储能风道，入口湍流强度一般在5%-10%之间。如果你不确定，设5%是个安全的起点。我见过有人设0.1%，结果算出来的流动几乎层流化，完全不对。
出口边界条件
出口通常用压力出口。注意：压力出口需要给定静压，不是总压。对于开口到大气的情况，静压设0 Pa（表压）即可。
这里有个坑：如果出口有回流，压力出口的收敛性会变差。我曾经算一个风道，出口有30%的回流，残差死活降不下去。后来我把出口延长了3倍管径，回流消失了，收敛也好了。
壁面边界条件
壁面分两种：绝热壁面和恒温壁面。对于储能系统，电池表面通常用恒热流边界条件（因为电池发热量已知），风道壁面用绝热或对流换热边界条件。
壁面函数的选择也很关键。如果你用SST k-ω模型，建议用增强壁面处理（Enhanced Wall Treatment），这样对y+的要求不那么苛刻，y+在1-5之间都能接受。
我的小技巧：在做风道优化时，我会在入口和出口各加一段延长段。入口延长5倍水力直径，出口延长10倍。这样能消除入口效应和出口回流对计算域的影响，结果更稳定。
3.4 知识体系总览
说了这么多，咱们用一张图把CFD基础的核心逻辑串起来。你照着这个框架去理解，就不会迷失在细节里。
这张图把CFD基础分成了三层：物理层（三大方程）、模型层（湍流模型）、应用层（边界条件）。你从底层往上走，每一步都踩实了，仿真结果自然靠谱。
好了，关于CFD基础就聊这么多。记住：方程是根，模型是枝，边界条件是叶。根深才能叶茂。下次你打开仿真软件时，先问问自己：这三个东西，我搞明白了吗？