
1424 H.-S. Dou et al. / International Journal of Thermal Sciences 47 (2008) 1422–1435
higher Re number and low curvature, the flow may transit to
turbulence earlier and yet does not violate Taylor’s criterion.
Recently, Dou [21,22] proposed a new energy gradient the-
ory to analyze flow instability and turbulent transition prob-
lems. In this theory, the critical condition for flow instability
depends both on the base flow and the disturbance which agrees
with the experimental observations. For a given disturbance, the
critical condition for flow instability and turbulent transition is
determined by the ratio (K ) of the gradient of total mechanical
energy in the transverse direction to the loss of total mechanical
energy in the streamwise direction. For a given flow geome-
try and fluid properties, when the maximum of K in the flow
field is larger than a critical value, it is expected that instability
would occur for some initial disturbances provided that the dis-
turbance energy is sufficiently large. For plane Poiseuille flow
(channel flow), Hagen–Poiseuille flow (pipe flow), and plane
Couette flow (simple shear flow), the findings based on the
theory are consistent with the experimental observations; for
the experimental determined critical condition, K
c
= 370–389
for all the above mentioned three types of flows below which
there is no occurrence of turbulence. In these comparisons, the
distribution of K was calculated for each flow and the value
of K
c
was obtained using the experimental data at critical con-
dition [21–23]. The theory also suggests the mechanism of
instability associated with an inflectional velocity profile for
viscous flows. The theory has been extended to curved flows
with similar derivations to parallel flows and three important
theorems have been obtained [24]. This theory has also been
employed to study the viscoelastic flows where the effect of
elastic force is dominating [25]. It should be mentioned that the
energy gradient theory is a semi-empirical theory since the crit-
ical value of K is observed and determined experimentally and
cannot be directly calculated from the theory so far. In this the-
ory, only the critical condition for the instability is sought after
and the detailed process of instability is not provided.
In this study, we apply the energy gradient theory to analyze
the Taylor–Couette flow between concentric rotating cylinders,
and aim to demonstrate that the mechanism of instability in
Taylor–Couette flow can be explained via the energy gradi-
ent concept. Through comparison with experiments, we show
that the energy gradient function K as a stability criterion is
sufficient to describe and characterize the flow instability in
Taylor–Couette flow. We also show that plane Couette flow can
be considered as just the limiting case of Taylor–Couette flow
when the curvature of the walls tends to zero. For flow between
concentric rotating cylinders, the flow instability may be in-
duced by rotation of the inner cylinder or the outer cylinder.
If it is induced by the former, a Taylor vortex cell pattern will
be formed when the critical condition is violated as in the exper-
iments; if it is induced by the latter, Taylor vortex cell pattern
will not occur and the flow may directly transit to turbulence
when the critical condition due to inertia force is reached as in
plane Couette flow [1–3,6]. In this study, only the critical con-
dition for the former situation is considered/treated.
2. Energy gradient theory revisited
Dou [21] proposed a mechanism with the aim to clarify the
phenomenon of transition from laminar flow to turbulence for
wall-bounded shear flows. In this mechanism, the whole flow
field is treated as an energy field. It is proposed that the gradi-
ent of total mechanical energy in the transverse direction of the
main flow and the total mechanical energy loss from viscous
friction in the streamwise direction dominate the instability
phenomena and hence the flow transition for a given distur-
bance. It is suggested that the energy gradient in the transverse
direction has the potential to amplify a velocity disturbance,
while the viscous friction loss in the streamwise direction can
resist and absorb this disturbance. The flow instability or the
transition to turbulence depends on the relative magnitude of
these two roles of energy gradient amplification and viscous
friction damping of the initial disturbance. In [22], more de-
tailed derivation has been given to exactly describe this mech-
anism, and this theory is termed as “energy gradient theory”.
Here, we give a short discussion for a better understanding of
the work presented in this study.
The equation of total mechanical energy for incompress-
ible flow by neglecting the gravitational energy can be written
as [21],
ρ
∂u
∂t
+∇
p +
1
2
ρu
2
= μ∇
2
u + ρ(u ×∇×u) (1)
For pressure driven flows, the derivatives of the total mechani-
cal energy in the transverse direction and the streamwise direc-
tion can be expressed, respectively, as [21–24],
∂E
∂n
=
∂(p + (1/2)ρu
2
)
∂n
= ρ(u × ω) ·
dn
|dn|
+ (μ∇
2
u) ·
dn
|dn|
= ρuω + (μ∇
2
u)
n
(2)
∂E
∂s
=
∂(p + (1/2)ρu
2
)
∂s
= ρ(u × ω) ·
ds
|ds|
+ (μ∇
2
u) ·
ds
|ds|
= (μ∇
2
u)
s
(3)
where ω =∇×u is the vorticity. Since there is no work input in
the pressure driven flows, the magnitude of the total mechanical
energy loss of unit volumetric fluid along the streamwise direc-
tion equals to the derivatives of the total mechanical energy in
the streamwise direction, that is
∂H
∂s
=−
∂E
∂s
(4)
For shear driven flows, the derivatives of the total mechan-
ical energy in the transverse direction is the same as Eq. (2).
The energy loss of unit volumetric fluid along the streamwise
direction equals to the derivatives of the total mechanical en-
ergy in the streamwise direction plus the work done to the fluid
by external,
∂H
∂s
=−
∂E
∂s
+
∂W
∂s
(5)
where W is the work done to the unit volume fluid by external.
For a given base of parallel flow, the fluid particles may move
in an oscillatory pattern in the streamwise direction if they are