From 79bed0056e5d74106093f09bf7d329b87a4edcbf Mon Sep 17 00:00:00 2001 From: mcgratta Date: Mon, 20 Nov 2023 13:46:07 -0500 Subject: [PATCH] FDS Verification: Add ht3d test case --- .../FDS_Verification_Guide.tex | 12 ++++++ .../FDS_verification_dataplot_inputs.csv | 1 + Verification/FDS_Cases.sh | 1 + .../ht3d_energy_conservation_8.fds | 37 +++++++++++++++++++ 4 files changed, 51 insertions(+) create mode 100644 Verification/Heat_Transfer/ht3d_energy_conservation_8.fds diff --git a/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex b/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex index 126fcb2d8ef..851af9faf9b 100644 --- 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 index 832322bb2e8..10f80460c14 100644 --- 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 index 5877bcdb4f9..90677462155 100755 --- 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 new file mode 100644 index 00000000000..65035e2bfa1 --- /dev/null +++ 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 /