Skip to content

Commit 4a095ed

Browse files
authored
Merge pull request #15590 from mcgratta/master
FDS Source: Adjust output of TWO_D and CYLINDRICAL
2 parents 6a17a7d + 32bf00c commit 4a095ed

File tree

9 files changed

+65
-51
lines changed

9 files changed

+65
-51
lines changed

Manuals/FDS_User_Guide/FDS_User_Guide.tex

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -991,12 +991,11 @@ \subsection{Basics}
991991
\subsection{Two-Dimensional and Axially-Symmetric Calculations}
992992
\label{info:2D}
993993

994-
The governing equations solved in FDS are written in terms of a three dimensional Cartesian coordinate system. However, a two dimensional Cartesian or two dimensional cylindrical (axially-symmetric) calculation can be performed by setting the \ct{J} in the \ct{IJK} triplet to 1 on the \ct{MESH} line. For axial symmetry, add \ct{CYLINDRICAL=T} to the \ct{MESH} line, and the coordinate $x$ is then interpreted as the radial coordinate $r$. If more than one mesh is used, all the meshes must be specified as 2-D or \ct{CYLINDRICAL}---you cannot mix 2-D, 3-D and cylindrical geometries. No boundary conditions should be set at the planes $y=\hbox{\tt YMIN=XB(3)}$ or $y=\hbox{\tt YMAX=XB(4)}$, nor at $r=\hbox{\tt XMIN=XB(1)}$ in an axially-symmetric calculation if $r=\hbox{XB(1)=0}$ (Note that \ct{XB(1)} does not have to be 0). For better visualizations, the difference between \ct{XB(4)} and \ct{XB(3)} should be small so that the Smokeview rendering appears to be in 2-D. An example of an axially-symmetric helium plume is given in Sec.~\ref{baroclinic_torque}.
995-
996-
When processing results for a \ct{CYLINDRICAL} simulation, note that integrated output quantities with the \ct{SPATIAL_STATISTIC} attribute apply only to the specified 2-D or cylindrical coordinates. Thus, the cylindrical coordinates define a cylindrical sector, like a slice of cake, even though Smokeview will not render it this way. The fully integrated quantity can be calculated by multiplying the reported value by $2 \pi \, \delta\theta$, where $\delta\theta$ is the difference between \ct{YMAX} and \ct{YMIN} in radians. The values chosen for \ct{YMAX} and \ct{YMIN} do not matter as long as the rendering in Smokeview is to your liking.
994+
The governing equations solved in FDS are written in terms of a three-dimensional Cartesian coordinate system. However, a two-dimensional Cartesian or two-dimensional cylindrical (axially-symmetric) calculation can be performed by setting the \ct{J} in the \ct{IJK} triplet to 1 on the \ct{MESH} line. For axial symmetry, add \ct{CYLINDRICAL=T} to the \ct{MESH} line, and the coordinate $x$ is then interpreted as the radial coordinate $r$. If more than one mesh is used, all the meshes must be specified as 2-D or \ct{CYLINDRICAL}---you cannot mix 2-D, 3-D and cylindrical geometries. No boundary conditions should be set at the planes $y=\hbox{\tt YMIN=XB(3)}$ or $y=\hbox{\tt YMAX=XB(4)}$, nor at $r=\hbox{\tt XMIN=XB(1)}$ in an axially-symmetric calculation if $r=\hbox{XB(1)=0}$ (Note that \ct{XB(1)} does not have to be 0). For better visualizations, the difference between \ct{XB(4)} and \ct{XB(3)} should be small so that the Smokeview rendering appears to be in 2-D. An example of an axially-symmetric helium plume is given in Sec.~\ref{baroclinic_torque}.
997995

998996
When performing solid phase heat transfer while using a 2-D \ct{CYLINDRICAL} coordinate system, you must designate \ct{GEOMETRY='CYLINDRICAL'} on a surface (\ct{SURF} line) that is facing radially outward (positive $r$ direction) or \ct{GEOMETRY='INNER CYLINDRICAL'} on a surface that is facing radially inward (negative $r$ direction). In the latter instance, you must also specify the \ct{INNER_RADIUS} (m) of the cylinder. For the outer cylindrical boundary, specify an \ct{INNER_RADIUS} if appropriate. Its default value is 0~m. Because your inward and outward facing boundaries might occur at various radii, you must create separate \ct{SURF} lines for each with the appropriate values of \ct{GEOMETRY} and \ct{INNER_RADIUS}. For an obstruction (\ct{OBST}), use \ct{SURF_ID6} to assign individual \ct{SURF ID}s to each of the six faces. Because this is a 2-D simulation, the third and fourth entries representing the \ct{SURF ID}s in the $y$ or angular direction can just be designated \ct{'INERT'}.
999997

998+
When processing results for a \ct{CYLINDRICAL} simulation, note that integrated output quantities with the \ct{SPATIAL_STATISTIC} attribute refer to the volume or surface area of the entire cylinder, not just the wedge. Smokeview renders the wedge as a 2-D slice. The values chosen for \ct{YMAX} and \ct{YMIN} do not matter as long as the rendering in Smokeview is to your liking. For a 2-D, non-cylindrical geometry, spatially integrated quanties shall be output in units of the quantity per unit meter.
1000999

10011000
\subsection{Multiple Meshes}
10021001
\label{info:multimesh}
@@ -9961,6 +9960,9 @@ \subsubsection{Limiting the Integration}
99619960
\end{lstlisting}
99629961
would output the total surface area in the volume \ct{XB} where the total heat flux exceeds 10~\unit{kW/m^2}.
99639962

9963+
\subsubsection{Two-Dimensional and Cylindrical Coordinate Systems}
9964+
9965+
If the computational domain is two-dimensional or cylindrical, some spatially-integrated quantities are adjusted to eliminate the dependence on the arbitrarily chosen $\delta y$ or $\delta \theta$. For a 2-D domain, a reported volume output will have units of m$^3$/m and an area output will have units of m$^2$/m. For a cylindrical domain, volume and area outputs are reported for the entire cylinder rather than the thin wedge on which the simulation is performed. This might cause confusion in cases where both the input parameters and simulation results involve volumetric or areal quantities; thus, it is good practice to perform a simple test case with a known result to verify that these adjustments have been performed properly.
99649966

99659967

99669968
\subsection{Temporally-Integrated Outputs}

Source/dump.f90

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -6741,6 +6741,12 @@ SUBROUTINE UPDATE_DEVICES_1(T,DT,NM)
67416741
IF (ICF>0) VOL = VOL*CUT_FACE(ICF)%ALPHA_CF
67426742
ENDIF
67436743

6744+
IF (CYLINDRICAL) THEN
6745+
VOL = TWOPI*VOL/DY(J)
6746+
ELSEIF (TWO_D) THEN
6747+
VOL = VOL/DY(J)
6748+
ENDIF
6749+
67446750
VALUE = GAS_PHASE_OUTPUT(T,DT,NM,I,J,K,DV%QUANTITY_INDEX(1),0,DV%Y_INDEX,DV%Z_INDEX,DV%ELEM_INDEX,&
67456751
DV%PART_CLASS_INDEX,DV%VELO_INDEX,DV%PIPE_INDEX,DV%PROP_INDEX,DV%REAC_INDEX,&
67466752
DV%MATL_INDEX)

Source/part.f90

Lines changed: 7 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1084,8 +1084,13 @@ SUBROUTINE INSERT_VOLUMETRIC_PARTICLES
10841084
ENDIF
10851085
EXIT GEOM_SELECT
10861086
ENDIF
1087-
INSERT_VOLUME = (X2-X1)*(Y2-Y1)*(Z2-Z1)
1088-
INPUT_VOLUME = (IN_X2-IN_X1)*(IN_Y2-IN_Y1)*(IN_Z2-IN_Z1)
1087+
IF (CYLINDRICAL) THEN
1088+
INSERT_VOLUME = 0.5_EB*(X2**2-X1**2)*(Y2-Y1)*(Z2-Z1)
1089+
INPUT_VOLUME = 0.5_EB*(IN_X2**2-IN_X1**2)*(IN_Y2-IN_Y1)*(IN_Z2-IN_Z1)
1090+
ELSE
1091+
INSERT_VOLUME = (X2-X1)*(Y2-Y1)*(Z2-Z1)
1092+
INPUT_VOLUME = (IN_X2-IN_X1)*(IN_Y2-IN_Y1)*(IN_Z2-IN_Z1)
1093+
ENDIF
10891094
CASE('CONE','CYLINDER')
10901095
X0 = 0.5_EB*(IN_X1+IN_X2)
10911096
Y0 = 0.5_EB*(IN_Y1+IN_Y2)

Utilities/Matlab/FDS_verification_dataplot_inputs.csv

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -49,12 +49,12 @@ d,box_burn_away4,Fires/box_burn_away4_git.txt,Fires/box_burn_away.csv,1,2,Time,M
4949
d,box_burn_away5,Fires/box_burn_away5_git.txt,Fires/box_burn_away.csv,1,2,Time,Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away5_devc.csv,2,3,Time,Mass fuel,FDS (fuel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away5),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away5,Relative Error,end,0.02,Fires,r*,r,TeX
5050
d,box_burn_away6,Fires/box_burn_away6_git.txt,Fires/box_burn_away.csv,1,2,Time,Mass_residue (kg)|Mass (kg),Ideal Fuel 1|Ideal Fuel 2,ko|ro,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away6_devc.csv,2,3,Time,Mass fuel 1|Mass fuel 2,FDS (fuel 1)|FDS (fuel 2),k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away6),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away6,Relative Error,end,0.02,Fires,r*,r,TeX
5151
d,box_burn_away7,Fires/box_burn_away7_git.txt,Fires/box_burn_away7.csv,1,2,Time,Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away7_devc.csv,2,3,Time,Mass fuel,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away7),Time (s),Mass (kg),0,20,1,0,0.04,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away7,Relative Error,end,0.02,Fires,r*,r,TeX
52-
d,box_burn_away_2D,Fires/box_burn_away_2D_git.txt,Fires/box_burn_away.csv,1,2,Time,Mass (kg)| Mass_end (kg),Ideal (fuel)|Ideal (solid),ko|ro,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away_2D_devc.csv,2,3,Time,Mass fuel|Mass solid,FDS (fuel)|FDS(solid),k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away_2D),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away_2D,Relative Error,end,0.02,Fires,r*,r,TeX
52+
d,box_burn_away_2D,Fires/box_burn_away_2D_git.txt,Fires/box_burn_away_2D.csv,1,2,Time,Mass (kg)| Mass_end (kg),Ideal (fuel)|Ideal (solid),ko|ro,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away_2D_devc.csv,2,3,Time,Mass fuel|Mass solid,FDS (fuel)|FDS(solid),k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away_2D),Time (s),Mass (kg),0,30,1,0,4.0,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away_2D,Relative Error,end,0.02,Fires,r*,r,TeX
5353
d,box_burn_away8,Fires/box_burn_away8_git.txt,Fires/box_burn_away.csv,1,2,Time,Mass4 (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away8_devc.csv,2,3,Time,Mass fuel,FDS (fuel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away8),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away8,Relative Error,end,0.02,Fires,r*,r,TeX
5454
d,box_burn_away9,Fires/box_burn_away9_git.txt,Fires/box_burn_away9.csv,1,2,Time,Mass,Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away9_devc.csv,2,3,Time,Mass fuel,FDS (fuel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away9),Time (s),Mass (kg),0,180,1,0,1,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away9,Relative Error,end,0.01,Fires,r*,r,TeX
5555
d,box_burn_away10,Fires/box_burn_away10_git.txt,Fires/box_burn_away10.csv,1,2,Time,Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away10_devc.csv,2,3,Time,Mass fuel,FDS (fuel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away10),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away10,Relative Error,end,0.02,Fires,r*,r,TeX
5656
d,box_burn_away11,Fires/box_burn_away11_git.txt,Fires/box_burn_away11.csv,1,2,Time,Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away11_devc.csv,2,3,Time,Mass fuel,FDS (fuel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away11),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away11,Relative Error,end,0.01,Fires,r*,r,TeX
57-
d,box_burn_away_2D_residue,Fires/box_burn_away_2D_residue_git.txt,Fires/box_burn_away.csv,1,2,Time,Mass_residue (kg)| Mass_end (kg),Ideal (fuel)|Ideal (solid),ko|ro,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away_2D_residue_devc.csv,2,3,Time,Mass fuel|Mass solid,FDS (fuel)|FDS(solid),k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away_2D_residue),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.91,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away_2D_residue,Relative Error,end,0.02,Fires,r*,r,TeX
57+
d,box_burn_away_2D_residue,Fires/box_burn_away_2D_residue_git.txt,Fires/box_burn_away_2D.csv,1,2,Time,Mass_residue (kg)| Mass_end (kg),Ideal (fuel)|Ideal (solid),ko|ro,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away_2D_residue_devc.csv,2,3,Time,Mass fuel|Mass solid,FDS (fuel)|FDS(solid),k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box_burn_away_2D_residue),Time (s),Mass (kg),0,30,1,0,4.0,1,no,0.05 0.91,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away_2D_residue,Relative Error,end,0.02,Fires,r*,r,TeX
5858
d,bucket_test_1,Sprinklers_and_Sprays/bucket_test_1_git.txt,Sprinklers_and_Sprays/bucket_test_1.csv,1,2,Time (s),Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/bucket_test_1_devc.csv,2,3,Time,Mass,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Accumulated Mass (bucket_test_1),Time (s),Mass (kg),0,40,1,0,110,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/bucket_test_1,Relative Error,end,0.02,Sprinklers and Sprays,bs,b,TeX
5959
d,bucket_test_2,Sprinklers_and_Sprays/bucket_test_2_git.txt,Sprinklers_and_Sprays/bucket_test_2.csv,1,2,Time (s),Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/bucket_test_2_devc.csv,2,3,Time,Mass,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Accumulated Mass (bucket_test_2),Time (s),Mass (kg),0,15,1,0,6,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/bucket_test_2,Relative Error,end,0.03,Sprinklers and Sprays,bs,b,TeX
6060
d,bucket_test_3,Sprinklers_and_Sprays/bucket_test_3_git.txt,Sprinklers_and_Sprays/bucket_test_3.csv,1,2,Time (s),Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/bucket_test_3_devc.csv,2,3,Time,Mass,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Accumulated Mass (bucket_test_3),Time (s),Mass (kg),0,10,1,0,0.01,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/bucket_test_3,Relative Error,end,0.02,Sprinklers and Sprays,bs,b,TeX
Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,3 +1,3 @@
11
"Time","div","Q_CONV","Q_DIFF","Q_SUM"
2-
1.5,-1.357e-8,1.16e-5,-1.16e-5,0
3-
2.0,-1.357e-8,1.16e-5,-1.16e-5,0
2+
1.5,-1.357e-5,1.16e-5,-1.16e-5,0
3+
2.0,-1.357e-5,1.16e-5,-1.16e-5,0
Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,3 @@
1+
Time,Mass (kg),Mass_residue (kg),Mass_end (kg)
2+
28,3.2,1.6,0.0
3+
30,3.2,1.6,0.0

Verification/Fires/box_burn_away_2D.fds

Lines changed: 8 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -1,17 +1,17 @@
11
&HEAD CHID='box_burn_away_2D', TITLE='Test BURN_AWAY feature' /
22

33
The FOAM box is evaporated away by the high thermal radiation
4-
from HOT surfaces. The mass of the box is 0.4^3 m3 * 20 kg/m3 = 1.28 kg.
4+
from HOT surfaces. The mass of the box is 0.4^2 m2 * 20 kg/m3 = 3.2 kg/m.
55
This should be compared to the final value of fuel density volume integral,
66
computed by the first DEVC.
77

88
2D version of box_burn_away
99

1010
The gas species is mixture fraction fuel.
1111

12-
&MESH IJK=10,1,10 XB=0.0,1.0,0.3,0.7,0.0,1.0 /
12+
&MESH IJK=10,1,10 XB=0.0,1.0,0.0,0.01,0.0,1.0 /
1313

14-
&TIME T_END=30. DT = 0.01/
14+
&TIME T_END=30., DT=0.01 /
1515

1616
&MATL ID = 'FOAM'
1717
HEAT_OF_REACTION = 800.
@@ -30,11 +30,11 @@ The gas species is mixture fraction fuel.
3030

3131
&REAC FUEL='METHANE',AUTO_IGNITION_TEMPERATURE=15000. /
3232

33-
&DUMP SMOKE3D=.FALSE., MASS_FILE=.TRUE. /
33+
&DUMP SMOKE3D=.FALSE., MASS_FILE=T /
3434

35-
&OBST XB=0.30,0.70,0.0,1.0,0.30,0.70, SURF_ID='FOAM SLAB', BULK_DENSITY=20., MATL_ID='FOAM' /
35+
&OBST XB=0.30,0.70,0.0,0.01,0.30,0.70, SURF_ID='FOAM SLAB', BULK_DENSITY=20., MATL_ID='FOAM' /
3636

37-
&SURF ID = 'HOT' TMP_FRONT = 1500., DEFAULT = .TRUE./
37+
&SURF ID='HOT', TMP_FRONT=1500., DEFAULT=T /
3838

3939
&BNDF QUANTITY='WALL TEMPERATURE' /
4040
&BNDF QUANTITY='BURNING RATE' /
@@ -44,9 +44,8 @@ The gas species is mixture fraction fuel.
4444
&SLCF PBY=0.5, QUANTITY='INTEGRATED INTENSITY' /
4545
&SLCF PBY=0.5, QUANTITY='MASS FRACTION',SPEC_ID='METHANE' /
4646

47-
&DEVC XB=0.0,1.0,0.3,0.7,0.0,1.0, QUANTITY = 'DENSITY',SPEC_ID='METHANE', SPATIAL_STATISTIC='VOLUME INTEGRAL' ID = 'Mass fuel'/
48-
&DEVC XB=0.0,1.0,0.3,0.7,0.0,1.0, QUANTITY = 'HRRPUV', SPATIAL_STATISTIC='VOLUME INTEGRAL', ID = 'HRR' /
49-
&DEVC XB=0.3,0.7,0.3,0.7,0.3,0.7, SPATIAL_STATISTIC='SURFACE INTEGRAL', QUANTITY = 'SURFACE DENSITY', ID = 'Mass solid' /
47+
&DEVC XB=0.0,1.0,0.0,0.01,0.0,1.0, QUANTITY='DENSITY', SPEC_ID='METHANE', SPATIAL_STATISTIC='VOLUME INTEGRAL', ID='Mass fuel'/
48+
&DEVC XB=0.3,0.7,0.0,0.01,0.7,0.7, IOR=3, SPATIAL_STATISTIC='SURFACE INTEGRAL', QUANTITY='SURFACE DENSITY', ID='Mass solid' /
5049

5150
&TAIL /
5251

Verification/Fires/box_burn_away_2D_residue.fds

Lines changed: 21 additions & 22 deletions
Original file line numberDiff line numberDiff line change
@@ -1,20 +1,20 @@
11
&HEAD CHID='box_burn_away_2D_residue', TITLE='Test BURN_AWAY feature' /
22

33
The FOAM box is evaporated away by the high thermal radiation
4-
from HOT surfaces. The mass of the box is 0.4^3 m3 * 20 kg/m3 = 1.28 kg.
4+
from HOT surfaces. The mass of the box is 0.4^2 m2 * 20 kg/m3 = 3.2 kg/m.
55
Half of the mass is converted to fuel and other half to residue.
66
Bulk density defines the combustible mass of the box:
7-
10 kg/m3 x 0.064 m3 = 0.64 kg.
7+
10 kg/m3 x 0.16 m2 = 1.6 kg/m.
88
This should be compared to the final value of fuel density volume integral,
99
computed by the first DEVC.
1010

1111
2D version of box_burn_away
1212

1313
The gas species is mixture fraction fuel.
1414

15-
&MESH IJK=10,1,10 XB=0.0,1.0,0.3,0.7,0.0,1.0 /
15+
&MESH IJK=10,1,10 XB=0.0,1.0,0.0,0.01,0.0,1.0 /
1616

17-
&TIME T_END=30. DT = 0.01/
17+
&TIME T_END=30. DT=0.01 /
1818

1919
&MATL ID = 'FOAM'
2020
HEAT_OF_REACTION = 100.
@@ -23,30 +23,30 @@ The gas species is mixture fraction fuel.
2323
DENSITY = 20.
2424
N_REACTIONS = 1
2525
NU_SPEC = 0.5
26-
NU_MATL = 0.5
26+
NU_MATL = 0.5
2727
SPEC_ID = 'METHANE'
28-
MATL_ID = 'RESIDUE'
28+
MATL_ID = 'RESIDUE'
2929
REFERENCE_TEMPERATURE= 100. /
3030

31-
&MATL ID = 'RESIDUE'
32-
CONDUCTIVITY = 10
33-
SPECIFIC_HEAT = 1.0
34-
DENSITY = 100
35-
EMISSIVITY = 1.0 /
31+
&MATL ID = 'RESIDUE'
32+
CONDUCTIVITY = 10
33+
SPECIFIC_HEAT = 1.0
34+
DENSITY = 100
35+
EMISSIVITY = 1.0 /
3636

37-
&SURF ID = 'FOAM SLAB'
38-
COLOR = 'TOMATO 3'
39-
VARIABLE_THICKNESS = T
40-
STRETCH_FACTOR = 1.0
41-
BURN_AWAY = T /
37+
&SURF ID = 'FOAM SLAB'
38+
COLOR = 'TOMATO 3'
39+
VARIABLE_THICKNESS = T
40+
STRETCH_FACTOR = 1.0
41+
BURN_AWAY = T /
4242

4343
&REAC FUEL='METHANE', AUTO_IGNITION_TEMPERATURE=15000. /
4444

45-
&DUMP SMOKE3D=.FALSE., MASS_FILE=.TRUE. /
45+
&DUMP SMOKE3D=.FALSE., MASS_FILE=T /
4646

47-
&OBST XB=0.30,0.70,0.0,1.0,0.30,0.70, SURF_ID='FOAM SLAB', BULK_DENSITY=10., MATL_ID='FOAM' /
47+
&OBST XB=0.30,0.70,0.0,0.01,0.30,0.70, SURF_ID='FOAM SLAB', BULK_DENSITY=10., MATL_ID='FOAM' /
4848

49-
&SURF ID = 'HOT' TMP_FRONT = 1100., DEFAULT = .TRUE./
49+
&SURF ID='HOT', TMP_FRONT=1100., DEFAULT=T /
5050

5151
&BNDF QUANTITY='WALL TEMPERATURE' /
5252
&BNDF QUANTITY='BURNING RATE' /
@@ -57,9 +57,8 @@ The gas species is mixture fraction fuel.
5757
&SLCF PBY=0.5, QUANTITY='DENSITY' /
5858
&SLCF PBY=0.5, QUANTITY='MASS FRACTION',SPEC_ID='METHANE' /
5959

60-
&DEVC XB=0.0,1.0,0.3,0.7,0.0,1.0, QUANTITY = 'DENSITY',SPEC_ID='METHANE', SPATIAL_STATISTIC='VOLUME INTEGRAL' ID = 'Mass fuel'/
61-
&DEVC XB=0.0,1.0,0.3,0.7,0.0,1.0, QUANTITY = 'HRRPUV', SPATIAL_STATISTIC='VOLUME INTEGRAL', ID = 'HRR' /
62-
&DEVC XB=0.3,0.7,0.3,0.7,0.3,0.7, SPATIAL_STATISTIC='SURFACE INTEGRAL', QUANTITY = 'SURFACE DENSITY' , ID = 'Mass solid'/
60+
&DEVC XB=0.0,1.0,0.0,0.01,0.0,1.0, QUANTITY='DENSITY', SPEC_ID='METHANE', SPATIAL_STATISTIC='VOLUME INTEGRAL', ID='Mass fuel' /
61+
&DEVC XB=0.3,0.7,0.0,0.01,0.7,0.7, IOR=3, SPATIAL_STATISTIC='SURFACE INTEGRAL', QUANTITY='SURFACE DENSITY', ID='Mass solid' /
6362

6463
&TAIL /
6564

0 commit comments

Comments
 (0)