diff --git a/Source/data.f90 b/Source/data.f90 index db2064946b0..2c5b823bbef 100644 --- 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 index b9e88e2db86..6da3385495c 100644 --- 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 index cd71c47d489..c52a5606ac3 100644 --- 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 index 7c2f34ab15b..91355aa347b 100644 --- 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