Merge pull request #13904 from mcgratta/master

FDS Verification: Remove CHECK_HT and clean up cases

Author: mcgratta

Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex

@@ -2252,7 +2252,7 @@ \section{Monin-Obukhov Similarity Profiles (\texorpdfstring{\ct{MO_velocity_prof
 \label{MO_velocity_profile_stable}
 \label{MO_velocity_profile_unstable}
 
-Atmospheric turbulence is affected by the stability of the boundary layer.  A stable layer (cool, heavy air at ground level) will suppress turbulence, while an unstable layer (warm, light air at ground level) will enhance turbulent mixing as buoyant plumes rise.  The theory governing these flows, Monin-Obukhov similarity theory, is discussed at some length in the FDS User Guide \cite{FDS_Users_Guide}.  Here we examine velocity profiles from a stable boundary layer and an unstable boundary layer.  When mean forcing is used for driving the wind field the Monin-Obukhov parameters determine the shape of the mean streamwise velocity profile.  These cases use a very tight \ct{DT_MEAN_FORCING_2} of 0.1 s in order to drive the flow field directly to the specified profile, therefore comfirming the target profile is being computed correctly in FDS.
+Atmospheric turbulence is affected by the stability of the boundary layer.  A stable layer (cool, heavy air at ground level) will suppress turbulence, while an unstable layer (warm, light air at ground level) will enhance turbulent mixing as buoyant plumes rise.  The theory governing these flows, Monin-Obukhov similarity theory, is discussed at some length in the FDS User Guide \cite{FDS_Users_Guide}.  Here we examine velocity profiles from a stable boundary layer and an unstable boundary layer. The Monin-Obukhov parameters determine the shape of the mean streamwise velocity profile.
 
 \begin{figure}[ht]
    \begin{tabular*}{\textwidth}{l@{\extracolsep{\fill}}r}
@@ -5941,11 +5941,11 @@ \subsection{Heating a Metal Sphere via Radiation and Convection}
 \label{particle_heating_convection}
 \label{particle_heating_radiation}
 
-A small metal sphere with mass, $m_{\rm s}=0.005$~kg, is suspended in a 1~m cube filled with $m_{\rm g}=0.318$~kg nitrogen with a specified specific heat, $c_p=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm g,i}=1073.15$~K. The metal has a specified specific heat, $c_{\rm s}=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm s,i}=293.15$~K. The walls of the box are adiabatic. In the first case, \ct{particle\_heating\_convection}, the sphere is heated via convection only and there is no radiation heat transfer. In the second case, \ct{particle\_heating\_radiation}, the sphere is heated via radiation only. The initial heat flux in both cases is approximately 75~kW/m$^2$, and the final temperature in both cases, $T_{\rm f}$, is found from solving an equation that equates the internal energy gained by the solid with the internal energy lost by the gas:
+Ten small metal spheres, each with mass, $m_{\rm s}=0.005$~kg, are suspended in a 1~m cube filled with $m_{\rm g}=0.318$~kg nitrogen with a specified specific heat, $c_p=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm g,i}=1073.15$~K. The metal has a specified specific heat, $c_{\rm s}=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm s,i}=293.15$~K. The walls of the box are adiabatic. In the first case, \ct{particle\_heating\_convection}, the spheres are heated via convection only and there is no radiation heat transfer. In the second case, \ct{particle\_heating\_radiation}, the spheres are heated via radiation only. The initial heat flux in both cases is approximately 75~kW/m$^2$, and the final temperature in both cases, $T_{\rm f}$, is found from solving an equation that equates the internal energy gained by the solid with the internal energy lost by the gas:
 \be
-   m_{\rm s} \, c_{\rm s} \, (T_{\rm f}-T_{\rm s,i}) = m_{\rm g} \, c_{v} \, (T_{\rm g,i}-T_{\rm f})  \quad ; \quad  c_v = c_p - \frac{R}{W} = 1 - \frac{8.3145}{28} \approx 0.703 \; \hbox{kJ/(kg}\cdot\hbox{K})
+   10 \, m_{\rm s} \, c_{\rm s} \, (T_{\rm f}-T_{\rm s,i}) = m_{\rm g} \, c_{v} \, (T_{\rm g,i}-T_{\rm f})  \quad ; \quad  c_v = c_p - \frac{R}{W} = 1 - \frac{8.3145}{28} \approx 0.703 \; \hbox{kJ/(kg}\cdot\hbox{K})
 \ee
-The final temperature, $T_{\rm f}=1056.1$~K or 782.9~$^\circ$C, as shown in Fig.~\ref{particle_heating_figs}.
+The final temperature, $T_{\rm f}=931.4$~K or 658.3~$^\circ$C, as shown in Fig.~\ref{particle_heating_figs}.
 \begin{figure}[ht]
 \includegraphics[height=2.2in]{SCRIPT_FIGURES/particle_heating_convection}
 \includegraphics[height=2.2in]{SCRIPT_FIGURES/particle_heating_radiation}

Utilities/Matlab/FDS_verification_dataplot_inputs.csv

@@ -486,10 +486,10 @@ d,particle_drag_U100_N1600,Sprinklers_and_Sprays/particle_drag_U100_N1600_git.tx
 d,particle_drag_U150_N1600,Sprinklers_and_Sprays/particle_drag_U150_N1600_git.txt,Sprinklers_and_Sprays/particle_drag_U150_N1600.csv,1,2,T,U,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_drag_U150_N1600_devc.csv,2,3,Time,U-VEL,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Gas phase velocity (particle\_drag\_F),Time (s),Velocity (m/s),0,2,1,0,200,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_drag_F,Absolute Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
 d,particle_drag_U10_N16,Sprinklers_and_Sprays/particle_drag_U10_N16_git.txt,Sprinklers_and_Sprays/particle_drag_U10_N16.csv,1,2,T,F,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_drag_U10_N16_devc.csv,2,3,Time,drag force,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Total drag force (particle\_drag\_A),Time (s),Force (N),0,100,1,0,1,-1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_drag_sum_A,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
 d,particle_flux,Sprinklers_and_Sprays/particle_flux_git.txt,Sprinklers_and_Sprays/particle_flux.csv,1,2,Time,mass,Expected (mass),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_flux_devc.csv,2,3,Time,mass,FDS (mass),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Particle Mass (particle\_flux),Time (s),Mass (kg),0,20,1,0,0.8,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/particle_flux,Relative Error,end,0.01,Sprinklers and Sprays,bs,b,TeX
-d,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
-f,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
-d,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
-f,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
+d,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
+f,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
+d,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
+f,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
 d,particle_isotropic_radi,Sprinklers_and_Sprays/particle_isotropic_radi_git.txt,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Delta e_ball,E\_w ball,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Q_rad ball,Q\_r ball,k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Enthalpy (particle\_isotropic\_radiation),Time (s),Enthalpy (kJ),0,50,1,0,12,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_isotropic_radi,Relative Error,end,0.015,Sprinklers and Sprays,ko,k,TeX
 f,particle_isotropic_radi,Sprinklers_and_Sprays/particle_isotropic_radi_git.txt,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Delta e_cyl,E\_w cyl,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Q_rad cyl,Q\_r cyl,r--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Enthalpy (particle\_isotropic\_radiation),Time (s),Enthalpy (kJ),0,50,1,0,12,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_isotropic_radi,Relative Error,end,0.015,Sprinklers and Sprays,ko,k,TeX
 f,particle_isotropic_radi,Sprinklers_and_Sprays/particle_isotropic_radi_git.txt,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Delta e_plate,E\_w plate,b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Q_rad plate,Q\_r plate,b--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Enthalpy (particle\_isotropic\_radiation),Time (s),Enthalpy (kJ),0,50,1,0,12,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_isotropic_radi,Relative Error,end,0.015,Sprinklers and Sprays,ko,k,TeX

Utilities/Matlab/scripts/atmospheric_boundary_layer.m

@@ -67,6 +67,11 @@
 
 T = T + (theta_0-T(12));
 
+ERROR = abs(u(end)-M2.data(end,2));
+if ERROR>2.
+   display(['Matlab Warning: atmospheric_boundary_layer Case ',num2str(i),' velocity out of tolerance. ERROR = ',num2str(ERROR),' m/s'])
+end
+
 plot(u,z,'ko'); hold on
 plot(M2.data(:,2),M2.data(:,1),'k-'); hold off
 set(gca,'FontName',Font_Name)
@@ -99,6 +104,11 @@
 set(gca,'Units',Plot_Units)
 set(gca,'Position',[Plot_X Plot_Y Plot_Width Plot_Height])
 
+ERROR = abs(T(end)-273.15-M2.data(end,3));
+if ERROR>1.0
+   display(['Matlab Warning: atmospheric_boundary_layer Case ',num2str(i),' temperature out of tolerance. ERROR = ',num2str(ERROR),' K'])
+end
+
 plot(T-273.15,z,'ko'); hold on
 plot(M2.data(:,3),M2.data(:,1),'k-'); hold off
 set(gca,'FontName',Font_Name)

Verification/Atmospheric_Effects/MO_velocity_profile_stable.fds

@@ -7,7 +7,7 @@
 
 &TIME T_END=10. /
 
-&MISC P_INF=90630, HUMIDITY=12.0, CHECK_HT=T /
+&MISC P_INF=90630, HUMIDITY=12.0 /
 &WIND L=69.4, U_STAR=0.369, Z_0=0.008, THETA_STAR=0.152, TMP_REF=30.8, Z_REF=1. /
 &RADI RADIATION=.FALSE. /
 

Verification/Complex_Geometry/geom_channel2.fds

@@ -15,7 +15,7 @@
 
 &RADI RADIATION=F /
 &PRES VELOCITY_TOLERANCE=1.E-6 /
-&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0., CHECK_HT=T /
+&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0. /
 &SURF ID='INLET', MASS_FLUX(1)=1., SPEC_ID(1)='TRACER', TAU_MF(1)=0./
 
 &VENT MB='XMIN', SURF_ID='INLET' /

Verification/Complex_Geometry/geom_channel_tmp.fds

@@ -18,7 +18,7 @@
 
 &RADI RADIATION=F /
 &PRES SOLVER='ULMAT'/
-&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0., CHECK_HT=T /
+&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0. /
 &SURF ID='INLET', FREE_SLIP=T, MASS_FRACTION(1)=1., SPEC_ID(1)='TRACER', VEL=-1, TAU_V=0., TMP_FRONT=600., RAMP_T='t1'/
 
 &RAMP ID='t1', T=0.,F=1./

Verification/Complex_Geometry/geom_channel_tmp2.fds

@@ -14,7 +14,7 @@
 
 &RADI RADIATION=F /
 &PRES VELOCITY_TOLERANCE=1.E-6 /
-&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0., CHECK_HT=T /
+&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0. /
 &SURF ID='INLET', TMP_FRONT=600., RAMP_T='t1'/
 
 &RAMP ID='t1', T=0.,F=1./

Verification/Complex_Geometry/geom_obst.fds

@@ -3,7 +3,6 @@
 &TIME T_END=1.0/
 &MESH IJK=15,15,15,XB=-1.0,1.0,-1.0,1.0,-1.0,1.0 /
 
-&MISC CHECK_HT=T /
 &REAC SOOT_YIELD=0.01,FUEL='PROPANE'/
 &SURF ID='BURNER',HRRPUA=600.0 / 
 

Verification/Heat_Transfer/back_wall_test_2.fds

@@ -1,11 +1,10 @@
-&HEAD CHID='back_wall_test_2', TITLE='Test 1-D heat transfer through rotated GEOM obstruction' /
+&HEAD CHID='back_wall_test_2_noht', TITLE='Test 1-D heat transfer through rotated GEOM obstruction' /
 
 &MESH XB=0,1,0,1,0,1, IJK=40,40,40 /
 &MESH XB=2,3,0,1,0,1, IJK=40,40,40 /
 &MESH XB=4,5,0,1,0,1, IJK=40,40,40 /
 
 &TIME T_END=20., DT=0.1 / 
-&MISC CFL_MAX=0.8, CHECK_HT=T / 
 &SURF ID='HOT', COLOR='RED', TMP_FRONT=1000., TAU_T=0., EMISSIVITY=1, HEAT_TRANSFER_COEFFICIENT=0 /
 &SURF ID='COLD',COLOR='BLUE',TMP_FRONT=20.,   TAU_T=0., EMISSIVITY=1, HEAT_TRANSFER_COEFFICIENT=0 /
 

Verification/Pressure_Solver/ulmat_2zone.fds

@@ -3,20 +3,11 @@
 &MESH IJK=10,10,5, XB=0,1,0,1,0,0.5/
 &MESH IJK=10,10,5, XB=0,1,0,1,0.5,1/
 
-#&MESH IJK=5,10,10, XB=0,0.5,0,1,0,1./
-#&MESH IJK=5,10,10, XB=0.5,1,0,1,0,1./
-
-#&MESH IJK=5,10,5, XB=0,0.5,0,1,0,0.5/
-#&MESH IJK=5,10,5, XB=0,0.5,0,1,0.5,1/
-#&MESH IJK=5,10,5, XB=0.5,1,0,1,0,0.5/
-#&MESH IJK=5,10,5, XB=0.5,1,0,1,0.5,1/
-
-
 &TIME T_END=10/
 
 &PRES CHECK_POISSON=T, SOLVER='ULMAT'/
 
-&MISC NOISE=F, CHECK_HT=T /
+&MISC NOISE=F /
 
 &RADI RADIATION=F/