turbo-design 1.3.2__py3-none-any.whl → 1.3.3__py3-none-any.whl

{turbo_design-1.3.2.dist-info → turbo_design-1.3.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: turbo-design
-Version: 1.3.2
+Version: 1.3.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.3.3.dist-info/RECORD

@@ -0,0 +1,33 @@
+turbodesign/__init__.py,sha256=N8Nu0I1vrlEHYJZZ7yhhg-FtbNbLyrgjy7hoicygUqg,325
+turbodesign/arrayfuncs.py,sha256=GHIlTyLfeNsNCQQoh5m1aXK3iyvU6RjdYQipkqpU_ps,519
+turbodesign/bladerow.py,sha256=FOs3I1dPHG0sPvfGH0PtZK9tyoTV9C8062_xGNZ4CrA,24974
+turbodesign/cantera_gas/co2.yaml,sha256=M2o_RzxV9B9rDkgkXJC-l3voKraFZguTZuKKt4F7S_c,887
+turbodesign/compressorspool.py,sha256=z8ZVczJ-EdZvIqqZArC6UdwC5MbU2SZh_MT2CGq5avY,2600
+turbodesign/coolant.py,sha256=evDtUFOYhfZKTVAsDWzxRzUX20gTvfjj0uybaWg4CsI,427
+turbodesign/enums.py,sha256=T4aFrGMnx9fgZkrExMBIY2IbCueA0utqr_afk-TUe9Q,1188
+turbodesign/inlet.py,sha256=L_IHE8zzm3uFMSow0iO-VxI71SnyQZwRhJMlJDcwZ4o,7329
+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
+turbodesign/loss/losstype.py,sha256=Vzj3UPk1_5IA3DtL56Eb5Z7LkYyMQ97DGZvgfki8ByA,708
+turbodesign/loss/turbine/TD2.py,sha256=evyTG4K6s56NGdjxrZFD3G7zsB8_2lJTBOxT3Ju2uG8,5207
+turbodesign/loss/turbine/__init__.py,sha256=d1UrbiIicO60FhDlOemo_irYW2asUaYoqlXUz-UG9zk,285
+turbodesign/loss/turbine/ainleymathieson.py,sha256=8UbugNmxI12I68HGctJe2lQpqjKZN1b40UOTzqU6hQQ,5861
+turbodesign/loss/turbine/craigcox.py,sha256=TrbNdpGdyD5kTRFKfNmbRDhCOjA0-GT5xoluLQpO1gA,9664
+turbodesign/loss/turbine/fixedefficiency.py,sha256=JJv4GUS5q9pvuf3akjBtlcHOykRYmRVFfDDrdWqH6S4,774
+turbodesign/loss/turbine/fixedpressureloss.py,sha256=jbm8dt8wtyN1RaOgpgYZsvwM0-0lv2Y73kXleyoy0zE,683
+turbodesign/loss/turbine/kackerokapuu.py,sha256=hUmZdyA1hAhN3_xbwBhka2dyk0Aq9rqq62S2ghuKSOk,5071
+turbodesign/loss/turbine/traupel.py,sha256=aZxFE9JDcCfi1gy-fqS1SxjjlhuDG-H-g3LPurty1M0,4164
+turbodesign/lossinterp.py,sha256=B2KEobp-nD9jwD6UINgBmTlH9kKyWg3UNvXxqfUsr-k,6198
+turbodesign/outlet.py,sha256=-0XpeEmTQpFjUQNYQxOp5WuHkF7cXfx59hxy7BmaxgY,2119
+turbodesign/passage.py,sha256=Ur6zZCzI9Y7HW0o5LrDq9yeatRain1mBLZBcTOzgjEU,11996
+turbodesign/radeq.py,sha256=NC4cRJG1Kz4n4zaLS5GFoIKBRSvTPuuj1dR1dgTa1mU,6475
+turbodesign/rotor.py,sha256=tHl9o5H4aQ6Etd4gqa8Ime1UK7k0de4GLt5Yb1sJdGs,1376
+turbodesign/solve_radeq.py,sha256=nLYlRtXElPgHaoUP9jwMulRmYKTJs_uQ1eCulk2V59I,1938
+turbodesign/spool.py,sha256=63Gc6U3-6zcHEJx8RDLN4kM12poo5E-M-nlbKvHm43I,14376
+turbodesign/stage.py,sha256=UP45sDKDLsAkO_WfDWJ6kqXU7cYKh_4QO01QZnSN1oQ,166
+turbodesign/td_math.py,sha256=czq-u1otamELybdBW5tP_05rzeTHo7lnrckzQkR4rG4,16183
+turbodesign/turbinespool.py,sha256=8ttt2DuExuQFt9SXraRWhaYHQDkgyh_B6AoJqdNgRlg,25046
+turbo_design-1.3.3.dist-info/METADATA,sha256=eJGhBnk9pjFtVMPNk6RiXg9DS5y1By2thmrniE1zxfo,734
+turbo_design-1.3.3.dist-info/WHEEL,sha256=IYZQI976HJqqOpQU6PHkJ8fb3tMNBFjg-Cn-pwAbaFM,88
+turbo_design-1.3.3.dist-info/RECORD,,

turbodesign/bladerow.py

@@ -14,6 +14,7 @@ from .passage import Passage
     
 
 class BladeRow:
+    id:int = 0
     stage_id:int = 0
     row_type: RowType = RowType.Stator
     loss_function:LossBaseClass

turbodesign/radeq.py

@@ -94,12 +94,12 @@ def radeq(row:BladeRow,upstream:BladeRow,downstream:BladeRow=None) -> BladeRow:
         # dP0_dr = float(interp1d(row.percent_hub_shroud, np.gradient(row.P0,row_radius))((r-row_radius[0])/(row_radius[-1]-row_radius[0]))) 
         
         C = (1 + np.tan(alpha)**2) * Vm**2/(2*Cp*T0)
-        if (C>1) & ((gamma/(gamma-1))<2):
-            raise Exception("Invalid value of C {C}, change reduce alpha or Vm")
+        if (C>1):
+            raise Exception(f"Invalid value of C {C:0.2f} which causes Vm to be nan.\nChange reduce alpha/beta for {row.row_type} {row.id}")
         B = (1-C)**(gamma/(gamma-1))
         A = -P0 * gamma/(gamma-1) * (1-C)**(1/(gamma-1)) * (1 + np.tan(alpha)**2)/(2*Cp)
         
-        eqn15_rhs = Vt**2/r - Vm**2/rm*np.sin(phi) - Vr*dVm_dm # right hand side of equation 15
+        eqn15_rhs = Vt**2/r - Vm**2/rm*np.cos(phi) - Vr*dVm_dm # right hand side of equation 15
         eqn15_rhs_simple = Vt**2/r # right hand side of equation 15 simplified for axial machines
         
         epsilon = 1e-10  # or another small threshold
@@ -107,7 +107,6 @@ def radeq(row:BladeRow,upstream:BladeRow,downstream:BladeRow=None) -> BladeRow:
             dVm_dr = T0/(2*Vm*A) * (rho*eqn15_rhs - B*dP0_dr) + Vm/(2*T0) * dT0_dr # Eqn 21
         else:
             dVm_dr = T0/(2*Vm*A) * (rho*eqn15_rhs_simple - B*dP0_dr) + Vm/(2*T0) * dT0_dr  # Eqn 21, simple 
-        
         ydot = np.array([dP0_dr,dT0_dr,dVm_dr])
 
         return ydot
@@ -133,18 +132,22 @@ def radeq(row:BladeRow,upstream:BladeRow,downstream:BladeRow=None) -> BladeRow:
     # P0_new = interp1d(hub_to_tip,res[:,0])(row_radius)
     # T0_new = interp1d(hub_to_tip,res[:,1])(row_radius)
     # Vm_new = interp1d(hub_to_tip,res[:,2])(row_radius)
-    
+    r_eval = row.r - mean_radius 
     # mean_radius_to_tip = np.linspace(0,tip_radius-mean_radius,len(row_radius)*5)
-    res1 = solve_ivp(ode_radeq_streamtube, t_span =[0, tip_radius-mean_radius], y0 = ics)
+    res1 = solve_ivp(ode_radeq_streamtube, t_span =[0, tip_radius-mean_radius], y0 = ics,
+                     t_eval=np.linspace(0,tip_radius-mean_radius,len(row_radius)*2))
     
     # mean_radius_to_hub = np.linspace(0,hub_radius-mean_radius,len(row_radius)*5)
-    res2 = solve_ivp(ode_radeq_streamtube, t_span = [hub_radius-mean_radius,0], y0 = ics)
+    res2 = solve_ivp(ode_radeq_streamtube, t_span = [0,hub_radius-mean_radius], y0 = ics,
+                     t_eval=np.linspace(0,hub_radius-mean_radius,len(row_radius)*2))
     
     mid_to_tip_vals = res1.y.transpose()
     mid_to_tip_r = res1.t + mean_radius
     mid_to_hub_vals = res2.y.transpose()
     mid_to_hub_r = res2.t + mean_radius
     mid_to_hub_vals = np.flipud(mid_to_hub_vals)
+    mid_to_hub_r = np.flipud(mid_to_hub_r)
+    
     hub_to_tip_vals = np.concatenate([mid_to_hub_vals[:-1,:],mid_to_tip_vals])
     
     r = np.concatenate([mid_to_hub_r[:-1], mid_to_tip_r])

turbodesign/spool.py

@@ -69,6 +69,7 @@ class Spool:
                 1. Conventional: Stator-Rotor-Stator-Rotor-etc. Set the RPM equally across all
                 2. Counter Rotation: We either use the RPM already persecribed for each blade row.
             '''    
+            self.blade_rows[i].id = i
             if (type(self.blade_rows[i]) != Inlet) and (type(self.blade_rows[i]) != Outlet):
                 self.blade_rows[i].rpm = rpm
                 self.blade_rows[i].axial_chord = self.blade_rows[i].location * self.passage.hub_length

turbodesign/td_math.py

@@ -132,7 +132,8 @@ def compute_power(row:BladeRow,upstream:BladeRow) -> None:
         a = np.sqrt(row.gamma*row.R*row.T_is)
         row.T0_is = row.T_is * (1+(row.gamma-1)/2*(row.V/a)**2)
         
-        row.power = row.massflow[-1] * row.Cp * (upstream.T0.mean() - row.T0.mean())
+        row.power = row.massflow[-1] * (row.Cp * (upstream.T0 - row.T0)).mean()
+        # row.power = sum(v * w for v, w in zip(row.power[1:], np.diff(row.massflow))) # Massflow weighted average 
         row.eta_static = row.power/ (row.massflow[-1]*row.Cp*(upstream.T0.mean()-row.T_is.mean()))
         row.eta_total = (upstream.T0.mean() - row.T0.mean()) / (upstream.T0.mean() - row.T0_is.mean())
         row.stage_loading = row.Cp*(upstream.T0.mean() - row.T0.mean())/row.U.mean()**2

turbo_design-1.3.2.dist-info/RECORD

@@ -1,55 +0,0 @@
-turbodesign/.DS_Store,sha256=t7V1NfD6bO4j43hG3uNmietQ86qxjKeImkJzP2ywlbo,8196
-turbodesign/__init__.py,sha256=N8Nu0I1vrlEHYJZZ7yhhg-FtbNbLyrgjy7hoicygUqg,325
-turbodesign/arrayfuncs.py,sha256=GHIlTyLfeNsNCQQoh5m1aXK3iyvU6RjdYQipkqpU_ps,519
-turbodesign/bladerow.py,sha256=v4tX27aOL7iuN0ILWYbtSbANWZ_7_LIxksu4ux1aQ2o,24959
-turbodesign/cantera_gas/co2.yaml,sha256=M2o_RzxV9B9rDkgkXJC-l3voKraFZguTZuKKt4F7S_c,887
-turbodesign/compressorspool.py,sha256=z8ZVczJ-EdZvIqqZArC6UdwC5MbU2SZh_MT2CGq5avY,2600
-turbodesign/coolant.py,sha256=evDtUFOYhfZKTVAsDWzxRzUX20gTvfjj0uybaWg4CsI,427
-turbodesign/enums.py,sha256=T4aFrGMnx9fgZkrExMBIY2IbCueA0utqr_afk-TUe9Q,1188
-turbodesign/inlet.py,sha256=L_IHE8zzm3uFMSow0iO-VxI71SnyQZwRhJMlJDcwZ4o,7329
-turbodesign/isentropic.py,sha256=YLTmNx69e_M5fTuLOmk_KvaQ6ABv0aZsityfBihJOmI,2003
-turbodesign/loss/.DS_Store,sha256=e_rC7wNyUbhwj0B0K6KYPxjmDluA2pn-nPFKtB93cfg,8196
-turbodesign/loss/__init__.py,sha256=ZJomzXa6ElduFruURukCrFwJQXMWS2aW8JSaVQ-M2r8,46
-turbodesign/loss/__pycache__/__init__.cpython-310.pyc,sha256=O-Uef93wkpiZg2LEZJlwSRHJmyWvzUVPdQhazmUDzqY,240
-turbodesign/loss/__pycache__/__init__.cpython-312.pyc,sha256=hb-aUXzBVDra2DhnX5EHnMxHhOGg5GtGTSYeDYOSV7k,242
-turbodesign/loss/__pycache__/losstype.cpython-310.pyc,sha256=Xng0JhBRPhBVq3ajlo2iXy0TU1X8eZpfAHEiJ3nuemg,1258
-turbodesign/loss/__pycache__/losstype.cpython-312.pyc,sha256=D2cC1pJyAsSeYpm4cl09qIQG_0Z-0lQpQ9Q9JdhbueI,1729
-turbodesign/loss/compressor/__init__.py,sha256=AbpHGcgLb-kRsJGnwFEktk7uzpZOCcBY74-YBdrKVGs,1
-turbodesign/loss/losstype.py,sha256=Vzj3UPk1_5IA3DtL56Eb5Z7LkYyMQ97DGZvgfki8ByA,708
-turbodesign/loss/turbine/TD2.py,sha256=evyTG4K6s56NGdjxrZFD3G7zsB8_2lJTBOxT3Ju2uG8,5207
-turbodesign/loss/turbine/__init__.py,sha256=d1UrbiIicO60FhDlOemo_irYW2asUaYoqlXUz-UG9zk,285
-turbodesign/loss/turbine/__pycache__/TD2.cpython-310.pyc,sha256=_xmRJL3-uKEX7wIfFNJwlKrTl5NEfFhqe--NsGHlbpA,4037
-turbodesign/loss/turbine/__pycache__/TD2.cpython-312.pyc,sha256=AKWBhyn5DqDolrhmJHYN1VfOz2qbtxSlBPfOo99GnCg,5657
-turbodesign/loss/turbine/__pycache__/__init__.cpython-310.pyc,sha256=08oqidbLwK0k2Bwe_mrK7TupB3h5UrEcCG-scobWCas,533
-turbodesign/loss/turbine/__pycache__/__init__.cpython-312.pyc,sha256=7_ZJqoesRMXAMEVhltiRTEdJBwCaX9AtOEhFcBQyABM,546
-turbodesign/loss/turbine/__pycache__/ainleymathieson.cpython-310.pyc,sha256=0aWYiEVGu3IfCHQ8b4-z3T3Ov7wdqDOHW1ERkMDNGy8,4887
-turbodesign/loss/turbine/__pycache__/ainleymathieson.cpython-312.pyc,sha256=FdgZSynbR0LiwunAj4UPVGkfbtZ9uoJDaeVtr25uelo,9397
-turbodesign/loss/turbine/__pycache__/craigcox.cpython-310.pyc,sha256=xsBf4UJSwijRq6r_WQRk7q9tYpS2id5XZJiF-ovMslg,6353
-turbodesign/loss/turbine/__pycache__/craigcox.cpython-312.pyc,sha256=nTzvA31LQKYQriHuYQNmHCjTMNpqphfpRL306ilbJzI,12163
-turbodesign/loss/turbine/__pycache__/fixedefficiency.cpython-310.pyc,sha256=D00hoTdPBqqm6J7m2iMMMb3fKsnZqdbwhXS6R8ezWxI,1263
-turbodesign/loss/turbine/__pycache__/fixedefficiency.cpython-312.pyc,sha256=CGw2WgkjJfKxvpUO0Qlw19w_6HnZqX47cRJIz763vgI,1498
-turbodesign/loss/turbine/__pycache__/fixedpressureloss.cpython-310.pyc,sha256=-cwVbUmXsqGnYlUccc1ZWpSl4X4ZO1-XJtQ78mps5N0,1191
-turbodesign/loss/turbine/__pycache__/fixedpressureloss.cpython-312.pyc,sha256=anr2ZSqXg5_AFQnTZYdRw2dTHK-wCpIDJeMCTGu_RMg,1370
-turbodesign/loss/turbine/__pycache__/kackerokapuu.cpython-310.pyc,sha256=Snq6kyzAbRQ_h9WSKD5sjQs_ADcUj-JHEVqVmPY9ixE,4432
-turbodesign/loss/turbine/__pycache__/kackerokapuu.cpython-312.pyc,sha256=ldnWoi8a1T2UE_iH-aw_L7M-Xgak4gYdPa02Hwdt5xc,8574
-turbodesign/loss/turbine/__pycache__/traupel.cpython-310.pyc,sha256=OUGOncOs84iMlPd5mGgtqiY95V2NU2JHyDXLF_kwKk8,3183
-turbodesign/loss/turbine/__pycache__/traupel.cpython-312.pyc,sha256=x_BnQ49DhhL_cMYzORlOLAoxd8GgGGmd-w4ec5kHDds,6445
-turbodesign/loss/turbine/ainleymathieson.py,sha256=8UbugNmxI12I68HGctJe2lQpqjKZN1b40UOTzqU6hQQ,5861
-turbodesign/loss/turbine/craigcox.py,sha256=TrbNdpGdyD5kTRFKfNmbRDhCOjA0-GT5xoluLQpO1gA,9664
-turbodesign/loss/turbine/fixedefficiency.py,sha256=JJv4GUS5q9pvuf3akjBtlcHOykRYmRVFfDDrdWqH6S4,774
-turbodesign/loss/turbine/fixedpressureloss.py,sha256=jbm8dt8wtyN1RaOgpgYZsvwM0-0lv2Y73kXleyoy0zE,683
-turbodesign/loss/turbine/kackerokapuu.py,sha256=hUmZdyA1hAhN3_xbwBhka2dyk0Aq9rqq62S2ghuKSOk,5071
-turbodesign/loss/turbine/traupel.py,sha256=aZxFE9JDcCfi1gy-fqS1SxjjlhuDG-H-g3LPurty1M0,4164
-turbodesign/lossinterp.py,sha256=B2KEobp-nD9jwD6UINgBmTlH9kKyWg3UNvXxqfUsr-k,6198
-turbodesign/outlet.py,sha256=-0XpeEmTQpFjUQNYQxOp5WuHkF7cXfx59hxy7BmaxgY,2119
-turbodesign/passage.py,sha256=Ur6zZCzI9Y7HW0o5LrDq9yeatRain1mBLZBcTOzgjEU,11996
-turbodesign/radeq.py,sha256=1_OOOBP_UglnnPuF5WKhMXTzfzXcuMDJ5A0Zsr_2Mf0,6201
-turbodesign/rotor.py,sha256=tHl9o5H4aQ6Etd4gqa8Ime1UK7k0de4GLt5Yb1sJdGs,1376
-turbodesign/solve_radeq.py,sha256=nLYlRtXElPgHaoUP9jwMulRmYKTJs_uQ1eCulk2V59I,1938
-turbodesign/spool.py,sha256=XKvYABvfigLulGLHH23X7k6aL6wvfuqgfLkKCMK9X8c,14338
-turbodesign/stage.py,sha256=UP45sDKDLsAkO_WfDWJ6kqXU7cYKh_4QO01QZnSN1oQ,166
-turbodesign/td_math.py,sha256=1hWf-B24jxiDfJHNUVOTC_0lEMZ8iu_4GNoKE4uxm6Y,16072
-turbodesign/turbinespool.py,sha256=8ttt2DuExuQFt9SXraRWhaYHQDkgyh_B6AoJqdNgRlg,25046
-turbo_design-1.3.2.dist-info/METADATA,sha256=8QYJqpja9gYWj8N5aKgfJcSj510Odrdir761wlhav0A,734
-turbo_design-1.3.2.dist-info/WHEEL,sha256=IYZQI976HJqqOpQU6PHkJ8fb3tMNBFjg-Cn-pwAbaFM,88
-turbo_design-1.3.2.dist-info/RECORD,,