## Strong anisotropy of superfluid 4 He counterflow turbulence

L. Biferale 1 , D. Khomenko , V. L'vov , A. Pomyalov , I. Procaccia 2 3 3 3 and G. Sahoo 4

1 Dept. of Physics, University of Rome, Tor Vergata, Roma, Italy 2 Laboratoire de Physique , D´ epartement de physique de l'ENS, ´ cole E Normale Sup´ erieure, PSL Research University, Sorbonne Universit´s, e CNRS, 75005 Paris, France Dept. of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot, Israel

3

4 Dept. of Mathematics and Statistics and Dept. of Physics, University of Helsinki, Finland

We report on a combined theoretical and numerical study of counterflow turbulence in superfluid 4 He in a wide range of parameters. The energy spectra of the velocity fluctuations of both the normal-fluid and superfluid components are strongly anisotropic. The angular dependence of the correlation between velocity fluctuations of the two components plays the key role. A selective energy dissipation intensifies as scales decrease, with the streamwise velocity fluctuations becoming dominant. Most of the flow energy is concentrated in a wavevector plane which is orthogonal to the direction of the counterflow. The phenomenon becomes more prominent at higher temperatures as the coupling between the components depends on the temperature and the direction with respect to the counterflow velocity.

## Introduction

the flow energy is concentrated in the wavevector plane orthogonal to the counterflow.

Below a critical temperature T λ ≈ 2 17 K liquid . 4 He behaves as a quantum fluid 1,2 , consisting of an inviscid superfluid, associated to the quantum ground state, and a gas of thermal excitations which make up the viscous normal fluid. Quantum mechanics 4 constrains the rotational motion of the superfluid in 4 He to discrete ˚ngstr¨m-width quantum vortex lines of fixed circulaA o tion. The thermal excitations scatter on a dense tangle of these vortices, thus inducing a mutual friction force between the normal fluid and the superfluid.

Turbulent superfluid Helium in a channel with a temperature gradient is a subject of extensive research for many decades 1-10 . In such a setting, so-called 'counterflow', the normal fluid flows from the hot end of the channel to the cold end while the superfluid flows in the opposite direction. Most attention was devoted so far to the measurement and the analysis of the density of vortex lines and to the mutual friction between the components.

Recent advances in the visualization techniques offer for the first time a direct access to the statistics of the velocity fluctuations of the normal fluid 11-13 and the superfluid 14-16 . It was shown that the large scale statistics of the normal fluid in the counterflow is very different 11-13 from the statistics of classical fluids. The theoretical analysis 17-19 highlighted the importance of correlations between the superfluid and the normal fluid components, which lead to the energy spectra of both components being steeper than their classical counterparts. Moreover, we have recently shown 19 that the direction of the mean relative velocity plays an important role; The correlation between the 4 He components decays slower for eddies stretched along the counterflow velocity. In contrast, the correlation of eddies, which are elongated in the orthogonal direction, decay faster, leading to their enhanced energy loss. As a result of this directionally preferred energy dissipation, the velocity fluctuation consist mostly of the streamwise components, while most of

Here we consider this phenomenon further and study its consequences in further detail. The paper is organized as follows: in the Sec. I we provide a sketch of a theory of counterflow turbulence with a stress on its anisotropy. In Sec. I A we introduce the basic set of coarse-grained equations for the counterflow. These are used for the theoretical analysis and for the numerical simulations. In Sec. I B we clarify how various approaches to the statistical description of the anisotropy energy surface are related. Next, in Sec. I C, we discuss the physical origin of the strong spectral anisotropy in counterflow turbulence. In Sec. II we present the results of the direct numerical simulations (DNS) of the two-fluid coarse-grained Eqs. (1a). The main conclusion is that the analytical predictions are confirmed. In the first subsection Sec. II A we discuss the simulation parameters and the numerical procedure. Next, in Sec. II B, we use standard statistical characteristics: one-dimensional (1D) energy spectra and cross-correlation functions, averaged over a spherical surface of radius k (i.e. over all directions of vector k ) to provide an overview of spectral properties of 4 He counterflow. We find that at the smallk regime, the normal-fluid and superfluid velocity components are indeed well correlated. As expected, mutual friction plays secondary role and the spherically-averaged spectra are similar to the spectra in the 4 He-coflow turbulence 20 , being only slightly steeper than the Kolmogorov-1941 (K41) spectra of classical hydrodynamic turbulence. On the other hand, at relatively large k the fluid components are practically uncorrelated; mutual friction provides a leading contribution to the energy dissipation and the counterflow spectra are similar to those in 3 He superfluid turbulence with the normal-fluid component at rest 21,22 . The spectra become strongly suppressed in comparison to K41 energy spectra.

The similarities between the inherently anisotropic counterflow energy spectra and the isotropic spectra in

the turbulent 4 He-coflow and 3 He are, however, superficial. To expose the differences, we discuss in Sec. II C the two-dimensional (2D) energy spectra which depend, besides the wavenumber k , upon the angle θ between the wave-vector k and the counterflow velocity U ns . Here we find that the spectra become more and more anisotropic with increasing k , being confined in k -space to a small range cos θ &lt; 0 1, i.e. . near the wavevectors plane which is orthogonal to U ns . This effect becomes stronger with increasing temperature. The tensor structure of the energy spectra, considered in Section II D, is also temperature dependent: the small-scale turbulent velocity fluctuations are dominated by only one vector component, parallel to U ns , becoming more so at higher temperature. Further, we compare several variants of differently averaged 1D spectra ( Sec. II E) and structure functions (Sec. II F) to expose other aspects of the spectral anisotropy in connection with possible experimental observations.

In Sec. III we summarize our findings: Counterflows exhibit strongly anisotropic energy distributions. The energy spectra are localized near a direction that is orthogonal to the counterflow. The phenomenon is similar to atmospheric turbulence with a strong stable stratification or to rotational turbulence 23-28 . On the other hand, the tensor structures of these two types of quasi2D turbulence are quite the opposite: in atmospheric turbulence the vertical component of the turbulent velocity fluctuations is suppressed by the stratification and only the two horizontal components are dominant 23 . In the counterflow turbulence the main contribution to the turbulent energy comes from one streamwise velocity projection, while the two cross-steam velocity projections are strongly suppressed. The observed phenomenon is mild at low temperatures and becomes more prominent as the temperature increases. We confirm, in agreement with Ref. 13 , that the structure functions of the turbulent velocities in the counterflow do not reflect in a quantitative manner the underlying energy spectra. However, the relative magnitude of the structure functions, measured in different directions, may qualitatively reflect the presence of the spectral anisotropy.

## I. QUALITATIVE ANALYSIS OF ANISOTROPIC COUNTERFLOW TURBULENCE

As we mentioned in the introduction, one of important properties of superfluid 4 He is the quantization of vorticity, which concentrates on vortex-lines of core radius a 0 ≈ 10 -8 cm with fixed circulation κ = h/M ≈ 10 -3 cm /s. 2 Here h is Planck's constant and M is the mass of the 4 He atom . 4 A complex tangle of these vortex lines with a typical inter-vortex distance 3 /lscript ∼ 10 -4 -10 -2 cm is a manifestation of superfluid turbulence 4 .

On the large scales this type of turbulence is commonly described by the two-fluid model. The density of 4 He ρ is modelled as a mixture of two fluid components: an inviscid superfluid and a viscous normal fluid, with respective densities ρ s and ρ n such that ρ = ρ s + ρ n . The fluid components are coupled by a mutual friction force 1,3,5,7-9 .

Large-scale turbulence in 4 He can be generated by various ways. In mechanically driven 4 He (so-called 'coflow'), the turbulent statistics is similar 10,20,30-32 to that of classical turbulence. In this case both components move in the same direction and the mutual friction force couples them almost at all scales. On the other hand, when a temperature gradient ∇ T is imposed in a channel closed at one end, the heat flux is carried away by the normal fluid with a mean velocity U n ∝ ∇ T , while the superfluid component flows in the opposite direction with the mean velocity U s . There is no net mass flow: ρ n U n + ρ s U s = 0. The counterflow velocity U ns = U n -U s creates a random vortex tangle with an energy spectrum, peaking at the intervortex scale /lscript and with a close to Gaussian statistics, as demonstrated experimentally in Refs. 14-16 and rationalized theoretically in Refs. 18,33 . At large enough U n , the laminar flow of the normal component become unstable, creating large scale turbulence with the energy spectrum dominated by 'the outer scale of turbulence' ∆ /greatermuch /lscript (e.g. about half-width of the channel). Although the particular mechanisms of the large scale superfluid motion generation are not known in details, recent indirect experimental evidence indicate 13,33 that the large-scale normal fluid motion gives rise to the superfluid turbulent motion due to the components' coupling by the mutual friction force.

## A. Coarse-grained equations for counterflow 4 He turbulence

Our approach 17,18,34 to large-scale counterflow turbulence is based on two Navier-Stokes equations (NSE) for the velocity fluctuations of the normal fluid and superfluid components U n ( r , t ) and U s ( r , t ). A complication arises from the fact that the counterflow is created in a channel and therefore is, in general, inhomogeneous. However, at large enough Reynolds numbers, the flow in the center of the channel can be approximated as almost space-homogeneous 35 . We therefore adopt a simplifying description with space homogeneity and stationarity. Further, we perform the standard Reynolds decomposition of the velocities into their mean and turbulent velocity fluctuations with zero mean:

$$\mathcal { U } _ { n } ( r, t ) \ & = \ U _ { n } + u _ { n } ( r, t ) \,, \ U _ { n } = \langle \mathcal { U } _ { n } ( r, t ) \rangle \, \ \ ( \text{la} ) \\ \mathcal { U } _ { s } ( r, t ) \ & = \ U _ { s } + u _ { s } ( r, t ) \,, \ U _ { s } = \langle \mathcal { U } _ { s } ( r, t ) \rangle \.$$

The mean velocities are taken below as externally prescribed parameters of the problem. Note that in the classical hydrodynamics, the Navier-Stokes equations are Galilean invariant and one can choose a reference system in which the constant mean velocity vanishes. In the twofluid counterflow, there is no such reference system and the mean velocities are necessarily present in the equa-

tions of motion for the turbulent velocity fluctuations:

$$[ \frac { \partial } { \partial t } + ( u _ { s } + U _ { s } ) \cdot \nabla ] u _ { s } - \frac { \nabla p _ { s } } { \rho _ { s } } \ = \nu _ { s } \, \Delta u _ { s } + f _ { n s } + \varphi _ { s } \,, ( 1 b ) \quad & \text{fuctua} \\ [ \frac { \partial } { \partial t } + ( u _ { n } + U _ { n } ) \cdot \nabla ] u _ { n } - \frac { \nabla p _ { n } } { \rho _ { n } } = \nu _ { n } \, \Delta u _ { n } - \frac { \rho _ { s } } { \rho _ { n } } f _ { n s } + \varphi _ { n } \,, \quad & \text{gested} \\ \text{form } f \\ \Omega _ { s } \, \text{was}$$

These equations are coupled by the mutual friction force 2,6,36 f ns . Here f ns ( r , t ) is a fluctuating (with zero mean) part of the total mutual friction force F ns ( r , t ):

$$\mathcal { F } _ { n s } ( r, t ) = F _ { n s } + f _ { n s } ( r, t ) \. \quad \ \ ( 1 c ) \ \text{me}$$

The pressures p n , p s in Eqs. (1b) are given by

$$p _ { n } & = \frac { \rho _ { n } } { \rho } \left [ p + \frac { \rho _ { s } u _ { n s } ^ { 2 } } { 2 } \right ], \, p _ { s } = \frac { \rho _ { s } } { \rho } \left [ p - \frac { \rho _ { n } u _ { n s } ^ { 2 } } { 2 } \right ], \, u _ { n s } = u _ { n } - u _ { s } \,. \quad \text{$\text{count}(\text{$\text{on exit}$} \\ ( 1 d ) \quad \text{$\text{$\text{used}$} \, of the normal $f$} \,f$} \,f$} \,.$$

The kinematic viscosity of the normal fluid component is ν n = η/ρ n , where η is the dynamical viscosity 29 of the normal 4 He component. The energy sink 34 in the equation for the superfluid component, with an effective superfluid viscosity 3 ν s , accounts for the energy dissipation at the intervortex scale /lscript due to vortex reconnections and energy transfer to Kelvin waves. The random forces ϕ s and ϕ n represent the forcing of the turbulent flow at large scales.

The physical origin of the mutual friction is the scattering of excitations that constitute the normal fluid on the vortex lines. Any motion of a vortex line relative to the normal fluid results 3 in a force per unit length of the line, which can be written as

$$& \overset { \dots } { \dots }, \, \dots \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overs ^ { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots} \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \oversets { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots  \, \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots  \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overs^{ \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \, \overset { \dots } { \dots } \\ & \text{Here} \, \, \text{$s^{\prime}$ is a unit vector along the length of the vertex $\text{ $\gamma_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \text{$i_{\alpha}$,} \quad \end{$$

Here s ′ is a unit vector along the length of the vortex, γ 0 and γ ′ 0 are some phenomenological parameters.

In order to estimate the coarse-grained mutual friction force f ns in Eqs. (1b) one needs to properly average the microscopic Eq.(2) for f . The result depends on the statistics of the quantum vortices that in turn depend on the particular turbulent flow properties. In particular, this procedure includes averaging of the force f over directions of the orientations s ′ in the vortex lines. In the relatively simple case of rotating turbulence, the vortices are oriented mostly along the axis of rotation. In this case 〈 s ′ 〉 may be directly related to the direction of the superfluid vorticity ω s : 〈 s ′ 〉 = ω s / ω | s | . The resulting equations were named the 'HVBK equations' 5,6 after Hall, Vinen, Bekarevich and Khalatnikov.

Clearly, the original HVBK equations are not applicable for the superfluid turbulence without global rotation, for which 〈 s ′ 〉 = 0. In this case, to obtain a coarsegrained representation one should average Eq.(2) over 'a physically small volume' of scale δ . This scale should be chosen to be much larger than /lscript , but still much smaller than the scale r of turbulent fluctuations under consideration, /lscript /lessmuch δ /lessmuch r . For such δ , the local line orientations s ′ and s ′ ⊗ s ′ in Eq. (2) can be considered as self-averaging in space as they are almost uncorrelated with the r -scale fluctuations u ns ( r ), which are treated as dynamical variables.

For co-flows a number of model expressions were suggested for the fluctuating part of the friction force in the form f ns ∝ Ω s u ns . Here the mutual friction frequency Ω was modeled as a dimensional estimate, assuming uns derlying Kolmogorov energy spectrum for the superfluid component. Examples of such models include Ω = s α ω | s | (e.g. in Refs. 32,37,38 ), Ω s α √ 〈 ω 2 s 〉 (e.g. in Refs. 17,20,39 ) and Ω s = ακ L (e.g. in Ref. 2,13,18,20 ). Here α is the dimensionless mutual friction parameter related 1 to γ 0 as αρ κ s = (1 + α 2 ) γ 0 , L is the vortex line density.

In counterflows the dynamics of the vortex tangle is dominated by the stretching of the vortex lines by the counterflow velocity and by their reconnections. Based on experiments in narrow slits, Gorter and Mellink 40 proposed to couple the equations of motion for the components' velocities by the mutual friction force of a phenomenological form F ns = Aρ ρ n s U 3 ns , where A is a temperature dependent constant. This form was later refined by Vinen 7 for homogeneous turbulence and an isotropic vortex tangle as F ns = G U ns , G = Aρ ρ n s U 2 ns . Taking into account the relation between the vortex line density and the counterflow velocity in the steady-state isotropic tangle, it can be further rewritten 2 as F ns = 2 3 ακ L U ns .

The tangle anisotropy with respect to the direction of U ns can be described by the Schwarz's indices 9 :

$$& \mathcal { U } _ { n } \text{ can be described by one suchwarz$$

Here L tot is the total vortex length in the whole vortex configuration C over which integrals are taken and ˆ r ‖ and ˆ r ⊥ are unit vectors in the directions parallel and perpendicular to U ns respectively. Using (3a) and (3b), the mutual friction force may be written as

$$F _ { n s } \, = \, \alpha \kappa \mathcal { L }$$

$$r _ { n s } \, & = \, \alpha \kappa \mathcal { L$$

Notably, the second non-dissipative term in Eq. (2) ∝ γ ′ 0 , vanishes by symmetry and does not contribute 9 to the averaged quantity F ns . The mean mutual friction force F ns , (3c), found earlier by Schwarz 9 , enters into the equations for the mean velocities U n and U s , which we do not discuss. The fluctuating part (3d) of the mutual friction force f ns ( r , t ) enters Eqs. (1b). The vector components of turbulent counterflow velocity fluctuations u ‖ ns and u ⊥ ns are oriented in the directions of ˆ r ‖ and ˆ r ⊥ respectively.

The definitions Eq. (3) do not take into account that the turbulent intensity of the normal fluid velocity √ 〈| u n | 2 〉 /U ns in the counterflow turbulence is not very small and can reach 12 values of about 0 25. . So, strictly speaking, in Eq. (3a) and (3b) the directions ˆ r ‖ and ˆ r ⊥ ,

should be taken along and orthogonal to the total counterflow velocity U ns . In addition, one should take into account the space-time dependence of L in Eq. (3d). However, since the leading contribution to √ 〈| u ns | 2 〉 originates from the outer scale of turbulence ∆ /greatermuch /lscript , for the motions of the scale /lscript this correction effectively adds to U ns and we can neglect the influence of the velocity fluctuations of the scale r on the vortex line density, replace an average of the product by the product of averages and consider L as constant. The equation (3d) may be identically rewritten as

$$\begin{array} {$$

Taking the numerical values of I ‖ and I ⊥ (see e.g. Tab. IV for T = 1 6 K in Ref. . 41 ) we see that ( I ⊥ -I ‖ ) / 2 ≈ 0 05. . This means that with a reasonable accuracy of about (5 -10)% we can neglect the anisotropy term in Eq. (4) and use the simple form f ns ( r , t ) = Ω s u ns ( r , t ) as a good approximation for f ns even for counterflow turbulence. Furthermore, in this paper we do not consider the dependence of L on the flow parameters and use the mutual friction frequency Ω s as a prescribed external control parameter. It can be estimated or measured for each particular flow conditions.

## B. Statistical characteristics of anisotropic turbulence

## 1. Velocity correlation function

A useful characterization of homogeneous superfluid 4 He turbulence is furnished by the three-dimensional (3D) correlation functions of the normal- and superfluid turbulent velocity fluctuations in k -representation:

$$\text{current velocity} \text{ communications in } \kappa \text{-representation}$$

Here δ ( k -k ′ ) is a 3D Dirac's delta function, v j ( k ) is the Fourier transform of u j ( r )

$$v _ { j } ( k ) \, = \, \int _ { \Delta _ { \Delta } } u _ { j } ( r ) \, \exp ( - i k \cdot r ) \, d r \,, \quad \ \ ( 5 c ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ =$$

$$u _ { j } ( r ) \, = \, \int v _ { j } ( k ) \, \exp ( - i k \cdot r ) \, \frac { d k } { ( 2 \pi ) ^ { 3 } } \,, \quad ( 5 \text{d} ) \quad \sqrt { k _ { \| } ^ { 2 } } _ { \text{relat} }$$

the indices α, β = { x, y, z } denote Cartesian coordinates, the subscripts ' i,j ' denote the normal ( i, j =n) or the superfuid ( i, j =s) fluid components and ∗ stands for complex conjugation. In the rest of the paper, we denote the trace of any tensor according to E ij ( k ) = ∑ α E αα ij ( k ) .

The correlation function E αβ ij ( k ) and the Fourier transform(5c) are defined such that the kinetic energy density per unite mass E j reads

$$\kappa _ { 1 } \, \quad \cdots \, \quad \cdots \, \quad \cdots \, \quad \cdots \\ E _ { j } = \frac { 1 } { 2 } \, \langle | u _ { j } ( r ) | ^ { 2 } \rangle = \frac { 1 } { 2 } \int \mathcal { E } _ { j j } ( k ) \frac { d ^ { 3 } k } { ( 2 \pi ) ^ { 3 } } \. \quad ( 5 e ) \\ \tau _ { 1 } \, \quad \cdots \, \quad \cdots \, \quad \cdots \, \quad \cdots \, \quad \cdots \,.$$

The dimensionality of the energy density is [ E jj ] =cm /s , 2 2 while the dimensionality of the 3D energy spectra is [ E jj ] =cm /s . 5 2

Due to the presence of a preferred direction (the counterflow velocity), the resulting turbulence has an axial symmetry around that direction. Accordingly, E αβ ij ( k ) depends only on two projections k ‖ and k ⊥ of the wavevector k : k ‖ ≡ U ns ( k · U ns ) /U 2 ns and k ⊥ ⊥ U ns , being independent of the angle ϕ in the ⊥ -plane, orthogonal to U ns : E αβ ij ( k ) ⇒ E αβ ij ( k , k ‖ ⊥ ). This allows us to define a two-dimensional (2D) object E αβ ij ( k , k ‖ ⊥ ) that still contains all the information about 2 nd -order statistics of the counterflow turbulence:

$$E _ { i j } ^ { \alpha \beta } ( k _ { \| }, k _ { \perp } ) \equiv \frac { k _ { \perp } } { 4 \pi ^ { 2 } } \mathcal { E } _ { i j } ^ { \alpha \beta } ( k _ { \| }, k _ { \perp } ) \. \quad \quad ( 6 \text{a} )$$

Another way to represent the same information is to introduce a polar angle θ = ∠ ( k U , ns ), to represent the wavevector length as k = √ k 2 ‖ + k 2 ⊥ and to define a 2D object ˜ E αβ ij ( k, θ ) in spherical coordinates:

$$\tilde { E } _ { i j } ^ { \alpha \beta } ( k, \theta ) \equiv \frac { k \mathcal { E } _ { i j } ^ { \alpha \beta } ( k ) } { 4 \pi ^ { 2 } } = \frac { k } { 4 \pi ^ { 2 } } \mathcal { E } _ { i j } ^ { \alpha \beta } ( k \cos \theta, k \sin \theta ). \ ( 6 b )$$

The dimensionality of 2D energy spectra and correlation functions is [ E αβ ij ] = [ ˜ E αβ ij ]=cm /s . 4 2

A more compact but less detailed information on the statistics of turbulence is provided by a set of onedimensional (1D) energy spectra. The most traditional are the 1D ' spherical ' energy spectra and the crosscorrelation function · E αβ ij ( k ), averaged over a spherical surface of radius k :

$$& \text{values} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \,  \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text` \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \,  \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \,\infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{array} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \. \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \, \infty \, \text{raw} \,$$

In the isotropic case, when E αβ ij depends only on k = √ k 2 ‖ + k 2 ⊥ , this representation simplifies to a well know relationship

$$\text{ites}, \quad \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text$$

Further information about the anisotropy of the 2 nd -order statistics is obtained by comparing the spherical

1D-spectra · E ij ( k ) with a set of 1D spectra averaged differently. A natural choice are spectra averaged over a cylinder of radius k ⊥ with the axis oriented along k ‖ . This results in the cylindrical 1D spectra

$$\circ E _ { i j } ^ { \alpha \beta } ( k _ { \perp } ) = \int E _ { i j } ^ { \alpha \beta } ( k _ { \| }, k _ { \perp } ) d k _ { \| }. \quad \quad ( 7 c ) \quad \underset { S _ { n } ^ { x x } | } { \text{line} }$$

Alternatively, one can average the 3D function E ij ( k , k ‖ ⊥ ) over a plane. Here we choose two planes:

- 1D-spectra ⊥ E αβ ij ( k ‖ ), averaged over a ⊥ -plane , orthogonal to U ns

$$\Phi _ { \quad } \Phi _ { i j } ^ { \alpha \beta } ( k _ { \| } ) = \int E _ { i j } ^ { \alpha \beta } ( k _ { \| }, k _ { \perp } ) d k _ { \perp } \. \quad \quad ( 7 d )$$

These spectra depend on the streamwise projection k ‖ of the wave-vector k .

- 1D-spectra ‖ E αβ ij ( k ⊥ ), averaged over the ‖ -plane , oriented along U ns . We chose for concreteness the plane ( k , k x z ), such that k ‖ = k x and the spectra depend on k y :

$$\lambda _ { \L } ^ { \L } E _ { i j } ^ { \alpha \beta } ( k _ { y } ) = \int \mathcal { E } _ { i j } ^ { \alpha \beta } \Big { ( } k _ { \L \| }, \sqrt { k _ { z } ^ { 2 } + k _ { y } ^ { 2 } } \, \Big { ) } \frac { d k _ { \L | } \, d k _ { z } } { 4 \pi ^ { 2 } } \. \quad ( 7 e ) \quad \text{In $th$grals}$$

Note that the 2D-, and 1D-energy spectra are defined such that the kinetic energy density per unit mass E j can be found as:

$$\text{be found as} \quad \cdots \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \\ 2 E _ { j } \ & = \ \int \mathcal { E } _ { j j } ( k ) \, \frac { d ^ { 3 } k } { ( 2 \pi ) ^ { 3 } } = \int E _ { j j } ( k _ { \| }, k _ { \perp } ) d k _ { \| } \, d k _ { \perp } \quad \text{$\quad $\quad $\quad $\quad $ \text{This i}$ for th} \\ & = \ \int \tilde { E } _ { j j } ^ { \alpha \beta } ( k, \theta ) \, k \, d k \, d \cos \theta = \int \bullet E _ { j j } ^ { \alpha \beta } ( k ) d k \quad ( 8 ) \quad \text{$\quad $\quad $\quad $\quad $\quad $\quad $ \text{grals} ( \quad \min$} \\ & = \ \int \stackrel { \perp } { \perp } E _ { j j } ^ { \alpha \beta } ( k _ { \| } ) d k _ { \| } = \int \| E _ { j j } ^ { \alpha \beta } ( k _ { y } ) d k _ { y } \. \\ \text{The two common at the time of $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $ \quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad \quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\quad $\ \end{cases}$$

The tensor structure of the energy spectra will be considered only in Sec. II D. In the rest of Sec. II we will restrict ourselves by discussing only scalar versions of the energy spectra, which are the traces of their tensorial counterparts.

## 2. Velocity structure functions

Another presentation of the statistics of turbulence is provided by the second-order velocity structure functions

$$\delta _ { R } u _ { j } ^ { \alpha } \, \equiv \, u _ { j } ^ { \alpha } ( R + r, t ) - u _ { j } ^ { \alpha } ( \, r, t ) \quad ( 9 \mathbf a )$$

The trace S j ( R ) ≡ ∑ α S αα j ( R ) is a measure of the kinetic energy of turbulent (normal or superfluid) velocity fluctuations on scale R . Recently, the streamwise normal velocity across a channel, v x ( y, t ) was measured using thin lines of the triplet-state He 2 molecular tracers created by a femptosecond-laser field ionization of He

$$o _ { R } u _ { j } ^ { - } \ & \equiv \ u _ { j } ^ { - } ( \mathfrak { H } + r, \tau ) - u _ { j } ^ { - } ( \ r, \tau ) \quad \text{(a)} \\ S _ { j } ^ { \alpha \beta } ( \mathfrak { R } ) \ & \equiv \ \left \langle \delta _ { R } u _ { j } ^ { \alpha } \, \delta _ { R } u _ { j } ^ { \beta } \right \rangle \. \quad \text{(9b)}$$

atoms 42 across the channel. This way, the transversal 2 nd -order structure functions 11,12 of the normal-fluid velocity differences S xx n ( R y ) were obtained. Similarly, one can use two or more tracer lines, separated in the streamline direction ̂ x , to measure the longitudinal structure function S xx n ( R x ) and even inclined structure function S xx n ( R ,R x y ).

Using the definition of the structure functions (9a) and the one-dimensional version of the inverse Fourier transform (5d) one gets

$$S _ { j } ^ { x x } ( R _ { y } ) \ = \ \underset { \infty } { 8 } \overset { \infty } { \underset { \infty } { \overset { \| } { \int } } } E _ { j } ^ { x x } ( k _ { y } ) \sin ^ { 2 } \frac { k _ { y } R _ { y } } { 2 } \, d k _ { y } \,, \ \ ( 1 0 \text{a} )$$

$$S _ { j } ^ { x x } ( R _ { x } ) \ = \ \sum _ { \substack { 0 \\ 0 } } ^ { \upsilon } E _ { n } ^ { x x } ( k _ { x } ) \sin ^ { 2 } \frac { k _ { x } R _ { x } } { 2 } \, d k _ { x } \. \ ( 1 0 b )$$

Analyzing the integrals (10) for the scale-invariant spectra E k ( ) ∝ k -m one concludes that they converge in the window of locality

$$1 < m < 3 \. \quad \quad ( 1 1 a )$$

In this window, the leading contribution to the integrals (10) comes from the region kR ∼ 1 and

$$S _ { j } ( R ) \in R ^ { n } \,, \ \ n = m - 1 \. \quad \ ( 1 1 b )$$

This is a well known relationship. For example, n = 2 3 / for the K41 spectrum with m = 5 3 [which satisfy (11a)]. / However for fast decaying spectra with m ≥ 3 the integrals (10a) diverge in the infrared region kR /lessmuch 1 with the main contribution coming from energy containing region k ∼ k 0 , giving

$$S _ { j } ^ { \alpha \alpha } ( R ) \in R ^ { 2 } \,. \quad \quad \ ( 1 1 c )$$

We see that connection between m and n for fastly decaying spectra with m &gt; 3 is lost. To recover it for m &gt; 3 we consider structure functions of the velocity second differences 43,44

$$\Delta _ { R } u _ { j } ^ { \alpha } \equiv u _ { j } ^ { \alpha } ( 2 R + r, t ) - 2 u _ { j } ^ { \alpha } ( R + r, t ) + u _ { j } ^ { \alpha } ( \ r, t ), ( 1 2 \tt a )$$

Now instead of Eqs. (10) we have

$$& \Delta _ { R } u _ { j } ^ { \omega } \equiv u _ { j } ^ { \omega } ( 2 R + r, t ) - 2 u _ { j } ^ { \omega } ( R + r, t ) + u _ { j } ^ { \omega } ( \ r, t ), ( 1 2 a ) \\ & \widetilde { S } _ { j } ^ { \alpha \beta } ( R ) \, \equiv \, \left \langle \Delta _ { R } u _ { j } ^ { \alpha } \Delta _ { R } u _ { j } ^ { \beta } \right \rangle \,. & ( 1 2 b ) \\ \text{ce is} \quad & \text{Now instead of Eqs. (10) we have}$$

$$\widetilde { S } _ { j } ^ { x x } ( R _ { y } ) \, = \, 3 2 \int _ { \substack { 0 \\ 0 \\ \infty } } ^ { \infty } \| E _ { j } ^ { x x } ( k _ { y } ) \sin ^ { 4 } \frac { k _ { y } R _ { y } } { 2 } \, d k _ { y } \,, \, \ ( 1 3 \tt a )$$

$$\widetilde { S } _ { j } ^ { x x } ( R _ { x } ) \, = \, 3 2 \int _ { 0 } ^ { \infty } \overset { \text{v} } { \overset { \text{v} } { \int } } E _ { j } ^ { x x } ( k _ { x } ) \sin ^ { 4 } \frac { k _ { x } R _ { x } } { 2 } \, d k _ { x } \,. \, ( 1 3 b ) \\.. \, \cdots \, \cdots \,. \, \cdots \,. \,.$$

Now these integrals converge in the extended window of locality

$$1 < m < 5 \. \quad \quad \ ( 1 4 \text{a} )$$

TABLE I: Parameters of simulations by columns: (# 1) Run #, (# 2); Temperature (K); (# 3 and # 4) Ratios of the superfluid and normal-fluid densities 29 , ρ /ρ s n and viscosities ν /ν s n ;(# 5) The mutual friction parameter α ; (# 6 and #7) The numerical values of the kinematic viscosity of the normal-fluid and superfluid components ˜ ν n and ˜ ; (# 8 and #9) The numerical values ν s of mutual friction frequency Ω and counterflow velocity V ; (# 10 and #11) the root mean square (rms) normal-fluid and superfluid turbulent velocity fluctuations v n T and v s T ; (# 12 and # 13) the normal-fluid and superfluid Reynolds numbers; (# 14) k × , Eq. (16d).

For details see Sect.II A. In all simulations: the number of collocation points along each axis is N = 256; the computational box size is L = 2 π , the range of forced wavenumbers k ˜ ϕ = [0 5 . , 1 5]. .

| 1   | 2    | 3       | 4       | 5    | 6       | 7       | 8   | 9   | 10            | 11            | 12   | 13   | 14   |
|-----|------|---------|---------|------|---------|---------|-----|-----|---------------|---------------|------|------|------|
| Run | T ,  | ρ s ρ n | ν s ν n | α    | ˜ ν n , | ˜ ν s , | Ω,  | V   | u n T =       | u s T =       | Re n | Re s | k ×  |
| #   | K    |         |         |      | × 10 3  | × 10 3  |     |     | √ 〈 u n 2 〉 | √ 〈 u s 2 〉 |      |      |      |
| 1   |      |         |         |      |         |         | 1   | 15  | 4.3           | 4.5           | 997  | 2595 | 0.34 |
| 2   | 1.65 | 4.2     | 0.49    | 0.11 | 3.0     | 1.38    | 20  | 15  | 4.2           | 4.2           | 1216 | 2667 | 6.88 |
| 3   |      |         |         |      |         |         | 1   | 0   | 3.6           | 3.6           | 1257 | 2953 | ∞    |
| 4   |      |         |         |      |         |         | 20  | 0   | 3.7           | 3.7           | 1316 | 2872 | ∞    |
| 5   |      |         |         |      |         |         | 1   | 15  | 3.4           | 3.5           | 908  | 994  | 0.18 |
| 6   | 1.85 | 1.75    | 1.07    | 0.18 | 3.0     | 1.85    | 20  | 15  | 3.5           | 3.5           | 1051 | 1056 | 3.57 |
| 7   |      |         |         |      |         |         | 1   | 0   | 3.3           | 3.3           | 1154 | 1239 | ∞    |
| 8   |      |         |         |      |         |         | 20  | 0   | 3.6           | 3.5           | 1179 | 1181 | ∞    |
| 9   |      |         |         |      |         |         | 1   | 15  | 4.3           | 4.2           | 1064 | 582  | 0.15 |
| 10  | 2.00 | 0.83    | 1.72    | 0.28 | 3.0     | 5.0     | 20  | 15  | 3.5           | 3.5           | 1153 | 664  | 1.5  |
| 11  |      |         |         |      |         |         | 1   | 0   | 3.3           | 3.3           | 1225 | 689  | ∞    |
| 12  |      |         |         |      |         |         | 20  | 0   | 3.6           | 3.5           | 1177 | 676  | ∞    |

In this window, the leading contribution to the integrals (13) again comes from the region kR ∼ 1 and, similarly to Eq.(11b), for the scale-invariant spectrum E k ( ) ∝ k -m we have

## C. Physical origin of the strong anisotropy of counterflow turbulence

$$\widetilde { S } _ { j } ^ { \alpha \alpha } ( R ) \in R ^ { n } \,, \quad n = m - 1 \,. \quad \ ( 1 4 b ) \quad \text{ anis} \\ \times 5 \, & \, t h e \, \text{into rata} \, ( 1 3 ) \, \text{divero} \, \text{in the infrared region}$$

For m&gt; 5 the integrals (13) diverge in the infrared region and

$$\widetilde { S } _ { j } ^ { \alpha \alpha } ( R ) \in R ^ { 4 } \,. \quad \quad ( 1 4 c ) \quad \begin{matrix} & \text{ial} \\ & \\ & \\ \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot & \cdot \\ \end{matrix}$$

It is worth noting that the relations (11b) and (14b) are valid in the limit of infinite inertial interval. For a finite inertial interval which is typical for the experimental conditions, the structure functions have a complicated functional dependence, mixing the inertial and viscous behavior and the original scaling of the energy spectra is reproduced over very short intervals of scales 13 . Nevertheless, when experimental conditions do not allow to measure the energy spectra directly, the structure functions remain the preferred tool to access the statistics of the velocity fluctuations. In a turbulent counterflow, where the energy spectra are not scale-invariant, the quantitative analysis of the structure functions may not be meaningful. Nevertheless a qualitative difference between structure functions measured along different directions may confirm the presence of the spectral anisotropy.

In a recent Letter 19 it was shown that the energy spectra in counterflow turbulence are expected to be strongly anisotropic. To keep the present paper self-contained, we repeat here some of that discussion and add further clarifications to the analysis.

We start with a balance equation 18,19 for the 2D energy spectrum ˜ ( E j k, θ ) in the counterflow turbulence with axial symmetry:

$$\underline { \tilde { E } _ { j } ( \iota } }$$

$$b ) \, \text{are} \quad \frac { \partial \tilde { E } _ { j } ( k, \theta, t ) } { \partial t } + \text{div} _ { k } [ \varepsilon _ { j } ( k ) ] = - \mathcal { D } _ { j } ^ { \text{mf} } ( k, \theta ) - \mathcal { D } _ { j } ^ { \nu } ( k, \theta ), ( 1 5 a ) \\ \text{finite} \quad \left ( \mathcal { D } _ { j } ^ { \text{mf} } ( k, \theta ) \, & = \, \Omega _ { j } [ \tilde { E } _ { j } ( k, \theta ) - \tilde { E } _ { n s } ( k, \theta ) ] \,, \quad \left ( 1 5 b \right ) \\ \text{func} \quad \left ( \mathcal { D } _ { j } ^ { \nu } ( k, \theta ) \, & = \, 2 \, \nu _ { j } k ^ { 2 } \tilde { E } _ { j } ( k, \theta ) \,. \quad \left ( 1 5 c \right )$$

Here div k [ ε j ( k )] is the transfer term due to inertial nonlinear effects. The terms on the right hand side describe the energy dissipation rate due to the mutual friction D mf j ( k, θ ) and due to the viscous effects D ν j ( k, θ ). To keep the presentation concise, we introduced a notation Ω n = Ω s ρ /ρ s n .

The origin of the energy spectra anisotropy in counterflow turbulence can be deduced from the form of the dissipation rate D mf j ( k, θ ) (15b). In this term, the crosscorrelation function ˜ E ns ( k, θ ) has the following form [cf.

FIG. 1: The spherically-averaged energy spectra and the cross-correlation functions for the counterflow and coflow at different temperatures. Panels (a),(b) and (c) - the K41-compensated energy spectra · E j ( k ) for the normal-fluid (solid lines) and superfluid components (dashed lines). Panels (d),(e) and (f) - the normalized cross-correlation functions · R k ( ). The four sets of lines in each panel correspond to V 0 Ω - green lines, 1 V 0 Ω 20 - brown lines, V 15 Ω - red lines, 1 V 15 Ω 20 - blue lines. The black lines, labeled 'Cl' in the panels (a),(b) and (c) correspond to the spectra of classical turbulence.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-65836c7f-a46d-47a9-a999-5174261ace39.jpg)

Eq.(13) in Ref. 17 ]:

$$\tilde { E } _ { n s } ( k, \theta ) = \frac { A B } { B ^ { 2 } + ( k \cdot U _ { n s } ) ^ { 2 } } \. \quad \ ( 1 6 \text{ \ \ } \text{(1.6)} \quad \text{ity} \quad \text{app}$$

Here A = Ω s ˜ E n ( k, θ )+Ω n ˜ ( E s k, θ ) and B can be written 18 as B = Ω ns = Ω n + Ω . s We further note 18,19 that when two components are highly correlated, the crosscorrelation may be accurately represented by the corresponding energy spectra. For wave numbers where the components are not correlated, ˜ E ns ( k, θ ) is small and the accuracy of its representation is less important. This allows us to decouple ˜ E ns ( k, θ ) in Eqs. (15) for each component as follows:

$$\tilde { E } _ { n s } ( k, \theta ) \ = \ \tilde { E } _ { j } ( k, \theta ) D ( k, \theta ) \,, \quad \quad ( 1 6 b ) \quad \text{case} _ { \text{the} \ \} }$$

$$\mathcal { L } _ { n s } ( \varkappa, \nu _ { J } ) \, = \, \mathcal { L } _ { j } ( \varkappa, \nu _ { J } ) \mathcal { L } ( \varkappa, \nu _ { J } ) \,, \quad \, \{ \i \i \i \rho \, \quad \text{the} \\ D ( k, \theta ) \, = \, 1 \Big / \Big [ 1 + \Big ( \frac { k U _ { n s } \cos \theta } { \Omega _ { n s } } \Big ) ^ { 2 } \Big ] \,. \quad ( 1 6 c ) \quad \text{cal} \, \quad \text{with} \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdots \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \end{cases}$$

Note that averaging Eq. (16c) over θ results in the equation for D k ( ), used in the theory of isotropic counterflow turbulence 18 :

$$\i u r o u e n c e \, \begin{matrix} & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & } & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ( & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & + & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & _ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & = & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ^ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & | & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & b & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

The function D k,θ ( ) in Eqs. (16) describes the level of decorrelation of the normal-fluid and superfluid velocity components by the counterflow velocity. Within the approximation (16b), D k,θ ( ) defines the rate of energy dissipation caused by mutual friction:

$$\mathcal { D } _ { j } ^ { \text{mf} } ( k, \theta ) = \Omega _ { j } \tilde { E } _ { j } \left [ 1 - D ( k, \theta ) \right ] \,. \quad \ \ ( 1 7 )$$

For small k or even for large k with k almost perpendicular to U ns (i.e cos θ /lessmuch 1), D k,θ ( ) /similarequal 1 , the normal and superfluid velocities are almost fully coupled and the rate of the energy dissipation D mf j ( k, θ ) /lessmuch Ω is small. In this j case, the role of mutual friction is minor and we expect the energy spectrum ˜ ( E j k, θ ) to be close to the classical prediction E K41 ( k ) ∝ k -5 / 3 . For large k and for k with cos θ ∼ 1, D k,θ ( ) /lessmuch 1, the velocity components are almost decoupled, D j ( k, θ ) /lessmuch 1 and the mutual-friction energy dissipation is maximal: D mf j ( k, θ ) ≈ Ω . j This situation is similar to that in 3 He with the normal-fluid component at rest, for which D mf s ( k, θ ) = Ω . s In this case, we can expect that the energy dissipation by mutual friction strongly suppresses the energy spectra, much below the K41 expectation E K41 ( k ) up to the level typical for the 3 He turbulence 21,22,36,45 . Next, we note that Ω ns = Ω +Ω s n = ρ Ω s /ρ n . Then Eq.(16c) may be rewrit-

FIG. 2: The angular dependencies of the energy spectra ˜ ( E j k, θ ) / · E j ( k ) and the cross-correlations R k, θ ( ), averaged in 3 wavenumber bands, for various flow condition. Panels (a),(b) and (c)- the 2D energy spectra ˜ ( E j k, θ ). The spectra of the normal-fluid are shown by solid lines and of the superfluid - by dashed lines. Panels (d)-(f)- the cross-correlations R k, θ ( ). The spectra and the cross-correlation for the coflow are shown by thin horizontal lines and marked V 0 Ω 20 , for the counterflow - by thick lines and marked V 15 Ω 20 . In all panels, red lines correspond to the averaging over wavenumber range 10 ≤ k &lt; 20 (labeled as k 10 ), green lines - to averaging over 20 ≤ k &lt; 60 (labeled as k 20 ) and blue lines - to the averaging over 60 ≤ k ≤ 80 (labeled as k 60 ). Note the log-linear scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-ce5460b4-a834-47cf-a063-ab95e916bb0a.jpg)

ten as

## II. STRONG ANISOTROPY OF ENERGY SPECTRA

## A. Simulation parameters and numerical procedure

The direct numerical simulation of the coupled Eqs. (1b) were carried out using a fully de-aliased pseudospectral code with a resolution of 256 3 collocation points in a triply periodic domain of size L = 2 π . The parameters of the simulations are summarized in Table I. To obtain the steady-state evolution, velocity fields of the normal and superfluid components are stirred by two independent isotropic random Gaussian forcings:

$$\langle \tilde { \varphi } _ { u } ( k, t ) \cdot \tilde { \varphi } _ { u } ^ { * } ( k ^ { \prime }, t ^ { \prime } ) \rangle = \Phi ( k ) \delta ( k - k ^ { \prime } ) \delta ( t - t ^ { \prime } ) \widehat { P } ( k ) \,, \, ( 1 9 ) \\ \mathrm e \, \text{en-\quad where } \widehat { P } ( k ) \, \text{is a projector assuring incompressibility and}$$

where ̂ P ( k ) is a projector assuring incompressibility and Φ( ) = Φ k 0 k -3 ; the forcing amplitude Φ 0 is nonzero only in a given band of Fourier modes: k ∈ [0 5 . , 1 5]. . Both components are forced with the same amplitude to allow direct comparison with simulations of the uncoupled equations. The time integration is performed using

$$D ( k, \theta ) = 1 \Big / \left [ 1 + \left ( \frac { \rho _ { n } k U _ { n s } \cos \theta } { \rho \Omega _ { s } } \right ) ^ { 2 } \right ]. \quad ( 1 8 ) \quad T |$$

At low temperatures ρ /ρ n is small and the velocity decorrelation is considerable only at large k , while at higher temperature ρ /ρ n ≈ 1 and the energy dissipation by mutual friction is effective at all scales.

Combining all these considerations, we expect the energy spectra ˜ ( E j k, cos θ ) to become more anisotropic with increasing k , with most of energy concentrated in the range of small cos θ , i.e. in the wavevector plane orthogonal to the counterflow velocity U ns . This effect is milder at low T and stronger at higher temperatures.

a 2-nd order Adams-Bashforth scheme with the viscous term exactly integrated. Simulations for all temperatures were carried out with the normal-fluid viscosity fixed at ˜ ν n = 0 003 and the value of ˜ . ν s is found using the known value of ratio ν /ν s n at each temperature.

To properly expose various aspects of the counterflow turbulence statistics we chose several sets of governing parameters for the simulations. Since the material parameters of 4 He are strongly temperature dependent 29 (see Table I, columns #3-5), we consider three temperatures, corresponding to an experimentally accessible range T = 1 65 1 85 and 2 0K. At low temperatures the . , . . superfluid component is dominant and has a lower viscosity, while at high T the density of the normal-fluid component is larger, while its kinematic viscosity is lower. At T = 1 85K, the densities and the viscosities of two . components are closely matched.

As was shown in previous studies 17,18,33 , the major role in the statistics of the counterflow in superfluid 4 He is played by the ratio of the mutual friction frequency and the counterflow Doppler frequency Ω ns / kU ( ns ) [cf. Eq.(16d)]. To explore various scenarios, we use one counterflow velocity U ns and two very different values of Ω . s To emphasis the importance of the flow anisotropy, we compare the results for the counterflow with simulations of coflow U ns = 0, keeping the rest of parameters unchanged. The detailed study of the statistics of the coflow was reported in Ref. 20 . Other parameters of the simulations were chosen based on dimensionless numbers: i)

the Reynolds numbers Re j = ∆ u ν k j 0 , ii) the turbulent in- tensity U ns ∆ u , and iii) the dimensionless cross-over scale q × = Ω ns k U 0 ns = k × k 0 . Here ∆ u is the root-mean-square

(rms) of the normal-fluid turbulent velocity fluctuations, k 0 = 1 is the outer scale of turbulence. The numerical values of the dimensionless counterflow velocity V = 15 and the mutual friction frequency Ω = 1 and 20 are listed in Table I, columns #8-9. In this way, for each temperature we have 4 runs, labeled below as V 0 Ω , 1 V 0 Ω 20 , V 15 Ω 1 and V 15 Ω 20 . Reference simulations of the uncoupled Eq. (1a) with Ω s = 0 , U ns = 0, representing classical hydrodynamic turbulence for the same parameters of the flow, are labeled as 'Cl'.

The correlation between components in the counterflow become gradually weaker with increasing cos θ for T = 1 65K, while for . T = 2 0K the normal fluid and . superfluid are essentially uncorrelated for cos θ /greaterorsimilar 0 1. .

## B. Spherical energy spectra · E jj ( k ) and cross-correlations · E ns ( k )

The energy spectra are influenced by a few competing factors: the viscous dissipation, the dissipation by mutual friction and the decoupling due to counterflow velocity. To find their relative importance, we first ig- nore the present angular dependence and consider the spherically averaged spectrum · E j ( k ) and the normalized cross-correlation function:

$$\bullet R ( k ) = \frac { 2 \bullet E _ { n s } ( k ) } { \bullet E _ { n n } ( k ) + \bullet E _ { s s } ( k ) } \,. \quad \quad ( 2 0 )$$

In the upper row of Fig.1 we plot the spectrum · E j ( k ), compensated by the classical scaling k 5 / 3 for different flow conditions. In the lower row we show the crosscorrelations Eq. (20). In this figure and in Figs. 2-7 the results for T = 1 65K are shown in the left column [panels . (a) and (d)], for T = 1 85K- in the middle column [pan-. els (b) and (e)] and for T = 2 0 K - in the right column . [panels (c) and (f)]. The effect of viscous dissipation is clearly seen in the spectra of the uncoupled components, corresponding to classical turbulence and marked 'Cl', black lines. The spectra almost coincide for T = 1 85K, . for which the viscosities are almost equal. The viscosity of the normal-fluid component (solid lines) is larger than for the superfluid (dashed lines) for T = 1 65K and . smaller for T = 2 0 K. .

Next, we add the coupling by the mutual friction force, creating a coflow (green and brown lines). The strongly coupled components ( V 0 Ω 20 , brown lines) are well correlated at all scales and move almost as one fluid. The corresponding spectra slightly differ only at the viscous scales. Note the additional dissipation due to mutual friction, leading to further suppression of the spectra compared to the uncoupled case, for T = 1 85 and 2 0K. . . At the lower temperature T = 1 65K the . energy exchange between components leads to stronger dissipation in the superfluid component and weaker dissipation in the normal-fluid component. For weaker coupling ( V 0 Ω , 1 green lines), the situation is completely different. The components are almost uncorrelated, especially at large k . The coupling between them is translated into very efficient dissipation by mutual friction, leading to spectra that are suppressed almost at all scales, especially at high temperature (see Ref. 20 for details).

In the presence of the counterflow velocity 17 , the two components are swept in opposite directions by the corresponding mean velocities. This leads to further decorrelation of the component's turbulent velocities, especially at small scales, for which the overlapping time is very short [cf. lines for V 15 Ω and 1 V 15 Ω 20 in Fig. 1(d-f)]. Even for the strong coupling, V 15 Ω 20 , blue lines, the velocities become progressively less correlated for all temperatures. The dissipation by mutual friction is very strong in this case, with both Ω and the velocity difference being large, leading to very strongly suppressed spectra, with · E nn ( k ) ≈ · E ss ( k ). At T = 1 65K there is still . some interval of scales with k /greaterorsimilar k 0 , for which the spectra are close to K41 scaling. The crossover scale agrees well with k × ≈ 7 for this case (see Table I, column #14). For higher temperatures this crossover scale become smaller and the classical-like behavior is not resolved. At weak coupling ( V 15 Ω , red lines), the velocities are essentially 1 uncorrelated and the spectra of the two components dif-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-ac1d4513-9f88-4f3e-b445-f54d0dc92465.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-bfcd4577-430a-455e-a7ab-da08e4e56de7.jpg)

FIG. 3: The superfluid component energy spectrum E s ( k || , k ⊥ ) in the counterflow. T = 2 0K. Note the difference in the magnitudes, shown by the color-bar range. .

(a) T = 1 65K, (b) . T = 1 85 K, (c) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-8ee07789-4e50-4af3-9f67-e073ed5ca148.jpg)

FIG. 4: Tensor decomposition of the 1D energy spectra K α j for the normal fluid (circles) and for the superfluid components (squares). The spectra for the coflow are labeled V 0 Ω 20 , for the counterflow V 15 Ω 20 . The K x j , K y j and K z j tensor components are annotated and denoted by red, blue and green lines, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-532b57a1-fd81-4e57-af53-5c9a91891752.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-0846f9d1-c15f-4513-a26a-240f269a3f92.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-ea1563da-5ae6-45f2-ac49-041e50509214.jpg)

fer and are very strongly suppressed, especially at high T .

## C. Angular dependence of 2D-energy spectra

The behavior of the spherically averaged energy spectra agrees well with the predictions of the theory 18 , based on the assumption of spectral isotropy. To explore the angular dependence of the energy spectra and the correlations ˜ E ij ( k, θ ) we plot in Fig. 2(a-c) the spectra ˜ ( E j k, θ ), normalized by the corresponding · E j ( k ) and in Fig. 2(d-f) the corresponding normalized cross-correlations ˜ ( R k, θ ). Given the discrete nature of the k -space in DNS, we further average them over 3 bands of wavenumbers. We do not account for the largest scales k ≈ k 0 which are influenced by the forcing and average the spectra and the cross-correlations over the k -ranges 10 ≤ k &lt; 20, 20 ≤ k &lt; 60 and 60 ≤ k ≤ 80. The corresponding lines are labeled as k 10 , k 20 and k 60 , respectively. Here we consider only strong coupling regime and plot the spectra and the cross-correlations for the coflow ( V 0 Ω 20 ) and for the counterflow ( V 15 Ω 20 ) .

for the counterflow, on the other hand, have a complicated form. Both the spectra and cross-correlation are largest for cos θ ≈ 0 and fall off very quickly with decreasing angle. The spectra decrease exponentially with cos θ , slower for small k (red lines for k 10 ) and faster for larger k (green and blue lines for k 20 and k 60 , respectively). This effect is stronger for the normal-fluid (super fluid) component at low temperatures (high temperature). Most of the energy is contained in the narrow range cos θ &lt; 0 1, . near the ⊥ -plane in the k -space, orthogonal to U ns .

The first observation is that the spectra and the crosscorrelation for the coflow are isotropic for all the conditions. The angular dependencies of ˜ ( E j k, θ ) and ˜ ( R k, θ )

To better quantify the angular energy distribution, we use the fact that the spectra ˜ ( E j k, θ ) have piecewise exponential dependence of cos θ , as is evident from Fig. 2(a)(c). We then estimate the cos θ -range, in which half of the total energy is contained, for different wavenumber bands. At T = 1 65K, for the small wavenumbers . k 10 band, this range is indeed cos( θ ) &lt; 0 1 for both the nor-. mal an superfluid components. With increasing temperature, this range decreases to cos θ &lt; 0 05 for . T = 1 85 K . and to cos θ &lt; 0 03 for the normal-fluid and cos . θ &lt; 0 025 . for the superfluid at T = 2 0 K. For the larger . k 20 band these values are 0 04 and 0 06 for the normal fluid and . . the superfluid, respectively, at T = 1 65K. For higher . temperatures, as well as for the high wavenumber k 60 band, about a half of the total energy is contained in a narrow cos θ &lt; 0 02 . -0 025 range for both components. .

௬

ෝ

ݕ

ෝ

ݖ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-3eb2ae86-cc98-4a7b-bb3f-154908527daf.jpg)

ୱ

ݑ

ୱ

௫

FIG. 5: Superfluid velocity u x s and u y s components 50 . T = 1 85 K, Ω = 20, . V = 15. The u z s component (not shown) is similar to the y component. The velocity magnitude is color coded, with red denoting positive and blue denoting negative values.

Indeed, the superfluid energy spectrum E ss ( k || , k ⊥ ), shown in Fig. 3, is strongly suppressed in k || direction, while decreasing slowly in the orthogonal direction, especially for T = 2 0 K. .

## D. Tensor structure of 1D energy spectra

Given such a strong anisotropy of the spectra in the counterflow, it is natural to expect that different components of the turbulent velocity fluctuations are excited to a different extent. In this section we consider the tensor structure of 1D-energy spectra · E αα jj ( k ) for α = x, y, z and clarify which components ( v x j along U ns or v y j , v z j , both orthogonal to U ns ) are most excited.

In Figs. 4 we plot the components of the spherical spectra for three temperatures as the ratios

$$\bullet K _ { j } ^ { \alpha } ( k ) \equiv 3 \bullet E _ { j j } ^ { \alpha \alpha } ( k ) / \bullet E _ { j j } ( k ) \. \quad \ \ ( 2 1 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The factor 3 was introduced to ensure that for the isotropic turbulence ∑ α = x,y,z · K α j ( k ) = 1.

Indeed, for the coflow (the almost horizontal lines, labeled V 0 Ω 20 ) all the velocity components are excited equally, except for the smallest wavenumbers. On the other hand, for the counterflow turbulence (lines labeled V 15 Ω 20 ) the contribution of the · K x j ( k ) component (shown by red lines) is dominant and monotonically increasing with k from the isotropic level · K x j ( k 0 ) ≈ 1 to the maximal possible level · K x j ( k ) ≈ 3. This means that the small-scale counterflow turbulence mainly consists of v x j ( k ) velocity fluctuations. The contribution of v y j and v z j fluctuations for k /greaterorsimilar 10 is negligible, especially at T = 2 0 K. .

Therefore, the counterflow turbulence represent a special kind of a quasi-2D turbulence, consisting mostly of the turbulent velocity fluctuations with only one streamwise projection u ‖ , which depends on the cross-stream coordinate r ⊥ : u ‖ ( r ⊥ , t ). This behavior is essentially different from other known types of quasi-2D turbulence, such as stably-stratified flow in the atmosphere 23-25 or rotational turbulence 26-28 , in which the leading contribution to the turbulent velocity field comes from the 2D velocity field u ⊥ that depends on r ⊥ : u ⊥ ( r ⊥ , t ). Such a u ‖ ( r ⊥ , t )-turbulence can be visually presented as narrow jets or thin sheets as illustrated in Fig. 5.

Note the difference with the strong acoustic turbulence. There the velocity field has tangential velocity breaks at the jets boundaries and the 1D energy spectrum E k ( ) ∝ k -2 . The energy spectra in the counterflow decay much faster. It means that the velocity fields at the jets boundaries are continuous together with some finite number of their derivatives. This is a consequence of the mutual friction that tends to smooth the velocity field.

## E. Comparison of 1D energy spectra and reconstruction of 3D spectra

The best way to study the anisotropy of hydrodynamic turbulence is to expand the statistical objects in the irreducible representations of the SO(3) symmetry group, see, e.g. Refs. 25,46-49 . In counterflow turbulence, an attempt to expand ˜ ( E j k, θ ) into a series with respect to Legendre polynomials,

$$\tilde { \ E } _ { j } ( k, \theta ) = \sum _ { \ell } \, \tilde { E } _ { j } ( k, \ell ) \, P _ { \ell } ( \cos \theta ) \,, \quad \ \ ( 2 2 )$$

and to study k -behavior of ˜ ( E j k, θ ), turned out to be ineffective. The very strong anisotropy of ˜ ( E j k, θ ) spectra required too many terms in the expansion(22) for an adequate reproduction of its angular dependence. Therefore, we choose another way to characterize the spectral anisotropy, which is more suitable in our case. We compare the normal-fluid and superfluid spherical · E j ( k ), cylinder ◦ E j ( k ⊥ ), and ‖ -, ⊥ -plane-averaged energy spectra ‖ E j ( k ⊥ ), ⊥ E j ( k ‖ ).

In the case of isotropy, all the four 1D energy spectra are proportional to each other

$$\L e ^ { \lambda } _ { i j } ( p ) \, \in \, ^ { \circ } E ^ { \alpha \beta } _ { i j } ( p ) \, \in \, ^ { \perp } E ^ { \alpha \beta } _ { i j } ( p ) \, \in \, ^ { \| } E ^ { \alpha \beta } _ { i j } ( p ) \,, \quad ( 2 3 )$$

differing only in numerical prefactors. Here p is the corresponding (dimensional, [ p ] =1/cm) wavenumber: p = k, k ⊥ , k ‖ or k y . By estimating contributions to the integrals in Eqs. (7a), (7c) and (7e) in the case of strong anisotropy (i.e coming from a narrow range with k ‖ /lessmuch k ⊥ ), one may show that the spectra are related as · E j ( p ) ≈ ◦ E j ( p ) ≈ C E ‖ j ( p ), where C is a numerical prefactor. This fact may explain the good agreement

FIG. 6: Various 1D energy spectra of the normal-fluid (solid lines) and the superfluid (dashed) components for the coflow (lines labeled V 0 Ω 20 ) and for the counterflow (lines labeled V 15 Ω 20 ). Different kinds of spectra are annotated in the figure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-0383b31e-de7c-4630-8770-ed4a54fefc7f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-72cccb1c-c8a9-4c5a-92d1-bc6844d05ee8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-b07415e8-09cd-4b6b-8ad3-0e65a0aa443f.jpg)

FIG. 7: The ⊥ -plane energy spectra ⊥ E j ( k ‖ ) of the normal-fluid (solid lines) and the superfluid (dashed) components for counterflow V 15 Ω 20 . The spectra are shown in a Log-linear scale in panel (a) and in Log-log scale in panel (b). Different temperatures are annotated in the figure. Note linear behavior at small k ‖ in panel (a) and at large k ‖ in panel (b) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-6ee71cab-40ef-4de9-807f-9021ac44c0d4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-60683df7-8400-402e-b43b-1b9f28a0c16a.jpg)

between the experimental spectra ‖ E n ( k ⊥ ), obtained in Ref. 13 and the prediction of the theory 18 for · E j ( k ). The integral in Eq. (7d) is different and the spectrum ⊥ E j ( k ‖ ) is expected to be confined to small k ‖ range.

These spectra, normalized by the energy density E n and compensated by the K41 factor p 5 / 3 , are shown in Fig. 6. The coflow spectra, appearing as almost horizontal lines, labeled V 0 Ω 20 , indeed differ by less than an order of magnitude for all T . The relation between various spectra for the counterlow is consistent with the above estimate, further confirming the strong spectral anisotropy. The degree to which the ⊥ E j ( p ) spectra, shown by green lines, are suppressed at different temperatures, agrees with the angular dependence of ˜ ( E j k, θ ), Fig. 2(a)-(c). While at T = 1 65 K the spectrum for the . k 10 -range at cos θ ≈ 1 is smaller by three orders of magnitude in the direction of the counterflow compared to the orthogonal plane, at T = 2 0 K this difference is almost ten orders of . magnitude. Accordingly, the ⊥ E j ( p ) spectrum at T = 2 0 . K is confined to less than a decade in p . To better quantify the steepness of the spectra we list in Table II the values of the ratios E j (10) /E j (1) for ⊥ E j ( p ) and || E j ( p ).

The analysis of the ( k, θ )-dependence of the ˜ ( E j k, θ ) energy spectra in Sec. II C showed that the overwhelming

TABLE II: The steepness of the energy spectra ⊥ E j ( k ) and ‖ E j ( k ) in the counterflow, characterized by the ratios ⊥ E j (10) / ⊥ E j (1) and ‖ E j (10) / ‖ E j (1).

|                            | T = 1 . 65K    | T = 1 . 85K    | T = 2 . 0K     |
|----------------------------|----------------|----------------|----------------|
| ‖ E n,s (10) / ‖ E n,s (1) | 10 - 3         | 2 . 0 × 10 - 4 | 2 . 5 × 10 - 4 |
| ⊥ E n (10) / ⊥ E n (1)     | 6 . 3 × 10 - 6 | 2 . 0 × 10 - 8 | 1 . 6 × 10 - 8 |
| ⊥ E s (10) / ⊥ E s (1)     | 2 . 5 × 10 - 5 | 3 . 2 × 10 - 8 | 1 . 6 × 10 - 9 |

part of the total turbulent energy is concentrated in the range of small cos θ , and consequently small k ‖ , say for k ‖ /lessorsimilar 10. For a semi-qualitative analysis of the 2D-spectra E j ( k , k ‖ ⊥ ) and ˜ ( E j k, θ ) in this range of k || , we assume a factorization

$$E ( k _ { \| }, k _ { \perp } ) \simeq f _ { 1 } ( k _ { \| } ) \, f _ { 2 } ( k _ { \perp } ) \. \quad \ \ ( 2 4 \tt a )$$

If so, using Eqs. (7c) and (7d), we can reconstruct the 2D energy spectra as follows

$$E _ { j } ( k _ { \| }, k _ { \perp } ) \simeq \frac { \overset { \circ } { E } _ { j } ( k _ { \perp } ) \, ^ { \perp } E _ { j } ( k _ { \| } ) } { E _ { j } } \,, \quad ( 2 4 b )$$

where E j is the energy density in the system given by Eq.(8).

FIG. 8: The velocity structure functions of the normal-fluid component S || n ( R ) (9b). The lines for the coflow are marked V 0 Ω 20 , for the counterflow V 15 Ω 20 . Various structure functions are marked in the figures. The dashed lines, labeled R 2 / 3 and the dot-dashed lines, marked R 2 , serve to guide the eye only.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-2b733a01-71b1-4403-818c-7827ab122230.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-3e6eb06a-e12c-4527-bcc5-552347201e96.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-d13ff9da-d3cd-4811-bc40-5b7fbceb5438.jpg)

FIG. 9: The structure functions S || j ( R || ) (a) and S || j ( R ⊥ ) (b) at T = 1 65 K. The lines for the counterflow with . V = 15 and different values of Ω are color-coded: Ω = 1 (red lines), Ω = 3 ( green lines), Ω = 20 ( blue lines). The structure functions for the normal-fluid component are shown by solid lines and for the superfluid component-by dashed lines. Structure functions for the classical turbulence ( V = 0 Ω = 0) are shown by black lines and labeled 'Cl'. ,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-a41c209b-250a-4f25-9f00-0efbe0e521d8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-bc16714f-a10d-4574-9dce-dda43e121674.jpg)

Furthermore, using Eqs. (6), we can also reconstruct the 2D spectra ˜ ( E k, θ ) from E k , k ( ‖ ⊥ ):

$$\tilde { E } _ { j } ( k, \theta ) \simeq \frac { \mathring { \mathrm E } _ { j } ( k \sin \theta ) \, ^ { \perp } E _ { j } ( k \cos \theta ) } { E _ { j } \sin \theta } \. \quad ( 2 4 c ) \quad \text{prop}$$

scale. At all temperatures the smallk ‖ behavior is exponential, while at larger k ‖ the spectra are consistent with the power-law behavior. Using this information, we propose the following form for the smallk ‖ -spectra:

In the range of small cos θ , where most of the turbulent energy is concentrated, Eq.(24c) can be simplified as follows: ˜ ( E j k, θ ) /similarequal ◦ E j ( k ) ⊥ E j ( k cos θ /E ) j . The θ -dependence of ˜ ( E j k, θ ) is therefore determined by ⊥ E j ( k cos θ ), i.e. cos θ appears in the combination k cos θ = k ‖ . This observation fully agrees with our theoretical prediction, that cos θ appears in the theory only via Eq. (16c) in the dimensionless factor kU ns cos θ/ Ω ns . We consider this agreement as an argument in a favor of the factorization assumption (24a) for small cos θ .

To take a closer look at ⊥ E j ( k ‖ ), we plot in Fig. 7 these spectra for different temperatures. To expose the functional dependence of the spectra, we use different scales in two panels: in the panel (a) the scale is Log-Linear, while in panel (b) the spectra are plotted in the Log-Log

$$\stackrel { \perp } { ^ { \perp } } E _ { j } ( k _ { \| } ) \simeq \frac { k _ { \| } \, \mathcal { E } _ { j } } { k _ { * } } \, \exp \left [ \, - \frac { k _ { \| } } { k _ { * } } \right ] \,. \quad \quad ( 2 5 )$$

It is tempting to relate the characteristic k ∗ to the crossover scale k × : k ∗ ∝ Ω ns /U ns = k × . Indeed, k ∗ estimated from Fig. 7(a) and k × (Table I, column # 14) have similar temperature trends. This gives additional support for factorization (24a) and for qualitative theoretical discussion of the problem in Sec. I C.

The observed steep power-law behavior of ⊥ E j ( k ‖ ) for larger k ‖ /greaterorsimilar 10, Fig. 7(b) with an apparent exponents m ≈ -7 may indicate a nonlocal energy transfer between largest and smallest scales, similar to the super-critical spectra in the superfluid 3 He.

FIG. 10: The flatness F || j ( R || ) (red lines) and F || j ( R ⊥ ) (blue lines) T = 1 85 K in the counterflow Ω . s = 20 , V = 15. The flatness for the normal-fluid component are shown by solid lines and for the superfluid component-by dashed lines. The corresponding flatnesses for the classical turbulence ( V = 0 Ω = 0) are shown by black and gray lines and labeled 'Cl'. ,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-c9414b67-93ad-4ec3-8c49-264f9e5d0b2d.jpg)

## F. The structure functions

The energy spectra ‖ E j ( k ⊥ ), ⊥ E j ( k ‖ ) may be translated into the corresponding structure functions, according to Eqs. (10) and (13). In Fig. 8 we show the structure functions (9b) S || n ( R ) ≡ S xx n ( R ) with the velocity differences taken in the direction of the counterflow R || ≡ R x and in the plane R ⊥ orthogonal to it. The structure functions (9b) for the coflow, shown in Fig. 8(a-c) as red and orange lines are similar to classical turbulence; at large scales they follow approximately R 2 / 3 scaling, gradually crossing over towards viscous R 2 behavior. The transition is very broad here, but the two apparent scaling ranges are evident. The cross-over scale increases with temperature. The structure functions, calculated along and across the counterflow direction, are similar, slightly differing mostly in the magnitude at all scales. The main difference from the uncoupled case (not shown) is the lower magnitude at all scales, reflecting the presence of addition energy dissipation by mutual friction. In the counterflow, the situation is different. Over most of the available range of scales, the structure functions, shown as dark and light blue lines in Fig. 8, appear to have an apparent scaling behavior close to R 2 , especially S || n ( R || ). The actual behavior depends on the flow conditions, in agreement with the results of Ref. 12 . The magnitudes of the structure functions are much lower than for the coflow. At the lower temperature T = 1 65K, . S || n ( R ⊥ ) has an overlap with the corresponding structure function in the coflow a large scales, which disappears with increasing temperature. The two types of the structure functions in the counterflow have significant difference in magnitude, with S || n ( R || ) being strongly suppressed. As it was suggested in Sec.I B 2, these structure functions do not quantitatively reflect the corresponding energy spectra, however the qualitative difference should be observable experimentally.

The influence of the coupling strength on the behavior of the structure functions is illustrated in Fig. 9 for T = 1 65K. Here, in addition to the weak coupling Ω = 1 . and the strong coupling Ω = 20 we consider also an intermediate coupling strength Ω = 3 4. . The structure functions for the classical turbulence are included for comparison. The general form is similar for all values of Ω, with all the structure functions in the counterflow being strongly suppressed compared to classical turbulence, especially S || j ( R || ). Note that at this temperature the structure functions of the normal fluid are more suppressed for weaker coupling, in accordance with the energy spectra in Fig. 1(a). For the transverse S || j ( R ⊥ ) the difference between the two fluid components is relatively small and the influence of the coupling strength is weak. This is consistent with the 2D energy spectra, shown in Fig. 3: the energy spectra in the transverse direction k ⊥ are weakly influenced by the mutual friction.

Additional information may be obtained from analysis of the flatness F || j ( R ) = P || j ( R / S ) [ || j ( R )] 2 , where P || j ( R ) = 〈 ( δ R u x j ) 4 ) 〉 is the forth-order structure function. In Fig. 10 we compare F || n ( R ⊥ ) and F || n ( R || ) with the flatness in the classical turbulence. In the transverse direction, F || n ( R ⊥ ) in the counterflow is growing towards small scales faster than in the classical turbulence at large scales. This indicate a moderate enhancement of intermittency at intermediate scales in this direction, in agreement with experimental results of Ref. 13 . In the longitudinal direction, similar to the structure functions, the flatness F || n ( R || ) is almost constant in the counterflow, leading to a much stronger discrepancy between the transverse and the longitudinal components that in the classical turbulence. This constant value reflects the behavior of the structure functions S || n ( R || ) ∝ R 2 || and P || n ( R || ) ∝ R 4 || over a wide range of scales that is a consequence of the energy spectra ⊥ E j ( k ‖ ) that fall off faster than k -3 .

The second difference structure functions ˜ S || j ( R ) (12b) are expected to better reflect the underlying spectra, at least for ‖ E k ( ⊥ ), since they have wider windows of locality up to k -5 . The steeper ⊥ E k ( ‖ ) result in S || ( R ‖ ) behaving as R 4 in most of the available range of scales. In principle, these exponents can fall within the windows of locality for the structure functions of the third difference (up to x ≤ 7) and of the fourth difference (up to x ≤ 9 0). . We do not discuss these objects due to the increasing difficulties in their measurements.

Having in mind possible experiments, we compare in Fig. 11 two types of the structure functions for the normal fluid in the counterflow. To allow a meaningful comparison we plot them normalized by the values at the largest R = R max and compensated by the corresponding vis-

FIG. 11: Comparison of the normalized compensated velocity structure functions † S || n ( R ⊥ ) ( /square ) and † S || n ( R ‖ ) ( © ) with † ˜ S || n ( R ⊥ )( /triangleright ) and † ˜ S || n ( R ‖ ) ( ♦ ) for the normal-fluid component in the counterflow.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-325625e4-9903-49f8-b059-35e6f8dc8046.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-173ed67d-8f93-4d9b-8cc9-fdf6259af74c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01144v1-figure-7c303f2e-03ba-4032-9c97-f44d7eaec6e7.jpg)

cous scaling

$$\begin{array} { r c l } \dag S _ { j } ^ { | | } ( R ) & = & S _ { j } ^ { | | } ( R ) / S _ { j } ^ { | | } ( R _ { \max } ) R ^ { - 2 } \,, \quad & ( 2 6 ) \quad \text{a n} \\ \dag \tilde { S } _ { j } ^ { | | } ( R ) & = & \tilde { S } _ { j } ^ { | | } ( R ) / \tilde { S } _ { j } ^ { | | } ( R _ { \max } ) R ^ { - 4 } \,. \quad & \text{tion} \\ \quad & \text{flow} \\ \end{array}$$

Indeed, the transition to the viscous behavior (the horizontal lines at small scales) occurs at smaller R for ˜ S || j ( R ) (marked by triangles and diamonds) than for S || j ( R ) (marked by squares and circles). As expected, the range of the condition-dependent apparent scaling at large scales also increases. In addition, the difference in the amplitudes of the structure functions in the longitudinal and transverse directions is much larger for ˜ S || ( R ), hopefully allowing more accurate detection of the anisotropy.

gated across the direction of the counterflow velocity. At the same time, most of the flow energy is confined to a narrow wavenumber plane, orthogonal to this direction, leading to a flow which is smooth along the counterflow direction and turbulent across it. Unlike rotational and atmospheric turbulence with stable stratification, in counterflow turbulence the streamwise velocity component plays the dominant role. This effect is progressively stronger at smaller scales and at higher temperatures. At low temperatures, the milder gradual increase of the small scale anisotropy is due to the smaller fraction of normal fluid and consequently weaker decorrelation. The structure functions of this anisotropic, non-scale invariant turbulent flow, do not allow to extract the quantitative information about the energy distribution over scales, but are expected to reveal strong differences between the directions along and orthogonal to the counterflow velocity.

## III. CONCLUSIONS

Both the theoretical considerations and the results of the numerical simulations presented indicate strong anisotropy in the energy distribution in counterflow turbulence. This is basically due to an angular dependence of the energy dissipation caused by the mutual friction force. It tends to suppress the velocity fluctuations elon-

1 R. J. Donnelly, Quantized Vortices in Hellium II (Cambridge 3 University Press, Cambridge, 1991).

2 Quantized Vortex Dynamics and Superfluid Turbulence , edited by C.F. Barenghi, R.J. Donnelly and W.F. Vinen, Lecture Notes in Physics 571 (Springer-Verlag, Berlin, 2001)

3 W. F. Vinen and J. J. Niemela, Quantum turbulence. J. Low Temp. Phys. 128 , 167 (2002).

4 R. P.Feynman, Application of quantum mechanics to liquid helium. Progress in Low Temperature Physics 1 , 17 (1955).

5 H. E. Hall and W. F. Vinen, The rotation of liquid helium

## Acknowledgments

L.B. Thanks Michele Buzzicotti for the data analysis and visualization. GS thanks AtMath collaboration at University of Helsinki.

II. I. Experiments on the propagation of second sound in uniformly rotating helium II. Proc. Roy. Soc. A 238 , 204 (1956).Donnel

6 I.L. Bekarevich, and I.M. Khalatnikov, Phenomenological Derivation of the Equations of Vortex Motion in He II, Sov. Phys. JETP 13 , 643 (1961).

7 W. F. Vinen, Mutual friction in a heat current in liquid helium II I. Experiments on steady heat currents, Proc. R. Soc. 240 , 114 (1957); Mutual friction in a heat current in liquid helium II. II. Experiments on transient effects, 240 , 128 (1957); Mutual friction in a heat current in liquid he-

- lium II III. Theory of the mutual friction, 242 , 493 (1957); Restricted access Mutual friction in a heat current in liquid helium. II. IV. Critical heat currents in wide channels, 243 , 400 (1958).
- 8 R. N. Hills and P. H. Roberts, Superfluid mechanics for a high density of vortex lines, Arch. Ration. Mech. Anal. 66 , 43 (1977).
- 9 K. W. Schwarz, Three-dimensional vortex dynamics in superfluid 4 He: Homogeneous superfluid turbulence, Phys. Rev. B 38 , 2398 (1988).
- 10 L. Skrbek and K. R. Sreenivasan, in Ten Chapters in Turbulence , edited by P. A. Davidson, Y. Kaneda, and K. R. Sreenivasan (Cambridge University Press, Cambridge, 2013), pp. 405-437.
- 11 A. Marakov, J. Gao, W. Guo, S. W. Van Sciver, G. G. Ihas, D. N. McKinsey, and W. F. Vinen. Visualization of the normal-fluid turbulence in counterflowing superfluid 4 He, Phys. Rev. B, 91 094503. (2015).
- 12 J. Gao, E. Varga, W. Guo and W. F. Vinen, Energy spectrum of thermal counterflow turbulence in superfluid Helium-4, Phys. Rev. B 96 , 094511 (2017).
- 13 S. Bao, W. Guo, V. S. L'vov, A. Pomyalov, Phys. Rev. B 98 , 174509 (2018).
- 14 M. La Mantia, L. Skrbek, Europhys. Lett. 105, 46002 (2014).
- 15 M. La Mantia, P. ˘ van˘ara, S c D. Duda, and L. Skrbek, Small-scale universality of particle dynamics in quantum turbulence, Phys. Rev B 94 , 184512 (2016).
- 16 M. La Mantia, Particle dynamics in wall-bounded thermal counterflow of superfluid helium, Physics of Fluids 29 , 065102 (2017);
- 17 D. Khomenko, V. S. L'vov, A. Pomyalov, and I. Procaccia, Counterflow induced decoupling in superfluid Turbulence. Phys. Rev. B 93 , 014516 (2016).
- 18 V. S. L'vov and A. Pomyalov, A theory of counterflow velocity dependence of superfluid 4 He turbulence statistics, Phys. Rev. B, 97 , 214513 (2018).
- 19 L. Biferale; D. Khomenko; V. L'vov; A. Pomyalov; I. Procaccia; G. Sahoo (2019). Superfluid Helium in ThreeDimensional Counterflow Differs Strongly from Classical Flows: Anisotropy on Small Scales. Physical Review Letters, 122 ,144501
- 20 L. Biferale, D. Khomenko, V.S. L'vov, A. Pomyalov, I. Procaccia, and G. Sahoo, Turbulent statistics and intermittency enhancement in coflowing superfluid 4 He, Phys. Rev.Fluids 3 , 024605 (2018).
- 21 L. Bou´ e, V.S. L'vov, A. Pomyalov, and I. Procaccia, Energy spectra of superfluid turbulence in 3 He, Phys. Rev. B 85 , 104502 (2012).
- 22 L. Biferale, D. Khomenko, V. L'vov, A. Pomyalov, I. Procaccia and G. Sahoo, Local and non-local energy spectra of superfuid 3 He turbulence, Phys. Rev. B. 95 , 184510 (2017).
- 23 A. Kumar, M. K. Verma and J. Sukhatme, Phenomenology of two-dimensional stably stratified turbulence under largescale forcing, J. of Turbulence, 18 , 219(2017).
- 24 A. Alexakis, L. Biferale, Phys. Rep. 767 769 -,1 (2018).
- 25 L. Biferale and I. Procaccia, Anisotropy in Turbulent Flows and in Turbulent Transport, Phys. Rep. 414 43, (2005)
- 26 L. Biferale, F. Bonaccorso, I.M. Mazzitelli, M.A.T. van Hinsberg, A.S. Lanotte, S. Musacchio, P. Perlekar, and F. Toschi. Phys. Rev. X 6 , 041036 (2016).
- 27 B. Gallet, A. Campagne, P.-P. Cortet, and F. Moisy, Phys. Fluids 26 , 035108 (2014).
- 28 B. Gallet, J. Fluid Mech. 783 , 412 (2015).
- 29 R. J. Donnelly, C. F. Barenghi , The Observed Properties of Liquid Helium at the Saturated Vapor Pressure, J. Phys. Chem. Ref. Data 27, 1217(1998).
- 30 C. F. Barenghi, V. S. L'vov, and P.-E. Roche, Experimental, numerical, and analytical velocity spectra in turbulent quantum fluid, Proc Natl Acad Sci USA 111 , 4683 (2014).
- 31 L. Skrbek, K.R. Sreenivasan Developed quantum turbulence and its decay. Phys Fluids 24 , 011301 (2012).
- 32 E. Rusaouen, B. Chabaud, J. Salort, Philippe-E. Roche. Intermittency of quantum turbulence with superfluid fractions from 0% to 96%. Physics of Fluids 29 , 105108 (2017).
- 33 S. Babuin, V.S. L'vov, A. Pomyalov, L. Skrbek, E. Varga, Coexistence and interplay of quantum and classic al turbulence in superfluid He-4: Phys. Rev. B, 94 , 174504 (2016).
- 34 L. Boue, V.S. L'vov, Y. Nagar, S.V. Nazarenko, A. Pomyalov, I. Procaccia, Energy and vorticity spectra in turbulent superfluid He-4 from T = 0 to T λ . Phys. Rev. B. 91 , 144501, (2015).
- 35 S. B. Pope, Turbulent Flows (Cambridge University Press, Cambridge, 2000).
- 36 V. S. L'vov, S. V. Nazarenko and G. E. Volovik, Energy spectra of developed superfluid turbulence, JETP Letters, 80 , 535 (2004).
- 37 J. Salort, B. Chabaud, E. Leveque, and P.E. Roche. Investigation of intermittency in superfluid turbulence. Jour. Phys. : Conf. Series, 318 , (2011).
- 38 J. Salort, P.E. Roche and Leveque, Mesoscale equipartition of kinetic energy in quantum turbulence EPL, 94 24001 (2011).
- 39 V. S. L'vov, S. V. Nazarenko, and L. Skrbek, Energy Spectra of Developed Turbulence in Helium Superfluids , Journal of Low Temperature Physics, 145 , 125 (2006).
- 40 C. J. Gorter and J. H. Mellink, Physica 15 , 285 (1949).
- 41 L. Kondaurova; V.S. L'vov; A.Pomyalov; I. Procaccia. Structure of a quantum vortex tangle in He-4 counterflow turbulence. Physical Review B. 89 014502 (2014).
- 42 D. N. McKinsey, C. R. Brome, J. S. Butterworth, S. N. Dzhosyuk, P. R. Huffman, C. E. H. Mattoni, J. M. Doyle, R. Golub, and K. Habicht, Habicht, Radiative decay of the metastable He 2 ( a 3 Σ ) molecule in liquid + u helium, Phys. Rev. A 59 , 200 (1999).
- 43 L. Biferale, M. Cencini, A. Lanotte, and D. Vergini, Inverse velocity statistics in two dimensions, Phys. Fluids 15 , 1012 (2003).
- 44 G.L. Eyink, Exctact results on stationary turbulence in 2D: Consequences of vorticity conservation, Physica D ( 91), 97 (1991).
- 45 V. S. L'vov and A. Pomyalov, Statistics of Quantum Turbulence in Superfluid He, J. Low Temp Phys 187 , 497 (2017).
- 46 V.S. L'vov and I. Procaccia, The universal scaling exponents of anisotropy in turbulence and their measurement. Physics of Fluids 8 , 2565 (1996).
- 47 I Arad, V.S. L'vov and I. Procaccia, Correlation functions in isotropic and anisotropic turbulence: The role of the symmetry group, Phys. Rev. E 59 , 6753 (1999).
- 48 I Arad, V.S. L'vov, and I. Procaccia, Anomalous scaling in anisotropic turbulence, Physica A. 288 , 280 (2000).
- 49 V.S. L'vov, I. Procaccia and V. Tiberkevich, Scaling exponents in anisotropic hydrodynamic turbulence, Phys. Rev. E 67 , 026312 (2003).
- 50 Please notice that in the Ref. 19 there is a typo referring to these data as u 2 s x , and u 2 s y , .
L. Biferale D. Khomenko V. S. L'vov A. Pomyalov I. Procaccia G. Sahoo
Strong anisotropy of superfluid $^4$He counterflow turbulence
We report on a combined theoretical and numerical study of counterflow
turbulence in superfluid $^{4}$He in a wide range of parameters. The energy
spectra of the velocity fluctuations of both the normal-fluid and superfluid
components are strongly anisotropic. The angular dependence of the correlation
between velocity fluctuations of the two components plays the key role. A
selective energy dissipation intensifies as scales decrease, with the
streamwise velocity fluctuations becoming dominant. Most of the flow energy is
concentrated in a wavevector plane which is orthogonal to the direction of the
counterflow. The phenomenon becomes more prominent at higher temperatures as
the coupling between the components depends on the temperature and the
direction with respect to the counterflow velocity.