mathematical model for fire spread - university of delhimaths.du.ac.in/events/iwm/talks/nita.pdf ·...
TRANSCRIPT
![Page 1: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/1.jpg)
Mathematical Model
for Fire Spread
Prof. Nita H. Shah
Department of Mathematics
Gujarat University
Ahmedabad
April 2, 2015 IWM2015 Delhi Uni. 1
![Page 2: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/2.jpg)
What is Fire?
• Heat is supplied from the fire to the
potential fuel, the surface is dehydrated
and further heating raises the surface
temperature until the fuel begins to
pyrolyze and release combustible
gases.
April 2, 2015 IWM2015 Delhi Uni. 2
![Page 3: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/3.jpg)
When the gas evolution rate from the
potential fuel is sufficient to support
combustion, the gas is ignited by the
flame and the fire advances to a new
position.
April 2, 2015 IWM2015 Delhi Uni. 3
![Page 4: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/4.jpg)
• Gradually, a constant rate of spread is
attained which is called “quasi-steady
state” wherein the fire advances at a
rate that is the average of all the
elemental rates.
April 2, 2015 IWM2015 Delhi Uni. 4
![Page 5: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/5.jpg)
A mathematical model to determine rate
of fire spread and its intensity is
formulated. The model is developed by
R.C. Rothermal (1972) and it is
considered as the basis in the National
Fire Danger Rating System.
April 2, 2015 IWM2015 Delhi Uni. 5
Objective
![Page 6: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/6.jpg)
The model considers physical and
chemical properties of the fuel and the
environmental conditions in which it is
expected to tingle.
April 2, 2015 IWM2015 Delhi Uni. 6
![Page 7: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/7.jpg)
• The physical properties incorporated are
1. fuel loading
2. fuel depth
3. fuel particle surface area to volume ratio
4. fuel particle moisture, and
5. the moisture content at which extinction is
expected.
April 2, 2015 IWM2015 Delhi Uni. 7
![Page 8: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/8.jpg)
• Environmental inputs are
1. mean wind velocity
2. slope of terrain.
April 2, 2015 IWM2015 Delhi Uni. 8
![Page 9: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/9.jpg)
Mathematical Model
April 2, 2015 IWM2015 Delhi Uni. 9
quasi-steady rate of spread, ft/min
Horizontal heat flux absorbed by a unit volume of the fuel
at the time of ignition B.t.u./ft2 min
Effective bulk density (the amount of fuel per unit volume of
the fuel bed raised to ignition ahead of the advancing fire),
lb/ft3
Heat of pre-ignition (the heat required to bring a unit weight
of fuel to ignition), B.t.u./lb
The gradient of the vertical intensity evaluated at a plane at
a constant depth, 𝑧𝑐, of the fuel bed, B.t.u./ft3 min
Horizontal coordinate
𝑧 Vertical coordinate
R
XigI
be
igQ
c
z
z
I
z
x
![Page 10: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/10.jpg)
• Using conservation of energy principle
to a unit volume of fuel ahead of an
advancing fire in a homogeneous fuel
bed (by Fransen (1971))
April 2, 2015 IWM2015 Delhi Uni. 10
0heat flux received
from the source(1)
heat required for
ignition by the potential fuel
c
zXig
z
be ig
II dx
zR
Q
![Page 11: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/11.jpg)
• The fuel-reaction zone interface is fixed
and the unit volume is moving at a
constant depth, zc, from x = - infinity
towards the interface at x = 0. The unit
volume ignites at the interface.
April 2, 2015 IWM2015 Delhi Uni. 11
Analysis
![Page 12: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/12.jpg)
• The heat required for ignition depends upon
(a)Ignition temperature
(b) moisture content of the fuel, and
(c)Amount of fuel involved in the ignition
process.
April 2, 2015 IWM2015 Delhi Uni. 12
![Page 13: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/13.jpg)
• The energy per unit mass required for
ignition is the heat of pre-ignition, .
April 2, 2015 IWM2015 Delhi Uni. 13
igQ
( , ), . . . / (2)
where
: ratio of fuel moisture to ovendry weight
: ignition temp.
ig f ig
f
ig
Q f M T B t u lb
M
T
![Page 14: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/14.jpg)
• The amount of fuel involved in the
ignition process is the effective bulk
density, .
• An effective heating number is ratio of
the effective bulk density to the actual
bulk density.
April 2, 2015 IWM2015 Delhi Uni. 14
be
(3)
be
b
![Page 15: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/15.jpg)
• The effective heating number is
dimensionless which will be nearly
unity for fine fuels and decrease to zero
as fuel size increases. Therefore,
April 2, 2015 IWM2015 Delhi Uni. 15
bulk density, fuel size (4) be f
![Page 16: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/16.jpg)
• Propagating Flux
It is given by
April 2, 2015 IWM2015 Delhi Uni. 16
0
gradient of the vertical flux integrated from
-infinity to the fire fron
horizontal
flu
5
x
t
c
z
p Xig
z
II I dx ( )
z
![Page 17: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/17.jpg)
April 2, 2015 IWM2015 Delhi Uni. 17
![Page 18: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/18.jpg)
• The figures indicate that the vertical
flux is more significant during wind-
driven and upslope fires because the
flame tilts over the potential fuel,
thereby, increasing radiation, causing
direct flame contact and convective
heat transfer to the potential fuel.
April 2, 2015 IWM2015 Delhi Uni. 18
![Page 19: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/19.jpg)
• Assume that the vertical flux is small
for no-wind fires and let . This is
the basic heat flux component to which
all additional effects of wind and slope
are related.
• With (3) – (5) in (1) and and R=R0
for the no-wind case, we get
April 2, 2015 IWM2015 Delhi Uni. 19
0
p pI I
0
p pI I
00
(6) p b igI R Q
![Page 20: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/20.jpg)
• The energy release rate of the fire front
is produced by burning gases released
from the organic matters in the fuels.
April 2, 2015 IWM2015 Delhi Uni. 20
Reaction Intensity
![Page 21: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/21.jpg)
• The rate of change of this organic
matter from a solid to a gas is a good
approximation of the subsequent heat
release rate of the fire.
April 2, 2015 IWM2015 Delhi Uni. 21
![Page 22: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/22.jpg)
• The heat release rate per unit area of
the front is called the reaction intensity
and is defined as
April 2, 2015 IWM2015 Delhi Uni. 22
Heat content of fmass loss
rate per unit area
(7
in
)
the fire front
l ue
R
dwI h
dt
![Page 23: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/23.jpg)
• The reaction intensity is a function of
fuel parameters such as particle size,
bulk density, moisture, and chemical
composition.
• The reaction intensity is the source of
the no-wind propagating flux
April 2, 2015 IWM2015 Delhi Uni. 23
0
(8)p RI f I
![Page 24: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/24.jpg)
• Wind and slope change the
propagating heat flux by exposing the
potential fuel to additional convective
and radiant heat.
April 2, 2015 IWM2015 Delhi Uni. 24
Effect of Wind & Slope
![Page 25: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/25.jpg)
April 2, 2015 IWM2015 Delhi Uni. 25
![Page 26: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/26.jpg)
• Let and 𝜙𝑠 represent the additional
propagating flux produced by wind and
slope. The total propagating flux is
• Approximate rate of spread (eq.(1)) becomes
April 2, 2015 IWM2015 Delhi Uni. 26
w
0
1 (9) p p w sI I
0
1(10)
p w s
b ig
IR
Q
![Page 27: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/27.jpg)
Specific heat of dry wood
Temperature range to ignition
Fuel moisture lb. water / lb. dry wood
Specific heat of water
Temperature range to boiling
Latent heat of vaporization
April 2, 2015 IWM2015 Delhi Uni. 27
Heat of Pre-ignition
pdC
igT
fM
pwC
bT
V
![Page 28: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/28.jpg)
• for cellulosic fuels is determined by
considering the change in specific heat
from ambient to ignition temperature
and the latent heat of vaporization of
the moisture as
April 2, 2015 IWM2015 Delhi Uni. 28
igQ
(11)ig pd ig f pw bQ C T M C T V
![Page 29: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/29.jpg)
• Taking temperature to ignition in the
range of 200c to 3200c and boiling
temperature to be @ 1000c, eq. (11)
becomes
April 2, 2015 IWM2015 Delhi Uni. 29
250 1116 (12) ig fQ M ,B.t.u. / lb
![Page 30: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/30.jpg)
• To determine the effective bulk density,
we need to compute the efficiency of
heating as a function of particle size.
April 2, 2015 IWM2015 Delhi Uni. 30
Effective Bulk Density
![Page 31: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/31.jpg)
• An exponential fit is given by
April 2, 2015 IWM2015 Delhi Uni. 31
1
138(13)
where is particle surface area to volume ratio, ft
exp
![Page 32: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/32.jpg)
• Rearrange eq.(7) as
• To solve (16), integrate over the
reaction zone depth D, and over the
limits of loading in the reaction zone.
April 2, 2015 IWM2015 Delhi Uni. 32
(15)
where the quassy-steady rate of spread
(16)
R
R
dw dxI h
dx dt
dxR,
dt
I dx Rhdw
x
w
![Page 33: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/33.jpg)
where
April 2, 2015 IWM2015 Delhi Uni. 33
0
( ) (18) r
n
wD
R R n r
w
I dx Rh dw I D Rh w w
D Reaction zone depth (front to rear)
Net initial fuel loading, lb/ft2
𝑤𝑟 Residue loading immediately after passage of the
reaction zone, lb/ft2
nw
![Page 34: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/34.jpg)
• The time taken for the fire front to travel
a distance equivalent to the depth of
one reaction zone is the reaction time,
Then eq. (18) becomes
(20)
April 2, 2015 IWM2015 Delhi Uni. 34
(19)R
D
R
n r
R
R
h w w I =
![Page 35: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/35.jpg)
• Next, we define a maximum reaction
intensity where there is no loading
residue left after the reaction zone is
passed and where the reaction time
remains unchanged.
April 2, 2015 IWM2015 Delhi Uni. 35
max (21) nR
R
hwI
![Page 36: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/36.jpg)
• The reaction zone efficiency is
In (23), the net fuel loading is given by
April 2, 2015 IWM2015 Delhi Uni. 36
max
(22)
(23)
n rR
R n
nR
R
w wI
I w
w hI
0
2
0
(24)1
where
ovendry fuel loading, lb/ft
fuel mineral content, lb minerals/lb dry fuel
n
T
T
ww
S
w
S
![Page 37: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/37.jpg)
• The reaction velocity denotes the
completeness and rate of fuel
consumption. It is defined as the ratio
of the reaction zone efficiency to the
reaction time.
April 2, 2015 IWM2015 Delhi Uni. 37
Reaction Velocity
*reaction velocity/min. (25)R
![Page 38: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/38.jpg)
1. Moisture content
2. Mineral content
3. Particle size, and
4. Fuel bed bulk density
April 2, 2015 IWM2015 Delhi Uni. 38
Effective Fuel Parameters for
Reaction Velocity
![Page 39: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/39.jpg)
Then
One can evaluate reaction velocity and the moisture
and mineral damping coeffiecients by experiments.
April 2, 2015 IWM2015 Delhi Uni. 39
Potential reaction velocity/min
Moisture damping coefficient having values
ranging from 1 to 0
Mineral damping coefficient having value
ranging from 1 to 0
M
s
(26)
(27)
then the reaction intensity is
M s
R n M sI w h
![Page 40: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/40.jpg)
• It is defined as
April 2, 2015 IWM2015 Delhi Uni. 40
Moisture Damping coefficient
max 0
(28)
f
RM
R M
I
I
![Page 41: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/41.jpg)
In adjacent fig., is fuel
moisture and 𝑴𝒙 is the
moisture content of the
fuel at which the fire will
not spread.
Moisture damping
coefficient accounts for
the decrease in intensity
caused by the
combination of fuel that
initially contained
moisture.
April 2, 2015 IWM2015 Delhi Uni. 41
fM
![Page 42: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/42.jpg)
• It was evaluated from thermo gravimetric
analysis (TGA) data of natural fuels by Philot
(1968). It is assumed that the ration of the
normalized decomposition rate will be same
as the normalized reaction intensity.
April 2, 2015 IWM2015 Delhi Uni. 42
Mineral Damping coefficient
![Page 43: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/43.jpg)
where 𝑺𝒆 is effective
mineral content.
April 2, 2015 IWM2015 Delhi Uni. 43
0.190.174 (30)s eS
![Page 44: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/44.jpg)
• Next, we need to consider two
parameters:
1. The reaction intensity – fuel bed
compactness, and
2. Fuel particle size
April 2, 2015 IWM2015 Delhi Uni. 44
Physical Fuel Parameters
![Page 45: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/45.jpg)
Assume that low values of fire intensity
and rate of spread occur at the two
extremes of compactness (loose and
dense).
In dense beds, this can be attributed to
low air-to-fuel ratio and to poor
penetration of the heat beyond the upper
surface of the fuel array. April 2, 2015 IWM2015 Delhi Uni. 45
![Page 46: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/46.jpg)
• In loose beds, low intensity and poor
spread are attributed to heat transfer
looses between particles and to lack of
fuels.
April 2, 2015 IWM2015 Delhi Uni. 46
![Page 47: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/47.jpg)
• Between these two extremes, there
must be an optimum best equilibrium
of air, fuel and heat transfer for both
maximum fire intensity and reaction
velocity.
April 2, 2015 IWM2015 Delhi Uni. 47
![Page 48: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/48.jpg)
• The compactness of the fuel bed is
quantified by the packing ration, which
is the ratio of the fuel array to fuel
particle density.
April 2, 2015 IWM2015 Delhi Uni. 48
3
3
packing ratio, (31)
where : fuel array bulk density, lb/ft
: fuel particle density, lb/ft
b
p
b
p
![Page 49: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/49.jpg)
• The surface area-to-volume ratio, is
used to quantify the particle size.
• For fuels that are long w.r.t. the
thickness
April 2, 2015 IWM2015 Delhi Uni. 49
4(32)
here denotes diameter of circular particles
or edge length of square particles.
d
d
![Page 50: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/50.jpg)
The reaction time,
is defined on the
derivative curve as the
time from initial mass
loss until the loss
stabilizes at a steady
rate.
During the reaction
time mass loss rate is
linear. Also the
duration of the
constant mass rate is
dependent on the
length of platform. April 2, 2015 IWM2015 Delhi Uni. 50
The reaction time
could be thought of as
the fire burned off.
R
![Page 51: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/51.jpg)
The mass loss rate is
The efficiency of fire is
April 2, 2015 IWM2015 Delhi Uni. 51
(33)
W : width of the platform.
n r
dmw w RW
dt
1(34)
n
dm
w RW dt
![Page 52: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/52.jpg)
• With efficiency and the reaction time,
the reaction velocity is
April 2, 2015 IWM2015 Delhi Uni. 52
1(35)
n R
dm
w RW dt
![Page 53: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/53.jpg)
The potential reaction
velocity to disassociate
the reaction velocity
from the effects of the
moisture and minerals
of the fuel is .
The reaction velocity must
drop to zero if there is no
fuel to support combustion.
April 2, 2015 IWM2015 Delhi Uni. 53
M s
Tightly packed fine fuels
have lower reaction
velocity than do longer
fuels at the same packing
ratio.
![Page 54: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/54.jpg)
The loss of reaction
intensity of the fine
fuel can be seen as
the difference in
flaming vigor.
April 2, 2015 IWM2015 Delhi Uni. 54
![Page 55: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/55.jpg)
To determine the
general equation as a
function of and 𝜎 is
given by
April 2, 2015 IWM2015 Delhi Uni. 55
1.5
1.5
0.8189
(36)495 0.0594
3.348 (37)op
![Page 56: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/56.jpg)
Combining these two eqs. gives
April 2, 2015 IWM2015 Delhi Uni. 56
max
0.1
exp 1 (38)
1with (39)
4.77 7.27
A
op op
A
A
![Page 57: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/57.jpg)
These will predict reaction velocity for
any combination of fuel particle size and
any packing ratio. The eqs. will predict
reasonable values when input
parameters are extrapolated. This will
help us to predict reaction intensity and,
subsequently, rate of spread over a wide
range of fuel and environmental
combinations. April 2, 2015 IWM2015 Delhi Uni. 57
![Page 58: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/58.jpg)
The no-wind propagating flux is
A ratio 𝝃 relating the propagating flux to
the reaction intensity is
April 2, 2015 IWM2015 Delhi Uni. 58
Propagating Flux
00( we know from (6))p b igI R Q
(41)
p o
R
I
I
![Page 59: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/59.jpg)
April 2, 2015 IWM2015 Delhi Uni. 59
![Page 60: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/60.jpg)
A respective fit is
April 2, 2015 IWM2015 Delhi Uni. 60
1exp (0.792 0.681 )( 0.1) (42)
192 0.259
![Page 61: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/61.jpg)
It is seen that
increases with
increase in , but at a
decreasing rate.
will attain a
maximum and then
decreases.
This is reasonable,
considering the fact
that the fuel array is
becoming so compact
that the intensity has
decreased.
April 2, 2015 IWM2015 Delhi Uni. 61
0pI
0pI
![Page 62: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/62.jpg)
Combining the heat source and heat sink
terms results into the no-wind rate of
spread equation as
April 2, 2015 IWM2015 Delhi Uni. 62
Rate of Spread
0 (43)R
b ig
IR
Q
![Page 63: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/63.jpg)
April 2, 2015 IWM2015 Delhi Uni. 63
Fire Flames
![Page 64: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/64.jpg)
Reference
1. Richard C. Rothermal (1972):
Intermountain Forest and Range
experiment station Forest service, Res.
Paper 115.
April 2, 2015 IWM2015 Delhi Uni. 64
![Page 65: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/65.jpg)
April 2, 2015 IWM2015 Delhi Uni. 65
![Page 66: Mathematical Model for Fire Spread - University of Delhimaths.du.ac.in/Events/IWM/talks/nita.pdf · Mathematical Model for Fire Spread Prof. Nita H. Shah Department of Mathematics](https://reader035.vdocument.in/reader035/viewer/2022070720/5ee0a82ead6a402d666bcf2b/html5/thumbnails/66.jpg)
April 2, 2015 IWM2015 Delhi Uni. 66