Correlazioni di altezza di fiamma nella Fire Safety Engineering
Raccolta strutturata delle principali correlazioni di altezza di fiamma usate nella Fire Safety Engineering, con formule, variabili, limiti di validità, note implementative e riferimenti bibliografici canonici.
Flame height correlations in Fire Safety Engineering
Context
This note organizes the main flame height correlations commonly referenced in Fire Safety Engineering for preliminary calculations, comparative assessments, and tool implementation.
The source material has been structured from the attached technical reference document, preserving formulas, variable definitions, typical outputs, validity ranges, and canonical bibliography.
General conventions
Q= heat release rate [kW]Q*= dimensionless fire Froude numberD= effective fire diameter [m]L= flame height [m]z= vertical coordinate from fire base [m]
Default ambient constants:
rho_inf = 1.2 kg/m³c_p = 1.005 kJ/(kg·K)T_inf = 293 Kg = 9.81 m/s²
Heskestad (1984, 1995)
Formula
L = -1.02 * D + 0.235 * Q^(2/5)
dove:
L = altezza media di fiamma (mean flame height) [m]
D = diametro del focolare [m]
Q = heat release rate [kW]
Forma adimensionale equivalente (utile per verifica):
(L + 1.02*D) / D = 0.235 * Q^(2/5) / D
ovvero: L/D + 1.02 = 0.235 * (Q / D^(5/2))^(2/5) * D^0 ...
In termini di Q*:
L/D = -1.02 + 3.7 * Q*^(2/5)
Variables
L : mean flame height [m]
Definita come la quota media del tip della fiamma visibile,
pari alla quota in cui la fiamma è presente il 50% del tempo.
D : diametro del focolare [m]
Per focolari non circolari usare il diametro equivalente: D = sqrt(4*A/pi)
Q : heat release rate convettivo + radiativo (totale) [kW]
Inputs and outputs
Input: D [m], Q [kW]
Output: L [m]
Output type: Mean flame height (altezza media): la fiamma è presente al 50% del tempo
a quota L. È la definizione più comune in fire engineering.
Validity and use
- Focolari circolari o con geometria per cui il diametro equivalente è
una buona approssimazione
- Q* da circa 0.001 a oltre 1000 (range sperimentale molto ampio)
- Non adatta per fiamme in crossflow (vento) o geometrie molto elongate
- Per Q* < ~1 (fiamme "tall and thin") la correlazione rimane valida
ma la precisione decresce leggermente
- Combustibili gassosi e liquidi (non validata specificamente per solidi
con strong pyrolysis front)
Notes
La correlazione di Heskestad è la più citata e usata in fire engineering.
È la base per la modellazione dei plume in NFPA 92, CIBSE Guide E, ecc.
Reference
Heskestad, G. (1984). Engineering Relations for Fire Plumes. Fire Safety Journal, 7(1), 25–32. [Versione aggiornata:] Heskestad, G. (1995). Fire Plumes. In: SFPE Handbook of Fire Protection Engineering, 2nd ed., Society of Fire Protection Engineers, Bethesda, MD. Section 2, Chapter 2.
Thomas (1963)
Formula
L / D = 0.0037 * (Q / (D^(5/2)))^0.61
oppure nella forma originale con portata massica:
L / D = C * (m_dot / (rho_inf * g^0.5 * D^(5/2)))^n
Forma più diffusa in letteratura (con Q):
L = 0.0037 * Q^0.61 * D^(-0.028)
(approssimazione; la forma esatta va ricavata dalla formulazione
adimensionale originale — vedere note)
Forma adimensionale originale di Thomas:
L / D = 55 * (m_dot / (rho_inf * sqrt(g*D) * D^2))^0.67
dove m_dot = portata massica di combustibile [kg/s]
Variables
L : mean flame height [m]
D : diametro del focolare [m]
m_dot : mass burning rate (portata massica di combustibile) [kg/s]
rho_inf : densità dell'aria ambiente [kg/m³] (default: 1.2 kg/m³)
g : accelerazione di gravità [m/s²] (= 9.81 m/s²)
Conversione Q → m_dot: m_dot = Q / (chi * deltaHc)
con chi = combustion efficiency (tipicamente 0.7–1.0)
e deltaHc = calore di combustione del combustibile [kJ/kg]
Inputs and outputs
Input: D [m], m_dot [kg/s], rho_inf [kg/m³], g [m/s²]
oppure D [m], Q [kW], deltaHc [kJ/kg], chi [-]
Output: L [m]
Output type: Mean flame height
Validity and use
- Correlazione sviluppata principalmente per legno e combustibili solidi
- Valida per focolari di grandi dimensioni (>1 m)
- Meno precisa per combustibili gassosi e focolari piccoli
- Non applicabile per fiamme in crossflow
Notes
La correlazione di Thomas è storicamente importante ma è stata in gran parte
sostituita dalla correlazione di Heskestad per applicazioni generali.
Rimane di riferimento per incendi di grandi dimensioni e per il confronto
con Heskestad nel regime di grande scala.
Reference
Thomas, P.H. (1963). The Size of Flames from Natural Fires. Ninth Symposium (International) on Combustion, The Combustion Institute, pp. 844–859.
Zukoski (1981)
Formula
L / D = 3.3 * Q*^(2/5)
con Q* = Q / (rho_inf * c_p * T_inf * sqrt(g) * D^(5/2))
oppure in forma esplicita:
L = 3.3 * D * [Q / (rho_inf * c_p * T_inf * sqrt(g) * D^(5/2))]^(2/5)
Variables
L : mean flame height [m]
D : diametro del focolare [m]
Q : heat release rate [kW]
Q* : fire Froude number adimensionale [-]
rho_inf : densità aria ambiente [kg/m³] (= 1.2 kg/m³)
c_p : calore specifico aria [kJ/(kg·K)] (= 1.005 kJ/(kg·K))
T_inf : temperatura ambiente [K] (= 293 K, ≈ 20°C)
g : accelerazione di gravità [m/s²] (= 9.81 m/s²)
Inputs and outputs
Input: D [m], Q [kW]
(rho_inf, c_p, T_inf, g sono costanti ambientali)
Output: L [m]
Output type: Mean flame height
Validity and use
- Valida per Q* > ~1 (fuochi di media-grande scala)
- Per Q* molto piccolo la correlazione sovrastima L
- Non considera il termine correttivo "-1.02*D" presente in Heskestad,
pertanto è meno accurata per focolari grandi a basso Q*
Reference
Zukoski, E.E., Kubota, T., and Cetegen, B. (1981). Entrainment in Fire Plumes. Fire Safety Journal, 3(3), 107–121.
NOTA: Alcuni autori attribuiscono la forma L/D = 3.3Q^(2/5) direttamente a Zukoski et al. (1981); altri la usano come approssimazione semplificata di Heskestad per Q* > 1. Verificare il contesto applicativo.
Mccaffrey (1979)
Formula
McCaffrey suddivide il plume in tre zone con leggi di potenza distinte:
ZONA 1 — Fiamma continua (continuous flaming): 0 < z/Q^(2/5) < 0.08
u = 6.84 * (z/Q^(2/5))^(1/2) [m/s] — velocità assiale
DeltaT / T_inf = 0.9 / (z/Q^(2/5))^1 — eccesso di temperatura relativo
ZONA 2 — Fiamma intermittente (intermittent flaming): 0.08 < z/Q^(2/5) < 0.20
u = 1.93 * (z/Q^(2/5))^0 [m/s]
DeltaT / T_inf = 0.9 / (z/Q^(2/5))^1 — continua
(Nota: in zona 2 la velocità è approssimativamente costante)
ZONA 3 — Plume (oltre la fiamma): z/Q^(2/5) > 0.20
u = 1.12 * (z/Q^(2/5))^(-1/3) [m/s]
DeltaT / T_inf = 0.9 / (z/Q^(2/5))^(5/3) (approssimazione)
ALTEZZA DI FIAMMA DERIVATA:
La transizione zona 1 → zona 2 definisce la "continuous flame height":
L_cont = 0.08 * Q^(2/5) [m]
La transizione zona 2 → zona 3 (tip della fiamma intermittente):
L_tip = 0.20 * Q^(2/5) [m]
NOTA: Queste altezze sono derivate dai limiti di zona, non sono formule
dirette per L; vanno usate consapevolmente.
Variables
z : quota verticale dalla base del focolare [m]
Q : heat release rate [kW]
u : velocità assiale del plume [m/s]
DeltaT : eccesso di temperatura rispetto all'ambiente [K]
T_inf : temperatura ambiente [K]
L_cont : altezza della fiamma continua [m]
L_tip : altezza del tip della fiamma intermittente [m]
Inputs and outputs
Input: z [m], Q [kW]
Output: u [m/s], DeltaT [K], L_cont [m], L_tip [m]
Output type: - Continuous flame height (L_cont): quota fino a cui la fiamma è
sempre presente (zona di combustione continua)
- Flame tip height / intermittent flame height (L_tip): quota massima
mediamente raggiunta dalla fiamma
Validity and use
- Sviluppata per focolari circolari di piccola-media dimensione
- Non include esplicitamente il diametro D: valida per D piccolo rispetto a z
- Non adatta per large pool fires o grandi focolari
- Combustibili gassosi (methano, propano) nella validazione originale
Notes
La correlazione di McCaffrey è fondamentale per la comprensione della
struttura del plume. In CFAST e altri modelli di zona viene usata per
il calcolo di temperatura e velocità assiale.
Reference
McCaffrey, B.J. (1979). Purely Buoyant Diffusion Flames: Some Experimental Results. NBSIR 79-1910, National Bureau of Standards (ora NIST), Washington, DC.
Delichatsios (1987)
Formula
L / D = 0.0287 * (T_f / T_inf)^(1/2) * (1/f_s)^(5/3) * Fr_f^(2/5) * D^0
Forma semplificata di uso comune:
L = C_D * Q^(2/5) * D^(1/5)
dove C_D dipende dalle proprietà del combustibile:
C_D = 0.211 * (1/f_s)^(5/3) * (T_f/T_s)^(1/2) * [coefficienti]
Forma adimensionale più esplicita:
L / D = -1.02 + [C_f * (Q_c / D^(5/2))^(2/5)]
(analoga a Heskestad ma con C_f dipendente dal combustibile)
Per combustibili "normali" (idrocarburi):
L = 0.235 * Q_c^(2/5) - 1.02 * D (coincide con Heskestad)
Forma completa con fr_s (stoichiometric mixture fraction):
L / D = f(Q*, f_s, T_ad)
dove T_ad = temperatura di fiamma adiabatica [K]
f_s = stoichiometric mixture fraction [-]
Variables
L : mean flame height [m]
D : diametro del focolare [m]
Q : heat release rate [kW]
Q_c : heat release rate convettivo [kW] (tipicamente 0.7*Q)
f_s : stoichiometric mixture fraction = Y_F,0 / (Y_F,0 + r)
con Y_F,0 = mass fraction combustibile alla sorgente
e r = stoichiometric air-to-fuel mass ratio [-]
T_f : temperatura di fiamma adiabatica [K]
T_inf : temperatura ambiente [K]
Inputs and outputs
Input: D [m], Q [kW] o Q_c [kW], f_s [-], T_f [K], T_inf [K]
Output: L [m]
Output type: Mean flame height (con correzione per le proprietà del combustibile)
Validity and use
- Applicabile a combustibili con diverso contenuto di soot e diverse
temperature di fiamma (H2, CO, idrocarburi, alcoli)
- Utile per confronto tra combustibili diversi nel medesimo scenario
- Richiede la conoscenza di f_s e T_f per il combustibile specifico
Reference
Delichatsios, M.A. (1987). Air Entrainment into Buoyant Jet Flames and Pool Fires. Combustion and Flame, 70(1), 33–46.
Babrauskas (1980)
Formula
L / D = 42 * (m_dot'' / (rho_inf * sqrt(g*D)))^0.61
oppure nella forma a due parametri:
L / D = C * (Q*)^n
con:
Q* = Q / (rho_inf * c_p * T_inf * sqrt(g) * D^(5/2))
C = 42, n = 0.61 (originale Babrauskas)
Forma esplicita in Q:
L = D * 42 * [Q / (rho_inf * c_p * T_inf * sqrt(g) * D^(5/2))]^0.61
Variables
L : mean flame height [m]
D : diametro del pool [m]
m_dot'': mass burning rate per unit area [kg/(m²·s)]
rho_inf: densità aria [kg/m³] (= 1.2 kg/m³)
g : accelerazione di gravità [m/s²] (= 9.81 m/s²)
Q : heat release rate [kW]
Inputs and outputs
Input: D [m], m_dot'' [kg/(m²·s)] oppure Q [kW]
Output: L [m]
Output type: Mean flame height per pool fire (fiamme da pozze di liquido combustibile)
Validity and use
- Specificamente sviluppata per pool fires (pozze di combustibile liquido)
- Range validazione: D da 0.1 m a oltre 10 m
- Combustibili: benzina, kerosene, JP-4 e altri idrocarburi liquidi
- Non adatta per fiamme in pressione o in crossflow
Reference
Babrauskas, V. (1980). Estimating Large Pool Fire Burning Rates. Fire Technology, 16(4), 251–261.
================================================================================
Flame Height Per Line Fires (Elongated Fire Sources)
Formula
Per fuoco lineare (line fire, rapporto L_f/W >> 1):
In regime dominato da bassa velocità di rilascio (weak plume regime):
H / W = 0.034 * Q_L'^(2/3)
In regime dominato da alto rilascio (strong plume regime):
H / W = 0.034 * Q_L'^(1/2) * Fr^(1/2)
con:
Q_L' = Q / (rho_inf * c_p * T_inf * g^(1/2) * W^(3/2)) [adimensionale]
Q' = Q per unità di lunghezza della linea di fuoco [kW/m]
W = larghezza della linea di fuoco [m]
H = altezza di fiamma [m]
Formula pratica comune (Delichatsios 1993):
H = 0.0148 * Q'^(2/3) per Q' < ~300 kW/m
H = 0.034 * Q'^(1/2) per Q' > ~300 kW/m
Variables
H : mean flame height [m]
W : larghezza della sorgente lineare [m]
Q' : heat release rate per unit length [kW/m]
Q_L': fire Froude number per line fire [-]
Inputs and outputs
Input: W [m], Q' [kW/m]
Output: H [m]
Output type: Mean flame height per sorgente lineare
Validity and use
- Rapporto L_f/W > 3–5 (sorgente effettivamente lineare)
- Non adatta per focolari quasi-circolari
Reference
Delichatsios, M.A. (1993). Transition from Momentum to Buoyancy-Controlled Turbulent Jet Diffusion Flames and Flame Height Relationships. Combustion and Flame, 92(4), 349–364.
Comparative summary
================================================================================
RIEPILOGO COMPARATIVO
================================================================================
Correlazione | Tipo output | Input principali | Regime
--------------------|--------------------------|-------------------|------------------
Heskestad (1984) | Mean flame height | D, Q | Generale (tutte scale)
Thomas (1963) | Mean flame height | D, m_dot | Fuochi grandi, legno
Zukoski (1981) | Mean flame height | D, Q | Q* > 1
McCaffrey (1979) | Cont./interm. flame ht. | z, Q | Piccola-media scala
Delichatsios (1987) | Mean flame height | D, Q, f_s, T_f | Combustibili diversi
Babrauskas (1980) | Mean flame height | D, m_dot'', Q | Pool fires
Thomas/Moorhouse | Mean flame height (wind) | D, Q, u_w | Crossflow
Hasemi (1984) | Mean flame height | D, Q, config. | Wall/corner fires
Delichatsios (1993) | Mean flame height | W, Q' | Line fires
Implementation notes
================================================================================
NOTE PER L'IMPLEMENTAZIONE PYTHON
================================================================================
1. Calcolo del diametro equivalente:
D_eq = (4 * A / pi) ** 0.5
dove A = area del focolare in m²
2. Calcolo del Fire Froude number Q*:
Q_star = Q / (rho_inf * c_p * T_inf * g**0.5 * D**2.5)
con Q in kW, D in m
3. Heskestad: per Q* molto piccolo (<0.001) o D molto grande, verificare
che L non risulti negativo; in tal caso L = 0 (fiamma ras terra).
4. McCaffrey: usare la variabile adimensionale xi = z / Q**(2/5)
per determinare in quale zona si trova la quota z.
5. Per la correlazione di Thomas è necessario convertire Q in m_dot:
m_dot = Q / (chi * deltaHc)
Valori tipici: chi = 1.0 (default), deltaHc da tabelle per combustibile.
Main references
================================================================================
RIFERIMENTI BIBLIOGRAFICI PRINCIPALI
================================================================================
[1] Heskestad, G. (1984). Engineering Relations for Fire Plumes.
Fire Safety Journal, 7(1), 25–32.
[2] Heskestad, G. (2016). Fire Plumes, Flame Height, and Air Entrainment.
In: SFPE Handbook of Fire Protection Engineering, 5th ed., Springer.
Chapter 13.
[3] Thomas, P.H. (1963). The Size of Flames from Natural Fires.
9th Symposium on Combustion, Combustion Institute, pp. 844–859.
[4] Zukoski, E.E., Kubota, T., and Cetegen, B. (1981).
Entrainment in Fire Plumes. Fire Safety Journal, 3(3), 107–121.
[5] McCaffrey, B.J. (1979). Purely Buoyant Diffusion Flames: Some
Experimental Results. NBSIR 79-1910, NBS, Washington, DC.
[6] Delichatsios, M.A. (1987). Air Entrainment into Buoyant Jet Flames
and Pool Fires. Combustion and Flame, 70(1), 33–46.
[7] Babrauskas, V. (1980). Estimating Large Pool Fire Burning Rates.
Fire Technology, 16(4), 251–261.
[8] Moorhouse, J. (1982). Scaling Criteria for Pool Fires Derived from
Large Scale Experiments. IChemE Symposium Series No. 71, pp. 165–179.
[9] Hasemi, Y. and Tokunaga, T. (1984). Flame Geometry Effects on the
Buoyant Plumes from Turbulent Diffusion Flames.
Combustion Science and Technology, 40(1–4), 1–17.
[10] Delichatsios, M.A. (1993). Transition from Momentum to Buoyancy-
Controlled Turbulent Jet Diffusion Flames and Flame Height Relationships.
Combustion and Flame, 92(4), 349–364.
[11] SFPE Handbook of Fire Protection Engineering, 5th edition (2016).
Springer, New York. (riferimento generale per tutte le correlazioni)
================================================================================
END OF DOCUMENT
================================================================================