Universality and Thermodynamics of Turbulence

TL;DR

This paper demonstrates a universal log-scale behavior in turbulence velocity structure functions across experiments and simulations, linking it to multifractal thermodynamics and phase transitions in velocity fields.

Contribution

It introduces the concept of log-universality in Eulerian velocity structure functions and connects it with multifractal thermodynamics and phase transitions.

Findings

Velocity structure functions follow a log-universal curve.

Log-universality relates to multifractal free energy properties.

A phase transition in velocity fields is associated with fluctuating dissipation scales.

Abstract

We investigate universality of the Eulerian velocity structure functions using velocity fields obtained from the stereoscopic particle image velocimetry (SPIV) technique in experiments and the direct numerical simulations (DNS) of the Navier-Stokes equations. We show that the numerical and experimental velocity structure functions up to order 9 follow a log-universality; we find that they collapse on a universal curve, if we use units that include logarithmic dependence on the Reynolds number. We then investigate the meaning and consequences of such log-universality, and show that it is connected with the properties of a "multifractal free energy", based on an analogy between multifractal and themodynamics. We show that in such a framework, the existence of a fluctuating dissipation scale is associated with a phase transition describing the relaminarisation of rough velocity fields with different Holder exponents. Such a phase transition has been already observed using the Lagrangian velocity structure functions, but was so far believed to be out of reach for the Eulerian data.

Tables

Table 1. Table 1: Parameters for the 5 experiments realized ( 𝐀 , 𝐁 , 𝐂 , 𝐃 ​ a ​ n ​ d ​ 𝐄 ) 𝐀 𝐁 𝐂 𝐃 𝑎 𝑛 𝑑 𝐄 (\mathbf{A},\mathbf{B},\mathbf{C},\mathbf{D}\ and\ \mathbf{E}) . F is the rotation frequency of the discs, Re refers to the Reynolds number based on the diameter of the tank, R λ subscript 𝑅 𝜆 R_{\lambda} is the Reynolds based on the Taylor micro-scale. η 𝜂 \eta gives the estimated Kolmogorov length according to the experiment and Δ ​ x Δ 𝑥 \Delta x refers to the spatial resolution of SPIV measurements. The last columns gives the number of frames used and the number of points over which are calculated the statistics. Except for ( E ), the Reynolds are much larger than those available with DNS. Table taken from Saw  et al. ( 2018 )

Case	F (Hz)	Glycerol part	Re	$R_{\lambda }$	$\eta$ (mm)	$\Delta x$	Frames	Symbol
A	5	0%	$3\times {10}^5$	$1,9\times {10}^3$	0.02	2.4	$3\times {10}^4$	$\circ$
B	5	0%	$3\times {10}^5$	$2,7\times {10}^3$	0.02	0.48	$3\times {10}^4$	$\square$
C	5	0%	$3\times {10}^5$	$2,5\times {10}^3$	0.02	0.24	$2\times {10}^4$	$◇$
D	1	0%	$4\times {10}^4$	$9,2\times {10}^2$	0.08	0.48	$1\times {10}^4$	$△$
E	1.2	59%	$6\times {10}^3$	$2,1\times {10}^2$	0.37	0.24	$3\times {10}^4$	$\star$

Table 2. Table 2: Parameters for the DNS. R λ subscript 𝑅 𝜆 R_{\lambda} is the Reynolds based on the Taylor micro-scale. η 𝜂 \eta is the Kolmogorov length. The third column gives resolution of the simulation through k max ​ η subscript 𝑘 𝜂 k_{\max}\eta ,where k max = N / 3 subscript 𝑘 𝑁 3 k_{\max}=N/3 is the maximum wavenumber. The fourth column gives the grid size; notice that the length of the box is 2 ​ π 2 𝜋 2\pi . Here, ℓ min subscript ℓ \ell_{\min} is the smallest scale available for the calculations of the wavelets. k ~ f subscript ~ 𝑘 𝑓 \tilde{k}_{f} is the forcing scale. The Sample columns gives the number of points (frames × \times gridsize) over which are calculated the statistics.

$R_{\lambda }$	$\eta$	$k_{\max }\eta$	$N_x\times N_y\times N_z$	${\ell }_{\min }/\eta$	$\widetilde{k_f}$	Samples	Symbol
25	0.079	3.35	${128}^3$	0.635	1	5000	$\star$
56	0.034	6.42	${256}^3$	0.31	1	105 000	$△$
80	0.020	1.68	${256}^3$	1.22	1	270 000	$\square$
90	0.017	5.70	${1024}^3$	0.36	1	10 000	$◇$
138	0.009	1.55	${512}^3$	1.37	1	12 000	$\circ$

Table 3. Table 3: Scaling exponents τ ​ ( p ) 𝜏 𝑝 \tau(p) and ζ ​ ( p ) 𝜁 𝑝 \zeta(p) found by the collapse method based on K62 universality for experimental data (subscript EXP) or numerical data (subscript DNS). The subscript SAW refers to the values obtained by Saw  et al. ( 2018 ) on the same set of experimental data, using velocity increments and Extended Self-Similarity technique Benzi  et al. ( 1993 ) . The exponents τ EXP ​ ( p ) subscript 𝜏 EXP 𝑝 \tau_{\text{EXP}}(p) (red square) and τ DNS subscript 𝜏 DNS \tau_{\text{DNS}} (blue circle) have been computed through a least square algorithm upon τ ​ ( p ) 𝜏 𝑝 \tau(p) , minimizing the scatter of the rescaled structure functions log [ ( L 0 η ) τ ​ ( p ) S ~ p ) ] \log\left[\left(\frac{L_{0}}{\eta}\right)^{\tau(p)}\tilde{S}_{p})\right] with respect to the line ( ℓ / η ) τ ​ ( p ) superscript ℓ 𝜂 𝜏 𝑝 (\ell/\eta)^{\tau(p)} . The corresponding ζ ​ ( p ) 𝜁 𝑝 \zeta(p) were inferred using the formula ζ ​ ( p ) = τ ​ ( p ) + ζ ​ ( 3 ) ​ p / 3 𝜁 𝑝 𝜏 𝑝 𝜁 3 𝑝 3 \zeta(p)=\tau(p)+\zeta(3)p/3 , where ζ ​ ( 3 ) 𝜁 3 \zeta(3) is computed in figure 4 -a.

exponent \ order	$p=1$	$p=2$	$p=3$	$p=4$	$p=5$	$p=6$	$p=7$	$p=8$	$p=9$
${\zeta }_{\text{SAW}}/{\zeta }_{\text{SAW}}\left(3\right)$	0.36	0.69	1	1.29	1.55	1.78	1.98	2.17	2.33
${\zeta }_{\text{DNS}}$	0.31	0.58	0.80	0.98	1.12	1.23	1.26	1.25	1.23
${\zeta }_{\text{EXP}}$	0.32	0.58	0.80	0.98	1.12	1.23	1.32	1.39	1.44
${\tau }_{\text{DNS}}$	0.04	0.05	0	-0.09	-0.21	-0.37	-0.61	-0.88	-1.17
${\tau }_{\text{EXP}}$	0.05	0.05	0	-0.09	-0.21	-0.36	-0.54	-0.74	-0.96

Table 4. Table 4: Summary of the analogy between the multifractal formalism of turbulence and thermodynamics.

	Thermodynamics	Turbulence
Temperature	$k_{B}T$	$1/p$
Energy	$E$	$\log \left({\mu }_{\ell }\right)$
Number of d.f.	$N$	$\log \left(\text{Re}\right)\equiv \log \left(L_0/\eta \right)/{\beta }_0$
Volume	$V$	$\log \left(\ell /\eta \right)$
Pressure	$P$	$\tau (p,\ell )$
Free energy	$F$	$\log \left({\tilde{S}}_{3p}\right)$

Equations

$$G_{ij}(\text{x},\ell )=\int \mathrm{d}\text{r}{\nabla }_j{\Phi }_{\ell }\left(\mathbf{r}\right)u_i(\mathbf{x}+\mathbf{r}),$$

$$\delta W(\text{u})(\text{x},\ell )=\ell \underset{ij}{\max }|G_{ij}(\text{x},\ell )|.$$

$$S_p(\ell )\equiv ⟨|{\delta }_{\ell }u{|}^p⟩\sim u_{\mathrm{K}}^p{\left(\frac{\ell }{\eta }\right)}^{p/3},$$

$${\tilde{S}}_p\equiv \frac{S_p}{S_3^{p/3}}=C_p,$$

$$S_p(\ell )\sim (ϵ\ell {)}^{p/3}{\left(\frac{\ell }{L}\right)}^{\mu p(3-p)},$$

$${\tilde{S}}_p\equiv \frac{S_p}{S_3^{p/3}}\sim A_p{\left(\frac{\ell }{L}\right)}^{\tau (p)},$$

$${\left(\frac{L_0}{\eta }\right)}^{\tau (p)}{\tilde{S}}_p\sim A_p{\left(\frac{\ell }{\eta }\right)}^{\tau (p)}.$$

$$h\left(\text{x},t,\ell \right)=\frac{\log \left(|{\delta }_{\ell }\text{u}(\text{x},t)|/u_0\right)}{\log (\ell /L)},$$

$$\mathbb{P}\left[\log ({\delta }_{\ell }u/u_0)=h\log (\ell /L)\right]\sim e^{\log (\ell /L)C(h)}={\left(\frac{\ell }{L}\right)}^{C(h)},$$

$$S_p(\ell )=⟨({\delta }_{\ell }u{)}^p⟩=\underset{h_{\min }}{\overset{h_{\max }}{\int }}{u}_0^p{\left(\frac{\ell }{L}\right)}^{ph+C(h)}\mathrm{d}h.$$

$$S_{p}(\ell)=(\epsilon\ell)^{p/3}\int\limits_{h_{\min}}^{h_{\max}}\Big{(}\frac{\ell}{L}\Big{)}^{p(h-1/3)+C(h)}\mathrm{d}h\sim(\epsilon\ell)^{p/3}\Big{(}\frac{\ell}{L}\Big{)}^{\tau(p)}.$$

$$S_{p}=(\epsilon\ell)^{p/3}\int\limits_{\eta_{h}\leq\ell}\Big{(}\frac{\ell}{L}\Big{)}^{p(h-1/3)+C(h)}\mathrm{d}h\sim(\epsilon\ell)^{p/3}\Big{(}\frac{\ell}{L}\Big{)}^{\tau(p,\ell)}.$$

$$\tilde{S}_{p}\equiv\frac{S_{p}}{S_{3}^{p/3}}=\int\limits_{-1+1/\theta}^{h_{\max}}\Big{(}\frac{\ell}{L}\Big{)}^{p(h-1/3)+C(h)}\mathrm{d}h\sim\exp\left(-\theta{\tau(p,\theta)\log(\text{Re})}\right),$$

$$\frac{\log ({\tilde{S}}_p)}{\log (\text{Re})}=-\tau (p,\theta )\theta ,\theta =\log (L/\ell )/\log (\text{Re}).$$

$$\gamma (\text{Re})\log \left(\frac{S_p}{A_p{u}_{\mathrm{K}}^p}\right)=G\left(p,\gamma (\text{Re})\log (\ell K_0/\eta )\right),$$

$$\beta (\text{Re})\left(\frac{\log ({\tilde{S}}_p/S_{0p})}{\log (L_0/\eta )}\right)=G\left(p,\beta (\text{Re})\frac{\log (\ell /\eta )}{\log (L_0/\eta )}\right),$$

$$S_p=⟨|\delta W(\text{u})(\text{x},\ell ){|}^p⟩.$$

$$\log \left(\frac{S_p}{u_{\mathrm{K}}^p}\right)=F\left(\log \left(\frac{\ell }{\eta }\right)\right).$$

$$\log \left[{\left(\frac{L_0}{\eta }\right)}^{\tau (p)}{\tilde{S}}_p\right]=F\left(\log \left(\frac{\ell }{\eta }\right)\right).$$

$$\beta (\text{Re})\left(\frac{\log ({\tilde{S}}_p/S_{0p})}{\log (L_0/\eta )}\right)=F\left(p,\beta (\text{Re})\frac{\log (\ell /\eta )}{\log (L_0/\eta )}\right),$$

$$C(h)=\underset{p}{\min }[ph+\zeta (p))].$$

$$C\Big{(}\frac{\mathrm{d}\,\zeta(p)}{\mathrm{d}\,p}\Big{|}_{p^{*}}\Big{)}=\zeta(p^{*})-p^{*}\frac{\mathrm{d}\,\zeta(p)}{\mathrm{d}\,p}\Big{|}_{p^{*}}.$$

$$C(h)=\frac{(h-a{)}^2}{2b},$$

$${\mu }_{\ell }=\frac{|\delta W_{\ell }{|}^3}{⟨|\delta W_{\ell }{|}^3⟩}.$$

$$\mathbb{P}\left[\log ({\mu }_{\ell })=E\log (\ell /\eta )\right]\sim e^{\log (\ell /\eta )S(E)},$$

$${\tilde{S}}_{3p}=\frac{S_{3p}}{S_3^p}=⟨e^{p\log ({\mu }_{\ell })}⟩,$$

$$\log ({\tilde{S}}_{3p})=F.$$

$$F(T,V,N)=NF(T,\frac{V}{N},1),$$

Full text

I Introduction

A well-known feature of any turbulent flow is the Kolmogorov-Richardson cascade by which energy is transferred from large to small length scales until the Kolmogorov length scales below which it is removed by viscous dissipation. This energy cascade is a non-linear and an out-of-equilibrium universal process. Moreover, the corresponding non-dimensional energy spectrum $E(k)/\epsilon^{2/3}\eta^{5/3}$ is an universal function of $k\eta$, where $\eta=(\nu^{3}/\epsilon)^{1/4}$ is the Kolmogorov length scale, $\epsilon$ the mean energy dissipation rate per unit mass, and $\nu$ the kinematic viscosity. However, there seem to be little dependences on the Reynolds number, boundary, isotropy or homogeneity conditions Dubrulle (2019). In facts, the energy spectrum is based upon a quantity, the velocity correlation, that is quadratic in velocity. Nevertheless, it is now well admitted that the universality does not carry over for statistical quantities that involve higher order moments. For example, the velocity structure functions of order $p$, given by $S_{p}(\ell)=\langle\|u(\textbf{x}+\textbf{r})-u(\textbf{x})\|^{p}\rangle_{\textbf{x},\|\textbf{r}\|=\ell}$ are not universal, at least when expressed in units of the Komogorov scale $\eta$ and velocity $u_{\mathrm{K}}=(\nu\epsilon)^{1/4}$(see below, section III.2 for an illustration).

The mechanism behind this universality breaking was identified by [P85], wherein a generalization of the Kolmogorov theory introduced, based on the hypothesis that a turbulent flow is multifractal. In this model, the velocity field is characterized locally by an exponent $h$, such that $|\delta_{\ell}u(\textbf{x})|\equiv\langle\|\textbf{u}(\textbf{x}+\textbf{r})-\textbf{u}(\textbf{x})\|\rangle_{\|\textbf{r}\|=\ell}\sim\ell^{h(\textbf{x})}$; here $h$ is a stochastic function that follows a large deviation property Eyink (2008) $\mathbb{P}\left(\log(|\delta_{\ell}u|/u_{0})=h\log\left(\ell/L_{0}\right)\right)\sim\left(\ell/L_{0}\right)^{C(h)}$, where $u_{0}$ (resp. $L_{0}$) is the caracteristic integral velocity (resp. length), and $C(h)$ is the multifractal spectrum. Velocity fields with $h<1$ are rough in the limit $\ell\to 0$. In real flows, any rough field with $h>-1$ can be regularized at sufficiently small scale (the ”viscous scale”) by viscosity. The first computation of such dissipative scale was performed by Paladin and Vulpiani [P87], who showed that it scales with viscosity like $\eta_{h}\sim\nu^{1/(1+h)}$, thereby generalizing the Kolmogorov scale, which corresponds to $h=1/3$. Such a dissipative scale fluctuates in space and time (along with $h$), resulting in non-universality for high order moments, at least when expressed in units of $\eta$ and $u_{\mathrm{K}}$.

A few years later, Frisch and Vergassola Frisch and Vergassola (1991) claimed that the universality of the energy spectrum can be recovered, if the fluctuations of the dissipative length scale are taken into account by introducing a new non-dimensionalisation procedure. The new prediction is that $\log\left(E(k)\epsilon^{-\frac{2}{3}}\eta^{-\frac{5}{3}}\right)/\log(\text{Re})$ should be a universal function of $\log(k\eta)/\log(\text{Re})$, where Re is the Reynolds number. This claim was examined by Gagne et al., later using data from the Modane wind tunnel experiments Castaing et al. (1993). They further suggested that the prediction can be extended to the velocity structure functions, so that, at any given $p$ $\log(S_{p}(\ell)/u_{\mathrm{K}}^{p})/\log(\text{Re})$ should be a universal function of $\log(\ell/\eta)/\log(\text{Re})$; they found good agreement for $p$ up to $6$. The velocity measurements, in the above experiments, were performed using hot wire anemometry, which provide access to only one component of velocity. To our knowledge, no further attempts have been made to check the claim with more realistic measurements.

The purpose of the present paper is to reexamine this claim; however, now using the velocity fields obtained from the Stereoscopic Particule Image Velocimetry (SPIV) in experiments and the direct numerical simulations (DNS) of the Navier-Stokes equations (NSE). We show that the numerical and experimental velocity structure functions up to order 9 follow a log-universality Castaing et al. (1993); they indeed collapse on a universal curve, if we use units that include $\log(\text{Re})$ dependence. We then investigate the meaning and consequences of such a log-universality, and show that it is connected with the properties of a ”multifractal free energy”, based on an analogy between multifractal and thermodynamics (see [M91] for summary). We show that in such a framework, the existence of a fluctuating dissipation length scale is associated with a phase transition describing the relaminarisation of rough velocity fields with different Hölder exponents.

II Experimental and numerical setup

II.1 Experimental facilities and parameters

We use experimental velocity field described in Saw et al. (2018). The radial, axial and azimuthal velocity are measured in a von Kármán flow, using Stereoscopic Particule Image Velocimetry technique at different resolutions $\Delta x$. The von Kármán flow is generated in a cylindrical tank of radius $R=$10\text{,}\mathrm{c}\mathrm{m}$$ through counter-rotation of two independent impellers with curved blades. The flow was maintained in a turbulent state at high Reynolds number by two independent impellers, rotating at frequency $F$. Figure 1 shows the sketch of the experimental setup. The five experiments were performed in conditions so that the non-dimensional mean energy dissipation per unit mass is constant. The viscosity was monitored using mixture of water and glycerol, so as to vary the Kolmogorov length $\eta$. Table 1 summarizes the different parameters; $R_{\lambda}=\lambda u^{\mathrm{rms}}/\nu$ is the Reynolds number based on the Taylor length scale $\lambda=\sqrt{\frac{\langle\textbf{u}^{2}\rangle}{\langle\nabla\textbf{u}^{2}\rangle}}$, the mean squared velocity $u^{\mathrm{rms}}$ and the viscosity $\nu$.

All velocity measurements are performed in a vertical plane that contains the rotation axis. The case A corresponds to measurements over the whole plane contained in between the two impellers, and extending from one side to the other side of the cylinder. Its resolution is 5 to 10 times coarser than similar measurements performed by zooming on a region centered around the symmetry point of the experiment (on the rotation axis, half way in between the two impellers), over a square window of size 2\text{\,}\mathrm{c}\mathrm{m}$\times$2\text{\,}\mathrm{c}\mathrm{m}.

II.2 Direct Numerical Simulation

The direct numerical simulations (DNS), based on pseudo-spectral method, were performed in order to compare with our experimental data. The DNS runs with $R_{\lambda}=25$, $R_{\lambda}=80$, $R_{\lambda}=90$ and $R_{\lambda}=138$ were performed using the NSE solver VIKSHOBHA Debue et al. (2018), whereas the run with $R_{\lambda}=56$ was carried out using another independent pseudo-spectral method based NSE solver. The velocity field $\bm{u}$ was computed on a $2\pi$ triply-periodic box.

Turbulent flow in a statistically steady state was obtained by using the Taylor-Green type external forcing in the NSE at wavenumber $k_{f}=1$ and amplitude $f_{0}=0.12$, the value of viscosity was varied in order to obtain different values of $R_{\lambda}$ (see Ref. Debue et al. (2018) for more details).

III Theoretical background

III.1 Velocity increments vs Wavelet Transform (WT) of velocity gradients

The classical theories of Kolmogorov Kolmogorov (1941); [K62] are based on the scaling properties of the velocity increment, defined as $\delta_{\ell}\textbf{u}={\bf u}({\bf x+r})-{\bf u}({\bf x})$, where $\ell=|{\bf r}|$ is the distance over which the increment is taken. As pointed out by [M91], a more natural tool to characterize the local scaling properties of the velocity field is the wavelet transform of the tensor $\partial_{j}u_{i}$, defined as:

[TABLE]

where $\Phi_{\ell}(\textbf{x})=\ell^{-3}\Phi(\textbf{x}/\ell)$ is a smooth function, non-negative with unit integral. In what follows, we choose a Gaussian function $\Phi(\textbf{x})=\exp(-\|\textbf{x}\|^{2}/2)/(2\pi)^{\frac{3}{2}}$ such that $\int\Phi(\textbf{r})\mathrm{d}\textbf{r}=1$. We then compute the wavelet velocity increments as

[TABLE]

This formulation is especially well suited for the analysis of the experimental velocity field, as it naturally allows to average out the noise. It has been verified that the wavelet based approach yields the same values for the scaling exponents as those computed from the velocity increments Debue et al. (2018).

III.2 K41 and K62 universality

In the first theory of Kolmogorov Kolmogorov (1941), the turbulence properties depend only on two parameters: the mean energy dissipation per unit mass $\epsilon$ and the viscosity $\nu$ . The only velocity and length unit that one can build using these quantities are the Kolmogorov length $\eta=(\nu^{3}/\epsilon)^{1/4}$ and velocity $u_{\mathrm{K}}=(\epsilon\nu)^{1/4}$. The structure functions are then self-similar in the inertial range $\eta\ll\ell\ll L_{0}$, where $L_{0}$ is the integral scale, and follow the universal scalings:

[TABLE]

which can also be recast into:

[TABLE]

where $C_{p}$ is a (non universal) constant.

This scaling is typical of a global scale symmetry solutions, and was criticized by Landau, who considered it incompatible with observed large fluctuations of the local energy dissipation. Kolmogorov then built a second theory (K62), in which fluctuations of energy dissipation were assumed to follow a log-normal statistics, and taken into account via an intermittency exponent $\mu$ and a new length scale $L$, thereby breaking the global scale invariance. The resulting velocity structure functions then follow the new scaling:

[TABLE]

which implies a new kind of universality involving the relative structure functions $\tilde{S}_{p}$ as:

[TABLE]

where $\tau(p)=\mu p(3-p)$ and $A_{p}$ is a constant. Such a formulation already predicts an interesting universality, if $L=L_{0}$, as we should have:

[TABLE]

Therefore, we should be able to collapse all structure functions, at different Reynolds number by plotting $(\frac{L_{0}}{\eta})^{\tau(p)}\tilde{S}_{p}$ as a function of $\frac{\ell}{\eta}$, given that $L_{0}/\eta\sim\text{Re}^{3/4}$. There is however no clear prediction about the value of $L$ and we show in the data analysis section that $L$ differs from $L_{0}$.

The relation (7) shows that $\log\left(\left(\frac{L_{0}}{\eta}\right)^{\tau(p)}\tilde{S}_{p}\right)$ is a linear function of $\log(\frac{\ell}{\eta})$. In principle, such universal scaling is not valid outside the inertial range, i.e. for example when $\ell<\eta$. To be more general than previously thought, it can however be shown using the multifractal formalism as first shown by Frisch and Vergassola (1991).

III.3 Multifractal and fluctuating dissipation length

For the multifractal (MFR) model, it is assumed that the turbulence is locally self-similar, so that there exists a scalar field $h(\textbf{x},\ell,t)$, such that

[TABLE]

for a range of scales in a suitable ”inertial range” $\eta_{h}\ll\ell\ll L$, where $L$ is a characteristic integral-length-scale, $\eta_{h}$ a cut-off length scale, and $u_{0}$ a characteristic large-scale velocity. This scale is a generalization of the Kolmogorov scale, and is defined as the scale where the local Reynolds number $\ell|\delta_{\ell}\bm{u}|/\nu$ is equal to $1$. Writing $\delta_{\ell}u=|\delta_{\ell}\textbf{u}|=u_{0}(\ell/L)^{h}$ leads to the expression of $\eta_{h}$ as a function of the global Reynolds number $\text{Re}=u_{0}L/\nu$ as $\eta_{h}\sim L\text{Re}^{-1/(1+h)}$. This scale thus appears as a fluctuating cut-off which depends on the scaling exponent and therefore on $\bm{x}$. This is the generalization of the Kolmogorov scale $\eta\sim\nu^{3/4}\equiv\eta_{\frac{1}{3}}$, and was first proposed in [P87]. Below $\eta_{h}$, the velocity field becomes laminar, and $|\delta_{\ell}\textbf{u}|\propto\ell$. When the velocity field is turbulent, $h\equiv\log(|\delta_{\ell}\textbf{u}|/u_{0})/\log(\ell/L)$ varies stochastically as a function of space and time. Also, if the turbulence is statistically homogeneous, stationary and isotropic, $h$ only depends on $\ell$, the scale magnitude. Therefore, formally, $h$ can be regarded as a continuous stochastic process labeled by $\log(\ell/L)$. By Kramer’s theorem (Touchette, 2009), one sees that as in the limit $\ell\to 0$, $\log(L/\ell)\to\infty$, we have

[TABLE]

where $C(h)$ is the rate function of $h$, also called multifractal spectrum. Formally, $C(h)$ can be interpreted as the co-dimension of the set where the local Hölder exponent at scale $\ell$ is equal to $h$. Using Gärtner-Elis theorem (Touchette, 2009), one can connect $C$ and the velocity structure functions as:

[TABLE]

To proceed further and make connection with previous section, we set $\epsilon=u_{0}^{3}/L$ so that $S_{p}(\ell)$ can now be written:

[TABLE]

This shows that $\tau(p)$ is the Legendre transform of the rate function $C(h+1/3)$, i.e. $\tau(p)=\min_{h}(p(h-1/3)+C(h))$, and equivalently, that $C(h)$ is the Legendre transform of $\tau(p)$. Because of this, it is necessarily convex. The set of points where it satisfies $C(h)\leq d$, represents the set of admissible or observable $h$, is therefore necessarily an interval, bounded by $-1\leq h_{\min}$ and $h_{\max}\leq 1$.

As noted by Frisch and Vergassola (1991), the scaling exponent $\zeta(p)=p/3+\tau(p)$ defined via Equation (11) is only constant in a range of scale where $\ell>\eta_{h}$ for any $h\in[h_{\min},h_{\max}]$. For small enough $\ell$, this condition is not met anymore, since as soon as $\ell<\eta_{h}$, all velocity fields corresponding to $h$ are ”regularized”, and do not contribute anymore to intermittency since they scale like $\ell$. This results in a slow dependence of $\zeta(p)$ with respect to the scale, which is obtained via the corrected formula:

[TABLE]

To understand the nature of the correction, we can compute the value of $h$ such that $\ell=\eta(h)$. It is simply: $h(\ell)=-1+\log(\text{Re})/\log(L/\eta)$. We note $\theta=\log(L/\ell)/\log(\text{Re})$. We can now rewrite equation (12) as:

[TABLE]

where $\tau(p,\theta)=\tau_{p}$ when $\theta\leq 1/(1+h_{\max})$ and $\tau(p,\theta)=p(\theta-1/3)+C(-1+1/\theta)$ when $1/(1+h_{\max})\leq\theta\leq 1/(1+h_{\min})$. As discussed by Frisch and Vergassola (1991), this implies a new form of universality that extends beyond the inertial range, into the so-called extended dissipative range, as;

[TABLE]

If the scale $L$ is constant and equal to $L_{0}$, the integral scale, then we have $\text{Re}=(L_{0}/\eta)^{4/3}$ and the multifractal universality implies that $\log(\tilde{S}_{p})/\log(L_{0}/\eta)$ is a function of $\log(\ell/\eta)/\log(L_{0}/\eta)$. When the function is linear, we thus recover the K62 universality. The multifractal universaility is thus a generalization of the K62 universality.

This form of universality is however not easy to test, as the scale $L$ is not known a priori, and may still depend on Re. In what follows, we demonstrate a new form of universality, that allows more freedom upon $L$ and encompass both K62 and multifractal universality.

III.4 General universality

Using the hypothesis that turbulence maximizes some energy transfer in the scale space, Castaing Castaing, B. (1989) suggested a new form of universality for the structure functions, that reads:

[TABLE]

where $A_{p}$ and $K_{0}$ are universal constant and $\beta$ and $G$ are general functions, $F$ being linear in the inertial range, $G(p,x)\sim\tau(p)x$. The validity of this universal scaling was checked by Gagne and Castaing Castaing et al. (1993) on data obtained from the velocity fields measured in a jet using hot wire anenometry. They found good collapse of the structure functions at different Taylor Reynolds $R_{\lambda}$, provided $\gamma(\text{Re})$ is constant at low Reynolds numbers and follows a law of the type: $\gamma(\text{Re})\sim\gamma_{0}/\log(R_{\lambda}/R_{*})$, where $R_{*}$ is a constant, whenever $R_{\lambda}>400$. Since we have $R_{\lambda}\sim\text{Re}^{1/2}$ and $(L_{0}/\eta)\sim\text{Re}^{3/4}$, we can rewrite equation (15) as:

[TABLE]

where $S_{0p}$ are some constants and $\beta$ and $F$ are general functions. Comparing with the K62 or MFR universality formulae (7) or (14), we see that formula (16) is a generalization of these two universality with $L=L_{0}$. It allows however more flexibility than K62 or MFR universality through the function $\beta(\text{Re})$, that is a new fitting function. We test these predictions in Section IV and provide a physical interpretation of (16) in Section V.

IV Check of universality using data analysis

The various universality are tested using the velocity structure functions based on the wavelet velocity increments Eq. (2), in order to minimize the noise in the experimental data. We define:

[TABLE]

We then apply this formula to both experimental data (Table 1) and numerical data (Table 2), to get wavelet velocity structure functions at various scales and Reynolds numbers.

IV.1 Check of K41 universality

The K41 universality (3) can be checked by plotting:

[TABLE]

This is shown in figure 2 for both experimental and numerical data. Obviously, the data do not collapse on a universal curve, meaning that K41 universality does not hold. This is well known, and is connected to intermittency effects Frisch (1996).

IV.2 Check of K62 universality

The K62 universality (7) can be checked by plotting:

[TABLE]

The collapse depends directly on $\tau(p)$, the intermittency exponents. Obtaining the best collapse of all curves is in fact a way to fit the best scaling exponents $\tau(p)$. We thus implemented a minimization algorithm that provides the values of $\tau(p)$ that minimized the distance between the curve and the line of slope $\tau(p)$. The values of $\tau(p)$ are reported in Table 3. The best collapse is shown on Figure 3-a for the DNS, and Figure 3-b for the experiment. The collapse is better for experiments than for the DNS. However, in both cases, there are significant differences in between points at different $R_{\lambda}$, at larger scales, showing that universality is not yet reached.

IV.3 Check of General Universality

We can now check the most general universality, by plotting:

[TABLE]

In this case, best collapse is obtained by fitting two families of parameters: $S_{0p}$, $\beta(\text{Re})$ that were obtained through a procedure of minimization. We take the DNS at $R_{\lambda}=138$ as the reference case, and find for both DNS and experiments, the values of $\beta(\text{Re})$ and $S_{0p}$ that best collapse the curves. The corresponding collapses are provided in figure 5. The collapse is good for any value of Re, except for the DNS at the lowest Reynolds number, which does not collapse in the far dissipative range.

IV.4 Function $\beta(\text{Re})$

Motivated by earlier findings by Castaing et al. (1993), we plot in figure 6 the value $1/\beta$ as a function of $R_{\lambda}$.

Our results are compatible with $1/\beta\sim\beta_{0}/\log(R_{\lambda})$, with $\beta_{0}\sim 4/3$ over the whole range of Reynolds number. For comparison, we provide also on figure 6 the values found by Gagne and Castaing Castaing et al. (1993) in jet of liquid Helium, shifted by an arbitrary factor to make our values coincide with them at large Reynolds number. This shift is motivated by the fact that $\beta(\text{Re})$ is determined up to a constant, depending upon the amplitude of the structure functions used as reference. At large Reynolds, our values are compatible with theirs. At low Reynolds, however, we do not observe the saturation of $1/\beta$ that is observed in the jet experiment of Castaing et al. (1993). An interpretation of the meaning of $\beta(\text{Re})$ will be provided in Section V.

IV.5 Scaling exponents

Our Collapse method enables us to obtain the scaling exponents of the structure functions $\zeta(p)$ by the following two methods:

i) Using the K62 universality, we get $\tau(p)$, and then $\zeta(p)=\zeta(3)p/3+\tau(p)$. These estimates still depend on the value of $\zeta(3)$, which is not provided by the K62 universality plot. To obtain it, we use a minimization procedure on both experimental $\log(S_{3}/u_{\mathrm{K}}^{3})$ from the one hand, and the numerical $\log(S_{3}/u_{\mathrm{K}}^{3})$ on the other hand (see figure 4-a), to compute $\zeta(3)$ as the value that minimizes the distance between the curve and a straight line of slope $\zeta(3)$. The values so obtained are reported in Table 3, and have been used to compute $\zeta(p)$ from $\tau(p)$.

ii) Using the general universality, we may also get $\tau_{p,\text{univ}}$ by a linear regression on the collapse curve. Note that since the data are collapsed, this provides a very good estimates of this quantity, with the lowest possible noise. In practice, we observe no significant differences with the two estimates; therefore, we only report the values obtained by following the first method.

The corresponding values are plotted in figure 4 and summarized in Table 3. Note that for both DNS and experiments, the value of $\zeta(3)$ is different from $1$, which is apparently incompatible with the famous Kolmogorov $4/5$th law, that predicts $\zeta(3)=1$. This is because we use absolute values of wavelet increments, while the Kolmogorov $4/5$th law uses signed values. We have checked that using unsigned values, we obtain a scaling that is closer to $1$, but with larger noise. Note also that when we consider the relative value $\zeta(p)/\zeta(3)$, we obtain values that are close to the values obtained Saw et al. (2018) on the same set of experimental data, using velocity increments and Extended Self-Similarity technique Benzi et al. (1993).

IV.6 Multifractal spectrum

From the values of $\zeta(p)$, one can get the multifractal spectrum $C(h)$ by performing the inverse Legendre transform:

[TABLE]

Practically, this amount to use the following formula:

[TABLE]

To estimate $C$, we thus first perform a polynomial interpolation of order 4 on $\zeta(p)$, then derivate the polynom to estimate $\frac{\mathrm{d}\,\zeta(p)}{\mathrm{d}\,p}$, thus get $C$ through equation (22). The result is provided in figure 6-b for both the DNS and the experiment.

The curve look like the portion of a parabola, corresponding to a log-normal statistics for the wavelet velocity increments. Specifically, fitting by the shape:

[TABLE]

we get $a=0.35$ and $b=0.045$. This parabolic fit also provides a good fit of the scaling exponents, as shown in figure 4 by performing Legendre transform of $C(h)$ given by equation (23).

V Thermodynamics and turbulence

V.1 Thermodynamical analogy

Multifractals obey a well-known thermodynamical analogy Bohr and Rand (1987); [M91]; Rinaldo et al. (1996) that will be useful to interpret and extend the general universality unraveled in the previous section. Indeed, consider the quantity:

[TABLE]

By definition $\mu_{\ell}$ is positive definite and $\langle\mu_{\ell}\rangle=1$ for any $\ell$. It is therefore the meaning of a scale dependent measure. It then also follows a large-deviation property as:

[TABLE]

where $S(E)$ is the large deviation function of $\log(\mu_{\ell})$ and has the meaning of an energy while $\log(\ell/\eta)$ has the meaning of a volume, and $\log(\mu_{\ell})/\log(\ell/\eta)$ is an energy density. Because of the definition of $\mu_{\ell}$, it is easy to see that $S$ is connected to $C$, the large deviation function of $|\delta W_{\ell}$. In fact, since in the inertial range where $\langle|\delta W_{\ell}|^{3}\rangle\sim\ell^{\zeta(3)}$, we have $S(E)=C(3h-\zeta(3))$. By definition, we also have:

[TABLE]

so that $Z=\tilde{S}_{3p}$ is the partition function associated to the variable $\log(\mu_{\ell})$, at the pseudo-inverse temperature $p=1/kT$. Taking the logarithm of $Z$, we then get the free energy $F$ as:

[TABLE]

By the Gärtner-Elis theorem, $F$ is the Legendre transform of the energy $S$: $F=\min_{E}(pE-S(E))$. The free energy a priori depends on the temperature; i.e. on $T=1/kp$, on the volume $V=\log(\ell/\eta)$ and on the number of degrees of freedom system $N$. If we identify $N=(1/\beta(\text{Re}))\log(L_{0}/\eta)$, we see that the general universality means:

[TABLE]

i.e. can be interpreted as extensivity of the free energy.

The thermodynamic analogy is thus meaningful and is summarized in Table 4. It can be used to derive interesting prospects.

V.2 Multifractal pressure and phase transition

Given our free energy, $F=\log(\tilde{S}_{3p})$, we can also compute the quantity conjugate to the volume, i.e. the multifractal pressure as: $P=\partial F/\partial V$. In the inertial range, where $\tilde{S}_{p}\sim\ell^{\tau(p)}$, we thus get $P=\tau(p)$, which only depends on the temperature. Outside the inertial range, $P$ has the meaning of a local scaling exponents also depends upon the scale, i.e., on the volume $V$ and on $N$ (Reynolds number). Using our universal functions derived in figure 5, we can then compute empirically the multifractal pressure $P$ and see how it varies as a function of $T$, $V$ and $N$. It is provided in figure 7 for $R_{\lambda}=25$ and $R_{\lambda}=56$, and in figure 8 for $R_{\lambda}=90$ and $R_{\lambda}=138$. We see that at low Reynolds number, the pressure decreases monotonically from the dissipative range, reaches a lowest points and then increases towards the largest scale. There is no clear flat plateau that would correspond to an ”inertial” range. In contrast, at higher Reynolds number, a plateau appears for $p=1$ to $p=4$ when going towards the largest scale, the value of the plateau corresponding to $\tau_{\text{DNS}}$. The plateau transforms into an inflection point for $p\geq 5$ making the derivative $\partial P/\partial V$ change sign. This is reminiscent of a phase transition occurring in the inertial range, with coexistence of two phases: one ”laminar” and one ”turbulent”. We interpret such a phase transition as the result of the coexistence of region of flows with different Hölder exponents, with areas where the flow has been relaminarized due to the action of viscosity, because of the random character of the dissipative scale (see below).

VI Conclusion

We have shown that a deep analogy exists between multifractal and classical thermodynamics. In this framework, one can derive from the usual velocity structure function an effective free energy that respects the classical extensivity properties, provided one uses a number of degrees of freedom (given by $N=1/\beta(\text{Re}))$ that scales like $\log(R_{\lambda}$). This number is much smaller than the classical $N\sim\text{Re}^{9/4}$ that is associated with the number of nodes needed to discretize the Navier-Stokes equation down to the Kolmogorov scale. It would be interesting to see whether this number is also associated with the dimension of a suitable ”attractor of turbulence”. Using the analogy, we also found the ”multifractal” equation of state of turbulence, by computing the multifractal pressure $P=\partial F/\partial V$. We found that for large enough $R_{\lambda}$ and $p$ (the temperature), the system obeys a phase transition, with coexistence of phase like in the vapor-liquid transition. We interpret this phase transition as the result of the coexistence of region of flows with different Hölder exponents, with areas where the flow has been relaminarized due to the action of viscosity, because of the random character of the dissipative scale. We note that this kind of phenomenon has already been observed in the context of Lagrangian velocity increments, using the local scaling exponent $\zeta(p,\tau)=\mathrm{d}(\log(S_{p}(\tau)))/\mathrm{d}(\log(\tau))$ Arneodo et al. (2008). The phase transition is then associated with the existence of a fluctuating dissipative time scale. It has further been shown that in a multifractal without fluctuating dissipative time scale, the local exponent decreases monotically from dissipative scale to large scale, implying a disappearance of the phase transition Biferale et al. (2004).