Merge pull request #12496 from drjfloyd/master

FDS Source: Issue #12481, fix Q_RAD_IN adjustment in part
firemodels · Feb 16, 2024 · 0edfabc · 0edfabc
2 parents 9f52309 + 4dac580
commit 0edfabc
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/Source/part.f90 b/Source/part.f90
@@ -3926,7 +3926,7 @@ SUBROUTINE PARTICLE_MASS_ENERGY_TRANSFER(T,DT,NM)
          B2%LP_CPUA(LPC%ARRAY_INDEX) = B2%LP_CPUA(LPC%ARRAY_INDEX) + &
                                        (1._EB-LPC%RUNNING_AVERAGE_FACTOR_WALL)*WGT*Q_CON_WALL/(B1%AREA*DT)
          IF (RADIATION) &
-         B1%Q_RAD_IN = (B1%AREA*DT*B1%Q_RAD_IN-WGT*DT*Q_DOT_RAD*B1%Q_RAD_IN/(B1%Q_RAD_IN+B1%Q_RAD_OUT+TWO_EPSILON_EB))/(B1%AREA*DT)
+         B1%Q_RAD_IN = MAX(0.,(B1%AREA*DT*B1%Q_RAD_IN-WGT*DT*Q_DOT_RAD*B1%Q_RAD_IN/(B1%Q_RAD_IN+B1%Q_RAD_OUT+TWO_EPSILON_EB))/(B1%AREA*DT))
       ENDIF
 
    ENDDO PARTICLE_LOOP