From ca0b6be8d99ef000079236b1f21d62b2dda9a394 Mon Sep 17 00:00:00 2001 From: Jason Floyd Date: Sun, 14 Dec 2025 09:45:21 -0500 Subject: [PATCH] FDS Source: Iterate BI-DIR for dynamic calibration --- Source/dump.f90 | 14 ++++++++++---- 1 file changed, 10 insertions(+), 4 deletions(-) diff --git a/Source/dump.f90 b/Source/dump.f90 index c3d4d40cc0..9de19a177d 100644 --- a/Source/dump.f90 +++ b/Source/dump.f90 @@ -7926,10 +7926,16 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z ! LJ AIR viscosity fit 100 K to 5000 K MU_G = 1.5205E-22_EB*PROBE_TMP**5 - 2.1417E-18_EB*PROBE_TMP**4 + 1.1402E-14_EB*PROBE_TMP**3 - & 2.9846E-11_EB*PROBE_TMP**2 + 5.9898E-8_EB*PROBE_TMP + 0.000002352_EB - RE_D = MIN(3800._EB,MAX(40._EB,PROBE_RHO*VEL*PY%PROBE_DIAMETER/MU_G)) - FAC = 1.533_EB-0.001366_EB*RE_D+0.000001688_EB*RE_D**2-0.0000000009706_EB*RE_D**3+& - 0.0000000000002555_EB*RE_D**4-2.484E-17_EB*RE_D**5 - GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*1._EB/FAC*SQRT(2*PROBE_DELTA_P/PROBE_RHO) + JJJ = 1 + BP_LOOP: DO + RE_D = MIN(3800._EB,MAX(40._EB,PROBE_RHO*VEL*PY%PROBE_DIAMETER/MU_G)) + FAC = 1.533_EB-0.001366_EB*RE_D+0.000001688_EB*RE_D**2-0.0000000009706_EB*RE_D**3+& + 0.0000000000002555_EB*RE_D**4-2.484E-17_EB*RE_D**5 + GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*1._EB/FAC*SQRT(2._EB*PROBE_DELTA_P/PROBE_RHO) + IF (JJJ > 9 .OR. ABS(VEL-GAS_PHASE_OUTPUT_RES)/(GAS_PHASE_OUTPUT_RES+TWO_EPSILON_EB) < 0.001_EB) EXIT BP_LOOP + VEL = 0.2_EB*VEL+0.8_EB*GAS_PHASE_OUTPUT_RES + JJJ = JJJ + 1 + ENDDO BP_LOOP ENDIF CASE(130) ! EXTINCTION