FDS Source: Use radiation-corrected TEMPERATURE for Smoke3D

mcgratta · Aug 10, 2023 · 38565bf · 38565bf
1 parent 32a4334
commit 38565bf
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/Source/dump.f90 b/Source/dump.f90
@@ -4104,6 +4104,7 @@ SUBROUTINE DUMP_SMOKE3D(T,DT,NM)
 REAL(FB) :: DXX,STIME
 REAL(EB), POINTER, DIMENSION(:,:,:) :: FF
 REAL(FB), ALLOCATABLE, DIMENSION(:) :: QQ_PACK
+REAL(EB) :: FR_C
 TYPE(SMOKE3D_TYPE), POINTER :: S3
 
 CALL POINT_TO_MESH(NM)
@@ -4130,6 +4131,13 @@ SUBROUTINE DUMP_SMOKE3D(T,DT,NM)
       ENDDO
    ENDDO
 
+   ! Adjust the temperature as it is used in the expression for the radiation source term
+
+   IF (S3%DISPLAY_TYPE=='TEMPERATURE' .AND. RTE_SOURCE_CORRECTION) THEN
+      FR_C = RTE_SOURCE_CORRECTION_FACTOR**0.25_EB
+      WHERE (CHI_R*Q>QR_CLIP) FF = (FF+TMPM)*FR_C - TMPM
+   ENDIF
+
    ! Interpolate data to cell nodes
 
    DO K=0,KBAR