歧管流量分配
理论与物理
概述
老师！歧管流量分配分析，在什么场景下使用呢？
歧管（分叉管、集管）是一种将流体从一根主管均匀分配到多根分支管路的部件。用于燃料电池堆、散热器、锅炉水管群、冷却水套等流量均匀性直接影响性能的场景。
控制方程
支配流量分配的物理原理是什么？
歧管各分支的流量分配，由主管内的静压分布和各分支流路阻力的平衡决定。基础是伯努利方程和连续性方程。
流向各分支 i 的流量，由分支点静压与分支出口压力之差驱动。
$C_d$ 是流量系数吧。会随分支形状变化吗？
是的。直角分支约为 $C_d \approx 0.6$〜$0.8$，光滑的喇叭口分支约为 $C_d \approx 0.9$〜$0.98$。
流量均匀度的定量指标
介绍几个评估分支流量均匀性的指标。
| 指标 | 定义 | 理想值 | ||
|---|---|---|---|---|
| 流量均匀性指数 $\gamma$ | $1 - \frac{1}{2n\bar{Q}}\sum | Q_i - \bar{Q} | $ | 1.0 |
| 分配不均因子 | $\frac{Q_{max} - Q_{min}}{\bar{Q}}$ | 0 | ||
| 标准差 $\sigma_Q$ | $\sqrt{\frac{1}{n}\sum(Q_i - \bar{Q})^2}$ | 0 | ||
| 变异系数 CV | $\sigma_Q / \bar{Q}$ | 0 |
燃料电池需要多高的流量均匀度？
PEFC（质子交换膜燃料电池）堆栈期望 CV < 5%。超过10%则电池单元间的温度/反应不均会加剧，导致堆栈性能下降。
U型 vs. Z型歧管
请讲解U型和Z型的区别。
入口和出口在同一侧的是U型（逆流），在相反侧的是Z型（平行流）。
| 配置 | 流量分布趋势 | 均匀性 |
|---|---|---|
| U型 | 两端分支流量大，中央少 | 稍有不均 |
| Z型 | 出口侧分支流量大 | 容易不均 |
| 混合型 | 取决于设计 | 优化空间大 |
这种趋势可以用 Bajura & Jones (1976) 的理论模型解释。主管内的静压，一方面因摩擦损失而下降，另一方面因流量从分支流出导致流速减小，动压恢复（Static Regain）会使静压上升。这两种效应的竞争决定了静压分布。
Static Regain 的思路和风道设计一样呢。
歧管流动理论的起源——燃料电池开发催生的均匀分配理论
歧管流量均匀分配的理论研究在1990年代燃料电池（PEMFC）开发热潮中迅速深化。由于燃料电池堆栈各单元氢气/空气的均匀供应对性能和耐久性至关重要，歧管形状优化成为重要课题。在结合了伯努利方程和动量守恒的 Hardy & Collins(1954) 管网理论基础上，Bajura & Jones(1976) 整理了Z型、U型歧管分配不均的理论公式。该理论预测表明，等截面U型歧管在分支数为10时，末端通道与中央通道的流量差可达30%以上，这成为现代CFD优化方向的基准。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才变成稳定水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非定常现象。那么定常分析是什么？就是“经过足够时间流动稳定后”的状态——即令此项为零。计算成本大幅降低，因此先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶丢进河里会怎样？会被水流带着往下游漂，对吧？这就是“对流”——流体运动搬运物质的效果。暖风的暖气能送到房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——此项包含“速度×速度”，因此是非线性的。也就是说，流速加快此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水。因为蜂蜜粘度（$\mu$）高，不易流动。粘度越大扩散项越强，流体运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。反之，Re数大的流动中对流占压倒性优势，扩散成为配角。
- 压力项 $-\nabla p$：推注射器的活塞，液体会从针头有力地射出，对吧？为什么？因为活塞侧高压，针头侧低压——这个压差产生了推动流体的力。大坝放水也是同样原理。天气图中等压线密集的地方会怎样？没错，会刮强风。“有压差的地方产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD的“压力”多为表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压/表压。
- 源项 $S_\phi$：受热的空气会上升——为什么？因为比周围轻（密度低），被浮力推上去。这个浮力作为源项添加到方程中。还有，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——冬天房间里开了暖气但暖空气不上升，得到这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3 时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq 近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用情况：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等方法）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件从体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压和绝对压力。可压缩分析使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C、水: 约998 kg/m³@20°C |
| 粘性系数 $\mu$ | Pa·s | 注意与运动粘性系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数值方法的细节
请讲解歧管CFD的具体实现。
网格策略
歧管中分支处的分离和涡流对流量分配影响很大，因此分支处的网格质量尤其重要。
| 区域 | 网格尺寸 | 备注 |
|---|---|---|
| 主管直管段 | D/20〜D/10 | 壁面棱柱层5层以上 |
| 分支汇合部 | D/40〜D/20 | 解析分离区域 |
| 分支管入口 | d/20〜d/10 | 影响流量系数 |
| 主管端部（闭端/开端） | D/30 | 驻点的压力恢复 |
需要特别细化分支处呢。
是的。分支拐角处会发生分离，形成缩流断面。如果不解析这个，就会高估流量系数，导致各分支流量预测精度下降。
边界条件
典型的边界条件设置：
典型的边界条件设置：
| 边界 | 条件 | 备注 |
|---|---|---|
| 主管入口 | 质量流量入口 | 指定总流量 |
| 各分支出口 | 压力出口 | 相同压力（大气开放时） |
| 壁面 | 无滑移、绝热 | 多为光滑壁面假设 |
把各分支出口设为相同的压力出口，流量就会自然分配吗？
是的。将各分支的出口压力设为相同值（例如表压0 Pa），CFD就会基于静压分布和流路阻力自动计算各分支的流量。这是歧管CFD的基本方法。
但是，如果各分支下游有不同的压力损失元件（例如燃料电池单元、散热器芯体），则需要在分支出口设置额外的阻力（如多孔介质跳跃等）。
湍流模型的选择
推荐的湍流模型是？
因为分支处的分离很重要，所以推荐 SST k-omega。k-epsilon 系倾向于低估分支处分离泡的尺寸。
求解器设置
| 参数 | 推荐设置 |
|---|---|
| 求解器 | 基于压力、定常 |
| 压力-速度耦合 | 耦合式（分支多时注重鲁棒性） |
| 对流格式 | 二阶迎风 |
| 梯度 | 基于单元的最小二乘法 |
| 收敛判定 | 残差 1e-5 + 全部分支流量监控 |
推荐耦合求解器，是因为分支多时压力-速度耦合会变难吗？
正是如此。SIMPLE 系在分支数超过10个时收敛可能会变慢。耦合求解器虽然消耗更多内存，但收敛鲁棒性更高。
计算结果评估
计算后需检查以下内容：
1. 通过 Report > Fluxes 确认各分支的质量流量
2. 入口流量与全部分支流量总和一致（质量守恒）
3. 绘制主管内静压分布图
4. 通过分支处速度矢量确认分离模式
5. 计算流量均匀性指数
歧管流量分配的数值方法——压力损失定律的离散化与收敛稳定性
歧管流量分配CFD分析中，同时求解从主管分出的众多通道时的数值稳定性是个难题。由于压力损失的非线性（ΔP ∝ V²），分支数增加时联立方程的条件数容易恶化。实践中有效的方法是：① 将 Hardy-Cross 法（压力平衡环路迭代）的一维网络代码与全三维CFD结合；② 在三维CFD中采用“混合网格”，即精细解析分支处、粗化主管。另外，当分支处发生湍流到层流的转换时（Re=500〜2300的转换区），可能出现定常解不收敛的情况，此时需要进行非定常分析或应用转换湍流模型（γ-Reθ）。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必须使用。
中心差分（Central Differencing）
二阶精度，但 Pe数 > 2 时会发生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
なった
詳しく
報告