混合層
理论与物理
概述
老师，混合层是什么？
当两个具有不同速度的平行流动汇合时，在其界面上形成的剪切层就是混合层（mixing layer）。它是最简单的自由剪切流之一，也是 Kelvin-Helmholtz 不稳定性的典型例子。
它会在什么场景中出现呢？
例如航空发动机排气与旁通流的界面、河流汇合处、大气中的锋面、从喷嘴喷出的射流外缘部等。它是直接影响混合与输运效率的重要流动。
基本参数
设两个流动的速度为 $U_1$（高速侧）和 $U_2$（低速侧），重要的参数如下。
- 速度比: $r = U_2 / U_1$
- 速度差: $\Delta U = U_1 - U_2$
- 对流速度: $U_c = (U_1 + U_2) / 2$
- 雷诺数: $Re = \Delta U \cdot \delta_\omega / \nu$（$\delta_\omega$ 为涡量厚度）
涡量厚度的定义是，
这是由速度剖面的最大梯度定义的混合层的代表性厚度。
Kelvin-Helmholtz 不稳定性
请讲解一下 K-H 不稳定性的理论。
从速度不连续面（vortex sheet）的线性稳定性分析出发。在不连续面上施加微小波状扰动 $\eta \propto e^{i(kx - \omega t)}$，则色散关系为，
由于虚部为正，所有波数 $k$ 下的扰动都会增长。也就是说速度不连续面在所有波数下都是不稳定的。增长率为 $\sigma = k \Delta U / 2$，波长越短增长越快。
所有波数都不稳定，意思是任何微小的扰动都会增长吗？
理论上如此，但实际上对于有限厚度的剪切层，粘性和厚度效应会抑制短波长的不稳定性。最被放大的波长约为 $\lambda \approx 7 \delta_\omega$，比它短的波会被稳定化。
自相似解与扩展率
混合层会向下游发展对吧。
是的。在足够下游处，混合层会以自相似的方式扩展。涡量厚度随 $x$ 成比例增大。
扩展率 $d\delta_\omega / dx$ 依赖于速度比。$r = 0$（单侧静止）情况下的实验值约为 $d\delta_\omega / dx \approx 0.16\text{--}0.18$。Brown & Roshko (1974) 的可视化实验是该领域的里程碑式研究。
自相似区域的速度剖面是怎样的？
用 $\eta = y / \delta_\omega$ 无量纲化后，
这里 $F$ 具有误差函数型的形状。在 Goertler 的解析解中，表示为 $F(\eta) = \frac{1}{2}[1 + \text{erf}(\sigma \eta)]$。
混合层理论的先驱——Kelvin勋爵与Helmholtz的业绩（1868年）
具有速度差的两个流体相接的“混合层”的不稳定性，由 Hermann von Helmholtz 和 Lord Kelvin（William Thomson）于1868年独立进行了理论解析。冠以两人之名的 Kelvin-Helmholtz（KH）不稳定性，是线性稳定性理论的先驱成果，它表明在界面张力为零的情况下，任何速度差都会产生不稳定性。有趣的是，该理论最初是在“大气科学”的语境下（卷云的形成机制）被讨论的。其向工程应用（射流·混合层设计）的展开，不得不等待半个世纪后实验流体力学的发展。现代 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 使用什么方法？
| 方法 | 适用场景 | 备注 |
|---|---|---|
| DNS | 基础研究。$Re_{\delta_\omega} < 10^4$ | 完全解析 K-H 涡的卷起与配对 |
| LES | 中～高 Re 数的混合过程 | SGS 模型的影响相对较小（自由剪切流） |
| RANS | 时间平均的扩展预测 | 涡的详细结构会丢失，但扩展率可以预测 |
| 时间混合层 | 涡动力学研究 | 在对流坐标系中模拟空间混合层。使用周期边界计算 |
时间混合层 vs 空间混合层
时间混合层是什么？
空间混合层沿流动方向发展，而时间混合层以均匀剪切为初始条件进行时间演化。计算域在流向和展向采用周期边界，对应于对空间混合层进行伽利略变换后的结果。
这是 DNS 中常用的方法。计算成本低，且易于对统计量进行空间平均。但是，无法模拟空间混合层的“来自上游扰动的传播”。
初始条件与扰动的设定
初始条件怎么设定？
时间混合层的 DNS/LES 中，
1. 基本流: $\bar{u}(y) = \frac{\Delta U}{2} \tanh(2y / \delta_{\omega,0})$
2. 2D扰动: 叠加最被放大的波长的模态。基本模态与次谐波（诱发涡配对）
3. 3D扰动: 添加展向模态（oblique modes）。三维转换所必需
4. 随机噪声: 宽频带扰动。再现自然的湍流转换
根据 Michalke (1964) 的线性稳定性分析，求得最被放大的模态的波数为 $k_{max} \delta_{\omega,0} / 2 \approx 0.4457$。
网格设计
网格该怎么划分？
这是时间混合层网格设计的要点。
- 流向 ($x$): 周期边界。长度至少为基本模态波长的 4 倍以上（为了追踪两次配对）
- 法向 ($y$): 取得足够宽以追踪混合层的增长。中心细密，远处稀疏（拉伸）
- 展向 ($z$): 周期边界。至少为三维结构波长的 2 倍。$L_z \geq 2\lambda_z$
- 分辨率: DNS 的话 $\Delta x \approx \Delta z \approx \delta_{\omega,0} / 10$。$\Delta y_{min} \approx \delta_{\omega,0} / 20$
对流项格式
对流项的格式用什么好？
混合层的 LES/DNS 中，低耗散格式是必须的。
- DNS: 中心差分（2 阶或 4 阶）。能量守恒格式最理想
- LES: Bounded Central Differencing（Fluent）、LUST（OpenFOAM）、中心差分+少量迎风
- RANS: Second Order Upwind 就足够
用迎风差分的话涡会消失掉对吧。
なった
詳しく
報告