Skip to content

Commit d8ae3b1

Browse files
committed
Merge branch 'master' into fl_mw_cor
2 parents 92299c9 + 9aeb558 commit d8ae3b1

File tree

10 files changed

+88
-69
lines changed

10 files changed

+88
-69
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: 23 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,7 @@
11

22
# Make Smokeview images for the FDS Manuals.
3-
# The script uses the command "smokeview" unless you provide an optional name.
3+
# The script gets the smokeview executable appropriate for your OS from the smv repository
4+
# unless you provide an optional full path to the executable as an argument.
45

56
import pandas as pd
67
import os
@@ -10,9 +11,9 @@
1011
import argparse
1112

1213
parser = argparse.ArgumentParser(description="A script to generate Smokeview images")
13-
parser.add_argument("message", nargs="?", default="smokeview", help="Optional smokeview name")
14+
parser.add_argument("message", nargs="?", default="null", help="Optional smokeview path")
1415
args = parser.parse_args()
15-
shell_command = args.message
16+
smokeview_path = args.message
1617

1718
os_name = platform.system()
1819

@@ -21,6 +22,7 @@
2122
raise FileNotFoundError("xvfb-run is not installed. Please install xvfb package.")
2223

2324
bindir = '../../../smv/Build/for_bundle'
25+
smvdir = '../../../smv/Build/smokeview/'
2426
outdir = '../../Verification/'
2527
original_dir = os.getcwd()
2628

@@ -29,26 +31,29 @@
2931
folder = df[0].values
3032
case = df[1].values
3133

32-
if os_name == "Linux":
33-
result = subprocess.run(["bash", "-i", "-c", f"alias {shell_command} 2>/dev/null | sed -E \"s/alias {shell_command}='(.*)'/\\1/\""],
34-
capture_output=True, text=True)
35-
smokeview_path=result.stdout.strip()
36-
else:
37-
smokeview_path = shutil.which(shell_command)
34+
if smokeview_path != "null":
35+
print("Using "+smokeview_path)
36+
elif os_name == "Linux":
37+
smokeview_path = smvdir + 'intel_linux_64/smokeview_linux_64'
38+
elif os_name == "Darwin":
39+
smokeview_path = smvdir + 'gnu_osx_64/smokeview_osx_64'
40+
elif os_name == "Windows":
41+
smokeview_path = smvdir + 'intel_win_64/smokeview_win_64'
3842

3943
for i in range(len(folder)):
4044
print('generating smokeview image ' + case[i])
4145
os.chdir(outdir + folder[i])
42-
if os_name == "Linux":
43-
try:
46+
try:
47+
if os_name == "Linux":
4448
subprocess.run(['xvfb-run','-w','10','-s','-fp /usr/share/X11/fonts/misc -screen 0 1280x1024x24','-a',smokeview_path,
45-
'-bindir',bindir,'-runscript', case[i] ], check=True, stdout=subprocess.DEVNULL, stderr=subprocess.DEVNULL)
46-
except subprocess.CalledProcessError as e:
47-
print(f"Smokeview failed with return code {e.returncode}")
48-
except FileNotFoundError:
49-
print(f"Smokeview executable not found: {smokeview_path}")
50-
else:
51-
subprocess.run([smokeview_path,'-bindir',bindir,'-runscript',case[i]], stdout=subprocess.DEVNULL, stderr=subprocess.DEVNULL)
49+
'-bindir',bindir,'-runscript', case[i] ], check=True, stdout=subprocess.DEVNULL, stderr=subprocess.DEVNULL)
50+
else:
51+
subprocess.run([smokeview_path,'-bindir',bindir,'-runscript',case[i]],
52+
check=True, stdout=subprocess.DEVNULL, stderr=subprocess.DEVNULL)
53+
except subprocess.CalledProcessError as e:
54+
print(f"Error: Smokeview failed with return code {e.returncode}")
55+
except FileNotFoundError:
56+
print(f"Error: Smokeview executable not found: {smokeview_path}")
5257

5358
os.chdir(original_dir)
5459

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

0 commit comments

Comments
 (0)