Merge pull request #13958 from rmcdermo/master

FDS Source: add output for BLOWING_CORRECTION
firemodels · Dec 30, 2024 · e838729 · e838729
2 parents 566b2a2 + 2afddbd
commit e838729
Show file tree

Hide file tree

Showing 4 changed files with 13 additions and 3 deletions.
diff --git a/Source/data.f90 b/Source/data.f90
@@ -1612,6 +1612,10 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 OUTPUT_QUANTITY(-81)%BNDF_APPROPRIATE = .FALSE.
 OUTPUT_QUANTITY(-81)%PROF_APPROPRIATE = .TRUE.
 
+OUTPUT_QUANTITY(-82)%NAME = 'BLOWING CORRECTION'
+OUTPUT_QUANTITY(-82)%UNITS= ''
+OUTPUT_QUANTITY(-82)%SHORT_NAME = 'bcor'
+
 ! Fire spread
 OUTPUT_QUANTITY(-90)%NAME = 'FIRE ARRIVAL TIME'
 OUTPUT_QUANTITY(-90)%UNITS = 's'

diff --git a/Source/dump.f90 b/Source/dump.f90
@@ -9262,13 +9262,13 @@ REAL(EB) FUNCTION SOLID_PHASE_OUTPUT(INDX,Y_INDEX,Z_INDEX,PART_INDEX,OPT_WALL_IN
       SOLID_PHASE_OUTPUT = Q_CON*0.001_EB
 
    CASE(77) ! CONVECTIVE HEAT TRANSFER REGIME
-      SOLID_PHASE_OUTPUT = 0
+      SOLID_PHASE_OUTPUT = 0._EB
       IF (SF%INCLUDE_BOUNDARY_PROP2_TYPE) SOLID_PHASE_OUTPUT = B2%HEAT_TRANSFER_REGIME
    CASE(78) ! SURFACE OXYGEN MASS FRACTION
-      SOLID_PHASE_OUTPUT = 0
+      SOLID_PHASE_OUTPUT = 0._EB
       IF (SF%INCLUDE_BOUNDARY_PROP2_TYPE) SOLID_PHASE_OUTPUT = B2%Y_O2_F
    CASE(79) ! SURFACE OXYGEN ITERATIONS
-      SOLID_PHASE_OUTPUT = 0
+      SOLID_PHASE_OUTPUT = 0._EB
       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
@@ -9284,6 +9284,9 @@ REAL(EB) FUNCTION SOLID_PHASE_OUTPUT(INDX,Y_INDEX,Z_INDEX,PART_INDEX,OPT_WALL_IN
             SOLID_PHASE_OUTPUT = B2%Y_O2_F*EXP(-MAX(0._EB,DEPTH-CHAR_FRONT)/(TWO_EPSILON_EB+ML%GAS_DIFFUSION_DEPTH(1)))
          ENDIF
       ENDIF
+   CASE(82) ! BLOWING CORRECTION
+      SOLID_PHASE_OUTPUT = 0._EB
+      IF (SF%INCLUDE_BOUNDARY_PROP2_TYPE) SOLID_PHASE_OUTPUT = B2%BLOWING_CORRECTION
    CASE(90) ! FIRE ARRIVAL TIME
       IF (PRESENT(OPT_WALL_INDEX)) THEN
          OUTPUT_INDEX = OPT_WALL_INDEX

diff --git a/Source/type.f90 b/Source/type.f90
@@ -358,6 +358,7 @@ MODULE TYPES
    REAL(EB) :: K_SUPPRESSION=0._EB   !< Suppression coefficent (m2/kg/s)
    REAL(EB) :: V_DEP=0._EB           !< Deposition velocity (m/s)
    REAL(EB) :: Y_O2_F=0._EB          !< Oxygen mass fraction at the surface
+   REAL(EB) :: BLOWING_CORRECTION=0._EB !< Ackermann blowing correction to heat transfer coefficient
 
    INTEGER :: SURF_INDEX=-1          !< Surface index
    INTEGER :: HEAT_TRANSFER_REGIME=0 !< 1=Forced convection, 2=Natural convection, 3=Impact convection, 4=Resolved

diff --git a/Source/wall.f90 b/Source/wall.f90
@@ -3665,6 +3665,7 @@ REAL(EB) FUNCTION HEAT_TRANSFER_COEFFICIENT(NMX,DELTA_N_TMP,H_FIXED,SFX,WALL_IND
 
 IF (SFX%BLOWING .AND. .NOT. SFX%BOUNDARY_FUEL_MODEL .AND. SIM_MODE /= DNS_MODE .AND. ALLOCATED(P1X%M_DOT_G_PP_ACTUAL)) THEN
    PHI = 0._EB
+   IF (SFX%INCLUDE_BOUNDARY_PROP2_TYPE) P2X%BLOWING_CORRECTION=0._EB
    ITMP = INT(TMP_FILM)
    DO I=1,N_TRACKED_SPECIES
       IF (ABS(P1X%M_DOT_G_PP_ACTUAL(I)) <= TWO_EPSILON_EB) CYCLE
@@ -3676,6 +3677,7 @@ REAL(EB) FUNCTION HEAT_TRANSFER_COEFFICIENT(NMX,DELTA_N_TMP,H_FIXED,SFX,WALL_IND
          HEAT_TRANSFER_COEFFICIENT = 0._EB
       ELSE
          HEAT_TRANSFER_COEFFICIENT = HEAT_TRANSFER_COEFFICIENT * PHI/(EXP(PHI)-1._EB)
+         IF (SFX%INCLUDE_BOUNDARY_PROP2_TYPE) P2X%BLOWING_CORRECTION = PHI/(EXP(PHI)-1._EB)
       ENDIF
    ENDIF
 ENDIF