turbo-design 1.1.1__tar.gz → 1.1.3__tar.gz

{turbo_design-1.1.1 → turbo_design-1.1.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: turbo-design
-Version: 1.1.1
+Version: 1.1.3
 Summary: TurboDesign is a library used to design turbines and compressors using radial equilibrium.
 Author: Paht Juangphanich
 Author-email: paht.juangphanich@nasa.gov

{turbo_design-1.1.1 → turbo_design-1.1.3}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "turbo-design"
-version = "1.1.1"
+version = "1.1.3"
 description = "TurboDesign is a library used to design turbines and compressors using radial equilibrium."
 authors = ["Paht Juangphanich <paht.juangphanich@nasa.gov>"]
 packages = [

{turbo_design-1.1.1 → turbo_design-1.1.3}/turbodesign/bladerow.py

@@ -486,6 +486,7 @@ class BladeRow:
             "euler_power":self.euler_power,
             "axial_chord":self.axial_chord,
             "aspect_ratio":self.aspect_ratio,
+            "num_blades":self.num_blades,
             "area": self.area,
             "mprime":self.mprime[-1],
             "Reynolds":self.Reynolds,

{turbo_design-1.1.1 → turbo_design-1.1.3}/turbodesign/inlet.py

@@ -54,7 +54,8 @@ class Inlet(BladeRow):
         self.beta2_metal = [0]
         self.P0_fun = interp1d(self.percent_hub_shroud,P0)
         self.T0_fun = interp1d(self.percent_hub_shroud,T0)
-        
+        self.mprime = [0]
+
         
     def initialize_velocity(self,passage:Passage,num_streamlines:int):
         """Initialize velocity calculations. Assumes streamlines and inclination angles have been calculated 
@@ -67,7 +68,6 @@ class Inlet(BladeRow):
 
         cutline,_,_ = passage.get_cutting_line(self.axial_location)
         self.x,self.r = cutline.get_point(np.linspace(0,1,num_streamlines))
-        
         for _ in range(10):
             T0_T = (1+(self.gamma-1)/2 * self.M**2)
 
@@ -115,7 +115,7 @@ class Inlet(BladeRow):
                 Area += 2*np.pi*C*(S/2*dx**2+self.r[j-1]*dx)
                 
         self.calculated_massflow = self.rho.mean()*self.Vm.mean() * Area
-    
+
 
     def get_total_pressure(self,percent_hub_shroud:Union[float,npt.NDArray]):
         """Returns the static pressure at a certain percent hub_shroud

{turbo_design-1.1.1 → turbo_design-1.1.3}/turbodesign/passage.py

@@ -184,7 +184,7 @@ class Passage:
         rshroud = self.rshroud(t_shroud)
         return line2D([xhub,rhub],[xshroud,rshroud]), t_hub, t_shroud
     
-    def get_xr_slice(self,t_span:float,axial_location:float):
+    def get_xr_slice(self,t_span:float,axial_location:Tuple[float,float]):
         """Returns the xr coordinates of a streamline, a line that is parallel to both hub and shroud
             
         Args:
@@ -194,7 +194,7 @@ class Passage:
         Returns:
             np.NDArray: _description_
         """
-        t_hub = np.linspace(0,axial_location,100)
+        t_hub = np.linspace(axial_location[0],axial_location[1],100)
         
         shroud_pts_cyl = np.vstack([self.xshroud(t_hub),self.rshroud(t_hub)]).transpose()
         hub_pts_cyl = np.vstack([self.xhub(t_hub),self.rhub(t_hub)]).transpose()
@@ -228,7 +228,9 @@ class Passage:
             cut,_,_ = self.get_cutting_line(p)
             x,r = cut.get_point(np.linspace(0,1,10))
             plt.plot(x,r,label=f'{p}',linestyle='dashed')
-            
+        
+        
+        plt.ylim([-self.rshroud_pts.max()*0.1, self.rshroud_pts.max()])
         plt.legend()
         plt.axis('scaled')
         plt.show()

{turbo_design-1.1.1 → turbo_design-1.1.3}/turbodesign/td_math.py

@@ -86,8 +86,10 @@ def compute_reynolds(rows:List[BladeRow],passage:Passage):
         rows (List[BladeRow]): Blade row to calculate the Reynolds number
         passage (Passage): Passage 
     """
-    for row in rows:
-        xr = passage.get_xr_slice(0.5,row.axial_location)
+    
+    for i in range(1,len(rows)):
+        row = rows[i]
+        xr = passage.get_xr_slice(0.5,[rows[i-1].axial_location,row.axial_location])
         dx = np.diff(xr[:,0])
         dr = np.diff(xr[:,1])
         c = np.sum(np.sqrt(dx**2+dr**2))
@@ -102,7 +104,8 @@ def compute_reynolds(rows:List[BladeRow],passage:Passage):
         mu = row.mu
         row.Reynolds = c*V*rho/mu
         row.mprime = mp
-        row.axial_chord = c
+        row.axial_chord = max(c,1E-12) # Axial chord
+        # row.num_blades = int(2*np.pi*row.r.mean() / row.pitch_to_chord * row.axial_chord)
 
 
     
@@ -192,7 +195,7 @@ def compute_quantities(row:BladeRow,upstream:BladeRow):
         else:
             row.P0 = upstream.P0
         row.T0 = upstream.T0 - T0_coolant_weighted_average(row)
-        row.T = row.T0 * (1+(row.gamma-1)/2*row.M**2)
+        row.T = row.T0 / (1+(row.gamma-1)/2*row.M**2)
         row.P = row.P0 * (row.T/row.T0)**((row.gamma)/(row.gamma-1))
         row.T0R = row.T + row.W**2 / (2*row.Cp)
         row.P0R = row.P*(row.T0R/row.T)**((row.gamma)/(row.gamma-1))