Skip to main content

Dynamic phase transition for the Taylor problem in the wide-gap case

Abstract

The main objective of this paper is to study the stability and the type of transition of the Taylor problem in the wide-gap case by using the averaging method, and we conclude that the stability of the Taylor problem in the wide-gap case is essentially the same with that in the case of the narrow-gap. The main technical tools are the spectral theory for linear and completely continuous fields, the dynamic bifurcation theory and the transition theory for incompressible flows, both developed by Ma and Wang (Bifurcation Theory and Applications, 2005; Stability and Bifurcation of Nonlinear Evolution Equations, 2007).

1 Introduction

In 1921, Taylor [1] observed and studied the stability of laminar flow, which is known as the Couette flow. Taylor studied the case, where the gap between two cylinders is smaller in comparison with the mean radius, and both of them rotate in the same direction. He found that when the Taylor number T is smaller than a critical value T c >0, called the critical Taylor number, the Couette flow is stable, and when the Taylor number crosses T c , the Couette flow breaks out into a cellular pattern which is radially symmetric.

Since Taylor’s work, there have been many studies on this kind of problem. Such as Couette [2] and Mallock [3] did a lot of experiments. Cloes [4] published the most comprehensive experimental results associated with the Taylor vortex and the secondary instability. Chandrasekhar [5], Walowit et al. [6], Drazin and Reid [7] studied the linear theories. Velte [8], Kirchgässner [9], Kirchgässner and Sorger [10], Yudovich [11], Ma and Wang [1214] studied the nonlinear theories. Especially, Ma and Wang established a new notion of bifurcation, called an attractor bifurcation, which was applied to the Taylor problem and obtained a series of fine results. This paper focuses on the Taylor problem in the wide-gap case. In this case, the radius of inner cylinder is small, while the radius of the outer cylinder is big. In addition to the same direction, the two cylinders could rot in the converse direction.

The main objective is to study the stability and the type of transition of the Taylor problem in the wide-gap case by using the averaging method and to compare with the Taylor problem in the narrow-gap case.

The main technical tools are the spectral theory for linear and completely continuous fields, the dynamic bifurcation theory and the transition theory for incompressible flows. These theories are directly applied to the Taylor problem in the wide-gap case.

The main conclusion is that the stability of the Taylor problem in the wide-gap case is essentially the same with that in the case of a narrow-gap. The main theorems are presented in Section 4, through which we can give pictures depicting the Couette flow stability in the wide-gap case, and compare with the Taylor problem in the narrow-gap case. In the later research, we intend to simulate the Taylor problem in the wide-gap case by using computer.

This paper is organized as follows. Section 2 introduces the governing equations for the Taylor problem. Section 3 studies the Taylor problem in the wide-gap case by using the averaging method, and establishes its mathematical frame. All the main theorems and the proofs are presented in Section 4.

2 Governing equations for the Taylor problem

The governing equations for the Taylor problem are the Navier-Stocks equations in the cylindrical coordinates (z,r,θ), which have the form

{ u z t + ( u ) u z = z ( p ρ ) + ν Δ u z , u r t + ( u ) u r u θ 2 r = r ( p ρ ) + ν ( Δ u r 2 r 2 u θ θ u r r 2 ) , u θ t + ( u ) u θ + u r u θ r = θ ( p ρ ) + ν ( Δ u θ + 2 r 2 u r θ u θ r 2 ) , ( r u r ) r + u θ θ + ( r u z ) z = 0 ,
(2.1)

where ν is the kinematic viscosity, ρ the density, u=( u z , u r , u θ ) the velocity field, p the pressure function, and

u = u z z + u r r + u θ r θ , Δ = 2 z 2 + 1 r r ( r r ) + 1 r 2 2 θ 2 .

The basic flow for (2.1) is the Couette flow, namely a steady state solution defined by

u z = u r = 0 , u θ = V ( r ) , p = ρ 1 r V 2 ( r ) d r , V ( r ) = a r + b r ,
(2.2)

where

a = Ω 1 η 2 1 μ / η 2 1 η 2 , b = Ω 1 r 1 2 ( 1 μ ) 1 η 2 , μ = Ω 2 Ω 1 , η = r 1 r 2 .

In order to investigate the stability of the Couette flow defined by (2.1), we need to consider the perturbed state

u z , u r , u θ +V(r),p+ρ 1 r V 2 (r)dr.

Assume that the perturbations are axisymmetric and independent of θ, from (2.1), we have:

{ u z t + ( u ˜ ) u z = ν Δ u z z ( p ρ ) , u r t + ( u ˜ ) u r u θ 2 r = ν ( Δ u r u r r 2 ) r ( p ρ ) + 2 V r u θ , u θ t + ( u ˜ ) u θ + u r u θ r = ν ( Δ u θ u θ r 2 ) ( V + V r ) u r , ( r u r ) r + ( r u z ) z = 0 ,
(2.3)

where

u ˜ = u z z + u r r , Δ = 2 z 2 + 2 r 2 + 1 r r .

The spatial domain for (2.3) is M=( r 1 , r 2 )×(0,L) R 2 , where L is the height of the field between the two cylinders. There are different physically sound boundary conditions.

In the radial direction, there are two kinds of boundary conditions:

Free boundary condition:

u z r =0, u r = u θ =0,at r= r 1 , r 2 .

Rigid boundary condition:

u z = u r = u θ =0,at r= r 1 , r 2 .

In the z direction, there are four kinds of boundary conditions:

Free slip boundary condition:

u z =0, u r z = U θ z =0,at z=0,L.

Dirichlet boundary condition (or rigid condition):

u z = u r = u θ =0,at z=0,L.

Free rigid boundary condition:

u z = 0 , u r z = u θ z = 0 , at  z = L , u z = u r = u θ = 0 , at  z = 0 .

Periodic condition:

u=( u z , u r , u θ ) is periodic in the z direction.

Using the transform:

{ x = l x , x = ( z , r , θ ) , t = l 2 t / ν , u = ν u / l , u = ( u z , u r , u θ ) , p = ρ ν 2 p / l 2 ,

where l is a certain length unit in (2.3).

We get the dimensionless form:

{ u z t = Δ u z p z ( u ˜ ) u z , u r t = ( Δ u r u r r 2 ) p r + 2 Ω 1 l 2 ν ( r 1 2 r 2 1 μ 1 η 2 η 2 μ 1 η 2 u θ ) + u θ 2 r ( u ˜ ) u r , u θ t = ( Δ u θ u θ r 2 ) + 2 Ω 1 l 2 ν η 2 μ 1 η 2 u r u r u θ r ( u ˜ ) u θ , ( r u r ) r + ( r u z ) z = 0 .
(2.4)

3 Averaging for the Taylor problem in the wide-gap case

In the case of a wide-gap, r 1 is small, while r 2 is big. The averaging method is that we let l= r 1 =1, and r 0 =( r 1 + r 2 )/2 replace variable r approximately. Since r x is big, r 1 and ν are small, we can ignore r n (n1), which are associated with the inner friction of fluid. In this paper, moreover, our analysis will be conducted using free boundary conditions in z and radial directions. Then from (2.4), we obtain the averaged governing equation for the case of the wide-gap:

{ u z t = Δ u z p z ( u ˜ ) u z , u r t = Δ u r p r + λ ( 1 r 0 2 κ ) u θ + u θ 2 r 0 ( u ˜ ) u r , u θ t = Δ u θ + λ κ u r u r u θ r 0 ( u ˜ ) u θ , u r r + u z z = 0 ,
(3.1)
{ u z r = 0 , u r = u θ = 0 , at  r = 1 , r 2 , u z = 0 , u r z = u θ z = 0 , at  z = 0 , L ,
(3.2)

where λ= T , T is the Taylor number, T and κ are defined by

T = 4 Ω 1 2 ( 1 μ 2 ) ν 2 ( 1 η 2 ) , κ = η 2 μ 1 μ , u ˜ = u z z + u r r , Δ = 2 z 2 + 2 r 2 ,

spatial domain for (3.1) is M=(1, r 2 )×(0,L) R 2 .

Let u ˜ =( u z , u r ),

H = { u = ( u ˜ , u θ ) L 2 ( M , R 3 ) | div u ˜ = 0 , u ˜ n | M = 0 } , H 1 = { u = ( u ˜ , u θ ) H 2 ( M , R 3 ) H | u  satisfies  ( 3.2 ) } .

By the Hodge decomposition theorem, L 2 (M, R 3 ) can be decomposed as follows.

L 2 ( M , R 3 ) = H G , H = { u = ( u ˜ , u θ ) L 2 ( M , R 3 ) | div u ˜ = 0 , u ˜ n | M = 0 } , G = { φ L 2 ( M , R 3 ) | φ H 1 ( M ) } , G H .

Here, H 1 (M) and H 2 (M, R 3 ) are the usual Sobolev spaces.

Then we can define an orthogonal projection, called Leray projection

P: L 2 ( M , R 3 ) H.

Let

L λ = A + λ B : H 1 H , G : H 1 H

be a mapping defined by

A ( u ) = P ( Δ u z , Δ u r , Δ u θ ) , B ( u ) = P ( 0 , ( 1 r 0 2 κ ) u θ , κ u r ) , G ( u ) = P ( ( u ˜ ) u z , ( u ˜ ) u r u θ 2 r 0 , ( u ˜ ) u θ + u θ u r r 0 ) ,
(3.3)

where u=( u z , u r , u θ ) H 1 , and the mapping P is the Leray projection.

Thus, (3.1) can be written in the abstract form:

d u d t = L λ u+Gu,u H 1 .
(3.4)

So far, the stability of the Couette flows of equation (2.1) equals to the stability of ( u z , u r , u θ )=0 of equation (3.4).

4 Main results and proofs

Let X 1 , X be two Banach spaces, X 1 X be a dense and compact inclusion. Consider the following nonlinear evolution equations

{ d u d t = L u + G u , u ( 0 ) = φ ,
(4.1)

where L: X 1 X is bounded linear operator, G: X 1 X is C r (r1) mapping.

Definition 4.1 [14]

Let ΩX be a bounded open set, we say that (4.1) is Lyapunov stable in Ω, if the solution of (4.1) u(t,φ)Ω, t0, for any initial point φΩ.

Definition 4.2 [14]

We say (4.1) is asymptotically Lyapunov stable in Ω, or u 0 is an asymptotically stable equilibrium point of (4.1) in Ω, if there exists a steady solution u 0 Ω, satisfying

lim t u ( t , φ ) u 0 X =0for any φΩ.

Lemma 4.3 [14]

Let L: X 1 x be a sectorial operator, G: X α X be a C r (r1) mapping for certain 0α<1, v be a steady solution of (4.1). If the spectral β(λ) of L+DG(v) satisfy Reβ(λ)<ε, for certain ε>0, then v is a locally asymptotically stable equilibrium point of (4.1), and it decays exponentially. Namely, there exist M>0, σ>0, δ>0, for the solution of (4.1) u(t,φ), u ( t , φ ) v X M e σ t (t0) hold true, as φv<δ.

Let

δ= 3 r 2 2 2 r 2 1 r 2 2 ( r 2 1 ) ( r 2 + 3 ) >0,as  r 2 > r 1 =1,

and as above

μ= Ω 2 Ω 1 ,η= r 1 r 2 .

Theorem 4.4 If μδ, or η 2 μ, then ( u z , u r , u θ )=0 is a locally asymptotically stable equilibrium point of (3.4), and it decays exponentially.

Theorem 4.5 If δ<μ< η 2 (μ0), we have the following conclusions for (3.4):

  1. (1)

    Equation (3.4) happens with a continuous transition at (u,λ)=(0, T c ), namely there is an attractor bifurcation, and it bifurcates exactly into two singular points v i (i=1,2), which attract two open subsets of U separately. U is the neighborhood of u=0.

  2. (2)

    The two singular points v i (i=1,2) can be written by:

    v 1 , 2 (λ)=± | β n 1 1 / σ | 1 2 ψ n 1 1 +o ( | β n 1 1 / σ | 1 2 ) ,

where σ= 1 4 β 2 n 1 2 a 1 4 b 1 2 , β n 1 1 is the first eigenvalue.

Remark 4.6 From the point of view of mathematics, Theorem 4.4 explains that there exists an open set Ω H 1 , which guarantees that if the initial point u 0 Ω, then the solution of (3.4) satisfies u ( t , u 0 ) H 0 (t).

Remark 4.7 From the point of view of physics, Theorem 4.4 explains that in the wide-gap case, the Couette flow of the Taylor problem is metastable. Namely, if the initial perturbation is in a certain range, the disturbed fluid will become the Couette flow in a short time. But if the initial perturbation is beyond that range above, the disturbed fluid will become another steady flow. The explanation is the same as that for Theorem 4.5.

Remark 4.8 Theorems 4.4-4.5 give the entire results of stability of the Taylor problem in the wide-gap case; see Figure 1. Σ 1 represents the unstable area, Σ 2 represents the stable area, T 0 0 and

T 0 = min ( n , l ) α 3 ν 2 ( 1 η 2 ) 2 4 a 2 η 2 ( 1 r 0 2 η 2 ) .
Figure 1
figure 1

Wide-gap Case.

Remark 4.9 [15]

In the narrow-gap case, let r 1 =3.55cm and r 2 =4.035cm, the results can be seen in Figure 2. Σ 1 represents the unstable area, Σ 2 represents the stable area.

Figure 2
figure 2

Narrow-gap Case.

Remark 4.10 By Figure 1 and Figure 2, we can conclude that the results of the stability of the Taylor problem in the wide-gap case, obtained by using the averaging method, are essentially the same as those of the Taylor problem in the narrow-gap case.

Proof of Theorem 4.4 Obviously, A and B are linear operators. According to [16], A is a homeomorphism, then it is a sectorial operator. According to the Sobolev compact embedding theorem [17] and the Leray projection, P is bounded, B is compact. So,

L λ =A+λB: H 1 H

is a sectorial operator, and also a linear completely continuous field.

Now,

Ω ( u ) u v + u 2 v d x u C 0 u L 2 v L 2 + u C 0 u L 2 v L 2 = u C 0 u W 1 , 2 v L 2 ,

we have

G u H u C 0 u W 1 , 2 .

According to [18], Theorem 7.2.2, the fraction spaces of L λ , H α =D( L λ α ), satisfy

u C 0 + u W 1 , 2 C u H α ,

where α> 3 4 . Then

G: H α H ( 3 4 < α < 1 )  is a bounded analytic mapping.

Consider the eigenvalue problem:

L λ ψ=β(λ)ψ,ψ=( u z , u r , u θ ) H 1 .

By (3.3), the abstract form can be referred to the following eigenvalue equations in H 1 :

Δ u z p z =β(λ) u z ,
(4.2)
Δ u r p r +λ ( 1 r 0 2 κ ) u θ =β(λ) u r ,
(4.3)
Δ u θ +λκ u r =β(λ) u θ ,
(4.4)
u r r + u z z =0,
(4.5)
u z r =0, u r = u θ =0,at r=1, r 2 ,
(4.6)
u z =0, u r z = u θ z =0,at z=0,L.
(4.7)

According to [5], we take the separation of variables as follows:

u z = d h ( z ) d t d R ( r ) d r , u r = a 2 h ( z ) R ( r ) , u θ = h ( z ) φ ( r ) .
(4.8)

We utilize (4.8) in (4.5) and (4.7) to obtain:

{ h R + a 2 h R = 0 , h ( 0 ) = h ( L ) = 0 { h = cos a z , a = n π L , n = 1 , 2 , .
(4.9)

We utilize (4.8), (4.9) in (4.4) to obtain:

h φ+h φ +λκ a 2 hR=β(λ)hφ

that equals to

( d 2 d r 2 a 2 ) φ=λκ a 2 R+β(λ)φ.
(4.10)

We utilize (4.8), (4.9) in (4.2), then differentiate with respect to r,

h R + h R ′′′′ 2 p r z =β(λ) h R ,

and utilize (4.8), (4.9) in (4.3), then differentiate with respect to z,

a 2 h R+ a 2 h R +λ ( 1 r 0 2 κ ) h φ 2 p z r =β(λ) a 2 h R,

the two equations above subtract to obtain:

( d 2 d r 2 a 2 ) 2 R=λ ( 1 r 0 2 κ ) φ+β(λ) ( d 2 d r 2 a 2 ) R.
(4.11)

Utilizing (4.8) in (4.6), we have

φ(1)=φ( r 2 )=R(1)=R( r 2 )= R (1)= R ( r 2 )=0,

then let

R=sinlπ ( r 1 r 2 1 ) ,φ= C n l sinlπ ( r 1 r 2 1 ) ,l=1,2,.

Combine (4.10), (4.11) to obtain:

{ α C n l = a 2 λ κ C n l β ( λ ) , α 2 = C n l ( 1 r 0 2 κ ) λ α β ( λ ) , α= a 2 + ( l π r 2 1 ) 2 ,
(4.12)

therefore, we get

{ C n l = 2 α λ κ a 2 ± Δ , β n l = α ± Δ 2 α ,
(4.13)

where

Δ=4α a 2 λ 2 κ ( 1 r 0 2 κ ) .

The corresponding eigenvectors of (4.2)-(4.7) are as follows:

ψ n l ={ n l π 2 L ( r 2 1 ) sin n π z L cos l π ( r 1 ) ( r 2 1 ) , ( n π L ) 2 cos n π z L sin l π ( r 1 ) ( r 2 1 ) , C n l cos n π z L sin l π ( r 1 ) ( r 2 1 ) .
(4.14)

Now, as (n,l)(0,0),

L λ ψ n l = β n l ( λ ) ψ n l , E 1 = { ψ n l | ψ n l  satisfies  ( 4.14 ) , ( n , l ) N × N } , B 1 = { β n l ( λ ) | β n l ( λ )  satisfies  ( 4.13 ) , ( n , l ) N × N } ,

as n=0, from (4.2)-(4.8),

L λ ψ 0 l = β 0 l (λ) ψ 0 l ,

where

ψ 0 l E 2 = { ψ 0 l | ψ 0 l = ( 0 , 0 , sin l π ( r 1 ) ( r 2 1 ) ) T , l = 1 , 2 , } , β 0 l B 2 = { β 0 l | β 0 l = ( l π r 2 1 ) 2 , l = 1 , 2 , } ,

as l=0, ψ n 0 =0.

Now, according to the Fourier expansion, B 1 B 2 are all the eigenvalues of (4.2)-(4.7), the corresponding eigenvectors E 1 E 2 form a complete basis in H 1 .

In view of the condition μ<δ, by simple calculation, we know that κ>0, κ r 0 2 1>0, then

Δ0,

thus,

Re β n l <ε(0<ε<α),

and

β 0 l = ( l π r 2 1 ) 2 <0.

Now, by Lemma 4.3, the proof of Theorem 4.4 is completed. □

Proof of Theorem 4.5 Utilizing the method in Theorem 4.4, we have the following results for the eigenvalue problem of L λ , here L λ is adjoint operator of L λ .

As (n,l)(0,0),

{ C n l = 2 α ( 1 r 0 2 κ ) α 2 λ ± Δ , β n l ( λ ) = α ± Δ 2 α ,
(4.15)

where

Δ=4α a 2 λ 2 κ ( 1 r 0 2 κ ) ,

the corresponding eigenvectors are as follows:

ψ n l ={ n l π 2 L ( r 2 1 ) sin n π z L cos l π ( r 1 ) ( r 2 1 ) , ( n π L ) 2 cos n π z L sin l π ( r 1 ) ( r 2 1 ) , C n l cos n π z L sin l π ( r 1 ) ( r 2 1 ) ,
(4.16)

satisfying

E 1 = { ψ n l | ψ n l  satisfies  ( 4.16 ) , ( n , l ) N × N } , B 1 = { β n l ( λ ) | β n l ( λ )  satisfies  ( 4.15 ) , ( n , l ) N × N } .

As n=0,

L λ ψ 0 l = β 0 l (λ) ψ 0 l ,

where

ψ 0 l E 2 = { ψ 0 l | ψ 0 l = ( 0 , 0 , sin l π ( r 1 ) ( r 2 1 ) ) T , l = 1 , 2 , } , β 0 l B 2 = { β 0 l | β 0 l = ( l π r 2 1 ) 2 , l = 1 , 2 , } .

As l=0, ψ n 0 =0.

Consider equation β n l (λ)=0, by (4.15), we have

α ( α + 1 r 0 2 ) 2 λ 2 a 2 κ ( 1 r 0 2 κ ) = 0 , λ 2 = α 3 a 2 κ ( 1 r 0 2 κ ) .

Now, we have the critical Taylor number

T c = min ( n , l ) N 2 α 3 a 2 κ ( 1 r 0 2 κ ) ,
(4.17)

through sample calculation, as (n,l)=( L 2 ( r 2 1 ) ,1), α 3 a 2 κ ( 1 r 0 2 κ ) reaches its minimal value. We assume that n+ 1 2 L 2 ( r 2 1 ) for n>0, then there exists only one couple ( n 1 ,l), making α 3 a 2 κ ( 1 r 0 2 κ ) reach the minimal value.

Therefore, we have

β n 1 l ( λ ) = { < 0 , if  λ < T c , = 0 , if  λ = T c , > 0 , if  λ > T c , β n l ( T c ) < 0 , ( n , l ) ( n 1 , 1 ) .
(4.18)

The reduced form of (3.4) into its center manifold in H is:

d x d t = β n 1 1 x+ 1 ψ n 1 1 , ψ n 1 1 H G ( ψ , ψ ) , ψ n 1 1 H ,
(4.19)

where ψ H 1 can be written by

ψ=x ψ n 1 1 +Φ,
(4.20)

Φ is a center manifold function. G(u,v) is a bilinear operator, satisfying

G(u,v)=P ( ( u ˜ ) v z , ( u ˜ ) v r u θ v θ r 0 , ( u ˜ ) v θ + u θ v r r 0 ) .
(4.21)

Through direct calculation, we have

G ( ψ n 1 1 , ψ n 1 1 ) , ψ J H ={ 0 , ψ J ψ 2 n 1 2 , a 1 2 b 1 c 11 c 11 L ( r 2 1 ) , ψ J = ψ 2 n 1 2 ,
(4.22)

where,

a 1 = n 1 π L , b 1 = π r 2 1 , c 11 = C n 1 1 , c 11 = C n 1 1 .

Here, we use the mark as follows:

o(2)=o ( x 2 ) +O ( | β n 1 l ( λ ) x 2 | ) ,o(3)=o ( x 3 ) +O ( | β n 1 l ( λ ) x 3 | ) .

Then the center manifold function can be written by

Φ= x 2 G ( ψ n 1 1 , ψ n 1 1 ) , ψ 2 n 1 2 H β 2 n 1 2 ψ 2 n 1 2 , ψ 2 n 1 2 H ψ 2 n 1 2 +o(2).
(4.23)

Direct computation yields that

G ( ψ n 1 1 , ψ 2 n 1 2 ) , ψ n 1 1 H = 2 a 1 2 b 1 c 11 c 11 L ( r 2 1 ) , G ( ψ 2 n 1 2 , ψ n 1 1 ) , ψ n 1 1 H = a 1 2 b 1 c 11 c 11 L ( r 2 1 ) , ψ n 1 1 , ψ n 1 1 H = 1 2 c 11 c 11 L ( r 2 1 ) , ψ 2 n 1 2 , ψ 2 n 1 2 H = 8 c 11 c 11 L ( r 2 1 ) .
(4.24)

Finally, we utilize (4.20), (4.22)-(4.24) in (4.19) to obtain

d x d t = β n 1 1 xσ x 3 +o(3),
(4.25)

where

σ= 1 4 β 2 n 1 2 a 1 4 b 1 2 >0.

Now, according to Theorem 6.9 in [14], the proof of Theorem 4.5 is completed. □

References

  1. Taylor GI: Experiments with rotating fluid. Proc. R. Soc. Lond. Ser. A 1921, 100: 114-121. 10.1098/rspa.1921.0075

    Article  Google Scholar 

  2. Couette MFA: Études sur le Frottement des Liquides. Ann. Chim. Phys. 1890, 21: 433-510.

    Google Scholar 

  3. Mallock A: Experiments on fluid viscosity. Philos. Trans. R. Soc. Lond. A 1896, 187: 41-56. 10.1098/rsta.1896.0003

    Article  Google Scholar 

  4. Cloes D: Transition in circular Couette flow. J. Fluid Mech. 1965, 21: 385-425. 10.1017/S0022112065000241

    Article  Google Scholar 

  5. Chandrasekhar S: Hydrodynamic and Hydromaganetic Stability. Dover, New York; 1981.

    Google Scholar 

  6. Walowit J, Tsao S, DiPrima RC: Stability of flow between arbitrarily spaced concentric cylindrical surfaces including the effect of a radial temperature gradient. J. Appl. Mech. 1964, 31: 585-593. 10.1115/1.3629718

    Article  MathSciNet  Google Scholar 

  7. Drazin P, Reid W: Hydromagnetic Stability. Cambridge University Press, Cambridge; 1981.

    Google Scholar 

  8. Velte W: Stabilität und Verzweigung stationärer Lösungen der Navier-Stokeschen Gleichungen beim Taylorproblem. Arch. Ration. Mech. Anal. 1966, 22: 1-14. 10.1007/BF00281240

    Article  MathSciNet  Google Scholar 

  9. Kirchgässner K: Die Instabilität der Strömung zwischen zwei rotierenden Zylindern gegenüber Taylor-Wirbeln für beliebige Spaltbreiten. Z. Angew. Math. Phys. 1961, 12: 14-30. 10.1007/BF01601104

    Article  MathSciNet  Google Scholar 

  10. Kirchgässner K, Sorger P: Branching analysis for the Taylor problem. Q. J. Mech. Appl. Math. 1969, 22: 183-209. 10.1093/qjmam/22.2.183

    Article  Google Scholar 

  11. Yudovich VI: Secondary flows and fluid instability between rotating cylinders. J. Appl. Math. Mech. 1966, 30: 1193-1199. 10.1016/0021-8928(66)90081-5

    Article  Google Scholar 

  12. Ma T, Wang S: Bifurcation Theory and Applications. World Scientific, Singapore; 2005.

    Google Scholar 

  13. Ma T, Wang S: Stability and bifurcation of the Taylor problem. Arch. Ration. Mech. Anal. 2006, 181: 149-176. 10.1007/s00205-006-0415-8

    Article  MathSciNet  Google Scholar 

  14. Ma T, Wang S: Stability and Bifurcation of Nonlinear Evolution Equations. Science Press, Beijing; 2007. (in Chinese)

    Google Scholar 

  15. Taylor GI: Stability of a viscous liquid contained between two rotating cylinders. Philos. Trans. R. Soc. Lond. A 1923, 223: 289-343. 10.1098/rsta.1923.0008

    Article  Google Scholar 

  16. Pazy A: Semigroups of Linear Operators and Applications to Partial Differential Equations. Springer, New York; 1983.

    Book  Google Scholar 

  17. Adams RA: Sobolev Spaces. Academic Press, New York; 1975.

    Google Scholar 

  18. Wang M: Semigroups and Evolution Equations. Science Press, Beijing; 2006. (in Chinese)

    Google Scholar 

Download references

Acknowledgements

The authors express their sincere thanks to the referees for helpful comments and suggestions which led to the improvement of the presentation and quality of the work.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zhibo Hou.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The authors typed, read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Hou, Z., Ma, T. Dynamic phase transition for the Taylor problem in the wide-gap case. Bound Value Probl 2013, 227 (2013). https://doi.org/10.1186/1687-2770-2013-227

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-227

Keywords