Compressor Stalls and Mobius Transformations 

Linear fractional transformations, also known as Mobius transformations, arise in many diverse contexts, both theoretical and practical. For example, in the field of gas turbine design it's been found that stall cells migrate around the face of a compressor following the pattern of an iterated Mobius transformation of a point in the complex plane. 

In 1984 F. K. Moore published a paper entitled "A Theory of Rotating Stall of Multistage Axial Compressors" in the Journal of Engineering for Gas Turbines and Power (vol 106, April 1984). In Part II of this paper Moore develops an equation for the axial and circumferential velocity disturbances in the airflow as a function of angular position q. Letting g(q) and h(q) denote these disturbances, respectively, Moore's equation is 

_{} 

where 
l = collection of terms analogous to oscillator mass 
m = parameter defining lag tendency outside compressor 
f = stall propagation speed coefficient 
K = pressure rise parameter at compressor inlet vanes 
y_{c}“ = second derivative of compressor characteristic (in 
absence of rotating stall) with respect to the average 
flow coefficient, V/U (where V = average axial flow 
speed and U = wheel speed at mean wheel diameter) 
d = constant performance increase (perhaps negative) due to 
rotating stall, equal to H  y, where H is the upstream 
total to downstream static pressurerise coefficient 

Moore comments that the disturbances g and h are periodic and have vanishing averages over a cycle: in the case of g because it is defined as a departure from the average axial velocity, and in the case of h because net circulation in the entrance flow is assumed to remain zero. He also notes that we need h and g such that h + ig is an analytic function of e^{i}^{q}. With this restriction, he states that he believes the following solution of equation (1) is unique 

_{} 

where n is the "wave number", i.e., the number of stall cells, and the constants A and h (as well as d) are determined by satisfying equation (1). Obviously if h = 0 the relation degenerates to a pure exponential so we will be concerned only with the cases when h is nonzero. Separating the real and imaginary parts of the above expression we have the individual velocity disturbance components: 

_{} 

Although not mentioned in Moore’s paper, it’s interesting that these equations can be mapped to the real and imaginary parts of a sequence of complex numbers z_{k} for k = 0, 1, 2, … generated by iterating a particular linear fractional (Mobius) transformation. To generate the velocity disturbance components for a given parameter h, set the initial complex number z_{0} equal to h/(1+h) and select an arbitrary circumferential step size Dq. Then we can iterate the Mobius transformation 

_{} 

and the real and imaginary parts of z_{k} = x_{k} + iy_{k} are proportional to the axial and circumferential velocity disturbance components at the angular position q = k Dq. Specifically, we have 

_{} 

To illustrate, consider the case h = 0.3. The figure below shows the velocity disturbance profiles as given by Moore's formulas (2), and then superimposed on those curves are the discrete values generated by iterations of the Mobius transformation (3) with a value of Dq corresponding to 12 degrees. 


Of course, the general Mobius transformation is of the form z → (az+b)/(cz+d), but any such transformation is "similar" (i.e., con conjugate) to a "bipolar" transformation 

_{} 

for which the kth iterate is simply 

_{} 

By the way, the squared trace of transformation (3) is given by 

_{} 

It's interesting to note that the condition for periodicity for iteration (4) is that w be a root of 1. This Mobius transformation has fixed points at 0 and 1, and the simple linear function that transforms the general Mobius transformation to this "bipolar" LFT maps the real axis to the "Riemann line" 1/2 + yi. (For a more detailed discussion of this topic, see Linear Fractional Transformations.) 

Although the general form of (4) allows the complex constant w to have a magnitude other than 1, Moore's formula assumes w = 1. Considering that the actual physical mechanism by which the stall cell propagates from one blade to another around the face of the compressor would seem to be more naturally modeled by the discrete transformations of the form (4) rather than the continuous relations of (2), we might consider the physical consequences of allowing the magnitude of w to vary slightly from 1. The figure below shows the result of setting w to 1.010. 


This shows that the velocity disturbance tends to dampen out over several cycles. On the other hand, if we reduce the magnitude of w to 0.996 we produce the results shown in the figure below. 


Thus, for a very slight reduction from unity in the magnitude of m we find that the velocity disturbance cycle becomes unstable. It would be interesting to determine whether actual instabilities in the propagation patterns exhibited by stall cells in real compressors conform to these profiles. 

Incidentally, to confirm that (2) actually is a solution of (1), we can substitute the expressions for h and g into (1) to arrive at the equation 

_{} 

where 
_{} 

In order for (5) to be satisfied for all q each of the coefficients must vanish. Solving the equation c_{2} = 0 for d gives 

_{} 

Substituting this expression for d into the equation c_{1} = 0 and solving for A, we get 

_{} 

Likewise we can substitute the expression (6) for d into the equation c_{0} = 0 and solve for A to give 

_{} 

These two expressions for A must be equal to each other, so we have the equation A_{(1)} – A_{(0)} = 0, which we can solve for l to give the relation 

_{} 

We can now substitute this expression for ln back into either the expression for A_{(0)} and A_{(1)} to give 

_{} 

Given values of m, f, y_{c}”, l, n, and K, we can compute h from equation (7), and then we can compute A from equation (8), and finally we can compute d from equation (6). By the way, Moore noted that K is actually a function of h so it isn’t really constant, but in this application the quantity 1–K is usually negligibly small, so he deleted the terms involving 1K. Setting 1K = 0 in equations (6), (7), (8), and also substituting for A in (6), we get 

_{} 

which agrees exactly with Moore’s equations (2.12), (2.13), and (2.14). Notice that our equations (6), (7), and (8) also match Moore’s equations (2.10), (2.8), and (2.9) respectively, except that there appears to be a typo in (2.10), and possibly an algebraic error in (2.8) and (2.9). Fortuitously, the errors in the latter two equations were in the terms involving 1K, so they dropped out of the simplified results. 
