The role of convective structures in the poloidal and toroidal rotation in tokamak
Abstract
The connection between the poloidal and the toroidal rotation of plasma in tokamak is important for the high confinement regimes, in particular in reactor regime. The sudden onset of closed convection structures in the poloidal section, due to the baroclinic production of vorticity, will sustain a fast increase of the poloidal velocity and a substantial effect on the toroidal rotation. However this is limited to the short time of the onset transition. In real plasma however there is random generation and suppression of convection cells and the sequence of these transient events can prove able to sustain the effect on the toroidal rotation. We formulate a simplified model which consists of a laminar sheared, regular, flow situated at the boundary of a region of driftwave turbulence. Vortical structures, randomly generated in this turbulent region are spontaneously advected toward the flow and are absorbed there, sustaining with their vortical content, the sheared flow. We examine this dynamics in the wavenumber space, using reasonable approximations. We derive a system of equations (which is a class of DevayStewartson system) and find that indeed, the vortices advected and absorbed into the layer can preserve its regular, poloidal, flow.
1 Introduction
Mixed regimes consisting of large scale flows (Hmode, Internal Transport Barriers) and turbulence are expected to be the current state in ITER. The rotation, either spontaneous or induced, will play a major role in the quality of the confinment. The efficiency of the sheared poloidal rotation to control the instabilities is considerably higher than that of the toroidal rotation but it is usually assumed that the poloidal rotation should be at the neoclassical level due to the damping induced by magnetic pumping. This is true if the drive of the poloidal rotation (besides the neoclassical one) relies on the Reynolds stress produced by a poloidally symmetric turbulence. Much higher drive of the poloidal rotation is however provided by flows associated with the convective structures that can be generated in the plasma crosssection beyond a threshold in the plasma pressure gradient. This drive overcomes the damping due to the magnetic pumping and the poloidal rotation is sustained. Cells of convection consisting of closed, large scale flows with unique sense of rotation can be spontaneously generated, triggered by streamers sustained by the baroclinic term able to generate vorticity. Similar to the ReyleighBenard first bifurcation (from purely conductive to convective regime), the onset is very fast and the drive exerted on the poloidal rotation leads to a fast time variation of the polarization radial electric field. This is sufficient to create a distinction between the phases (first half, second half) of bounce on a banana of trapped ions and, implicitely, leads to acceleration in the toroidal direction. This phenomenological model has been proposed recently [1] and the estimations are compatible with the observed effect of reversal of the toroidal rotation in tokamak. There are several aspects that still must be investigated and in particular the nature of suddenimpulse of the generation of the poloidal flow as envelope of the peripheral velocity of the convective structures. It is only on a very short time interval (onset of convection) that the time variation of the radial electric field can produce a toroidal effect on the bananas. If the convection cells become saturated and stationary then the efficiency of sustaining toroidal acceleration of bananas decreases rapidly to zero. The process must be repeated if it is to sustain its effect on the toroidal rotation.
On the other hand, it is highly improbable that the structure consisting of cells of convection breaking the azimuthal symmetry can exist as a stationary state. This is because the convection in the cells is much more efficient as radial transport of energy than any diffusive process sustained by small scale turbulence and the equilibrium profiles would be modified, suppressing the origin of convection. It is more realistic to regard the generation of convective structures as a stochastic process consisting of a random sequence of transient events, taking place on a range of spatial scales between scale drift vortices up to cells of convection with dimensions comparable with the minor radius . We should study the way these transitory structures create the impulsive increase of the poloidal velocity such as, on the average, to sustain the toroidal effect on bananas. A marginal stability regime can be established, where the external sources sustain the gradient of the pressure, the convection rolls are formed transitorily and the excess of the gradients relative to a threshold is transformed into pulses of heat and momentum in the radial direction. This problem requires a treatment based on the stochastic nature of the generation and destruction of the transversal convection cells, seen as a random transitory process. Such a treatment simply does not exist, not even for the famous RayleighBenard system.
This problem is important per se, beyond the phenomenological model for the reversal of the toroidal rotation. This process has been identified in experiments on a linear machine [2] and explained in the context of residual stress able to sustain the coherent flow. The sequence consisting of tilt, stretch and absorbtion of drift turbulent eddies was found to contribute significantly to the Reynolds stress. Coalescence of strongly tilted, almost poloidally oriented, eddies and formation of radially periodic and poloidally closed () flow was proposed in [3] as transition to rotation given by the stable solution of the nonlinear equation for the electrostatic potential at the edge (the FlierlPetviashvili equation). The absorbtion of a coherent convection cell inside a layer of sheared rotation is mentioned by Drake et al. [4]. It has also been discussed in similar terms in Refs.[5], [6].
The structure of the paper is as follows. We formulate a model which is relevant for the physical problem of sustainement of poloidal rotation by stochastic rise and decay of vortices on different scales. We estimate the average density of high amplitude vortices that are generated inside a volume of plasma (in our case it is question of an area within the poloidal crosssection since the third  toroidal  axis is considered irrelevant). Then we remind the known mechanisms on motion of a vortex on a background of gradient of vorticity, in order to see how random elements of vorticity of only one sign will join the flow layer where they are absorbed. Finally, we examine this dynamics in the wavenumber space and derive a model system of equations reflecting the interaction between the laminar flow and the incoming drops of vorticity. This is actually the Davey Stewartson system and indeed, the class of solutions called “solitoffs” show that the final state exhibits unidirectional (poloidal) alignment of the resulting flow. This confirms that the elements of vorticity coming from the turbulent region into the laminar flow do not perturb but are conforming themselves to the geometry of the flow, i.e. in the poloidal direction.
2 Formulation of the model
The effect of advection of impulselike convection structures (streamers, transient convection rolls of all sizes from down to the small scale of high amplitude robust drift vortices) to the layer of sheared laminar flow and their role in sustaining the laminar flow is a process that has aquired recognition for its importance in creating a coherent flow pattern. In the physics of the atmosphere and in fluid physics this is a well known phenomenon. It has been examined within the model of potential vorticity mixing [7] for the generation of an atmospheric jet, the equivalent of sheared laminar flow. Recent numerical simulation indicate the interesting and unexpected random switching of the flow between the two geometries: laminar flow and respectively a periodic chain of positive and negative vortices [8] which indicates that these two states are indeed neighbor in the functional space of the system’s configurations. The examination of the problem in real space is particularly difficult since it requires to analytically describe unequivocally the emergence of a coherent structure out of turbulence. One of the possiblity would be to derive a timeevolution equation for the effective number of degrees of freedom that are involved in the motion, as it results from a Karhunen  Loewe analysis. Such an approach is under examination but with the reserve that it may not really become a quantitative instrument.
Alternatively we can examine the system in the spectral space. The wave kinetic equation is the standard instrument to describe the dynamics of the energy density on elementary spectral intervals (seen as pseudoparticles), induced by nonlinear processes of wave turbulence [9]. We must adapt this approach to our particular problem, inspired by similar adaptations that have been made in several cases in the past. The most common models for dissipative equations of motion are the timedependent Landau  Ginzburg models. The model function is , the order parameter of the superfluid. In the absence of dissipation the general form of the equation is [10]
where is a real constant and the hamiltonian is
The total number of particles must be constant, as well as the energy const. The time dependent LG, is universal and occurs whenever there is a twophase system and a possible phase transition between them, with a function representing the order parameter. Of course, such a situation is subiacent to our more general and more complex problem. We take this as a suggestion to simplify our problem and reduce our objectives with the promise that they would become more accessible: we will try to prove that the advection of the transient convection structures to the laminar flow layer will not perturb the geometry of the flow, which remains mainly poloidally oriented. Then the conservation of the vorticity will guarantee that the content of vorticity of the individual structure joining the laminar sheared flow will be absorbed into the new, coalesced, flow and further will be redistributed to become homogeneous along the direction of flow. Or, the vorticity of a laminar flow is simply the derivative of the unidirectional velocity with respect to the transverse coordinate. The increase of the content of vorticity in the region of the laminar flow simply means that the total momentum (in case of a slab geometry) or angular momentum (if the flow is confined in a poloidal layer) is increasing. We then have a drive that may sustain the flow against the magnetic pumping. All this now requires to prove that absorbtion of the transient convection structures preserves the laminar flow and merge their vorticity content into the sheared flow.
In connection with our problem of stochastic stationarity of sustainment of poloidal flow by intermittent rise and decay of convection cells we construct a simple model to be investigated. We consider a layer of laminar sheared rotation where we assume that the instabilities are suppressed (this is an ideal representation of the zonal flow or of the mode layer). The layer is connex to a region where drifttype turbulence exists. The width of the layer is much smaller than the region of turbulence. In the driftturbulent region there is stochastic generation of convective motions and of strong vortical structures. The physical effect is different in the two limits: for large convective events there is an effect which we see as a jetlike impulse acting on the laminar layer; for small and robust vortices acceding the sheared flow, the effect consists of a forcing applied to the layers of the flow that have the same direction of flow as the rotation in the vortex while the layers where the two directions are opposite are decelerated. This enhances the local gradient of the velocity (i.e. the shear). Since the decelerated layers are in contact with the bulk, turbulent, plasma, the angular momentum is transferred by viscosity. While the thin laminar layer is accelerated, the large mass of turbulent plasma will only get a low rate of compensating rotation. This ensures the conservation of angular momentum. Although large and small vortices interact in different ways with the sheared flow layer, we will attempt a simplified treatment that retains the common aspects.
The shear rate (first derivative of the velocity of flow with respect to the coordinate transversal to the layer) is assumed nonuniform, which is equivalent to a gradient of vorticity (involving the second order derivative). The vorticity starts from zero at the boundary with the turbulent region and increases with the distance from it. Since the background has a gradient of vorticity a vortex will move in a direction which depends on the sign of its vorticity relative to the one of the background: the clumps (positive sign circulation vortices) are ascending the gradient of vorticity, toward the maximum of vorticity in the rotation layer and the holes (negative sign circulation vortices) are moving toward the minimum of vorticity, equivalently  they are repelled from the sheared layer [11]. In short, vortices of only one sign are joining the layer. The clump vortices will be absorbed by the sheared velocity layer with the transformation of their vorticity content into a local enhancement of the shear of the flow. This is a source of momentum which sustains the sheared flow, in particular sustains the sheared poloidal rotation against the damping due to the magnetic pumping. However it may appear as a spotaneous separation of ordered motion out of turbulence, and would pose a problem of entropy decrease. We have to remember that the vortices and the directed motions belonging to convection events are actually supported, via the baroclinic terms, by the equilibrium gradients (which are externally sustained). Similar physical processes are probably behind the formation of the Internal Transport Barriers and the formation of the mode rotation layer. A last comment is in order, in connection with what appears to be a spontaneous separation of the vorticities of the opposite signs. This is indeed what happens in the ideal fluid at relaxation as result of the inverse cascade, and is shown by experiments and numerical simulation [12] and also derived analytically [13]. Similar separation of the two signs of vorticity has been seen for plasma [14] and atmosphere, also supported by analytical derivation [15]. The result is a final state consisting of a pair of vortices of opposite signs that have collected all positive and respectively all negative vorticities in the field (a dipole). We now note that the dipolar state can be realized as either two distinct vortices of opposite signs, or, in a circular geometry, as a central part with vorticity of one sign surrounded by an annular region of vorticity of the opposite sign. Transitions are possible between the two ”dipolar” geometries of vorticity separation and the change is very fast. We only mention this facts in support of our phenomenological representation of the sustainment of the sheared rotation, but we will not develop this subject here.
We take the layer of plasma flow in the ( longitudinal, streamwise, poloidal in tokamak) direction with sheared velocity varying in the ( transversal on the layer of flow, radial in tokamak) direction. It is described by a scalar function which is the complex amplitude of slow spatial variation of the envelope of the eddies of drift turbulence. From the components of the velocity are derived as where is the versor of the magnetic field. We also define a scalar function which represents the field associated to the incoming vortices that will interact and will be absorbed into the sheared, laminar, flow, transfering to it their vorticity content. The process in which a localized vortex is peeledoff and absorbed into the sheared flow may be seen as a transfer of energy in the spectrum from the small spatial scales involved in the structure of the vortex toward larger spatial scales of the sheared flow, with much weaker variations on the poloidal direction. This can be seen as a propagation in space.
We must first estimate what is the number of high amplitude, robust vortices that are generated by the driftturbulence in a unit area. This problem is treated in the Appendix 1 where it is shown that the density of strong vortices is where is the typical size of the drift eddies normalized to .
3 How these vortices ascend to the flow layer?
The answer is: due to the existence of the background gradient of vorticity. This has been clearly described for nonneutral plasma in [11] and is also known in the physics of fluids [17], the physics of the atmosphere [16] and in astrophysics (accretion disks) [18].
Since the background has a gradient of vorticity a vortex will move in a direction which depends on the sign of its vorticity relative to the one of the background:

the prograde vortices are moving toward the maximum of vorticity and

the retrograde vortices are moving toward the minimum of vorticity.
The prograde vortices eventually are absorbed by the sheared layer and they contribute with their vorticity content to the momentum of the background flow. This is a source of momentum which sustains the sheared flow, in particular sustains the sheared poloidal rotation against the damping due to the magnetic pumping.
A clump of vorticity (positive circulation ) goes up the gradient of background vorticity. For a sheared velocity flow with poloidal velocity increasing linearly with radius, with a coefficient
but in a layer which has background vorticity along , the radial motion of the positive vortex (position ) is
where [19]. Schecter and Dubin [11] find . The position of the vortex is greater but comparable to the typical length , . The function is bounded by .
The negativecirculation vortices, called holes are moving in the opposite direction therefore they remain in the turbulent region and are repelled by the layer of sheared flow with global conservation of angular momentum. This means that only vortices of a certain sign ascend to the flow and are absorbed in it. In the following we will analyse this process as an interaction between the two fields: and the field representing the incoming vortices.
4 Dynamics in the  space of the interaction between the laminar flow field and the field of the vortical cells
In order to derive a dynamic equation for the scalar streamfunction of the sheared flow we remind the properties resulting from multiple spacetime analysis of a turbulent field, like the drift waves interacting with a coherent, sheared flow. The typical equation in the hydrodynamic description is RayleighKuo or barotropic equation and the multiple spacetime scale analysis leads to the Nonlinear Schrodinger Equation for the amplitude of the envelope of the oscillatory waves excited in the fluid [20]. The general structure of this equation is
(1) 
The interesting term is the nonlinearity (the last term) which may be seen as arising from a selfinteraction potential (see the next Section). The later potential results naturally from the multiple scale analysis and it is the effect on the envelope amplitude of the nonlinearity consisting of convection of the vorticity by its own velocity field. We must work with an amplitude in the spectral space, which we see as the envelope of oscillations representing propagation of perturbation in the spectrum, this propagation being a reflection of what is the process of the incoming vortex dissolving itself and being absorbed by the flow,  in the physical space. Looking for an analytical model for in the space we will use the structure which is suggested above.
We expect to have a diffusion in space and the same selflimitting nonlinear term
(2) 
and we have added a source that represents the interaction of the basic flow with the vortices that are randomly generated in the drift turbulence and are joining the flow.
The addition of “drops” of vorticity is represented by the space scalar function
(3) 
Several space scales are present in the structure of the function in real space. The reason is that we must consider that the vortices that are continuously generated in the turbulent region have two characteristics:

the space extension of individual vortices coming into the layer has a wide range of values. The range extends from vortices of the dimension (the natural result of excitations sustained by the ionpolarization drift nonlinearity) up to convective events representing streamers with closed (rolltype) geometry and spatial extension given by the inverse of for of only few units.

the spatial distribution of the place where these vortical structures are reaching the layer of sheared flow is random over practically all the circumference of the poloidal flow.
We will draw conclusions about the space profile of based on the contributions coming in from various scales in the real space. The most simple description neglects the random aspect of the content of and simply ennumerates the possible spatial scales. For this we adopt a discrete representation of the dimensions of vortices, starting from a minimal size of the order of up to a low fraction of the circumference of the poloidal rotation layer, this last scale corresponding to large convection rolltype events.
The above elementary assumption, i.e. the discrete number of space extension of vortices, implies that the function presents higher amplitudes around wavenumbers which correspond to the inverses of the spatial scales. The simplifying assumption which takes uniform spatial distribution of the events of arrival of vortices (of any scale) to the layer implies periodicities on various spatial scales.
We consider separately the direction (poloidal). One periodicity can be at the level of a smallest vortex, . This means that a sequence of maxima of appear with spatial periodicity and this corresponds to the addition of smallest vortices into the flow layer. Further, on a longer space scale there is another periodicity, . This means that the sequence on the lower scale (with much denser spatial granulations of vorticity, ) is modulated on a longer space scale with a periodicity . This is because larger vortices (but not yet convection rolls) are absorbed into the flow layer. Large vortices carry smaller vortices. These considerations can be extended to other, intermediate, levels of periodicity. Finally we can have the largest (accessible) spatial scale where a convective event occurs as a result of a stochastic event, or a streamer sustained by the baroclinic term and having the closeup property.
We note that the periodicity in the poloidal direction cannot occur without a corresponding periodicity in the radial direction, as far as the vortices and convective events are concerned. It is reasonable to assume a linear dependence : small scale periodicities in the poloidal direction (small scale robust drift vortices) occur on a similar radial scale, i.e. again small periodicity in the radial direction. With larger wavelengths, we expect to involve also larger radial wavelengths since the physical process involves more fluid motion in the closedup convection cell.
This suggests a proportionality between the periodicities on the poloidal direction, , and the periodicities on the radial, , direction, in the structure of . This means that the Fourier components along (e.g. at ) are proportional with the Fourier components along (at ) with the same proportionality coefficient on all the spectral interval. The fact that there is (assumed) proportionality on the periodicities on and on , requires for analytical description a hyperbolic operator
(4) 
where the D ’Alamebrtian operator is the signature of the correlated periodicities on the poloidal and radial directions of the field representing the incoming drift vortices  up to convection rolls, generated in the drift turbulence region, moving against the gradient of vorticity, and finally absorbed into the rotation layer. The ”velocity ” is the ratio of the periods on the two space directions and we take it . The centers of the vortices reach points that belong to the lines . The general form of the solution of the homogeneous equation is
(5) 
When is ”free” it is double periodic and it verifies the above equation. However the destruction (peelingoff) of the incoming vortices and their absorbtion through interaction with the background flow, , is the loss of the double periodicity. The departure relative to the space double periodicity (the ”free” state of ) feeds the flow field such as to enhance it in spectral regions comparable to those where the vortices are localised, i.e. higher . Basically, if the free is localised on close to (equivalent to an uniform poloidal flow), the source induced by is a charge which acts through a Poisson equation on the ”potential” given by the amplitude . We should recall that we want to represent the physical process in which incoming vortices transfer via absorbtion their vortical content to the sheared flow, sustaining and/or enhancing the poloidal, , rotation. This is mainly manifested as the evolution of an initial oscillation generated by on the spectrum with collection of the energy in a region close to , since this corresponds in real space to the uniformization of the flow in the poloidal direction, after an event of absorbtion of a vortex. Then, ignoring the less significant periodicity in direction, corresponding to periodic perturbation propagating in the (radial) direction, the interaction may be represented as
(6) 
The constant is a measure of the permitivity of the equivalent electrostatic problem. Using the specific terminology measures the decrease of the effectiveness of the effect of the charge (right hand side of the Eq.(6)) on the potential due to the “polarization”. Essentially, the amount of destruction of the incoming vortex, measured by the departure from the “free” state by the right hand side of the Eq.(6) is transferred to the background flow by first exciting elementary waves in real space, while is only an envelope. The perturbation of the sheared flow (via the left hand side of the Eq.(6)) results from the nonlinear interaction of these excited waves, i.e. by convecting the modified vorticity with the background velocity flow. This is the nonlinear term of the righthand side of Eq.(2). Then it is reasonable to assume .
Returning to the equation for we specify the interaction between the two fields and in the simplest way
(7) 
Then
(8) 
to which we add the equation for ,
(9) 
These two equations are known as DaveyStewartson system and our case is DSI [21]. It is exactly integrable and several analytical solutions are available, in terms of Riemann theta functions or Jacobi elliptic functions. Of particular importance for our physical problem is the long wave limit solutions, as obtained by Chow using the Hirota method. We reproduce here his result [21], with details given in Appendix B.
(10) 
These solutions are called solitoffs and they decay in all directions except a preferred one. They are ”semiinfinite solitary waves”. For example, for and the solution shows a peak of around which in physical terms means quasiuniform flow along the poloidal direction. This is compatible with what we expect from this analytical model: the turbulence advects “drops” of vorticity into the rotation layer and this does not destroy the flow but sustains it. It is actually a discrete Reynolds stress, or a process which is reversed version of the usual KelvinHelmholtz instability.
5 Discussion
The Figure (1) represents the amplitude of the solution given in Eq.(10). We now want to understand in more detail what was the origin of the final spectral structure (i.e. the solution written above) in our previous analysis, which was mostly qualitative and based on several assumptions.
We note that the solution looks compatible with our idea about the buildup and/or sustainment of a sheared laminar flow out of random convection of momentum from the deep region where turbulence is active. We first refer to the term . Essentially this term is the quantitative translation of the physical fact that the flow, after a perturbation (which later we will specify as ”absorbing” an incoming vortex), will try to reestablish the uniformity of the flow by spreading the new perturbing vorticity as shear of the flow. Although this is a very approximative image we expect that this term in the action functional (i.e. in the integral of the Lagrangian density) to penalize departures of the flow configurations from uniformity of the flow. The departures are measured as the square of the difference between the actual magnitude of and a reference uniform configuration. Then the system moves in the function space in a potential with two extrema (symmetric relativ to ), so that in its Lagrangian there is a term of nonlinear selfinteraction of the scalar field , of the general type
with . This means that the stationary states will try to attain one of the equilbria given by the extremum of . In our simplified case the ”reference configuration” that the system takes as ground state above which it penalizes any departure is simply . This means that we represent the simultaneous wiping out of the local perturbation and the suppression of the shear, without having any reason for the latter, except for tractability. The potential makes the uniformity ( in this simplified description) an attractive state, if no source term is present. With the nonlinearity becomes the same as in the Nonlinear Schrodinger Equation, as we have mentioned above, in the construction of the model.
On the other hand we see that in the presence of a source the solution for fixed, greater than zero, , is not exactly at the ”poloidal uniform flow” but shifted, the spectrum is peaked in at a small but finite value. This means that our description actually obtains longwave poloidal oscillations but not exactly symmetric poloidal flow. This shift is due to the average of the spatial extensions of the incoming perturbations and is dependent on the scale factor (’wavelength’) . The two parameters and must be seen as representative space scales along the poloidal respectively radial directions.
We have simply neglected the structure of the spectrum on (radial). In the solution it occurs with a uniform space contributions for positive wavenumbers (the function rises rapidly after a space scale and remains constant after that for all ). This means that the structure of the flow on the transversal direction (a cross section along the radius) contains all possible motions and modulations, without however placing emphasis on any of them. It becomes effective the assumption that the convections and vortices are coming in all sizes, and we have neglected the second derivative of the flow amplitude to in Eq.(6), compatible with the approximate independence of on .
We note that the identical reproduction of the localized maxium close to at periodic intervals along axis is the result of the periodicity
of the Jacobi elliptic function , where
For slightly elongated perturbations along the direction of the flow (poloidal)
we have
which means that the periodicity does not appear to be relevant in our case.
6 Conclusion
In conclusion, this simple analytical analysis supports the idea that the randomly generated vortical motions, including the convection rolls that rise and decay transiently in the plane transversal to the magnetic field, are able to sustain poloidal rotation. When the time scale is very short, as for convection rolls driven by the baroclinic effect, the poloidal flow arising as envelope of the rolls has fast time variation and the neoclassical polarization effect acts effectively on the toroidal rotation. This process accompanies every event consisting of configuration of a closed convection pattern. Then, even if these events are random, their statistical average is an effective source of poloidal flow.
Acknowledgement 1
Work supported partially by the Contracts BS2 and BS14 of the Association EURATOM  MEdC Romania. The views presented here do not necessarly represents those of the European Commission.
Appendix A Appendix. Estimation of the density of high amplitude vortices generated in the driftwave turbulence
In this Appendix we discuss the problem of estimating the average number of strong vortices that are generated in a twodimensional region of driftturbulent plasma. Being a complex problem, it will be the subject of a separate work. We limit ourselves here to mention the main steps of this calculation, revealing the necessity to link apparently different systems.
We adopt the definition that a strong, robust vortex of high amplitude is one which is a physical realisation of the mathematical singular vortex, the latter being characterised by . This approximation is useful since we will now look for the statistically averaged density of singular vortices in a region of turbulence.
The following connections can be made:

The drift wave turbulence in the hydrodynamic regime and at spatial scales that are much larger than is governed by the cuasithreedimensional HasegawaMima equation. However, for strong magnetic field, there is not too much difference to ignore the intrinsic length, the Larmor radius, and to consider the approximate description given by the scalefree hydrodynamic model, i.e. the Euler equation

It is known that this equation is equivalent to a model of pointlike vortices interacting in plane by a selfgenerated longrange (Coulombian) potential as known from prestigious theoretical work Kirchhoff, Onsager, Montgomery, etc. (see [22]).

The dynamics of the pointlike vortices has been mapped onto a field theoretical model consisting of a complex scalar (matter) field , , and a gauge field (which mediates the interaction) [13]. The Lagrangian density consists of the kinetic term for , the ChernSimons term for and a nonlinear interaction. The field variables are in the algebra due to the vortical nature of the physical elements
with the ladder generators in the cartan algebra. Introducing , , the vorticity is
and in the relaxed states (selfdual) we have
here is the streamfunction of the flow, with , .

We have noted that this equation is exactly the equation describing the surfaces in real space with the property of constant mean curvature const (CMC). If the curvatures of the surface are and the CMC surfaces have

The fact that the selfduality states and the CMC surfaces are governed by the same equation sinhPoisson has a deep significance. Extending the parallel between these two system at states which are not at selfduality or at CMC, (i.e. the dynamical states preceeding the asymptotic relaxed states) we have established the following mapping

It then results that the umbilic points of a surface in real space, i.e. a point where
is a point where
and this means that
at these points.

Now we see that looking for the points where there is a very high value of the vorticity (almost singular vortices ) means to look for the umbilic points () of a surface.

We actually consider a turbulent region where the strong, cuasisingular vortices appear at random. Equivalently, we must consider a surface in which has fluctuating shape. The statistical properties of the random strong vortical structures arising in turbulence is the same as the statistical properties of the set of umbilic points on the fluctuating surface. On the other hand the statistics of the fluctuations of the surface is assumed Gaussian.

We then have to find the statistical properties of the set of umbilic points on a surface whose shape fluctuates with the Gaussian statistics.
This problem is solved in Ref. [23]. For a Gaussian correlation of the surface heigth , the density of umbilic points is
and this is also the density of high amplitude robust drifttype, vortices generated by turbulence. Here is the linear dimension of a typical eddy of the drift wave turbulence, normalized to .
Appendix B Appendix. Details about parameters of the exact solutions
The solution written in Eqs.(10) and (4) have been obtained by Chow [21] using the Hirota method. The parameters are defined with the expression of written in terms of the Riemann functions , the general solution. The connection with the original expression in terms of two ”spatial” coordinates consists of replacing, for our particular case
Then, one has
The variables are pure imaginary parameters of the nome which defines the expression of the theta functions as series in powers of . The relations between the space scaling factors of the and variables and the two nome variables is
The coefficient is
It is possible to replace the functions with Jacobi elliptic functions and so new notations are introduced:
and are the complete elliptic integrals of modulus . We see that the new (stretched) spatial variables are and . Since in this solution and are regarded as spatial coordinates, their coefficients and are called ”wavenumbers”. In our use of this solution the situation is reversed, and which means that are spatial wavlengths along these directions.
The restriction exists: and the long wavelength limit is
leading to
For definitions and properties of the Riemann functions see Abramovitz and Stegun [24] or NIST Digital Library of Mathematical Functions [25].
References
 [1] F. Spineanu and M. Vlad, Nucl. Fusion 52 (2012) 114019.
 [2] M. Xu, R.G. Tynan, P.H. Diamond et al. Phys. Rev.Lett. 107 (2011) 055003.
 [3] F. Spineanu, M. Vlad, K. Itoh, S.I Itoh, http://arxiv.org/pdf/physics/0311138.
 [4] J.F. Drake, J.M. Finn, P. Guzdar et al. Phys. Fluids B4 (1992) 488.
 [5] R.Z. Sagdeev, V.D. Shapiro, V.I. Shevchenko, Sov.J.Plasma.Phys. 4(3) (1978) 306.
 [6] G.I. Soloviev, V.D. Shapiro, R.C.V. Somerville, B. Shkoller, J. Atmos. Sci. 53 (1996) 2671.
 [7] R.K. Scott and D.G. Dritschell, J. Fluid Mech. 711 (2012) 576.
 [8] F. Bouchet and E. Simonnet, Phys. Rev. Lett. 102 (2009) 094501.
 [9] V.E. Zakharov, V.S. L’vov, G. Falkovich, Kolmogorov Spectra of Turbulence I, Spinger Verlag, Berlin, 1997.
 [10] P.C. Hohenberg, B.I. Halperin, Rev. Mod. Phys. 49 91977) 435.
 [11] D.A. Schecter and D.H.E. Dubin, Phys.Fluids 13 (2001) 1704.
 [12] D. Montgomery, W.H. Mathaeus, W.T. Stribling, D. Martinez and S. Oughton, Phys. Fluids A4, (1992) 3.
 [13] F. Spineanu, M. Vlad, Phys. Rev.E 67, (2003) 046309.
 [14] R. Kinney, J. C. McWilliams and T. Tajima, Phys. Plasmas 2 (1995) 3623.
 [15] F. Spineanu, M. Vlad, Phys. Rev. Lett. 94 (2005) 025001 and http://arxiv.org/pdf/physics/0501020.
 [16] B. Wang, X. Li and L. Wu, Journal of Atmospheric Sciences 54 (1997) 1462.
 [17] P.S. Marcus, J. Fluid Mech. 215 (1990) 393.
 [18] H. Lin, J.A. Barranco, P.S. Marcus, Center for Turbulence Research, Annual Research Briefs 2002, Stanford University.
 [19] D.A. Schecter, D.H.E. Dubin, Phys. Rev. Lett. 83 (1999) 2191.
 [20] F. Spineanu and M. Vlad, Phys. Rev. Lett. 89 (2002) 185001, and: http://arxiv.org/pdf/physics/0204050.
 [21] K.W. Chow, Wave Motions 35 (2002) 71.
 [22] R.H. Kraichnan and D.C. Montgomery, Rep.Prog.Phys.43 (1980) 547.
 [23] M.V. Berry and J.H. Hannay, J. Phys.A: Math. Gen., 10 (1977) 1809.
 [24] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, National Bureau of Standards Applied Mathematics Series 55, (Tenth printing) 1972, Washington D.C.
 [25] NIST Digital Library of Mathematical functions, http://dlmf.nist.gov.