Skip to content

Commit c9655cc

Browse files
mcgrattamcgratta
authored
Merge pull request #14096 from mcgratta/master
FDS Source: Reduce QR_CLIP
2 parents 9d4904d + 24cff7b commit c9655cc

File tree

2 files changed

+24
-17
lines changed

2 files changed

+24
-17
lines changed

‎Manuals/FDS_User_Guide/FDS_User_Guide.tex‎

Lines changed: 23 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -2875,6 +2875,12 @@ \subsection{Solid Phase Gas Transport}
28752875

28762876
Suppose, for example, that the solid is composed of a layer of insulation on top of a layer of plastic on top of a layer of steel. The steel is impermeable. By setting \ct{LAYER_DIVIDE=2.0} to the \ct{SURF} line whose first layer is the insulation directs the vapors generated by the insulation and plastic to be driven out of the exterior surface of the insulation. Similarly, for the \ct{SURF} line that is applied to the steel, specifying \ct{LAYER_DIVIDE=0.0} would indicate that no fuel vapors are to escape the steel surface. Note that in this instance, the sum of the values of \ct{LAYER_DIVIDE} is not equal to the number of layers, but this is not a problem because the layer of steel does not generate any gases.
28772877

2878+
The solid phase output \ct{QUANTITY='PYROLYSIS DEPTH'} can check if the \ct{LAYER_DIVIDE} is being applied correctly. For example, the input line
2879+
\begin{lstlisting}
2880+
&BNDF QUANTITY='PYROLYSIS DEPTH', CELL_CENTERED=T /
2881+
\end{lstlisting}
2882+
shows the distance from the surface where the pyrolyates are directed towards that surface, as opposed to the opposite side of the solid.
2883+
28782884

28792885

28802886
\subsection{Reaction Rates}
@@ -11465,45 +11471,46 @@ \section{Solid Phase Output Quantities}
1146511471
\ct{BACK WALL TEMPERATURE} & Section~\ref{info:BACK} & $^\circ$C & B,D \\ \hline
1146611472
\ct{BLOWING CORRECTION} & Section~\ref{info:blowing} & & B,D \\ \hline
1146711473
\ct{BURNING RATE} & Mass loss rate of fuel & \unit{kg/(m^2.s)} & B,D \\ \hline
11468-
\ct{CONDENSATION HEAT FLUX} & Section~\ref{info:condensation} & \unit{kW/m^2} & B,D \\ \hline
11469-
\ct{CONVECTIVE HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11470-
\ct{CONVECTIVE HEAT FLUX GAUGE} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11471-
{\scriptsize\tt CONVECTIVE HEAT TRANSFER REGIME} & Section \ref{info:convection} & & B,D \\ \hline
11472-
\ct{CPUA}$^2$ & Section~\ref{bucket_test_1} & \unit{kW/m^2} & B,D \\ \hline
11473-
\ct{CPUA_Z}$^1$ & Section~\ref{bucket_test_1} & \unit{kW/m^2} & B,D \\ \hline
11474+
\ct{CONDENSATION HEAT FLUX} & Section~\ref{info:condensation} & \unit{kW/m^2} & B,D \\ \hline
11475+
\ct{CONVECTIVE HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11476+
\ct{CONVECTIVE HEAT FLUX GAUGE} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11477+
\ct{CONVECTIVE HEAT TRANSFER REGIME} & Section \ref{info:convection} & & B,D \\ \hline
11478+
\ct{CPUA}$^2$ & Section~\ref{bucket_test_1} & \unit{kW/m^2} & B,D \\ \hline
11479+
\ct{CPUA_Z}$^1$ & Section~\ref{bucket_test_1} & \unit{kW/m^2} & B,D \\ \hline
1147411480
\ct{DEPOSITION VELOCITY} & Section~\ref{info:deposition} & m/s & B,D \\ \hline
1147511481
\ct{FRICTION VELOCITY} & Section~\ref{info:yplus} & m/s & B,D \\ \hline
11476-
\ct{GAUGE HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11477-
\ct{ENTHALPY FLUX WALL} & Section~\ref{info:enthalpy_flux} & \unit{kW/m^2} & B,D \\ \hline
11478-
\ct{TOTAL HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11482+
\ct{GAUGE HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11483+
\ct{ENTHALPY FLUX WALL} & Section~\ref{info:enthalpy_flux} & \unit{kW/m^2} & B,D \\ \hline
11484+
\ct{TOTAL HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
1147911485
\ct{EMISSIVITY} & Surface emissivity (usually constant) & & B,D \\ \hline
1148011486
\ct{FIRE ARRIVAL TIME} & Section \ref{info:fire_spread_output} & \si{s} & B,D \\ \hline
1148111487
\ct{FIRE RESIDENCE TIME} & Section \ref{info:fire_spread_output} & \si{s} & B,D \\ \hline
1148211488
\ct{GAS DENSITY} & Gas Density near wall & \si{kg/m^3} & B,D \\ \hline
1148311489
\ct{GAS TEMPERATURE} & Gas Temperature near wall & $^\circ$C & B,D \\ \hline
1148411490
\ct{HEAT TRANSFER COEFFICIENT} & Section \ref{info:convection} & \si{W/(m^2.K)} & B,D \\ \hline
11485-
\ct{HRRPUA} & $\dq''$ & \unit{kW/m^2} & B,D \\ \hline
11486-
\ct{INCIDENT HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11491+
\ct{HRRPUA} & $\dq''$ & \unit{kW/m^2} & B,D \\ \hline
11492+
\ct{INCIDENT HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
1148711493
\ct{INSIDE WALL TEMPERATURE} & Section~\ref{info:DEPTH} & $^\circ$C & D,Pr \\ \hline
1148811494
\ct{INSIDE WALL DEPTH} & Section~\ref{info:DEPTH} & m & D,Pr \\ \hline
1148911495
\ct{MASS FLUX}$^{1,4}$ & Section~\ref{info:wallflux} & \unit{kg/(m^2.s)} & B,D \\ \hline
1149011496
\ct{MASS FLUX WALL}$^1$ & Section~\ref{info:wallflux} & \unit{kg/(m^2.s)} & B,D \\ \hline
1149111497
\ct{MPUA}$^2$ & Section~\ref{bucket_test_1} & kg/m$^2$ & B,D \\ \hline
11492-
\ct{MPUA_Z}$^1$ & Section~\ref{bucket_test_1} & kg/m$^2$ & B,D \\ \hline
11498+
\ct{MPUA_Z}$^1$ & Section~\ref{bucket_test_1} & kg/m$^2$ & B,D \\ \hline
1149311499
\ct{NORMAL VELOCITY} & Wall normal velocity & m/s & B,D \\ \hline
1149411500
\ct{NORMALIZED HEATING RATE} & Section~\ref{info:material_components} & W/g & D \\ \hline
1149511501
\ct{NORMALIZED HEAT RELEASE RATE} & Section~\ref{info:material_components} & W/g & D \\ \hline
1149611502
\ct{NORMALIZED MASS}$^4$ & Section~\ref{info:material_components} & & D \\ \hline
1149711503
\ct{NORMALIZED MASS LOSS RATE}$^4$ & Section~\ref{info:material_components} & 1/s & D \\ \hline
1149811504
\ct{OXIDATIVE HRRPUA} & Section~\ref{sec:reaction_rates} & \unit{kW/m^2} & B,D \\ \hline
1149911505
\ct{PRESSURE COEFFICIENT} & Section~\ref{info:pressure_coefficient} & & B,D \\ \hline
11500-
\ct{RADIATIVE HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11501-
\ct{RADIOMETER} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11502-
\ct{REFERENCE_HEAT_FLUX} & Section~\ref{info:scaled_burning} & \unit{kW/m^2} & B,D \\ \hline
11506+
\ct{PYROLYSIS DEPTH} & Section~\ref{info:LAYER_DIVIDE} & m & B,D \\ \hline
11507+
\ct{RADIATIVE HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11508+
\ct{RADIOMETER} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
11509+
\ct{REFERENCE_HEAT_FLUX} & Section~\ref{info:scaled_burning} & \unit{kW/m^2} & B,D \\ \hline
1150311510
\ct{SOLID CONDUCTIVITY}$^4$ & Section~\ref{info:DEPTH} & \si{W/(m.K)} & D,Pr \\ \hline
1150411511
\ct{SOLID DENSITY}$^4$ & Section~\ref{info:DEPTH} & kg/m$^3$ & D,Pr \\ \hline
1150511512
\ct{SOLID ENTHALPY}$^4$ & Section~\ref{info:DEPTH} & kJ/m$^3$ & D,Pr \\ \hline
11506-
\ct{SOLID MASS FRACTION}$^5$ & Section~\ref{info:DEPTH} & \si{kg/kg} & D,Pr \\ \hline
11513+
\ct{SOLID MASS FRACTION}$^5$ & Section~\ref{info:DEPTH} & \si{kg/kg} & D,Pr \\ \hline
1150711514
\ct{SOLID SPECIFIC HEAT}$^4$ & Section~\ref{info:DEPTH} & \si{kJ/(kg.K)} & D,Pr \\ \hline
1150811515
\ct{SUBSTEPS} & Section~\ref{info:solid_phase_stability} & & B,D \\ \hline
1150911516
\ct{SURFACE DENSITY}$^4$ & Section~\ref{info:material_components} & kg/m$^2$ & B,D \\ \hline

‎Source/cons.f90‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -471,7 +471,7 @@ MODULE GLOBAL_CONSTANTS
471471
REAL(EB) :: RTE_SOURCE_CORRECTION_FACTOR=1._EB !< Multiplicative factor used in correcting RTE source term
472472
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$
473473
REAL(EB) :: KFST4_SUM=0._EB !< \f$ \sum_{ijk} 4 \kappa_{ijk} \sigma T_{ijk}^4 V_{ijk} \f$
474-
REAL(EB) :: QR_CLIP=10._EB !< Lower bound of \f$ \chi_{\rm r} \dot{q}_{ijk}''' \f$ below which no source correction is made
474+
REAL(EB) :: QR_CLIP=1._EB !< Lower bound of \f$ \chi_{\rm r} \dot{q}_{ijk}''' \f$ below which no source correction is made
475475
REAL(EB) :: C_MAX=100._EB !< Maximum value of RAD_Q_SUM/KFST4_SUM
476476
REAL(EB) :: C_MIN=0.1_EB !< Minimum value of RAD_Q_SUM/KFST4_SUM
477477

0 commit comments

Comments
 (0)