Rotating Gas Pipe (subsonic applications)

In the present section a rotating gas pipe with a varying cross section and friction is considered. Although the gas pipe Fanno is a special case of the rotating gas pipe, its governing equations constitute a singular limit to the equations presented here. Therefore, for a gas pipe without rotation and with constant cross section the equations here do not apply. The equivalent of Equation (

) now reads ([26], Table 10.2 on page 515):

\frac{dM^2}{ M^2} = \left [ \frac{ 1 + \frac{\kappa -1}{2} M^2 }{1-M^2} \right ] \left [ -2 \frac{dA}{A} - \frac{r \omega^2}{c_p T_t}\left( \frac{\kappa+1}{\kappa -1} \right) dx + \frac{\lambda dx}{D} \kappa M^2 \right ],

(99)

where r is the shortest distance from the rotational axis, \omega is the rotational speed and A is the local cross section of the pipe. Assuming that the radius R of the pipe varies linearly along its length 0 <= x <= L:

R(x)=\frac{(L-x)R_1 + x R_2}{L},

(100)

\frac{dA}{A} = \frac{2(R_2-R_1)dx}{(L-x)R_1 + x R_2}.

(101)

Taking for r, R, D and T_t the mean of their values at the end of the pipe one obtains for the second term in Equation (

) [\alpha + \beta M^2]dx where

\alpha = \frac{-8(D_2-D_1)}{L(D_1+D_2)} - ( r_1+r_2 ) \frac{\omega^2}{c_p ({T_t}_1+{T_t}_2)} \left ( \frac{\kappa+1}{\kappa-1} \right )

(102)

\beta = \frac{2 \lambda \kappa }{D_1 + D_2} .

(103)

\frac{dZ}{Z} = \left [ \frac{1+\frac{\kappa-1}{2}Z }{1-Z} \right ] (\alpha + \beta Z) dx,

(104)

\frac{a}{Z} + \frac{b}{\alpha+\beta Z} + \frac{c}{1+\frac{\kappa-1}{2}Z } =dx,

(105)

a = \frac{1}{\alpha },

(106)

b= \frac{2(\alpha+\beta) \beta }{\alpha[\alpha(\kappa-1)-2 \beta]}

(107)

c=\frac{-(1+\kappa)(1-\kappa)}{2[\alpha(1-\kappa)+2 \beta]} .

(108)

From the above equations one notices that for a non-rotating pipe with constant cross section \alpha=0 and a and b become undeterminate. Therefore, although the gas pipe Fanno is a special case, the present formulas cannot be used for this element type. Integrating Equation (

) leads to:

f:= a \ln \frac{Z2}{Z1} + \frac{b}{\beta } \ln \frac{\alpha + \beta Z_2}{\alpha + \beta Z_1} + \frac{2c}{\kappa-1} \ln \left( \frac{1 + \frac{\kappa-1}{2}Z_2 }{1 + \frac{\kappa-1}{2}Z_1} \right ) = L.

(109)

\frac{\partial f}{\partial M_1} = - \left [ \frac{a}{Z_1} + \frac{b}{(\alpha + \beta Z_1)} + \frac{c}{\left ( 1 + \frac{\kappa -1}{2} Z_1 \right) } \right] 2 M_1

(110)

\frac{\partial f}{\partial M_2} = \left [ \frac{a}{Z_2} + \frac{b}{(\alpha + \beta Z_2)} + \frac{c}{\left ( 1 + \frac{\kappa -1}{2} Z_2 \right) } \right] 2 M_2.

(111)

Focussing on the subsonic range, one has 0 \le Z_1,Z_2 \le 1. Therefore, the only term in Equation (

) which may cause problems is the second term. This is because \alpha and \beta do not necessarily have the same sign, therefore the logarithm may be undefined, i.e. the function \alpha + \beta Z may have a zero in between the ends of the pipe. This boils down to the condition (cf. Equation (

)) that in part of the element the Mach number is increasing and in part decreasing.

In general, convergence of a pipe and friction leads to increasing Mach numbers, divergence and centrifugal forces to decreasing Mach numbers. Sonic conditions should be avoided during the calculation. Especially if sonic conditions are observed at the end of a converged calculation, the result may not be correct.

Although the rotating pipe is adiabatic, i.e. no heat is transported to the envoronment, the total temperature changes due to conversion of rotational energy into heat or vice versa. Centrifugal motion leads to a total temperature increase, centripetal motion to a decrease. The change in total temperature amounts to [26]:

dT_t = \frac{r \omega^2}{c_p} dx.

(112)

T_t - {T_t}_1 = \frac{\omega^2 }{c_p} \left [ r_1 + \left ( \frac{r_2-r_1}{2} \right ) \frac{x}{L} \right ] x.

(113)

Evaluating this expression for x=L yields the total temperature increase across the pipe. In order to estimate the total pressure increase (e.g. to arrive at sensible initial conditions) one can again use the formulas in [26] (discarding the friction effect):

\frac{dp_t}{p_t} = \frac{\kappa }{\kappa-1} \frac{r \omega^2}{c_p T_t} dx.

(114)

Substituting a linear relationship for r and the result just derived for T_t leads to:

\displaystyle \frac{dp_t}{p_t}	\displaystyle = \left( \frac{\kappa }{\kappa-1}\right) \frac{\omega^2}{c_p} \fr... ...ft [ r_1 + \left( \frac{r_2-r_1}{2} \right ) \frac{x}{L} \right ] x \right \} }	(115)
	\displaystyle = \left ( \frac{ \kappa }{\kappa-1}\right) \frac{2\left [ x + \fr... ...2 + \frac{2Lr_1}{r_2-r_1}x + \frac{2Lc_p {T_t}_1}{\omega^2(r_2-r_1)} \right ] }	(116)
	\displaystyle = \left ( \frac{ \kappa }{\kappa-1}\right) d \ln \left [ x^2 + \frac{2Lr_1}{r_2-r_1}x + \frac{2Lc_p {T_t}_1}{\omega^2(r_2-r_1)} \right ].	(117)

\frac{{p_t}_2}{{p_t}_1} = \left [ 1 + \frac{L \omega^2}{c_p {T_t}_1} \left ( \frac{r_1+r_2}{2} \right ) \right ] ^ {\left (\frac{\kappa }{\kappa-1} \right)} .

(118)

It is important to notice that the rotating gas pipe is to be used in the relative (rotational) system (since the centrifugal force only exists in the rotational system). If used in the absolute system it has to be preceded by an absolute to relative element and followed by a relative to absolute element.

The rotating gas pipe is described by the following parameters (to be specified in that order on the line beneath the *FLUID SECTION, TYPE=ROTATING GAS PIPE card):