## Lack of self-similarity in transverse velocity increments and circulation statistics in two-dimensional turbulence

Nicol´s a P. M¨ uller 1, 2 and Giorgio Krstulovic 2

1

Laboratoire de Physique de l' ´ cole E normale sup´ erieure, ENS, Universit´ PSL, e CNRS, Sorbonne Universit´ e, Universit´ Paris e Cit´, e F-75005 Paris, France 2 Universit´ Cˆte d'Azur, Observatoire de la Cˆte d'Azur, CNRS, Laboratoire Lagrange, e o o Boulevard de l'Observatoire CS 34229 - F 06304 NICE Cedex 4, France

(Dated: November 27, 2024)

We study the statistics of longitudinal and transverse structure functions, as well as velocity circulation, in the inverse energy cascade of two-dimensional turbulence. Using direct numerical simulations of the incompressible Navier-Stokes equations, we show that transverse structure functions exhibit an anomalous scaling, in contrast to the self-similar behavior of longitudinal ones. We derive an analytical relation that shows that the scaling exponents of transverse structure functions and velocity circulation are related in two-dimensional turbulence.

Two-dimensional (2D) turbulence is a paradigmatic system in fluid dynamics as it strongly differs from threedimensional (3D) turbulence. One of the most remarkable phenomenon which takes place in 2D turbulent flows is the inverse energy cascade, a process in which energy is transferred towards large scales leading to the formation of large-scale coherent structures [1-3]. This mechanism is relevant in geophysical flows such as atmospheres and oceans, which exhibit quasi-2D properties due to the suppression of motion in one direction induced by confinement, stratification, rotation, or other effects [4-6]. A prime example of this phenomenon is Jupiter's giant red spot [7].

Unlike the direct energy cascade in three dimensions [8], it has been widely accepted that the 2D inverse energy cascade is not intermittent. More precisely, most of the attention has been devoted to the longitudinal velocity increments δu ∥ r = ( u x ( + r ) -u x ( )) · ˆ, with r u the velocity field, x the space coordinate and r the increment. This quantity has been shown to be self-similar for a given length scale r within the inertial range [9, 10]. It follows that the longitudinal structure functions (LSF) of order p

$$S _ { p } ^ { \| } ( r ) = \langle | \delta u _ { r } ^ { \| } | ^ { p } \rangle \sim r ^ { C _ { p } ^ { \| } }$$

exhibit scaling exponents ζ ∥ p that follow a self-similar dependence p/ 3, consistent with Kraichnan-Leith-Batchelor (KLB) theory [1-3]. The angular brackets ⟨ ⟩ . indicate averaging over space and time. Note that this is in stark contrast to 3D homogeneous and isotropic turbulence, where strong fluctuations in the velocity field lead to deviations from the Kolmogorov self-similar theory (K41) predicting ζ ∥ p = p/ 3. Deviations from this scaling are typically described using multifractal theories [8].

In recent years, the study of the intermittent behavior in turbulent flows has taken a new direction through the analysis of velocity circulation, a conserved quantity in ideal fluids defined as the line integral of the velocity field around a closed loop C r of linear length scale r

$$\Gamma _ { r } = \oint _ { \mathcal { C } _ { r } } u \cdot \text{d} l.$$

High-resolution direct numerical simulations (DNS) of 3D turbulence showed an intermittent behavior in the statistics of velocity circulation, with scaling exponents λ p of a different nature than velocity increments [11]. That is, the naive relation between the scaling exponents of structure functions and circulation moments, expressed as λ p = ζ p + , is p not satisfied. This behavior was also observed in quantum turbulence [12-14], where the nature of vortices drastically differs from that of classical fluids, motivating the development of new intermittency theories based on circulation [15-17]. A recent experimental study in quasi-2D turbulence showed that the statistics of velocity circulation is intermittent [18], and was later confirmed in numerical simulations of 2D classical and quantum turbulence [19]. This striking observation is in strong contrast with the established self-similar picture of 2D turbulent flows. An explanation for the surprising difference between fluctuations of circulation and velocity increments remains elusive.

A priori, the relation between increments and circulation can be understood when writing explicitly the velocity circulation around a square loop of opposite vertices ( x , y 0 0 ) and ( x , y 1 1 ) = ( x 0 + r, y 0 + ) r

$$\Gamma _ { r } \ = \ \int _ { x _ { 0 } } ^ { x _ { 0 } + r } \left [ u _ { x } ( x, y _ { 0 } + r ) - u _ { x } ( x, y _ { 0 } ) \right ] d x + \int _ { y _ { 0 } } ^ { y _ { 0 } + r } \left [ u _ { y } ( x _ { 0 } + r, y ) - u _ { y } ( x _ { 0 }, y ) \right ] d y.$$

Wenote that the integrands correspond to the transverse velocity increments in different directions [20]. The transverse scaling exponents ζ ⊥ p are defined from the transverse structure functions (TSFs) as

$$S _ { p } ^ { \perp } ( r ) = \langle | \delta u _ { r } ^ { \perp } | ^ { p } \rangle \sim r ^ { \zeta _ { p } ^ { \perp } },$$

with δ u ⊥ r = δ u r -δu r ∥ r ˆ the transverse velocity increments. It follows that a refined dimensional guess for the circulation scaling exponent would be λ p = ζ ⊥ p + , but this relation does not hold in 3D turbulence. p

Scaling exponents of LSFs and TSFs in 3D homogeneous and isotropic turbulence are known to differ, the latter being more intermittent [21]. Furthermore, high-resolution DNS showed that the transverse scaling exponents ζ ⊥ p saturate for high-order moments [22]. In 2D, the TSFs have received much less attention.

In this Letter, we study the statistics of velocity increments and circulation in the inverse energy cascade of 2D turbulence by means of DNS. We look independently at the longitudinal, transverse and circulation structure functions and do a general comparison of their scaling properties. We also derive an analytical relation between the circulation and the transverse scaling exponents.

We describe the dynamics of a 2D flow using the incompressible Navier-Stokes equations for the scalar vorticity field ω = ( ∇ × u ) · ˆ z

$$\partial _ { t } \omega + \{ \omega, \psi \} = \nu \nabla ^ { 2 } \omega - \alpha \omega + f,$$

where the vorticity and the stream function satisfy the relation ω = -∇ 2 ψ , with ν the kinematic viscosity and f an external forcing. The Poisson brackets are defined as { ω, ψ } = ∂ ω∂ ψ x y -∂ ω∂ ψ y x , and the term -αω is a linear frictional damping that dissipates energy at large scales. We analyze data from two different DNS, one kindly provided by Guido Boffetta (RUN-A) with N x = N y = N = 2048 linear collocation points, and a second one at a larger resolution performed by our group (RUN-B) with N = 6144 [19]. In both cases, Eq. (5) is solved in a 2D periodic square of side 2 π using a pseudo-spectral method and evolved using a second-order Adams-Bashforth (RUN-A) or Runge-Kutta (RUN-B) scheme. For RUN-A, the forcing is a random white noise in time with a Gaussian correlation in space, which satisfies ⟨ f ( x , t ) f (0 , t ′ ) ⟩ = F l 0 2 f δ t ( -t ′ ) exp[ -x / 2 (2 ℓ 2 f )] with ℓ f the forcing scale [10]. For RUN-B, we use a constant in time Gaussian forcing in Fourier space ˆ ( f k ) = exp[ - | ( k | -k f ) 2 / (2∆ k 2 )], with k f = ℓ -1 f the forcing wave number and ∆ k its width. We choose the parameters α and k f in order to maximize the width of the inertial range. Due to the implementation of different forcings in each simulation, the statistical properties of the 2D turbulent flow slightly differ at the forcing length scale, but become negligible in the inertial range. The turbulent properties of the flow are averaged in time once the system reaches a statistically steady state. We average 28 fields for RUN-A and 453 fields for RUN-B. The parameters of each DNS are shown in Table I.

The energy spectra for both runs develop a clear k -5 / 3 scaling law between the large-scale dissipation wave number k α ∼ ℓ -1 α with ℓ α ≃ ϵ 1 / 2 α -3 / 2 , and the forcing wave number k f (see Fig. 1). In this same range of scales, the energy flux

$$\Pi ( k ) = \langle u ^ { < k } \cdot [ u \cdot \nabla u ] \rangle$$

takes negative values, signature of the energy cascading towards large scales. The superscript &lt;k indicates that the velocity field considers only wave numbers such that k &lt; | k | and ϵ = k -1 f ∑ k f k =1 Π( ) is the mean energy flux. k The formation of a condensate at small wave numbers is avoided by the damping term. Note that the energy flux is not exactly flat in the inertial range due to the use of a linear damping. It has been observed that employing an hypoviscosity term flattens the energy flux in the inertial range [24]. It is expected that in the limit of an infinite scale separation ℓ α /ℓ f →∞ , the flux becomes constant in a larger range of scales. Figure 1 (c) shows the energy spectra compensated by the KLB prediction Cϵ 2 / 3 k -5 / 3 with C = 6 [23].

For the analysis of structure functions, we compute the velocity increments along lines in both the x and y directions. The top row of Fig. 2 shows the LSF defined as in Eq. (1) for even-order moments up to p = 10 for both DNS. The

TABLE I. Parameters of the 2D numerical simulations. N is the number of linear collocation points, α the linear friction coefficient, ν the viscosity, L/L f the ratio between the system size and the forcing scale, ϵ the mean energy flux of the inverse cascade, ℓ I = ∫ k -1 E k ( )d k/ ∫ E k ( )d k the integral scale, and R I = ( ℓ /L I f ) 2 / 3 a Reynolds number based on the integral and forcing scales.

| RUNS   |    N | α        | ν          |   L/L f |     ϵ |   ℓ I /L |   R I |
|--------|------|----------|------------|---------|-------|----------|-------|
| RUN-A  | 2048 | 0 . 02   | 10 - 5     |     256 | 0.001 |    0.135 | 10.61 |
| RUN-B  | 6144 | 0 . 0328 | 8 × 10 - 6 |     768 | 0.033 |    0.272 | 35.2  |

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-ca7c1564-d9b5-4c08-a366-c8404021bac9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-c99fe15e-2a25-46c0-8652-76b4ebe5f27f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-9941cc50-e9ec-4a40-bbdc-1ab3c520e7c3.jpg)

FIG. 1. (a) Energy spectra, (b) energy flux, and (c) compensated energy spectra for both DNS. The wave number is normalized by the forcing wave number k f , and the flux is normalized by its mean value ϵ . For the compensated spectra we used the constant C = 6 [23].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-a628a8bd-3415-4e51-9b0b-df8aaa6fea34.jpg)

f

FIG. 2. Longitudinal structure functions in numerical simulations of 2D turbulence in the inverse energy cascade for moments up to p = 10 for (left column) RUN-A and (right column) RUN-B. Top panels show the structure functions as a function of r , and the bottom panels show the ESS method. The shaded area indicates the inertial range. The insets show the local slopes ζ ∥ p = dlog S / ∥ p dlog r . Dashed horizontal lines in the insets show the self-similar scaling ζ ss p = p/ 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-2b486b14-98dc-4d33-9bf6-173da7ce5b59.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-96bb5c2c-74e5-4351-88e9-a10bd877feac.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-2b51dde4-94a8-4ee3-a240-dae12265580f.jpg)

insets display the local slopes, defined as the logarithmic derivative ζ ∥ p ( r ) = d log S / ∥ p dlog r . Note that for large scales r &gt; ℓ α , the flow is decorrelated and the structure functions become flat, while for small scales r &lt; ℓ f they follow the viscous scaling S p ∼ r 2 p . In the inertial range, the LSFs present logarithmic corrections which affect the scaling properties of the flow [9]. The deviations from a pure power-law scaling are stronger for RUN-A than for RUN-B, probably due to the different resolutions used in each dataset, as well as the peculiarities of the forcing. To reduce these systematic deviations, we employ the extended self-similarity (ESS) method with respect to the third-order moment S ∥ 3 ( r ) = ⟨| δu ∥ r | 3 ⟩ [25]. The shaded area indicates the inertial range ℓ f &lt; r &lt; ℓ α defined from the ESS method. Note that, due to the implementation of different forcings, the moments at the injection scale clearly differ from each other. Nevertheless, this effect in the inertial range is suppressed employing the ESS method, as shown in the bottom row of Fig. 2. The new local slopes defined as d log S / p dlog S 3 display clear scaling properties that follow the

well-known self-similar behavior ζ ∥ p = p/ 3 [10].

FIG. 3. Transverse structure functions up to p = 10 in numerical simulations of 2D turbulence in the inverse energy cascade for (left column) RUN-A and (right column) RUN-B. Top panels show the structure functions as a function of r , and the bottom panels show the ESS method. The shaded area indicates the inertial range. The insets show the local slopes ζ ⊥ p ( r ) = dlog S ⊥ p / dlog r . Dashed horizontal lines in the insets show the self-similar scaling ζ p = p/ 3.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-573f9546-8fae-40d6-a5df-b6926ccce865.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-80af6d87-114f-4330-84ca-26f9b397ebd7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-78cb3115-ecc3-403e-836a-08e6740f61eb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-01148fc4-29cc-48eb-aa1f-c34c741e07e6.jpg)

Following the same procedure, we compute the TSFs. Figure 3 shows the TSFs for both datasets. Again, the structure functions display logarithmic corrections in the inertial range that we reduce by employing the ESS method with respect to S ⊥ 3 . Strikingly, the local slopes deviate from the self-similar prediction, suggesting that fluctuations of transverse increments are stronger than longitudinal ones. This behavior is somehow similar to the case of 3D turbulence, where the ζ ⊥ p saturate for p ≥ 10 [22]. We do not observe this saturation in 2D turbulence.

We now focus on the circulation statistics to see if there is a relation with the intermittency of TSFs. We compute the velocity circulation around squared planar loops of size r over the whole space. The calculations are performed in Fourier space to preserve the spectral precision from the simulations [12]. The circulation moments ⟨| Γ r | p ⟩ are shown in Fig. 4. The local slopes also exhibit poor scaling properties, which are improved by employing the ESS method with respect to the second-order moment ⟨| Γ r | 2 ⟩ . The insets of the bottom panels in Fig. 4 show the deviations of the local slopes with respect to the self-similar prediction 4 p/ 3.

Figure 5 summarizes the results of this work. The left panels show the longitudinal and transverse scaling exponents obtained from both databases using the ESS method, exhibiting a self-similar behavior for the former and an anomalous scaling for the latter. The error bars indicate the maximum and minimum values of the local slopes in the inertial range. The right panel shows the circulation scaling exponents λ p obtained from the ESS method, and we compare them with the transverse exponents via ζ ⊥ p + . p We also include experimental data extracted from Zhu et al. [18], obtained without employing the ESS method. The insets show the anomalous deviations of the scaling exponents with respect to the self-similar predictions ζ ss p = p/ 3 and λ ss p = 4 p/ 3, i.e.

$$\frac { \zeta _ { p } ^ { \| } - \zeta _ { p } ^ { \text{ss} } } { \zeta _ { p } ^ { \text{ss} } }, \ \frac { \zeta _ { p } ^ { \perp } - \zeta _ { p } ^ { \text{ss} } } { \zeta _ { p } ^ { \text{ss} } }, \ \text{ and } \frac { \lambda _ { p } - \lambda _ { p } ^ { \text{ss} } } { \lambda _ { p } ^ { \text{ss} } }.$$

Despite the clear differences observed in the forcing scales between the two numerical simulations, the scaling exponents obtained from the ESS method in the inertial range are statistically equivalent. The different trend observed in the inset of Fig.5a might be due to statistical errors and it is less than 3%. Remarkably, both datasets lead to the same transverse anomalous scaling. Moreover, circulation scaling exponents measured in experiments coincide well with the

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-d441f213-3262-4d8e-b6b0-e1c91e65f18b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-4fdf2b1e-75b9-4841-9272-bc6ef5019e71.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-f3d09d35-adc4-4456-895f-6a0ff5fe17c8.jpg)

〈

|

|

〉

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-04dd94c3-a605-4b01-8390-53803f04dcae.jpg)

〈

|

|

〉

FIG. 4. Circulation moments up to p = 10 in numerical simulations of 2D turbulence in the inverse energy cascade for (left column) RUN-A and (right column) RUN-B. Top panels show the circulation moments as a function of r , and the bottom panels show the ESS method. The shaded area indicates the inertial range. Dashed horizontal lines in the insets show the self-similar scaling λ p = 4 p/ 3. The top insets show the local slopes λ p ( r ) = d log Γ ⟨ p r ⟩ / dlog r , and the bottom insets show the deviations from the self-similar scaling.

ones measured in both simulations without using the ESS method. These results show the robustness of the turbulent properties of the inverse cascade in 2D flows, at least when the ESS method is used.

Starting from Eq. (3) relating circulation and transverse velocity increments, it is possible to show using the H¨lder o inequality that (see Supplemental Material [20])

$$\lambda _ { p } \leq \zeta _ { p } ^ { \perp } + p.$$

This inequality implies that, if the TSFs are intermittent, i.e. ζ ⊥ p &lt; p/ 3 for p &gt; 3, then the circulation moments also are, as it follows that λ p &lt; 4 p/ 3 for p &gt; 3. This relation is obeyed in our simulations, as shown in Fig. 5.c. We remark that the derivation of Eq. (8) is valid only in the limit r/L f ≫ 1 for the inverse energy cascade. In 3D turbulence, using data extracted from Iyer et al. [11] and Iyer et al. [22], we observe that the inequality still holds, however, the demonstration of Eq. (8) needs to be modified.

To conclude, we analyzed the statistical properties of two-dimensional turbulence from direct numerical simulations of the incompressible Navier-Stokes equations of two different databases. By forcing at small scales, a stationary inverse energy cascade is generated which exhibits well-known features, such as the k -5 / 3 scaling law for the energy spectrum and a negative energy flux. The study of high-order moments structure functions confirmed the well-known self-similar behavior of longitudinal velocity increments but revealed anomalous deviations for the transverse ones. We argue that the intermittent behavior of this quantity explains the anomalous exponents of velocity circulation through the inequality (8). For a future work we propose to investigate the inverse energy cascade in a dual cascade regime to study how non-local effects of the enstrophy cascade affect large-scale fluctuations. We also propose to study more realistic geophysical quasi-2D flows such as rotating or stratified turbulence. The characterization of structure functions and circulation moments in these systems could bring some insights into anisotropic effects and could be useful for atmospheric modeling.

FIG. 5. Scaling exponents in the inverse energy cascade of 2D turbulence. Panels (a) and (b) show the longitudinal and transverse scaling exponents, respectively, for both numerical simulations up to p = 10. Panel (c) shows the circulation exponents λ p compared to the transverse exponents ζ ⊥ p + . p The insets show the anomalous deviations with respect to the self-similar prediction p/ 3 for velocity increments and 4 p/ 3 for circulation defined in Eq. (7). Error bars are defined as the maximum and minimum values of the local slopes, which are not shown in the insets. Black crosses show experimental scaling exponents extracted from [18].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-947172d1-d67a-4376-80ff-f10de319d80d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-c51db2fb-f57f-414a-8cb3-5e5ec699b0dc.jpg)

## ACKNOWLEDGMENTS

We are grateful to Guido Boffetta for fruitful scientific discussions and for kindly sharing his 2D turbulence data. This work was supported by the Agence Nationale de la Recherche through the project GIANTE ANR-18-CE30-002001. GK was also funded by the Simons Foundation project 'Collaboration in Wave Turbulence' (award ID 651471) Computations were carried out at the M´socentre SIGAMM hosted at the Observatoire de la Cˆte d'Azur. e o

## Appendix: Relation between circulation and transverse scaling exponents

We consider the circulation Γ r ( s 0 ) around a squared loop of size r , with one corner of the loop placed at s 0 = ( x , y 0 0 ). It follows that

$$& \left ( | \Gamma _ { r } | ^ { p } \right ) = \frac { 1 } { V } \int | \Gamma _ { r } ( s _ { 0 } ) | ^ { p } d s _ { 0 } \leq \frac { 1 } { V } \int \left [ \int _ { x _ { 0 } } ^ { z _ { 0 } + r } | u _ { x } ( x, y _ { 0 } + r ) - u _ { x } ( x, y _ { 0 } ) | d x + \int _ { y _ { 0 } } ^ { y _ { 0 } + r } | u _ { y } ( z _ { 0 } + r, y ) - u _ { y } ( z _ { 0 }, y ) | d y \right ] ^ { p } d s _ { 0 }, \\ & \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod:                                                                                                                                                                                                        
\quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod:                                                                                                                                                 
\quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{ mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod:                                                                                                                                                                      
\quad \text{mod: } \quad \text{mod: } \quad \text{mod: } \quad \text{mod: }                                                                                                                                                                                                        \quad \text{mod:                                                                                                                                                                                                
\quad \text{mod:- } \quad \text{mod:                                                                                                                                                                                               
\quad \text{mod:-                                                                                                                                                                                                        }                                                                                                                                                                                                       }
\quad \text{mod:-                                                                                                                                                                                                }        }                                                                                                                                                                                                }       
\quad \text{mod:-                                                                                                                                                                                         }               }                                                                                                                                                                                         }              
\quad \text{mod:-                                                                                                                                                                                  }                      }                                                                                                                                                                                  }                     
\quad \text{mod:-                                                                                                                                                                           }                             }                                                                                                                                                                           }                            
\quad \text{mod:-$$

where we applied the triangular inequality several times. Following a similar procedure to the one in Iyer et al. [11], we now apply the H¨lder inequality for each of the two integrals in the right hand side. o It leads to

$$\cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \\ & \langle \Gamma _ { r } | ^ { p } \rangle \leq \frac { 1 } { V } \int \int _ { x _ { 0 } } \left [ r ^ { 1 / q } \left ( \int _ { x _ { 0 } } ^ { x _ { 0 } + r } | u _ { x } ( x, y _ { 0 } + r ) - u _ { x } ( x, y _ { 0 } ) | ^ { p } d x \right ) ^ { 1 / p } + r ^ { 1 / q } \left ( \int _ { y _ { 0 } } ^ { y _ { 0 } + r } | u _ { y } ( x _ { 0 } + r, y ) - u _ { p } ( x _ { 0 }, y ) | ^ { p } d y \right ) ^ { 1 / p } \right ] \, d s _ { 0 } \\ \text{with $n$ and $a$ as of $i$ in $n^{-1} \cdot n^{-1} - 1$ $f$ for $n$} \, \alpha \, \beta \, \cdots \, 1$$

with p and q satisfying p -1 + q -1 = 1 for p, q &gt; 1.

For a sufficiently large Reynolds numbers, assuming homogeneity, isotropy, and at a fixed r in the inertial range, we can approximate each inner integral by r ⟨| δu ⊥ r | p ⟩ = rS ⊥ p ( r ). Fig. 6 shows the validity of this approximation in the inertial range for RUN-A. The outer integral cancels out and we obtain

$$\langle | \Gamma _ { r } | ^ { p } \rangle & \leq 2 ^ { p } r ^ { p / q } r ^ { p / p } ( S _ { p } ^ { \perp } ) ^ { p / p } = 2 ^ { p } r ^ { p } S _ { p } ^ { \perp } ( r ).$$

Finally, we use the fact that the circulation moments and TSFs follow the scaling properties ⟨| Γ r | p ⟩ ∼ ( r/L f ) λ p and S ⊥ p ∼ ( r/L f ) ζ ⊥ p , with L f the forcing scale. For the inertial range of the inverse energy cascade in two-dimensional turbulence, we take the limit r/L f ≫ 1, so we obtain an inequality for the scaling exponents

$$\lambda _ { p } \leq \zeta _ { p } ^ { \perp } + p.$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00631v2-figure-ac99da28-c4c4-486a-b531-5a3f54e18056.jpg)

r

FIG. 6. Validation of the approximation performed between Eqs.(10) and (11) for p = 2, with u = u x and v = u y . The angular brackets ⟨ ⟩ . indicate averaging in space.

- [1] R. H. Kraichnan, Inertial Ranges in Two-Dimensional Turbulence, Physics of Fluids 10 , 1417 (1967).
- [2] C. E. Leith, Diffusion Approximation for Two-Dimensional Turbulence, Physics of Fluids 11 , 671 (1968).
- [3] G. K. Batchelor, Computation of the Energy Spectrum in Homogeneous Two-Dimensional Turbulence, Physics of Fluids 12 , II (1969).
- [4] P. A. Davidson, Turbulence in Rotating, Stratified and Electrically Conducting Fluids (Cambridge University Press, Cambridge, 2013).
- [5] A. Alexakis and L. Biferale, Cascades and transitions in turbulent flows, Physics Reports 767-769 , 1 (2018).
- [6] R. Marino, P. D. Mininni, D. Rosenberg, and A. Pouquet, Inverse cascades in rotating stratified turbulence: Fast growth of large scales, EPL (Europhysics Letters) 102 , 44006 (2013).
- [7] R. M. B. Young and P. L. Read, Forward and inverse kinetic energy cascades in Jupiter's turbulent weather layer, Nature Physics 13 , 1135 (2017).
- [8] U. Frisch, Turbulence: The Legacy of A.N. Kolmogorov , 1st ed. (Cambridge University Press, Cambridge, 1995).
- [9] J. Paret and P. Tabeling, Intermittency in the two-dimensional inverse cascade of energy: Experimental observations, Physics of Fluids 10 , 3126 (1998).
- [10] G. Boffetta, A. Celani, and M. Vergassola, Inverse energy cascade in two-dimensional turbulence: Deviations from Gaussian behavior, Physical Review E 61 , R29 (2000).
- [11] K. P. Iyer, K. R. Sreenivasan, and P. K. Yeung, Circulation in High Reynolds Number Isotropic Turbulence is a Bifractal, Physical Review X 9 , 041006 (2019).
- [12] N. P. M¨ller, J. I. Polanco, and G. Krstulovic, Intermittency of Velocity Circulation in Quantum Turbulence, Physical u Review X 11 , 011053 (2021).
- [13] J. I. Polanco, N. P. M¨ller, and G. Krstulovic, Vortex clustering, polarisation and circulation intermittency in classical u and quantum turbulence, Nature Communications 12 , 7090 (2021).
- [14] N. P. M¨ uller, Y. Tang, W. Guo, and G. Krstulovic, Velocity circulation intermittency in finite-temperature turbulent superfluid helium, Physical Review Fluids 7 , 104604 (2022).
- [15] A. Migdal, Clebsch confinement and instantons in turbulence, International Journal of Modern Physics A 35 , 2030018 (2020).
- [16] G. B. Apolin´rio, L. Moriconi, R. M. Pereira, and V. J. Valad˜o, Vortex gas modeling of turbulent circulation statistics, a a Physical Review E 102 , 041102(R) (2020).
- [17] L. Moriconi, R. M. Pereira, and V. J. Valad˜o, Vortex polarization a and circulation statistics in isotropic turbulence, Physical Review E 109 , 045106 (2024).
- [18] H.-Y. Zhu, J.-H. Xie, and K.-Q. Xia, Circulation in Quasi-2D Turbulence: Experimental Observation of the Area Rule and Bifractality, Physical Review Letters 130 , 214001 (2023).

- [19] N. P. M¨ller and G. Krstulovic, Exploring the Equivalence between Two-Dimensional Classical and Quantum Turbulence u through Velocity Circulation Statistics, Physical Review Letters 132 , 094002 (2024).
- [20] See Supplemental Material at URL-will-be-inserted-by-publisher for details on the analytical derivation of Eq. (8).
- [21] S. Chen, K. R. Sreenivasan, M. Nelkin, and N. Cao, Refined Similarity Hypothesis for Transverse Structure Functions in Fluid Turbulence, Physical Review Letters 79 , 2253 (1997).
- [22] K. P. Iyer, K. R. Sreenivasan, and P. K. Yeung, Scaling exponents saturate in three-dimensional isotropic turbulence, Physical Review Fluids 5 , 054605 (2020).
- [23] G. Boffetta and R. E. Ecke, Two-Dimensional Turbulence, Annual Review of Fluid Mechanics 44 , 427 (2012).
- [24] G. Boffetta, Energy and enstrophy fluxes in the double cascade of two-dimensional turbulence, Journal of Fluid Mechanics 589 , 253 (2007).
- [25] R. Benzi, S. Ciliberto, R. Tripiccione, C. Baudet, F. Massaioli, and S. Succi, Extended self-similarity in turbulent flows, Physical Review E 48 , R29 (1993).
Nicolás Pablo Müller Giorgio Krstulovic
Lack of self-similarity in transverse velocity increments and circulation statistics in two-dimensional turbulence
We study the statistics of longitudinal and transverse structure functions,
as well as velocity circulation in the inverse energy cascade of
two-dimensional turbulence. By means of direct numerical simulations of the
incompressible Navier-Stokes equations, we show that transverse structure
functions exhibit an anomalous scaling, in contrast to the self-similar
behavior of longitudinal ones. We derive an analytical relation that shows that
the scaling exponents of transverse structure functions and velocity
circulation are related in two-dimensional turbulence.