Merge pull request #14618 from cxp484/FireX

FireX FDS Source: Add Thickened flame model for laminar case

Author: cxp484

Source/cons.f90

@@ -955,4 +955,7 @@ MODULE CHEMCONS
 REAL(EB), ALLOCATABLE, DIMENSION(:,:) :: CVODE_SUBSTEP_DATA
 INTEGER :: TOTAL_SUBSTEPS_TAKEN
 
+! FOR THICKENED FLAME MODEL
+REAL(EB) :: FLAME_THICK_FACTOR=1.0_EB
+
 END MODULE CHEMCONS

Source/divg.f90

@@ -30,6 +30,7 @@ SUBROUTINE DIVERGENCE_PART_1(T,DT,NM)
 USE COMPLEX_GEOMETRY, ONLY : CC_CGSC, CC_UNKZ, CC_SOLID, CC_CUTCFE
 USE CC_SCALARS, ONLY : ADD_CUTCELL_PSUM,ADD_LINKEDCELL_PSUM,SET_EXIMDIFFLX_3D,SET_EXIMRHOHSLIM_3D,&
                        SET_EXIMRHOZZLIM_3D,CC_DIVERGENCE_PART_1,CC_VELOCITY_FLUX,CFACE_PREDICT_NORMAL_VELOCITY
+USE CHEMCONS, ONLY : FLAME_THICK_FACTOR
 
 INTEGER, INTENT(IN) :: NM
 REAL(EB), INTENT(IN) :: T,DT
@@ -149,6 +150,10 @@ SUBROUTINE DIVERGENCE_PART_1(T,DT,NM)
          ENDIF
       ENDIF
 
+      IF (FLAME_THICK_FACTOR > 1._EB) THEN
+         RHO_D = RHO_D*FLAME_THICK_FACTOR
+      ENDIF
+
       ! Manufactured solution
 
       IF (PERIODIC_TEST==7) RHO_D = DIFF_MMS

Source/fire.f90

@@ -413,6 +413,8 @@ SUBROUTINE SET_SPECIES_SOURCE_TERM_CELL(DT, I, J, K, ZZ_NEW, ZZ_OLD, REAC_SOURCE
 
 USE PHYSICAL_FUNCTIONS, ONLY: GET_SPECIFIC_GAS_CONSTANT,GET_SPECIFIC_HEAT, &
                               GET_SENSIBLE_ENTHALPY_Z,IS_REALIZABLE
+USE CHEMCONS, ONLY : FLAME_THICK_FACTOR
+
 REAL(EB), INTENT(IN) :: DT
 REAL(EB), INTENT(INOUT) :: ZZ_NEW(N_TRACKED_SPECIES), ZZ_OLD(N_TRACKED_SPECIES)
 REAL(EB), INTENT(IN) :: REAC_SOURCE_TERM_TMP(N_TRACKED_SPECIES), Q_REAC_TMP(N_REACTIONS)
@@ -436,7 +438,12 @@ SUBROUTINE SET_SPECIES_SOURCE_TERM_CELL(DT, I, J, K, ZZ_NEW, ZZ_OLD, REAC_SOURCE
       RETURN
    ENDIF
 ENDIF
+
 DZZ = ZZ_NEW - ZZ_OLD
+IF (FLAME_THICK_FACTOR > 1._EB) THEN
+   DZZ = DZZ/FLAME_THICK_FACTOR
+ENDIF
+
 ! Update RSUM and ZZ
 DZZ_IF: IF ( ANY(ABS(DZZ) > DZZ_CLIP) ) THEN
    IF (ABS(Q(I,J,K)) > TWO_EPSILON_EB) Q_EXISTS = .TRUE.

Source/read.f90

@@ -4604,10 +4604,10 @@ END SUBROUTINE PROC_SPEC_2
 !> \brief Read the COMBustion namelist line
 
 SUBROUTINE READ_COMB
-USE CHEMCONS, ONLY: ODE_MIN_ATOL, EQUIV_RATIO_CHECK, MIN_EQUIV_RATIO, MAX_EQUIV_RATIO, DO_CHEM_LOAD_BALANCE
+USE CHEMCONS, ONLY: ODE_MIN_ATOL, EQUIV_RATIO_CHECK, MIN_EQUIV_RATIO, MAX_EQUIV_RATIO, DO_CHEM_LOAD_BALANCE,FLAME_THICK_FACTOR
 REAL(EB) :: ODE_REL_ERROR
 NAMELIST /COMB/ CHECK_REALIZABILITY,COMPUTE_ADIABATIC_FLAME_TEMPERATURE, DO_CHEM_LOAD_BALANCE, EQUIV_RATIO_CHECK, &
-                EXTINCTION_MODEL,FINITE_RATE_MIN_TEMP, FIXED_MIX_TIME,FREE_BURN_TEMPERATURE,&
+                EXTINCTION_MODEL,FINITE_RATE_MIN_TEMP, FIXED_MIX_TIME,FLAME_THICK_FACTOR,FREE_BURN_TEMPERATURE,&
                 INITIAL_UNMIXED_FRACTION, MAX_CHEMISTRY_SUBSTEPS, MAX_EQUIV_RATIO, MIN_EQUIV_RATIO, &
                 N_FIXED_CHEMISTRY_SUBSTEPS, ODE_MIN_ATOL,ODE_REL_ERROR,ODE_SOLVER,SUPPRESSION,TAU_CHEM, &
                 TAU_FLAME,ZZ_MIN_GLOBAL

Source/wall.f90

@@ -435,6 +435,7 @@ SUBROUTINE SURFACE_HEAT_TRANSFER(NM,T,SF,BC,B1,WALL_INDEX,CFACE_INDEX,PARTICLE_I
 USE COMPLEX_GEOMETRY, ONLY : CC_CGSC, CC_SOLID
 USE PHYSICAL_FUNCTIONS, ONLY : GET_SPECIFIC_GAS_CONSTANT,GET_VISCOSITY,GET_MOLECULAR_WEIGHT
 USE DEVICE_VARIABLES, ONLY : PROPERTY,PROPERTY_TYPE
+USE CHEMCONS, ONLY : FLAME_THICK_FACTOR
 
 REAL(EB), INTENT(IN) :: T
 INTEGER, INTENT(IN) :: NM
@@ -942,6 +943,9 @@ SUBROUTINE SURFACE_HEAT_TRANSFER(NM,T,SF,BC,B1,WALL_INDEX,CFACE_INDEX,PARTICLE_I
                               ENDIF
                               RHO_D = 0.5_EB*( RHO_OTHER*D_Z_OTHER + B1%RHO_G*D_Z_G ) + RHO_D_TURB
                         END SELECT MODE_SELECT
+                        IF (FLAME_THICK_FACTOR > 1._EB) THEN
+                           RHO_D = RHO_D*FLAME_THICK_FACTOR
+                        ENDIF
                         SELECT CASE(BC%IOR)
                            CASE( 1)
                               ARO = MM%DY(JJO)*MM%DZ(KKO)/(DY(BC%JJ)*DZ(BC%KK))
@@ -1011,6 +1015,7 @@ END SUBROUTINE SURFACE_HEAT_TRANSFER
 !> \details Calculate the diffusion coefficient, rho*D (kg/m/s)
 
 SUBROUTINE CALCULATE_RHO_D_F(B1,BC,WALL_INDEX,CFACE_INDEX)
+USE CHEMCONS, ONLY : FLAME_THICK_FACTOR
 
 INTEGER, INTENT(IN), OPTIONAL :: WALL_INDEX,CFACE_INDEX
 REAL(EB) :: MU_G
@@ -1044,7 +1049,9 @@ SUBROUTINE CALCULATE_RHO_D_F(B1,BC,WALL_INDEX,CFACE_INDEX)
          B1%RHO_D_F(N) = B1%RHO_F*D_Z(ITMP,N)
       ENDDO
 END SELECT
-
+IF (FLAME_THICK_FACTOR > 1._EB) THEN
+   B1%RHO_D_F = B1%RHO_D_F*FLAME_THICK_FACTOR
+ENDIF
 END SUBROUTINE CALCULATE_RHO_D_F