Merge pull request #12242 from mcgratta/master

FDS Verification: Add ht3d test case
firemodels · Nov 20, 2023 · aa0207a · aa0207a
2 parents b41f54f + 79bed00
commit aa0207a
Show file tree

Hide file tree

Showing 4 changed files with 51 additions and 0 deletions.
diff --git a/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex b/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex
@@ -4465,6 +4465,18 @@ \subsection{Energy Conservation in a 3-D Solid (\texorpdfstring{\textct{ht3d\_en
 \label{fig:ht3d_energy_conservation_4}
 \end{figure}
 
+\FloatBarrier
+
+Figure~\ref{fig:ht3d_energy_conservation_8} compares the integrated net heat flux versus the internal enthalpy in a case where a vertical steel column is exposed to a radiant heater. The column is modeled in FDS using thin obstructions to represent a 10~mm thick web and two 18~mm thick flanges. One flange is subjected to the heat, which then passes through the web. 
+
+\begin{figure}[ht]
+\centering
+\includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_energy_conservation_8}
+\caption[Additional \textct{ht3d\_energy\_conservation} test case 8]{Comparison of the integrated net heat flux versus the internal enthalpy for a steel column.}
+\label{fig:ht3d_energy_conservation_8}
+\end{figure}
+
+
 \FloatBarrier
 
 \subsection{Heat Transfer to a Steel Beam (\texorpdfstring{\textct{ht3d\_beam\_heating}}{ht3d\_beam\_heating})}

diff --git a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
@@ -269,6 +269,7 @@ d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_4_git.txt,Heat
 d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_5_git.txt,Heat_Transfer/ht3d_energy_conservation_5_devc.csv,2,3,Time,H1,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_5_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation\_5),Time (s),Enthalpy (kJ),0,100,1,0,4,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_5,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX
 d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_6_git.txt,Heat_Transfer/ht3d_energy_conservation_6_devc.csv,2,3,Time,H1,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_6_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation\_6),Time (s),Enthalpy (kJ),0,100,1,0,4,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_6,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX
 d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_7_git.txt,Heat_Transfer/ht3d_energy_conservation_7_devc.csv,2,3,Time,H1,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_7_devc.csv,2,3,Time,Q_net_1+Q_net_2,Integrated Heat Flux,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation\_7),Time (s),Enthalpy (kJ),0,300,1,0,30,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_7,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX
+d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_8_git.txt,Heat_Transfer/ht3d_energy_conservation_8_devc.csv,2,3,Time,H1,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_8_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation\_8),Time (s),Enthalpy (kJ),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_8,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX
 d,ht3d_ibeam,Heat_Transfer/ht3d_ibeam_git.txt,Heat_Transfer/ht3d_ibeam_FEM_results.csv,2,3,Time,Ts_x195_40|Ts_x145_30|Ts_x095_20|Ts_x025_40|Ts_x195_1|Ts_x025_1,FEM 1|FEM 2|FEM 3|FEM 4|FEM 5|FEM 6,ro|k^|bd|gsq|mv|c>,0,100000,,0,100000,-1.00E+09,1.00E+09,20,Heat_Transfer/ht3d_ibeam_devc.csv,2,3,Time,TS_x195-40|TS_x145-30|TS_x095-20|TS_x025-40|TS_x195-01|TS_x025-01,FDS 1|FDS 2|FDS 3|FDS 4|FDS 5|FDS 6,r-|k-|b-|g-|m-|c-,0,100000,,0,100000,-1.00E+09,1.00E+09,20,HT3D I-beam Surface Temperature (ht3d\_ibeam),Time (s),Temperature (°C),0,3600,1,0,1000,1,no,0.05 0.90,EastOutside,,1.2,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_ibeam_TS,Relative Error,end,8.00E-02,Heat Transfer,r^,r,TeX
 d,ht3d_mass_conservation,Heat_Transfer/ht3d_mass_conservation_git.txt,Heat_Transfer/ht3d_mass_conservation.csv,1,2,Time,Mass,Exact (Mass),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_mass_conservation_mass.csv,2,3,Time,WOOD MOISTURE,FDS (WOOD MOISTURE),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Mass Balance  (ht3d\_mass\_conservation),Time (s),Mass (kg),0,180,1,0,0.3,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_mass_conservation,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX
 d,ht3d_mass_conservation,Heat_Transfer/ht3d_mass_conservation_2_git.txt,Heat_Transfer/ht3d_mass_conservation.csv,1,2,Time,Mass,Exact (Mass),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_mass_conservation_2_mass.csv,2,3,Time,WOOD MOISTURE,FDS (WOOD MOISTURE),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Mass Balance  (ht3d\_mass\_conservation\_2),Time (s),Mass (kg),0,180,1,0,0.3,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_mass_conservation_2,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX

diff --git a/Verification/FDS_Cases.sh b/Verification/FDS_Cases.sh
@@ -256,6 +256,7 @@ $QFDS -p 8 -d Heat_Transfer ht3d_energy_conservation_4.fds
 $QFDS -p 8 -d Heat_Transfer ht3d_energy_conservation_5.fds
 $QFDS -p 8 -d Heat_Transfer ht3d_energy_conservation_6.fds
 $QFDS -p 8 -d Heat_Transfer ht3d_energy_conservation_7.fds
+$QFDS -p 36 -d Heat_Transfer ht3d_energy_conservation_8.fds
 $QFDS -d Heat_Transfer ht3d_ibeam.fds
 $QFDS -d Heat_Transfer ht3d_mass_conservation.fds
 $QFDS -d Heat_Transfer ht3d_mass_conservation_2.fds

diff --git a/Verification/Heat_Transfer/ht3d_energy_conservation_8.fds b/Verification/Heat_Transfer/ht3d_energy_conservation_8.fds
@@ -0,0 +1,37 @@
+&HEAD CHID='ht3d_energy_conservation_8', TITLE='Heated Column' /
+
+&MESH IJK=5,12,30, XB=-2.30,-2.05,-0.90,-0.30,0.00,1.50, MULT_ID='mesh' /
+&MULT ID='mesh', DX=0.25, DY=0.60, DZ=1.50, I_UPPER=5, J_UPPER=2, K_UPPER=1 /
+
+&TIME T_END=900. /                                   
+
+&MISC NEIGHBOR_SEPARATION_DISTANCE=3 / 
+
+&SURF ID='HOT', TMP_FRONT=500, COLOR='RED', HEAT_TRANSFER_COEFFICIENT=0 /
+&OBST XB=-1.40, -1.20, -0.10, 0.10, 00.40, 01.00, COLOR='RED', SURF_ID6='HOT','INERT','INERT','INERT','INERT','INERT' /
+
+&OBST XB=-2.000, -1.982, -0.150, 0.150, 00.05, 02.50, SURF_ID='Column' / 
+&OBST XB=-2.000, -1.500, -0.005, 0.005, 00.05, 02.50, SURF_ID='Column' / 
+&OBST XB=-1.518, -1.500, -0.150, 0.150, 00.05, 02.50, SURF_ID='Column' / 
+
+&SURF ID='Column', MATL_ID='STEEL', HT3D=T, CELL_SIZE=0.025, COLOR='BLACK' /
+&MATL ID='STEEL', SPECIFIC_HEAT=0.45, CONDUCTIVITY=50, DENSITY=7850, EMISSIVITY=0.7 /
+
+&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' /
+
+&BNDF QUANTITY='WALL_TEMPERATURE', CELL_CENTERED=T /
+
+&DUMP DT_DEVC=20 /
+
+'WALL ENTHALPY' is the energy (kJ) of the volume of solid bounded by the surface cell. The CONVERSION_FACTOR is intended to
+cancel out the cell area 0.05 m x 0.05 m
+
+&DEVC XB=-2.01,-1.49,-0.16,0.16,0.04,2.51, QUANTITY='WALL ENTHALPY', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='H1', IOR=-1, TIME_AVERAGED=F, RELATIVE=T, CONVERSION_FACTOR=400, SURF_ID='Column' /
+
+&DEVC XB=-2.01,-1.49,-0.16,0.16,0.04,2.51, QUANTITY='NET HEAT FLUX', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='Q_net', SURF_ID='Column' /
+
+&TAIL /