_{1}

The stability of a heat-conducting flow due to the pumping of a fluid around the annulus of horizontal porous cylinders is studied. The basic flow is under the action of radial flow and a radial temperature gradient. The objects of investigations are different regimes and bifurcations which may arise in this flow.

Due to the technological and engineering importance, the investigation of the stability of curved flows has been a source of interest for several years.

Experimental studies of main flow in curved pipes carried out by Eustice [

The stability investigation of Dean flow was continued by many authors (see [

Important results were obtained also in the study of Dean flow stability of complex fluids. [

In these papers, investigations are mainly carried out in the linear approximation, which allows one to judge the first loss of stability of main flow. It was shown that the first instability of the basic steady flow can be not only a monotonic instability leading to the Dean vortices but also can be oscillatory instability resulting in a time-periodic secondary regime.

For practical purposes, it is important to consider the transitions to more complex regimes, the appearance of which precedes the development of higher instability and leads to turbulence.

For Dean flow stability, the numerical and experimental investigations of transitions to turbulence were carried out, for instance, in [

Our aim is to study possible transitions to complex modes in a heat-coducting flow between horizontal porous cylinders with radial flow and radial temperature gradient. The linear analysis of that problem for small axisymmetric perturbations has been performed in [

We consider the steady flow in an annulus between the fixed horizontal heated cylinders to be maintained by a constant azimuthal pressure gradient in the presence of a radial flow through the porous cylinders and of a radial temperature gradient. We denote the radii and temperature of the inner and outer cylinders by R 1 , R 2 and T 1 , T 2 , respectively. A constant flow through the annulus is

maintained by a suitable pumping circuit. It is assumed that the fluid inflow through one cylinder is equal to the fluid outflow through the other. The system that will be considered here is shown in

In the cylindrical coordinates ( r , φ , z ) with the z-axis coinciding with that of the cylinders, the Navier-Stokes equations as wall as those for heat transfer, continuity and state are the form [

∂ v ′ r ∂ t + ( V ′ , ∇ ) v ′ r − v ′ φ 2 r + 1 ρ ′ ∂ Π ′ ∂ r = ν ( Δ v ′ r − v ′ r r 2 − 2 r 2 ∂ v ′ φ ∂ φ ) , ∂ v ′ φ ∂ t + ( V ′ , ∇ ) v ′ φ + v ′ φ v ′ r r + 1 r ρ ′ ∂ Π ′ ∂ φ = ν ( Δ v ′ φ − v ′ φ r 2 + 2 r 2 ∂ v ′ r ∂ φ ) , ∂ v ′ z ∂ t + ( V ′ , ∇ ) v z ′ + 1 ρ ′ ∂ Π ′ ∂ z = ν Δ v ′ z , ∂ T ′ ∂ t + ( V ′ , ∇ ) T ′ = χ Δ T ′ , ∂ ρ ′ ∂ t + d i v ( ρ ′ V ′ ) = 0 , ρ ′ = ρ 0 ( 1 − β ( T ′ − T 0 ) ) , (1)

where V ′ = { v ′ r , v ′ φ , v ′ z } is the velocity vector; T ′ is temperature; Π ′ is pressure; ρ ′ is the fluid density; t is time; ν , χ , β are, respectively, the coefficients of kinematic viscosity, thermal diffusivity and thermal expansion, ρ 0 ; and T 0 denote constant density and temperature, respectively.

Operators Δ and ∇ in the cylindrical coordinates are of the form

Δ = ∂ 2 ∂ r 2 + 1 r ∂ ∂ r + 1 r 2 ∂ 2 ∂ φ 2 + ∂ 2 ∂ z 2 , ∇ = { ∂ ∂ r , 1 r ∂ ∂ φ , ∂ ∂ z } .

Equation (1) are considered under the condition

∫ 0 2 π ∫ 1 R 2 R 1 ρ ′ v ′ z r d r d φ = 0 , (2)

which ensures the absence of the fluid discharge through the cavity cross-section of the cylinders.

Let us assume that on the surface of the cylinders the following boundary conditions

v ′ r = U 0 , v ′ φ = 0 , v ′ z = 0 , T ′ = T 1 ( r = R 1 ) , v ′ r = U 0 R , v ′ φ = 0 , v ′ z = 0 , T ′ = T 2 ( r = R 2 ) , (3)

are fulfilled, where R = R 2 R 1 .

Problem (1)-(3) possesses a group of symmetries G = S O ( 2 ) ∗ O ( 2 ) which are invariant with respect to rotations L φ δ near the cylinder axis, to shifts L φ h along that axis, and also to inversion J . These factors act on the vector-function F ( v r , v φ , v z , T ) according to the rules

( L φ δ F ) ( t , r , φ , z ) = F ( t , r , φ + δ , z ) , ( L z h F ) ( t , r , φ , z ) = F ( t , r , φ , z + h ) , ( J F ) ( t , r , φ , z ) = { v r ( t , r , φ , − z ) , v φ ( t , r , φ , − z ) , − v z ( t , r , φ , − z ) , T ( t , r , φ , − z ) } (4)

for any real δ and h .

As is easily seen, problem (1)-(3) admits an exact solution

V 0 = { u 0 ( r ) , v 0 ( r ) , 0 } , T 0 = c 1 + c 2 r ϰ P r , u 0 ( r ) = R 1 U 0 r , v 0 ( r ) = ( K ϰ ( a r ϰ + 1 + b r − r ) , ϰ ≠ − 2 , K 2 ( a 1 ln r + b 1 r ) , ϰ = − 2 , ∂ Π 0 ∂ r = ρ ( u 0 2 + v 0 2 ) r , (5)

where

K = 1 2 ρ ν ( ∂ Π 0 ∂ θ ) 0 , a = R 2 − 1 ( R ϰ + 2 − 1 ) R 1 ϰ , a 1 = R 1 2 ( R 2 − 1 ) ln R ,

b = R 2 2 ( R ϰ − 1 ) R ϰ + 2 − 1 , b 1 = − R 1 2 ln R 2 − R 2 2 ln R 1 ln R ,

c 1 = T 1 R Pr ϰ − T 2 R ϰ Pr − 1 , c 2 = T 2 − T 1 R 1 ϰ Pr − 1

and ϰ = U 0 R 1 ν is the radial Reynolds number, P r = ν χ is the Prandtl number.

The radial flow is inward for ϰ < 0 (converging flow) and outward for ϰ > 0 (diverging flow).

The flow with the velocity vector V 0 , temperature T 0 and pressure Π 0 will be called the Dean stationary flow of a heat-conducting fluid with a radial temperature gradient and a radial flow.

A solution of problem (1)-(3) is sought in the form

V ′ = V 0 + V ( v r , v φ , v z ) , T ′ = T 0 + T , Π ′ = Π 0 + Π . (6)

Taking into account that the main stationary flow is a rotating shear flow, we denote rotation shear S by V m d , where V m is an average velocity in the azimuthal

direction, d = R 2 − R 1 is a gap width between cylinders. Introducing dimensionless variables for time, velocity, temperature and pressure by S , R 2 , S R 2 , T 2 − T 1 , ν ρ ′ S in the system (1), for the vector-functions F = { v r , v φ , v z , T } and F 1 = { u r , u φ , u z , T 1 } , we obtain the following nonlinear problem of finding perturbations V , T and Π :

∂ F ∂ t + N F − 1 Re M F + 1 Re ∇ 1 Π = − L ( F , F ) , ( ∇ 1 , r F ) = 0 , F | r = 1 , R = 0 , (7)

where

M F = { Δ 1 v r − 1 − ϰ r 2 v r − 2 r 2 ∂ v φ ∂ φ , Δ 1 v φ − 1 + ϰ r 2 v φ + 2 r 2 ∂ v r ∂ φ , Δ 1 v z , 1 Pr Δ 1 T } ,

N F = ω 1 ∂ F ∂ φ + { Ra ω 2 T − 2 ω 1 v φ , − g 1 v r , 0 , g 2 Pr v r } ,

L ( F , F 1 ) = { ( F , ∇ 1 ) u r − v φ u φ r , ( F , ∇ 1 ) u φ + v r u φ r , ( F , ∇ 1 ) u z , ( F , ∇ 1 ) T 1 } ,

Δ 1 = ∂ 2 ∂ r 2 + 1 − ϰ r ∂ ∂ r + 1 r 2 ∂ 2 ∂ φ 2 + ∂ 2 ∂ z 2 , ∇ 1 = { ∂ ∂ r , 1 r ∂ ∂ φ , ∂ ∂ z ,0 } ,

Re = V m d ν istheazimuthalReynoldsnumber ,

Ra = β ( T 2 − T 1 ) 2 istheRayleighnumber ,

ω 1 = V m g ( r ) r , ω 2 = ω 1 2 r , V m = K ϰ ⋅ d R 2 D 4 ,

g ( r ) = d R 2 D 1 r ϰ + 2 + D 2 − r 2 r D 4 ,

g 1 ( r ) = d v 0 d r + v 0 r = K ϰ ⋅ d R 2 D 1 ( ϰ + 2 ) r ϰ − 2 D 4 ,

g 2 ( r ) = ϰ Pr 2 ⋅ R ϰ Pr R ϰ Pr − 1 r ϰ Pr − 1 ,

D 1 = ( R 2 − 1 ) R ϰ R ϰ + 2 − 1 , D 2 = 1 − D 1 ,

D 3 = R 2 − 1 R 2 , D 4 = D 2 ln R − ϰ D 3 2 ( ϰ + 2 ) .

Note, that D e = Re ( d R 2 ) 1 2 is the Dean number.

Problem (7) is written in terms of the Boussinesq approximation, which is based on the assumption that the thermal expansion coefficient is small [

The theoretical and experimental studies have shown that after the loss of stability of Dean flow there occured secondary modes via either axisymmetric or nonaxisymmetric disturbances as Dean vortices and oscillatory modes in the form of traveling waves. To study the transition to complex regimes of special attention are the points of intersection of neutral curves, corresponding to the two above-mentioned kinds of the secondary flows, since at these points with a high probability may appear various regimes, including the complex one.

Let ( Ra 0 , Re 0 ) be the point lying on the plane of parameters ( Ra , Re ) and corresponding to the intersection of the neutral curves corresponding to the monotonic ( m = 0 ) axisymmetric and oscillatory nonaxisymmetric loss of stability of main flow (3). In this paper we will explore the complex modes that may appear in the vicinity of intersection points of neutral curves corresponding to the vortex and azimuthal wave bifurcations. Under the definite values of parameters of the problem, the neutral curves may be nonintersecting that indicates that under the corresponding values of parameters of the problem we cannot expect the appearance of complex regimes.

To construct neutral curves, we assume that the perturbations V , temperature T and pressure Π are infinitely small. Neglecting in (7) the nonlinear terms and seeking for a solution of that problem in the case of axisymmetric monotone perturbations, we obtain the following spectral problem

( M − Re N ) Φ 0 = ∇ 1 p 0 , ( ∇ 1 , r Φ 0 ) = 0 , Φ 0 | r = 1 , R = 0 , (8)

where

Φ 0 = { u 0 ( r ) , v 0 ( r ) , i w 0 ( r ) , τ 0 ( r ) } e i α z , p 0 = q 0 ( r ) e i α z . (9)

By the substitution of (9) into (8), after dividing the variables in (8) we obtain the spectral problem for the ordinary differential equations

( L − 1 − ϰ r 2 − α 2 ) u 0 = d q 0 d r − 2 Re ω 1 v 0 + Re Ra ω 2 τ 0 , ( L − 1 + ϰ r 2 − α 2 ) v 0 = − Re g 1 u 0 ,

( L − α 2 ) w 0 = α q 0 , ( L − α 2 ) τ 0 = Re g 2 u 0 , d u 0 d r + u 0 r − α w 0 = 0 , u 0 = v 0 = w 0 = τ 0 = 0 ( r = 1 , R ) , (10)

where L = d 2 d r 2 + 1 − ϰ r d d r .

The neutral curves corresponding to the bifurcation of azimuthal waves origination are found by solving the spectral problem

( M − Re N − i c Re ) Φ 1 = ∇ 1 p 1 , ( ∇ 1 , r Φ 1 ) = 0 , Φ 1 | r = 1 , R = 0. (11)

A solution of problem (11) is sought in the form

Φ 1 = { u 1 ( r ) , v 1 ( r ) , w 1 ( r ) , τ 1 ( r ) } e − i ( m φ + α z ) , p 1 = q 1 ( r ) e − i ( m φ + α z ) , (12)

where c is the unknown frequency of autooscillations, m is the azimuthal wave number. Substitution of (12) into (11) after dividing the variables in (11), we get the spectral problem for the ordinary differential equations:

[ L 1 − 1 − ϰ r 2 − α 2 − i Re ( c − m ω 1 ) ] u 1 = d q 1 d r − 2 Re ω 1 v 1 + Re Ra ω 2 τ 1 − 2 i m r 2 v 1 , [ L 1 − 1 + ϰ r 2 − α 2 − i Re ( c − m ω 1 ) ] v 1 = − i m r q 1 − Re g 1 ( r ) u 1 + 2 i m r 2 u 1 , [ L 1 − α 2 − i Re ( c − m ω 1 ) ] w 1 = − i α q 1 , [ L 1 − α 2 − i Re Pr ( c − m ω 1 ) ] τ 1 = Re g 2 u 1 , d u 1 d r + u 1 r − i m r v − i α w 1 = 0 , u 1 = v 1 = w 1 = τ 1 = 0 ( r = 1 , R ) , (13)

where

L 1 = d 2 d r 2 + 1 − ϰ r d d r − m 2 r 2 .

Problems of eigenvalues (10) and (13) have been solved by the shooting method for fixed ϰ , α , R , m , Pr , Ra . Thus, for the fixed values of these parameters we established the dependence of the critical value of the azimuthal Reynolds number Re and the neutral mode frequency c corresponding to the bifurcation of vortices and azimuthal waves origination on a Rayleigh number Ra. Further, using the Newton method, we minimize the difference between the obtained critical values of Re. This allows us to calculate with sufficient exactness the values Ra 0 , Re 0 and c 0 corresponding to the point of intersection of neutral curves.

The calculations in this paper were performed for the case R = 2 (radius of the outer cylinders is two times greater than that of the inner ones), m = 0 , 1 , for various values of axial wave number α , Pr = 0.71 and for small absolute values ϰ ( 3 ≤ ϰ ≤ 3 ). The results of calculations are presented in

As our calculations show, these intersections of neutral curves take place especially when temperature of the outer cylinder is higher than that of the inner one ( Ra 0 > 0 ) for sufficiently large values of the wave axial number and when

ϰ | α = 8 | α = 5 | α = 4 | ||||||
---|---|---|---|---|---|---|---|---|---|

Ra 0 | Re 0 | c 0 | Ra 0 | Re 0 | c 0 | Ra 0 | Re 0 | c 0 | |

3 | 3.627 | 20.9 | 2.795 | 0.714 | 68.445 | 2.473 | 0.873 | 87.84 | 2.49 |

2 | 3.14 | 22.07 | 2.792 | 0.667 | 82.88 | 2.458 | 0.809 | 81.245 | 2.472 |

1.5 | 2.52 | 24.45 | 2.7826 | 0.65 | 59.86 | 2.456 | 0.778 | 77.421 | 2.4709 |

0.5 | 2.084 | 26.938 | 2.763 | 0.64 | 53.09 | 2.473 | 0.75 | 68.464 | 2.492 |

0.2 | 1.675 | 30.2 | 2.727 | 0.65 | 51.488 | 2.48 | 0.7524 | 65.438 | 2.507 |

−0.3 | - | - | - | 0.7 | 46.25 | 2.52 | 0.772 | 60.164 | 2.541 |

−0.5 | - | - | - | 0.72 | 44.73 | 2.525 | 0.785 | 58.028 | 2.558 |

−1 | - | - | - | 0.84 | 39.73 | 2.593 | 0.835 | 52.86 | 2.606 |

−1.5 | - | - | - | 1.036 | 34.038 | 2.6678 | 0.902 | 48.268 | 2.6604 |

−3 | - | - | - | 1.784 | 25.479 | 2.877 | 1.157 | 39.24 | 2.824 |

the liquid moves through the inner cylinder to the outer one. This indicates that short-wave axially directed perturbations generating vortices and azimuthal waves are interacting, and hence one may expect the appearance of complex regimes. If values of the axial wave numbers decrease, i.e., for the corresponding long-wave perturbations there are also an intersections of neutral curves for the converging flow.

If temperature of the inner cylinder exceeds that of the outer one ( Ra 0 < 0 ) , neutral curves do not intersect, and thus it is difficult to expect the occurrence of complex modes.

To investigate the secondary flows and the appearance of high instabilities in the flow (3) we use the nonlinear theory of bifurcation of hydrodynamic flows with cylindrical symmetry (see [

This theory allows one to find various liquid motion regimes existing in the vicinity of the points of intersection of neutral curves corresponding to the two types of secondary instability-vortices and azimuthal waves in the Dean flow. This theory has been applied in a rather wide class of problems, such, for example, as the Couette flow [

The problems are reduced to the investigation of nonlinear dynamical autonomous systems of the amplitude equation, which is a generalization of Landaus amplitude equation.

The S O ( 2 ) ∗ O ( 2 ) symmetry enables one to reduce the six-dimensional amplitude system with unknown complex amplitudes η 0 , η , η 1 to the four-dimensional motor subsystem for the modules of the amplitudes with free parameters σ , μ (the damping decrements of the monotonic and oscillatory perturbations, respectively). To the equilibria of this system, which is lying on the invariant subspaces, there correspond the motions of a fluid, which has the concrete physical nature [

As our calculations shows, depending on the parameter values of the problem, the motor subsystem may have no equilibria or may have equilibria of the above mentioned types. It is found, that transition schemes may turn out to be both rather complicated and absolute trivial.

We present here the scheme of equilibria bifurcations of the motor subsystem, which we consider most interesting and allowing us to judge about the transitions characteristic of the considered system of small absolute values of radial Reynolds number ϰ .

Consider in detail the case for ϰ = 3 (dividing flow), α = 4 (perturbations are π 2 periodical in the axial direction), Ra 0 = 0.873 (temperature of the outer cylinder is higher than that of the inner one).

For σ < 0 and μ r < 0 ( μ r is the real part of damping decrement of oscillatory perturbations μ ) the equilibria fail to take place, but for μ r > 0 of all the equilibria of the motor subsystem there exist only unstable main stationary flow and unstable purely azimuthal waves.

For σ = 0 and μ r < 0 there exist only stable main stationary flow,but for μ r > 0 there are unstable main stationary flow and vortices.

For σ > 0 and μ r < 0 the unstable main stationary flow and vortices exists fon any values of free parameters; vortices are stable for μ r < μ r 2 and unstable for μ r > μ r 2 .

For σ > 0 , μ r > 0 the transitions associated with bifurcations are graphically shown in

The single lines show J-symmetric equilibria, the double lines indicate J-connected pair of equilibria. Stable equilibria are drawn by solid lines and unstable equilibria by dotted lines. The circles are the points at which the motor subsystem cycles bifurcate.

For μ r = 0 from main flow (3) there simultaneowly bifurcate unstable pure azimuthal waves and unstable J-connected pair of spiral waves. Spiral waves are stable for μ r > μ r 4 , for μ r = μ r 4 from them bifurcates an unstable connected pair of general equilibria. This pair exists in the range μ r 4 < μ r < μ r 8 and disappear for μ r = μ r 8 merging with unstable mixed azimuthal waves. On the other hand at the point μ r = μ r 9 mixed azimuthal waves merging with purely azimuthal waves.

For μ r = μ r 3 , of a set of limit symmetric cycles of the motor subsystem in the form of number eight are branching off from mixed azimuthal waves and exists

in the range μ r 3 < μ r < μ r 5 . They were calculated by using the software prepared in [

Consider the case ϰ = − 0.5 (converging flow), α = 4 (perturbations are π / 2 periodical in the axial direction), the Rayleigh number Ra 0 = 0.785 (temperature of the outer cylinder is higher than that of the inner) for σ > 0 . The transitions associated with the bifurcation are graphically shown in

In comparison to the case under consideration for diverging flow, the picture of bifurcations for converging flow becomes more complicated.

A characteristic feature of the transitions in this case are the existence of several stable equilibria in the range changes μ r , for example, μ r 2 < μ r < μ r 4 . For the values of a free parameter belonging to such ranges, in the experiments, the hysteresis phenomena are observed [

In this case, in contrast to the case ϰ > 0 for μ r = μ r 4 , from the mixed azimuthal waves branch off the nonsymmetric limit stable cycles similar to the deformed number eight, they lose their stability because of bifurcations of the period duplications at the point μ r = 5.3 and from it branches off a J-connected pair of stable II―revolution nonsymmetric cycles which exist for 5.3 < μ r < 5.7 . A further increase of μ r leads to the formation of another series of duplications of J-connected pairs of cycles which for μ r = 5.9 results in to the generation of J-connected pair of chaotic attractors; the form of their phase trajectories differs slightly from that of the trajectories of the generated cycles; changes gradually its form and disappears for μ r = 6.2 . It can be calculated in a very narrow range of variation of the parameter.In the same way, as well as with vertical cylinders (see [

In paper under consideration, we have examined possible transitions to complex modes in a heat-conducting flow produced by pumping a fluid around the annulus of a system of fixed porous horizontal cylinders. The liquid is under the action of a radial flow and a radial temperature gradient.

1) Neutral curves corresponding to vortices and azimuthal waves are calculated.

2) It is established that for diverging flow with small absolute values of radial Reynolds numbers and for sufficiently large values of the wave axial number, the intersections of neutral curves take place especially when temperature of the outer cylinder is higher than that of the inner one ( Ra 0 > 0 ) . This indicates that the short-wave axially directed perturbations generating vortices and azimuthal waves are interacting, and hence one may expect the appearance of complex regimes. If values of the axial wave numbers decrease, i.e., for the corresponding long-wave perturbations there are also intersections of neutral curves for the converging flow.

If temperature of the inner cylinder exceeds that of the outer one ( Ra 0 < 0 ) , neutral curves do not intersect, and thus it is difficult to expect the occurrence of complex modes.

1) Equilibria of a motor subsystem which are lying on the invariant subspaces are calculated. They correspond to stationary or periodic regimes of the fluid motion.

2) Equilibria of a general state are calculated. They corresponds to a quasiperiodic two-frequency solution of the amplitude system.

3) It is shown that for a converging flow there are several stable equilibria on the range of change μ r . For the values of a free parameter belonging to such ranges, in the experiments, hysteresis phenomena are observed [

4) It is established that for the certain values of the problem, there exist stable or unstable limit cycles of the motor subsystem.

5) It is shown that for converging flow, when temperature of the outer cylinder is higher than that of the inner one, the sequential duplications of a period of stable limit cycles lead to the generation of chaotic attractors and consequently as well as with vertical cylinders [

Shapakidze, L. (2017) On the Bifurcations of Dean Flow between Porous Horizontal Cylinders with a Radial Flow and a Radial Temperature Gradient. Journal of Applied Mathematics and Physics, 5, 1725-1738. https://doi.org/10.4236/jamp.2017.59146