FDS Source: add OXIDATIVE HRRPUA solid phase output. Issue #13664

Author: ericvmueller

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -2868,7 +2868,7 @@ \subsection{Reaction Rates}
 \be
 X_{\rm O_2}(x) = X_{\rm O_2,f}\exp(-x/L_{\rm g})
 \ee
-where $L_{\rm g}$ is the gas diffusion length scale and $X_{\rm O_2,f}$ is computed to satisfy simultaneous mass transfer and solid phase reaction rates (see FDS Tech Guide \cite{FDS_Math_Guide}). $n_{{\rm O_2},ij}$ is prescribed under the name {\ct N\_O2(j)} on the {\ct MATL} line of the $i$th material. It is zero by default. $L_{\rm g}$ is prescribed under the name {\ct GAS\_DIFFUSION\_DEPTH(j)}, and it is 0.001 m by default.
+where $L_{\rm g}$ is the gas diffusion length scale and $X_{\rm O_2,f}$ is computed to satisfy simultaneous mass transfer and solid phase reaction rates (see FDS Tech Guide \cite{FDS_Math_Guide}). $n_{{\rm O_2},ij}$ is prescribed under the name {\ct N\_O2(j)} on the {\ct MATL} line of the $i$th material. It is zero by default. $L_{\rm g}$ is prescribed under the name {\ct GAS\_DIFFUSION\_DEPTH(j)}, and it is 0.001 m by default. The heat release rate per unit area associated with oxidation reactions can be output using the solid phase {\ct QUANTITY} called {\ct 'OXIDATIVE HRRPUA'}. This is an integral in depth of the contributions from all reactions and all materials for a given solid surface.
 
 You can specity {\ct MAX\_REACTION\_RATE(j)} (\si{kg/m^3.s}) to limit the reaction rate, $r_{ij}$, below a specified value.
 
@@ -11472,6 +11472,7 @@ \section{Solid Phase Output Quantities}
 {\ct NORMALIZED HEAT RELEASE RATE}              & Section~\ref{info:material_components}        & W/g                & D            \\ \hline
 {\ct NORMALIZED MASS}$^4$                       & Section~\ref{info:material_components}        &                    & D            \\ \hline
 {\ct NORMALIZED MASS LOSS RATE}$^4$             & Section~\ref{info:material_components}        & 1/s                & D            \\ \hline
+{\ct OXIDATIVE HRRPUA}                          & Section~\ref{sec:reaction_rates}              & \unit{kW/m^2}      & B,D          \\ \hline
 {\ct PRESSURE COEFFICIENT}                      & Section~\ref{info:pressure_coefficient}       &                    & B,D          \\ \hline
 {\ct RADIATIVE HEAT FLUX}                       & Section~\ref{info:heat_flux}                  & \unit{kW/m^2}           & B,D          \\ \hline
 {\ct RADIOMETER}                                & Section~\ref{info:heat_flux}                  & \unit{kW/m^2}           & B,D          \\ \hline

Source/data.f90

@@ -1604,6 +1604,10 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 OUTPUT_QUANTITY(-79)%SHORT_NAME = 'O2_iter'
 OUTPUT_QUANTITY(-79)%BNDF_APPROPRIATE = .TRUE.
 
+OUTPUT_QUANTITY(-80)%NAME = 'OXIDATIVE HRRPUA'
+OUTPUT_QUANTITY(-80)%UNITS= 'kW/m2'
+OUTPUT_QUANTITY(-80)%SHORT_NAME = 'hrrpua_O2'
+
 ! Condensation
 OUTPUT_QUANTITY(-100)%NAME = 'CONDENSATION HEAT FLUX'
 OUTPUT_QUANTITY(-100)%UNITS= 'kW/m2'

Source/dump.f90

@@ -9252,6 +9252,8 @@ REAL(EB) FUNCTION SOLID_PHASE_OUTPUT(INDX,Y_INDEX,Z_INDEX,PART_INDEX,OPT_WALL_IN
    CASE(79) ! SURFACE OXYGEN ITERATIONS
       SOLID_PHASE_OUTPUT = 0
       IF (SF%INCLUDE_BOUNDARY_PROP2_TYPE) SOLID_PHASE_OUTPUT = B2%Y_O2_ITER
+   CASE(80) ! OXIDATIVE HRRPUA
+      SOLID_PHASE_OUTPUT = B1%Q_DOT_O2_PP*0.001_EB
 
    CASE(100) ! CONDENSATION HEAT FLUX
       SOLID_PHASE_OUTPUT = B1%Q_CONDENSE * 0.001_EB

Verification/WUI/char_oxidation_1.fds

@@ -79,6 +79,7 @@
 &VENT MB='ZMAX',SURF_ID='OPEN' /
 
 &DEVC XB=-0.5, 0.5,-0.5,0.5,0.0,3.0, QUANTITY='MPUV', PART_ID='foliage part', ID='mass foliage', SPATIAL_STATISTIC='VOLUME INTEGRAL'  /
+&DEVC XB=-0.5, 0.5,-0.5,0.5,0.0,3.0, QUANTITY='OXIDATIVE HRRPUA', PART_ID='foliage part', ID='char hrr', SPATIAL_STATISTIC='SURFACE INTEGRAL'  /
 
 &TAIL /