turbo-design 1.1.0__tar.gz → 1.1.2__tar.gz

{turbo_design-1.1.0 → turbo_design-1.1.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: turbo-design
-Version: 1.1.0
+Version: 1.1.2
 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.0 → turbo_design-1.1.2}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "turbo-design"
-version = "1.1.0"
+version = "1.1.2"
 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.0 → turbo_design-1.1.2}/turbodesign/bladerow.py

@@ -97,6 +97,8 @@ class BladeRow:
     phi:npt.NDArray = field(default_factory=lambda: np.array([0]))                      # Inclination angle x,r plane. AY td2.f
     rm: npt.NDArray = field(default_factory=lambda: np.array([0]))                      # Curvature
     incli_curve_radii: npt.NDArray = field(default_factory=lambda: np.array([0]))       # radius at which curvature was evaluated
+    mprime:npt.NDArray = field(default_factory=lambda: np.array([0]))                   # Mprime distance
+    axial_chord:float = 0
     
     Yp: float = 0                   # Pressure loss
     power:float = 0                 # Watts 
@@ -104,7 +106,8 @@ class BladeRow:
     P0_P:float = 0                  # Total to Static Pressure Ratio 
     Power_Type:PowerType
     euler_power:float = 0
-
+    Reynolds:float = 0
+    
     # Used for loss calculations
     _blade_to_blade_gap:float = 0.025 # Gap between blade in terms of percent chord.
     
@@ -118,7 +121,7 @@ class BladeRow:
     _tip_clearance:float = 0 # Clearance as a percentage of span or blade height
 
     _inlet_to_outlet_pratio = [0.06,0.95]
-        
+    
     @property
     def inlet_to_outlet_pratio(self) -> Tuple[float,float]:
         """This is what is varied by the optimization. 
@@ -474,6 +477,7 @@ class BladeRow:
             "P":self.P.tolist(),
             "T":self.T.tolist(),
             "rho":self.rho.tolist(),
+            "mu":self.mu,
             "Yp":self.Yp,
             "Power":self.power,
             "P0_P": self.P0_P,
@@ -482,7 +486,11 @@ class BladeRow:
             "euler_power":self.euler_power,
             "axial_chord":self.axial_chord,
             "aspect_ratio":self.aspect_ratio,
-            "area": self.area
+            "num_blades":self.num_blades,
+            "area": self.area,
+            "mprime":self.mprime[-1],
+            "Reynolds":self.Reynolds,
+            "axial_chord":self.axial_chord
         }
 
         return data

{turbo_design-1.1.0 → turbo_design-1.1.2}/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.0 → turbo_design-1.1.2}/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.0 → turbo_design-1.1.2}/turbodesign/td_math.py

@@ -5,6 +5,7 @@ import numpy.typing as npt
 from .bladerow import BladeRow, compute_gas_constants
 from .enums import RowType, LossType
 from scipy.integrate import trapezoid
+from .passage import Passage
 
 def T0_coolant_weighted_average(row:BladeRow) -> float:
     """Calculate the new weighted Total Temperature array considering coolant
@@ -78,6 +79,36 @@ def compute_massflow(row:BladeRow) -> None:
     row.calculated_massflow = massflow[-1]
     row.area = total_area
 
+def compute_reynolds(rows:List[BladeRow],passage:Passage):
+    """Calculates the Reynolds Number 
+
+    Args:
+        rows (List[BladeRow]): Blade row to calculate the Reynolds number
+        passage (Passage): Passage 
+    """
+    
+    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))
+        mp = [2/(xr[i,1]+xr[i-1,1])*np.sqrt(dr[i-1]**2 + dx[i-1]**2) for i in range(1,len(xr[:,1]))]
+        mp = np.hstack([[0],np.cumsum(mp)])
+        
+        if row.row_type == RowType.Rotor:
+            V = row.W.mean()
+        else:
+            V = row.V.mean()
+        rho = row.rho.mean()
+        mu = row.mu
+        row.Reynolds = c*V*rho/mu
+        row.mprime = mp
+        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)
+
+
+    
 def compute_power(row:BladeRow,upstream:BladeRow) -> None:
     """Calculates the power
 
@@ -169,56 +200,6 @@ def compute_quantities(row:BladeRow,upstream:BladeRow):
         row.T0R = row.T + row.W**2 / (2*row.Cp)
         row.P0R = row.P*(row.T0R/row.T)**((row.gamma)/(row.gamma-1))
    
-def compute_quantities_power(row:BladeRow,upstream:BladeRow):
-    """Calculation of all quantites after radial equilibrium has been solved assuming we know the power at the exit
-        
-    Note:
-        Radial Equilibrium gives P0, T0, Vm. This code assumes the loss either enthalpy or pressure loss has already been calculated 
-
-        compute_velocity has been called so we know W, Wt, V, Vt, U, M, M_rel
-
-        Static Pressure and Temperature should come from Total Temperature and Pressure + Velocity 
-
-    Args:
-        row (BladeRow): current blade row. All quantities are at exit
-        upstream (BladeRow): upstream blade row. All quantities are at exit
-
-    """
-    if row.row_type == RowType.Rotor:
-        Cp_avg = (row.Cp+upstream.Cp)/2
-        row.T0R = upstream.T0R - T0_coolant_weighted_average(row) - (upstream.U**2-row.U**2)/(2*Cp_avg)
-
-        # Factor in T0R_drop. Convert T0R drop to absolute terms
-        T_drop = (upstream.T0R - row.T0R) - row.W**2/(2*row.Cp) # row.T0R contains the drop
-        T0_drop = T_drop*(1+(row.gamma-1)/2*row.M**2)
-
-        # Adjust Total Temperature to match power
-        T0 = upstream.T0 - row.power/row.eta_total/(row.total_massflow*row.Cp) + T0_drop
-
-        if row.loss_function == LossType.Pressure: 
-            row.P0R = upstream.P0R - row.Yp*(upstream.P0R-row.P)
-            row.T0 = T0
-            row.T = row.T0/(1+(row.gamma-1)/2*row.M**2) 
-            row.P = row.P0R*(row.T/row.T0R)**(row.gamma/(row.gamma-1))
-            
-        elif row.loss_function == LossType.Enthalpy:                
-            row.P0 = row.P*(T0/row.T)**(row.gamma/(row.gamma-1))
-            row.T = T0 - row.W**2/(2*row.Cp) + T_drop
-            row.T0 = T0
-            row.P = row.P0*(row.T0/row.T)**(row.gamma/(row.gamma-1))
-            row.P0R = row.P * (row.T0R/row.T)**((row.gamma)/(row.gamma-1))
-
-    elif row.row_type == RowType.Stator:        
-        row.T0 = upstream.T0 - T0_coolant_weighted_average(row)
-        if row.loss_function == LossType.Pressure:
-            row.P0 = upstream.P0 - row.Yp*(upstream.P0-row.P)
-        else:
-            row.P0 = upstream.P0
-        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))
-
 def stator_calc(row:BladeRow,upstream:BladeRow,downstream:BladeRow=None,calculate_vm:bool=True):
     """Given P0, T0, P, alpha2 of stator calculate all other quantities
 

{turbo_design-1.1.0 → turbo_design-1.1.2}/turbodesign/turbinespool.py

@@ -8,7 +8,7 @@ from .passage import Passage
 from scipy.interpolate import interp1d
 import numpy as np
 import numpy.typing as npt
-from .td_math import inlet_calc,rotor_calc, stator_calc, compute_massflow, compute_power, compute_gas_constants
+from .td_math import inlet_calc,rotor_calc, stator_calc, compute_massflow, compute_power, compute_gas_constants, compute_reynolds
 from .solve_radeq import adjust_streamlines, radeq
 from scipy.optimize import minimize_scalar, minimize, fmin_slsqp
 from .inlet import Inlet
@@ -136,7 +136,7 @@ class TurbineSpool(Spool):
             self.__match_massflow()
         elif self.massflow_constraint == MassflowConstraint.BalanceMassFlow:
             self.__balance_massflow()
-    
+
     
     def __match_massflow(self):
         """ Matches the massflow between streamtubes by changing exit angles. Doesn't use radial equilibrium.
@@ -171,10 +171,11 @@ class TurbineSpool(Spool):
                             row.alpha2[0] = 1/(len(row.alpha2)-1)*row.alpha2[1:].sum()
                     upstream = compute_gas_constants(upstream)
                     row = compute_gas_constants(row)
-                    
+                
                     
             # Step 3: Adjust streamlines to evenly divide massflow
             adjust_streamlines(self.blade_rows,self.passage)
+        compute_reynolds(self.blade_rows,self.passage)
             
     def __balance_massflow(self):
         """ Balances the massflow between rows. Use radial equilibrium.
@@ -284,7 +285,10 @@ class TurbineSpool(Spool):
             self.blade_rows[-1].P = self.blade_rows[-1].get_static_pressure(self.blade_rows[-1].percent_hub_shroud)
         err = calculate_error(self.blade_rows[:-1])
         print(f"Massflow convergenced error:{err}")
-            
+        
+        # calculate Reynolds number
+        compute_reynolds(self.blade_rows,self.passage)
+        
         # finetune = True
         # if finetune:        
         #     print('Finetune static pressure between stages')
@@ -453,7 +457,7 @@ def calculate_massflows(blade_rows:List[BladeRow],calculate_vm:bool=False):
                     row = compute_gas_constants(row)
                     compute_massflow(row)
                     compute_power(row,upstream)
-
+    
 def massflow_loss_function(exit_angle:float,index:int,row:BladeRow,upstream:BladeRow,downstream:BladeRow=None):
     """Finds the blade exit angles that balance the massflow throughout the stage