turbo-design 1.0.9__tar.gz → 1.0.11__tar.gz

{turbo_design-1.0.9 → turbo_design-1.0.11}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.3
 Name: turbo-design
-Version: 1.0.9
+Version: 1.0.11
 Summary: TurboDesign is a library used to design turbines and compressors using radial equilibrium.
 Author: Paht Juangphanich
 Author-email: paht.juangphanich@nasa.gov
@@ -10,6 +10,7 @@ Classifier: Programming Language :: Python :: 3.9
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
 Requires-Dist: cantera
 Requires-Dist: findiff
 Requires-Dist: matplotlib

{turbo_design-1.0.9 → turbo_design-1.0.11}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "turbo-design"
-version = "1.0.9"
+version = "1.0.11"
 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.0.9 → turbo_design-1.0.11}/turbodesign/passage.py

@@ -59,7 +59,7 @@ class Passage:
         """Gets the streamline at a certain percent radius 
 
         Args:
-            t_radial (float): _description_
+            t_radial (float): percent between hub and shroud
 
         Returns:
             Tuple containing:
@@ -186,7 +186,7 @@ class Passage:
         
         plt.figure(num=1,clear=True,dpi=150,figsize=(15,10))
         plt.plot(self.xhub_pts,self.rhub_pts,label='hub',linestyle='solid',linewidth=2,color='black')
-        plt.plot(self.xshroud_pts,self.rshroud_pts,label='hub',linestyle='solid',linewidth=2,color='black')
+        plt.plot(self.xshroud_pts,self.rshroud_pts,label='shroud',linestyle='solid',linewidth=2,color='black')
         for p in percent_axial:
             cut,_,_ = self.get_cutting_line(p)
             x,r = cut.get_point(np.linspace(0,1,10))

{turbo_design-1.0.9 → turbo_design-1.0.11}/turbodesign/solve_radeq.py

@@ -34,4 +34,5 @@ def adjust_streamlines(blade_rows:List[BladeRow],passage:Passage):
             phi, rm, r = passage.streamline_curvature(x_streamline,r_streamline)
             row.phi[i] = float(interp1d(t_streamline,phi)(row.axial_location))
             row.rm[i] = float(interp1d(t_streamline,rm)(row.axial_location))
-            row.r[i] = float(interp1d(t_streamline,r)(row.axial_location))
+            row.r[i] = float(interp1d(t_streamline,r)(row.axial_location))
+            row.x[i] = float(interp1d(t_streamline,x_streamline)(row.axial_location))

{turbo_design-1.0.9 → turbo_design-1.0.11}/turbodesign/spool.py

@@ -28,6 +28,7 @@ class Spool:
     _fluid: Solution = Solution("air.yaml")
 
     massflow_constraint: MassflowConstraint
+    _adjust_streamlines:bool = True
     
     # Inlet Conditions 
     def __init__(self,passage:Passage,
@@ -86,7 +87,26 @@ class Spool:
             new_gas (ct.Solution.Solution): New gas mixture
         """
         self._fluid = newFluid
+        for i in range(len(self.blade_rows)):
+            self.blade_rows[i].fluid = self.fluid
+
+    @property
+    def adjust_streamlines(self) -> bool:
+        """Variable that determines whether or not to try and adjust the streamlines to balance massflow
+        This can be a source of failure. Setting it to false may help with debugging. 
 
+        Returns:
+            bool: True = streamlines are adjusted after calculation; False = No changes in streamlines
+        """
+        return self._adjust_streamlines
+    
+    @adjust_streamlines.setter
+    def adjust_streamlines(self,val:bool):
+        """Variable that determines whether or not to try and adjust the streamlines to balance massflow
+        This can be a source of failure. Setting it to false may help with debugging. 
+        """
+        self._adjust_streamlines = val
+    
 
     def set_blade_row_rpm(self,index:int,rpm:float):
         """sets the rpm of a particular blade row"""    
@@ -229,7 +249,7 @@ class Spool:
             x_start = 0
             y_max = 0; y_min = 0
             plt.figure(num=1,clear=True)
-            for i in range(1,len(self.blade_rows)):
+            for i in range(1,len(self.blade_rows)-1): # Don't plot inlet or outlet 
                 row = self.blade_rows[i]
                 x_end = x_start + row.Vx.mean()
                 dx = x_end - x_start

{turbo_design-1.0.9 → turbo_design-1.0.11}/turbodesign/td_math.py

@@ -4,6 +4,7 @@ import math
 import numpy.typing as npt
 from .bladerow import BladeRow, compute_gas_constants
 from .enums import RowType, LossType
+from scipy.integrate import trapezoid
 
 def T0_coolant_weighted_average(row:BladeRow) -> float:
     """Calculate the new weighted Total Temperature array considering coolant
@@ -99,7 +100,7 @@ def compute_power(row:BladeRow,upstream:BladeRow) -> None:
         row.eta_static = row.power/ (row.massflow[-1]*row.Cp*(upstream.T0.mean()-row.T0_is.mean()))
         row.eta_total = row.power / (row.massflow[-1]*row.Cp * (upstream.T0.mean()-row.T0_is.mean()))
         row.stage_loading = row.Cp*(upstream.T0.mean() - row.T0.mean())/row.U.mean()**2
-        row.euler_power = (row.massflow[-1]*((row.U+upstream.U)/2).mean()*(upstream.Vt-row.Vt).mean())
+        row.euler_power = row.massflow[-1]* (upstream.U*upstream.Vt - row.U*row.Vt).mean()
     
 def compute_quantities(row:BladeRow,upstream:BladeRow):
     """Calculation of all quantites after radial equilibrium has been solved assuming we know the static pressure at the exit.
@@ -360,8 +361,10 @@ def inlet_calc(row:BladeRow):
         row (BladeRow): _description_
     """
     area = row.Vm.copy()*0
+    # Estimate the density
     row.T = row.T0
     row.P = row.P0 
+    row.rho = row.P/(row.T*row.R)
     total_area = 0 
     for _ in range(5): # Lets converge the Mach and Total and Static pressures
         for j in range(1,len(row.percent_hub_shroud)):
@@ -381,7 +384,7 @@ def inlet_calc(row:BladeRow):
         row.Vr = row.Vm*np.sin(row.phi)
         row.Vt = row.Vm*np.cos(row.phi)*np.tan(row.alpha2)
         row.V = np.sqrt(row.Vt**2+row.Vm**2)
-           
+        # Fine tune the Temperature and Pressure and density
         row.M = row.V/np.sqrt(row.gamma*row.R*row.T)
         row.T = row.T0 * 1/(1+(row.gamma-1)/2*row.M**2)
         row.P = row.P0 * (row.T/row.T0)**(row.gamma/(row.gamma-1))

{turbo_design-1.0.9 → turbo_design-1.0.11}/turbodesign/turbinespool.py

@@ -239,7 +239,7 @@ class TurbineSpool(Spool):
                 else:
                     for i in range(1,len(blade_rows)-1):
                         for j in range(self.num_streamlines):
-                            blade_rows[i].P[j] = P[i][j]*x0[(i-1)*self.num_streamlines+j]    # x0 size = num_streamlines -1 
+                            blade_rows[i].P[j] = P[j]*x0[(i-1)*self.num_streamlines+j]    # x0 size = num_streamlines -1 
                     
                 calculate_massflows(blade_rows,True)
                 print(x0)
@@ -247,7 +247,7 @@ class TurbineSpool(Spool):
             except:
                 for i in range(1,len(blade_rows)-1):
                     for j in range(self.num_streamlines):
-                        blade_rows[i].P[j] = P[i][j]
+                        blade_rows[i].P[j] = P[j]
                 calculate_massflows(blade_rows,True)
                 return 10
         # Break apart the rows to stages
@@ -271,11 +271,15 @@ class TurbineSpool(Spool):
         self.blade_rows[0].massflow = np.linspace(0,1,self.num_streamlines)*self.blade_rows[1].total_massflow_no_coolant
         inlet_calc(self.blade_rows[0]) # adjust the inlet to match massflow 
         
-        for _ in range(2):
-            adjust_streamlines(self.blade_rows[:-1],self.passage)
+        if self.adjust_streamlines:
+            for _ in range(2):
+                adjust_streamlines(self.blade_rows[:-1],self.passage)
+                self.blade_rows[-1].transfer_quantities(self.blade_rows[-2])
+                self.blade_rows[-1].P = self.blade_rows[-1].get_static_pressure(self.blade_rows[-1].percent_hub_shroud)
+                balance_massflows(x,self.blade_rows[:-1],self.blade_rows[0].P0,self.blade_rows[-1].P) 
+        else:
             self.blade_rows[-1].transfer_quantities(self.blade_rows[-2])
             self.blade_rows[-1].P = self.blade_rows[-1].get_static_pressure(self.blade_rows[-1].percent_hub_shroud)
-            balance_massflows(x,self.blade_rows[:-1],self.blade_rows[0].P0,self.blade_rows[-1].P) 
         err = calculate_error(self.blade_rows[:-1])
         print(f"Massflow convergenced error:{err}")
             
@@ -340,7 +344,13 @@ class TurbineSpool(Spool):
                 t,x,r = self.passage.get_streamline(p)
                 x_streamline[j,indx] = float(interp1d(t,x)(row.axial_location))
                 r_streamline[j,indx] = float(interp1d(t,r)(row.axial_location))
-                
+        
+        Pratio_Total_Total = np.mean(self.blade_rows[0].P0 / self.blade_rows[-2].P0)
+        Pratio_Total_Static = np.mean(self.blade_rows[0].P0 / self.blade_rows[-2].P)
+        FlowFunction = np.mean(massflow)*np.sqrt(self.blade_rows[0].T0)*self.blade_rows[0].P0/1000 # kg sqrt(K)/(sec kPa)
+        CorrectedSpeed = self.rpm * np.pi/30 / np.sqrt(self.blade_rows[0].T0.mean())   # rad/s * 1/sqrt(K)
+        EnergyFunction = (self.blade_rows[0].T0 - self.blade_rows[-2].T0) * 0.5* (self.blade_rows[0].Cp + self.blade_rows[-2].Cp) / self.blade_rows[0].T0 # J/(KgK)
+        EnergyFunction = np.mean(EnergyFunction)
         data = {            
             "blade_rows": blade_rows,
             "massflow":np.mean(massflow),
@@ -355,7 +365,12 @@ class TurbineSpool(Spool):
             "total-total_efficiency":total_total_efficiency,
             "total-static_efficiency":total_static_efficiency,
             "stage_loading":stage_loading,
-            "degree_of_reaction":degree_of_reaction
+            "degree_of_reaction":degree_of_reaction,
+            "Pratio_Total_Total":Pratio_Total_Total,
+            "Pratio_Total_Static":Pratio_Total_Static,
+            "FlowFunction":FlowFunction,
+            "CorrectedSpeed":CorrectedSpeed,
+            "EnergyFunction":EnergyFunction
         }
         # Dump all the Python objects into a single JSON file.
         class NumpyEncoder(json.JSONEncoder):