Size 1 MW axial turbine

------------------------------------------------------------------------
See also: AxialTurbineDesign, DisplayLatexTable
------------------------------------------------------------------------
%--------------------------------------------------------------------------
%   Copyright (c) 2018 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------

tI      = 300;
t0      = 1589;
atmToPa = 101325;
power   = 0.5; % MW
cP      = 424; % He/Xe mixture
mDot    = power*1e6/(cP*(t0-tI))
[~, l] = AxialTurbineDesign( 'name',          '1 MW Engine',...
                             'power',         power,...
                             'r tip',         0.1,...
                             'omega',         40000*pi/30,...
                             'gamma',         1.659,...
                             'heat capacity', cP,...
                             'mass flow',     mDot,...
                             't0 inlet',    1589,...
                             'p0 inlet',    2*atmToPa);

DisplayLatexTable(l);


%--------------------------------------

% $Id: ca093fb9f18fa576cc1a8c96a69887213ce9fa00 $
mDot =

      0.91485


rH =

     0.086965      0.08117     0.063668     0.045604


sigma =

   6.7262e+07   9.4156e+07   1.6412e+08    2.186e+08


aBlade =

      0.23255

                        Name         1 MW Engine          
                    Reaction                0.50          
                       Power                0.50       MW 
                       $r_t$               0.100        m 
                       $r_h$               0.087        m 
                       $r_m$               0.093        m 
                    $\Omega$             4188.79    rad/s 
                   $\dot{m}$                0.91     kg/s 
                      $\phi$                0.40        	 
                      $\psi$                1.22        	 
                    N stages                   3          
                         $u$               391.6      m/s 
                       $t_0$              1589.0    deg-K 
                       $p_0$            202650.0       Pa 
                    $\eta_P$               0.771          
                      $\eta$               0.832          
                    $\gamma$                1.66          
                       $A_x$                0.01    m$^2$ 
                   $\zeta_N$                0.00          
                   $\zeta_R$                0.00          
                      $\eta$                0.83          
                    $\rho_b$             4800.00 kg/m$^3$ 
                    $\rho_s$             4800.00 kg/m$^3$ 
                    $\rho_h$             4800.00 kg/m$^3$ 
                   $\sigmaY$              940.00      MPA 
                    $\alpha$ [15.07 70.02 15.07]      deg 
                     $\beta$       [15.07 70.02]      deg 
                       $m_s$               11.12       kg 
                       $m_b$               14.48       kg 
                       $m_h$                6.77       kg 
                       $m_t$               32.37       kg 
                       $L_s$                0.00        m 
                  $N blades$              433.00          
                   $spacing$                0.00          
$Final Stage Pressure Ratio$                0.15