fixed bug in metal temperature calc with HX

MIT-LAE · Jun 11, 2024 · 98b029f · 98b029f
1 parent dc5e481
commit 98b029f
Showing 1 changed file with 8 additions and 8 deletions.
diff --git a/src/engine/tfoper.jl b/src/engine/tfoper.jl
@@ -907,9 +907,16 @@ function tfoper!(gee, M0, T0, p0, a0, Tref, pref,
                   Trr_Tb = Trr_gmi4 * (gmi4_Rt4 * Rt4_Tb + gmi4_cpt4 * cpt4_Tb)
 
                   if (icool == 1)
+                        # Heat exchanger to cool turbine cooling air
+                        ht_tc = ht3 + Δh_TurbC #Specific enthalpy of turbine cooling air
+                        Tt_tc, Tttc_httc, _ = gas_tsetd(alpha, nair, ht_tc, Tt3) #Temperature of turbine cooling air
+
+                        httc_ht3 = 1.0
+                        Tttc_Tt3 = Tttc_httc * httc_ht3 / Tt3_ht3
+
                         #------ epsrow(.) is assumed to be passed in.. calculate Tmrow(.)
                         Tmrow_copy = Tmcalc(ncrowx, ncrow,
-                              Tt3, Tb, dTstrk, Trrat,
+                              Tt_tc, Tb, dTstrk, Trrat,
                               efilm, tfilm, StA, epsrow)
                         Tmrow[:] = Tmrow_copy[:]
 
@@ -928,13 +935,6 @@ function tfoper!(gee, M0, T0, p0, a0, Tref, pref,
                         fc_Tb = 0.0
                         fc_Mi = fc_fo * fo_Mi
 
-                        # Heat exchanger to cool turbine cooling air
-                        ht_tc = ht3 + Δh_TurbC #Specific enthalpy of turbine cooling air
-                        Tt_tc, Tttc_httc, _ = gas_tsetd(alpha, nair, ht_tc, Tt3) #Temperature of turbine cooling air
-
-                        httc_ht3 = 1.0
-                        Tttc_Tt3 = Tttc_httc * httc_ht3 / Tt3_ht3
-
                   else
                         # Heat exchanger to cool turbine cooling air
                         ht_tc = ht3 + Δh_TurbC #Specific enthalpy of turbine cooling air