Coull Hodson Design PDF
Coull Hodson Design PDF
Flow Angles
The flow coefficient sets the axial velocity: Fig. 4 Geometry estimation using a parabolic camberline
s
w_ x turning airfoils. It is first assumed that the camber-line of a blade
Vx / U / (7) can be approximated by a parabolic curve, which was shown to be
m_ w
reasonable by Horlock [12]:
Together with the flow angles, the inter-blade velocities can now ycamber A x2 B x (10)
be calculated:
The basic approach is illustrated in Fig. 4. The camber-line of the
Vi Vx = cos ai (8) blade is assumed to align with the flow angles at inlet and exit,
Wi Vx = cos bi (9) uniquely setting the constants A and B in Eq. (10). The length of
the camber-line b and the true chord C may then be calculated
where the index i represents each inter-blade gap (Fig. 2). Figure from the known axial chord Cx . The length of the blade suction
3(a) shows the variation of the exit velocity from the stator row surface S0 will largely depend on the camber-line of the blade,
(normalized by the datum value), which is largest for designs with and so the following simple relationship was proposed:
high flow coefficient and low stage loading.
Compressible flow relationships are then used to determine the S0 1:15 b (11)
Mach numbers, and hence the total and static values of pressure,
temperature and density at the inlet and exit of each blade row. To The factor of 1.15 is based on a survey of LP turbine blades
correctly calculate the total pressures throughout the stage, the from the literature. Figure 3(c) shows the ratio of this estimated
loss coefficients must be known: the calculation is therefore iter- surface length to the axial chord over the design space. In reality
ated to arrive at the correct result. the thickness distribution of the turbine blade will also make a
Together with the specified mass flow rate, the calculated den- contribution to the surface length, but Eq. (11) is sufficiently accu-
sity determines the flow area, which expands through the machine rate for the current purposes.
to maintain constant axial velocity. The flow area determines the For profile losses, the appropriate Reynolds number is based on
span: Fig. 3(b) shows the variation over the design space. The the suction surface length and exit velocity ([5]). Figure 3(d) dem-
span is largest for designs with high stage loading and low flow onstrates that this parameter is large for highly turning designs;
coefficients, which have low axial velocity. Together with the the variation in suction surface length (Fig. 3(c)) therefore domi-
specified aspect ratios (Table 1), the axial chords of the stators nates the variation in exit flow velocity (Fig. 3(b)) in setting the
and rotors are thus determined. Reynolds number.
It is necessary at this point to make some estimates of the blade
geometry, since some loss correlations rely on true chord, 3.5 Surface Velocity Distributions, Pitch and Loss
camber-line length and surface lengths. The relationships between Prediction. One of the profile loss models considered in this pa-
these geometric parameters will be very different for high and low per is based on the preliminary design correlation previously
Fig. 3 Calculated flow and geometry parameters across the design space: (a) stator exit velocity; (b) mean span; (c) ratio of
S0 =Cx ; (d) Reynolds number
m_ passage DVh
Zw (18)
hCx P01 P2
Fig. 7 Predicted efficiency for constant Zw 1:10 using the Fig. 8 Diffusion factor for constant Zw 1:10 using the models
models of Refs. [6] and [7] of Refs. [6] and [7]
separated and transitional flows [13]. Figure 10(b) also shows a high /) is far more rapid than either the Coull and Hodson or the
weak tendency for loss to increase at high flow coefficients, which simple Denton methods suggest, and therefore does not appear to
is evident as a slight curling of the contours in the bottom right of be realistic.
the chart. This effect is driven by variations in Reynolds number The profile lost efficiency predicted using the correlations of
(Fig. 3(d)) and exit dynamic pressure (from Fig. 3(a))2. Craig and Cox [9] is presented in Fig. 10(f). Although the loss
Figure 10(c) shows the profile lost efficiency predicted by the increases for high-turning designs, it also increases strongly as the
method of Ainley and Mathieson [2], who interpolated between flow coefficient rises, which does not match the simple Denton
loss charts for nozzle blades with axial inlet flow (i.e., a1 0 analysis in Fig. 10(a). (The kink in the data is due to extrapolation
for a stator), and impulse blades (a1 a2 ) for different flow errors from Craig and Coxs charts.)
angles and pitch-to-chord ratios. A correction is then made for the This comparison shows that the model of Coull and Hodson [5]
trailing edge thickness. The overall trend of Fig. 10(c) is some- gives the most realistic predictions of profile loss across the design
what similar to Fig. 10(a) and 10(b) except that high losses are space. This method also has higher fidelity than the others, since it
predicted for low turning designs (high / and low w), which is depends on the shape of the blade velocity distributions rather
caused by extrapolating the method to the high pitch-to-chord than just the pitch.
ratios of these designs.3 The predicted profile losses for designs
with high turning (low / and high w) are very high. In fact, even
if the secondary losses are zero, this model predicts lower effi- 5.2 The Influence of Flow Angles - Secondary Loss
ciency than was observed by Smith for such designs (Fig. 1). Models. Figure 11(a) shows the predicted lost efficiency due to
Figure 10(d) and 10(e) show two later modifications of the Ain- the secondary loss model of Craig and Cox [9]. This plot shows
ley and Mathieson method. Dunham and Came [3] accounted for that secondary losses are largely responsible for the observed drop
Reynolds number and Mach number variations, which makes only in efficiency at high flow coefficients. Together with the increase
small changes to the predictions (Fig. 10(d)). This model therefore in profile loss with the stage loading coefficient (Fig. 10(b)), the
suffers from the same problems as the original method. Kacker overall trend of the Smith chart is reproduced (Fig. 9). The varia-
and Okapuu [4] suggested multiplying the Dunham and Came pre- tion across the design space in Fig. 11(a) is due to a combination
dictions by 2/3 (Fig. 10(e)), which produces more reasonable lev- of competing effects. The Craig and Cox correlations are pre-
els of efficiency. However the increase in losses for high-turning sented graphically rather than as analytical expressions, but the
designs (high w and low /) and low-turning designs (low w and trends of Fig. 11(a) are driven by variation in the following pa-
rameters, in approximately decreasing order of influence:
2
(1) Velocity ratio through the blade row. Secondary losses tend
For a high speed datum turbine, this effect is slightly magnified due to the high
Mach numbers in this region of the Smith chart.
to be higher when there is a low overall acceleration
3
Note that these designs have low flow turning, so for a constant circulation they through the stage. Designs with high flow coefficient and
must have a high pitch-to-chord ratio. low stage loading therefore have high secondary losses.
(2) Exit dynamic pressure. Pressure losses scale with the The details of the Ainley and Mathieson secondary loss model
dynamic pressure, which is high for designs with high flow are given in Ref. [14]. They chose the following form for the sec-
coefficient and low stage loading (see Fig. 3(b)). ondary loss coefficient, based on the work of Carter [15]:
(3) Aspect ratio. Secondary losses increase as the aspect ratio
is decreased. A constant aspect ratio h=Cx has been speci- CL 2 cos2 a2
Ys k (28)
fied in this study, but the Craig and Cox correlation is sensi- s=C cos3 am
tive to the ratio of span to camberline h=b, which decreases
for highly cambered blades, driving a slight increase in sec- where k is a function of the effective inter-row flow areas and the
ondary losses towards the top left of the Smith chart. hub-to-tip ratio. The vector-mean flow direction am is defined as:
(4) Blade pitch. Increasing the blade pitch causes the secondary am tan1 tan a1 tan a2 (29)
flow structures to develop to a larger size (since they are
less restricted by the blade surfaces) and penetrate further The lift coefficient CL is analogous to that for an external wing
towards the midspan region. Craig and Cox use the pitch- section, and is defined using the vector-mean velocity:
to-camberline ratio, which tends to be larger for designs s
with low stage loading coefficient. CL 2 tan a1 tan a2 cos am (30)
(5) Lift coefficient. Highly loaded blades exhibit a stronger C
pressure difference between the suction and pressure surfa-
ces; this pressure difference drives the overturning of the The square of the lift coefficient in Eq. (28) drives the large
endwall flow. According to the Craig and Cox definition, increases in secondary loss at low flow coefficient and high stage
designs with a higher stage loading coefficient have higher loading (Fig. 11(b)). A similar effect was noted in point [5] above
lift coefficients, and so tend to have higher secondary for the Craig and Cox method, but here the influence of lift coeffi-
losses. cient dominates over the other effects.4
(6) Reynolds number. The Craig and Cox model assumes mod- One problem with the Ainley and Mathieson method is that the
est changes in secondary loss with Reynolds number, definition of lift coefficient CL (Eq. (30)) is not appropriate for
approximately following a turbulent trend ( Re0:2 ). moderate or highly cambered blades, as they themselves pointed
out in Ref. [14]. For example, the vector-mean velocity for a high-
The predicted secondary lost efficiency using the method of
turning impulse blade (a1 a2 ) is simply the axial velocity Vx ,
Ainley and Mathieson [3] is presented in Fig. 11(b). The observed
which is much smaller than the actual flow velocities in the cascade.
trend is very different to the Craig and Cox result, and very high
losses are predicted for designs with high stage loading coeffi-
cient. One is left wondering why the two methods give such dif- 4
It should be noted that Ainley and Mathieson account for the machine hub-to-tip
ferent trends, and which is more representative of reality. ratio rather than the aspect ratio.
The reference velocity is therefore artificially low, and thus CL is ultra-high-lift blade designs for LP turbines. Higher lift designs
artificially high for such high-turning designs. Furthermore, Eq. (28) require fewer blades and therefore offer a weight and cost saving
was formulated for relatively low turning blades and it is therefore of for the engine, but tend to have lower efficiency. This section
questionable validity. In comparison, Craig and Cox used results examines the trade-off between lift and efficiency.
from a much wider range of blade designs to formulate their model. Calculations are shown in Fig. 12 for designs with the same
Figure 11(c) shows the predictions according to Dunham and flow angles as the datum turbine (/ 0:9, w 2) but varying
Came [3], who updated the Ainley and Mathieson secondary loss lift. Figure 12(a) shows the lost efficiency due to profile loss
model to account for aspect ratio. Predictions using the subse- alone. (Noting that these designs have velocity distributions of the
quent modification by Kacker and Okapuu [9] are presented in style shown in Fig. 5, the diffusion factor for each design is also
Fig. 11(d). Both methods produce predict lower losses than the indicated on the horizontal axes.) The method of Denton [13] pre-
original method, which is arguably more realistic, but still predict dicts the lowest loss, which is unsurprising since it assumes fully
the wrong efficiency trend across the design space. turbulent boundary layers. The Coull and Hodson [6] model indi-
Predictions using the correlation of Traupel [7] are presented in cates a modest increase in profile loss with increasing lift, consist-
Fig. 11(e). This method successfully captures the loss of effi- ent with such forward-loaded designs. However, there is a slight
ciency at high flow coefficient and low stage loading, and the pre- drop off in loss for designs with low lift (DF < 0:20), which is
dictions are very similar to those of Craig and Cox. (The very unlikely to be physical. Such designs are close to the limit of the
slight kinks at high stage loading coefficient are a result of extrap- high-lift experimental data used to formulate this method ([5]). It
olating from Traupels graphs.) should also be noted that the Coull and Hodson method does not
Figure 11(f) shows the predictions of the recent secondary loss account for the base pressure loss on the turbine trailing edge,
correlation of Benner et al. [17], assuming that the displacement which can be significantly higher for lower lift designs. Figure 12
thickness of the inlet endwall boundary layer is d 0:05h. This shows that the method of Craig and Cox [9] predicts a rapid rise
method correctly captures the increase in loss at high flow coeffi- in loss above Co 0:77, probably due to extrapolating beyond
cients. Although promising, it is not clear how to select an appro- their experimental design space. The Ainley and Mathieson [3]
priate value of d at the preliminary stages of design, which has a and Dunham and Came [4] methods both predict very high profile
large influence on the predictions. losses. For this particular set of flow angles, the Kacker and Oka-
Considering the trends across the design space, the secondary puu [5] correlation predicts similar values to the Coull and Hod-
loss models of Craig and Cox and Traupel are the most reasonable son method.
of those considered. This analysis says nothing of the absolute ac- The secondary lost efficiency is presented in Fig. 12(b) for the
curacy of these correlations, which may only be determined by models considered in the previous section. Of these, the method
direct comparison of predictions and test data. It should also be of Ainley and Mathieson [3] predicts very high losses in line with
noted that both methods are inherently approximate as they take Fig. 11. The models of Dunham and Came [3], Kacker and Oka-
no account of the shape of the surface velocity distributions, puu [5] and Benner et al. [17] (assuming an inlet endwall bound-
which are known to have a significant impact on the secondary ary layer with d 0:05h) predict that secondary loss stays
losses. For example, Gier et al. [11] compared two high-lift LP almost constant with increasing lift. These predictions are con-
turbine cascades with similar flow turning and loading. The veloc- trary to experimental experience (e.g., Ref. [11]). The Craig and
ity peak on the T162 design was relatively far forward on the sur- Cox [9] and Traupel [7] models predict a rise in secondary flow
face, and this exhibited a 30% higher secondary loss coefficient losses with lift, but the increase is significantly less than that
than the aft loaded T161. The mean-line secondary loss models observed by Gier et al. [11] and Praisner et al. [20] for high-lift
considered here lack the fidelity to capture this variation, and will LP turbine blades. Despite this inadequacy, these methods give
therefore be prone to (significant) errors. There is therefore con- the most realistic trends of the models considered.
siderable scope to develop improved methods in the future.
Fig. 13 Predicted lost efficiency with varying Reynolds number for the datum turbine:
(a) profile loss; (b) secondary loss