Merge remote-tracking branch 'firemodels/master' into FireX

Author: cxp484

Manuals/Bibliography/FDS_general.bib

@@ -1838,7 +1838,7 @@ @TECHREPORT{Farouki:1981
 
 @PHDTHESIS{Fernandez:2018,
   author       = {Sebastian Ferreyro Fernandez},
-  title        = {ADVANCED SOOT AND RADIATION MODELS FOR LAMINAR AND TURBULENT FLAMES},
+  title        = {Advanced Soot And Radiation Models For Laminar And Turbulent Flames},
   year         = {2018},
   school       = {The Pennsylvania State University}
 }

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -2271,7 +2271,7 @@ \subsubsection{Default Convective Heat Transfer Model}
 \end{center}
 \end{table}
 
-The length scale, $L$, is specified by \ct{CONVECTION_LENGTH_SCALE} on the \ct{SURF} line. By default, it is 1~m for plates, and the diameter of a sphere or cylinder.
+The length scale, $L$, is specified by \ct{CONVECTION_LENGTH_SCALE} on the \ct{SURF} line. By default, it is 1~m for plates, and the diameter of a sphere or cylinder. The velocity used to compute Re can be output with the solid phase quantity \ct{'TANGENTIAL VELOCITY'}.
 
 \subsubsection{Output for Convective Heat Transfer Regime}
 \label{info:convection_regime}
@@ -11741,6 +11741,7 @@ \section{Solid Phase Output Quantities}
 \ct{SUBSTEPS}                                  & Section~\ref{info:solid_phase_stability}      &                    & B,D          \\ \hline
 \ct{SURFACE DENSITY}$^4$                       & Section~\ref{info:material_components}        & kg/m$^2$           & B,D          \\ \hline
 \ct{SURFACE DEPOSITION}$^1$                    & Section~\ref{info:deposition}                 & kg/m$^2$           & B,D          \\ \hline
+\ct{TANGENTIAL VELOCITY}                       & Section~\ref{info:convection}                 & m/s                & B,D          \\ \hline
 \ct{TOTAL MASS FLUX WALL}$^1$                  & Section~\ref{info:wallflux}                   & kg/s/m$^2$         & B,D          \\ \hline
 \ct{VELOCITY ERROR}                            & Section~\ref{info:PRES}                       & m/s                & B,D          \\ \hline
 \ct{VISCOUS STRESS WALL}                       & Section~\ref{info:distributed_forces}         & Pa                 & B,D          \\ \hline

…de/SCRIPT_FIGURES/Sandia_Jets/.gitignore …S/Sandia_Jets_Pools_Fireballs/.gitignoreManuals/FDS_Validation_Guide/SCRIPT_FIGURES/Sandia_Jets/.gitignore renamed to Manuals/FDS_Validation_Guide/SCRIPT_FIGURES/Sandia_Jets_Pools_Fireballs/.gitignore Manuals/FDS_Validation_Guide/SCRIPT_FIGURES/Sandia_Jets/.gitignore renamed to Manuals/FDS_Validation_Guide/SCRIPT_FIGURES/Sandia_Jets_Pools_Fireballs/.gitignore

@@ -1752,6 +1752,10 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 OUTPUT_QUANTITY(-100)%SPEC_ID_REQUIRED = .FALSE.
 OUTPUT_QUANTITY(-100)%BNDF_APPROPRIATE = .TRUE.
 
+OUTPUT_QUANTITY(-101)%NAME = 'TANGENTIAL VELOCITY'
+OUTPUT_QUANTITY(-101)%UNITS= 'm/s'
+OUTPUT_QUANTITY(-101)%SHORT_NAME = 'u_tang'
+
 END SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 
 

Source/data.f90

@@ -10703,6 +10703,9 @@ REAL(EB) FUNCTION SOLID_PHASE_OUTPUT(INDX,Y_INDEX,Z_INDEX,PART_INDEX,OPT_WALL_IN
    CASE(100) ! CONDENSATION HEAT FLUX
       SOLID_PHASE_OUTPUT = B1%Q_CONDENSE * 0.001_EB
 
+   CASE(101) ! TANGENTIAL VELOCITY
+      IF (ASSOCIATED(B1)) SOLID_PHASE_OUTPUT = B1%U_TANG
+
 END SELECT SOLID_PHASE_SELECT
 
 END FUNCTION SOLID_PHASE_OUTPUT

Source/dump.f90

@@ -1188,73 +1188,83 @@ SUBROUTINE CALC_AFT_REAC_AND_PROD(ZZ,ZZ_REAC,ZZ_PROD)
 
 REAL(EB), INTENT(IN) :: ZZ(N_TRACKED_SPECIES)
 REAL(EB), INTENT(OUT) :: ZZ_REAC(N_TRACKED_SPECIES),ZZ_PROD(N_TRACKED_SPECIES)
+REAL(EB) :: ZZ_EQUIV(N_TRACKED_SPECIES)
 REAL(EB) :: EQUIV, X,Y,Z,A,B,C,D,E, SUM_ZZ
 REAL(EB) :: ZZ_MASS, ZZ_FUEL_O2_MASS
 REAL(EB) :: REAC_MASS, PROD_MASS, FUEL_MASS, O2_MASS, CO2_MASS, H2O_MASS
 INTEGER :: NS
-
-CALL CALC_EQUIV_RATIO(ZZ(1:N_TRACKED_SPECIES), EQUIV)
-IF (SIM_MODE /= DNS_MODE .AND. USE_MIXED_ZN_AFT_TMP) THEN
-   EQUIV = MAX(1._EB,EQUIV)
-ENDIF
+REAL(EB), PARAMETER:: EQUIV_TOL=1.E-4_EB
 
 ! Adjust for diluents
 ZZ_FUEL_O2_MASS = 0._EB
 ZZ_MASS = SUM(ZZ(1:N_TRACKED_SPECIES))
+ZZ_EQUIV = 0.0_EB
 ZZ_REAC(1:N_TRACKED_SPECIES) = ZZ(1:N_TRACKED_SPECIES)
 ZZ_PROD(1:N_TRACKED_SPECIES) = ZZ(1:N_TRACKED_SPECIES)
 DO NS=1, N_TRACKED_SPECIES
-   IF (SPECIES_MIXTURE(NS)%OXR > TWO_EPSILON_EB) THEN ! FUEL
+   IF (SPECIES_MIXTURE(NS)%OXR > TWO_EPSILON_EB .AND. &
+          ABS(SPECIES_MIXTURE(NS)%EQUIV - 1.0_EB) > EQUIV_TOL) THEN ! FUEL
       ZZ_REAC(NS) = 0._EB
       ZZ_PROD(NS) = 0._EB
       ZZ_FUEL_O2_MASS = ZZ_FUEL_O2_MASS + ZZ(NS)
+      ZZ_EQUIV(NS) = ZZ(NS)
    ENDIF   
 ENDDO
 ZZ_REAC(I_O2) = 0._EB
 ZZ_PROD(I_O2) = 0._EB
+ZZ_EQUIV(I_O2) = ZZ(I_O2)
 ZZ_FUEL_O2_MASS = ZZ_FUEL_O2_MASS + ZZ(I_O2)
 
-! Based on CxHyOz + aO2 = bCO2 + cH2O + dCxHyOz + eO2 
-X=SPECIES_MIXTURE(I_FUEL)%ATOMS(6) !C
-Y=SPECIES_MIXTURE(I_FUEL)%ATOMS(1) !H
-Z=SPECIES_MIXTURE(I_FUEL)%ATOMS(8) !O
-A=0.5_EB*((2._EB*X+0.5_EB*Y)/EQUIV -Z) ! a is a function of equivalence ratio
-IF (ABS(EQUIV - 1.0_EB) < TWO_EPSILON_EB) THEN ! Stoich
-   B = X
-   C = 0.5_EB*Y
-   D = 0._EB ! No fuel
-   E = 0._EB ! No O2
-ELSEIF (EQUIV > 1) THEN ! Rich
-   D = (2._EB*X+0.5_EB*Y-2._EB*A-Z)/(2._EB*X+0.5_EB*Y-Z)
-   B = (1._EB-D)*X
-   C = 0.5_EB*(1._EB-D)*Y
-   E = 0._EB ! No O2
-ELSE !EQUIV < 1, Lean
-   B = X
-   C = 0.5_EB*Y
-   D = 0._EB ! No fuel
-   E = A + 0.5*Z - B - 0.5_EB*C
-ENDIF
+IF (ZZ_FUEL_O2_MASS > TWO_EPSILON_EB) THEN
+   SUM_ZZ=SUM(ZZ_EQUIV)
+   ZZ_EQUIV = ZZ_EQUIV/(SUM_ZZ + TWO_EPSILON_EB)
+   CALL CALC_EQUIV_RATIO(ZZ_EQUIV(1:N_TRACKED_SPECIES), EQUIV)
+   IF (SIM_MODE /= DNS_MODE .AND. USE_MIXED_ZN_AFT_TMP) THEN
+      EQUIV = MAX(1._EB,EQUIV)
+   ENDIF
 
-! Setup reactants
-FUEL_MASS = 1._EB*SPECIES_MIXTURE(I_FUEL)%MW
-O2_MASS = A*SPECIES_MIXTURE(I_O2)%MW
-REAC_MASS = FUEL_MASS+O2_MASS
-ZZ_REAC(I_FUEL)=FUEL_MASS/REAC_MASS*ZZ_FUEL_O2_MASS
-ZZ_REAC(I_O2)=O2_MASS/REAC_MASS*ZZ_FUEL_O2_MASS
-SUM_ZZ=SUM(ZZ_REAC)
-ZZ_REAC = ZZ_REAC/SUM_ZZ
-
-! Setup products
-CO2_MASS = B*SPECIES_MIXTURE(I_CO2)%MW
-H2O_MASS = C*SPECIES_MIXTURE(I_H2O)%MW
-FUEL_MASS = D*SPECIES_MIXTURE(I_FUEL)%MW
-O2_MASS = E*SPECIES_MIXTURE(I_O2)%MW
-PROD_MASS = CO2_MASS+H2O_MASS+FUEL_MASS+O2_MASS
-ZZ_PROD(I_CO2)=ZZ_PROD(I_CO2)+CO2_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
-ZZ_PROD(I_H2O)=ZZ_PROD(I_H2O)+H2O_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
-ZZ_PROD(I_FUEL)=FUEL_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
-ZZ_PROD(I_O2)=O2_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
+   ! Based on CxHyOz + aO2 = bCO2 + cH2O + dCxHyOz + eO2 
+   X=SPECIES_MIXTURE(I_FUEL)%ATOMS(6) !C
+   Y=SPECIES_MIXTURE(I_FUEL)%ATOMS(1) !H
+   Z=SPECIES_MIXTURE(I_FUEL)%ATOMS(8) !O
+   A=0.5_EB*((2._EB*X+0.5_EB*Y)/EQUIV -Z) ! a is a function of equivalence ratio
+   IF (ABS(EQUIV - 1.0_EB) < EQUIV_TOL) THEN ! Stoich
+      B = X
+      C = 0.5_EB*Y
+      D = 0._EB ! No fuel
+      E = 0._EB ! No O2
+   ELSEIF (EQUIV > 1) THEN ! Rich
+      D = (2._EB*X+0.5_EB*Y-2._EB*A-Z)/(2._EB*X+0.5_EB*Y-Z)
+      B = (1._EB-D)*X
+      C = 0.5_EB*(1._EB-D)*Y
+      E = 0._EB ! No O2
+   ELSE !EQUIV < 1, Lean
+      B = X
+      C = 0.5_EB*Y
+      D = 0._EB ! No fuel
+      E = A + 0.5*Z - B - 0.5_EB*C
+   ENDIF
+
+   ! Setup reactants
+   FUEL_MASS = 1._EB*SPECIES_MIXTURE(I_FUEL)%MW
+   O2_MASS = A*SPECIES_MIXTURE(I_O2)%MW
+   REAC_MASS = FUEL_MASS+O2_MASS
+   ZZ_REAC(I_FUEL)=FUEL_MASS/REAC_MASS*ZZ_FUEL_O2_MASS
+   ZZ_REAC(I_O2)=O2_MASS/REAC_MASS*ZZ_FUEL_O2_MASS
+   SUM_ZZ=SUM(ZZ_REAC)
+   ZZ_REAC = ZZ_REAC/SUM_ZZ
+
+   ! Setup products
+   CO2_MASS = B*SPECIES_MIXTURE(I_CO2)%MW
+   H2O_MASS = C*SPECIES_MIXTURE(I_H2O)%MW
+   FUEL_MASS = D*SPECIES_MIXTURE(I_FUEL)%MW
+   O2_MASS = E*SPECIES_MIXTURE(I_O2)%MW
+   PROD_MASS = CO2_MASS+H2O_MASS+FUEL_MASS+O2_MASS
+   ZZ_PROD(I_CO2)=ZZ_PROD(I_CO2)+CO2_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
+   ZZ_PROD(I_H2O)=ZZ_PROD(I_H2O)+H2O_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
+   ZZ_PROD(I_FUEL)=FUEL_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
+   ZZ_PROD(I_O2)=O2_MASS/PROD_MASS*ZZ_FUEL_O2_MASS
+ENDIF
 SUM_ZZ=SUM(ZZ_PROD)
 ZZ_PROD = ZZ_PROD/SUM_ZZ
 

Source/fire.f90

@@ -3572,6 +3572,7 @@ SUBROUTINE READ_SPEC
                                 * ELEMENT(8)%MASS / SPECIES_MIXTURE(N1)%MW
 
    SPECIES_MIXTURE(N1)%OXA = SPECIES_MIXTURE(N1)%ATOMS(8) * ELEMENT(8)%MASS / SPECIES_MIXTURE(N1)%MW
+   SPECIES_MIXTURE(N1)%EQUIV = SPECIES_MIXTURE(N1)%OXR/(SPECIES_MIXTURE(N1)%OXA+TWO_EPSILON_EB)
 ENDDO
 
 !If SOOT_OXIDATION is enabled make sure to save accumulation on the wall and check that it is an AEROSOL species

Source/read.f90

@@ -650,6 +650,7 @@ MODULE TYPES
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: R50
    REAL(EB) :: OXR                  !< Required oxygen for complete combustion (gm/gm-species)
    REAL(EB) :: OXA                  !< Available oxygen for combustion (gm/gm-species)
+   REAL(EB) :: EQUIV                !< Species equiv ratio OXR/OXA
 
 
 END TYPE SPECIES_MIXTURE_TYPE