Merge pull request #14636 from mcgratta/master

FDS Source: Add Q_RAD_IN for fixed gas temp; add new verification cases

Author: mcgratta

Manuals/Bibliography/FDS_general.bib

@@ -5820,6 +5820,13 @@ @INPROCEEDINGS{Semmes:IAFSS11
   note         = {Accepted}
 }
 
+@MANUAL{SFPE_S.02,
+  title        = {{SFPE Standard S.02 -- Calculation Methods to Predict the Thermal Performance of Structures \& Fire Resistive Assemblies}},
+  organization = {Society of Fire Protection Engineers},
+  address      = {Gaithersburg, Maryland},
+  year         = {2007}
+}
+
 @ARTICLE{Sarwar:2013,
     author={M. Sarwar and K. A. M. Moinuddin and G. R. Thorpe},
     title={Large Eddy Simulation of Flow over a Backward Facing Step using {Fire Dynamics Simulator (FDS)}},

Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex

@@ -5016,6 +5016,36 @@ \subsection{3-D Heat Conduction in a Layered Solid (\texorpdfstring{\ct{ht3d\_pi
 \label{fig:ht3d_pile}
 \end{figure}
 
+\clearpage
+
+\section{SFPE Verification Cases}
+
+The Society of Fire Protection Engineers (SFPE) has developed a standard entitled {\em S.02 -- Calculation Methods to Predict the Thermal Performance of Structures \& Fire Resistive Assemblies}~\cite{SFPE_S.02} that contains an appendix with verification cases to benchmark basic heat transfer calculations. This section contains several of these cases.
+
+\subsection{Case 6: 2-D Heat Transfer with Cooling by Convection}
+\label{SFPE_Case_6}
+
+A 2~m by 2~m square column ($k=1$~W/(m~K), $\rho=1$~kg/m$^3$, $c=0.001$~kJ/(kg~K), $\epsilon=0$) with an initial temperature of 1000~°C cools via convection only. Assuming that $h=1$~W/(m$^2$~K) and the surrounding air temperature is 0~°C, calculate the temperature at the center of the column as a function of time.
+
+This is an unusual set of parameters, but nevertheless this case is used to test a multi-dimensional heat transfer solver. The results are shown in Fig.~\ref{fig:SFPE_Case_6_7}.
+
+\subsection{Case 7: 2-D Heat Transfer by Convection and Radiation}
+\label{SFPE_Case_7}
+
+A 0.2 m by 0.2 m square column ($k=1$~W/(m~K), $\rho=2400$~kg/m$^3$, $c=1$~kJ/(kg~K), $\epsilon=0.8$) is heated according to the ISO~834 time-temperature curve
+\be
+   T(t) = T_\infty + 345 \, \ln \left( 8t/60 + 1 \right)
+\ee
+where the time, $t$, is in seconds. Assuming that $h=10$~W/(m$^2$~K) and that the initial temperature is $T_\infty=273$~K, calculate the temperature at the column center, corner and middle side surface as a function of time. The results are shown in Fig.~\ref{fig:SFPE_Case_6_7}.
+
+\begin{figure}[ht]
+\includegraphics[height=2.2in]{SCRIPT_FIGURES/SFPE_Case_6}
+\includegraphics[height=2.2in]{SCRIPT_FIGURES/SFPE_Case_7}
+\caption[The SFPE heat transfer verification cases 6 and 7]{Comparisons of FDS 2-D heat transfer calculations with solutions generated by a finite-element solver.}
+\label{fig:SFPE_Case_6_7}
+\end{figure}
+
+
 
 
 \chapter{Pyrolysis}

Source/wall.f90

@@ -385,6 +385,7 @@ SUBROUTINE NEAR_SURFACE_GAS_VARIABLES(T,SF,BC,B1,LP,TW,WALL_INDEX,PARTICLE_INDEX
 
    IF (SF%TMP_GAS_FRONT > 0._EB) THEN
       B1%TMP_G = TMPA + EVALUATE_RAMP(T-T_BEGIN,SF%RAMP(TIME_TGF)%INDEX)*(SF%TMP_GAS_FRONT-TMPA)
+      B1%Q_RAD_IN = B1%EMISSIVITY*SIGMA*B1%TMP_G**4
    ELSE
       B1%TMP_G =  TMP(BC%IIG,BC%JJG,BC%KKG)
    ENDIF

Utilities/Matlab/FDS_verification_dataplot_inputs.csv

@@ -612,6 +612,8 @@ d,realizable_mass_fractions,Flowfields/realizable_mass_fractions_git.txt,Flowfie
 d,rms_example,Controls/rms_example_git.txt,Controls/rms_example.csv,1,2,Time,mean|sd,Expected Mean|Expected RMS,k--|r--,0,100000,,5,20,-1.00E+09,1.00E+09,0,Controls/rms_example_devc.csv,2,3,Time,u|u_rms,FDS (u)|FDS (u\_rms),k-|r-,0,100000,,5,20,-1.00E+09,1.00E+09,0,rms\_example,Time (s),Velocity (m/s),0,20,1,0,2.5,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/rms_example,Relative Error,mean,4.00E-02,Controls,mx,m,TeX
 d,screen_drag_1,Sprinklers_and_Sprays/screen_drag_1_git.txt,Sprinklers_and_Sprays/screen_drag_1.csv,1,1,Time,Exact 0.1|Exact 0.4|Exact 0.8,Exact 0.1|Exact 0.4|Exact 0.8,ko|ro|bo,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/screen_drag_1_devc.csv,2,3,Time,FDS 0.1|FDS 0.4|FDS 0.8,FDS 0.1|FDS 0.4|FDS 0.8,k-|r-|b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Screen Pressure Drop,Time (s),Pressure Drop (Pa),0,4,1,0,150,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/screen_drag_1,Relative Error,end,0.03,Sprinklers and Sprays,bs,b,TeX
 d,screen_drag_2,Sprinklers_and_Sprays/screen_drag_2_git.txt,Sprinklers_and_Sprays/screen_drag_2.csv,1,1,Time,Exact 0.1|Exact 0.4|Exact 0.8,Exact 0.1|Exact 0.4|Exact 0.8,ko|ro|bo,0,100000,,6,10,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/screen_drag_2_devc.csv,2,3,Time,FDS 0.1|FDS 0.4|FDS 0.8,FDS 0.1|FDS 0.4|FDS 0.8,k-|r-|b-,0,100000,,6,10,-1.00E+09,1.00E+09,0,Screen Pressure Drop,Time (s),Pressure Drop (Pa),0,10,1,0,1500,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/screen_drag_2,Relative Error,mean,0.08,Sprinklers and Sprays,bs,b,TeX
+d,SFPE_Case_6,Heat_Transfer/SFPE_Case_6_git.txt,Heat_Transfer/SFPE_Case_6.csv,1,2,Time,Temperature,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/SFPE_Case_6_devc.csv,2,3,Time,T,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,SFPE Case 6,Time (s),Temperature (°C),0,1,1,0,1200,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/SFPE_Case_6,Relative Error,end,0.01,Heat Transfer,bs,b,TeX
+d,SFPE_Case_7,Heat_Transfer/SFPE_Case_7_git.txt,Heat_Transfer/SFPE_Case_7.csv,1,2,Time,Center|Side|Corner,Exact Center|Exact Side|Exact Corner,ko|ro|go,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/SFPE_Case_7_devc.csv,2,3,Time,Center|Side|Corner,FDS Center|FDS Side|FDS Corner,k-|r-|g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,SFPE Case 7,Time (min),Temperature (°C),0,180,60,0,1200,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/SFPE_Case_7,Relative Error,end,0.01,Heat Transfer,bs,b,TeX
 d,shrink_swell,Pyrolysis/shrink_swell_git.txt,Pyrolysis/shrink_swell.csv,1,2,Time,M_1 (kg/m2),Expected,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/shrink_swell_devc.csv,2,3,Time,M_1|M_2|M_3,Shrink 1|Shrink 2|Shrink 3,k--|b-|r-.,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Shrink Mass Conservation (shrink\_swell),Time (s),Surface density (kg/m²),0,15,1,0,1,1,no,0.05 0.90,South,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/shrink_mass_conservation,Relative Error,end_1_1,0.01,Pyrolysis,rx,r,TeX
 d,shrink_swell,Pyrolysis/shrink_swell_git.txt,Pyrolysis/shrink_swell.csv,1,2,Time,M_4 (kg/m2),Expected,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/shrink_swell_devc.csv,2,3,Time,M_4|M_5|M_6,Swell 1|Swell 2|Swell 3,k--|b-|r-.,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Swell Mass Conservation (shrink\_swell),Time (s),Surface density (kg/m²),0,15,1,0,2,1,no,0.05 0.90,South,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/swell_mass_conservation,Relative Error,end_1_1,0.01,Pyrolysis,rx,r,TeX
 d,shrink_swell,Pyrolysis/shrink_swell_git.txt,Pyrolysis/shrink_swell.csv,1,2,Time,L_1 (m)|L_2 (m)|L_3 (m),Expected 1|Expected 2|Expected 3,ko|bd|rv,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/shrink_swell_devc.csv,2,3,Time,L_1|L_2|L_3,Shrink 1|Shrink 2|Shrink 3,k--|b-|r-.,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Shrinking Thickness (shrink\_swell),Time (s),Thickness (m),0,15,1,0,0.0015,1,no,0.05 0.90,South,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/shrink_thickness,Relative Error,end_1_1,0.01,Pyrolysis,rx,r,TeX
@@ -805,4 +807,4 @@ f,anca-couce,Pyrolysis/anca-couce-fig2_2p5K_git.txt,Pyrolysis/Anca-Couce-fig1.cs
 f,anca-couce,Pyrolysis/anca-couce-fig2_2p5K_git.txt,Pyrolysis/Anca-Couce-fig1.csv,1,2,2.5K/min TMP,2.5K/min MLR,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/anca-couce-fig2_2p5K_cat_tga.csv,2,3,Temp,LIGNIN MLR,LIGNIN,k:,0,100000,,0,100000,-1.00E+09,1.00E+09,0,0% O_2 (3C Wood Pyrolysis),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,3.20E-03,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/anca_couce_fig2_2p5K,N/A,max,0.1,pine wood TGA,kd,k,TeX
 s,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
 g,sphere_leak,Complex_Geometry/sphere_leak_git.txt,Complex_Geometry/sphere_leak.csv,1,2,Time,Pressure,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Complex_Geometry/sphere_leak_devc.csv,2,3,Time,Pressure,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Rise (sphere\_leak),Time (s),Pressure (Pa),0,100,1,0,5000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/sphere_leak,Relative Error,max,0.05,Pressure Effects,k+,k,TeX
-d,sphere_radiate,Complex_Geometry/sphere_radiate_git.txt,Complex_Geometry/sphere_radiate.csv,1,2,Time,HF,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Complex_Geometry/sphere_radiate_devc.csv,2,3,Time,HF1,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat Flux (sphere\_radiate),Time (s),Heat Flux (kW/m²),0,0.01,1,0,8,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/sphere_radiate,Relative Error,max,0.07,Radiation,bs,b,TeX
+d,sphere_radiate,Complex_Geometry/sphere_radiate_git.txt,Complex_Geometry/sphere_radiate.csv,1,2,Time,HF,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Complex_Geometry/sphere_radiate_devc.csv,2,3,Time,HF1,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat Flux (sphere\_radiate),Time (s),Heat Flux (kW/m²),0,0.01,1,0,8,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/sphere_radiate,Relative Error,max,0.07,Radiation,bs,b,TeX

Verification/FDS_Cases.sh

@@ -324,6 +324,8 @@ $QFDS -p 4 -d Heat_Transfer ht3d_sphere_48.fds
 $QFDS -p 64 -d Heat_Transfer ht3d_sphere_96.fds
 $QFDS -p 4 -d Heat_Transfer back_wall_test.fds
 $QFDS -p 3 -d Heat_Transfer back_wall_test_2.fds
+$QFDS -d Heat_Transfer SFPE_Case_6.fds
+$QFDS -d Heat_Transfer SFPE_Case_7.fds
 $QFDS -d Heat_Transfer thermocouple_time_constant.fds
 
 $QFDS -d HVAC ashrae7_fixed_flow.fds

Verification/Heat_Transfer/SFPE_Case_6.csv

@@ -0,0 +1,9 @@
+Time,Temperature
+0.0,1000
+0.1,986.4
+0.2,903.8
+0.4,690.2
+0.6,514.7
+0.8,382.7
+1.0,284.5
+

Verification/Heat_Transfer/SFPE_Case_6.fds

@@ -0,0 +1,37 @@
+&HEAD CHID='SFPE_Case_6', TITLE='Verification Case 6, SFPE Guide to Thermal Performance of Resistive Assemblies'/
+
+&MESH IJK=24,24,3, XB=-1.2,1.2,-1.2,1.2,-0.1,0.2 /
+
+&TIME T_END=1, DT=0.01, WALL_INCREMENT=1 /
+
+&MISC SOLID_PHASE_ONLY=T, TMPA=1000. /
+
+&RADI RADIATION=F /
+
+&OBST XB=-1.0, 1.0,-1.0, 1.0,0.0,0.1, SURF_IDS='ADIABATIC','STUFF','ADIABATIC', MATL_ID='STUFF' /
+
+&SURF ID='STUFF', COLOR='GRAY', HEAT_TRANSFER_COEFFICIENT=1.0, HT3D=T, CELL_SIZE=0.1, TMP_GAS_FRONT=0. /
+&SURF ID='ADIABATIC', HEAT_TRANSFER_COEFFICIENT=0.0, HT3D=T, CELL_SIZE=0.1 /
+
+&MATL ID            = 'STUFF'
+      SPECIFIC_HEAT = 0.001
+      CONDUCTIVITY  = 1.0
+      DENSITY       = 1.0 /
+
+&BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=T/
+&BNDF QUANTITY='NET HEAT FLUX', CELL_CENTERED=T/
+
+&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' /
+
+&DUMP DT_DEVC=0.03 /
+
+&SLCF PBZ=0.005, QUANTITY='TEMPERATURE', CELL_CENTERED=T /
+
+&DEVC ID='T', XYZ=0,0,0.01, QUANTITY='TEMPERATURE' /
+
+&TAIL /

Verification/Heat_Transfer/SFPE_Case_7.csv

@@ -0,0 +1,8 @@
+Time,Center,Side,Corner
+0,0,0,0
+1800,9,721,809
+3600,127,873,921
+5400,315,952,984
+7200,492,1005,1028
+9000,640,1045,1062
+10800,757,1077,1089