Merge pull request #12356 from drjfloyd/master

FDS Source: prevent HTC overflow with BLOWING
firemodels · Jan 10, 2024 · c886744 · c886744
2 parents 525e92e + 14bb7ee
commit c886744
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/Source/wall.f90 b/Source/wall.f90
@@ -3268,7 +3268,11 @@ REAL(EB) FUNCTION HEAT_TRANSFER_COEFFICIENT(NMX,DELTA_N_TMP,H_FIXED,SURF_INDEX_I
    ENDDO
    IF (ABS(PHI)>TWO_EPSILON_EB .AND. ABS(HEAT_TRANSFER_COEFFICIENT)>=TWO_EPSILON_EB) THEN
       PHI = PHI / HEAT_TRANSFER_COEFFICIENT
-      HEAT_TRANSFER_COEFFICIENT = HEAT_TRANSFER_COEFFICIENT * PHI/(EXP(PHI)-1._EB)
+      IF (PHI > 100._EB) THEN
+         HEAT_TRANSFER_COEFFICIENT = 0._EB
+      ELSE
+         HEAT_TRANSFER_COEFFICIENT = HEAT_TRANSFER_COEFFICIENT * PHI/(EXP(PHI)-1._EB)
+      ENDIF
    ENDIF
 ENDIF