1. Thermodynamic Isentropic Expansion of Real Gases
Unlike isenthalpic throttling across a Joule-Thomson valve where enthalpy is conserved (\(\Delta h = 0\)), a turboexpander converts fluid enthalpy directly into mechanical shaft rotation. For an ideal gas undergoing reversible adiabatic (isentropic) expansion:
\[ T_{2s} = T_1 \left( \frac{P_2}{P_1} \right)^{\frac{\gamma - 1}{\gamma}} \]
For real hydrocarbon gas mixtures at high pressures (\(40\text{--}90\text{ barg}\)), real-gas compressibility factor \(Z\) alters the expansion path. The isentropic discharge temperature \(T_{2s}\) and actual discharge temperature \(T_2\) are evaluated as:
\[ T_{2s} = T_1 \left( \frac{P_2}{P_1} \right)^{\frac{Z_{\text{avg}} (\gamma - 1)}{\gamma}} \qquad T_2 = T_1 - \eta_{\text{exp}} (T_1 - T_{2s}) \]
2. Specific Isentropic Head & Shaft Work Generation
The specific isentropic head \(H_{\text{isen}}\) developed across the expander stage is calculated by integrating the real-gas \(V dP\) work:
\[ H_{\text{isen}} = Z_{\text{avg}} \frac{R}{MW} T_1 \left( \frac{\gamma}{\gamma - 1} \right) \left[ 1 - \left( \frac{P_2}{P_1} \right)^{\frac{\gamma - 1}{\gamma}} \right] \quad [\text{kJ/kg}] \]
The gross shaft power delivered to the rotor is:
\[ W_{\text{exp}} = \dot{m} \cdot H_{\text{isen}} \cdot \eta_{\text{exp}} \quad [\text{kW}] \]
3. Direct-Coupled Booster Compressor Power Balance
The expander turbine and centrifugal booster compressor share a single rigid shaft supported by active magnetic bearings (AMB) or oil film bearings. Mechanical bearing and windage losses are accounted for by mechanical efficiency \(\eta_{\text{mech}} \approx 98\%\):
\[ W_{\text{booster}} = W_{\text{exp}} \cdot \eta_{\text{mech}} = \dot{m}_{\text{residue}} \cdot \Delta h_{\text{comp}} \]
The recompression pressure ratio achieved across the single-stage booster wheel is solved from the polytropic compression head equation:
\[ \frac{P_{\text{boost,out}}}{P_{\text{boost,in}}} = \left[ 1 + \frac{\Delta h_{\text{comp}} \cdot (\gamma - 1)}{Z_{\text{avg}} R_{\text{gas}} T_{\text{suction}} \gamma} \right]^{\frac{\gamma}{\gamma - 1}} \]
4. Balje Dimensionless Specific Speed (\(N_s\))
Radial inflow expander aerodynamic peak efficiency (\(84\text{--}88\%\)) requires operating within the optimal Balje specific speed window (\(N_s = 45\text{ to }120\)):
\[ N_s = \frac{N \sqrt{Q_2}}{H_{\text{isen}}^{0.75}} \]
where \(N\) is rotational speed (\(\text{RPM}\)), \(Q_2\) is discharge volumetric flow (\(\text{m}^3/\text{s}\)), and \(H_{\text{isen}}\) is head in meters of fluid column.