turbo-design 1.0.9__py3-none-any.whl → 1.0.10__py3-none-any.whl

{turbo_design-1.0.9.dist-info → turbo_design-1.0.10.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: turbo-design
-Version: 1.0.9
+Version: 1.0.10
 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.0.9.dist-info → turbo_design-1.0.10.dist-info}/RECORD

@@ -20,14 +20,14 @@ turbodesign/loss/turbine/kackerokapuu.py,sha256=hUmZdyA1hAhN3_xbwBhka2dyk0Aq9rqq
 turbodesign/loss/turbine/traupel.py,sha256=aZxFE9JDcCfi1gy-fqS1SxjjlhuDG-H-g3LPurty1M0,4164
 turbodesign/lossinterp.py,sha256=B2KEobp-nD9jwD6UINgBmTlH9kKyWg3UNvXxqfUsr-k,6198
 turbodesign/outlet.py,sha256=SwTwoL6XnWlsrE5A6wLxu3qXkydDDUsdKKeBbJm4mrQ,2069
-turbodesign/passage.py,sha256=1Erl03ISwV4fTwRjgh2cCTdNpUpT1O4zrZh8xO-1zSU,7577
+turbodesign/passage.py,sha256=9cRSCZIS4LqlIcO8mOeePN42hmcc4g8tFDZ3F86WM7E,7594
 turbodesign/radeq.py,sha256=4unpFad6hvgaYwYLBNI8o9bszyvoNiizYjMqRX0Mjcc,10012
 turbodesign/rotor.py,sha256=tHl9o5H4aQ6Etd4gqa8Ime1UK7k0de4GLt5Yb1sJdGs,1376
-turbodesign/solve_radeq.py,sha256=tdVbRAUbXaq8WCOjNUkTLQzND6orRUSEC2je3M_ly5g,1800
-turbodesign/spool.py,sha256=-svy5RbNZrrHocq6Fyejz7RA_jS7QwA-BFPx87TV6yI,13048
+turbodesign/solve_radeq.py,sha256=2VGzFSVkyN0rTPFHaCQiD3v166cdEs39x7x7RuzLgmw,1886
+turbodesign/spool.py,sha256=FImXUwB4bJSAFroae-Y0fL5T7hBCiHg17Uo2zRvcFnA,13864
 turbodesign/stage.py,sha256=UP45sDKDLsAkO_WfDWJ6kqXU7cYKh_4QO01QZnSN1oQ,166
-turbodesign/td_math.py,sha256=36HkPCtRADlnUv1GHawy9GODMZU2LN9WBscSxMtO--c,17260
-turbodesign/turbinespool.py,sha256=SmG4g5jcemAat1tbir2bpk8L2FlQ9x03IKE0LOomX0Y,24781
-turbo_design-1.0.9.dist-info/METADATA,sha256=Jzk6ZlOf-BWdBE_AeGxTMBNERUwqwdB9a_gceqja3ko,683
-turbo_design-1.0.9.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
-turbo_design-1.0.9.dist-info/RECORD,,
+turbodesign/td_math.py,sha256=28cGzmTMs2evC0vMEp0oN6UmSJDHsKfXw51XJSMxmo8,17370
+turbodesign/turbinespool.py,sha256=s9cZpK9245pjUYI5hOo81SsO8UqG-TVUJlSN7SQu7yw,25877
+turbo_design-1.0.10.dist-info/METADATA,sha256=jQ5YdRxxq3NcYV-B7beGfmgTGTsGhMzE8PC8Y3W23QU,684
+turbo_design-1.0.10.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
+turbo_design-1.0.10.dist-info/RECORD,,

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:

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))

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,
@@ -87,6 +88,23 @@ class Spool:
         """
         self._fluid = newFluid
 
+    @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 +247,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

turbodesign/td_math.py

@@ -360,8 +360,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 +383,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))

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):