Skip to content

Commit d4ff074

Browse files
committed
Merge remote-tracking branch 'firemodels/master' into FireX
2 parents bb3cb0b + a5a4fc7 commit d4ff074

File tree

52 files changed

+11095
-4806
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

52 files changed

+11095
-4806
lines changed

Manuals/Bibliography/FDS_general.bib

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -6602,6 +6602,14 @@ @TECHREPORT{Vettori:2
66026602
month = {December},
66036603
}
66046604

6605+
@INPROCEEDINGS{Vigne:2009,
6606+
author = {G. Vigne and J. J\"{o}nsson},
6607+
year = {2009},
6608+
month = {October},
6609+
title = {{Experimental Research---Large-Scale Tunnel Fire Tests and the use of CFD modelling to predict Heat Flux and Thermal Behavior}},
6610+
booktitle = {Proceedings of the International Congress on Fire Protection and Life Safety in Buildings and Transportation Systems, Santander, Spain}
6611+
}
6612+
66056613
@ARTICLE{Vreman:2004,
66066614
author = {Vreman, B.},
66076615
title = {{An eddy-viscosity subgrid-scale model for turbulent shear flow: Algebraic theory and applications}},

Manuals/FDS_Validation_Guide/Experiment_Chapter.tex

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -10,27 +10,27 @@ \chapter{Description of Experiments}
1010
\section{Aalto Woods Experiments}
1111
\label{Aalto_Woods_Description}
1212

13-
Rinta-Paavola et al. developed a material model for Norway spruce and Scots pine woods to predict char front progress and heat release in burning timber~\cite{Rinta-Paavola:2023}. The material models are calibrated using temperature, mass loss rate and heat release rate measurements of wood specimens in a cone calorimeter. The model is intended for a number of different fire scenarios, such as fully developed fire with only little oxygen in the compartment, or the decay period after flaming with exposed char still smoldering. Cone calorimeter measurements were conducted under flaming, inert, and and smoldering conditions to replicate these scenarios. The estimated material models are able to predict the experiments in the scale of cone calorimeter.
13+
Rinta-Paavola et al.~developed a material model for Norway spruce and Scots pine woods to predict char front progress and heat release in burning timber~\cite{Rinta-Paavola:2023}. The material models are calibrated using temperature, mass loss rate and heat release rate measurements of wood specimens in a cone calorimeter. The model is intended for fully-developed compartment fires or the subsequent decay period with exposed, smoldering char. Cone calorimeter measurements were conducted under flaming, inert, and smoldering conditions to replicate these scenarios. The material models are able to predict the bench-scale experiments.
1414

15-
The material model for spruce developed using cone calorimeter measurements is tested in prediction of heat release rates in two previously existing experiments: an ISO 9705 room corner test~\cite{Sundstrom:1998} (test code M12 in the original source) and a modified room corner test originally reported in Finnish in~\cite{Hietaniemi:2001}. English description of the modified test is available in~\cite{Hietaniemi:1}.
15+
The material model for spruce developed using cone calorimeter measurements is evaluated for its prediction of heat release rates in two previous experiments: an ISO~9705 room corner test~\cite{Sundstrom:1998} (test code M12 in the original source) and a modified room corner test originally reported in Finnish in~\cite{Hietaniemi:2001}. An English description of the modified test is available in Ref.~\cite{Hietaniemi:1}.
1616

17-
In a room corner test, the walls and the ceiling of a 2.4~m by 3.6~m by 2.4~m (width, length, height) test compartment are lined with the tested material. The 2.4~m wide front wall has a 0.8~m wide, 2.0~m high door opening in the middle, and the back corner adjacent to the wall opposite to the door has a gas burner with an output of 100~kW for the first ten minutes of the test. If flashover does not occur during this time, the burner heat release is increased to 300~kW for another ten minutes. The test is terminated after 20 minutes or when flashover occurs.
17+
In a room corner test, the walls and ceiling of a 2.4~m by 3.6~m by 2.4~m tall compartment are lined with the sample material. The 2.4~m wide front wall has a 0.8~m wide, 2.0~m high door opening in the middle, and the back corner adjacent to the wall opposite to the door has a gas burner with an output of 100~kW for 10~min. If flashover does not occur during this time, the burner heat release is increased to 300~kW for another 10~min. The test is terminated after 20~min or when flashover occurs.
1818

19-
In the test M12~\cite{Sundstrom:1998} the walls and the ceiling are lined with 10~mm thick spruce timber, whereas in the modified test~\cite{Hietaniemi:2001} the walls and ceiling only in the proximity of the gas burner was lined with 22~mm thick spruce timber, elsewhere being exposed concrete. In both of the studied tests the flashover occurred within the first ten minutes, so the burner heat release was 100~kW over the entire duration of the test.
19+
In the test M12~\cite{Sundstrom:1998}, the walls and the ceiling are lined with 10~mm thick spruce timber, whereas in the modified test~\cite{Hietaniemi:2001} the walls and ceiling only in the proximity of the gas burner were lined with 22~mm thick spruce timber, elsewhere being exposed concrete. In both tests, flashover occurred within the first 10~min.
2020

2121
\subsubsection{Modeling Notes}
2222

23-
In the simulations of flaming cone calorimeter experiments, the computational domain is 20~cm in length and width and 30~cm in height, with a grid resolution of 1~cm. For the inert and smoldering experiments the domain is 30~cm by 30~cm by 30~cm with the resolution of 10~cm, as the gaseous phase was not directly modeled. In the simulations of room corner tests, the computational domain continued for 1.5~m outside after the 20~cm thick front wall. In the corner where the burner was located, the grid resolution was 2~cm over a volume that extends for 0.5~m from the corner at x-, and y directions, and spans the entire height from floor to the ceiling. Elsewhere, a resolution of 10~cm was used.
23+
In the simulations of flaming cone calorimeter experiments, the computational domain is 20~cm in length and width and 30~cm in height, with a grid resolution of 1~cm. For the inert and smoldering experiments, the domain is 30~cm by 30~cm by 30~cm with a resolution of 10~cm because the gaseous phase is not directly modeled. In the simulations of room corner tests, the computational domain is extended 1.5~m outside the 20~cm thick front wall. In the corner where the burner is located, the grid resolution is 2~cm over a volume that extends 0.5~m from the corner and spans the entire height from floor to ceiling. Elsewhere, a resolution of 10~cm is used.
2424

25-
Neither~\cite{Sundstrom:1998} nor~\cite{Hietaniemi:2001} reported the initial moisture content of the timbers, so a moisture mass fraction of 0.09 similar to the cone calorimeter specimens in~\cite{Rinta-Paavola:2023}, was assumed in the beginning of room corner tests. This is a relevant moisture content for timber in normal room air.
25+
Neither Ref.~\cite{Sundstrom:1998} nor Ref.~\cite{Hietaniemi:2001} reports the initial moisture content of the timbers; thus, a moisture mass fraction of 0.09 is assumed, similar to the cone calorimeter specimens in Ref.~\cite{Rinta-Paavola:2023}.
2626

27-
The parameters in the original research of Rinta-Paavola et al.~\cite{Rinta-Paavola:2023} were estimated for FDS 6.7.7, but in FDS 6.7.8 and later the reaction rate equation for solid pyrolysis is different, possibly requiring scaling of the pre-exponential factor when parameters are carried from a FDS version with the old reaction rate definition to a newer one. See~\cite{Rinta-Paavola:2023} for greater detail. Table~\ref{Aalto_Woods_Properties} presents appropriately scaled parameters that yield results consistent to the original simulations in~\cite{Rinta-Paavola:2023} when initial moisture mass fraction is 0.09. If a different initial moisture content is used, re-scaling of the frequency factor would be necessary to gain most accurate possible agreement with the model presented in~\cite{Rinta-Paavola:2023}.
27+
The parameters in the original research of Rinta-Paavola~et~al.~\cite{Rinta-Paavola:2023} were estimated for FDS 6.7.7, but in FDS 6.7.8 and later the reaction rate equation for solid pyrolysis is different, possibly requiring scaling of the pre-exponential factor when parameters are carried from a FDS version with the old reaction rate definition to a newer one. See Ref.~\cite{Rinta-Paavola:2023} for details. Table~\ref{Aalto_Woods_Properties} presents appropriately scaled parameters that yield results consistent with the original simulations in Ref.~\cite{Rinta-Paavola:2023} when the initial moisture mass fraction is 0.09. If a different initial moisture content is used, re-scaling of the frequency factor would be necessary to obtain the most accurate possible agreement with the model presented in Ref.~\cite{Rinta-Paavola:2023}.
2828

2929

3030
\section{ArupFire Tunnel Fire Experiments}
3131
\label{ArupFire_Tunnel_Fire_Description}
3232

33-
Gabriele Vigne and Jimmy J\"{o}nsson of ArupFire conducted a series of fire experiments within a tunnel with a 50~m$^2$ cross section. The tunnel is located in La Ribera del Folgoso, Spain. It is approximately 6.5~m high, 8~m wide and 300~m long. Five replicate tests were conducted using a 1~m by 2~m steel pan filled with heptane on water. Near-ceiling temperatures were measured 2~m, 4~m, 6~m and 8~m from the plume centerline. The peak heat release rate was approximately 5.3~MW.
33+
Gabriele Vigne and Jimmy J\"{o}nsson~\cite{Vigne:2009} conducted a series of fire experiments within a tunnel with a 50~m$^2$ cross section. The tunnel is located in La Ribera del Folgoso, Spain. It is approximately 6.5~m high, 8~m wide and 300~m long. Five replicate tests were conducted using a 1~m by 2~m steel pan filled with heptane on water. Near-ceiling temperatures were measured 2~m, 4~m, 6~m and 8~m from the plume centerline. The peak heat release rate was approximately 5.3~MW.
3434

3535

3636
\section{ATF Corridors Experiments}

Source/radi.f90

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -4443,6 +4443,7 @@ SUBROUTINE RADIATION_FVM
44434443
IF (WC%BOUNDARY_TYPE/=SOLID_BOUNDARY) CYCLE
44444444
B1 => BOUNDARY_PROP1(WC%B1_INDEX)
44454445
SF => SURFACE(WC%SURF_INDEX)
4446+
IF (SF%SKIP_INRAD) INRAD_W(IW) = 0._EB
44464447
IF (SF%EXTERNAL_FLUX > TWO_EPSILON_EB) THEN
44474448
IF (ABS(T_BEGIN) <= SPACING(B1%T_IGN)) THEN
44484449
TSI = T
@@ -4461,6 +4462,7 @@ SUBROUTINE RADIATION_FVM
44614462
IF (CFACE(ICF)%BOUNDARY_TYPE==NULL_BOUNDARY) CYCLE
44624463
B1 => BOUNDARY_PROP1(CFA%B1_INDEX)
44634464
SF => SURFACE(CFA%SURF_INDEX)
4465+
IF (SF%SKIP_INRAD) INRAD_W(ICF) = 0._EB
44644466
IF (SF%EXTERNAL_FLUX > TWO_EPSILON_EB) THEN
44654467
IF (ABS(T_BEGIN) <= SPACING(B1%T_IGN)) THEN
44664468
TSI = T

Source/read.f90

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -8021,7 +8021,7 @@ SUBROUTINE READ_SURF(QUICK_READ)
80218021
VEG_LSET_WIND_EXP,VEG_LSET_BETA,VEG_LSET_HT,VEG_LSET_SIGMA,VEG_LSET_ROS_00, &
80228022
VEG_LSET_M1,VEG_LSET_M10,VEG_LSET_M100,VEG_LSET_MLW,VEG_LSET_MLH,VEG_LSET_SURF_LOAD,VEG_LSET_FIREBASE_TIME,&
80238023
VEG_LSET_CHAR_FRACTION,VEG_LSET_WIND_HEIGHT,VEL_PART,INIT_PER_AREA,TIME_STEP_FACTOR
8024-
LOGICAL :: DEFAULT,VEG_LSET_SPREAD,VEG_LSET_TAN2,TGA_ANALYSIS,VARIABLE_THICKNESS,HT3D,THERM_THICK,VEG_LSET_ROS_FIXED
8024+
LOGICAL :: DEFAULT,VEG_LSET_SPREAD,VEG_LSET_TAN2,TGA_ANALYSIS,VARIABLE_THICKNESS,HT3D,THERM_THICK,VEG_LSET_ROS_FIXED,SKIP_INRAD
80258025
LOGICAL, ALLOCATABLE, DIMENSION(:) :: DUPLICATE
80268026
! Ember generating variables
80278027
REAL(EB) :: EMBER_GENERATION_HEIGHT(2),EMBER_IGNITION_POWER_MEAN,EMBER_IGNITION_POWER_SIGMA,EMBER_TRACKING_RATIO,EMBER_YIELD
@@ -8050,7 +8050,7 @@ SUBROUTINE READ_SURF(QUICK_READ)
80508050
RAMP_V,RAMP_V_X,RAMP_V_Y,RAMP_V_Z,&
80518051
REFERENCE_HEAT_FLUX,REFERENCE_HEAT_FLUX_TIME_INTERVAL,MINIMUM_SCALING_HEAT_FLUX,MAXIMUM_SCALING_HEAT_FLUX,&
80528052
REFERENCE_THICKNESS,REMESH_RATIO,&
8053-
RGB,ROUGHNESS,SHAPE_FACTOR,SPEC_ID,&
8053+
RGB,ROUGHNESS,SHAPE_FACTOR,SKIP_INRAD,SPEC_ID,&
80548054
SPREAD_RATE,STRETCH_FACTOR,SURFACE_VOLUME_RATIO,&
80558055
TAU_EF,TAU_MF,TAU_PART,TAU_Q,TAU_T,TAU_V,TEXTURE_HEIGHT,TEXTURE_MAP,TEXTURE_WIDTH,&
80568056
TGA_ANALYSIS,TGA_CONVERSION_FACTOR,TGA_DT,TGA_DUMP,TGA_FINAL_TEMPERATURE,TGA_HEATING_RATE,THICKNESS,&
@@ -8577,6 +8577,7 @@ SUBROUTINE READ_SURF(QUICK_READ)
85778577
SF%NO_SLIP = NO_SLIP
85788578
SF%FYI = FYI
85798579
SF%EXTERNAL_FLUX = 1000._EB*EXTERNAL_FLUX
8580+
SF%SKIP_INRAD = SKIP_INRAD
85808581
SF%HORIZONTAL = HORIZONTAL
85818582
SF%INIT_IDS = INIT_IDS
85828583
IF (INIT_IDS(1)/='null') INIT_INVOKED_BY_SURF = .TRUE.
@@ -9266,6 +9267,7 @@ SUBROUTINE SET_SURF_DEFAULTS
92669267
ROUGHNESS = -1._EB !4.5E-5_EB ! meters, commercial steel
92679268
ENDIF
92689269
SHAPE_FACTOR = 0.25_EB
9270+
SKIP_INRAD = .FALSE.
92699271
SPEC_ID = 'null'
92709272
SPREAD_RATE = -1._EB
92719273
STRETCH_FACTOR = 2._EB

Source/type.f90

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -998,6 +998,7 @@ MODULE TYPES
998998
LOGICAL :: INERT_Q_REF !< Treat REFERENCE_HEAT_FLUX as an inert atmosphere test
999999
LOGICAL :: ALLOW_UNDERSIDE_PARTICLES=.FALSE. !< Allow droplets to move along downward facing surfaces
10001000
LOGICAL :: ALLOW_SURFACE_PARTICLES=.TRUE. !< Allow particles to live on a solid surface
1001+
LOGICAL :: SKIP_INRAD = .FALSE. !< Only apply external flux to a surface
10011002
INTEGER :: GEOMETRY,BACKING,PROFILE,HEAT_TRANSFER_MODEL=0,NEAR_WALL_TURB_MODEL=5
10021003
CHARACTER(LABEL_LENGTH) :: PART_ID
10031004
CHARACTER(LABEL_LENGTH) :: ID,TEXTURE_MAP,LEAK_PATH_ID(2)

0 commit comments

Comments
 (0)