Merge pull request #74 from fastresearchgroup/Nuss

Updated Nusselt Number

LoadedTACO/src/NusseltNumber.py

@@ -0,0 +1,58 @@
+ # -*- coding: utf-8 -*-
+"""
+Created on Sun Dec. 10 10:57:35 2023
+
+@author: skyet
+"""
+import LoadedTACO.src.HT9Props as clad
+import LoadedTACO.src.HexDhCal as Geom
+import LoadedTACO.src.HegNu as Nu
+import LoadedTACO.src.TempBulkCal as TempBulk
+import TACOCAT_Read_In_File as TCinput
+import LoadedTACO.src.Geometry_Value as Geometry
+from scipy.integrate import trapz
+from scipy.integrate import quad
+from LoadedTACO.src.Fuel_Props import Fuel_props
+from LoadedTACO.src.Geometry_Value import Core_Geometry
+from LoadedTACO.src.Coolant_Value import Coolant
+
+#References: 
+# 1.)https://www.sciencedirect.com/science/article/pii/S073519332030347X
+#2.)https://www.thermal-engineering.org/what-is-dittus-boelter-equation-definition/
+Uinlet = TCinput.Uinlet #average inlet velocity in a subchannel - m/s
+Coolant_Type = TCinput.Coolant
+Geometry_Type = TCinput.Geometry
+Pr = Coolant[Coolant_Type]["Cp"]*Coolant[Coolant_Type]["nu"]*Coolant[Coolant_Type]["rho"]/Coolant[Coolant_Type]["k"] #Prandtl Number Calculation
+Re = (Uinlet*Core_Geometry[Geometry_Type]["InnerHydraulicDiameter"]/Coolant[Coolant_Type]["nu"])
+b=0.4 #This is for heating. May need to revise in the future
+f = 0.316*Re**(-0.25) #Re should be Re1. This is to test Nusselt Number Capabilities. Needs to be revised in the future
+Dh=Core_Geometry[Geometry_Type]["InnerHydraulicDiameter"]
+L=Geometry.z
+Prw=Pr #We assume that Prandlt Number at the wall is the same as the Prandlt number at the middle of the channel
+# Defining Nusselt Number functions
+
+
+def Nu_DB(Re,Pr,b):
+    Nu=0.023*(Re**0.8)*(Pr**b)
+    return Nu
+    #Nu_DB (1,5,7)
+    #print (Nu)
+
+
+def Nu_Gn(f,Re,Pr,Dh,L,Prw):
+    Nu=(((f/8)*(Re-1000)*Pr)/(1+12.7*((f/8)**(1/2))*((Pr**(2/3))-1)))*((1+((Dh/L)**(2/3)))*((Pr/Prw)**0.11))
+    return Nu
+   # Nu_Gn (9,11,13,15,17,19)
+    #print (Nu)
+# Preparing Functions for dictionary
+Nu_dittus = Nu_DB(Re,Pr,b)
+Nu_Gniel=Nu_Gn(f,Re,Pr,Dh,L,Prw)
+
+Nus_DB={"NuCorrelation":Nu_dittus,}
+Nus_Gn={"NuCorrelation":Nu_Gniel,}
+#Dictionary
+Nusselt={
+    "DittusBoelter":Nus_DB,
+    "Gnilenski":Nus_Gn,
+}
+

TACOCAT.py

@@ -13,8 +13,8 @@
 from LoadedTACO.src.Geometry_Value import Core_Geometry
 from LoadedTACO.src.Coolant_Value import Coolant
 from LoadedTACO.src.Flux_Profiles import Fluxes
-
-#Assumptions
+from LoadedTACO.src.NusseltNumber import Nusselt
+ #Assumptions
 #1. The core thermal production is assumed to set after heat deposition
 #(i.e. gamma isn't relevant)
 
@@ -37,10 +37,11 @@
 Qth = TCinput.Qth
 Tbulkin = TCinput.Tbulkin #Bulk Temperature of coolant at the Inlet - C
 Uinlet = TCinput.Uinlet #average inlet velocity in a subchannel - m/s
-
+Nusselt_Type=TCinput.Nusselt
 #Coolant Parameters
 Cp = Coolant[Coolant_Type]["Cp"]
 rho = Coolant[Coolant_Type]["rho"]
+Nu=Nusselt[Nusselt_Type]["NuCorrelation"]
 
 #Geometry Parameters
 z = Geometry.z
@@ -70,8 +71,9 @@
 
 # Coolant Calculations
 mdot = Uinlet*Coolant[Coolant_Type]["rho"]*Core_Geometry[Geometry_Type]["CoolantFlowArea"] # Mass flow rate for the fluid - kg/s
-Pe = (Uinlet*Core_Geometry[Geometry_Type]["InnerHydraulicDiameter"]/Coolant[Coolant_Type]["nu"])*Pr # Peclet Number for Fluid
-Nu = Nu.Nu(Geometry.PtoD,Pe)
+Re = (Uinlet*Core_Geometry[Geometry_Type]["InnerHydraulicDiameter"]/Coolant[Coolant_Type]["nu"])
+Pe = Re*Pr # Peclet Number for Fluid
+Nu = Nusselt[Nu_correlation]["NuCorrelation"]
 h = Nu*Coolant[Coolant_Type]["k"]/Core_Geometry[Geometry_Type]["InnerHydraulicDiameter"] #Heat Transfer Coefficient for Rod Bundles - W/m^2 - C
 
 #Core Temperature Calculations
@@ -204,4 +206,4 @@ def HotFBulkTemp(Tbulkin,FluxPro,z,NFuel,qlinHotF,Cp,Uinlet,rho,A_xs):
 plt.grid()
 k = k + 1
 
-plt.show()
+plt.show()

TACOCAT_Read_In_File.py

@@ -4,7 +4,7 @@
 Reactor_Title = "Placeholder"
 ## Print Logicals for saving plots and data files 
 # (0 - save, 1 or higher - do not save)
-print_logic = 1
+print_logic = 0
 data_logic = 1
 
 #Provide chemical composition of fuel. Fuel types include: U_Zr10, UC, UO2, PuO2, ThO2, UN, U3Si2, MOX, U
@@ -19,6 +19,9 @@
 #Provide flux profile. Flux profile includes: Flatline, Linear, Exponential, Chopped_Cosine
 Flux = "Chopped_Cosine"
 
+#Provide desired Nusselt Number Corrilation. Corrilations include:Nus_DB and Nus_Gn
+Nusselt="Nus_DB"
+
 Hc = 0.35 #Active Height of Core is 2m
 Qth = 3*10**6 #Core Thermal Production - W
 #Mdot = #Total mass flow rate for the reactor. This will then be divided for subchannels.