turbo-design 1.1.1__py3-none-any.whl → 1.1.3__py3-none-any.whl

{turbo_design-1.1.1.dist-info → turbo_design-1.1.3.dist-info}/METADATA

@@ -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.dist-info → turbo_design-1.1.3.dist-info}/RECORD

@@ -1,11 +1,11 @@
 turbodesign/__init__.py,sha256=N8Nu0I1vrlEHYJZZ7yhhg-FtbNbLyrgjy7hoicygUqg,325
 turbodesign/arrayfuncs.py,sha256=GHIlTyLfeNsNCQQoh5m1aXK3iyvU6RjdYQipkqpU_ps,519
-turbodesign/bladerow.py,sha256=Ie3TDxaLzOSM38dQ0xNyjFEwHWbHXMq1FcksTVEvkqc,24540
+turbodesign/bladerow.py,sha256=I0KYtw7c5cy32vxq1D54TrIkfsdbRlDWsnoCIADtCUA,24582
 turbodesign/cantera_gas/co2.yaml,sha256=M2o_RzxV9B9rDkgkXJC-l3voKraFZguTZuKKt4F7S_c,887
 turbodesign/compressorspool.py,sha256=z8ZVczJ-EdZvIqqZArC6UdwC5MbU2SZh_MT2CGq5avY,2600
 turbodesign/coolant.py,sha256=nuj7ryUXuBhuOPt2jHjHwtEM_LVno7A8Wm14P0lJn-8,436
 turbodesign/enums.py,sha256=T4aFrGMnx9fgZkrExMBIY2IbCueA0utqr_afk-TUe9Q,1188
-turbodesign/inlet.py,sha256=q-gzmA58SU0Z-y4Mwv7_fybJJ5Jn8t_w-6obG95IvpU,6317
+turbodesign/inlet.py,sha256=QeSvFi35JeOPevjQrz1TmAhoWkScXT9BRA_5oRTnDEE,6322
 turbodesign/isentropic.py,sha256=YLTmNx69e_M5fTuLOmk_KvaQ6ABv0aZsityfBihJOmI,2003
 turbodesign/loss/__init__.py,sha256=ZJomzXa6ElduFruURukCrFwJQXMWS2aW8JSaVQ-M2r8,46
 turbodesign/loss/compressor/__init__.py,sha256=AbpHGcgLb-kRsJGnwFEktk7uzpZOCcBY74-YBdrKVGs,1
@@ -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=eZGaB5Gq9Frm80TTDJaLm1TxJV2-_fHHIeqjPuvtb8E,9550
+turbodesign/passage.py,sha256=a3DMc-QdlFEdjOEL1cockcTUFu_1iT1PSDrGo1qKN9s,9659
 turbodesign/radeq.py,sha256=RCsiQzcHvR5dozK6XG0eysAhXYewF0BY_qX0f1knURo,10091
 turbodesign/rotor.py,sha256=tHl9o5H4aQ6Etd4gqa8Ime1UK7k0de4GLt5Yb1sJdGs,1376
 turbodesign/solve_radeq.py,sha256=2VGzFSVkyN0rTPFHaCQiD3v166cdEs39x7x7RuzLgmw,1886
 turbodesign/spool.py,sha256=h4m--HR5g8o01Kkhl4Mbvx7WORQez3YpncDE7XxbNaw,14172
 turbodesign/stage.py,sha256=UP45sDKDLsAkO_WfDWJ6kqXU7cYKh_4QO01QZnSN1oQ,166
-turbodesign/td_math.py,sha256=En8g8AcldaZN9qmG33wrt2EuTqKfSQ8UO2PWod6a0IE,15970
+turbodesign/td_math.py,sha256=GjJdH9Mym1_vMSDT7lxilzSLVHBM5-Z5UE0wfsH2hr4,16153
 turbodesign/turbinespool.py,sha256=U_TMv5XNTlRWP9NJWTaaFWfNeTEGEuSVWt3HDi8QyFI,26124
-turbo_design-1.1.1.dist-info/METADATA,sha256=C1hmGnuG5wJahAn1UDbHQPQCwfUeGK4Mx3XgwGucUHg,734
-turbo_design-1.1.1.dist-info/WHEEL,sha256=IYZQI976HJqqOpQU6PHkJ8fb3tMNBFjg-Cn-pwAbaFM,88
-turbo_design-1.1.1.dist-info/RECORD,,
+turbo_design-1.1.3.dist-info/METADATA,sha256=F1sI0aDeah1bWEhX57xeQTvt5aMMwHR9bkWyV7VFa-M,734
+turbo_design-1.1.3.dist-info/WHEEL,sha256=IYZQI976HJqqOpQU6PHkJ8fb3tMNBFjg-Cn-pwAbaFM,88
+turbo_design-1.1.3.dist-info/RECORD,,

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,

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

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

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