From 9420b1f0bce87d47f4dc0e52be3b0683095a02a9 Mon Sep 17 00:00:00 2001 From: mcgratta Date: Fri, 13 Dec 2024 16:51:21 -0500 Subject: [PATCH] FDS Source: #13885. Fix bug at Atmospheric interpolated boundary --- Source/wall.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/Source/wall.f90 b/Source/wall.f90 index 689fe356957..fb5493d2408 100644 --- a/Source/wall.f90 +++ b/Source/wall.f90 @@ -645,10 +645,11 @@ SUBROUTINE SURFACE_HEAT_TRANSFER(NM,T,SF,BC,B1,WALL_INDEX,CFACE_INDEX,PARTICLE_I ! interp or extrap RHO_OTHER for jump in vertical grid resolution, linear in temperature to match heat flux in divg PBAR_G = PBAR_P(BC%KKG,B1%PRESSURE_ZONE) PBAR_OTHER = EVALUATE_RAMP(MM%ZC(EWC%KKO_MIN),I_RAMP_P0_Z) - ZZ_GET(1:N_TRACKED_SPECIES) = MAX(0._EB,MIN(1._EB,ZZP(BC%II,BC%JJ,BC%KK,1:N_TOTAL_SCALARS)/RHOP(BC%II,BC%JJ,BC%KK))) + ZZ_GET(1:N_TRACKED_SPECIES) = MAX(0._EB,MIN(1._EB,ZZP(BC%II,BC%JJ,BC%KK,1:N_TOTAL_SCALARS))) CALL GET_SPECIFIC_GAS_CONSTANT(ZZ_GET,RSUM(BC%II,BC%JJ,BC%KK)) DENOM = PBAR_G/B1%RHO_G/RSUM_G + DDO*(PBAR_OTHER/RHO_OTHER/RSUM(BC%II,BC%JJ,BC%KK) - PBAR_G/B1%RHO_G/RSUM_G) RHOP(BC%II,BC%JJ,BC%KK) = PBAR_P(BC%KK,B1%PRESSURE_ZONE)/RSUM(BC%II,BC%JJ,BC%KK)/DENOM + TMP(BC%II,BC%JJ,BC%KK) = PBAR_P(BC%KK,B1%PRESSURE_ZONE)/(RSUM(BC%II,BC%JJ,BC%KK)*RHOP(BC%II,BC%JJ,BC%KK)) ENDIF ! Correction for variable molecular weights