Merge pull request #12472 from mcgratta/master

FDS Source: Minor fix to HT3D thin obstruction noding

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -2208,7 +2208,7 @@ \subsubsection{Output for Convective Heat Transfer Regime}
 \subsubsection{Specified Convective Heat Transfer Coefficient}
 
 To specify the convective heat transfer coefficient, set it to a constant using \\
-{\ct HEAT\_TRANSFER\_COEFFICIENT} on the {\ct SURF} line in units of \si{W/(m$^2$.K)} with optional time dependent ramp using {\ct RAMP\_HEAT\_TRANSFER\_COEFFICIENT}. If the back side of the solid obstruction faces the exterior of the computational domain and the solid conducts heat, the heat transfer coefficient of the back side may be specified using {\ct HEAT\_TRANSFER\_COEFFICIENT\_BACK} with optional time dependent ramp ramp using {\ct RAMP\_HEAT\_TRANSFER\_COEFFICIENT\_BACK}. This back side condition is appropriate for a {\ct SURF} line with {\ct BACKING='VOID'} or {\ct BACKING='EXPOSED'}. 
+{\ct HEAT\_TRANSFER\_COEFFICIENT} on the {\ct SURF} line in units of \si{W/(m$^2$.K)} with optional time dependent ramp using {\ct RAMP\_HEAT\_TRANSFER\_COEFFICIENT}. If the back side of the solid obstruction faces the exterior of the computational domain and the solid conducts heat, the heat transfer coefficient of the back side may be specified using {\ct HEAT\_TRANSFER\_COEFFICIENT\_BACK} with optional time dependent ramp ramp using {\ct RAMP\_HEAT\_TRANSFER\_COEFFICIENT\_BACK}. This back side condition is appropriate for a {\ct SURF} line with {\ct BACKING='VOID'} or {\ct BACKING='EXPOSED'}.
 
 \subsubsection{Specifying the Heat Flux at a Solid Surface}
 \label{info:net_and_convective_heat_flux}
@@ -2535,7 +2535,7 @@ \subsubsection{Example: Unprotected Structural Steel}
 \subsection{Surface Linings}
 \label{info:linings}
 
-In some applications, you might want to perform 3-D heat conduction within collections of abutting obstructions, each with a {\ct MATL\_ID}. The obstructions might have some outer lining, like fabric on sofa cushions, for example, and this lining is thinner than a single gas phase grid cell. In other words, this lining material cannot be represented by an {\ct OBST}. However, you can specify one or more layers of multi-component materials using {\ct MATL\_ID}, {\ct MATL\_MASS\_FRACTION}, and {\ct THICKNESS} on the {\ct SURF} line(s) that are assigned to the {\ct OBST}s that also have specified {\ct MATL\_ID}s and {\ct MATL\_MASS\_FRACTION}s. 
+In some applications, you might want to perform 3-D heat conduction within collections of abutting obstructions, each with a {\ct MATL\_ID}. The obstructions might have some outer lining, like fabric on sofa cushions, for example, and this lining is thinner than a single gas phase grid cell. In other words, this lining material cannot be represented by an {\ct OBST}. However, you can specify one or more layers of multi-component materials using {\ct MATL\_ID}, {\ct MATL\_MASS\_FRACTION}, and {\ct THICKNESS} on the {\ct SURF} line(s) that are assigned to the {\ct OBST}s that also have specified {\ct MATL\_ID}s and {\ct MATL\_MASS\_FRACTION}s.
 
 The way linings are handled depends on the size of the obstruction. If the obstruction is a flat plate whose depth is less than one-half a grid cell and is represented as a zero-cell thick sheet, then any layer materials are assumed to {\em add} to the thickness of the underlying plate. For example, if a plate of steel is 1~cm thick and it is assumed to be coated on both sides by layers of insulation that are 2~cm thick, the assembly is still treated as a thin sheet from the standpoint of the gas phase computation, but the depth is taken as 5~cm for the purpose of the heat transfer computation.
 
@@ -2557,12 +2557,10 @@ \subsubsection{Example: Insulated Structural Steel}
 
 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 applied to each side via a {\ct SURF} line. 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.
 \begin{lstlisting}
-&OBST XB=0.12,1.88,0.51,0.53,0.21,0.79, SURF_ID6='STEEL SLAB','STEEL SLAB',
-      'WEB INSULATION','WEB INSULATION','STEEL SLAB','STEEL SLAB', MATL_ID='STEEL' /
-&SURF ID='STEEL SLAB', HT3D=T, COLOR='BLACK', CELL_SIZE=0.1 /
-&SURF ID='WEB INSULATION', MATL_ID='STUFF', COLOR='BEIGE', THICKNESS=0.01, HT3D=T /
+&OBST XB=..., SURF_ID='INSULATION', MATL_ID='STEEL', CELL_SIZE=0.1 /
+&SURF ID='INSULATION', MATL_ID='STUFF', COLOR='BEIGE', THICKNESS=0.01, HT3D=T /
 \end{lstlisting}
-Note that the {\ct CELL\_SIZE} refers to heat transfer in the lateral direction. The heat transfer in the direction normal to the insulation and thin steel plate will be gridded according to the {\ct SURF} {\ct 'WEB INSULATION'}. 
+Note that the {\ct CELL\_SIZE} refers to heat transfer in the lateral direction. The heat transfer in the direction normal to the insulation and thin steel plate will be gridded according to parameters on the {\ct SURF} line {\ct 'INSULATION'}.
 
 
 

Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex

@@ -4507,28 +4507,26 @@ \subsection{Energy Conservation in a 3-D Solid (\texorpdfstring{\textct{ht3d\_en
 
 \FloatBarrier
 
-\subsection{Heat Transfer to a Steel Beam (\texorpdfstring{\textct{ht3d\_beam\_heating}}{ht3d\_beam\_heating})}
+\subsection{Heat Transfer to an I-Beam (\texorpdfstring{\textct{ht3d\_beam\_heating}}{ht3d\_beam\_heating})}
 \label{ht3d_beam_heating}
 
-A steel I-beam 1.76~m long, with flanges that are 0.22~m wide and 0.03~m thick and a web that is 0.6~m deep and 0.02~m thick, is heated within a sealed compartment that is filled with a gas that has a volume $V_{\rm g}=0.8$~m$^3$, initial temperature $T_{\rm g,0}=1073.15$~K, pressure $p_0=101.325$~kPa, molecular weight $W=0.028$~kg/mol, specific heat $c_{\rm g}=1$~kJ/(kg$\cdot$K), and density $\rho_{\rm g}=p_0 W/(RT_{\rm g,0})\approx0.318$~kg/m$^3$. $R\approx0.008314$~kJ/(mol$\cdot$K) is the universal gas constant. The steel has a density $\rho_{\rm s}=7500$~kg/m$^3$, specific heat $c_{\rm s}=0.5$~kJ/(kg$\cdot$K), and volume $V_{\rm s}=0.041536$~kg/m$^3$, and initial temperature $T_{\rm s,0}=293.15$~K. Initially, the gas and steel have internal energies given by:
+An aluminum I-beam $L=1.76$~m long and $d=0.58$~m deep (i.e.~in height), with flanges that are $b_{\rm f}=0.22$~m wide and $t_{\rm f}=0.005$~m thick and a web that is $t_{\rm w}=0.006$~m thick, is heated within a sealed compartment with volume $V_{\rm c}=0.8$~m$^3$, initial temperature $T_{\rm g,0}=1073.15$~K, pressure $p_0=101.325$~kPa, molecular weight $W=0.028$~kg/mol, specific heat $c_{\rm g}=1$~kJ/(kg$\cdot$K), and density $\rho_{\rm g}=p_0 W/(RT_{\rm g,0})\approx0.318$~kg/m$^3$. $R\approx0.008314$~kJ/(mol$\cdot$K) is the universal gas constant. The aluminum has a density $\rho_{\rm a}=2700$~kg/m$^3$, specific heat $c_{\rm a}=0.9$~kJ/(kg$\cdot$K), and initial temperature $T_{\rm a,0}=293.15$~K. Initially, the gas and aluminum have internal energies given by:
 \be
-   U_{\rm g} = \rho_{\rm g} V_{\rm g} c_{\rm g} T_{\rm g,0} - p_0 V_{\rm g} \approx 191.95 \; \hbox{kJ} \quad ; \quad
-   U_{\rm s} = \rho_{\rm s} V_{\rm s} c_{\rm s} T_{\rm s,0} \approx 45661.04 \; \hbox{kJ}
+   U_{\rm g} = \rho_{\rm g} V_{\rm g} c_{\rm g} T_{\rm g,0} - p_0 V_{\rm g} \approx 189.56 \; \hbox{kJ} \quad ; \quad
+   U_{\rm a} = \rho_{\rm a} V_{\rm a} c_{\rm a} T_{\rm a,0} \approx 7046.0 \; \hbox{kJ}
 \ee
 After several minutes, the equilibrium temperature is expected to be:
 \be
-   T = \frac{U_{\rm g}+U_{\rm s}}{\rho_{\rm s} V_{\rm s} c_{\rm s} + \rho_{\rm g} V_{\rm g} c_{\rm g} - n R} \approx 294.045 \; \hbox{K} \; \equiv 20.895 \; ^\circ \hbox{C}
+   T = \frac{U_{\rm g}+U_{\rm a}}{\rho_{\rm a} V_{\rm a} c_{\rm a} + \rho_{\rm g} V_{\rm g} c_{\rm g} - n R} \approx 298.84 \; \hbox{K} \; \equiv 25.69 \; ^\circ \hbox{C}
 \ee
-where $n=p_0 V_{\rm g}/(RT_{\rm g,0}) \approx 9.0847$ is the number of moles of gas.
+where $n=p_0 V_{\rm g}/(RT_{\rm g,0}) \approx 8.973$ is the number of moles of gas. The left hand plot in Fig.~\ref{fig:ht3d_beam_heating} displays the average surface temperature of the I-beam.
 
-The left hand plot in Fig.~\ref{fig:ht3d_beam_heating} displays the average surface temperature of the steel beam. Its computed value is approximately 4~\% lower than the exact value because the overlap of the web and flanges is not accounted for with the current 3-D heating algorithm in FDS.
-
-In a second case, a steel plate that is 1.76~m long, 0.58~m wide, and 0.02~m thick is coated with a 1~cm thick layer of thermal insulation (density $\rho_{\rm i}=20$~kg/m$^3$, specific heat $c_{\rm i}=1.0$~kJ/(kg$\cdot$K)) and exposed to the same thermal environment. The expected equilibrium temperature is 21.808~$^\circ$C. The right hand plot in Fig.~\ref{fig:ht3d_beam_heating} compares the computed and exact temperatures.
+In a second case, the I-beam is coated with a 1~cm thick layer of thermal insulation (density $\rho_{\rm i}=20$~kg/m$^3$, specific heat $c_{\rm i}=1.0$~kJ/(kg$\cdot$K)) and exposed to the same thermal environment. The expected equilibrium temperature is 25.62~$^\circ$C. The right hand plot in Fig.~\ref{fig:ht3d_beam_heating} compares the computed and exact temperatures.
 
 \begin{figure}[ht]
 \includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_beam_heating_1}
 \includegraphics[height=2.2in]{SCRIPT_FIGURES/ht3d_beam_heating_2}
-\caption[The \textct{ht3d\_beam\_heating} test cases]{(Left) Average surface temperature of a steel beam immersed in hot gas. (Right) Average temperature of an insulated steel plate exposed to the same conditions.}
+\caption[The \textct{ht3d\_beam\_heating} test cases]{(Left) Average surface temperature of an I-beam immersed in hot gas. (Right) Average temperature of an insulated I-beam exposed to the same conditions.}
 \label{fig:ht3d_beam_heating}
 \end{figure}
 

Source/init.f90

@@ -1992,8 +1992,8 @@ SUBROUTINE INITIALIZE_HT3D_WALL_CELLS(NM)
       DO IOR=1,3
          IF (ABS(BC%IOR)==IOR) CYCLE
          IF (.NOT.IOR_AVOID(-IOR) .AND. .NOT.IOR_AVOID(IOR)) THEN
-            WRITE(LU_ERR,'(A,I0,A,I0)') 'ERROR(424): HT3D thin solid must have at least one face exposed, Mesh=',NM,&
-                                        ', IOR=',IOR
+            WRITE(LU_ERR,'(7(A,I0))') 'ERROR(424): HT3D thin solid must have at least one face exposed in direction ',IOR,&
+                                      ': Mesh=',NM,', IOR=',BC%IOR,', IIG=',BC%IIG,', JJG=',BC%JJG,', KKG=',BC%KKG,', I=',I
             STOP_STATUS = SETUP_STOP
             RETURN
          ENDIF
@@ -4187,6 +4187,10 @@ SUBROUTINE FIND_THIN_WALL_BACK_INDEX(NM,ITW)
          CASE(2) ; ONE_D%LAYER_THICKNESS(1) = OB%UNDIVIDED_INPUT_LENGTH(2)
          CASE(3) ; ONE_D%LAYER_THICKNESS(1) = OB%UNDIVIDED_INPUT_LENGTH(3)
       END SELECT
+      IF (OB%CELL_SIZE>0._EB) THEN
+         ONE_D%CELL_SIZE(1) = OB%CELL_SIZE
+         ONE_D%STRETCH_FACTOR(1) = 1._EB
+      ENDIF
       EXIT FIND_BACK_THIN_WALL_CELL
    ENDIF
 

Utilities/Matlab/FDS_verification_dataplot_inputs.csv

@@ -258,8 +258,8 @@ d,hrrpuv_reac_single,Species/hrrpuv_reac_single_git.txt,Species/hrrpuv_reac_sing
 d,hrrpuv_reac_soot,Species/hrrpuv_reac_soot_git.txt,Species/hrrpuv_reac_soot_devc.csv,2,3,Time,hrrpuv|q1,hrrpuv|q1,ko|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Species/hrrpuv_reac_soot_devc.csv,2,3,Time,q2|q1+q2,q2|q1+q2,g-|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Reaction Test (hrrpuv\_reac\_soot),Time (s),Heat Release Rate (kW/m³),0,5,1,0,2000,1,no,0.05 0.90,West,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/hrrpuv_reac_soot,Relative Error,end_1_2,1.00E-05,Species,kd,k,TeX
 d,ht1d_pile,Heat_Transfer/ht1d_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/ht1d_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  (ht1d\_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/ht1d_pile,Relative Error,slope,1.00E-02,Heat Transfer,r^,r,TeX
 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_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,10,60,20,30,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_beam_heating_1,Relative Error,end,4.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,10,60,20,30,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ht3d_beam_heating_2,Relative Error,end,4.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

Verification/Heat_Transfer/ht3d_beam_heating.py

@@ -0,0 +1,37 @@
+L       = 1.76
+b_f     = 0.22
+t_f     = 0.005
+t_w     = 0.006
+t_i     = 0.01
+d       = 0.58
+rho_a   = 2700.
+rho_i   = 20.
+V_a     = L*(b_f*t_f*2+(d-2*t_f)*t_w)
+V_c     = 0.8
+V_i     = (L+2*t_i)*( (b_f+2*t_i)*(t_f+2*t_i)*2 + (d+2*t_i-2*(t_f+2*t_i))*(t_w+2*t_i) ) - V_a
+p_0     = 101.325
+R       = 0.008314
+W       = 0.028
+T_g_0   = 1073.15
+rho_g   = p_0*W/(R*T_g_0)
+c_g     = 1.
+c_a     = 0.9
+c_i     = 1.
+T_a_0   = 293.15
+V_g_1   = V_c - V_a
+n_1     = p_0*V_g_1/(R*T_g_0)
+U_g     = rho_g*V_g_1*c_g*T_g_0 - p_0*V_g_1
+U_a     = rho_a*V_a*c_a*T_a_0
+U_i     = rho_i*V_i*c_i*T_a_0
+T_1     = (U_g+U_a)/(rho_a*V_a*c_a + rho_g*V_g_1*c_g - n_1*R)
+V_g_2   = V_c - V_a - V_i
+n_2     = p_0*V_g_2/(R*T_g_0)
+T_2     = (U_g+U_a+U_i)/(rho_i*V_i*c_i + rho_a*V_a*c_a + rho_g*V_g_2*c_g - n_2*R)
+print('n_1 is ' + str(n_1) + ' mol')
+print('n_2 is ' + str(n_2) + ' mol')
+print('U_g is ' + str(U_g) + ' kJ')
+print('U_a is ' + str(U_a) + ' kJ')
+print('U_i is ' + str(U_i) + ' kJ')
+print('T_1 is ' + str(T_1-273.15) + ' C')
+print('T_2 is ' + str(T_2-273.15) + ' C')
+

Verification/Heat_Transfer/ht3d_beam_heating_1.csv

@@ -1,4 +1,4 @@
-Time,U,T
-240,-139.3549,20.89468
-270,-139.3549,20.89468
-300,-139.3549,20.89468
+Time,T
+540,25.690
+570,25.690
+600,25.690

Verification/Heat_Transfer/ht3d_beam_heating_1.fds

@@ -8,24 +8,23 @@
 
 &RADI RADIATION=F /
 
-&MISC STRATIFICATION=F, GVEC=0,0,0 /
-
-&OBST XB=0.12,1.88,0.51,0.53,0.21,0.79, SURF_ID='STEEL SLAB' /
-&OBST XB=0.12,1.88,0.41,0.63,0.76,0.79, SURF_ID='STEEL SLAB' /
-&OBST XB=0.12,1.88,0.41,0.63,0.21,0.24, SURF_ID='STEEL SLAB' /
+&OBST XB=0.12,1.88,0.510,0.516,0.210,0.790, MATL_ID='ALUMINUM', SURF_ID='ALUMINUM SLAB', CELL_SIZE=0.1 /
+&OBST XB=0.12,1.88,0.403,0.623,0.785,0.790, MATL_ID='ALUMINUM', SURF_ID='ALUMINUM SLAB', CELL_SIZE=0.1 /
+&OBST XB=0.12,1.88,0.403,0.623,0.210,0.215, MATL_ID='ALUMINUM', SURF_ID='ALUMINUM SLAB', CELL_SIZE=0.1 /
 
 &INIT XB=0.00,2.00,0.30,0.70,0.00,1.00, TEMPERATURE=800 /
 
-&SURF ID='STEEL SLAB', MATL_ID='STEEL', HT3D=T, COLOR='BLACK', CELL_SIZE=0.1, HEAT_TRANSFER_COEFFICIENT=20 /
+&SURF ID='ALUMINUM SLAB', HT3D=T, COLOR='BLACK' /
 &SURF ID='WALL', COLOR='GRAY', ADIABATIC=T, DEFAULT=T /
 
-&MATL ID='STEEL', DENSITY=7500, SPECIFIC_HEAT=0.5, CONDUCTIVITY=50 /
+&MATL ID='ALUMINUM', DENSITY=2700, SPECIFIC_HEAT=0.9, CONDUCTIVITY=240 /
 
 &BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=T /
 &SLCF PBY=0.5, QUANTITY='TEMPERATURE', CELL_CENTERED=T /
 
+&DUMP DT_DEVC=5. /
+
 &DEVC ID='temp', QUANTITY='TEMPERATURE', SPATIAL_STATISTIC='MEAN', XB=0.0,2.0,0.3,0.7,0.0,1.0 /
 &DEVC ID='wall temp', QUANTITY='WALL TEMPERATURE', SPATIAL_STATISTIC='MEAN', XB=0.05,1.95,0.35,0.65,0.05,0.95 /
-&DEVC ID='U', QUANTITY='INTERNAL ENERGY', SPATIAL_STATISTIC='VOLUME INTEGRAL', XB=0.0,2.0,0.3,0.7,0.0,1.0, RELATIVE=T /
 
 &TAIL /

Verification/Heat_Transfer/ht3d_beam_heating_2.csv

@@ -1,4 +1,4 @@
-Time,U,T
-3000,-139.1699,21.80814
-3300,-139.1699,21.80814
-3600,-139.1699,21.80814
+Time,T
+540,25.621
+570,25.621
+600,25.621