Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2009 / Article
Special Issue

Modeling Experimental Nonlinear Dynamics and Chaotic Scenarios

View this Special Issue

Review Article | Open Access

Volume 2009 |Article ID 526945 |

J. A. P. Aranha, K. P. Burr, I. C. Barbeiro, I. Korkischko, J. R. Meneghini, "Flow Around a Slender Circular Cylinder: A Case Study on Distributed Hopf Bifurcation", Mathematical Problems in Engineering, vol. 2009, Article ID 526945, 16 pages, 2009.

Flow Around a Slender Circular Cylinder: A Case Study on Distributed Hopf Bifurcation

Academic Editor: JosƩ Roberto Castilho Piqueira
Received17 Dec 2008
Accepted06 Jan 2009
Published21 Apr 2009


This paper presents a short overview of the flow around a slender circular cylinder, the purpose being to place it within the frame of the distributed Hopf bifurcation problems described by the Ginzburg-Landau equation (GLE). In particular, the chaotic behavior superposed to a well tuned harmonic oscillation observed in the range Re > 270, with Re being the Reynolds number, is related to the defect-chaos regime of the GLE. Apparently new results, related to a Kolmogorov like length scale and the rms of the response amplitude, are derived in this defect-chaos regime and further related to the experimental rms of the lift coefficient measured in the range Re > 270.

1. Flow Around a Circular Cylinder: An Overview

Let a cylinder with a circular cross section in the plane š±=(š‘„,š‘¦) exposed to an incident flow š‘ˆš¢; if d is the circle diameter and šœˆ the kinematic fluid viscosity, the Reynolds number is defined by Re=š‘ˆš‘‘/šœˆ. Assuming a unit system where šœŒ=š‘ˆ=š‘‘=1, with šœŒ being the fluid density, the two-dimensional (2D) velocity and pressure fields, respectively {š®(š±,š‘”);š‘(š±,š‘”)}, satisfy the Navier-Stokes equations, and the boundary conditions (šœ•š‘‰š‘: cylinder cross section)

Equations (1.1) and (1.2) has a steady solution š®š‘ (š±) that is however stable only in the range Re<46.5; for Re>46.5 a limit cycle solution, oscillating with the Strouhal frequency šœ”š‘ ā‰ˆš‘ˆ/š‘‘, is observed. This limit cycle is stable in the 2D contextā€”namely, if the perturbation is restricted to the plane š±=(š‘„,š‘¦)ā€”in a large range of Reynolds numbers and Figure 1 shows typical flow visualizations in the steady (Re = 40) and limit cycles regimes (Re = 102;161).

The periodic flow in the limit cycle regime can be expanded in its harmonic components by Fourier series decomposition; namely, if š®(š±,š‘”)=š‘¢(š±,š‘”)š¢+š‘£(š±,š‘”)š£ is the flow field then

The time average š®š‘œ(š±)=š‘¢š‘œ(š±)š¢+š‘£š‘œ(š±)š£ of š®(š±,š‘”) is a flow field symmetric with respect to the š‘„-axis, with š‘¢š‘œ(š±) being an even function of š‘¦(š‘¢š‘œ(š‘„,š‘¦)=š‘¢š‘œ(š‘„,āˆ’š‘¦)) and š‘£š‘œ(š±) an odd one (š‘£š‘œ(š‘„,š‘¦)=āˆ’š‘£š‘œ(š‘„,āˆ’š‘¦)); the first harmonic š®1(š±)=š‘¢1(š±)š¢+š‘£1(š±)š£ is an anti-symmetric field, with {š‘¢1(š‘„,š‘¦)=āˆ’š‘¢1(š‘„,āˆ’š‘¦);š‘£1(š‘„,š‘¦)=š‘£1(š‘„,āˆ’š‘¦)} and, as a rule, one can show for a circular cylinder that the even harmonics are symmetric and the odd ones anti-symmetric: the in-line force (drag) depends thus only on the even modes while the transverse force (lift) depends only on the odd modes. Figure 2 displays, for Re=100, the functions {(š‘¢1,š‘(š±),š‘£1,š‘(š±));(š‘¢2,š‘(š±),š‘£2,š‘(š±))} obtained both from numerical simulation of the 2D flow and from PIV measurement of an actual flow experiment. Both plots are visually very similar and confirm the symmetry/anti-symmetry behavior quoted above; furthermore, the agreement between them is also quantitative: defining the normalized mode amplitude by the ratio š“š‘›(Re) =max|š‘¢š‘›(š±)|/max|š‘¢š‘œ(š±)|, Figure 3 displays the function š“š‘›(Re) determined numerically in the range 60 ā‰¤ Re ā‰¤ 600 and also the same value obtained experimentally at Re = 100.

It must be observed also that these numerical results indicate a hierarchy between the modes amplitudes and a square root singularity near š‘š‘1ā‰…46.5, namely

typical of a Hopf supercritical bifurcation. This point will be elaborated in the following; in fact, if one writes (1.3) in the complex form, placing (1.5) into (1.1) and separating the harmonic parcels {exp (inšœ”š‘ š‘”); š‘›=0,1,2,ā€¦} one obtains the sequence of problems, the first one, that determines š®š‘œ(š±), being nonlinear and the remaining ones linear, as usual in an asymptotic expansion. In fact, if higher order terms in the ā€œsmall parameterā€ šœ€ are disregarded, one may express, to leading order, the field š®1(š±) in the form with {šœ†av=š‘–šœ”š‘ ;šžav(š±)} being the eigenvalueā€”eigenvector of the homogeneous problem defined in (1.7): this is consistent with a numerical result due to Barkley [2], stating that the averaged flow š®š‘œ(š±) is marginally stable (Real šœ†av=0) with respect to 2D perturbation.

But this is not enough for the present purpose: the final goal is to solve the tri-dimensional (3D) problem for a slender cylinder having to solve basically the 2D cross section problem: besides the obvious economy in the degrees of freedom needed in the numerical computation, the 2D flow is well organized (laminar) while the 3D one is chaotic (turbulent), as it will be seen later in this paper. As discussed in Aranha [3], the flow around a slender cylinder can be asymptotically approximated by the Ginzburg-Landau equation but one needs then, first of all, to express the harmonic mode š®1(š±)ā‹…exp(š‘–šœ”š‘ š‘”) in the form š‘Ž(š‘”)ā‹…šž(š±), as in (1.7), with š‘Ž(š‘”) satisfying Landauā€™s equation

The hope is that such š‘Ž(š‘”)ā‹…šž(š±), with the related eigenvalue-eigenvector {šœ†(Re)=šœŽ + š‘–šœ”;šž(š±;Re)}, coalesce with the standard Hopf bifurcation in the limit šœ€ā†’0 (Reā†’š‘š‘1), while recovering the Fourier series expansion (1.5) when šœ€2 is ā€œsmallā€ but finite in the range Reā‰«š‘š‘1. This double requisite obliges to look for a basic stationary flow that is neither the steady solution š®š‘ (š±), useless far from the bifurcation, nor the averaged flow š®š‘œ(š±), always marginally stable and so unsuited to describe a Hopf bifurcation.

A clue is given by the following observation: the steady state solution š®š‘ (š±;Re), that becomes unstable at Re=š‘š‘1, satisfies the homogeneous 0th-order equation (1.6); with šŸš‘œ(š±) defined in (1.6), if one considers instead of š®š‘ (š±;Re) the field the stability of this flow coalesce with the one related to the steady state solution š®š‘ (š±;Re) with an error of order šœ€4 in the limit šœŽā†’0(Reā†’š‘š‘1): in Hopf bifurcation one has šœŽ(Re)=š›¼ā‹…(1āˆ’š‘š‘1/Re)ā‰ˆšœ€2 when Reā†’š‘š‘1 and this relation is recovered if the field defined in (1.9) is used, instead of the standard steady state solution š®š‘ (š±;Re), as the basic field.

At š‘š‘1 one has strictly š®š‘ (š±;š‘š‘1)ā‰”š®š‘œ(š±;š‘š‘1); however, as Re increases the steady state solution š®š‘ (š±;Re) presents a bulbous region in the wake, similar to the one indicated in Figure 1 but with a length increasing linearly with Re: the difference between š®š‘ (š±;Re) and š®š‘œ(š±;Re) becomes enormous in the range Re > 200, in despite of the fact that the forcing term šŸš‘œ(š±) in the problem that defines š®š‘œ(š±;Re) be small, of order šœ€2. This apparent paradox is in fact due to an extreme susceptibility of the steady flow š®š‘ (š±;Re) to the influence of ā€œsmall forcesā€, either applied directly, as šŸš‘œ(š±), or else indirectly, as the constraint forces that appear on the outer contour of the finite domain used in the numerical computation; for example, to determine numerically š®š‘ (š±) at Re = 600 with reasonable accuracy one needs to discretize a circle with radius 1000d: only then the ā€œsmall constraint forcesā€ in the outer circle becomes small enough to not impair convergence. In the other hand, the presence of the small forcing term šŸš‘œ(š±) seems to regularize the problem, since then the time average field š®š‘œ(š±;Re) is robust: it can be easily determined numerically, without any major concern about the region size to be discretized, and it also changes weakly with the Reynolds number.

In the asymptotic solution that leads to Landau's equation (1.8) terms of order šœ€4 are ignored and the fields {Ģ‚š®š‘œ;Ģ‚š‘š‘œ} can be thus determined by solving the regular linear system where the term (š›æš®ā‹…āˆ‡)š›æš®ā‰…š’Ŗ(šœ€4) was disregarded; notice that the linear operator (1.10) is regular since its eigenvalue šœ† with largest real part is given by {Realšœ†=0; Imag šœ†=šœ”š‘ ā‰…š’Ŗ(1)}, see (1.7).

The eigenvalue-eigenvector {šœ†(Re)=šœŽ+š‘–šœ”;šž(š±;Re)} corresponding to the basic flow defined in (1.9) is solution of the problem and since {Ģ‚š®š‘œāˆ’š®š‘œ;Ģ‚š‘š‘œāˆ’š‘š‘œ}ā‰…š’Ŗ(šœ€2), comparing (1.7) to (1.11) one obtains

Observing that the harmonic components {š®š‘›(š±);š‘›>0} in (1.5) tend to zero as ā€–š±ā€–ā†’āˆž, and so it does the functions {šŸš‘œ(š±);Ī”šŸš‘œ(š±)}, one considers now the solution of satisfying the same boundary conditions (1.2). The steady state solution of (1.13), defined in (1.9), becomes unstable for Re>š‘š‘1, the only unstable mode being given by {šœ†(Re)=šœŽ+š‘–šœ”;šž(š±;Re)}, solution of (1.11); the solution of (1.13) in the unstable range Re>š‘š‘1 can be thus expressed by means of the standard asymptotic series where (āˆ—) stands for the complex conjugate of the expression in the left and, as usual, the mode amplitude š‘Ž(š‘”) is assumed to change slowly in time, the slow time being proportional to the amplitude square.

By placing (1.14) into (1.13) and separating terms of like orders in šœ€, a sequence of linear problems is obtained, allowing to compute the fields {š®š›¼š›½(š±);š‘š›¼š›½(š±)}. Details will be omitted here but two points must be commented. First, the operator that determines š®31(š±) is exactly the one defined in (1.11) and it is thus singular: the solvability condition (Fredholm alternative) of this š®31-equation leads to Landau's equation (1.8); second, for future reference, the field š®20(š±) is solution of the equation where the relation (Ģ‚š®š‘œāˆ’š®š‘œ)ā‹…|š‘Ž(š‘”)|2š®20ā‰…š’Ŗ(šœ€4) was used.

The 2D systems (1.1) and (1.13) have both the same singularity at š‘š‘1ā‰…46.5 and are both regular in ā€œall rangeā€ Re>š‘š‘1, a result numerically confirmed by Henderson [4] up to Re = 1000; since one system differ from the other only by a forcing term of order šœ€4, one should have asymptotically a result that will be explored next. In fact, recalling that {Ģ‚š®š‘œāˆ’š®š‘œ;Ģ‚š‘š‘œāˆ’š‘š‘œ}ā‰…š’Ŗ(šœ€2) one has, with the help of (1.16),

Two results can be derived directly from the latter equality (see also (1.9) and (1.15)), and thus it follows from (1.14) that the asymptotic solution of (1.13), based on the Landauā€™s equation (1.8), recovers the observed 2D periodic (limit cycle) solution (1.5) in the range Reā‰«š‘š‘1, with an error of order šœ€4. Or in short: Landau's equation (1.8), strictly valid in a close neighborhood of a Hopf supercritical bifurcation, can be extended in the present flow problem to the range Reā‰«š‘š‘1, where a neat periodic solution persists.

Finally, once the 2D numerical solution š®(š±,š‘”) is determined and its harmonic components {š®š‘›(š±);š‘›=0,1,2,3} are computed, the unstable mode {šœ†=šœŽ+š‘–šœ”;šž(š±);š‘š‘’(š±)} and the coefficients {šœ‡;š‘3} of Landauā€™s equation can be easily estimated, as elaborated below. In fact, using the approximations and the normalization of the mode šž(š±), multiplying (1.11) by šžāˆ—(š±), integrating by parts and using āˆ‡ā‹…šžāˆ—=0, one obtains with āˆ‡šžāˆ¶āˆ‡šžāˆ—=āˆ‡š‘’š‘„ā‹…āˆ‡š‘’āˆ—š‘„+āˆ‡š‘’š‘¦ā‹…āˆ‡š‘’āˆ—š‘¦. Using again the approximation šžavā‰…š®1(š±)/|š‘Ž|ā‰…šž, with the same error šœ€2 in (1.7), the following identity can be derived, and subtracting one expression from the other, while using Ģ‚š®š‘œ=š®š‘œāˆ’š›æš®, one obtains with š›æš®(š±;Re)ā‰…š’Ŗ(šœ€2) being solution of the regular linear system (1.10).

Notice that (1.22) reaffirms, as it should, the estimated orders in (1.12) and placing the š‘Ž(š‘”) defined in (1.18) into Landau's equation (1.8) one has

Summarizing: through the 2D simulation one obtains {š®(š±,š‘”);š‘(š±,š‘”)} and from the Fourier expansion in the harmonics of the observed frequency šœ”š‘  one determines the averaged flow {š®š‘œ(š±;Re);š‘š‘œ(š±;Re)} and the first harmonic {š®1(š±;Re);š‘1(š±;Re)} defined in (1.5). Solving the linear system (1.10) the field š›æš®(š±;Re) can be computed and so the coefficients of Landau's equation using (1.19), (1.22), and (1.23). The gain in this extra computation is certainly marginal in the context of the 2D problem; however, as it will be discussed in the following sections, Landau's equation is the basis of the 3D Ginzburg-Landau equation (GLE) and with it one can possible predict an asymptotic approximation of the 3D behavior without having to resort to a 3D numerical computation of the flow field. In this context, the proposed approximation is similar to existing ā€œslender body theoriesā€ in applied mechanics, as for example the Lifting Line Theory in the Aerodynamics of slender wings: in all of them one takes profit of the body slenderness to correct asymptotically the 2D solution. But before one addresses this 3D extension of Landauā€™s equation it is worth to mention some general features of the actual 3D flow around a slender cylinder.

2. Features of the 3D Flow Around a Slender Cylinder

The 2D flow around a slender cylinder is unstable with respect to 3D-perturbation for Re > 190 and this is instability, known experimentally for a long time, has only recently been verified theoretically in a comprehensive numerical study done by Henderson [4]. The plot of the Strouhal number St = š‘“š‘ š‘‘/š‘ˆ(š‘“š‘ =šœ”š‘ /2šœ‹) as a function of Re, see Figure 4, portrays this instability in a very clear way and Henderson [4] has shown that the bifurcation at š‘š‘2ā‰…190 is subcritical while a second one at š‘š‘3ā‰…260 is supercritical. The curve St(Re) presents a hysteretic behavior in the range 180<Reā‰¤260, where two competing solutions, corresponding to two distinct attractors, can appear depending on the initial conditions; as usual in a subcritical bifurcation, the presence of the two attractors can be detected a little before the critical point š‘š‘2ā‰…190.

In the range Re > 260 the Strouhal frequency changes weakly with Re and the flow pattern presents a well tuned frequency immersed in a chaotic (turbulent) background. In this range of Reynolds numbers the most conspicuous experimental result is, certainly, the ā€œlift crisisā€ observed by Norberg [5] and briefly commented below.

The sectional transverse (lift) force š‘™(š‘§,š‘”) was measured by Norberg [5] in the range 250 < Re < 10ā€‰000 and the rms of the lift coefficient šœš‘™(š‘§,š‘”)=š‘™(š‘§,š‘”)/((1/2)šœŒš‘ˆ2š‘‘) was plotted as a function of Re, the obtained result being shown in Figure 5. The ā€œlift crisisā€ corresponds to the sharp drop of rms šœš‘™ at Re ā‰ˆ 260, reaching a minimum at Re ā‰ˆ 1000 of about 20% of the 2D value and there remaining up to Re ā‰ˆ 5000, where the value ofā€‰ā€‰rms šœš‘™ starts a slow recovering.

The behavior is similar to the well known ā€œdrag crisisā€ in the range 105< Re < 106, although even sharper, and it should be also related to the chaotic (turbulent) flow observed when Re > 260. The main purpose in the present paper is to indicate that the Ginzburg-Landau Equation (GLE) has the potential ability to recover Norbergā€™s ā€œlift crisisā€ and this point will be addressed next.

3. Ginzburg-Landau Equation in the Defect Chaos Regime

The 2D unstable mode š‘Ž(š‘”)ā‹…šž(š±)exp(š‘–šœ”š‘”) is triggered by a random perturbation distributed along the cylinderā€™s span and one should expect, as a consequence, a certain phase-lag of this mode in the z-direction: the amplitude a must then change with the span coordinate š‘§, namely, š‘Ž=š‘Ž(š‘§,š‘”). The z-dependence of the mode amplitude should modify the 2D Landauā€™s equation by a parcel proportional to a z-derivative of a and observing that there is no preferred z-direction this derivative should be even in z: the obvious choice here is to take the second derivative šœ•š‘§š‘§š‘Ž. This is perhaps the simplest argument to introduce the Ginzburg-Landau Equation (GLE), as done by Ginzburg in 1950 in his joint study with Landau on super-conductivity, see Ginzburg [6]. It was introduced then as a phenomenological model, namely, as an equation that emulates the overall behavior of an observed phenomenon, and as such has been used in Physics, see Aranson and Kramer [7], to analyze a class of problems related to a distributed Hopf bifurcation; the flow around a slender cylinder is just an example of it.

In this context, the GLE was first proposed as a phenomenological model by AbarĆØde and Monkewitz [8] and studied by Monkewitz and co-authors in several papers; particularly interesting is the work by Monkewitz et al. [9] where some subtle aspects of the flow are theoretically predicted and confirmed experimentally. These works were restricted, however, to the range Re < 160, within the stable range of the 2D periodic flow, and the purpose here is to extend it to the unstable regime Re >š‘š‘2ā‰… 190.

Normalizing time, space and amplitude by using {š‘”ā†šœŽš‘”; š‘§ā†(šœŽ/š›¾)1/2š‘§; š‘Žā†(šœŽ/šœ‡)1/2š‘Ž} the same equation (3.1) is obtained with šœŽ=š›¾=šœ‡=1: the behavior of the GLE depends only on the dispersion coefficients {š‘1;š‘3} and it is not difficult to show, via Fourier Transform of the perturbed equation, that the 2D solution becomes unstable with respect to 3D perturbation when š‘1ā‹…š‘3>1; incidentally, this stability condition is usually called the Benjamin-Feir condition, in honour of a stability study in water waves done by these authors, see Benjamin and Feir [11]. In Figure 6 it is shown the results of a comprehensive numerical study done by Shraiman et al. [10] in the unstable region of the dispersion plane (š‘1,š‘3). It discloses, at first, two distinct chaotic regimes: one, very mild, called ā€œphase chaosā€, is characterized by a ā€œturbulenceā€ superposed on the uniform 2D phase and restricted to a small strip on the unstable region š‘1ā‹…š‘3>1; the other, very energetic and covering the remaining of the unstable region, called ā€œdefect chaosā€, is related directly to the amplitude size |š‘Ž|: the ā€œdefectsā€ are the points in time-space plane (š‘§,š‘”) where the amplitude is null and the iso-phases either stop or bifurcate at them, see the plots of the iso-phases in the detached figures in Figure 6. Shraiman et al. [10] also observed a thin strip, coined bi-chaotic, close to the threshold curve š‘1ā‹…š‘3=1 and in fact penetrating a little into the stable region š‘1ā‹…š‘3<1, where the GLE has two chaotic attractors.

It seems then that the GLE, with recognized predictive ability in the stable range (Re < 190 or š‘1ā‹…š‘3<1), may be useful also in the unstable range (Re > 190 or š‘1ā‹…š‘3>1) since, as in the flow problem, it presents a bi-chaotic behavior in the vicinity of the threshold point (Re ā‰ˆ 190 or š‘1ā‹…š‘3ā‰ˆ1) and a chaotic one when Re ā‰« 190 or š‘1ā‹…š‘3ā‰«1. The difficulty here is first of all operational, since it seems awkward to adjust the parameters of the phenomenological GLE to the empirical data of the now chaotic flow, and also conceptual in some sense, once it is understood that GLE can model the problem just in the vicinity of Hopf bifurcation but not far from it, although Monkewitz et al. [9] used GLE to model properly the flow problem at a Re almost three times larger than critical Reynolds š‘š‘1ā‰… 46.5.

However, as seen in the first section, Landauā€™s equation can be extended far beyond bifurcation and the GLE can be obtained as an asymptotic approximation of the 3D flow related to (1.13), allowing one to determine the curve (š‘3(Re); š‘1(Re)) representing the flow problem. This computation was not done yet, however, only scarce results are available; meanwhile, it seems interesting to check whether or not GLE has the potential ability to recover the main features of the observed 3D flow. For example, one certainly should expect that the curve (š‘3(Re); š‘1(Re)) crosses the Benjamin-Feir curve š‘1ā‹…š‘3 = 1 at Re ā‰ˆ 190, penetrating after the defect chaos regime trough the bi-chaotic region, corresponding to the hysteretic behavior in the range 180 < Re ā‰¤ 260 observed in Figure 4; as Re rises above 260 the curve (š‘1(Re); š‘2(Re)) must go even deeper into the defect chaos regime and, in particular, Norberg's lift crisis must be predicted if the GLE approximation is consistent. But the transverse (lift) force is due to the odd harmonics, and so it is proportional to š‘Ž(š‘§,š‘”): the sharp drop in rms šœš‘™ must be related, in the GLE context, to a sharp drop in rms |š‘Ž(š‘§,š‘”)|. The purpose here is to discuss this point while revealing some interesting aspects of the GLE in the defect chaos regime, that may have an interest in itself.

As in a turbulent flow regime, the chaotic solution in the ā€œdefect chaosā€ regime is characterized by a cascade of length scales š‘˜āˆ’1 limited below by a ā€œKolmogorov scaleā€ (š‘˜(kol))āˆ’1, where the dissipated power, proportional to š›¾ā‹…|šœ•š‘§š‘Ž|2, is of order of the power given by the instability, proportional to šœŽā‹…|š‘Ž|2; it follows that

Equation (3.1) was integrated in the region 0ā‰¤š‘§ā‰¤š‘™=1000 in the time interval 18000ā‰¤š‘”ā‰¤20000, using the periodic boundary condition {š‘Ž(0,š‘”)=š‘Ž(š‘™,š‘”);šœ•š‘§š‘Ž(0,š‘”)=šœ•š‘§š‘Ž(š‘™,š‘”)}. Figure 7 shows the wavenumber spectrum š‘†(š‘˜) of š‘Ž(š‘§,š‘”) for pairs of values (š‘3;š‘1) and it is clear that the energy is almost exhausted in the region š‘˜>š‘˜kolā‰ˆ1. This behavior was observed in all numerical experiments in the grid {1ā‰¤š‘3ā‰¤20;1ā‰¤š‘1ā‰¤20;š‘3ā‹…š‘1>1}.

The intensity of the response can be also estimated by the wavenumber spectrum integral, and in Figure 8(a) the values of š¼[š‘†(š‘˜)] and rms |š‘Ž(š‘§š‘š,š‘”)|, š‘§š‘š=(1/2)š‘™, are plotted again for several points in the ā€œdispersion planeā€ (š‘3;š‘1). The almost exact agreement between the two plots indicates that the random signal š‘Ž(š‘§,š‘”) is weakly stationary, namely

Figure 8(a) shows that the rms of |š‘Ž(š‘§,š‘”)| decreases monotonically with š‘3, kept š‘1 constant, but when š‘3 is constant it increases with š‘1, also monotonically in the range š‘3>4. This behavior can be inferred from an identity of the GLE. In fact, if (3.1) is multiplied by š‘Žāˆ— and integrated in the interval 0ā‰¤š‘§ā‰¤š‘™, one obtains, after using the periodicity of the boundary conditions and the weak stationary condition (3.4), the identities

Now, if š‘Ž(š‘§,š‘”)=|š‘Ž(š‘§,š‘”)|ā‹…exp[š‘–šœ‘(š‘§,š‘”)] and introducing the average frequency šœ”š‘š by the expression one obtains from (3.5)

This relation was obtained under the weak stationary assumption (3.4) and it seems reasonable to assume that the intensity of rms |š‘Ž(š‘§,š‘”)| can be gauged by šœ’; notice, in particular, that šœ’ is monotonically increasing with š‘1 when šœ”š‘š/š‘3<1 and decreasing with š‘3 increasing, in accordance to the observed in Figures 8(a), 8(b) for rms |š‘Ž(š‘§,š‘”)|. From the relation šœ’ā‰…š‘Ÿš‘šš‘ |š‘Ž(š‘§,š‘”)| it follows also the asymptotic relations

The expression (ii) in (3.8) can be related to the Kolmogorov scale (3.2). In fact, lets recall, first of all, a standard result: by assuming an harmonic wave solution š‘Ž(š‘§,š‘”)=|š‘Ž|ā‹…exp[š‘–(š‘˜ā‹…š‘§+šœ”ā‹…š‘”)] of the GLE (3.1) one obtains the dispersion relation, depending on the ā€œdispersion coefficientsā€ (š‘1;š‘3). For a ā€œrandom waveā€ one may take (š‘Ÿš‘šš‘ š‘Ž(š‘§,š‘”))2 in the place of |š‘Ž|2 in (3.9) and if š‘˜=š‘˜(kol)ā‰ˆ1 one obtains, with the help of (ii) in (3.8), šœ”ā‰ˆšœ”š‘š, or in short: the averaged frequency defined in (3.6) is the ā€œKolmogorov frequency scaleā€ of the random signal š‘Ž(š‘§,š‘”) in the limit š‘3ā†’āˆž, kept š‘1 constant; in this limit šœ”š‘š tends to a bounded value šœ”āˆž(š‘1), see Figure 8(b). The data of Figure 8(a) confirm, in the limit {š‘3ā†’āˆž;š‘1=1}, the asymptotic behavior š‘Ÿš‘šš‘ š‘Ž(š‘§,š‘”)ā‰…šœŒā‹…š‘3āˆ’1/2 with šœŒā‰ˆ(š‘1+šœ”š‘š)1/2; in reality, šœŒā‰…1.49 from the data of Figure 8(a) while (š‘1+šœ”š‘š)1/2ā‰… 2.17 from Figure 8(b).

One expects then that š‘Ÿš‘šš‘ |š‘Ž(š‘§,š‘”)| diminishes monotonically with increasing š‘3, a result confirmed by the direct evaluation of š‘Ÿš‘šš‘ |š‘Ž(š‘§,š‘”)| in the dispersion plane (š‘3,š‘1), see Figure 9; notice that expression (i) in (3.8) is also recovered, a result consistent with the ā€œphase chaosā€ regime identified in Figure 6.

In the flow problem, the linear dispersion coefficient š‘1 is not expected to change too much with Re but š‘3, defined by the ratio (šœ”š‘ āˆ’šœ”)/šœŽ, see (1.23), apparently does: the difference (šœ”š‘ āˆ’šœ”) is small but fairly constant while Ļƒ appears to drop sharply for Re above 100, the ratio (šœ”š‘ āˆ’šœ”)/šœŽ=š‘3 becoming very large then: as it was seen, if š‘3ā‰«1 then š‘Ÿš‘šš‘ |š‘Ž|ā‰Ŗ1 and thus š‘Ÿš‘šš‘ šœš‘™ā‰Ŗ1. This result must be confirmed by a more refined numerical solution but it indicates, anyway, the ability of the GLE to predict Norberg's ā€œlift crisis.ā€ Or, in other words, if the actual curve (š‘3(Re); š‘1(Re)) in fact penetrates the unstable range through the bi-chaotic region at Reā‰ˆš‘š‘2ā‰…190 and š‘3(Re) increases rapidly with Re, then the GLE, together with the asymptotic expansion (1.14), defines in fact a reduced Navier-Stokes Equation for the flow around a slender cylinder, the practical importance of it being commented below.

4. Conclusion

In this paper the possibility to solve asymptotically the flow around a slender cylinder using a 2D computation and the Ginzburg-Landau equation to obtain the 3D correction was elaborated, stressing the regime above Re ā‰… 190, where three dimensionality has a marked influence. Although one must wait more refined numeric results to reach a definitive conclusion, the qualitative behavior of GLE in the range š‘3ā‹…š‘1>1 matches very well the most important qualitative features of the flow around a slender cylinder in the range Re > 190; as already discussed in Monkewitz et al. [9], the matching between both is impressive in the range Re < 160.

A practical problem where the present study may be relevant is related to the fatigue analysis of ā€œrisersā€ (vertical ducts) in the offshore oil production systems, essential to assure the safe operation of these systems during its projected life: risers are exposed to ocean currents and oscillate transversally in the elastic modes with natural frequencies close to the flowā€™s Strouhal frequency, the related cyclic stress causing fatigue of the material. The use of a 3D Navier-Stokes code to obtain practical answers is, however, completely out of question in the present stage of development, not only due to computer time needed, but also for the lack of confidence in the numerical results of the enormous discrete system related to it. The reduced Navier-Stokes equation, represented by the GLE, opens an opportunity to a feasible and relatively cheap computation: in it, the complex coupling between the incoming flow and the riser's elasticity can be represented by a coupled set of equationsā€”one of them being the (extended) GLE, the other representing the riserā€™s elastic behaviorā€”both depending only on the space variable along the riser's span, turning the discrete model orders of magnitude smaller. This is the main motivation to study this problem at NDF.


The authors acknowledge the financial support from FINEP-CTPetro, FAPESP, PETROBRAS and CNPq.


  1. M. van Dyke, An Album of Fluid Mechanics, The Parabolic Press, Stanford, Calif, USA, 1982.
  2. D. Barkley, ā€œLinear analysis of the cylinder wake mean flow,ā€ Europhysics Letters, vol. 75, no. 5, pp. 750ā€“756, 2006. View at: Publisher Site | Google Scholar | MathSciNet
  3. J. A. P. Aranha, ā€œWeak three dimensionality of a flow around a slender cylinder: the Ginzburg-Landau equation,ā€ Journal of the Brazilian Society of Mechanical Sciences and Engineering, vol. 26, no. 4, pp. 355ā€“367, 2004. View at: Publisher Site | Google Scholar
  4. R. D. Henderson, ā€œNonlinear dynamics and pattern formation in turbulent wake transition,ā€ The Journal of Fluid Mechanics, vol. 352, pp. 65ā€“112, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  5. C. Norberg, ā€œFluctuating lift on a circular cylinder: review and new measurements,ā€ Journal of Fluids and Structures, vol. 17, no. 1, pp. 57ā€“96, 2003. View at: Publisher Site | Google Scholar
  6. V. L. Ginzburg, On Superconductivity and Superfluidity, Nobel Lecture, Springer, Berlin, Germany, 2003.
  7. I. S. Aranson and L. Kramer, ā€œThe world of the complex Ginzburg-Landau equation,ā€ Reviews of Modern Physics, vol. 74, no. 1, pp. 99ā€“143, 2002. View at: Publisher Site | Google Scholar | MathSciNet
  8. P. AbarĆØde and P. A. Monkewitz, ā€œA model for the formation of oblique shedding and ā€œchevronā€ patterns in cylinder wake,ā€ Physics of Fluids A, vol. 4, no. 4, pp. 744ā€“756, 1992. View at: Publisher Site | Google Scholar
  9. P. A. Monkewitz, C. H. K. Williamson, and G. D. Miller, ā€œPhase dynamics of KĆ”rmĆ”n vortices in cylinder wakes,ā€ Physics of Fluids, vol. 8, no. 1, pp. 91ā€“96, 1996. View at: Publisher Site | Google Scholar | MathSciNet
  10. B. I. Shraiman, A. Pumir, W. van Saarlos, P. C. Hohenberg, H. ChatĆ©, and M. Holen, ā€œSpatiotemporal chaos in the one-dimensional complex Ginzburg-Landau equation,ā€ Physica D, vol. 57, no. 3-4, pp. 241ā€“248, 1992. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  11. T. B. Benjamin and J. E. Feir, ā€œThe desintegration of wave trains on deep waterā€”part I: theory,ā€ The Journal of Fluid Mechanics, vol. 27, no. 3, pp. 417ā€“430, 1967. View at: Publisher Site | Google Scholar

Copyright © 2009 J. A. P. Aranha et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles