Merge pull request #15594 from rmcdermo/mw_cor_2

FDS Source: IMPORTANT, add back FLUX_LIMITER_MW_CORRECTION with 2 ghost cells as default, remove single coefficient scheme

Author: rmcdermo

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -9139,7 +9139,7 @@ \section{Flux Limiters}
 Central differencing                & \ct{'CENTRAL'}  \\
 Godunov                             & \ct{'GODUNOV'}  \\
 Superbee (VLES, SVLES default)      & \ct{'SUPERBEE'}  \\
-CHARM (DNS, LES default)            & \ct{'CHARM'}  \\ 
+CHARM (DNS, LES default)            & \ct{'CHARM'}  \\
 \hline
 \end{tabular}
 \end{table}
@@ -14426,8 +14426,8 @@ \chapter{Error Codes}
 610 \> \ct{HOLE ... Cannot overlap HOLEs with a DEVC or CTRL_ID.} \> Section~\ref{info:HOLE}  \\
 611 \> \ct{OBST ... has a BULK_DENSITY but zero volume.} \> Section~\ref{info:BURN_AWAY}  \\
 612 \> \ct{OBST ... must have a volume to be assigned HT3D.} \> Section~\ref{checkerboard}  \\
-614 \> \ct{OBST_ID ... cannot have a SURF with NODE_ID} \> Section~\ref{info:hvac_geom}  \\    
-615 \> \ct{OBST ... cannot overlap OBST ...}            \> Section~\ref{info:OVERLAY}  \\    
+614 \> \ct{OBST_ID ... cannot have a SURF with NODE_ID} \> Section~\ref{info:hvac_geom}  \\
+615 \> \ct{OBST ... cannot overlap OBST ...}            \> Section~\ref{info:OVERLAY}  \\
     \>                                                  \>                          \\
 701 \> \ct{problem with GEOM, local SURF_ID index ... out of bounds.} \> Section~\ref{info:GEOM_Basics}  \\
 702 \> \ct{problem with GEOM, SURF_IDS not defined properly.} \> Section~\ref{info:first_geom}  \\

Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex

@@ -5141,7 +5141,7 @@ \subsection{Case 8: 2-D Heat Transfer with Temperature-Dependent Conductivity}
 \FloatBarrier
 
 
-\subsection{Case 9: 2-D Heat Transfer in a Composite Section with Temperature-Dependent Conductivity}
+\subsection{Case 9: 2-D Heat Transfer in a Composite Section, Variable Conductivity}
 \label{SFPE_Case_9}
 
 A hollow square metal tube ($\rho=7850$~kg/m$^3$, $c=0.6$~kJ/(kg~K), $\epsilon=0.8$) is filled with an insulation material ($k=0.05$~W/(m~K), $\rho=50$~kg/m$^3$, $c=1$~kJ/(kg~K)). The thermal conductivity of the metal tube varies linearly with temperature such that its value is 54.7~W/(m~K) at 0~°C, 27.3~W/(m~K) at 800~°C, and 27.3~W/(m~K) at 1200~°C. The tube walls are 0.5~mm thick, and the exterior width of the assembly is 0.201~m. The surrounding air temperature is 1000~°C, and the initial temperature of the assembly is 0~°C. Assuming that the heating is by convection and radiation, and that the heat transfer coefficient is 10~W/(m$^2$~K), calculate the temperature at the center of the tube as a function of time (Fig.~\ref{fig:SFPE_Case_9}).

Source/cons.f90

@@ -278,6 +278,7 @@ MODULE GLOBAL_CONSTANTS
 LOGICAL :: STORE_FIRE_RESIDENCE=.FALSE.             !< Flag for tracking residence time of spreading fire front over a surface
 LOGICAL :: STORE_LS_SPREAD_RATE=.FALSE.             !< Flag for outputting local level set spread rate magnitude
 LOGICAL :: TEST_NEW_CHAR_MODEL=.FALSE.              !< Flag to envoke new char model
+LOGICAL :: FLUX_LIMITER_MW_CORRECTION=.TRUE.        !< Flag for MW correction ensure consistent equation of state at face
 
 INTEGER, ALLOCATABLE, DIMENSION(:) :: CHANGE_TIME_STEP_INDEX      !< Flag to indicate if a mesh needs to change time step
 INTEGER, ALLOCATABLE, DIMENSION(:) :: SETUP_PRESSURE_ZONES_INDEX  !< Flag to indicate if a mesh needs to keep searching for ZONEs

Source/data.f90

@@ -1360,42 +1360,6 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 OUTPUT_QUANTITY(553)%UNITS  = 'm/s'
 OUTPUT_QUANTITY(553)%SHORT_NAME  = 'V_LS'
 
-OUTPUT_QUANTITY(560)%NAME  = 'BFX'
-OUTPUT_QUANTITY(560)%UNITS  = ''
-OUTPUT_QUANTITY(560)%SHORT_NAME  = 'bfx'
-OUTPUT_QUANTITY(560)%CELL_POSITION = CELL_FACE
-OUTPUT_QUANTITY(560)%IOR = 1
-
-OUTPUT_QUANTITY(561)%NAME  = 'BFY'
-OUTPUT_QUANTITY(561)%UNITS  = ''
-OUTPUT_QUANTITY(561)%SHORT_NAME  = 'bfy'
-OUTPUT_QUANTITY(561)%CELL_POSITION = CELL_FACE
-OUTPUT_QUANTITY(561)%IOR = 2
-
-OUTPUT_QUANTITY(562)%NAME  = 'BFZ'
-OUTPUT_QUANTITY(562)%UNITS  = ''
-OUTPUT_QUANTITY(562)%SHORT_NAME  = 'bfz'
-OUTPUT_QUANTITY(562)%CELL_POSITION = CELL_FACE
-OUTPUT_QUANTITY(562)%IOR = 3
-
-OUTPUT_QUANTITY(563)%NAME  = 'BFX MINUS'
-OUTPUT_QUANTITY(563)%UNITS  = ''
-OUTPUT_QUANTITY(563)%SHORT_NAME  = 'bfx-'
-OUTPUT_QUANTITY(563)%CELL_POSITION = CELL_FACE
-OUTPUT_QUANTITY(563)%IOR = 1
-
-OUTPUT_QUANTITY(564)%NAME  = 'BFY MINUS'
-OUTPUT_QUANTITY(564)%UNITS  = ''
-OUTPUT_QUANTITY(564)%SHORT_NAME  = 'bfy-'
-OUTPUT_QUANTITY(564)%CELL_POSITION = CELL_FACE
-OUTPUT_QUANTITY(564)%IOR = 2
-
-OUTPUT_QUANTITY(565)%NAME  = 'BFZ MINUS'
-OUTPUT_QUANTITY(565)%UNITS  = ''
-OUTPUT_QUANTITY(565)%SHORT_NAME  = 'bfz-'
-OUTPUT_QUANTITY(565)%CELL_POSITION = CELL_FACE
-OUTPUT_QUANTITY(565)%IOR = 3
-
 ! Boundary Quantities (Negative indices)
 
 OUTPUT_QUANTITY(-1)%NAME = 'RADIATIVE HEAT FLUX'