firemodels · mcgratta · Jan 23, 2025 · Jan 23, 2025 · Jan 23, 2025
diff --git a/Manuals/FDS_Technical_Reference_Guide/Radiation_Chapter.tex b/Manuals/FDS_Technical_Reference_Guide/Radiation_Chapter.tex
@@ -335,9 +335,9 @@ \subsection{Radiation Contribution to Energy Equation}
 
 \subsection{Correction of the Emission Source Term}
 
-In calculations of limited spatial resolution, the source term, $I_{\rm b}$, defined in Eq.~(\ref{emission_source_term}) requires special treatment in the flaming region of the fire. Typical FDS calculations use grid cells that are tens of centimeters in size, and consequently the computed temperatures constitute a bulk average for a given grid cell and are considerably lower than the maximum temperature in a diffusion flame. Because of its fourth-power dependence on the temperature, the source term must be modeled in those grid cells where combustion occurs. Elsewhere, the computed temperature is used directly to compute the source term. It is assumed that this ``flaming region'' is where the local, nominal radiative loss is greater than a specified lower bound, $\chi_{\rm r} \dq'''>10$~kW/m$^3$. In this region, the global radiative fraction model is used. The emission source term is multiplied by a corrective factor, $C$:
+In calculations of limited spatial resolution, the source term, $I_{\rm b}$, defined in Eq.~(\ref{emission_source_term}) requires special treatment in the flaming region of the fire. Typical FDS calculations use grid cells that are tens of centimeters in size, and consequently the computed temperatures constitute a bulk average for a given grid cell and are considerably lower than the maximum temperature in a diffusion flame. Because of its fourth-power dependence on the temperature, the source term must be modeled in those grid cells where combustion occurs. Elsewhere, the computed temperature is used directly to compute the source term. It is assumed that this ``flaming region'' is where the local, nominal radiative loss is greater than a specified lower bound, $\chi_{\rm r} \dq'''>1$~kW/m$^3$. In this region, the global radiative fraction model is used. The emission source term is multiplied by a corrective factor, $C$:
 \be I_{\rm b,f}(\bx) = C \, \frac{\sigma \, T(\bx)^4}{\pi}   \quad ; \quad
-    C = \min \left( 100 \; , \; \max \left[0.1 \; , \; \frac{\sum_{\chi_{\rm r}\dq'''_{ijk}>10} \left( \chi_{\rm r} \, \dq'''_{ijk} + \kappa_{ijk} \, U_{ijk} \right) \, V_{ijk}}{\sum_{\chi_{\rm r} \dq'''_{ijk}>10} \left( 4 \, \kappa_{ijk} \, \sigma \, T^4_{ijk} \right) \, V_{ijk}} \right] \right) \label{corrected_emission_source_term}
+    C = \min \left( 100 \; , \; \max \left[0.1 \; , \; \frac{\sum_{\chi_{\rm r}\dq'''_{ijk}>1} \left( \chi_{\rm r} \, \dq'''_{ijk} + \kappa_{ijk} \, U_{ijk} \right) \, V_{ijk}}{\sum_{\chi_{\rm r} \dq'''_{ijk}>1} \left( 4 \, \kappa_{ijk} \, \sigma \, T^4_{ijk} \right) \, V_{ijk}} \right] \right) \label{corrected_emission_source_term}
 \ee
 When the source term defined in Eq.~(\ref{corrected_emission_source_term}) is substituted into Eq.~(\ref{net_emission}), the net radiative emission from the flaming region becomes the desired fraction of the total heat release rate. Note that this correction factor is bounded below by 0.1 and above by 100. These bounds are somewhat arbitrary, meant to prevent spurious behavior at the start of a simulation.
 

diff --git a/Manuals/FDS_User_Guide/FDS_User_Guide.tex b/Manuals/FDS_User_Guide/FDS_User_Guide.tex
@@ -5756,15 +5756,15 @@ \subsection{Radiation Option 2. Optically-Thin Limit; Specified Radiative Fracti
 \subsection[Radiation Option 3. Optically-Thick; Specified Radiative Fraction]{Radiation Option 3. Optically-Thick; Specified Radiative Fraction (LES Default)}
 \label{info:RTE_Source_Correction}
 
-In its normal operation, the RTE transfers energy from hot, emitting gases, like flames, to colder, absorbing gases like water vapor or soot particulate. The absorption coefficient, $\kappa$, computed using RadCal, governs both the emission and absorption of thermal radiation. Because flame temperatures are not well-resolved for typically large-scale fire simulations, the source term in the RTE is adjusted in grid cells for which the radiative fraction, $\chi_{\rm r}$, times the local heat release rate per unit volume, $\dot{q}'''$, is greater than 10~kW/m$^3$
+In its normal operation, the RTE transfers energy from hot, emitting gases, like flames, to colder, absorbing gases like water vapor or soot particulate. The absorption coefficient, $\kappa$, computed using RadCal, governs both the emission and absorption of thermal radiation. Because flame temperatures are not well-resolved for typically large-scale fire simulations, the source term in the RTE is adjusted in grid cells for which the radiative fraction, $\chi_{\rm r}$, times the local heat release rate per unit volume, $\dot{q}'''$, is greater than 1~kW/m$^3$
 \be
-   \chi_{\rm r} \, \dot{q}''' > 10 \; \hbox{kW/m}^3 \label{clip}
+   \chi_{\rm r} \, \dot{q}''' > 1 \; \hbox{kW/m}^3 \label{clip}
 \ee
 The adjustment ensures that the net radiative emission from the combusting region (i.e. the fire) is the specified \ct{RADIATIVE_FRACTION} multiplied by the total combustion energy generated in this region. Elsewhere, hot and cold gases emit and absorb thermal radiation according to their bulk temperature and radiative properties, in particular the absorption coefficient, $\kappa$.
 
 The net radiative loss from the computational domain is reported as a function of time in the column \ct{Q_RADI} in the output file \ct{CHID_hrr.csv}. The absolute value of \ct{Q_RADI} divided by the total heat release rate, \ct{HRR}, is usually not exactly equal to the specified \ct{RADIATIVE_FRACTION}. The reason for this is that the specified radiative fraction of the fire's energy can be reabsorbed by colder combustion products such as smoke and water vapor, thereby decreasing the absolute value of \ct{Q_RADI}. Or, hot layer smoke and combustion products can heat up and emit thermal radiation, adding to the absolute value of \ct{Q_RADI}.
 
-The correction factor that is applied to the RTE source term in the region defined by Eq.~(\ref{clip}) by default is bound between 1 and 100, meaning that the correction factor only increases the net radiative output of the combusting region, if necessary, to achieve the desired \ct{RADIATIVE_FRACTION}. However, you can change the default behavior of the correction as follows. First, you can force the RTE source term to be modified in all grid cells by changing the 10 in Eq.~(\ref{clip}) to -1 via \ct{QR_CLIP} on the \ct{RADI} line, in which case the solver will apply the radiative fraction to the entire domain, not just the cells where combustion occurs. This will essentially force the net radiative loss from the entire domain to obey the \ct{RADIATIVE_FRACTION}. Second, you can allow the RTE source term to increase or decrease in value to achieve the desired \ct{RADIATIVE_FRACTION} by changing the lower limit of the correction factor, \ct{C_MIN}, from its default value of 1 to, say, 0.5 on the \ct{RADI} line. The corresponding parameter, \ct{C_MAX}, limits the correction factor to 100. The time-varying correction factor can be output using a device with quantity \ct{RTE SOURCE CORRECTION FACTOR}. Since \ct{RTE SOURCE CORRECTION FACTOR} is a global value, the position of the device does not matter.
+The correction factor that is applied to the RTE source term in the region defined by Eq.~(\ref{clip}) by default is bound between 1 and 100, meaning that the correction factor only increases the net radiative output of the combusting region, if necessary, to achieve the desired \ct{RADIATIVE_FRACTION}. However, you can change the default behavior of the correction as follows. First, you can force the RTE source term to be modified in all grid cells by changing the 1 in Eq.~(\ref{clip}) to -1 via \ct{QR_CLIP} on the \ct{RADI} line, in which case the solver will apply the radiative fraction to the entire domain, not just the cells where combustion occurs. This will essentially force the net radiative loss from the entire domain to obey the \ct{RADIATIVE_FRACTION}. Second, you can allow the RTE source term to increase or decrease in value to achieve the desired \ct{RADIATIVE_FRACTION} by changing the lower limit of the correction factor, \ct{C_MIN}, from its default value of 1 to, say, 0.5 on the \ct{RADI} line. The corresponding parameter, \ct{C_MAX}, limits the correction factor to 100. The time-varying correction factor can be output using a device with quantity \ct{RTE SOURCE CORRECTION FACTOR}. Since \ct{RTE SOURCE CORRECTION FACTOR} is a global value, the position of the device does not matter.
 
 \subsection[Radiation Option 4. Optically-Thick; Unspecified Radiative Fraction]{Radiation Option 4. Optically-Thick; Unspecified Radiative Fraction (DNS Default)}
 
@@ -13104,7 +13104,7 @@ \section{\texorpdfstring{{\tt RADI}}{RADI} (Radiation Parameters)}
 \ct{NUMBER_RADIATION_ANGLES}      & Integer      & Section~\ref{info:RADI_Resolution}   &            & 100            \\ \hline
 \ct{OPTICALLY_THIN}                & Logical      & Section~\ref{info:CHI_R}             &            & \ct{F}  \\ \hline
 \ct{PATH_LENGTH }                  & Real         & Section~\ref{info:RadCal}            &   m        & 0.1            \\ \hline
-\ct{QR_CLIP}                       & Real         & Section~\ref{info:CHI_R}             & kW/m$^3$   & 10             \\ \hline
+\ct{QR_CLIP}                       & Real         & Section~\ref{info:CHI_R}             & kW/m$^3$   & 1              \\ \hline
 \ct{RADIATION}                      & Logical      & Section~\ref{info:radiation_off}     &            & \ct{T}   \\ \hline
 \ct{RADIATION_ITERATIONS}          & Integer      & Section~\ref{info:RADI_Resolution}   &            & 1              \\ \hline
 \ct{RADTMP}                         & Real         & Section~\ref{info:RADI_Two_Phase}    & $^\circ$C   & 900            \\ \hline

diff --git a/Source/cons.f90 b/Source/cons.f90
@@ -471,7 +471,7 @@ MODULE GLOBAL_CONSTANTS
 REAL(EB) :: RTE_SOURCE_CORRECTION_FACTOR=1._EB   !< Multiplicative factor used in correcting RTE source term
 REAL(EB) :: RAD_Q_SUM=0._EB   !< \f$ \sum_{ijk} \left( \chi_{\rm r} \dot{q}_{ijk}''' + \kappa_{ijk} U_{ijk} \right) V_{ijk} \f$
 REAL(EB) :: KFST4_SUM=0._EB   !< \f$ \sum_{ijk} 4 \kappa_{ijk} \sigma T_{ijk}^4 V_{ijk} \f$
-REAL(EB) :: QR_CLIP=10._EB    !< Lower bound of \f$ \chi_{\rm r} \dot{q}_{ijk}''' \f$ below which no source correction is made
+REAL(EB) :: QR_CLIP=1._EB     !< Lower bound of \f$ \chi_{\rm r} \dot{q}_{ijk}''' \f$ below which no source correction is made
 REAL(EB) :: C_MAX=100._EB     !< Maximum value of RAD_Q_SUM/KFST4_SUM
 REAL(EB) :: C_MIN=0.1_EB      !< Minimum value of RAD_Q_SUM/KFST4_SUM
 

diff --git a/Source/main.f90 b/Source/main.f90
@@ -1517,7 +1517,7 @@ END SUBROUTINE PRESSURE_ITERATION_SCHEME
 
 SUBROUTINE CALCULATE_RTE_SOURCE_CORRECTION_FACTOR
 
-REAL(EB), PARAMETER :: WGT=0.5_EB
+REAL(EB), PARAMETER :: WGT=0.9_EB
 REAL(EB) :: TNOW
 
 TNOW = CURRENT_TIME()

diff --git a/Source/radi.f90 b/Source/radi.f90
@@ -3475,7 +3475,7 @@ SUBROUTINE RADIATION_FVM
 REAL(EB) :: XID,YJD,ZKD,KAPPA_PART_SINGLE,DLF,DLA(3),TSI,TMP_EXTERIOR,TEMP_ORIENTATION(3)
 REAL(EB), ALLOCATABLE, DIMENSION(:) :: ZZ_GET
 INTEGER :: IID,JJD,KKD,IP
-LOGICAL :: UPDATE_INTENSITY, UPDATE_QRW2
+LOGICAL :: UPDATE_INTENSITY
 REAL(EB), POINTER, DIMENSION(:,:,:) :: IL,UIIOLD,KAPPA_PART,KFST4_PART,EXTCOE,SCAEFF,SCAEFF_G,IL_UP
 REAL(EB), POINTER, DIMENSION(:)     :: OUTRAD_W,INRAD_W,OUTRAD_F,INRAD_F,IL_F
 TYPE (OMESH_TYPE), POINTER :: M2
@@ -3574,8 +3574,6 @@ SUBROUTINE RADIATION_FVM
    ENDDO
 ENDIF
 
-UPDATE_QRW2 = .FALSE.
-
 ! Loop over spectral bands
 
 BAND_LOOP: DO IBND = 1,NUMBER_SPECTRAL_BANDS