喷流（射流）
喷流（射流）的理论基础
概述
先生，喷流说白了就是从喷嘴喷出来的流动，对吧？
没错。喷流（jet）是从喷嘴或孔口射入周围流体中的流动。产业应用范围广泛。喷气发动机排气、焊枪、空调出风口、化工厂混合器、喷墨打印机等。
流体力学上，喷流是自由剪切流（free shear flow）的典型例子，与混合层和尾流并列为乱流的基本研究对象。
喷流的分类
喷流也有分类吗？
按几何形状分类是这样的。
| 种类 | 形状 | 自相似区速度衰减 | 扩展率 |
|---|---|---|---|
| 轴对称圆形喷流 | 圆形喷嘴 | $u_c / U_0 \propto (x/D)^{-1}$ | $\delta / x \approx 0.10$ |
| 平面喷流 | 狭缝喷嘴 | $u_c / U_0 \propto (x/h)^{-1/2}$ | $\delta / x \approx 0.11$ |
| 矩形喷流 | 矩形喷嘴 | 近处为平面喷流状，远处为轴对称性 | 与长宽比相关 |
轴对称的衰减速度更快呢。
对的。轴对称喷流由于周向各方向的夹带（entrain）作用，使动量扩散更快。
喷流的区域结构
让我们从上游到下游整理圆形喷流的结构。
1. 势核区（Potential Core） ($0 < x < x_c$)：喷嘴出口速度 $U_0$ 在中心轴保持不变。$x_c \approx 4\text{--}6D$
2. 过渡区 ($x_c < x < 20D$ 左右)：中心速度开始衰减
3. 自相似区 ($x > 20\text{--}30D$)：速度轮廓形成自相似形式
势核长度取决于入口乱流强度。乱流强度越高，势核越短。
自相似解
请告诉我自相似解的具体形式。
在轴对称喷流的自相似区中，时间平均速度轮廓具有以下形式。
中心速度的衰减由动量守恒推导。
其中 $B_u \approx 5.8\text{--}6.2$ 是实验常数，$x_0$ 是虚拟原点。在高斯轮廓假设下，
$B_u$ 的值因为研究者不同而略有差异，对吧。
这是因为初始条件（喷嘴出口边界层厚度、乱流强度、速度轮廓形状）的差异。Hussein et al. (1994) 的精密测量报告 $B_u = 5.8$，Panchapakesan & Lumley (1993) 报告 $B_u = 6.06$。
动量守恒
喷流中动量是守恒的，对吗？
周围为静止流体时，轴向动量通量保持不变。
将此关系与自相似轮廓假设相结合，得到 $u_c \propto x^{-1}$ 和 $\delta \propto x$。
喷流理论的确立——从Prandtl混合长理论到乱流喷流
自由喷流（Free Jet）的理论分析基于Prandtl(1925)的混合长理论发展而来。圆形自由喷流中心轴速度Uc随距喷嘴出口距离x的衰减规律为Uc ∝ x⁻¹，半值半径约每个出口径增加0.1倍，满足相似律。20世纪50-60年代，Tolmien、Görtler等得出了严格的解析解，之后Wygnanski & Fiedler(1969)的精密实验证实了乱流喷流的自相似性。这一自相似性（Self-Similarity）的发现成为现代RANS模型参数调校的基准，k-ε模型中Cμ=0.09的取值就是基于这一实验数据的历史。
喷流（射流）的数值计算方法
数值方法的选择
喷流的CFD采用什么样的手法？
喷流是自由剪切流，不需要壁面分辨，所以与LES的配合很好。
| 方法 | 应用场景 | 备注 |
|---|---|---|
| RANS ($k$-$\varepsilon$) | 时间平均扩展率预测 | 注意圆形喷流异常（round jet anomaly） |
| RANS (SST $k$-$\omega$) | 一般工程计算 | 比 $k$-$\varepsilon$ 更准确地预测喷流扩展 |
| LES | 喷流噪声、混合过程详解 | 入口条件的设置很关键 |
| DNS | 低Re喷流基础研究 | Re < $10^4$ 左右才可行 |
圆形喷流异常
什么是圆形喷流异常？
标准 $k$-$\varepsilon$ 模型对平面喷流的扩展率预测良好，但对轴对称喷流的扩展率预测约过大 $40\%$。这源于 $C_{\varepsilon 1}$ 常数的问题，平面喷流和轴对称喷流不能用同一常数。
对策如下：
- 将 $C_{\varepsilon 1}$ 从 $1.44$ 改为 $1.60$（Pope修正）
- 采用 SST $k$-$\omega$ 模型（改善了喷流扩展预测）
- 使用 Realizable $k$-$\varepsilon$（$C_{\mu}$ 变为可变，改善了喷流行为）
入口条件的设置
喷嘴出口的速度分布怎样设置？
用LES解喷流时，入口条件对结果影响很大。
- 均匀速度轮廓（top-hat）：最简单但不现实。喷嘴出口没有边界层，初始剪切层发展会改变
- 管流轮廓：$u(r) = U_c (1 - (r/R)^n)$。常用 $n=7$（乱流1/7幂律）
- 包含喷嘴内部的计算：最准确。直接求解喷嘴内边界层发展
乱流变动的注入也很重要。方法包括：
- 合成乱流生成法（SEM: Synthetic Eddy Method）：Jarrin et al. (2006)
- 循环法：从喷嘴内断面循环数据
- 数字滤波法：Klein et al. (2003)
仅指定乱流强度还不够啊。
对RANS来说，在入口指定 $k$ 和 $\varepsilon$（或 $\omega$）就够了。但LES需要在入口给定有空间和时间相关性的变动速度场，否则，非物理的适应区会变长，势核长度会偏离。
网格设计
喷流的网格有什么要注意的？
以下要点很重要。
- 喷嘴出口附近的剪切层：需要喷嘴唇厚度的1/10以下的网格。为了分辨剪切层初期不稳定性
- 轴向区域长度：如果想看到自相似区，需要 $30D$ 以上。噪声分析需要 $50D$ 以上
- 径向：喷流边界外侧也要保留充分的区域（$10D$ 以上）
- 夹带边界：侧面边界设为压力条件（允许夹带）。流速固定不行
如果侧面是墙的话，流入被阻止所以夹带会受阻，对吧。
正是。侧面的压力条件设置不当会在喷嘴近处产生非物理的低压区，影响喷流扩展。
喷流CFD的边界条件设置——入口乱流强度设置错误导致的误差
喷流CFD分析中容易被忽视的是入口边界条件，特别是乱流强度（TI）和乱流尺度（L）的设置。如果乱流强度设置过大，喷流核的衰减比实验快，过小则势核区会过度延伸。在没有实验数据参考的情况下，喷流直径D的5%左右的乱流强度和D的10%左右的尺度是工程初始估计值。喷嘴内流速轮廓是完全展开流还是不是，也会很大影响扩散特性——"均匀轮廓"假设要与实验条件（是否有收缩喷嘴）配合选择。没有事前实验数据的情况下，细化入口条件是喷流CFD精度改进的第一步。
喷流（射流）的工程应用
分析流程
请教我喷流CFD分析的步骤。
典型流程如下。
1. 目的明确化：时间平均扩展预测（RANS）还是混合过程、噪声的详细分析（LES）
2. 喷嘴形状定义：圆形、矩形、同轴二重（coaxial）、收缩喷嘴等
3. 计算区域设计：轴向 $30\text{--}50D$，径向 $10\text{--}15D$
4. 网格生成：剪切层部分细密处理。结构化网格（O-H型）在精度上有优势
5. 边界条件：入口轮廓、侧面出口压力条件、有coflow时设定外部流
6. 计算实行与统计：LES时，扣除初期瞬态，统计 $50\text{--}100$ 流动通过次数以上
7. 验证：$u_c(x)$、半值宽 $r_{1/2}(x)$、Reynolds应力轮廓与文献对比
验证用基准
应与什么样的实验数据比较？
轴对称圆形喷流的代表基准数据是这些。
| 研究者 | Re | $B_u$ | $r_{1/2}/x$ | 备注 |
|---|---|---|---|---|
| Panchapakesan & Lumley (1993) | 11,000 | 6.06 | 0.096 | HWA计测，参考数据的标准 |
| Hussein, Capp & George (1994) | 95,500 | 5.8 | 0.094 | LDA+HWA，高Re |
| Burattini et al. (2005) | 50,000 | 5.9 | 0.095 | PIV计测 |
$r_{1/2}$ 是半值宽，也就是速度为中心速度一半处的半径，对吧。
对。$r_{1/2} / x \approx 0.094\text{--}0.096$ 是一个非常稳健的值，不太依赖Re或喷嘴条件。如果合不上，说明网格或乱流模型有问题。
常见设置错误
喷流计算容易陷入什么误区？
列举几个代表性的吧。
| 问题 | 原因 | 对策 |
|---|---|---|
| 扩展率过大 | 标准 $k$-$\varepsilon$ 的圆形喷流异常 | 改为 Realizable $k$-$\varepsilon$ 或 SST |
| 势核过长 | LES入口没有乱流变动 | 用SEM或数字滤波法注入变动 |
| 喷流偏向 | 网格非对称性 | 采用圆柱坐标网格或足够细的非结构网格 |
| 夹带不足 | 侧面边界为墙面条件 | 侧面改为压力inlet/outlet |
| 统计不收敛 | 抽样时间不足 | 确保 flow-through time 100次以上 |
flow-through time 是什么？
$T_{ft} = L / U_0$。$L$ 是计算区域的轴向长度。相当于喷流通过一次计算区域的时间。统计的独立样本数要基于这个倍数来估计计算时间。
工厂排气口的喷流扩散——CFD揭示了烟囱高度设计的秘诀
工厂烟囱或排气口喷出射流的扩散分析是环保评估的必要项目。自由喷流的中心轴速度随距离x的衰减遵循1/x律，半值宽约每直径D增加0.1倍，符合高斯分布理论值。实际工程CFD案例中，相邻建筑的尾流导致的"下洗"引发了意想不到的浓度上升，超过了地方政府的规制值。CFD分析表明，增加排气速度比提高烟囱高度更有效率地改善扩散，这一发现被用于优化设备改造的成本。这类案例已被多家企业采纳。
喷流（射流）的软件对比
工具别喷流分析功能
适合喷流分析的CFD工具是哪些？
按喷流分析相关功能比较。
| 工具 | 乱流模型 | 噪声分析 | 合成乱流入口 |
|---|---|---|---|
| Ansys Fluent | 所有RANS、LES (WALE, Dynamic), DES | FW-H声学类比 | 涡旋法、SEM (UDF) |
| STAR-CCM+ | 同上 + IDDES | FW-H、宽频噪声 | 合成乱流生成器 |
| OpenFOAM | 全模型自由实现 | FW-H (libAcoustics) | turbulentDFSEMInlet |
| Nek5000 | DNS/LES（谱元法） | 直接计算 | 循环法 |
喷流噪声分析工作流
喷流噪声的计算怎样进行？
有两个方法。
直接法（CAA: 计算气动声学）：直接求解声波传播。需要以声波波长对网格进行离散，计算成本极高
混合法：近场用LES求解，远场声传播用Ffowcs Williams-Hawkings (FW-H) 方程计算
Fluent中启用声学模型 → FW-H，在围住喷流的封闭面（可透过）收集数据，就可计算远场声压。
OpenFOAM 中的喷流LES
用OpenFOAM做喷流LES的设置请教一下。
pimpleFoam + WALE SGS模型是标准配置。
```
constant/turbulenceProperties:
simulationType LES;
LES { LESModel WALE; delta cubeRootVol; }
system/fvSchemes:
ddt: backward;
div(phi,U): Gauss LUST grad(U); // 75%中心差分 + 25%风上
0/U (入口):
type turbulentDFSEMInlet;
delta 2.0;
nCellsEddy 5;
mapMethod nearestCell;
```
LUST格式是第一次听说。
Linear Upwind Stabilized Transport的缩写，是中心差分和线性风上的混合格式。LES中纯中心差分容易不稳定，所以常用这种混合格式。
选型指南
到底应该选哪个工具？
看目的而定。
- RANS时间平均混合预测：Fluent或STAR-CCM+。GUI设置方便
- LES详细涡结构分析：OpenFOAM（自由度高）。Nek5000（高精度DNS/LES）
- 喷流噪声：Fluent的FW-H最成熟。STAR-CCM+也支持。OpenFOAM有libAcoustics库
- 燃烧喷流：Fluent的flamelet/PDF模型、OpenFOAM的reactingFoam、或Cantera耦合
喷流CFD商用工具选择——OpenFOAM和Fluent的精度与手工成本权衡
对于自由喷流、冲击喷流的CFD分析，乱流模型的选择和工具的实现质量影响结果很大。ANSYS公司的Fluent的RSM（Reynolds应力模型）实现成熟，能相对好地再现喷流的各向异性乱流。OpenFOAM的LES求解器丰富，适合喷流混合的瞬间结构分析——但前提是需要超算环保和HPC运维技能。StarCCM+的DES（Detached Eddy Simulation）的RANS-LES切换边界可通过GUI直观确认。小规模设计检查通常用RANS/Fluent，研究用途的高精度分析用LES/OpenFOAM，已成为多数实务工程师的事实标准。
喷流（射流）的前沿研究
喷流不稳定性理论
喷流的涡结构是由什么不稳定性决定的？
喷流剪切层有多种不稳定性模式。
- Kelvin-Helmholtz不稳定性：剪切层的卷起。最初出现的不稳定性。波长与剪切层的动量厚度 $\theta$ 成正比
- 首选模式：整个喷流的柱状不稳定性。$\mathrm{St}_D = fD/U_0 \approx 0.3\text{--}0.5$
- 方位角模式：$m = 0$（轴对称）、$m = \pm 1$（螺旋形）等。$m = 0$ 和 $m = 1$ 最易放大
Crow & Champagne (1971) 报告喷流的首选模式为 $\mathrm{St}_D \approx 0.30$。这对应于喷流柱的固有不稳定性。
大尺度涡结构与噪声
喷流噪声与涡结构的关系是怎样的？
亚音速喷流的噪声有两个成分。
1. 细尺度乱流噪声：小尺度乱流的声音。向全方向辐射。高频成分
2. 大尺度结构噪声：大尺度波动结构（不稳定波）的声音。向下游方向（浅角度）辐射。峰值周频率 $\mathrm{St}_D \approx 0.2$
后者也叫Mach波辐射。当对流速度超过周围音速时会形成Mach锥。在超音速喷流中起主导作用。
喷流控制
有什么方法可以促进喷流混合或降低噪声？
介绍几个已经实用化的方法。
- 锯齿形喷嘴：喷嘴后端有锯齿突起。破坏大尺度结构来降低噪声约 $2\text{--}3$ dB。实际航空发动机中使用
- 微喷流注入：从喷嘴唇周围喷射小喷流。控制剪切层发展
- 等离子体驱动：用DBD等离子体向流动注入动量
- 突起物（tabs）：在喷嘴内壁设置突起。产生轴向涡增加混合
DNS/LES 的最前沿
喷流大规模计算有什么成果？
此外，利用物理信息神经网络（PINN）从有限的测量数据重建整个喷流流场的尝试也在进行。Raissi等人的研究是先驱。
10亿网格真是太多了。在什么样的计算机上运行？
数千到数万块GPU的大规模并行计算。Charles（Cascade Technologies的非结构网格求解器）针对GPU优化。学术代码中nekRS（Nek5000的GPU版）也备受关注。
超音速喷流的冲击波晶胞结构——尖啸音与CFD分析的局限
超音速喷流（不完全膨胀喷流）中形成交替出现压缩波和膨胀波的"冲击波晶胞（Shock Cell）"结构。这种结构与剪切层的相互作用产生离散的"尖啸音（Screech Tone）"，是航空发动机和风洞试验的噪声问题。频率可用Powell公式按喷嘴径和出口马赫数预测，但CFD（RANS）很难精密再现尖啸音的准确频率和声压。目前最好的方法是LES+压缩可压声学（CAA）的组合，但计算成本是普通分析的100倍以上。
喷流（射流）的故障排除
常见故障
喷流计算有哪些常见问题？
逐个介绍。
1. 扩展率与实验不符
症状：喷流半值宽 $r_{1/2}$ 的梯度与文献值（$\approx 0.094$）偏差 $20\%$ 以上。
原因和对策：
- 使用了标准 $k$-$\varepsilon$ → 改用Realizable $k$-$\varepsilon$ 或SST $k$-$\omega$（圆形喷流异常）
- 计算区域太短 → 自相似区（$x > 20D$）要包含在内
- 喷嘴出口条件不合适 → 用再现实验出口边界层厚度的入口轮廓
2. 势核过长（LES）
实验中 $x_c \approx 5D$，但LES得到 $8D$ 左右。
原因：入口没有给定现实的乱流变动。从层流剪切层的Kelvin-Helmholtz涡卷起延迟。
对策：
- 使用合成乱流入口条件（SEM、DFSEM）
- 计算域包含喷嘴内部（从上游 $3\text{--}5D$ 处从喷嘴内部开始求解）
- 调整入口乱流强度到实验值（通常 $1\text{--}5\%$）
3. 喷流非物理性偏向
原因和对策：
| 原因 | 对策 |
|---|---|
| 网格的非对称性 | 采用绕喷流轴对称的网格。O-H型结构网格最好 |
| 出口边界条件的反射 | 出口远离处理（$> 30D$）。采用对流流出条件 |
| 侧面边界的影响 | 径向区域扩大到 $> 10D$ |
| coflow的不均一性 | 均一化外部流入口条件 |
4. 夹带边界的问题
从侧面流入的流体处理困难。
喷流夹带周围流体，所以侧面边界必须允许自然流入。
- 压力inlet/outlet：指定全压和静压。允许反流
- Fluent：侧面设为Pressure Outlet，"反流方向"改为"垂直于边界"
- OpenFOAM：
pressureInletOutletVelocity+totalPressure组合 - NG：速度inlet（均匀流 $U = 0$）会阻碍流入
5. 统计量收敛困难
喷流的统计量（特别是Reynolds应力）收敛很慢。
目安：
- 一阶统计（平均速度）：约 $50 T_{ft}$ 收敛
- 二阶统计（Reynolds应力）：需 $100\text{--}200 T_{ft}$
- 三阶及以上（斜度等）：需 $500 T_{ft}$ 以上
需要很长的计算时间呢。
对的。加快统计收敛的方法是结合方位角方向的平均（轴对称喷流），或多个快照的集合平均。
喷流CFD改变乱流模型时结果大幅变化——标准k-ε的"怪癖"
自由喷流的CFD分析中，"改变乱流模型后中心速度衰减率变化很大"的经历不少工程师都有过。标准k-ε模型对轴对称喷流的扩散有20～30%的过高评估，这从很久前就为人所知（Pope,1978）。这称为圆形喷流异常（Round Jet Anomaly），源于标准k-ε的模型常数Cε1=1.44是针对平板喷流优化造成的。为了解决这个问题，开发了Pope修正（涡伸长项追加）、RNGk-ε、Realizable k-ε等，但根本解决还没到位。喷流CFD务必进行多模型结果比较，并与实验值对照，这是必须的。
相关主题
价值
更详细
报告