Mesh resolution - Revision history

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

2026-11-02T16:38:44Z

Add 4 watermarked SVG illustrations from the Expert visual plan

← Older revision		Revision as of 17:38, 2 November 2026
Line 22:		Line 22:

	The grid also limits the geometry. In the fire tests of the National Institute of Standards and Technology (NIST) for the World Trade Center investigation, a 10 cm grid described the test compartment to within about 5 cm.<ref name="M5F-28">McGrattan et al. (2005), p. 28.</ref> On a 50 cm grid, a 20 cm shift of the top of a window changed the depth of the hot upper layer by the same amount. The upper layer temperature then changed by roughly 100 °C.<ref name="M5F-28" /> Small objects, such as cable trays and floor trusses, may not fit on a coarse grid at all.<ref name="N7-5-2">McGrattan (2007), p. 5-2.</ref><ref name="M5F-27">McGrattan et al. (2005), p. 27.</ref>		The grid also limits the geometry. In the fire tests of the National Institute of Standards and Technology (NIST) for the World Trade Center investigation, a 10 cm grid described the test compartment to within about 5 cm.<ref name="M5F-28">McGrattan et al. (2005), p. 28.</ref> On a 50 cm grid, a 20 cm shift of the top of a window changed the depth of the hot upper layer by the same amount. The upper layer temperature then changed by roughly 100 °C.<ref name="M5F-28" /> Small objects, such as cable trays and floor trusses, may not fit on a coarse grid at all.<ref name="N7-5-2">McGrattan (2007), p. 5-2.</ref><ref name="M5F-27">McGrattan et al. (2005), p. 27.</ref>
			[[File:fse-mesh-resolution-v4.svg\|thumb\|upright=1.5\|alt=A wall with a window drawn on a fine and on a coarse grid; on the coarse grid the window edges jump to the nearest cell lines.\|A window in a wall on a 0.1 m grid (left) and a 0.5 m grid (right): the opening snaps to the nearest cell faces. Based on McGrattan et al. (2005), p. 28, and McGrattan et al. (2026), p. 6.]]

	== Discretisation error and computing cost ==		== Discretisation error and computing cost ==
Line 31:		Line 32:

	Halving the cell size therefore also halves the time step. Together, these effects increase the computing time by a factor of 16.<ref name="N7-5-1" /> For LES, keeping the CFL number close to 1 also keeps the time step consistent with the cell size. If a faster run is needed, it is better to use a coarser grid than to take much larger time steps on a fine grid.<ref name="FDS-47" />		Halving the cell size therefore also halves the time step. Together, these effects increase the computing time by a factor of 16.<ref name="N7-5-1" /> For LES, keeping the CFL number close to 1 also keeps the time step consistent with the cell size. If a faster run is needed, it is better to use a coarser grid than to take much larger time steps on a fine grid.<ref name="FDS-47" />
			[[File:fse-mesh-resolution-v2.svg\|thumb\|upright=1.5\|alt=Log chart: as cells shrink from 0.4 m to 0.025 m, relative error falls from 16 to 0.06 while relative computing time rises from 0.004 to 256.\|Relative discretisation error and computing time of a second-order model compared with a 0.1 m grid. Based on McGrattan (2007), p. 5-1.]]

	A validation study of FDS for nuclear power plants shows this trade-off in practice. A plume temperature in one test was calculated on grids of 10 cm, 7.5 cm and 5 cm. The 10 cm run took a few hours on a single 2.4 GHz processor, while the 5 cm run took a few days.<ref name="N7-5-1" />		A validation study of FDS for nuclear power plants shows this trade-off in practice. A plume temperature in one test was calculated on grids of 10 cm, 7.5 cm and 5 cm. The 10 cm run took a few hours on a single 2.4 GHz processor, while the 5 cm run took a few days.<ref name="N7-5-1" />
Line 80:		Line 82:

	The validation report for nuclear power plants states that past experience has shown a ratio of 5 to 10 to produce favourable results at a moderate computational cost.<ref name="N7-6-5" /> The NIST World Trade Center investigation used a ratio of 10 when it chose a 10 cm grid for its laboratory fire tests.<ref name="M5F-19">McGrattan et al. (2005), p. 19.</ref> A later study compared FDS flame heights with an empirical correlation for <math>Q^</math> values from 0.1 to 10,000. It used ratios of 5, 10 and 20. Except at the lowest values of <math>Q^</math>, the predictions fell within the uncertainty of the correlation itself.<ref name="NS1-6-26" /> The table lists the ratios used in the earlier validation study.<ref name="N7-6-4" />		The validation report for nuclear power plants states that past experience has shown a ratio of 5 to 10 to produce favourable results at a moderate computational cost.<ref name="N7-6-5" /> The NIST World Trade Center investigation used a ratio of 10 when it chose a 10 cm grid for its laboratory fire tests.<ref name="M5F-19">McGrattan et al. (2005), p. 19.</ref> A later study compared FDS flame heights with an empirical correlation for <math>Q^</math> values from 0.1 to 10,000. It used ratios of 5, 10 and 20. Except at the lowest values of <math>Q^</math>, the predictions fell within the uncertainty of the correlation itself.<ref name="NS1-6-26" /> The table lists the ratios used in the earlier validation study.<ref name="N7-6-4" />
			[[File:fse-mesh-resolution-v1.svg\|thumb\|upright=1.5\|alt=Three grids over the same fire: the fire spans about 3 cells of 0.2 m, 6 cells of 0.1 m and 12 cells of 0.05 m.\|The same 340 kW fire (characteristic diameter about 0.62 m) on grids of 0.2 m, 0.1 m and 0.05 m cells. Based on McGrattan (2007), p. 6-5, and Ministry of Business, Innovation and Employment (2013), p. 219.]]

	{\| class="wikitable"		{\| class="wikitable"
Line 98:		Line 101:

	The NBS Multi-Room tests show that the ratio alone does not decide the outcome. Only four cells spanned the burner, yet the predicted temperatures in all three rooms were within the experimental uncertainty. The ceiling was low, so the plume was less important for the average room conditions.<ref name="N7-6-5" />		The NBS Multi-Room tests show that the ratio alone does not decide the outcome. Only four cells spanned the burner, yet the predicted temperatures in all three rooms were within the experimental uncertainty. The ceiling was low, so the plume was less important for the average room conditions.<ref name="N7-6-5" />
			[[File:fse-mesh-resolution-v3.svg\|thumb\|upright=1.5\|alt=Chart of D* divided by cell size against heat release rate for cells of 0.05, 0.1, 0.2 and 0.4 m, with the band from 5 to 10 shaded.\|Ratio D*/δx for fires of 100 kW to 5000 kW and four cell sizes; the shaded band marks the range of 5 to 10. Based on McGrattan (2007), p. 6-5.]]

	=== Worked example ===		=== Worked example ===

Mark: FSE wiki: wording corrections after fact-check (source of the 5 to 10 rule of thumb; translation of section 11.4) (written by Claude Opus 5.5)

2026-11-01T22:33:01Z

FSE wiki: wording corrections after fact-check (source of the 5 to 10 rule of thumb; translation of section 11.4) (written by Claude Opus 5.5)

← Older revision		Revision as of 23:33, 1 November 2026
Line 79:		Line 79:
	The expressions describe buoyant fires. A value of <math>Q^*</math> above about 2.5 indicates a high-momentum jet fire rather than a buoyant fire.<ref name="CK-3-205">Cox and Kumar (2002), p. 3-205.</ref> The ratio is a guide for choosing a grid, not a measure of accuracy in itself.<ref name="N7-6-5" />		The expressions describe buoyant fires. A value of <math>Q^*</math> above about 2.5 indicates a high-momentum jet fire rather than a buoyant fire.<ref name="CK-3-205">Cox and Kumar (2002), p. 3-205.</ref> The ratio is a guide for choosing a grid, not a measure of accuracy in itself.<ref name="N7-6-5" />

	The validation ~~study~~ for nuclear power plants ~~found~~ that a ratio of 5 to 10 ~~usually produces~~ favourable results at a moderate computational cost.<ref name="N7-6-5" /> The NIST World Trade Center investigation used a ratio of 10 when it chose a 10 cm grid for its laboratory fire tests.<ref name="M5F-19">McGrattan et al. (2005), p. 19.</ref> A later study compared FDS flame heights with an empirical correlation for <math>Q^</math> values from 0.1 to 10,000. It used ratios of 5, 10 and 20. Except at the lowest values of <math>Q^</math>, the predictions fell within the uncertainty of the correlation itself.<ref name="NS1-6-26" /> The table lists the ratios used in the earlier validation study.<ref name="N7-6-4" />		The validation report for nuclear power plants states that past experience has shown a ratio of 5 to 10 to produce favourable results at a moderate computational cost.<ref name="N7-6-5" /> The NIST World Trade Center investigation used a ratio of 10 when it chose a 10 cm grid for its laboratory fire tests.<ref name="M5F-19">McGrattan et al. (2005), p. 19.</ref> A later study compared FDS flame heights with an empirical correlation for <math>Q^</math> values from 0.1 to 10,000. It used ratios of 5, 10 and 20. Except at the lowest values of <math>Q^</math>, the predictions fell within the uncertainty of the correlation itself.<ref name="NS1-6-26" /> The table lists the ratios used in the earlier validation study.<ref name="N7-6-4" />

	{\| class="wikitable"		{\| class="wikitable"
Line 137:		Line 137:
	A third misconception is that good agreement on a coarse grid proves that the grid is adequate. NIST obtained encouraging results when it simulated a laboratory test on its coarse World Trade Center grid. The agreement occurred because a stable, uniform hot layer formed, so the conservation of energy controlled the temperatures. NIST stressed that coarsening a grid to that extent does not, as a rule, give such good results.<ref name="M5F-27" />		A third misconception is that good agreement on a coarse grid proves that the grid is adequate. NIST obtained encouraging results when it simulated a laboratory test on its coarse World Trade Center grid. The agreement occurred because a stable, uniform hot layer formed, so the conservation of energy controlled the temperatures. NIST stressed that coarsening a grid to that extent does not, as a rule, give such good results.<ref name="M5F-27" />

	Even a careful grid study leaves some uncertainty. The NUREG-1824 validation reports treat the grid size of a CFD model as part of the intrinsic uncertainty of the model, which they did not try to quantify.<ref name="N2-1-8">Hamins and McGrattan (2007), p. 1-8.</ref> Coarse cells give no reliable information at scales smaller than the cells.<ref name="VR-11.4" /> It is also easy to produce impressive images of a flow, so the sensitivity of the solution to further mesh refinement has to be shown, not assumed.<ref name="CK-3-206" />		Even a careful grid study leaves some uncertainty. The NUREG-1824 validation reports treat the grid size of a CFD model as part of the intrinsic uncertainty of the model, which they did not try to quantify.<ref name="N2-1-8">Hamins and McGrattan (2007), p. 1-8.</ref> Coarse cells give little or no reliable information at scales smaller than the cells.<ref name="VR-11.4" /> It is also easy to produce impressive images of a flow, so the sensitivity of the solution to further mesh refinement has to be shown, not assumed.<ref name="CK-3-206" />

	== See also ==		== See also ==

Mark: FSE wiki: align footnotes and references with the canonical reference list (written by Claude Opus 5.5)

2026-11-01T22:10:00Z

FSE wiki: align footnotes and references with the canonical reference list (written by Claude Opus 5.5)

Show changes

Mark: New FSE wiki article: Mesh resolution (written by Claude Opus 5.5)

2026-11-01T21:58:11Z

New FSE wiki article: Mesh resolution (written by Claude Opus 5.5)

New page

'''Mesh resolution''', also called grid resolution, is the fineness of the computational grid in a [[computational fluid dynamics]] (CFD) model. The model divides the space into many small cells, which together form the grid, and solves the flow equations in every cell.<ref name="VR-11.1">Van Rossum (2027), §11.1.</ref><ref name="VR-11.3">Van Rossum (2027), §11.3.</ref> Mesh resolution is usually expressed as the cell size, written δx. In the [[Fire Dynamics Simulator]] (FDS), the cell size also sets the size of the smallest flow structures that the model calculates directly.<ref name="FDS-6">McGrattan et al. (2026), p. 6.</ref><ref name="FDS-47">McGrattan et al. (2026), p. 47.</ref>

Smaller cells show more detail, but they do not automatically give a better or more reliable result.<ref name="VR-11.3" /> They also cost much more computing time. Halving the cell size in FDS reduces the numerical error by a factor of four, but it increases the computing time by a factor of sixteen.<ref name="N7-5-1">U.S. Nuclear Regulatory Commission and Electric Power Research Institute (2007b), p. 5-1.</ref> Engineers therefore choose a cell size with a rule of thumb based on the size of the fire, and then test that choice with a grid sensitivity study.<ref name="N7-6-5">U.S. Nuclear Regulatory Commission and Electric Power Research Institute (2007b), p. 6-5.</ref> The report of a CFD study should explain why the chosen grid is fine enough, which sensitivity check was carried out and which uncertainty remains.<ref name="VR-11.3" />

== Cells, grids and filter width ==
A CFD model breaks the three-dimensional space of a fire scenario into cells, and the collection of cells is called the computational grid or mesh.<ref name="VR-11.3" /> General-purpose CFD software can often use complex cell shapes, such as tetrahedra, to follow complex geometry. Fire-specific software such as FDS uses only rectangular cells, so every object is reduced to rectangular blocks. Rectangular cells are simpler and faster, but they represent fire and flow effects around complex shapes less accurately.<ref name="VR-11.3" />

FDS prefers uniform meshes, in which all cells have the same size. The only numerical parameters that the user chooses are the three dimensions of the grid. Walls and other obstructions "snap" to the nearest cell faces, so the geometry is only as detailed as the grid allows.<ref name="FDS-6" /> Several meshes can be combined and run in parallel on different processors.<ref name="FDS-6" /> This allows a fine mesh near the fire and coarser meshes in areas where little happens, for example in a large car park.<ref name="VR-11.3" />

The cell size has a second meaning in FDS, which is a [[large eddy simulation]] (LES) model. LES computes the average values of mass, momentum and energy in each cell. The effect of the smaller, unresolved motions inside a cell is described by submodels. The width of the averaging "filter" is linked directly to the cell size.<ref name="FDS-6" /> For cells that are not perfect cubes, FDS uses the cube root of the cell volume:<ref name="FDS-27">McGrattan et al. (2026), p. 27.</ref>

<math display="block">\Delta = \left( \delta x \, \delta y \, \delta z \right)^{1/3}</math>

A coarser grid therefore means that more of the turbulence and combustion has to be modelled rather than calculated.<ref name="FDS-28" /> The cell shape matters as well. Good practice is to keep the cell proportions as close as possible to 1:1:1, so that cells are not strongly stretched in one direction.<ref name="VR-11.3" />

== Why cell size matters ==
CFD models solve an approximate form of the conservation equations of mass, momentum and energy. Replacing derivatives by differences between neighbouring cells introduces an error that depends on the size of the cells.<ref name="N7-5-1" /> In some numerical schemes, part of this error acts like extra diffusion, so sharp features in the flow are smoothed out. This is called numerical diffusion.<ref name="CK-3-202">Cox and Kumar (2002), p. 3-202.</ref>

A grid can only describe features that are larger than its cells. The smallest turbulent eddies in a fire are about one millimetre in size. Calculating them directly, which is called direct numerical simulation (DNS), is impractical except for special research flames.<ref name="FDS-27" /><ref name="FDS-28">McGrattan et al. (2026), p. 28.</ref> Resolving the chemical reaction zone of a fire would likewise require cells smaller than 1 mm.<ref name="CK-3-194">Cox and Kumar (2002), p. 3-194.</ref> For domains between metres and kilometres in size, the affordable cell size in most LES fire calculations ranges from centimetres to metres.<ref name="FDS-28" />

In large fire simulations, cells are typically tens of centimetres wide. The thin flame sheet then lies inside a single cell, and the computed cell-average temperature can be much lower than the temperature in the flame itself.<ref name="FDS-12">McGrattan et al. (2026), p. 12.</ref> Measured time-averaged gas temperatures in fires seldom exceed about 1300 °C. Because CFD averages temperatures over each cell, it is unlikely to predict such peak values.<ref name="CK-3-206" />

The grid also limits the geometry. In the fire tests of the National Institute of Standards and Technology (NIST) for the World Trade Center investigation, a 10 cm grid described the test compartment to within about 5 cm.<ref name="M5F-28">McGrattan et al. (2005), p. 28.</ref> On a 50 cm grid, a 20 cm shift of the top of a window changed the depth of the hot upper layer by the same amount. The upper layer temperature then changed by roughly 100 °C.<ref name="M5F-28" /> Small objects, such as cable trays and floor trusses, may not fit on a coarse grid at all.<ref name="N7-5-2">U.S. Nuclear Regulatory Commission and Electric Power Research Institute (2007b), p. 5-2.</ref><ref name="M5F-27">McGrattan et al. (2005), p. 27.</ref>

== Discretisation error and computing cost ==
FDS uses second-order accurate approximations in both time and space. The discretisation error is therefore proportional to the square of the cell size and of the time step.<ref name="N7-5-1" /> In theory, halving the cell size reduces this error by a factor of four. The price is high, because the number of cells grows by a factor of two in each of the three directions.<ref name="N7-5-1" />

The time step must shrink as well, because the model advances in small steps of time. The size of each step is limited by the cell size and by the flow speed in the cell.<ref name="FDS-47" /> The Courant-Friedrichs-Lewy (CFL) condition states that a fluid element should not travel through more than one cell in a single time step:<ref name="FDS-47" />

<math display="block">\mathrm{CFL} = \frac{\delta t \, \lVert \mathbf{u} \rVert}{\delta x} \approx 1</math>

Halving the cell size therefore also halves the time step. Together, these effects increase the computing time by a factor of 16.<ref name="N7-5-1" /> For LES, keeping the CFL number close to 1 also keeps the time step consistent with the cell size. If a faster run is needed, it is better to use a coarser grid than to take much larger time steps on a fine grid.<ref name="FDS-47" />

A validation study of FDS for nuclear power plants shows this trade-off in practice. A plume temperature in one test was calculated on grids of 10 cm, 7.5 cm and 5 cm. The 10 cm run took a few hours on a single 2.4 GHz processor, while the 5 cm run took a few days.<ref name="N7-5-1" />

The 5 cm prediction was noticeably better, because the entrainment of air into the plume was described more faithfully. Further refinement made no noticeable difference to the result.<ref name="N7-5-1" /> The point of diminishing returns depends on the application and on the quantity of interest. A hot gas layer temperature usually needs a less fine grid than the heat flux to a target close to the fire.<ref name="N7-5-1" />

The cost explains why very large problems use coarse cells. The NIST simulations of the fires in the World Trade Center towers used cells of 0.5 m by 0.5 m by 0.4 m, with about 125,000 cells per floor. One simulation of 105 minutes of fire in one tower took about a week on eight computers.<ref name="NIST-122">National Institute of Standards and Technology (2005), p. 122.</ref> By the factor of 16 described above, halving the cell size in all three directions would have raised this to about one million cells per floor and roughly sixteen weeks of computing time.<ref name="N7-5-1" />

== Characteristic fire diameter and the D*/δx ratio ==
Whether a cell is "small" depends on the size of the fire. The physical diameter of a fire is not always well defined: a burning room has no clear diameter, whereas a round pan of burning liquid has an obvious one.<ref name="N7-6-5" /> The [[dimensionless heat release rate]] <math>Q^*</math> compares the heat release rate of a fire with its physical diameter <math>D</math>:<ref name="N7-6-4">U.S. Nuclear Regulatory Commission and Electric Power Research Institute (2007b), p. 6-4.</ref>

<math display="block">Q^* = \frac{\dot{Q}}{\rho_\infty c_p T_\infty \sqrt{g D}\, D^2}</math>

A large value of <math>Q^*</math> describes a fire with a high energy output for its size, such as an oil well blowout. A low value describes a fire such as a brush fire, and most common accidental fires have values of about 1.<ref name="N7-6-4" /> Because the physical diameter is often unclear, CFD practice uses the characteristic fire diameter <math>D^*</math> instead:<ref name="N7-6-5" />

<math display="block">D^* = \left( \frac{\dot{Q}}{\rho_\infty c_p T_\infty \sqrt{g}} \right)^{2/5}</math>

When <math>Q^*</math> is about 1, <math>D^*</math> is comparable to the physical diameter of the fire. The ratio <math>D^*/\delta x</math> gives the number of cells that span this characteristic diameter.<ref name="NS1-6-26">U.S. Nuclear Regulatory Commission and Electric Power Research Institute (2016), pp. 6-26–6-27.</ref> The greater the ratio, the more of the fire dynamics are resolved directly and the more accurate the simulation becomes.<ref name="N7-6-5" />

{| class="wikitable"
! Symbol !! Description !! Unit
|-
| <math>\Delta</math> || LES filter width || m
|-
| <math>\delta x, \delta y, \delta z</math> || cell size in the three coordinate directions || m
|-
| <math>\delta t</math> || time step || s
|-
| <math>\lVert \mathbf{u} \rVert</math> || local flow speed || m/s
|-
| <math>\dot{Q}</math> || heat release rate of the fire || kW
|-
| <math>Q^*</math> || dimensionless heat release rate || –
|-
| <math>D</math> || physical diameter of the fire || m
|-
| <math>D^*</math> || characteristic fire diameter || m
|-
| <math>\rho_\infty</math> || density of the ambient air || kg/m³
|-
| <math>c_p</math> || specific heat of air || kJ/(kg·K)
|-
| <math>T_\infty</math> || ambient temperature || K
|-
| <math>g</math> || acceleration due to gravity || m/s²
|}

The expressions describe buoyant fires. A value of <math>Q^*</math> above about 2.5 indicates a high-momentum jet fire rather than a buoyant fire.<ref name="CK-3-205">Cox and Kumar (2002), p. 3-205.</ref> The ratio is a guide for choosing a grid, not a measure of accuracy in itself.<ref name="N7-6-5" />

The validation study for nuclear power plants found that a ratio of 5 to 10 usually produces favourable results at a moderate computational cost.<ref name="N7-6-5" /> The NIST World Trade Center investigation used a ratio of 10 when it chose a 10 cm grid for its laboratory fire tests.<ref name="M5F-19">McGrattan et al. (2005), p. 19.</ref> A later study compared FDS flame heights with an empirical correlation for <math>Q^*</math> values from 0.1 to 10,000. It used ratios of 5, 10 and 20. Except at the lowest values of <math>Q^*</math>, the predictions fell within the uncertainty of the correlation itself.<ref name="NS1-6-26" /> The table lists the ratios used in the earlier validation study.<ref name="N7-6-4" />

{| class="wikitable"
! Test series !! Heat release rate (kW) !! <math>D^*</math> (m) !! <math>\delta x</math> (m) !! <math>D^*/\delta x</math>
|-
| ICFMP BE #2 || 1800–3600 || 1.2–1.6 || 0.13 || 9–12
|-
| ICFMP BE #3 || 400–2300 || 0.7–1.3 || 0.10 || 7–13
|-
| ICFMP BE #4 || 3500 || 1.6 || 0.10 || 16
|-
| ICFMP BE #5 || 400 || 0.7 || 0.10 || 7
|-
| FM/SNL || 500 || 0.7 || 0.05 || 14
|-
| NBS Multi-Room || 100 || 0.4 || 0.10 || 4
|}

The NBS Multi-Room tests show that the ratio alone does not decide the outcome. Only four cells spanned the burner, yet the predicted temperatures in all three rooms were within the experimental uncertainty. The ceiling was low, so the plume was less important for the average room conditions.<ref name="N7-6-5" />

=== Worked example ===
A worked example in the New Zealand commentary on Verification Method C/VM2 (2013) models a hotel corridor with FDS. The fire grows to 340 kW, when sprinklers activate, and the guest room, corridor and lift lobby use cells of 0.1 m. The stairway, further from the fire, uses cells of 0.2 m.<ref name="CVM2C-219">Ministry of Business, Innovation and Employment (2013), p. 219.</ref> The ambient values are taken as <math>T_\infty</math> = 293 K, <math>\rho_\infty</math> = 1.2 kg/m³, <math>c_p</math> = 1.0 kJ/(kg·K) and <math>g</math> = 9.81 m/s².<ref name="VR-3.2">Van Rossum (2027), §3.2.</ref><ref name="VR-6.7">Van Rossum (2027), §6.7.</ref>

The denominator of the expression for <math>D^*</math> is 1.2 × 1.0 × 293 × √9.81, which is about 1101 kW/m<sup>5/2</sup>.<ref name="N7-6-5" /> The characteristic diameter is then (340 / 1101)<sup>2/5</sup>, or about 0.62 m. With 0.1 m cells, the ratio <math>D^*/\delta x</math> is about 6.2, inside the range of 5 to 10.

A ratio between 5 and 10 would need cells between about 0.062 m and 0.12 m. Halving the cells to 0.05 m would raise the ratio to about 12, at roughly 16 times the computing cost. The stairway cells of 0.2 m give a ratio of only about 3.1. The commentary uses these coarser cells because the stairway is further from the fire.<ref name="CVM2C-219" />

The same ambient values reproduce the table above. A 500 kW fire gives <math>D^*</math> = 0.73 m, which the validation study lists as 0.7 m. A 3500 kW fire gives 1.59 m, listed as 1.6 m.<ref name="N7-6-4" />

== Grid sensitivity studies ==
Determining which cell size to use for a given calculation is known as a grid sensitivity study.<ref name="N7-5-1" /> The idea is to run the same scenario on successively finer grids and to compare the results.<ref name="CK-3-206">Cox and Kumar (2002), p. 3-206.</ref> The aim is the coarsest grid that still gives reliable results.<ref name="VR-11.3" /> In the validation study for nuclear power plants, each grid was refined until the results no longer changed from run to run.<ref name="N7-6-5" />

A chapter on CFD in the ''SFPE Handbook of Fire Protection Engineering'' gives a checklist for such studies. The chapter deals mainly with Reynolds-averaged models rather than with LES.<ref name="CK-3-194" /> The grid study should preferably use at least three grid resolutions, ideally by doubling the number of cells in all directions. The vertical spacing needs special attention in buoyant flows.<ref name="CK-3-206" />

The mesh should be fine at the fire source and at openings. The aspect ratio of the cells should typically not exceed 50.<ref name="CK-3-206" /> The time step must suit the grid, and global mass and heat balances should be better than 95 %.<ref name="CK-3-206" />

A converged solution is not the same as a grid-independent solution. An iterative solver can reach very small residuals on a grid that is still too coarse. The solution must also be shown to be independent of the grid.<ref name="CK-3-202" /><ref name="CK-3-203">Cox and Kumar (2002), p. 3-203.</ref> In the terms of [[verification and validation]], a grid study is a verification task. It checks whether the model equations are solved with adequate accuracy, not whether the physics in the model is correct.<ref name="CK-3-203" />

The result of a grid study depends on the quantity examined. In one validation test, a 20 cm grid still predicted the average hot gas layer temperature reasonably well.<ref name="N7-5-2" /> However, the cable trays were too small to include, and the door and vent were described less accurately. The compartment pressure therefore had a greater error.<ref name="N7-5-2" /> Coarse grids can give reasonable values for quantities that follow directly from the conservation of mass and energy, such as average temperatures and pressures. Other results become less reliable or cannot be obtained at all.<ref name="N7-5-2" />

NIST tested its World Trade Center grid with office workstation fires on cells of 25 cm and 50 cm.<ref name="M5F-37">McGrattan et al. (2005), p. 37.</ref> Both grids gave similar results, because the burning rate depended mainly on heat from the hot upper layer.<ref name="M5F-39">McGrattan et al. (2005), p. 39.</ref> The report warns that this is not guaranteed for other applications. A fire that is small compared with its compartment needs a much finer grid, because coarse grids smear out steep gradients.<ref name="M5F-39" />

== Mesh resolution in design practice ==
Design models often combine meshes of different sizes, with a fine mesh near the fire and coarser meshes elsewhere.<ref name="VR-11.3" /> In the validation study, the two test series with a small fire in a large space used a fine grid in the plume, where the entrainment of fresh air matters most. Coarse grids were used elsewhere.<ref name="N7-6-5" />

Building regulations sometimes have to allow for the grid. In New Zealand, C/VM2 (second edition, 2025) lets the leakage through walls be modelled in CFD as two vents, one at floor level and one at ceiling level. If the required leakage area is smaller than one cell, it may be increased to fit the grid. The combined modelled area must not exceed 5 times the calculated leakage area.<ref name="CVM2-16">Ministry of Business, Innovation and Employment (2025), p. 16.</ref> The 2013 commentary example does this for a 10 mm gap around a door. It models the gap as two vents of 0.1 m by 0.1 m, which fit the 0.1 m mesh.<ref name="CVM2C-219" />

The same edition of C/VM2 asks for the design fire in a CFD model to be a burning area sized from a heat release rate per square metre. For most buildings, it allows 500 to 1000 kW/m², a range that accommodates different heat release rates and mesh sizes. A zone model uses a single value of 250 kW/m² instead.<ref name="CVM2-20">Ministry of Business, Innovation and Employment (2025), p. 20.</ref>

The quantity of interest also guides the grid near particular features. Predicting the activation of a heat detector or sprinkler requires a realistic ceiling jet, and therefore a reasonably fine mesh just below the ceiling.<ref name="CK-3-207">Cox and Kumar (2002), p. 3-207.</ref> Openings need enough cells as well. In one NIST test, a window was spanned by seven cells of 10 cm. The model predicted the inflow of air well, but it could not resolve the details of the outflow velocity, which varied by at least 3 m/s over the window height.<ref name="M5F-26">McGrattan et al. (2005), p. 26.</ref>

== Limits and common misconceptions ==
A frequent misconception is that a finer grid always gives a better answer. Finer cells give a more specific result, but this does not by itself make the result more reliable.<ref name="VR-11.3" /> The combustion, turbulence and water spray models in CFD remain approximations, and the modeller still has to prescribe the fire scenario and its heat release rate.<ref name="VR-11.4">Van Rossum (2027), §11.4.</ref>

A second misconception is that the rule of thumb for <math>D^*/\delta x</math> replaces a grid study. The validation study states explicitly that it does not.<ref name="N7-6-5" /> The ratio describes the fire plume, while other features, such as openings, targets or a ceiling jet, may need their own resolution.<ref name="N7-5-2" /><ref name="CK-3-207" />

A third misconception is that good agreement on a coarse grid proves that the grid is adequate. NIST obtained encouraging results when it simulated a laboratory test on its coarse World Trade Center grid. The agreement occurred because a stable, uniform hot layer formed, so the conservation of energy controlled the temperatures. NIST stressed that coarsening a grid to that extent does not, as a rule, give such good results.<ref name="M5F-27" />

Even a careful grid study leaves some uncertainty. The NUREG-1824 validation reports treat the grid size of a CFD model as part of the intrinsic uncertainty of the model, which they did not try to quantify.<ref name="N2-1-8">U.S. Nuclear Regulatory Commission and Electric Power Research Institute (2007a), p. 1-8.</ref> Coarse cells give no reliable information at scales smaller than the cells.<ref name="VR-11.4" /> It is also easy to produce impressive images of a flow, so the sensitivity of the solution to further mesh refinement has to be shown, not assumed.<ref name="CK-3-206" />

== See also ==
* [[Fire modelling]]
* [[Zone model]]
* [[Sensitivity analysis]]
* [[Finite element analysis]]
* [[Hand calculations]]

== Notes ==
<references />

== References ==
* Cox, G., & Kumar, S. (2002). Modeling enclosure fires using CFD. 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-194–3-218). National Fire Protection Association.
* 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.
* McGrattan, K., Hostikka, S., Floyd, J., McDermott, R., Vanella, M., Mueller, E., & Paul, C. (2026). ''Fire Dynamics Simulator technical reference guide: Volume 1. Mathematical model'' (6th ed., NIST Special Publication 1018-1). National Institute of Standards and Technology. https://doi.org/10.6028/NIST.SP.1018
* Ministry of Business, Innovation and Employment. (2013). ''Commentary for Building Code clauses C1–C6 and Verification Method C/VM2''.
* Ministry of Business, Innovation and Employment. (2025). ''C Protection from fire: Verification Method C/VM2. Framework for fire safety design'' (2nd ed.).
* National Institute of Standards and Technology. (2005). ''Final report on the collapse of the World Trade Center towers'' (NIST NCSTAR 1).
* U.S. Nuclear Regulatory Commission, & Electric Power Research Institute. (2007a). ''Verification and validation of selected fire models for nuclear power plant applications: Volume 2. Experimental uncertainty'' (NUREG-1824, EPRI 1011999).
* U.S. Nuclear Regulatory Commission, & Electric Power Research Institute. (2007b). ''Verification and validation of selected fire models for nuclear power plant applications: Volume 7. Fire Dynamics Simulator (FDS)'' (NUREG-1824, EPRI 1011999).
* U.S. Nuclear Regulatory Commission, & Electric Power Research Institute. (2016). ''Verification and validation of selected fire models for nuclear power plant applications: Supplement 1'' (NUREG-1824 Supplement 1, EPRI 3002002182).
* Van Rossum, E. (2027). ''Introduction to Fire Safety Engineering''. Fire Safety Academy.

[[Category:Fire Modelling]]