turbo-design 1.3.0__py3-none-any.whl → 1.3.1__py3-none-any.whl

{turbo_design-1.3.0.dist-info → turbo_design-1.3.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
-Metadata-Version: 2.3
+Metadata-Version: 2.1
 Name: turbo-design
-Version: 1.3.0
+Version: 1.3.1
 Summary: TurboDesign is a library used to design turbines and compressors using radial equilibrium.
 Author: Paht Juangphanich
 Author-email: paht.juangphanich@nasa.gov
@@ -10,7 +10,6 @@ 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.3.0.dist-info → turbo_design-1.3.1.dist-info}/RECORD

@@ -21,13 +21,13 @@ turbodesign/loss/turbine/traupel.py,sha256=aZxFE9JDcCfi1gy-fqS1SxjjlhuDG-H-g3LPu
 turbodesign/lossinterp.py,sha256=B2KEobp-nD9jwD6UINgBmTlH9kKyWg3UNvXxqfUsr-k,6198
 turbodesign/outlet.py,sha256=-0XpeEmTQpFjUQNYQxOp5WuHkF7cXfx59hxy7BmaxgY,2119
 turbodesign/passage.py,sha256=Ur6zZCzI9Y7HW0o5LrDq9yeatRain1mBLZBcTOzgjEU,11996
-turbodesign/radeq.py,sha256=3XuSf_KtOqddJclcxh9v1V2YWMvbWfd0YxcHGkic7hg,11361
+turbodesign/radeq.py,sha256=hmHNmP5Fj0qgwq0o_j_bSyC4ktL3XFuC3QaebH5Ju_M,11361
 turbodesign/rotor.py,sha256=tHl9o5H4aQ6Etd4gqa8Ime1UK7k0de4GLt5Yb1sJdGs,1376
 turbodesign/solve_radeq.py,sha256=mGVBHmc20pAXEDKGzyiU3BVUsmY9gJl7zhbtHZqxARg,1956
 turbodesign/spool.py,sha256=XKvYABvfigLulGLHH23X7k6aL6wvfuqgfLkKCMK9X8c,14338
 turbodesign/stage.py,sha256=UP45sDKDLsAkO_WfDWJ6kqXU7cYKh_4QO01QZnSN1oQ,166
 turbodesign/td_math.py,sha256=NMKgHkc36LCrr4bZxjxx7TcA4ztaRps1cIsm1tucMrk,16069
-turbodesign/turbinespool.py,sha256=I5e1l8SsvwsQZZkOv7ric83-uNvT9MExwCimJEAhaI8,25194
-turbo_design-1.3.0.dist-info/METADATA,sha256=qSAB19ivwujv2mpnh8XwkwdSuVNIbiYNMYozn6ccPW0,734
-turbo_design-1.3.0.dist-info/WHEEL,sha256=IYZQI976HJqqOpQU6PHkJ8fb3tMNBFjg-Cn-pwAbaFM,88
-turbo_design-1.3.0.dist-info/RECORD,,
+turbodesign/turbinespool.py,sha256=VmoyLH-TqaX9uXF5xWhwNriLPKZ5trjLiRg8NdV9MiE,25268
+turbo_design-1.3.1.dist-info/METADATA,sha256=3b9pATZ2DBpozyTFfEHxS5YxjY-pCCfwgYEQ5NGnjkA,683
+turbo_design-1.3.1.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
+turbo_design-1.3.1.dist-info/RECORD,,

{turbo_design-1.3.0.dist-info → turbo_design-1.3.1.dist-info}/WHEEL

@@ -1,4 +1,4 @@
 Wheel-Version: 1.0
-Generator: poetry-core 2.0.1
+Generator: poetry-core 1.9.0
 Root-Is-Purelib: true
 Tag: py3-none-any

turbodesign/radeq.py

@@ -102,13 +102,13 @@ def radeq(row:BladeRow,upstream:BladeRow,downstream:BladeRow=None) -> BladeRow:
         if (C>1) & ((gamma/(gamma-1))<2):
             raise Exception("Invalid value of C {C}, change reduce alpha or Vm")
         B = (1-C)**(gamma/(gamma-1))      
-        A = P0 * gamma/(gamma-1) * (1-C)**(1/(gamma-1))
+        A = P0 * gamma/(gamma-1) * (1-C)**(1/(gamma-1)) * (1+np.tan(alpha)**2)/(2*Cp)
         
         epsilon = 1e-10  # or another small threshold
         if abs(rm) > epsilon:
-            dVm_dr = Cp*T0/(Vm*A*(1+np.tan(alpha)**2)) * (rho*(Vt**2/r + Vm**2/rm * np.cos(phi) - Vr*dVm_dm) - B*dP0_dr) + Vm/(2*T0) *dT0_dr  # Eqn 6
+            dVm_dr = Vm/(2*T0**2) * dT0_dr - rho/(2*A*Vm) * (Vt**2/r + Vm/rm*np.cos(phi)-Vr*dVm_dm) + B/(2*A*Vm) * dP0_dr  # Eqn 21
         else:
-            dVm_dr = Cp*T0/(Vm*A*(1+np.tan(alpha)**2)) * (rho*(Vt**2/r - Vr*dVm_dm) - B*dP0_dr) + Vm/(2*T0) *dT0_dr  # Eqn 6
+            dVm_dr = Vm/(2*T0**2) * dT0_dr - rho/(2*A*Vm) * (Vt**2/r) + B/(2*A*Vm) * dP0_dr  # Eqn 21, modified 
         
         ydot = np.array([dP0_dr,dT0_dr,dVm_dr])
 

turbodesign/turbinespool.py

@@ -278,7 +278,7 @@ class TurbineSpool(Spool):
         
         print(f"Looping to converge massflow")
         past_err = -100; loop_iter = 0; err = 0.001
-        while np.abs((err-past_err)/err)>0.05:
+        while (np.abs((err-past_err)/err)>0.05) and loop_iter<10:
             if len(outlet_P) == 1:
                 # x = balance_massflows(0.658,self.blade_rows[:-1],self.blade_rows[0].P0,self.blade_rows[-1].P)
                 res = minimize_scalar(fun=balance_massflows,args=(self.blade_rows,self.blade_rows[0].P0,self.blade_rows[-1].P),bounds=outlet_P[0],tol=0.001,options={'disp': True},method='bounded')
@@ -302,7 +302,9 @@ class TurbineSpool(Spool):
             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)
             
+            past_err = err
             err = self.__massflow_std__(self.blade_rows)
+            loop_iter += 1 
             print(f"Loop {loop_iter} massflow convergenced error:{err}")
         
         # calculate Reynolds number