MQH correlation - Revision history

Jan: Add 4 watermarked SVG illustrations from the Expert visual plan

2026-11-02T16:18:59Z

Add 4 watermarked SVG illustrations from the Expert visual plan

← Older revision		Revision as of 17:18, 2 November 2026
Line 13:		Line 13:
	Van Rossum describes the same idea as an energy balance on the smoke layer. The convective part of the heat release rate enters the layer, and energy leaves with the outflowing smoke and by heat transfer to the walls and ceiling. Whatever remains changes the temperature of the layer.<ref name="VR-3.7">Van Rossum (2027), §3.7.</ref> For a quasi-steady fire, McCaffrey, Quintiere and Harkleroad state that the heat released by the fire equals the heat carried out with the hot gas plus the heat lost to the boundaries:<ref name="WT-3-175" />		Van Rossum describes the same idea as an energy balance on the smoke layer. The convective part of the heat release rate enters the layer, and energy leaves with the outflowing smoke and by heat transfer to the walls and ceiling. Whatever remains changes the temperature of the layer.<ref name="VR-3.7">Van Rossum (2027), §3.7.</ref> For a quasi-steady fire, McCaffrey, Quintiere and Harkleroad state that the heat released by the fire equals the heat carried out with the hot gas plus the heat lost to the boundaries:<ref name="WT-3-175" />
	<math display="block">\dot{Q} = \dot{m}_g c_p \left( T_g - T_\infty \right) + \dot{q}_\mathrm{loss}, \qquad \dot{q}_\mathrm{loss} = h_k A_T \left( T_g - T_\infty \right)</math>		<math display="block">\dot{Q} = \dot{m}_g c_p \left( T_g - T_\infty \right) + \dot{q}_\mathrm{loss}, \qquad \dot{q}_\mathrm{loss} = h_k A_T \left( T_g - T_\infty \right)</math>
			[[File:fse-mqh-correlation-v1.svg\|thumb\|upright=1.5\|alt=Cross-section of a room fire with a hot upper layer, hot gas leaving through the top of a door, air entering below and heat flowing into walls.\|Energy balance on the hot upper layer: the heat released by the fire leaves with the hot gas through the opening or is lost to the walls, ceiling and floor. Based on Walton and Thomas (2002), pp. 3-172 and 3-175.]]

	The heat loss is written as an effective heat transfer coefficient multiplied by the total area of the enclosing surfaces and the temperature rise.<ref name="WT-3-175" /> Solving for the temperature rise gives two dimensionless groups. The first compares the heat released with the heat that the outflowing gas can carry away. The second compares the heat lost to the boundaries with the heat carried away by the gas.<ref name="WT-3-176" /> The outflow depends on the height of the neutral plane in the opening, which in turn depends mainly on the gas temperature, the heat release rate and the geometry. The outflow is therefore replaced by the term <math>\sqrt{g}\rho_\infty A_0 \sqrt{H_0}</math>, and the remaining effects are absorbed into the fitted constants.<ref name="WT-3-176" />		The heat loss is written as an effective heat transfer coefficient multiplied by the total area of the enclosing surfaces and the temperature rise.<ref name="WT-3-175" /> Solving for the temperature rise gives two dimensionless groups. The first compares the heat released with the heat that the outflowing gas can carry away. The second compares the heat lost to the boundaries with the heat carried away by the gas.<ref name="WT-3-176" /> The outflow depends on the height of the neutral plane in the opening, which in turn depends mainly on the gas temperature, the heat release rate and the geometry. The outflow is therefore replaced by the term <math>\sqrt{g}\rho_\infty A_0 \sqrt{H_0}</math>, and the remaining effects are absorbed into the fitted constants.<ref name="WT-3-176" />
Line 67:		Line 68:

	The [[thermal inertia]] of the lining is the product <math>k\rho c</math>. NUREG-1805 lists typical values, for example 0.18 (kW/(m²·K))²·s for gypsum board and 2.9 (kW/(m²·K))²·s for concrete.<ref name="IS-2-11">Iqbal and Salley (2004), p. 2-11.</ref> A lining with a low thermal inertia gives a low coefficient, and the same fire then produces a hotter layer. If walls and ceiling are made of different materials, an area-weighted average of <math>h_k</math> is used.<ref name="WT-3-176" />		The [[thermal inertia]] of the lining is the product <math>k\rho c</math>. NUREG-1805 lists typical values, for example 0.18 (kW/(m²·K))²·s for gypsum board and 2.9 (kW/(m²·K))²·s for concrete.<ref name="IS-2-11">Iqbal and Salley (2004), p. 2-11.</ref> A lining with a low thermal inertia gives a low coefficient, and the same fire then produces a hotter layer. If walls and ceiling are made of different materials, an area-weighted average of <math>h_k</math> is used.<ref name="WT-3-176" />
			[[File:fse-mqh-correlation-v3.svg\|thumb\|upright=1.5\|alt=Line chart: heat transfer coefficient of a gypsum plaster lining falls with time, then drops to a constant value at 161 seconds.\|Effective heat transfer coefficient of 16 mm gypsum plaster against time; the thermal penetration time is about 161 s. Based on Walton and Thomas (2002), pp. 3-176–3-177.]]

	== Worked example ==		== Worked example ==
Line 74:		Line 76:

	The full correlation then gives temperature rises of 227 K after 10 s, 306 K after 60 s and 455 K after 600 s. The handbook reports 227 K, 307 K and 453 K, because it rounds the intermediate terms.<ref name="WT-3-177" /> After ten minutes the layer is at about 480 °C. This is below the 500 °C often used as a flashover criterion and within the 600 °C range of the data. The 750 kW fire is also well below the stoichiometric heat release rate for this door, 1,500 × <math>A_0 \sqrt{H_0}</math> ≈ 2,170 kW, which is the most that the incoming air can burn completely.<ref name="WT-3-184">Walton and Thomas (2002), p. 3-184.</ref>		The full correlation then gives temperature rises of 227 K after 10 s, 306 K after 60 s and 455 K after 600 s. The handbook reports 227 K, 307 K and 453 K, because it rounds the intermediate terms.<ref name="WT-3-177" /> After ten minutes the layer is at about 480 °C. This is below the 500 °C often used as a flashover criterion and within the 600 °C range of the data. The 750 kW fire is also well below the stoichiometric heat release rate for this door, 1,500 × <math>A_0 \sqrt{H_0}</math> ≈ 2,170 kW, which is the most that the incoming air can burn completely.<ref name="WT-3-184">Walton and Thomas (2002), p. 3-184.</ref>
			[[File:fse-mqh-correlation-v2.svg\|thumb\|upright=1.5\|alt=Line chart: upper layer temperature rises with heat release rate, faster for warm walls after 10 minutes than after 60 seconds.\|Upper layer temperature from the MQH correlation for the example room, 60 s after ignition and after the lining has heated through. Based on Walton and Thomas (2002), pp. 3-176–3-177.]]

	== Heat release rate for flashover ==		== Heat release rate for flashover ==
Line 81:		Line 84:

	A simpler method by Babrauskas gives about 1,090 kW for the same room, using 750 × <math>A_0 \sqrt{H_0}</math>.<ref name="WT-3-185" /> Babrauskas compared the methods and found similar results over the range of room sizes of most interest. The MQH method departs from the others for small rooms, and all the methods represent the test data conservatively.<ref name="WT-3-186" /> NUREG-1805 adds that the minimum heat release rate for flashover grows with the size of the compartment. With too little ventilation flashover cannot occur, while too much ventilation dilutes and cools the smoke.<ref name="IS-13-5">Iqbal and Salley (2004), p. 13-5.</ref>		A simpler method by Babrauskas gives about 1,090 kW for the same room, using 750 × <math>A_0 \sqrt{H_0}</math>.<ref name="WT-3-185" /> Babrauskas compared the methods and found similar results over the range of room sizes of most interest. The MQH method departs from the others for small rooms, and all the methods represent the test data conservatively.<ref name="WT-3-186" /> NUREG-1805 adds that the minimum heat release rate for flashover grows with the size of the compartment. With too little ventilation flashover cannot occur, while too much ventilation dilutes and cools the smoke.<ref name="IS-13-5">Iqbal and Salley (2004), p. 13-5.</ref>
			[[File:fse-mqh-correlation-v4.svg\|thumb\|upright=1.5\|alt=Bar chart for one room: design fire 750 kW, flashover 860 kW by MQH and 1,087 kW by Babrauskas, stoichiometric limit 2,173 kW.\|Heat release rates for the example room: the design fire, the rate for flashover by two methods and the stoichiometric limit of the door. Based on Walton and Thomas (2002), pp. 3-184–3-186.]]

	The result depends on the chosen [[flashover criteria]]. NUREG-1805 uses a temperature rise of 500 °C, sometimes 600 °C, or a radiant heat flux of 20 kW/m² at the floor.<ref name="IS-13-3">Iqbal and Salley (2004), pp. 13-3–13-4.</ref> Van Rossum links a layer temperature of 500 °C to a radiation of about 20 kW/m² at the floor, at which many materials pyrolyse rapidly.<ref name="VR-4.7">Van Rossum (2027), §4.7.</ref> In Australia, the data sheets of the Fire Safety Verification Method cite the SFPE engineering standard S.01. This standard bases its approach on the MQH and Thomas methods and uses 600 °C or 20 kW/m² as flashover criteria in fire models.<ref name="ABCB-12">SFPE (2011), as cited in Australian Building Codes Board (2020a), Data Sheet B1, p. 12.</ref>		The result depends on the chosen [[flashover criteria]]. NUREG-1805 uses a temperature rise of 500 °C, sometimes 600 °C, or a radiant heat flux of 20 kW/m² at the floor.<ref name="IS-13-3">Iqbal and Salley (2004), pp. 13-3–13-4.</ref> Van Rossum links a layer temperature of 500 °C to a radiation of about 20 kW/m² at the floor, at which many materials pyrolyse rapidly.<ref name="VR-4.7">Van Rossum (2027), §4.7.</ref> In Australia, the data sheets of the Fire Safety Verification Method cite the SFPE engineering standard S.01. This standard bases its approach on the MQH and Thomas methods and uses 600 °C or 20 kW/m² as flashover criteria in fire models.<ref name="ABCB-12">SFPE (2011), as cited in Australian Building Codes Board (2020a), Data Sheet B1, p. 12.</ref>

Mark: FSE wiki: corrections after independent fact-check (written by Claude Opus 5.5)

2026-11-01T23:22:30Z

FSE wiki: corrections after independent fact-check (written by Claude Opus 5.5)

← Older revision		Revision as of 00:22, 2 November 2026
Line 1:		Line 1:
	The '''MQH correlation''' is a simple formula for estimating the average temperature of the hot smoke layer in a room fire with a door or window open to the outside. It is named after McCaffrey, Quintiere and Harkleroad, who developed it in 1981 by fitting a simple energy balance to data from about 100 experimental fires.<ref name="IS-2-8">Iqbal and Salley (2004), p. 2-8.</ref> The formula gives the rise of the upper layer temperature above ambient from the [[heat release rate]] of the fire, the size of the opening and the heat lost to the walls, ceiling and floor.<ref name="WT-3-176">Walton and Thomas (2002), p. 3-176.</ref> It belongs to the methods for compartment fire temperatures that are most widely accepted in fire protection engineering.<ref name="WT-3-175">Walton and Thomas (2002), p. 3-175.</ref><ref name="IS-2-8" />		The '''MQH correlation''' is a simple formula for estimating the average temperature of the hot smoke layer in a room fire with a door or window open to the outside. It is named after McCaffrey, Quintiere and Harkleroad, who developed it in 1981 by fitting a simple energy balance to data from about 100 experimental fires.<ref name="IS-2-8">Iqbal and Salley (2004), p. 2-8.</ref> The formula gives the rise of the upper layer temperature above ambient from the [[heat release rate]] of the fire, the size of the opening and the heat lost to the walls, ceiling and floor.<ref name="WT-3-176">Walton and Thomas (2002), p. 3-176.</ref> It belongs to the methods for compartment fire temperatures that are most widely accepted in fire protection engineering.<ref name="WT-3-175">Walton and Thomas (2002), p. 3-175.</ref><ref name="IS-2-8" />

	The correlation is one of the standard [[hand calculations]] of fire safety engineering.<ref name="WT-3-175" /> It is used to estimate the gas temperature and, in a rearranged form, the heat release rate needed for flashover.<ref name="WT-3-185">Walton and Thomas (2002), p. 3-185.</ref> It is valid for upper layer temperatures up to about 600 °C and predicts room averages, not local temperatures near the fire.<ref name="WT-3-176" /><ref name="WT-3-175" /> In a validation study against full-scale experiments, it overpredicted the layer temperature by about 17 % on average.<ref name="NRC-6-57">Stroup et al. (2016), p. 6-57.</ref>		The correlation is one of the standard [[hand calculations]] of fire safety engineering.<ref name="WT-3-175" /> It is used to estimate the gas temperature and, in a rearranged form, the heat release rate needed for flashover.<ref name="WT-3-185">Walton and Thomas (2002), p. 3-185.</ref> It is valid for upper layer temperatures up to about 600 °C and predicts room averages, not local temperatures near the fire.<ref name="WT-3-176" /><ref name="WT-3-175" /> In a validation study against full-scale experiments, it overpredicted the rise of the layer temperature above ambient by about 17 % on average.<ref name="NRC-6-57">Stroup et al. (2016), p. 6-57.</ref>

	== Origin and data ==		== Origin and data ==
Line 22:		Line 22:
	<math display="block">\Delta T_g = 480 \left( \frac{\dot{Q}}{\sqrt{g}\, c_p \rho_\infty T_\infty A_0 \sqrt{H_0}} \right)^{2/3} \left( \frac{h_k A_T}{\sqrt{g}\, c_p \rho_\infty A_0 \sqrt{H_0}} \right)^{-1/3}</math>		<math display="block">\Delta T_g = 480 \left( \frac{\dot{Q}}{\sqrt{g}\, c_p \rho_\infty T_\infty A_0 \sqrt{H_0}} \right)^{2/3} \left( \frac{h_k A_T}{\sqrt{g}\, c_p \rho_\infty A_0 \sqrt{H_0}} \right)^{-1/3}</math>

	With standard ambient values of <math>g</math> = 9.8 m/s², <math>c_p</math> = 1.05 kJ/(kg·K), <math>\rho_\infty</math> = 1.2 kg/m³ and <math>T_\infty</math> = 295 K, the formula reduces to the form ~~used in most hand calculations~~:<ref name="WT-3-176" /><ref name="IS-2-9">Iqbal and Salley (2004), p. 2-9.</ref>		With standard ambient values of <math>g</math> = 9.8 m/s², <math>c_p</math> = 1.05 kJ/(kg·K), <math>\rho_\infty</math> = 1.2 kg/m³ and <math>T_\infty</math> = 295 K, the formula reduces to the simpler form that NUREG-1805 also uses:<ref name="WT-3-176" /><ref name="IS-2-9">Iqbal and Salley (2004), p. 2-9.</ref>
	<math display="block">\Delta T_g = 6.85 \left( \frac{\dot{Q}^2}{A_0 \sqrt{H_0}\, h_k A_T} \right)^{1/3}</math>		<math display="block">\Delta T_g = 6.85 \left( \frac{\dot{Q}^2}{A_0 \sqrt{H_0}\, h_k A_T} \right)^{1/3}</math>

Line 85:		Line 85:

	== Accuracy and validation ==		== Accuracy and validation ==
	The NRC and the Electric Power Research Institute (EPRI) compared the correlation with full-scale experiments in naturally ventilated rooms. All results were within or above the experimental uncertainty, with no sign of systematic underprediction.<ref name="NRC-6-6">Stroup et al. (2016), p. 6-6.</ref> ~~The~~ bias factor, which shows the average overprediction or underprediction, was 1.17 with a relative standard deviation of 0.15.<ref name="NRC-6-1">Stroup et al. (2016), p. 6-1.</ref><ref name="NRC-6-57" /> For comparison, the zone model [[CFAST]] scored 1.21 and 0.38, and the [[Fire Dynamics Simulator]] (FDS) scored 1.02 and 0.07.<ref name="NRC-6-57" />		The NRC and the Electric Power Research Institute (EPRI) compared the correlation with full-scale experiments in naturally ventilated rooms. All results were within or above the experimental uncertainty, with no sign of systematic underprediction.<ref name="NRC-6-6">Stroup et al. (2016), p. 6-6.</ref> For the rise of the layer temperature, the bias factor, which shows the average overprediction or underprediction, was 1.17 with a relative standard deviation of 0.15.<ref name="NRC-6-1">Stroup et al. (2016), p. 6-1.</ref><ref name="NRC-6-57" /> For comparison, the zone model [[CFAST]] scored 1.21 and 0.38, and the [[Fire Dynamics Simulator]] (FDS) scored 1.02 and 0.07.<ref name="NRC-6-57" />

	The power 2/3 also governs how errors in the input carry through. A relative error in the heat release rate gives two thirds of that relative error in the temperature rise, so 7.5 % uncertainty in the heat release rate gives about 5 % in temperature.<ref name="NRC-3-35">Stroup et al. (2016), p. 3-35.</ref> In the investigation of the World Trade Center fires, the National Institute of Standards and Technology (NIST) used the same reasoning. A 5 % uncertainty in the measured heat release rate corresponded to 3.3 % in the temperature rise, or roughly ±20 °C at 600 °C. The observed differences between simulations and measurements were 5 to 20 °C.<ref name="NIST-24">McGrattan et al. (2005), p. 24.</ref>		The power 2/3 also governs how errors in the input carry through. A relative error in the heat release rate gives two thirds of that relative error in the temperature rise, so 7.5 % uncertainty in the heat release rate gives about 5 % in temperature.<ref name="NRC-3-35">Stroup et al. (2016), p. 3-35.</ref> In the investigation of the World Trade Center fires, the National Institute of Standards and Technology (NIST) used the same reasoning. A 5 % uncertainty in the measured heat release rate corresponded to 3.3 % in the temperature rise, or roughly ±20 °C at 600 °C. The observed differences between simulations and measurements were 5 to 20 °C.<ref name="NIST-24">McGrattan et al. (2005), pp. 24–25.</ref>

	NIST also noted that the correlation rests on test data below 600 °C. In a test with a layer temperature of about 750 °C, a 10 % higher heat release rate raised the simulated temperature by 9 %. The correlation implies 6.7 %, and NIST attributed the difference to its simplifications and its limited range.<ref name="NIST-28">McGrattan et al. (2005), p. 28.</ref>		NIST also noted that the correlation rests on test data below 600 °C. In a test with a layer temperature of about 750 °C, a 10 % higher heat release rate raised the simulated temperature by 9 %. The correlation implies 6.7 %, and NIST attributed the difference to its simplifications and its limited range.<ref name="NIST-28">McGrattan et al. (2005), p. 28.</ref>

Mark: New FSE wiki article: MQH correlation (written by Claude Opus 5.5)

2026-11-01T22:12:11Z

New FSE wiki article: MQH correlation (written by Claude Opus 5.5)

New page

The '''MQH correlation''' is a simple formula for estimating the average temperature of the hot smoke layer in a room fire with a door or window open to the outside. It is named after McCaffrey, Quintiere and Harkleroad, who developed it in 1981 by fitting a simple energy balance to data from about 100 experimental fires.<ref name="IS-2-8">Iqbal and Salley (2004), p. 2-8.</ref> The formula gives the rise of the upper layer temperature above ambient from the [[heat release rate]] of the fire, the size of the opening and the heat lost to the walls, ceiling and floor.<ref name="WT-3-176">Walton and Thomas (2002), p. 3-176.</ref> It belongs to the methods for compartment fire temperatures that are most widely accepted in fire protection engineering.<ref name="WT-3-175">Walton and Thomas (2002), p. 3-175.</ref><ref name="IS-2-8" />

The correlation is one of the standard [[hand calculations]] of fire safety engineering.<ref name="WT-3-175" /> It is used to estimate the gas temperature and, in a rearranged form, the heat release rate needed for flashover.<ref name="WT-3-185">Walton and Thomas (2002), p. 3-185.</ref> It is valid for upper layer temperatures up to about 600 °C and predicts room averages, not local temperatures near the fire.<ref name="WT-3-176" /><ref name="WT-3-175" /> In a validation study against full-scale experiments, it overpredicted the layer temperature by about 17 % on average.<ref name="NRC-6-57">Stroup et al. (2016), p. 6-57.</ref>

== Origin and data ==
McCaffrey, Quintiere and Harkleroad published the correlation in 1981, and Walton and Thomas later reported it in the ''SFPE Handbook of Fire Protection Engineering''.<ref name="IS-2-8" /> The constants were fitted to data from over 100 experimental fires, both steady and transient. The fuels included cellulosic materials, synthetic polymers and gaseous hydrocarbons.<ref name="WT-3-176" /> The test compartments were 0.3 to 2.7 m high, with floor areas of 0.14 to 12.0 m² and a variety of window and door sizes.<ref name="WT-3-176" />

NUREG-1805, a set of fire hazard calculation methods of the United States Nuclear Regulatory Commission (NRC), describes the data as 100 fires from eight test series. The rooms were of conventional size, and the temperature rises ranged from 20 to 600 °C.<ref name="IS-2-8" /> In these tests the fire was placed away from the walls, and most of the fuel was near the centre of the room.<ref name="IS-2-8" /><ref name="WT-3-176" /> NUREG-1805 summarises the main trend: the larger the heat release rate and the smaller the opening, the higher the upper layer temperature.<ref name="IS-2-8" />

== Energy balance ==
The correlation rests on the two-layer picture of a room fire. The fire plume carries hot gases to the ceiling, where they form a [[hot gas layer]] that deepens with time. A fairly sharp interface separates this layer from the cooler air below. Once the layer reaches the top of a door or window, hot gas flows out through the upper part of the opening and outside air flows in below. Walton and Thomas call this description a two-layer or [[zone model]].<ref name="WT-3-172">Walton and Thomas (2002), p. 3-172.</ref>

Van Rossum describes the same idea as an energy balance on the smoke layer. The convective part of the heat release rate enters the layer, and energy leaves with the outflowing smoke and by heat transfer to the walls and ceiling. Whatever remains changes the temperature of the layer.<ref name="VR-3.7">Van Rossum (2027), §3.7.</ref> For a quasi-steady fire, McCaffrey, Quintiere and Harkleroad state that the heat released by the fire equals the heat carried out with the hot gas plus the heat lost to the boundaries:<ref name="WT-3-175" />
<math display="block">\dot{Q} = \dot{m}_g c_p \left( T_g - T_\infty \right) + \dot{q}_\mathrm{loss}, \qquad \dot{q}_\mathrm{loss} = h_k A_T \left( T_g - T_\infty \right)</math>

The heat loss is written as an effective heat transfer coefficient multiplied by the total area of the enclosing surfaces and the temperature rise.<ref name="WT-3-175" /> Solving for the temperature rise gives two dimensionless groups. The first compares the heat released with the heat that the outflowing gas can carry away. The second compares the heat lost to the boundaries with the heat carried away by the gas.<ref name="WT-3-176" /> The outflow depends on the height of the neutral plane in the opening, which in turn depends mainly on the gas temperature, the heat release rate and the geometry. The outflow is therefore replaced by the term <math>\sqrt{g}\rho_\infty A_0 \sqrt{H_0}</math>, and the remaining effects are absorbed into the fitted constants.<ref name="WT-3-176" />

Measurements show how the heat divides between these terms. In experimental wood crib fires in a compartment with a floor area of 29 m², 47 to 65 % of the heat left with the hot gases through the windows. Another 15 to 30 % went into the walls, ceiling and floor, and the rest returned to the fuel or radiated out through the windows. With smaller windows, a larger share of the heat went into the enclosing surfaces.<ref name="WT-3-175" />

== The correlation ==
Fitting the two groups to the test data gave a power law with a constant of 480 and exponents of 2/3 and −1/3:<ref name="WT-3-176" />
<math display="block">\Delta T_g = 480 \left( \frac{\dot{Q}}{\sqrt{g}\, c_p \rho_\infty T_\infty A_0 \sqrt{H_0}} \right)^{2/3} \left( \frac{h_k A_T}{\sqrt{g}\, c_p \rho_\infty A_0 \sqrt{H_0}} \right)^{-1/3}</math>

With standard ambient values of <math>g</math> = 9.8 m/s², <math>c_p</math> = 1.05 kJ/(kg·K), <math>\rho_\infty</math> = 1.2 kg/m³ and <math>T_\infty</math> = 295 K, the formula reduces to the form used in most hand calculations:<ref name="WT-3-176" /><ref name="IS-2-9">Iqbal and Salley (2004), p. 2-9.</ref>
<math display="block">\Delta T_g = 6.85 \left( \frac{\dot{Q}^2}{A_0 \sqrt{H_0}\, h_k A_T} \right)^{1/3}</math>

{| class="wikitable"
! Symbol !! Description !! Unit
|-
| <math>\Delta T_g</math> || rise of the upper layer temperature above ambient, <math>T_g - T_\infty</math> || K
|-
| <math>\dot{Q}</math> || heat release rate of the fire || kW
|-
| <math>A_0</math> || area of the opening || m²
|-
| <math>H_0</math> || height of the opening || m
|-
| <math>A_T</math> || total area of the enclosing surfaces, without the openings || m²
|-
| <math>h_k</math> || effective heat transfer coefficient of the boundaries || kW/(m²·K)
|-
| <math>g</math> || acceleration due to gravity || m/s²
|-
| <math>c_p</math> || specific heat of the gas || kJ/(kg·K)
|-
| <math>\rho_\infty</math> || ambient air density || kg/m³
|-
| <math>T_\infty</math> || ambient temperature || K
|-
| <math>k</math> || thermal conductivity of the lining || kW/(m·K)
|-
| <math>\rho</math>, <math>c</math> || density and specific heat of the lining || kg/m³, kJ/(kg·K)
|-
| <math>\delta</math> || thickness of the lining || m
|-
| <math>t</math>, <math>t_p</math> || time since ignition, thermal penetration time || s
|}

The term <math>A_0 \sqrt{H_0}</math> is the [[ventilation factor]] of the opening.<ref name="WT-3-174">Walton and Thomas (2002), p. 3-174.</ref> For several openings, the values of <math>A_0 \sqrt{H_0}</math> are added together. The area <math>A_T</math> covers the walls, ceiling and floor, with the area of the openings subtracted.<ref name="IS-2-9" /> Because the temperature rise grows with the heat release rate to the power 2/3, doubling the heat release rate raises it by a factor of about 1.6. Halving the ventilation factor raises it by a factor of about 1.26, the cube root of 2.

== Heat losses to the boundaries ==
The effective heat transfer coefficient depends on the lining and on how long the fire has burnt. Heat first soaks into the lining, and only later flows through it to the back face. The time this takes is the thermal penetration time:<ref name="WT-3-176" /><ref name="IS-2-10">Iqbal and Salley (2004), p. 2-10.</ref>
<math display="block">t_p = \frac{\rho c}{k} \left( \frac{\delta}{2} \right)^2</math>

After this time, heat conduction through the lining is close to steady, and the coefficient is the conductivity divided by the thickness. Before it, the lining behaves like a thick solid that keeps absorbing heat, and the coefficient falls with time:<ref name="WT-3-176" /><ref name="IS-2-10" />
<math display="block">h_k = \frac{k}{\delta} \quad \text{for } t > t_p, \qquad h_k = \left( \frac{k \rho c}{t} \right)^{1/2} \quad \text{for } t \le t_p</math>

The [[thermal inertia]] of the lining is the product <math>k\rho c</math>. NUREG-1805 lists typical values, for example 0.18 (kW/(m²·K))²·s for gypsum board and 2.9 (kW/(m²·K))²·s for concrete.<ref name="IS-2-11">Iqbal and Salley (2004), p. 2-11.</ref> A lining with a low thermal inertia gives a low coefficient, and the same fire then produces a hotter layer. If walls and ceiling are made of different materials, an area-weighted average of <math>h_k</math> is used.<ref name="WT-3-176" />

== Worked example ==
Walton and Thomas apply the method to a room of 3 m by 3 m with a height of 2.4 m. The room has a door 1.8 m high and 0.6 m wide, and a steady fire of 750 kW.<ref name="WT-3-177">Walton and Thomas (2002), p. 3-177.</ref> The lining is 16 mm of gypsum plaster on metal lath with a density of 1,440 kg/m³. Its conductivity is 0.48 × 10⁻³ kW/(m·K) and its specific heat 0.84 kJ/(kg·K). The ambient air is at 27 °C (300 K), with a density of 1.18 kg/m³ and a specific heat of 1.0 kJ/(kg·K).<ref name="WT-3-177" />

The door has an area of 1.08 m², which gives a ventilation factor of 1.08 × √1.8 ≈ 1.45 m<sup>5/2</sup>. The enclosing surfaces add up to 28.8 m² of wall, 9 m² of floor and 9 m² of ceiling, minus the door, which gives 45.72 m². The thermal penetration time is (1,440 × 0.84 / 0.00048) × 0.008², or about 161 s.<ref name="WT-3-177" /> At 10 s and 60 s the lining is still absorbing heat, so <math>h_k</math> is 0.24 and 0.098 kW/(m²·K). At 600 s the steady value applies: 0.00048 / 0.016 = 0.030 kW/(m²·K).<ref name="WT-3-177" />

The full correlation then gives temperature rises of 227 K after 10 s, 306 K after 60 s and 455 K after 600 s. The handbook reports 227 K, 307 K and 453 K, because it rounds the intermediate terms.<ref name="WT-3-177" /> After ten minutes the layer is at about 480 °C. This is below the 500 °C often used as a flashover criterion and within the 600 °C range of the data. The 750 kW fire is also well below the stoichiometric heat release rate for this door, 1,500 × <math>A_0 \sqrt{H_0}</math> ≈ 2,170 kW, which is the most that the incoming air can burn completely.<ref name="WT-3-184">Walton and Thomas (2002), p. 3-184.</ref>

== Heat release rate for flashover ==
Rearranging the correlation gives the heat release rate needed to reach a chosen layer temperature. Walton and Thomas take a temperature rise of 500 °C above an ambient of 295 K as the condition for flashover, with standard air properties.<ref name="WT-3-185" /> This gives:
<math display="block">\dot{Q}_\mathrm{fo} = 610 \left( h_k A_T A_0 \sqrt{H_0} \right)^{1/2}</math>
where <math>\dot{Q}_\mathrm{fo}</math> is the heat release rate for flashover in kilowatts. Starting from the simplified form gives a coefficient of 620 instead, and either value is acceptable within the accuracy of the expression.<ref name="WT-3-186">Walton and Thomas (2002), p. 3-186.</ref> For the example room with <math>h_k</math> = 0.03 kW/(m²·K), the result is about 860 kW.<ref name="WT-3-186" />

A simpler method by Babrauskas gives about 1,090 kW for the same room, using 750 × <math>A_0 \sqrt{H_0}</math>.<ref name="WT-3-185" /> Babrauskas compared the methods and found similar results over the range of room sizes of most interest. The MQH method departs from the others for small rooms, and all the methods represent the test data conservatively.<ref name="WT-3-186" /> NUREG-1805 adds that the minimum heat release rate for flashover grows with the size of the compartment. With too little ventilation flashover cannot occur, while too much ventilation dilutes and cools the smoke.<ref name="IS-13-5">Iqbal and Salley (2004), p. 13-5.</ref>

The result depends on the chosen [[flashover criteria]]. NUREG-1805 uses a temperature rise of 500 °C, sometimes 600 °C, or a radiant heat flux of 20 kW/m² at the floor.<ref name="IS-13-3">Iqbal and Salley (2004), pp. 13-3–13-4.</ref> Van Rossum links a layer temperature of 500 °C to a radiation of about 20 kW/m² at the floor, at which many materials pyrolyse rapidly.<ref name="VR-4.7">Van Rossum (2027), §4.7.</ref> In Australia, the data sheets of the Fire Safety Verification Method cite the SFPE engineering standard S.01. This standard bases its approach on the MQH and Thomas methods and uses 600 °C or 20 kW/m² as flashover criteria in fire models.<ref name="ABCB-12">SFPE (2011), as cited in Australian Building Codes Board (2020a), Data Sheet B1, p. 12.</ref>

== Accuracy and validation ==
The NRC and the Electric Power Research Institute (EPRI) compared the correlation with full-scale experiments in naturally ventilated rooms. All results were within or above the experimental uncertainty, with no sign of systematic underprediction.<ref name="NRC-6-6">Stroup et al. (2016), p. 6-6.</ref> The bias factor, which shows the average overprediction or underprediction, was 1.17 with a relative standard deviation of 0.15.<ref name="NRC-6-1">Stroup et al. (2016), p. 6-1.</ref><ref name="NRC-6-57" /> For comparison, the zone model [[CFAST]] scored 1.21 and 0.38, and the [[Fire Dynamics Simulator]] (FDS) scored 1.02 and 0.07.<ref name="NRC-6-57" />

The power 2/3 also governs how errors in the input carry through. A relative error in the heat release rate gives two thirds of that relative error in the temperature rise, so 7.5 % uncertainty in the heat release rate gives about 5 % in temperature.<ref name="NRC-3-35">Stroup et al. (2016), p. 3-35.</ref> In the investigation of the World Trade Center fires, the National Institute of Standards and Technology (NIST) used the same reasoning. A 5 % uncertainty in the measured heat release rate corresponded to 3.3 % in the temperature rise, or roughly ±20 °C at 600 °C. The observed differences between simulations and measurements were 5 to 20 °C.<ref name="NIST-24">McGrattan et al. (2005), p. 24.</ref>

NIST also noted that the correlation rests on test data below 600 °C. In a test with a layer temperature of about 750 °C, a 10 % higher heat release rate raised the simulated temperature by 9 %. The correlation implies 6.7 %, and NIST attributed the difference to its simplifications and its limited range.<ref name="NIST-28">McGrattan et al. (2005), p. 28.</ref>

== Limitations ==
Walton and Thomas list the limits given by McCaffrey and co-workers. The correlation holds up to an upper layer temperature of about 600 °C. It applies to steady and growing fires, as long as heat conduction into the walls is the main transient process.<ref name="WT-3-176" /> It does not apply to rapidly growing fires in large enclosures, where much of the growth happens before smoke leaves the room. The heat release rate and the thermal properties of the lining must be known, and the data contain few ventilation-controlled fires and no combustible linings.<ref name="WT-3-176" />

NUREG-1805 adds several conditions. The correlation predicts average temperatures and should not be used for the activation of detectors or sprinklers, or for heating by direct flame contact.<ref name="IS-2-23">Iqbal and Salley (2004), p. 2-23.</ref> It treats heat flow into the boundaries as one-dimensional and uses one heat transfer coefficient for the whole room. It is not valid with the coefficient 6.85 for a fire against a wall or in a corner, and it does not allow an opening in the ceiling.<ref name="IS-2-23" /> The validation study also notes that the correlation assumes a single material for ceiling, walls and floor.<ref name="NRC-6-6" />

The size of the room is less critical than it may seem. NUREG-1805 advises caution for large compartments, but it also reports Quintiere's view that the correlation works for any room size because it is dimensionless.<ref name="IS-2-23" /><ref name="IS-2-24">Iqbal and Salley (2004), p. 2-24.</ref> According to the same report, Quintiere also gave an upper limit for the heat release rate, based on about 3,000 kJ per kilogram and the air inflow term 0.5 <math>A_0 \sqrt{H_0}</math>.<ref name="IS-2-24" /> Walton and Thomas use the same two numbers for the heat released per kilogram of air and for the air inflow of a room near flashover.<ref name="WT-3-184" /> Rooms with several connected spaces, ceiling vents or forced ventilation need other correlations or a computer fire model.<ref name="NRC-6-6" /><ref name="WT-3-175" />

=== Related correlations ===
Foote, Pagni and Alvares (FPA) applied the same energy balance to rooms with mechanical ventilation. Their correlation uses the measured ventilation flow instead of the flow through an opening, with a coefficient of 0.63 and exponents of 0.72 and −0.36.<ref name="WT-3-177" /> It was fitted to well-ventilated tests in a compartment of 6 m by 4 m and 4.5 m high, with ventilation rates of 110 to 325 g/s.<ref name="WT-3-177" /> For closed rooms, NUREG-1805 gives the method of Beyler, which assumes that the room leaks enough to prevent a pressure build-up.<ref name="IS-2-12">Iqbal and Salley (2004), p. 2-12.</ref> In the validation study of the NRC and EPRI, the FPA correlation had a bias factor of 1.29 and the Beyler method 1.04.<ref name="NRC-6-57" />

== See also ==
* [[Compartment fire]]
* [[Ventilation-controlled fire]]
* [[Smoke filling]]
* [[Vent flow]]
* [[Fire-induced pressure]]

== Notes ==
<references />

== References ==
* Australian Building Codes Board. (2020a). ''Data sheets: Fire Safety Verification Method'' (Handbook annex, Print version 1.1).
* Iqbal, N., & Salley, M. H. (2004). ''Fire Dynamics Tools (FDTs): Quantitative fire hazard analysis methods for the U.S. Nuclear Regulatory Commission Fire Protection Inspection Program'' (NUREG-1805). U.S. Nuclear Regulatory Commission, Office of Nuclear Reactor Regulation.
* McGrattan, K. B., Bouldin, C., & Forney, G. P. (2005). ''Computer simulation of the fires in the World Trade Center towers'' (NIST NCSTAR 1-5F). National Institute of Standards and Technology.
* Stroup, D., Lindeman, A., McGrattan, K., Peacock, R., Overholt, K., Joglar, F., LeStrange, S., & Montanez, S. (2016). ''Verification and validation of selected fire models for nuclear power plant applications: Supplement 1'' (NUREG-1824, Supplement 1; EPRI 3002002182). U.S. Nuclear Regulatory Commission, Office of Nuclear Regulatory Research; Electric Power Research Institute.
* Van Rossum, E. (2027). ''Introduction to Fire Safety Engineering''. Fire Safety Academy.
* Walton, W. D., & Thomas, P. H. (2002). Estimating temperatures in compartment fires. In P. J. DiNenno, D. Drysdale, C. L. Beyler, W. D. Walton, R. L. P. Custer, J. R. Hall Jr., & J. M. Watts Jr. (Eds.), ''SFPE handbook of fire protection engineering'' (3rd ed., pp. 3-171–3-188). National Fire Protection Association; Society of Fire Protection Engineers.

[[Category:Enclosure Fires]]