No Arabic abstract
We study the evolution of kinetic and magnetic energy spectra in magnetohydrodynamic flows in the presence of strong cross helicity. For forced turbulence, we find weak inverse transfer of kinetic energy toward the smallest wavenumber. This is plausibly explained by the finiteness of scale separation between the injection wavenumber and the smallest wavenumber of the domain, which here is a factor of 15. In the decaying case, there is a slight increase at the smallest wavenumber, which is probably explained by the dominance of kinetic energy over magnetic energy at the smallest wavenumbers. Within a range of wavenumbers covering almost an order of magnitude the decay is purely exponential, which is argued to be a consequence of a suppression of nonlinearity due to the presence of strong cross helicity.
Hydromagnetic turbulence produced during phase transitions in the early universe can be a powerful source of stochastic gravitational waves (GWs). GWs can be modelled by the linearised spatial part of the Einstein equations sourced by the Reynolds and Maxwell stresses. We have implemented two different GW solvers into the {sc Pencil Code} -- a code which uses a third order timestep and sixth order finite differences. Using direct numerical integration of the GW equations, we study the appearance of a numerical degradation of the GW amplitude at the highest wavenumbers, which depends on the length of the timestep -- even when the Courant--Friedrichs--Lewy condition is ten times below the stability limit. This degradation leads to a numerical error, which is found to scale with the third power of the timestep. A similar degradation is not seen in the magnetic and velocity fields. To mitigate numerical degradation effects, we alternatively use the exact solution of the GW equations under the assumption that the source is constant between subsequent timesteps. This allows us to use a much longer timestep, which cuts the computational cost by a factor of about ten.
Numerical simulations are made for forced turbulence at a sequence of increasing values of Reynolds number, R, keeping fixed a strongly stable, volume-mean density stratification. At smaller values of R, the turbulent velocity is mainly horizontal, and the momentum balance is approximately cyclostrophic and hydrostatic. This is a regime dominated by so-called pancake vortices, with only a weak excitation of internal gravity waves and large values of the local Richardson number, Ri, everywhere. At higher values of R there are successive transitions to (a) overturning motions with local reversals in the density stratification and small or negative values of Ri; (b) growth of a horizontally uniform vertical shear flow component; and (c) growth of a large-scale vertical flow component. Throughout these transitions, pancake vortices continue to dominate the large-scale part of the turbulence, and the gravity wave component remains weak except at small scales.
The previously reported non-equilibrium dissipation law is investigated in turbulent flows generated by various regular and fractal square grids. The flows are documented in terms of various turbulent profiles which reveal their differences. In spite of significant inhomogeneity and anisotropy differences, the new non-equilibrium dissipation law is observed in all these flows. Various transverse and longitudinal integral scales are measured and used to define the dissipation coefficient $C_{varepsilon}$. It is found that the new non-equilibrium dissipation law is not an artefact of a particular choice of the integral scale and that the usual equilibrium dissipation law can actually coexist with the non-equilibrium law in different regions of the same flow.
We investigate non-equilibrium turbulence where the non-dimensionalised dissipation coefficient $C_{varepsilon}$ scales as $C_{varepsilon} sim Re_{M}^{m}/Re_{ell}^{n}$ with $mapprox 1 approx n$ ($Re_M$ and $Re_{ell}$ are global/inlet and local Reynolds numbers respectively) by measuring the downstream evolution of the scale-by-scale energy transfer, dissipation, advection, production and transport in the lee of a square-mesh grid and compare with a region of equilibrium turbulence (i.e. where $C_{varepsilon}approx mathrm{constant}$). These are the main terms of the inhomogeneous, anisotropic version of the von K{a}rm{a}n-Howarth-Monin equation. It is shown in the grid-generated turbulence studied here that, even in the presence of non-negligible turbulence production and transport, production and transport are large-scale phenomena that do not contribute to the scale-by-scale balance for scales smaller than about a third of the integral-length scale, $ell$, and therefore do not affect the energy transfer to the small-scales. In both the non-equilibrium and the equilibrium decay regions, the peak of the scale-by-scale energy transfer scales as $(overline{u^2})^{3/2}/ell$ ($overline{u^2}$ is the variance of the longitudinal fluctuating velocity). In the non-equilibrium case this scaling implies an imbalance between the energy transfer to the small scales and the dissipation. This imbalance is reflected on the small-scale advection which becomes larger in proportion to the maximum energy transfer as the turbulence decays whereas it stays proportionally constant in the further downstream equilibrium region where $C_{varepsilon} approx mathrm{constant}$ even though $Re_{ell}$ is lower.
We present results from an ensemble of 50 runs of two-dimensional hydrodynamic turbulence with spatial resolution of 2048^2 grid points, and from an ensemble of 10 runs with 4096^2 grid points. All runs in each ensemble have random initial conditions with same initial integral scale, energy, enstrophy, and Reynolds number. When both ensemble- and time-averaged, inverse energy cascade behavior is observed, even in the absence of external mechanical forcing: the energy spectrum at scales larger than the characteristic scale of the flow follows a k^(-5/3) law, with negative flux, together with a k^(-3) law at smaller scales, and a positive flux of enstrophy. The source of energy for this behavior comes from the modal energy around the energy containing scale at t=0. The results shed some light into connections between decaying and forced turbulence, and into recent controversies in experimental studies of two-dimensional and magnetohydrodynamic turbulent flows.