From 500884bfcdbc3ba2d019fd973a5e1d9ba4a0628e Mon Sep 17 00:00:00 2001 From: mcgratta Date: Wed, 1 Nov 2023 13:40:42 -0400 Subject: [PATCH] FDS Source: Allow different BACKINGs for HT3d --- Manuals/FDS_User_Guide/FDS_User_Guide.tex | 8 ++-- .../FDS_Verification_Guide.tex | 5 ++- Source/init.f90 | 3 +- .../FDS_verification_dataplot_inputs.csv | 17 ++++---- Verification/FDS_Cases.sh | 1 + .../ht3d_energy_conservation_7.fds | 39 +++++++++++++++++++ 6 files changed, 58 insertions(+), 15 deletions(-) create mode 100644 Verification/Heat_Transfer/ht3d_energy_conservation_7.fds diff --git a/Manuals/FDS_User_Guide/FDS_User_Guide.tex b/Manuals/FDS_User_Guide/FDS_User_Guide.tex index 5c23e9b1247..0fd9079cf69 100644 --- a/Manuals/FDS_User_Guide/FDS_User_Guide.tex +++ b/Manuals/FDS_User_Guide/FDS_User_Guide.tex @@ -2443,7 +2443,7 @@ \subsection{3-D Heat Conduction (Beta)} \subsubsection{Limitations} \begin{enumerate} -\item {\ct HT3D} cannot be applied to an exterior boundary; it must be applied to an {\ct OBST} that is completely within the computational domain. +\item {\ct HT3D} cannot be applied to an exterior boundary; it must be applied to an {\ct OBST} for which at least one face in each coordinate direction is exposed. \item Avoid contact between 3-D and 1-D solids. If two sides of a 3-D solid touch 1-D solids, there will be no lateral heat conduction computed in that particular direction. \item If your 3-D obstruction extends beyond meshes that abut, add the parameter {\ct NEIGHBOR\_SEPARATION\_DISTANCE} to the {\ct MISC} line. Any mesh within this distance of another mesh will share geometry information for use in the 3-D heat conduction calculation. \item By default, the interior nodes are clustered near the surface and stretched out deeper within the solid. If you want to maintain uniform spacing, set {\ct CELL\_SIZE} on the {\ct SURF} line to indicate the desired interior node spacing. The {\ct CELL\_SIZE} is typically the same as the gas phase cells. @@ -2478,7 +2478,7 @@ \subsubsection{Example: Steel Assembly} The quantity {\ct 'WALL ENTHALPY'} is the total energy within the volume carved out by a wall cell from front to back. This is why the {\ct IOR} is included because the integration need only be carried out for upward facing wall cells. Also, because the unit associated with this quantity is kJ rather than kJ/m$^2$, the {\ct CONVERSION\_FACTOR} is needed to cancel out the area of each wall cell (0.01~m by 0.01~m) after the summation. The parameter {\ct RELATIVE=T} means that only the change in enthalpy is desired. . -\subsubsection{Example: Uncoated Structural Steel} +\subsubsection{Example: Unprotected Structural Steel} Structural steel members like I-beams are difficult to model for two reasons. First, the underlying grid must be reasonably well-resolved to capture the cross-sectional shape, and second, the lateral heat transfer can be important. At best, I-beams are typically modeled as a collection of thin (i.e. zero cell thick) obstructions representing the web and flanges. \begin{lstlisting} @@ -2502,7 +2502,7 @@ \subsubsection{Example: Uncoated Structural Steel} \end{figure} -\subsubsection{Example: Coated Structural Steel} +\subsubsection{Example: Insulated Structural Steel} \label{info:normal_direction_only} An insulated steel beam can be modeled as a hybrid of 1-D and 3-D objects. The heat conducted through the insulation varies mainly in the direction normal to the surface, while the heat conducted along the steel beam varies mainly in the lateral direction. The following lines provide an example where a steel plate is coated with insulation. The steel plate is entered as a relatively thin obstruction, and the insulation is entered as a slightly thicker obstruction with a {\ct SURF} line that indicates that heat is to be conducted in the {\ct NORMAL\_DIRECTION\_ONLY} even though it is declared to be {\ct HT3D}. This latter parameter ensures that the thin steel obstruction and the insulation obstruction are included in the calculation of solid overlap volumes. Note that all dimensions are exact to ensure that the overlap volumes are computed properly. @@ -2930,7 +2930,7 @@ \subsection{The Heat of Reaction} \subsection{Liquid Fuels} \label{info:liquid_fuels} -The evaporation rate of a liquid fuel is analogous to the convective heating rate in that the evaporation rate is a function of a mass transfer coefficient, $h_{\rm m}$, much like thermal convection is a function of the heat transfer coefficient, $h$, discussed in Section~\ref{info:convection}\footnote{As with the convective heat transfer coefficient, there is an option to specify a fixed {\ct MASS\_TRANSFER\_COEFFICIENT} (m/s) on the {\ct SURF} line that describes a liquid pool.}. The FDS Technical Reference Guide~\cite{FDS_Tech_Guide} provides further details on how the evaporation rate is computed. +The evaporation rate of a liquid fuel is analogous to the convective heating rate in that the evaporation rate is a function of a mass transfer coefficient, $h_{\rm m}$, much like thermal convection is a function of the heat transfer coefficient, $h$, discussed in Section~\ref{info:convection}\footnote{As with the convective heat transfer coefficient, there is an option to specify a fixed {\ct MASS\_TRANSFER\_COEFFICIENT} (m/s) on the {\ct SURF} line that describes a liquid pool.}. The FDS Technical Reference Guide~\cite{FDS_Tech_Guide} provides further details on how the evaporation rate is computed. The properties of a liquid fuel are given on the {\ct MATL} line: \begin{lstlisting} diff --git a/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex b/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex index a38e4ffa8d2..126fcb2d8ef 100644 --- a/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex +++ b/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex @@ -4458,9 +4458,10 @@ \subsection{Energy Conservation in a 3-D Solid (\texorpdfstring{\textct{ht3d\_en \begin{tabular*}{\textwidth}{l@{\extracolsep{\fill}}r} \includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_energy_conservation_4} & \includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_energy_conservation_5} \\ -\multicolumn{2}{c}{\includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_energy_conservation_6}} +\includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_energy_conservation_6} & +\includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_energy_conservation_7} \end{tabular*} -\caption[Additional \textct{ht3d\_energy\_conservation} test cases, 4, 5, and 6]{Comparison of the integrated net heat flux versus the internal enthalpy for a homogenous solid block of insulation material (upper left), a block of steel with a single layer of insulation (upper right), and a block with multiple layers of multi-component insulation (bottom).} +\caption[Additional \textct{ht3d\_energy\_conservation} test cases, 4, 5, 6, and 7]{Comparison of the integrated net heat flux versus the internal enthalpy for a homogenous solid block of insulation material (upper left), a block of steel with a single layer of insulation (upper right), a block with multiple layers of multi-component insulation (bottom left), a block of steel that abuts an open boundary (bottom right).} \label{fig:ht3d_energy_conservation_4} \end{figure} diff --git a/Source/init.f90 b/Source/init.f90 index ffe6d80ddbd..0a49821d286 100644 --- a/Source/init.f90 +++ b/Source/init.f90 @@ -3620,7 +3620,7 @@ SUBROUTINE FIND_WALL_BACK_INDEX(NM,IW) WC => M%WALL(IW) SF => SURFACE(WC%SURF_INDEX) -IF (SF%THERMAL_BC_INDEX/=THERMALLY_THICK .OR. SF%BACKING/=EXPOSED) RETURN +IF (SF%THERMAL_BC_INDEX/=THERMALLY_THICK) RETURN BC => M%BOUNDARY_COORD(WC%BC_INDEX) ICG = M%CELL_INDEX(BC%IIG,BC%JJG,BC%KKG) IF (M%CELL(ICG)%SOLID) RETURN @@ -3718,6 +3718,7 @@ SUBROUTINE FIND_WALL_BACK_INDEX(NM,IW) ! Determine if the back face is found IF ((.NOT.OM%CELL(IC)%SOLID .AND. OM%CELL(IC)%WALL_INDEX(IOR)>0) .OR. NOM==0) THEN ! the back wall face is found + IF (NOM>0 .AND. SF%BACKING/=EXPOSED) RETURN ! No need to assign back cell information for anything but exposed backing ONE_D%BACK_INDEX = OM%CELL(IC)%WALL_INDEX(IOR) ONE_D%BACK_MESH = NOM ONE_D%BACK_SURF = OM%CELL(IC)%SURF_INDEX(IOR) diff --git a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv index 464236461a8..5b925ca06d1 100644 --- a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv +++ b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv @@ -259,14 +259,15 @@ d,ht1d_pile,Heat_Transfer/ht1d_pile_git.txt,Heat_Transfer/ht1d_pile.csv,2,3,Dept d,ht3d_pile,Heat_Transfer/ht3d_pile_git.txt,Heat_Transfer/ht1d_pile.csv,2,3,Depth,Temperature,Exact,k-,0,100000,,0.05,0.13,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_pile_prof_1.csv,2,3,Depth,Left,FDS,ko,0,100000,,0.07,0.11,-1.00E+09,1.00E+09,20,Heat Transfer (ht3d\_pile),Depth (m),Temperature (°C),0,0.18,1,0,600,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_pile,Relative Error,slope,1.00E-02,Heat Transfer,r^,r,TeX d,ht3d_beam_heating,Heat_Transfer/ht3d_beam_heating_1_git.txt,Heat_Transfer/ht3d_beam_heating_1.csv,1,2,Time,T,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,20,Heat_Transfer/ht3d_beam_heating_1_devc.csv,2,3,Time,wall temp,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,20,Heat Transfer (ht3d\_beam\_heating\_1),Time (s),Beam Temperature (°C),0,5,60,20,24,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_beam_heating_1,Relative Error,end,6.00E-02,Heat Transfer,r^,r,TeX d,ht3d_beam_heating,Heat_Transfer/ht3d_beam_heating_2_git.txt,Heat_Transfer/ht3d_beam_heating_2.csv,1,2,Time,T,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,20,Heat_Transfer/ht3d_beam_heating_2_devc.csv,2,3,Time,wall temp,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,20,Heat Transfer (ht3d\_beam\_heating\_2),Time (s),Beam Temperature (°C),0,60,60,20,24,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_beam_heating_2,Relative Error,end,5.00E-02,Heat Transfer,r^,r,TeX -d,ht3d_demo,Heat_Transfer/ht3d_demo_git.txt,Heat_Transfer/ht3d_demo_devc.csv,2,3,Time,H,Enthalpy,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_demo_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_demo),Time (s),Enthalpy (kJ),0,100,1,0,10,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/ht3d_demo,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX -d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_git.txt,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,E3D,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation),Time (s),Enthalpy (kJ),0,10,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX -d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_git.txt,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,E3D,Enthalpy,k-o,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,r-o,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation),Time (s),Enthalpy (kJ),9.5,10,1,92,97,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_fine,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX -d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_2_git.txt,Heat_Transfer/ht3d_energy_conservation_2_devc.csv,2,3,Time,E3D,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_2_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\_2),Time (s),Enthalpy (kJ),0,10,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_2,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX -d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_3_git.txt,Heat_Transfer/ht3d_energy_conservation_3_devc.csv,2,3,Time,E3D,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_3_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\_3),Time (s),Enthalpy (kJ),0,10,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_3,Relative Error,end,3.00E-02,Heat Transfer,r^,r,TeX -d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_4_git.txt,Heat_Transfer/ht3d_energy_conservation_4_devc.csv,2,3,Time,H1,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_4_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\_4),Time (s),Enthalpy (kJ),0,100,1,0,3,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_4,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX -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,3,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,3,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_demo,Heat_Transfer/ht3d_demo_git.txt,Heat_Transfer/ht3d_demo_devc.csv,2,3,Time,H,Enthalpy,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_demo_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_demo),Time (s),Enthalpy (kJ),0,100,1,0,10,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/ht3d_demo,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX +d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_git.txt,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,E3D,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation),Time (s),Enthalpy (kJ),0,10,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX +d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_git.txt,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,E3D,Enthalpy,k-o,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_devc.csv,2,3,Time,Q_net,Integrated Heat Flux,r-o,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Energy Balance (ht3d\_energy\_conservation),Time (s),Enthalpy (kJ),9.5,10,1,92,97,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_fine,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX +d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_2_git.txt,Heat_Transfer/ht3d_energy_conservation_2_devc.csv,2,3,Time,E3D,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_2_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\_2),Time (s),Enthalpy (kJ),0,10,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_2,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX +d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_3_git.txt,Heat_Transfer/ht3d_energy_conservation_3_devc.csv,2,3,Time,E3D,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_3_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\_3),Time (s),Enthalpy (kJ),0,10,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_energy_conservation_3,Relative Error,end,3.00E-02,Heat Transfer,r^,r,TeX +d,ht3d_energy_conservation,Heat_Transfer/ht3d_energy_conservation_4_git.txt,Heat_Transfer/ht3d_energy_conservation_4_devc.csv,2,3,Time,H1,Enthalpy,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/ht3d_energy_conservation_4_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\_4),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_4,Relative Error,end,1.00E-02,Heat Transfer,r^,r,TeX +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_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 fff6e65224f..7381bf4d360 100755 --- a/Verification/FDS_Cases.sh +++ b/Verification/FDS_Cases.sh @@ -254,6 +254,7 @@ $QFDS -d Heat_Transfer ht3d_energy_conservation_3.fds $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 -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_7.fds b/Verification/Heat_Transfer/ht3d_energy_conservation_7.fds new file mode 100644 index 00000000000..576f14bda80 --- /dev/null +++ b/Verification/Heat_Transfer/ht3d_energy_conservation_7.fds @@ -0,0 +1,39 @@ +&HEAD CHID='ht3d_energy_conservation_7' / + +&MESH IJK=25,25,8, XB=-0.25,0.00,-0.25,0.00,-0.05,0.03, MULT_ID='mesh' / +&MULT ID='mesh', DX=0.25, DY=0.25, DZ=0.08, I_UPPER=1, J_UPPER=1, K_UPPER=1 / + +&TIME T_END=300 / + +&VENT DB='XMIN', SURF_ID='OPEN' / +&VENT DB='XMAX', SURF_ID='OPEN' / +&VENT DB='YMIN', SURF_ID='OPEN' / +&VENT DB='YMAX', SURF_ID='OPEN' / +&VENT DB='ZMIN', SURF_ID='OPEN' / +&VENT DB='ZMAX', SURF_ID='OPEN' / + +&OBST XB=-0.20, 0.20,-0.20, 0.20,-0.05, 0.05, SURF_IDS='SLAB TOP','SLAB','SLAB', MATL_ID='STEEL' / +&OBST XB=-0.01, 0.01,-0.01, 0.01, 0.06, 0.07, SURF_ID='HOT' / + +&SURF ID='HOT', TMP_FRONT=1000, COLOR='RED' / + +&SURF ID='SLAB', HT3D=T, COLOR='BEIGE', MATL_ID='STEEL', CELL_SIZE=0.01 / +&SURF ID='SLAB TOP', HT3D=T, COLOR='BEIGE', MATL_ID='STEEL', CELL_SIZE=0.01, BACKING='INSULATED' / + +&MATL ID='STEEL', DENSITY=7500, SPECIFIC_HEAT=0.5, CONDUCTIVITY=50. / + +&BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=T / + +&SLCF PBY=0.001, QUANTITY='TEMPERATURE', CELL_CENTERED=T / + +&DUMP DT_DEVC=12. / + +'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.01 m x 0.01 m + +&DEVC XB=-0.25,0.25,-0.25,0.25,-0.05,0.11, QUANTITY='WALL ENTHALPY', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='H1', IOR=-1, TIME_AVERAGED=F, RELATIVE=T, CONVERSION_FACTOR=10000, SURF_ID='SLAB' / + +&DEVC XB=-0.25,0.25,-0.25,0.25,-0.05,0.11, QUANTITY='NET HEAT FLUX', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='Q_net_1', SURF_ID='SLAB' / +&DEVC XB=-0.25,0.25,-0.25,0.25,-0.05,0.11, QUANTITY='NET HEAT FLUX', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='Q_net_2', SURF_ID='SLAB TOP' / + +&TAIL /