Skip to content

add test cases comparing 3D smoke visualization using SMOKE3D_VERSION=1 and SMOKE3D_VERSION=2 #2676

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
gforney merged 3 commits into from
Oct 8, 2026
Merged

add test cases comparing 3D smoke visualization using SMOKE3D_VERSION=1 and SMOKE3D_VERSION=2 #2676

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
4162d34
smv vv: add test cases comparing versions 1 and 2 of smoke3d files
gforney Oct 8, 2026
5bdefcc
smv vv: add s3d version test cases to vv case lsit
gforney Oct 8, 2026
8a2b29e
smv vg: add SMOKE3D_VERSION test cases to verification guide
gforney Oct 8, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
27 changes: 25 additions & 2 deletions Manuals/SMV_Verification_Guide/SMV_Verification_Guide.tex
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1061,8 +1061,8 @@ \section{Vector Slices}

\section{3D Smoke}

Figure \ref{figsmoketest} presents images verifying the display of 3D smoke and fire
(heat release rate per unit volume or HRRPUV). A series of 3D smoke images are drawn
Figure \ref{figsmoketest} shows images verifying the display of 3D smoke and fire
(heat release rate per unit volume or HRRPUV). Three 3D smoke images are drawn
at \SI{1}{s}, \SI{10}{s} and \SI{30}{s}. The images contain semi-transparent slices
derived from both soot density and HRRPUV data.

Expand All @@ -1080,6 +1080,29 @@ \section{3D Smoke}
\label{figsmoketest}%
\end{figure}

Figure \ref{figsmokeversion} shows images verifying the FDS
{\tt \&DUMP SMOKE3D\_VERSION}\ keyword.
Images in the top row were generated using {\tt SMOKE3D\_VERSION=1}.
Images in the bottom row were generated using {\tt SMOKE3D\_VERSION=2}.
The top and bottom sets of images should be the same.

\begin{figure}[bph]
\begin{center}
\begin{tabular}{ccc}
\includegraphics[height=\figheightD]{SCRIPT_FIGURES/plume_s3d1_smoke_01}&
\includegraphics[height=\figheightD]{SCRIPT_FIGURES/plume_s3d1_smoke_10}&
\includegraphics[height=\figheightD]{SCRIPT_FIGURES/plume_s3d1_smoke_30}\\
\includegraphics[height=\figheightD]{SCRIPT_FIGURES/plume_s3d2_smoke_01}&
\includegraphics[height=\figheightD]{SCRIPT_FIGURES/plume_s3d2_smoke_10}&
\includegraphics[height=\figheightD]{SCRIPT_FIGURES/plume_s3d2_smoke_30}
\end{tabular}
\end{center}
\caption
[A test comparing 3D smoke version 1 and 2 format outputs at \SI{1}{s}, \SI{10}{s} and \SI{30}{s}]
{A test comparing 3D smoke version 1 and 2 format outputs at \SI{1}{s}, \SI{10}{s} and \SI{30}{s}}
\label{figsmokeversion}%
\end{figure}

\clearpage

\section{Plot3D}
Expand Down
155 changes: 155 additions & 0 deletions Verification/Visualization/plume_s3d1.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,155 @@
&HEAD CHID='plume_s3d1',TITLE='Plume whirl case' /

same as plume5a except there is a blockage in the middle of the scene to block the flow
The purpose of this case is to demonstrate the curved flow (via streak lines) that results.

&MESH IJK=16,16,32, XB=0.0,1.6,0.0,1.6,0.0,3.2 /

&MISC IBLANK_SMV=.TRUE. /

&DUMP NFRAMES=400 DT_PL3D=8.0, DT_SL3D=0.1, SMOKE3D_VERSION=1 /

&INIT XB=0.2,1.4,0.2,1.4,0.5,2.2 TEMPERATURE=600.0 /

&TIME T_END=40. / Total simulation time

&MATL ID = 'FABRIC'
FYI = 'Properties completely fabricated'
SPECIFIC_HEAT = 1.0
CONDUCTIVITY = 0.1
DENSITY = 100.0
N_REACTIONS = 1
NU_SPEC = 1.
SPEC_ID = 'PROPANE'
REFERENCE_TEMPERATURE = 350.
HEAT_OF_REACTION = 3000.
HEAT_OF_COMBUSTION = 15000. /

&MATL ID = 'FOAM'
FYI = 'Properties completely fabricated'
SPECIFIC_HEAT = 1.0
CONDUCTIVITY = 0.05
DENSITY = 40.0
N_REACTIONS = 1
NU_SPEC = 1.
SPEC_ID = 'PROPANE'
REFERENCE_TEMPERATURE = 350.
HEAT_OF_REACTION = 1500.
HEAT_OF_COMBUSTION = 30000. /

&SURF ID = 'UPHOLSTERY_LOWER'
FYI = 'Properties completely fabricated'
RGB = 151,96,88
BURN_AWAY = .FALSE.
MATL_ID(1:2,1) = 'FABRIC','FOAM'
THICKNESS(1:2) = 0.002,0.1
/

&SURF ID = 'UPHOLSTERY_UPPER'
FYI = 'Properties completely fabricated'
RGB = 151,96,88
BURN_AWAY = .FALSE.
TMP_FRONT = 600.0
/
&REAC SOOT_YIELD=0.01,FUEL='PROPANE'/
&SURF ID='BURNER',HRRPUA=600.0,PART_ID='tracers' / Ignition source

&VENT XB=0.5,1.1,0.5,1.1,0.1,0.1,SURF_ID='BURNER' / fire source on kitchen stove

&OBST XB=0.5,1.1,0.5,1.1,0.0,0.1 /
&OBST XB=0.3,1.3,0.3,1.3,0.4,0.8 SURF_ID='UPHOLSTERY_LOWER'/
&HOLE XB=0.6,1.0,0.2,0.8,0.3,0.9 /
&OBST XB=0.3,1.3,0.3,1.3,1.2,1.6 SURF_ID='UPHOLSTERY_UPPER' /

&VENT MB='XMIN',SURF_ID='OPEN'/
&VENT MB='XMAX',SURF_ID='OPEN'/
&VENT MB='YMIN',SURF_ID='OPEN'/
&VENT MB='YMAX',SURF_ID='OPEN'/
&VENT MB='ZMAX',SURF_ID='OPEN'/

&ISOF QUANTITY='TEMPERATURE',VALUE(1)=100.0 / Show 3D contours of temperature at 100 C
&ISOF QUANTITY='TEMPERATURE',VALUE(1)=200.0 / Show 3D contours of temperature at 200 C
&ISOF QUANTITY='TEMPERATURE',VALUE(1)=620.0 / Show 3D contours of temperature at 620 C
&ISOF QUANTITY='TEMPERATURE',VALUE(1:2)=200.0,400.0 / Show 3D contours of temperature at 200 C
&ISOF QUANTITY='TEMPERATURE' , VALUE(1)=400, QUANTITY2='HRRPUV' /

&PART ID='tracers',MASSLESS=.TRUE.,
QUANTITIES(1:3)='PARTICLE U','PARTICLE V','PARTICLE W'
SAMPLING_FACTOR=10 / Description of massless tracer particles. Apply at a
solid surface with the PART_ID='tracers'

X slices

&SLCF PBX=0.8,QUANTITY='TEMPERATURE',VECTOR=.TRUE. / Add vector slices colored by temperature
&SLCF PBX=0.8,QUANTITY='TEMPERATURE',CELL_CENTERED=.TRUE. /

Y slices

&SLCF PBY=0.8,QUANTITY='VISIBILITY' /
&SLCF PBY=0.8,QUANTITY='TEMPERATURE',VECTOR=.TRUE. /
&SLCF PBY=0.8,QUANTITY='VOLUME FRACTION' SPEC_ID='PROPANE' /
&SLCF PBY=0.8,QUANTITY='VOLUME FRACTION' SPEC_ID='OXYGEN' /
&SLCF PBY=0.8,QUANTITY='VOLUME FRACTION' SPEC_ID='CARBON DIOXIDE' /
&SLCF PBY=0.8,QUANTITY='VOLUME FRACTION' SPEC_ID='CARBON MONOXIDE' /
&SLCF PBY=0.8,QUANTITY='TEMPERATURE',CELL_CENTERED=.TRUE. /


Z slices

&SLCF PBZ=0.4,QUANTITY='TEMPERATURE',VECTOR=.TRUE. /
&SLCF PBZ=0.4,QUANTITY='TEMPERATURE',CELL_CENTERED=.TRUE. /
&SLCF PBZ=1.6,QUANTITY='TEMPERATURE',VECTOR=.TRUE. /
&SLCF PBZ=1.6,QUANTITY='TEMPERATURE',CELL_CENTERED=.TRUE. /
&SLCF PBZ=3.2,QUANTITY='TEMPERATURE',VECTOR=.TRUE. /
&SLCF PBZ=3.2,QUANTITY='TEMPERATURE',CELL_CENTERED=.TRUE. /

3d slices

&SLCF XB=0.0,1.6,0.0,1.6,0.0,3.2,QUANTITY='TEMPERATURE' ,VECTOR=.TRUE. /
&SLCF XB=0.0,1.6,0.0,1.6,0.0,3.2, QUANTITY='TEMPERATURE',CELL_CENTERED=.TRUE. / 3D slice
&SLCF XB=0.0,1.6,0.0,1.6,0.0,3.2,QUANTITY='VOLUME FRACTION' SPEC_ID='CARBON DIOXIDE' /
&SLCF XB=0.0,1.6,0.0,1.6,0.0,3.2,QUANTITY='VOLUME FRACTION' SPEC_ID='CARBON MONOXIDE' /
&SLCF XB=0.0,1.6,0.0,1.6,0.0,3.2,QUANTITY='VOLUME FRACTION' SPEC_ID='OXYGEN' /
&SLCF XB=0.0,1.6,0.0,1.6,0.0,3.2, QUANTITY='DENSITY',SPEC_ID='SOOT',CELL_CENTERED=.TRUE. / 3D slice

&BNDF QUANTITY='GAUGE HEAT FLUX' / Common surface quantities. Good for monitoring fire spread.
&BNDF QUANTITY='BURNING RATE' /
&BNDF QUANTITY='WALL TEMPERATURE' /
&BNDF QUANTITY='WALL TEMPERATURE' CELL_CENTERED=.TRUE. /

&DEVC XYZ=1.2,0.8,0.0 QUANTITY='U-VELOCITY', ID='U-VEL_1' /
&DEVC XYZ=1.2,0.8,0.6 QUANTITY='U-VELOCITY', ID='U-VEL_2' /
&DEVC XYZ=1.2,0.8,1.2 QUANTITY='U-VELOCITY', ID='U-VEL_3' /
&DEVC XYZ=1.2,0.8,1.8 QUANTITY='U-VELOCITY', ID='U-VEL_4' /
&DEVC XYZ=1.2,0.8,2.4 QUANTITY='U-VELOCITY', ID='U-VEL_5' /
&DEVC XYZ=1.2,0.8,3.0 QUANTITY='U-VELOCITY', ID='U-VEL_6' /

&DEVC XYZ=1.2,0.8,0.0 QUANTITY='V-VELOCITY', ID='V-VEL_1' /
&DEVC XYZ=1.2,0.8,0.6 QUANTITY='V-VELOCITY', ID='V-VEL_2' /
&DEVC XYZ=1.2,0.8,1.2 QUANTITY='V-VELOCITY', ID='V-VEL_3' /
&DEVC XYZ=1.2,0.8,1.8 QUANTITY='V-VELOCITY', ID='V-VEL_4' /
&DEVC XYZ=1.2,0.8,2.4 QUANTITY='V-VELOCITY', ID='V-VEL_5' /
&DEVC XYZ=1.2,0.8,3.0 QUANTITY='V-VELOCITY', ID='V-VEL_6' /

&DEVC XYZ=1.2,0.8,0.0 QUANTITY='W-VELOCITY', ID='W-VEL_1' /
&DEVC XYZ=1.2,0.8,0.6 QUANTITY='W-VELOCITY', ID='W-VEL_2' /
&DEVC XYZ=1.2,0.8,1.2 QUANTITY='W-VELOCITY', ID='W-VEL_3' /
&DEVC XYZ=1.2,0.8,1.8 QUANTITY='W-VELOCITY', ID='W-VEL_4' /
&DEVC XYZ=1.2,0.8,2.4 QUANTITY='W-VELOCITY', ID='W-VEL_5' /
&DEVC XYZ=1.2,0.8,3.0 QUANTITY='W-VELOCITY', ID='W-VEL_6' /

&DEVC XYZ=1.2,0.8,0.0 QUANTITY='TEMPERATURE', ID='TEMP_1' /
&DEVC XYZ=1.2,0.8,0.6 QUANTITY='TEMPERATURE', ID='TEMP_2' /
&DEVC XYZ=1.2,0.8,1.2 QUANTITY='TEMPERATURE', ID='TEMP_3' /
&DEVC XYZ=1.2,0.8,1.8 QUANTITY='TEMPERATURE', ID='TEMP_4' /
&DEVC XYZ=1.2,0.8,2.4 QUANTITY='TEMPERATURE', ID='TEMP_5' /
&DEVC XYZ=1.2,0.8,3.0 QUANTITY='TEMPERATURE', ID='TEMP_6' /

&DEVC XYZ=1.2,0.8,0.0 QUANTITY='VISIBILITY', ID='VIS_1' /
&DEVC XYZ=1.2,0.8,0.6 QUANTITY='VISIBILITY', ID='VIS_2' /
&DEVC XYZ=1.2,0.8,1.2 QUANTITY='VISIBILITY', ID='VIS_3' /
&DEVC XYZ=1.2,0.8,1.8 QUANTITY='VISIBILITY', ID='VIS_4' /
&DEVC XYZ=1.2,0.8,2.4 QUANTITY='VISIBILITY', ID='VIS_5' /
&DEVC XYZ=1.2,0.8,3.0 QUANTITY='VISIBILITY', ID='VIS_6' /

&TAIL /
Loading
Loading