| name | jet-impingement |
| description | Jet impingement heat transfer — stagnation Nusselt number, Martin correlation, H/D ratio, array impingement, spent flow management, crossflow degradation, confined vs. unconfined, turbine blade leading edge impingement, ANSYS CFD validation, pressure drop, pumping power. |
| metadata | {"priority":7,"promptSignals":{"phrases":["jet impingement","impingement cooling","impingement heat transfer","jet array cooling","impingement Nusselt","stagnation point heat transfer"],"minScore":3}} |
Jet Impingement Heat Transfer — Complete Skill
Fundamentals
Impingement Regions
Three zones of jet impingement:
- Free jet region: jet exits nozzle; potential core (length ≈ 6–8 D) + developing region
- Stagnation zone: jet deflects parallel to surface; thin boundary layer; maximum heat transfer
- Wall jet region: flow spreads radially; boundary layer grows; decreasing heat transfer
Stagnation point Nusselt number:
Nu_stag = h_stag × D / k [dimensionless; h_stag = peak heat transfer coefficient; D = nozzle diameter; k = fluid conductivity]
Martin Correlation (Single Round Jet)
Martin (1977) — most widely used for engineering:
Nu(r, H, D, Re, Pr) = G × F × Re^0.667 × Pr^0.42 [valid: 2,000 ≤ Re ≤ 400,000; 2 ≤ H/D ≤ 12; 0.004 ≤ f ≤ 0.04]
G — radial distribution function:
G = (1 + (H/(0.6D))^6)^(-0.05) × (2√f × (1 - 2.2√f)) / (1 + 0.2(H/D - 6)√f)
F — radial function from stagnation:
At stagnation point (r = 0): F = F₀ (maximum)
F₀ = 1 - [1.1 × r/(D × f^(-1/2))]^n where n varies with H/D
Simplified stagnation Nusselt (Jambunathan correlation):
Nu_0 = C × Re^m × Pr^0.4 [C = 0.5–1.0; m = 0.5–0.7 depending on H/D]
For H/D = 6 (optimal): Nu_0 = 0.7 × Re^0.5 × Pr^0.4 [approximately]
Reynolds number for jet:
Re = ρ × V_j × D / μ = 4 × ṁ / (π × D × μ) [V_j = jet exit velocity; ṁ = mass flow rate per nozzle]
Single Round Jet — Stagnation and Radial Distribution
Average Nu over target area (circular zone radius r₀ = 2D):
Nu_avg ≈ 0.5 × Nu_0 [rough estimate; detailed integration needed from Martin]
H/D ratio effects:
H/D = 2–4: close impingement; short potential core; high Re sensitivity; risk of jet recirculation
H/D = 6–8: optimal for most applications; maximum Nu₀; potential core just reaches surface
H/D > 12: jet decelerates; Nu decreases; turbulence dissipates; velocity uniformity better (preferable for arrays)
Optimal H/D = 6: maximum stagnation heat transfer at given Re and pumping power
Nozzle area ratio (f):
f = A_nozzle / A_target [fraction of target covered by nozzle area; f = π/4 × D²/p² for square array pitch p]
Optimal f ≈ 0.01–0.04 for maximum heat transfer from Martin correlation
Array Impingement
Inline and Staggered Arrays
Array configuration parameters:
S_L = streamwise pitch between rows [mm]
S_T = transverse pitch between nozzles [mm]
Open area ratio: f = (π/4 × D²) / (S_L × S_T)
Crossflow effect (inline array with impingement on bottom):
Spent flow must exit — if channeled in one direction → crossflow velocity increases with streamwise position → degrades Nu for downstream jets
Crossflow momentum ratio:
G_c/G_j = (number of upstream jets × ṁ_jet) / (ṁ_jet × √(S_L/D))
Corrected Nu for crossflow (Florschuetz et al.):
Nu_i / Nu_0 = 1 - C_n × (G_c / G_j)^1.2 [C_n ≈ 0.6–0.8 depending on geometry; Nu_0 = stagnation Nu without crossflow]
Downstream jets in adverse crossflow: Nu degradation up to 30–50% for long arrays
Design strategies to mitigate crossflow:
- Staggered outflow (alternate spent flow direction between rows)
- Staggered nozzle layout (crossflow distributes more uniformly)
- Exit holes between jets (spent flow exits between impingement zones)
- Confined channel with exit ports every 4–6 jets
Slot Jets (2D Jets)
Slot jet Nu (Gardon-Akfirat):
Nu = C × Re^n [Re = ρ × V_j × W / μ; W = slot width]
n = 0.625 (for inline slot), 0.695 (turbulent)
Advantageous for elongated targets (turbine trailing edge cooling)
Confined vs. Unconfined Impingement
Unconfined (free jet impinging on surface):
Spent flow exits freely in all radial directions; crossflow not an issue for single jet
Nu lower than confined due to less turbulence augmentation from wall
Confined (nozzle plate close to target surface):
Spent flow must pass between nozzle plate and target → channel flow + impingement
H/D < 4 → significant confinement; recirculation between nozzle and target
Confinement increases turbulence → augments Nu near stagnation (up to 20–30%)
Can cause pressure oscillations and noise at small H/D (< 2)
Turbine Blade Leading Edge Impingement
Application: internal jet array cools turbine leading edge (highest heat load zone)
Jet exits from insert holes → impinges on leading edge inner surface → exits via showerhead film holes
Leading edge thermal load:
q'' = h_external × (T_gas - T_blade) + h_internal × (T_blade - T_coolant) [net heat balance]
Thermal effectiveness: φ = (T_gas - T_blade) / (T_gas - T_coolant_in) → want φ near 1
Insert impingement holes:
Typically D = 0.5–2 mm; H/D = 3–6; 6–12 holes per leading edge row
Material: nickel superalloy insert (same as blade or cast separately)
Crossflow in leading edge cavity: spent coolant exits via downstream holes → limited crossflow for first few jets
SENB geometry for impingement in curved surface: use effective D and radial distance r/D corrected for curvature
Pressure Drop and Pumping Power
Nozzle pressure drop:
ΔP_nozzle = (1/2) × ρ × V_j² × (1/C_d² - 1) [C_d ≈ 0.61 for sharp-edged orifice; lower for rounded entry]
For C_d = 0.61: ΔP_nozzle ≈ 1.69 × (1/2) × ρ × V_j²
Total pressure drop:
ΔP_total = ΔP_supply + ΔP_nozzle + ΔP_exit [supply plenum + nozzle + discharge]
Pumping power:
P_pump = Q_flow × ΔP_total / η_pump [W; Q_flow = total volumetric flow; η_pump = 0.5–0.8]
Figure of merit (FOM):
FOM = Q_heat / P_pump = h_avg × A_target × ΔT / (Q_flow × ΔP_total) [dimensionless]
Higher FOM → better thermal-hydraulic performance; use to compare configurations
Performance factor:
J = Nu_avg / Re^0.5 [Colburn-like j-factor; higher J = more efficient heat transfer per momentum]
CFD Validation Approach
Turbulence models for impingement:
Standard k-ε: overestimates Nu at stagnation (known issue)
Realizable k-ε / SST k-ω: improved stagnation prediction; ±10–20% accuracy
RSM (Reynolds Stress Model): best accuracy for complex impingement; high cost
SST k-ω with γ-Reθ transition model: for low Re or transition effects
Mesh requirements:
Y+ = 1–5 for wall-resolved; 5 integration points in viscous sublayer
Nozzle region: structured mesh; 15–20 cells across D
Impingement zone: refined mesh; Δr < D/20 in radial direction
Validation procedure:
Compare Nu distribution with Martin correlation or experimental data; target ±15% at stagnation; ±25% in wall jet region
Design Example
Problem: cool a 50 × 50 mm electronics chip with a water jet array; q'' = 100 W/cm² = 10⁶ W/m²; T_max = 85°C; T_water = 25°C
Required h: h = q'' / (T_max - T_water) = 10⁶ / 60 = 16,667 W/(m²·K)
Required Nu for D = 2 mm nozzle, water k = 0.6 W/(m·K):
Nu = h × D / k = 16,667 × 0.002 / 0.6 = 55.6
From Martin correlation, H/D = 6:
55.6 = 0.7 × Re^0.5 × Pr^0.4 [Pr_water ≈ 6.1 at 25°C]
Re^0.5 = 55.6 / (0.7 × 6.1^0.4) = 55.6 / (0.7 × 2.10) = 37.8 → Re = 1,429 (laminar edge)
Need higher Re → increase D to 3 mm or use higher velocity
Standards and References
| Source | Scope |
|---|
| Martin (1977) ASME J. Heat Transfer | Fundamental correlation for round and slot jets |
| Florschuetz et al. (1981) NASA TM | Array impingement with crossflow |
| Jambunathan et al. (1992) Int. J. Heat Fluid Flow | Review of round jet impingement |
| ASME PTC 30.1 | Air-cooled heat exchangers (includes jet principles) |
| NASA TM-2003-212601 | Turbine blade impingement cooling design |
| ASTM E1969 | Standard practice for impingement testing (erosion context) |
Output
Provide: jet geometry (D [mm], H/D, array pitch S_L × S_T [mm]), fluid properties (k [W/m·K], μ [Pa·s], Pr), jet velocity V_j [m/s] and Reynolds number Re, Martin correlation (stagnation Nu_0, average Nu_avg), heat transfer coefficient h [W/(m²·K)], heat flux capacity at given ΔT [W/cm²], crossflow parameter G_c/G_j and Nu degradation [%] for array configuration, nozzle pressure drop ΔP [kPa], total pumping power [W], thermal-hydraulic figure of merit FOM, comparison with CFD validation (turbulence model, mesh Y+, ±% error vs. correlation), and applicable reference (Martin 1977, Florschuetz 1981, ASME PTC 30.1).