碰撞射流热传递

碰撞射流（jet impingement）是指从喷嘴喷出的射流碰撞到壁面的流动，在碰撞点（stagnation point）附近可获得非常高的热传递系数。可实现通常管内强制对流的2～3倍热传递率，广泛应用于燃气轮机叶片内部冷却、钢板淬火、电子部件冷却等领域。

流动分为三个区域。(1) 自由射流区域（free jet region）：从喷嘴到壁面间存在势流核心的区域。(2) 碰撞域（impingement zone）：壁面附近流动急剧改变方向的区域。(3) 壁面射流区域（wall jet region）：沿壁面呈放射状扩展的区域。

有的。Martin（1977）提出的相关式被广泛应用。单个圆形喷嘴的停滞点努塞尔数为

超过势流核心长度（通常4～6D）后，射流会扩散导致碰撞时速度下降，努塞尔数随之减小。反之当 $H/D < 4$ 时，在封闭几何体中会出现cross-flow效应（使用过的射流干扰新射流），性能会下降。

正确。配列射流（array impingement）中射流间干涉和cross-flow影响变得至关重要。孔间距 $S/D$ 越小，面积平均努塞尔数越高，但cross-flow对下游喷孔冷却性能的削弱也越明显。通常在 $S/D = 4$～$8$ 范围内设计。Florschuetz等（1981）的相关式是配列射流的标准参考数据。

完全同意。碰撞射流是乱流模型的标准基准问题，许多模型都会过度预测碰撞域的努塞尔数。主要原因是对碰撞域乱流能量生成项的过度评估。

具体而言，标准k-ε模型会将停滞点努塞尔数过度预测为实验值的2倍以上。v2f模型和SST k-ω有明显改善，但通常仍有20～30%的过度预测。基于ω方程的模型总体上优于k-ε系列。

经文献调查，v2f模型（Fluent中称为 $\overline{v^2}$-$f$ 模型）在碰撞射流中预测最接近实验值。其次是SST k-ω。但v2f仅在Fluent和OpenFOAM（v2f湍流模型）中可用，STAR-CCM+没有标准实装。

碰撞域需要 $y^+ < 1$ 来确保壁面法向精度，壁面平行方向相对喷嘴直径 $D$ 需 $\Delta r / D \approx 0.02$～$0.05$ 的分辨率。从喷嘴出口到壁面间也要布置足量的单元来解析射流发展。典型的二维轴对称计算需5万～20万单元，三维配列射流则需数百万单元。

单个圆形喷嘴可用轴对称（轴对称）模型，计算效率很高。Fluent和STAR-CCM+都提供轴对称求解器。但喷嘴为矩形或存在cross-flow时必须用三维。

研究阶段广泛应用。射流剪切层会产生Kelvin-Helmholtz不稳定性，碰撞域的大尺度涡结构导致努塞尔数时间变化。LES可直接解析涡结构，使时间平均努塞尔数分布比RANS更接近实验。DES（分离涡模拟）和SBES也是有用的中间选择。

高密度数据中心和新一代功率半导体中，微喷嘴液体碰撞冷却在研究中。直径0.5mm喷嘴配列可在芯片表面直接喷射，达到 $Nu \sim 100$～$500$ 的高性能。

首先用单喷嘴轴对称问题验证乱流模型。Baughn & Shimizu（1989）和Cooper等（1993）的实验数据是标准验证数据集。对停滞点努塞尔数和径向努塞尔数分布都要比对。

配列射流用Florschuetz等的实验数据比对。CFD结果需评估面积平均努塞尔数、喷嘴中心线上的努塞尔数分布以及喷嘴间最小努塞尔数。

(1) 首先做网格收敛性分析排除数值误差。(2) 尝试不同乱流模型（v2f、SST k-ω、RSM）。(3) 还是不行则改变壁面处理方式（壁面函数 vs 低雷数）。(4) 最后手段是做LES或DES的非定常计算。实务中通常以"全工况20%误差以内"为目标。

以定常轴对称单喷嘴为例。(1) 在Mesh Setup中选择Axisymmetric模型。(2) 乱流模型选Viscous Models > k-omega > SST。如可能则试v2f模型（Viscous Models > Transition > v2f）。(3) 启用Enhanced Wall Treatment。(4) 打开Energy方程。

边界条件：喷嘴入口设velocity inlet（指定流速和温度），壁面设constant heat flux或constant temperature。出口设pressure outlet。若无约束壁，则在上方和侧面设pressure outlet。

从零速度场开始会收敛缓慢。沿射流方向给定均匀速度作初值会加快收敛。也可用Hybrid Initialization自动设置。

Physics Models中选Steady + Segregated Flow + Segregated Fluid Temperature + Turbulent + SST k-ω。确认All y+ Wall Treatment已启用。碰撞射流对壁面附近网格要求高，用Prism Layer Mesher确保总厚度和层数充分。

STAR-CCM+标准配置中无v2f模型，碰撞域过度预测的疑虑时可试Elliptic Blending k-epsilon（EB k-ε）模型。其壁面近处处理与v2f相当。

求解器用simpleFoam或buoyantSimpleFoam。在constant/turbulenceProperties中指定kOmegaSST。v2f模型在OpenFOAM v2212以后的v2f class中可用。壁面边界条件用nutUSpaldingWallFunction、kqRWallFunction、omegaWallFunction。

用blockMeshDict创建楔形网格（角度约5度），前后面设wedge边界条件。或用extrudeMeshDict的axisymmetric选项。中心轴不用empty（二维）而用axis边界条件。

旋转通过科里奥利力和离心浮力两种方式影响碰撞射流的热传递。旋转数 $Ro = \Omega D / u_j$ 达0.1以上时，努塞尔数分布相对非旋转状态明显不同。压力面侧科里奥利力将喷流压向壁面，努塞尔数增加；吸力面侧则相反，努塞尔数下降。

用旋转坐标系（Rotating Reference Frame / Moving Reference Frame）。Fluent的Frame Motion或sliding mesh、STAR-CCM+的Rotating Physics、OpenFOAM的SRFSimpleFoam都可设置。旋转坐标系中科里奥利力和离心力以源项形式加入动量方程。

伴沸碰撞射流更复杂。核沸腾区气泡生成与离脱大幅增强热传递，$h$ 可达单相流的10倍以上。临界热流密度（CHF）超限后转为膜沸腾，传热急剧恶化。

用VOF法或欧拉多相流模型追踪气液界面，结合蒸发-冷凝模型可近似求解。Fluent的Evaporation-Condensation Model和STAR-CCM+的Boiling Model都基于Lee模型处理蒸发/冷凝。但气泡核生成的直接求解计算成本太大，实务中多用Wall Boiling Model（RPI模型）这类壁面沸腾模型。

IBM和Intel的研究团队发表的微喷嘴冷却使用直径 $D = 50$～$500$μm、$H/D = 2$～$5$ 的配列，目标除热>1000 W/cm²。CFD中层流-乱流转捩很重要，$Re_D = 500$～$5000$ 范围须覆盖。该区域Transition SST模型有用。

硅基的通孔（TSV）加工成喷嘴的"集成微喷嘴冷却"也在研究，作为三维集成电路（3D-IC）冷却方案备受关注。这类微细结构的CFD还需考虑壁面粗糙度和滑移速度（slip velocity）影响。

这是碰撞射流CFD最常见的问题。原因是乱流模型对碰撞域乱流能量生成的过度评估。对策按优先顺序排列如下。

(1) 改变乱流模型：标准k-ε → SST k-ω → v2f（Fluent）/ EB k-ε（STAR-CCM+）。(2) 使用Kato-Launder修正：Fluent中乱流模型Options启用Production Limiter。STAR-CCM+中打开Realizability修正。(3) 试用RSM（雷诺应力模型）：计算成本增加但直接处理碰撞域各向异性。

标准乱流能量生成项 $P_k = \mu_t S^2$ 改为 $P_k = \mu_t S \Omega$（$S$ 为应变速率，$\Omega$ 为涡度）。碰撞域 $S$ 很大但 $\Omega$ 很小，故生成项被抑制。简单但效果显著。

cross-flow（使用过的射流积存干扰新射流）未被正确建模的可能。检查计算域出口边界条件。出口太近会导致cross-flow积存不足。RANS也易低估cross-flow对喷流的偏转。配列射流的精确预测需用LES。

碰撞射流壁面近处速度、温度梯度陡峭，网格依赖性强。应用Richardson外推法做三水准网格收敛性分析，直到GCI < 5% 为止。特别是碰撞点附近壁面平行方向网格常不足。碰撞域应确保 $\Delta r / D < 0.03$。

$y^+$ 远小于1（$y^+ < 0.1$）时，ω方程的壁面边界值会很大，数值不稳定。Fluent中将ω的欠松弛因子降至0.6。STAR-CCM+的All y+ Wall Treatment会根据 $y^+$ 值自动调整，相对稳定。