Non Linear Vortex Structures in Stratified Fluid Driven by Small-Scale Helical Force

In this work, we consider the effect of a small-scale helical driving force on fluid with a stable temperature gradient with Reynolds number . At first glance, this system does not have any instability. However, we show that a large scale vortex instability appears in the fluid despite its stable stratification. In a non-linear mode this instability becomes saturated and gives a large number of stationary spiral vortex structures. Among these structures there is a stationary helical soliton and a kink of the new type. The theory is built on the rigorous asymptotical method of multi-scale development. 1 R


Introduction
The importance of the generation processes of large-scale coherent vortex structures in hydrodynamics is well known.A large-scale vortex means a vortex which is generated by a much smaller scale force or in turbulence with a characteristic scale much smaller than a vortex scale.When these coherent structures appear in smallscale turbulence, they play a key role in transfer processes (see for instance [1]).Numerical and laboratory experiments [2][3][4][5][6][7], confirm the existence of coherent vortex structures, especially for two-dimensional or quasi two-dimensional turbulence [7][8][9].Notably, they are well observed in geophysical hydrodynamics like various cyclones in the planet's atmospheres [10,11].Sometimes the appearance of large scale vortex structures is accompanied by the inverse cascade of energy both in the three-dimensional case (AKA-effect [12]), and in quasitwo-dimensional cases as well [3][4][5][6]8,9].It may be said that the inverse cascade itself is also one of the mechanisms of the generation of large-scale structures [4,13].The generation of large scale slow movements by small scale external forces in a rotating stratified fluid was also studied numerically in works [14,15].One of the important large scale instabilities in non compressible fluid is the AKA-effect (Anisotropic Kinetic Alpha effect) which was found in work of Frisch, She and Sulem [16].In this work, a large scale instability appears under the impact of small scale force in which parity is broken (with zero helicity).In the following work [17] the inverse cascade of energy and the non linear mode of instability saturation were studied.Despite the fact that the broken parity is a more general notion than helicity, it is the helicity rot 0   v v which is the widespread mechanism of broken parity in hydrodynamical flow.For instance, the turbulence becomes helical when rotation and stratification are taken into account [18][19][20].Therefore one may consider the small-scale helical force the parametrization of this turbulence.The injection of a helical external force into hydrodynamic systems was considered in many works ( [21][22][23][24]), and as a result it was understood that a small-scale turbulence able to generate large-scale perturbations cannot be simply homogeneous, isotropic and helical [25], but must have additional special properties.In some cases, the existence of largescale instability was shown (vortex dynamo or hydrodynamic  -effect).(In the magneto hydrodynamics of conductive fluid the  -effect is well known [26] ).In particular, in work [22] it is shown that large-scale instability exists in convective systems with small-scale helical turbulence.These works as well as the results of numerical modelling are described in detail in review [27], which are focused essentially on the possible application of these results to the issue of tropical cyclone origination.
In this work we consider the theory of large scale vortex generation in stratified fluid under the impact of small scale helical force.Let us suppose that there is a stable temperature stratification in fluid.To this fluid with the Reynolds number let us apply a smallscale, helical, external force.This force will maintain in the fluid small-scale helical fluctuations of velocity field We consider the fluid as being boundless.At first glance there are no instabilities at all in this system.However, we show in this work that despite stable stratification, a large-scale vortex instability appears in the fluid which leads to the generation of large-scale vortex structures.The theory of this instability is built rigourously using the method of asymptotical multi-scale development similar to what was done in work of Frisch, She and Sulem for the theory of AKA-effect [16].But the equations which we solve differ considerably from equations in work [16].In addition to linear theory, we also develop and study in details the non-linear theory of this instability saturation.We devote special attention to stationary, non-linear, periodical vortex structures which appear as a result of the saturation of found instability.Among these structures, there is a spiral vortex soliton and kink of the new type.In order to distinguish our instability from others in stratified fluid we consider the case of stable stratification.Nevertheless our theory permits the examination of unstable stratification as well by means of substitution However, in this case we have to consider that the usual convective instability is eliminated and the Raleigh number is reasonably small. .

Ra Ra  
Our work is arranged as follows: in Section 2 we set forth the formulation of the problem and equations for stable stratification in Boussinesq approximation; in Section 3 we examine the principal scheme of multi scale development and we give secular equations.In Section 4 we describe external force properties and calculate the Reinolds stress.In Section 5 we discuss the non-linear stage of the instability and its saturation.We study the equations of non-linear instability and its stationary solutions.It is shown that due to the hamiltonian nature of these equations a large number of stationary vortex structures of spiral type appear.We also demonstrate that there are solutions in the form of the spiral soliton and the kink of new type.The obtained results are discussed in the conclusion in Section 6.

Main Equations and Formulation of the Problem
Let us consider the equations for the motion of non compressible fluid with a constant temperature gradient in the Boussinesq approximation: -is the unit vector in the direction of axis ,  -is the thermal expansion coefficient, The buoyancy force and the external force 0 0 are taken into account in Euler Equation (1).Let us note the force tic amplitude of external force.We designate characteristic velocity, which is engendered by external force as -is Rayleigh number on the scale 0  .Further for the purpose of simplification we will consider the case Pr 1  .We pass to the new temperature T T R  , and obtain: Although we essentially pay attention to stable temperature stratification, unstable stratification can also be considered in the frame of this scheme.We use dimensionless writing of the equation more typically for the problem of convection.
We will consider as a small parameter of asymptotical development the Reynolds number  .The parameter will be considered neither big nor small, without any impact on development scheme (i.e.outside of the scheme parameters).

Ra
Let us examine the following formulation of the problem.We consider the external force as being small and of high frequency.This force leads to small scale fluctuations in velocity and temperature against a background of equilibrium.After averaging, these quickly oscillating fluctuations vanish.Nevertheless, due to small non-linear interactions in some orders of perturbation theory, non zero terms can occur after averaging.This means that they are not oscillatory, that is to say large scale.From a formal point of view these terms are secular, i.e. create conditions for the solvability of the large scale asymptotic development.So, to find and study the solvability equations i.e. the equations for large scale perturbations is actually the purpose of this work.Let us designate further the small scale variables as , and large scale ones as , and derivatives of large scale variables misunderstanding occurs between the temperature and the large scale time since time is argument and temperature is function).To construct a multi scale asymptotic development we follow the method which is proposed in work [16].We could start by establishing linear theory for instability development and after that pass to non linear theory, but as the non linear theory is technically less bulky, so we construct the non linear theory directly and then consider the linear limit.

The Multi-Scale Asymptotical Development
Let us search the solution for Equations ( 3) and ( 4) in the following form: First of all, we develop space and time derivatives in Equations ( 3) and ( 4) into asymptotical series of the form: .
Substituting these expressions into the initial Equations ( 3) and ( 4) and gathering together the terms of the same order, we obtain the equations of multi scale asymptotical development and write down the obtained equations up to order inclusive.Let us present the algebraical structure of the asymptotical development of the Equations ( 3) and (4) for the non linear theory (we will not write indices because they can be restored trivially at any moment).In the order there is only the equation:

) X
In the order 2 R  we have the equation : In the order 1 R  we get a system of equations: The system of Equations ( 12) and ( 13) gives secular terms: In zero order we have the following system of equations: Copyright © 2013 SciRes.
These equations give one secular equation: Consider the equations of the first approximation R: From this system of equations follow the secular equations:     The secular Equations ( 22)-( 24), are clearly obviously satisfied for velocity field geometry: In the second order , we obtain equations: It is easy to see that in the order there are no secular terms.

R
Let us come now to the most important order .In this order we obtain equations: From this we get the main secular equation: In these equations we do not write the law index   Besides there are secular equations: The Equations ( 33)-( 35) are satisfied in the previous geometry: , 0 , and, Const.
There is also an equation to find the pressure 3 P  : 3 .
These formulae show that when one knows the velocity it is possible to restore temperature and pressure.

Calculations of the Reynolds Stresses
It is clear that the essential equation for finding the non linear alpha-effect is Equation (31).In order to obtain these equations in the closed form we need to calculate the Reynolds stresses . First of all we have to calculate the fields of zero approximation 0 from the asymptotical development in zero order we have the equations: Let us introduce the operator : Using the operator , we write down Equations ( 38) and (39) in the form: Eliminating the temperature and pressure from Equation (41) we obtain: Here is the projection operator:

Dividing this equation by
, we can write it in the form: where is the operator: It is easy to make sure by a direct check that the inverse operator has the form: Consequently the expression for the velocity and temperature , takes the form: In order to use these formulae we have to specify in explicit form the helical external force 0 j F .The most simple and natural way is to specify the external force as deterministic.(Certainly, it is possible to specify the external force in a statistical way with specifying random field correlators, but this leads to more bulky calculations).As it is known helicity means that 0 0 rot 0  F F .Let us specify the force like so: where or It is evident that , where  -is the single pseudo scalar, i.e. helicity is equal to: The formulae (50) and (52) allow us to easily make intermediate calculations, but in the final formulae we obviously shall take 0 0 0 , , f k  as equal to one, since external force is dimensionless and depends only on dimensionless space and time arguments.The force (50) is physically simple and can be realized in laboratory experiments and in numerical simulation.
It is easy to write down the force (50) in complex form.It is evident that: where vectors A and B has the form: and 1 2 ,   D are given by formulae (52).The effect of the operator 0 on proper function has obviously the form: From this it is evident that: From the formulae (48) and (54), follows that the field is composed of four terms: where As was stated earlier, in scalar operators one can take Here we introduced the following notations: . Taking into consideration these formulae we can write down the velocities in the form: Copyright © 2013 SciRes.

OJFD
In order to calculate the Reynolds stresses we have first of all to calculate the expression: 0 0 01 01 03 03 2 Re .
Taking into account the formula (65 ), we obtain: Similarly taking into account formula (66), we obtain: It is clear that the components and are of interest.To begin with we consid o ponents of the tensor

T AA AA D A
The first bracket in the ( 71) is equal to zer why: o, which is Now consider the component : as far as .Consider now the component : The fir formula ( 74) is equal to zero, n: st bracket in the the Taking into account: The components take the form: Now, when we have these tensors components, w obtain closed equations for velocity.

Large-Scale Instability and Non Linear Vortex Structures
Let us write down in the explicit form the equations for non linear instability:


It is easy to see that with small values of the variables the Equations ( 82) and ( 83) are reduced to linear describe the linear stage of instability. 1  ons designates the single pseudo scalar because expressi , e Equations ( 84) and (85) differ from equation of ffect [16] by the coefficient are components of Th AKA-e rotW .
 only, but the im- portant difference between our equations and equations of AKA e ce of th y -eff ct is the presen e Ra leigh number in  coefficient.As a result: we obtain the instability increment 2 , Here new variables are introduced -are integration constants.The system of Equations ( 89) and (90) can be write down in the hamiltonian form: Here the variable plays the role of time hamiltonian and the H has t e form: h where function

 
U x has the form:   we are mostly interested in a boundless problem, it should be noted that thick closed lines correspond to the non linear structures which are also the solutions of the boundary problem with a rigid boundary: L where is the period over phase trajectory, which gets out th , from poi and gets back to the nt with The space st ed in      also solutions for the boundary problem with free boundaries.Thereby in the hamiltonian scheme which we consider there are three kinds of solutions: periodical waves, solitons and solutions moving to infinity.The last ones are not of interest from the point of view of the problem of large cale instability stabilizatio In conclusion it should be remembered that the system of the n. , W W 1 P and contributes to th ation for temperature (32).Closure of this eq made in much the same way as closure for locity.

n of Results
In this work, it is shown that in fluid with stable stratification a large scale instability appears under the action of small scale helical force.The result of the instability is the generation of vortex structures of the Beltrami type.These vortices have the characteristic vertical dimension e equuation is ve Nevertheless this equation is secondary and here we do not give the result of this closure.real space, it describes one spiral turn of the velocity vector field around the axe z.The soliton and kink are also the solutions for the boundary problem with free boundaries.

. Conclusions and Discussio 6
Let us return to the formulatio external helical force 0 f is given in the explicit fo order to make calculations more transparent.Str speaking, its explicit form is not very important for the existence of  -effect itself.It is only necessary that 0 0 rot f f  .The external force could be chosen statistically with specifying the correlator: It is fundamental that the last term   , G r  (helicity) in this co lator is not equal to zero, otherwise the rre  - effect is absent.Nevertheless the statistical method is more bulky since it requires us to specify the functions , , A B G and calculations of rather complicated integrals.
If we specify the external force dynamically, the averaging over fast oscillations is performed easily.In conclusion, it should be noted that te perature stratification is necessary for the existence of the instability.Previously it was supposed that this stratification was stable.However, the formulae for the large scale instability also admit the transition to an unstable fluid stratification, i.e. allow the substitution .Ra Ra   But one has to rethat the number Ra has to be sufficiently small so that the usual convective instability should not appear in the system.
Reynolds number on the scale 0 describe a spiral vortex structure (circularly polaed plane wave) with an am e which increases exponentially with time.These waves are sometimes called Beltrami runaways since for them there is no usual hydrodynamical linear instability vanishes but the non If the external force has a zero h linear remains.elicity, then the  -term vanishes in accordance with the Reynolds stress tensor [25].Helicity is taken into ure itself.If perature nishes general theorem of tem account in the external force struct the gradient va , then it is evident that the  -term also ationa s is o vanishes.It is clear that with the increasing 1 2, W W of the non linear terms decrease and the instability becomes saturated.As a result of the development and stabilization of instability, non linear vortex helical structures appear.The study of the form of these st structure f interest.For that purpose we take (82), these equations over z , obtain: z

Figure 1 .
Figure 1.Phase picture of the dynamical system with Ra = 2, C 1 = C 2 = 0.The bold line shows the phase trajectory which comes out of the point (1,1) and after the "time" Z = L comes back to the same point.This trajectory presents the stationary solution of the boundary problem with the rigid boundaries in the layer whose thickness is L = z.

Figure 2 .
Figure 2. Phase picture of the dynamical system with Ra = 3, C 1 = 0, C 2 = 0.The bold line shows the trajectory which corresponds to the stationary solution of the boundary problem with rigid boundaries with z = 0 and z = L.

Figure 5 . 9 . 10 .
Figure 5. Example of a vortex structure with Ra = 1.6, C 1 = 0, C 2 = 0. C 2 = 0.1 is shown in Figure 9.As above, the periodic vortex structures correspond to closed trajectories around elliptical points.Localized solutions (solitons) corre- to the separatrix on Figure 9. Since the separatrix connects two different hyperbolic points the soliton has n z   igure 6. F Phase picture of a dynamical system with Ra = 2, C 1 = 0, C 2 = 0.

Figure 7 .
Figure 7. Helical vortex structure with Ra = 2, C 1 = 0, C 2 = 0.This structure corresponds to the closed trajectory around the elliptical point in Figure 6.

Figure 9 .Figure 10 . 9 .
Figure 9. Phase picture of the dynamical system with Ra = 2, C 1 = 0.1, C 2 = 0.1.One can see the appea yperbolical and two elliptical points.rance of two h dimension much bigger than the vertical one.Since the vertical component of the z W e stratifi ha tion atio is equal to zero in the main approximation cation is stable, then the found instability ve any relation to convection.The structure of which describes the instability in linear n is the same as the equation of and th does not the equa approxim  -effector more precisely as the equation of AKA-effect.But as opposed to the AKA-effect, Rayleigh number.This means the inst ility vanishes in the non stratified fluid.As a resu stability generates plane spiral waves with circular rization (Beltrami runaway).With increases in amp de, the instability and its stabilization are described b n linear theory.With the instability has an essentially non linear r from the very begi ng.Stationary equations to be hamiltonian that is why they are a rich set of cal spiral vortex struc Notwithstanding that, in this work was essentially paid to the bounee problem, it should noted that so e perioout to be solutions for the boundary problem wit ing t correspond to the separatrix on the phase plane.T soliton of the new type.