BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alexandros Alexakis

Alexandros Alexakis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynphysics.plasm-phnlin.CDastro-ph.SRcond-mat.stat-mechphysics.ao-phastro-ph.EPmath-phmath.MPnlin.AO

Catalog footprint

What is connected

38works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2026arXiv

Ensemble of Fixed Points in Multi-branch Shell Models of Turbulent Cascades

Stationary solutions of a shell model of turbulence defined on a dyadic tree topology are studied. Each node's amplitude is expressed as the product of amplitude multipliers associated with its ancestors, providing a recursive representation of the cascade process. A geometrical rule governs the tree growth, and we prove the existence of a continuum of fixed points, including the Kolmogorov solution, that sustain a strictly forward energy cascade. Sampling along randomly chosen branches defines a homogeneous Markov chain, enabling a stochastic characterization of extended self-similarity and intermittency through the spectral properties of the associated Feynman-Kac operators. Numerical simulations confirm the theoretical predictions, showing that multi-branch shell models offer a minimal yet physically rich framework for exploring the complexity of nonlinear energy transfer across scales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bistability of the large-scale dynamics in quasi-two-dimensional turbulence

In many geophysical and astrophysical flows, suppression of fluctuations along one direction of the flow drives a quasi-2D upscale flux of kinetic energy, leading to the formation of strong vortex condensates at the largest scales. Recent studies have shown that the transition towards this condensate state is hysteretic, giving rise to a limited bistable range in which both the condensate state as well as the regular 3D state can exist at the same parameter values. In this work, we use direct numerical simulations of thin-layer flow to investigate whether this bistable range survives as the domain size and turbulence intensity are increased. By studying the time scales at which rare transitions occur from one state into the other, we find that the bistable range grows as the box size and/or Reynolds number Re are increased, showing that the bistability is neither a finite-size nor a finite-Re effect. We furthermore predict a crossover from a bimodal regime at low box size, low Re to a regime of pure hysteresis at high box size, high Re, in which any transition from one state to the other is prohibited at any finite time scale.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Inducing intermittency in the inverse cascade of two dimensional turbulence by a fractal forcing

We demonstrate that like in the forward cascade of three dimensional turbulence that displays intermittency (lack of self-similarity) due to the concentration of energy dissipation in a small set of fractal dimension less than three, the inverse cascade of two-dimensional turbulence can also display lack of self-similarity and intermittency if the energy injection is constrained in a fractal set of dimension less than two. A series of numerical simulations of two dimensional turbulence are examined, using different forcing functions of the same forcing length-scale but different fractal dimension $D$ that varies from the classical $D=2$ case to the point vortex case $D=0$. It is shown that as the fractal dimension of the forcing is decreased from $D=2$, the self-similarity is lost and intermittency appears, with the scaling of the different structure functions $\langle |δu_\|^p| \rangle\propto r^{ζ_p}$ differs from the dimensional analysis prediction $ζ_p=p/3$. The present model thus provides a unique example that intermittency is controlled and can thus shed light and provide test beds for multi-fractal models of turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local fluxes in MHD Turbulence

Using highly resolved direct numerical simulations we examine the statistical properties of the local energy flux rate $Π_\ell(x)$ towards small scales for three isotropic turbulent magnetohydrodynamic flows, which differ in strength and structure of the magnetic field. We analyse the cascade process both in the kinetic and magnetic energy, disentangling the different flux contributions to the overall energy dynamics. The results show that the probability distribution of the local energy flux develops long tails related to extreme events, similar to the hydrodynamic case. The different terms of the energy flux display different properties and show sensitivity on the type of the flow examined. We further examine the joint pdf between the local energy flux and the gradients of the involved fields. The results point out a correlation with the magnetic field gradients, showing however a dispersion much stronger than what is observed in hydrodynamic flows. Finally, it is also shown that the local energy flux shows some dependence on the local amplitude of the magnetic field. The present results have implications for subgrid scale models that we discuss.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

$λ$-Navier-Stokes turbulence

We investigate numerically the model proposed in Sahoo et al [Phys. Rev. Lett. 118, 164501, (2017)] where a parameter $λ$ is introduced in the Navier-Stokes equations such that the weight of homochiral to heterochiral interactions is varied while preserving all original scaling symmetries and inviscid invariants. Decreasing the value of $λ$ leads to a change in the direction of the energy cascade at a critical value $λ_c \sim 0.3$. In this work, we perform numerical simulations at varying $λ$ in the forward energy cascade range and at changing the Reynolds number $\mathrm{Re}$. We show that for a fixed injection rate, as $λ\to λ_c$, the kinetic energy diverges with a scaling law $\mathcal{E} \propto (λ-λ_c)^{-2/3}$. The energy spectrum is shown to display a larger bottleneck as $λ$ is decreased. The forward heterochiral flux and the inverse homochiral flux both increase in amplitude as $λ_c$ is approached while keeping their difference fixed and equal to the injection rate. As a result, very close to $λ_c$ a stationary state is reached where the two opposite fluxes are of much higher amplitude than the mean flux and large fluctuations are observed. Furthermore, we show that intermittency as $λ_c$ is approached is reduced. The possibility of obtaining a statistical description of regular Navier-Stokes turbulence as an expansion around this newly found critical point is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Do extreme events trigger turbulence decay? - a numerical study of turbulence decay time in pipe flows

Turbulence locally created in laminar pipe flows shows sudden decay or splitting after a stochastic waiting time. In laboratory experiments, the mean waiting time was observed to increase double-exponentially as the Reynolds number (Re) approaches its critical value. To understand the origin of this double-exponential increase, we perform many pipe flow direct numerical simulations of the Navier-Stokes equations, and measure the cumulative histogram of the maximum axial vorticity field over the pipe (turbulence intensity). In the domain where the turbulence intensity is not small, we observe that the histogram is well-approximated by the Gumbel extreme-value distribution. The smallest turbulence intensity in this domain roughly corresponds to the transition value between the locally stable turbulence and a meta-stable (edge) state. Studying the Re dependence of the fitting parameters in this distribution, we derive that the time scale of the transition between these two states increases double-exponentially as Re approaches its critical value. On the contrary, in smaller turbulence intensities below this domain, we observe that the distribution is not sensible to the change of Re. This means that the decay time of the meta-stable state (to the laminar state) is stochastic but Re-independent in average. Our observation suggests that the conjecture made by Goldenfeld et al. to derive the double-exponential increase of turbulence decay time is approximately satisfied in the range of Re we studied. We also discuss using another extreme-value distribution, Fréchet distribution, instead of the Gumbel distribution to approximate the histogram of the turbulence intensify, which reveals interesting perspectives.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Energy cascades in rapidly rotating and stratified turbulence within elongated domains

We study forced, rapildy rotating and stably stratified turbulence in an elongated domain using an asymptotic expansion at simultaneously low Rossby number $Ro\ll1$ and large domain height compared to the energy injection scale, $h=H/\ell_{in}\gg1$. The resulting equations depend on the parameter $λ=(h Ro )^{-1}$ and the Froude number $Fr$. An extensive set of direct numerical simulations (DNS) is performed to explore the parameter space $(λ,Fr)$. We show that a forward energy cascade occurs in one region of this space, and a split energy cascade outside it. At weak stratification (large $Fr$), an inverse cascade is observed for sufficiently large $λ$. At strong stratification (small $Fr$) the flow becomes approximately hydrostatic and an inverse cascade is always observed. For both weak and strong stratification, we present theoretical arguments supporting the observed energy cascade phenomenology. Our results shed light on an asymptotic region in the phase diagram of rotating and stratified turbulence, which is difficult to attain by brute-force DNS.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Geometric microcanonical theory of two-dimensional Truncated Euler flows

This paper presents a geometric microcanonical ensemble perspective on two-dimensional Truncated Euler flows, which contain a finite number of (Fourier) modes and conserve energy and enstrophy. We explicitly perform phase space volume integrals over shells of constant energy and enstrophy. Two applications are considered. In a first part, we determine the average energy spectrum for highly condensed flow configurations and show that the result is consistent with Kraichnan's canonical ensemble description, despite the fact that no thermodynamic limit is invoked. In a second part, we compute the probability density for the largest-scale mode of a free-slip flow in a square, which displays reversals. We test the results against numerical simulations of a minimal model and find excellent agreement with the microcanonical theory, unlike the canonical theory, which fails to describe the bimodal statistics. This article is part of the theme issue "Mathematical problems in physical fluid dynamics".

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Symmetry breaking in a turbulent environment

In this work we investigate symmetry breaking in the presence of a turbulent environment. The transition from a symmetric state to a symmetry-breaking state is demonstrated using two examples: (i) the transition of a two-dimensional flow to a three dimensional flow as the fluid layer thickness is varied and (ii) the dynamo instability in a thin layer flow as the magnetic Reynolds number is varied. We show that these examples have similar critical exponents that differ from the mean-field predictions. The critical behavior can be related to the multiplicative nature of the fluctuations and can be predicted in certain limits using results from the statistical properties of random interfaces. Our results indicate the possibility of existence of a new class of out-of-equilibrium phase transition controlled by the multiplicative noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Critical transition in fast-rotating turbulence within highly elongated domains

We study rapidly rotating turbulent flows in a highly elongated domain using an asymptotic expansion at simultaneously low Rossby number $Ro\ll 1$ and large domain height compared to the energy injection scale, $h=H/\ell_{in}\gg 1$. We solve the resulting equations using an extensive set of direct numerical simulations for different parameter regimes. As the parameter $λ= (h Ro)^{-1}$ is increased beyond a threshold $λ_c$, a transition is observed from a state without an inverse energy cascade to a state with an inverse energy cascade. } For large Reynolds number and large horizontal box size, we provide evidence for criticality of the transition in terms of the large-scale energy dissipation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the local energy flux of turbulent flows

We investigate the local energy flux rate $Π_\ell(\bf x)$ towards small scales in isotropic turbulent flows using direct numerical simulations and applying different low-pass filters. Two different filters are examined, a sharp Fourier filter and a Gaussian filter. The probability density function (pdf) of the local energy flux is calculated for the different filters and for different filtering scales. It is shown that the local energy flux is a largely fluctuating quantity taking both negative and positive values and this is more pronounced for the sharp filter. The variance, the skewness and the kurtosis of these fluctuations are shown to increase as the filtering scale is decreased. Furthermore we calculate the joint pdf of $Π_\ell(\bf x)$ with the local filtered strain rate $S_\ell$ and the enstrophy $Ω_\ell$. The flux shows a good correlation with the strain but not with the enstrophy. It is shown that its conditional mean value scales like $<{Π_\ell}>_S \propto \ell^2 S_\ell^{3} $ in support to the Smagorinsky eddy viscosity model. Nonetheless strong fluctuations exist around this value that also need to be modeled. We discuss the implications of our results for subgrid scale models, and propose new modelling directions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Role of the forcing dimensionality in thin-layer turbulent energy cascades

We investigate the transition from forward to inverse energy cascade in turbulent flows in thin layers, varying the functional form of the forcing and the thickness of the layer. We show that, as the forcing function becomes more three-dimensional, the inverse cascade is suppressed and the critical height hc, where the transition occurs, is decreased. We study the dependence of this critical height on a parameter r which measures the dimensionality of the forcing and thus construct a phase space diagram in the parameter space r-h. We discuss the effect of Reynolds number and domain size.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Energy fluxes in quasi-equilibrium flows

We examine the relation between the absolute equilibrium state of the spectrally truncated Euler equations (TEE) predicted by Kraichnan (1973) to the forced and dissipated flows of the spectrally truncated Navier-Stokes (TNS) equations. In both of these idealized systems a finite number of Fourier modes is kept contained inside a sphere of radius $k_{max}$ but while the first conserves energy in the second energy is injected by a body-force $\bf{f}$ and dissipated by the viscosity $ν$. For the TNS system stochastically forced with energy injection rate $\mathcal{I}_\mathcal{E}$ we show, using an asymptotic expansion of the Fokker-Planck equation, that in the limit of small $k_{max}η$ (where $η=(ν^3/\mathcal{I}_\mathcal{E})^{1/4}$ the Kolmogorov lengthscale) the flow approaches the absolute equilibrium solution of Kraichnan with such an effective "temperature" so that there is a balance between the energy injection and the energy dissipation rate. Numerical simulations verify the predictions of the model for small values of $k_{max}η$. For intermediate values of $k_{max}η$ a transition from the quasi-equilibrium "thermal" state to Kolmogorov turbulence is observed. If the forcing is applied at small scales and the dissipation acts only at large scales then the equipartition spectrum appears at all scales for all values of $ν$. In both cases a finite forward or inverse flux is present even for the cases where the flow is close to the equilibrium state solutions. However, unlike the classical turbulence where an energy cascade develops with a mean energy flux that is large compared to its fluctuations, the quasi-equilibrium state has a mean flux of energy that is subdominant to the large flux fluctuations observed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Turbulent cascade, bottleneck and thermalized spectrum in hyperviscous flows

In many simulations of turbulent flows the viscous forces $ν\nabla^2 {\bf u}$ are replaced by a hyper-viscous term $-ν_p(-\nabla^2)^{p}{\bf u}$ in order to suppress the effect of viscosity at the large scales. In this work we examine the effect of hyper-viscosity on decaying turbulence for values of $p$ ranging from $p=1$ (regular viscosity) up to $p=100$. Our study is based on direct numerical simulations of the Taylor-Green vortex for resolutions from $512^3$ to $2048^3$. Our results demonstrate that the evolution of the total energy $E$ and the energy dissipation $ε$ remain almost unaffected by the order of the hyper-viscosity used. However, as the order of the hyper-viscosity is increased ,the energy spectrum develops a more pronounced bottleneck that contaminates the inertial range. At the largest values of $p$ examined, the spectrum at the bottleneck range has a positive power-law behavior $E(k)\propto k^α$ with the power-law exponent $α$ approaching the value obtained in flows at thermal equilibrium $α=2$. This agrees with the prediction of Frisch et al. [Phys. Rev. Lett. 101, 144501 (2008)] who suggested that at high values of $p$, the flow should behave like the truncated Euler equations (TEE). Nonetheless, despite the thermalization of the spectrum, the flow retains a finite dissipation rate up to the examined order, which disagrees with the predictions of the TEE system implying suppression of energy dissipation. We reconcile the two apparently contradictory results, predicting the value of $p$ for which the hyper-viscous Navier-Stokes goes over to the TEE system and we discuss why thermalization appears at smaller values of $p$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fluctuations of electrical conductivity: a new source for astrophysical magnetic fields

We consider the generation of magnetic field by the flow of a fluid for which the electrical conductivity is nonuniform. A new amplification mechanism is found which leads to dynamo action for flows much simpler than those considered so far. In particular, the fluctuations of the electrical conductivity provide a way to bypass anti-dynamo theorems. For astrophysical objects, we show through three-dimensional global numerical simulations that the temperature-driven fluctuations of the electrical conductivity can amplify an otherwise decaying large scale equatorial dipolar field. This effect could play a role for the generation of the unusually tilted magnetic field of the iced giants Neptune and Uranus.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large-scale instabilities of helical flows

Large-scale hydrodynamic instabilities of periodic helical flows are investigated using $3$D Floquet numerical computations. A minimal three-modes analytical model that reproduce and explains some of the full Floquet results is derived. The growth-rate $σ$ of the most unstable modes (at small scale, low Reynolds number $Re$ and small wavenumber $q$) is found to scale differently in the presence or absence of anisotropic kinetic alpha (\AKA{}) effect. When an $AKA$ effect is present the scaling $σ\propto q\; Re\,$ predicted by the $AKA$ effect theory [U. Frisch, Z. S. She, and P. L. Sulem, Physica D: Nonlinear Phenomena 28, 382 (1987)] is recovered for $Re\ll 1$ as expected (with most of the energy of the unstable mode concentrated in the large scales). However, as $Re$ increases, the growth-rate is found to saturate and most of the energy is found at small scales. In the absence of \AKA{} effect, it is found that flows can still have large-scale instabilities, but with a negative eddy-viscosity scaling $σ\propto ν(b Re^2-1) q^2$. The instability appears only above a critical value of the Reynolds number $Re^{c}$. For values of $Re$ above a second critical value $Re^{c}_{S}$ beyond which small-scale instabilities are present, the growth-rate becomes independent of $q$ and the energy of the perturbation at large scales decreases with scale separation. A simple two-modes model is derived that well describes the behaviors of energy concentration and growth-rates of various unstable flows. In the non-linear regime (at moderate values of $Re$) and in the presence of scale separation, the forcing scale and the largest scales of the system are found to be the most dominant energetically.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The fate of alpha dynamos at large $Rm$

At the heart of today's solar magnetic field evolution models lies the alpha dynamo description. In this work, we investigate the fate of alpha-dynamos as the magnetic Reynolds number $Rm$ is increased. Using Floquet theory, we are able to precisely quantify mean field effects like the alpha and beta effect (i) by rigorously distinguishing dynamo modes that involve large scale components from the ones that only involve small scales, and by (ii) providing a way to investigate arbitrary large scale separations with minimal computational cost. We apply this framework to helical and non-helical flows as well as to random flows with short correlation time. Our results determine that the alpha-description is valid for $Rm$ smaller than a critical value $Rm_c$ at which small scale dynamo instability starts. When $Rm$ is above $Rm_c$ the dynamo ceases to follow the mean field description and the growth rate of the large scale modes becomes independent of the scale separation while the energy in the large scale modes is inversely proportional to the square of the scale separation. The results in this second regime do not depend on the presence of helicity. Thus alpha-type modeling for solar and stellar models needs to be reevaluated and new directions for mean field modeling are proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal lengthscale for a turbulent dynamo

We demonstrate that there is an optimal forcing length scale for low Prandtl number dynamo flows, that can significantly reduce the required energy injection rate. The investigation is based on simulations of the induction equation in a periodic box of size $2πL$. The flows considered are turbulent ABC flows forced at different forcing wavenumbers $k_f$ simulated using a subgrid turbulent model. The critical magnetic Reynolds number $Rm_c^T$ decreases as the forcing wavenumber $k_f$ increases from the smallest allowed $k_{min}=1/L$. At large $k_f$ on the other hand, $Rm_c^T$ increases with the forcing wavenumber as $Rm_c^T \propto \sqrt{ k_f}$ in agreement with mean-field scaling prediction. At $k_f L\simeq 4$ an optimal wavenumber is reached where $Rm_c^T$ obtains its minimum value. At this optimal wavenumber $Rm_c^T$ is smaller by more than a factor of ten than the case forced in $k_f=1$. This leads to a reduction of the energy injection rate by three orders of magnitude when compared to the case that the system is forced in the largest scales and thus provides a new strategy for the design of a fully turbulent experimental dynamo.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statistical equilibria of large scales in dissipative hydrodynamic turbulence

We present a numerical study of the statistical properties of three-dimensional dissipative turbulent flows at scales larger than the forcing scale. Our results indicate that the large scale flow can be described to a large degree by the truncated Euler equations with the predictions of the zero flux solutions given by absolute equilibrium theory, both for helical and non-helical flows. Thus, the functional shape of the large scale spectra can be predicted provided that scales sufficiently larger than the forcing length scale but also sufficiently smaller than the box size are examined. Deviations from the predictions of absolute equilibrium are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Transition from 2D HD to 2D MHD turbulence

We investigate the critical transition from an inverse cascade of energy to a forward energy cascade in a two-dimensional magneto-hydrodynamic flow as the ratio of magnetic to mechanical forcing amplitude is varied. It is found that the critical transition is the result of two competing processes. The first process is due to hydrodynamic interactions, cascades the energy to the large scales. The second process couples small scale magnetic fields to large scale flows transferring the energy back to the small scales via a non-local mechanism. At marginality the two cascades are both present and cancel each other. The phase space diagram of the transition is sketched.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the edge of an inverse cascade

We demonstrate that systems with a parameter-controlled inverse cascade can exhibit critical behavior for which at the critical value of the control parameter the inverse cascade stops. In the vicinity of such a critical point standard phenomenological estimates for the energy balance will fail since the energy flux towards large length scales becomes zero. We demonstrate these concepts using the computationally tractable model of two-dimensional magneto-hydrodynamics in a periodic box. In the absence of any external magnetic forcing the system reduces to hydrodynamic fluid turbulence with an inverse energy cascade. In the presence of strong magnetic forcing the system behaves as 2D magneto-hydrodynamic turbulence with forward energy cascade. As the amplitude of the magnetic forcing is varied a critical value is met for which the energy flux towards the large scales becomes zero. Close to this point the energy flux scales as a power law with the departure from the critical point and the normalized amplitude of the fluctuations diverges. Similar behavior is observed for the flux of the square vector potential for which no inverse flux is observed for weak magnetic forcing, while a finite inverse flux is observed for magnetic forcing above the critical point. We conjecture that this behavior is generic for systems of variable inverse cascade.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rotating Taylor-Green Flow

The steady state of a forced Taylor-Green flow is investigated in a rotating frame of reference. The investigation involves the results of 184 numerical simulations for different Reynolds number $\mathrm{Re}_{_{F}}$ and Rossby number $\mathrm{Ro}_{_F}$. The large number of examined runs allows a systematic study that enables the mapping of the different behaviors observed to the parameter space ($\mathrm{Re}_{_{F}},\mathrm{Ro}_{_{F}}$), and the examination of different limiting procedures for approaching the large $\mathrm{Re}_{_{F}}$ small $\mathrm{Ro}_{_F}$ limit. Four distinctly different states were identified: {\it laminar, intermittent bursts, quasi-2D condensates, and weakly rotating turbulence}. These four different states are separated by power-law boundaries $\mathrm{Ro}_{_F} \propto \mathrm{Re}_{_{F}}^{-γ}$ in the small $\mathrm{Ro}_{_F}$ limit. In this limit, the predictions of asymptotic expansions can be directly compared to the results of the direct numerical simulations. While the first order expansion is in good agreement with the results of the linear stability theory, it fails to reproduce the dynamical behavior of the quasi-2D part of the flow in the nonlinear regime, indicating that higher order terms in the expansion need to be taken in to account. The large number of simulations allows also to investigate the scaling that relates the amplitude of the fluctuations with the energy dissipation rate and the control parameters of the system for the different states of the flow. Different scaling was observed for different states of the flow, that are discussed in detail. The present results clearly demonstrate that the limits small Rossby and large Reynolds do not commute and it is important to specify the order in which they are taken.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-organisation and non-linear dynamics in driven magnetohydrodynamic turbulent flows

Magnetohydrodynamic turbulent flows driven by random mechanical and electromagnetic external forces of zero helicities are investigated by means of direct numerical simulations. It is shown that despite the absence of helicities in the forcing, the system is attracted to self-organized helical states that exhibit laminar behaviour despite the large value of the Reynolds numbers examined. We demonstrate that the correlation time of the external forces is controlling the time spent on these states, i.e. for short correlation times the system remains in the turbulent state while as the correlation time is increased the system spends more and more time in the self-organised states. As a result, time averaged statistics can significantly be affected by the time spent on these states. These results have important theoretical implications for the understanding of the suppression of non-linearities in plasma fusion devises as well as in astrophysical observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The role of dissipation in flexural wave turbulence: from experimental spectrum to Kolmogorov-Zakharov spectrum

The Weak Turbulence Theory has been applied to waves in thin elastic plates obeying the Föppl-Von Kármán dynamical equations. Subsequent experiments have shown a strong discrepancy between the theoretical predictions and the measurements. Both the dynamical equations and the Weak Turbulence Theory treatment require some restrictive hypotheses. Here a direct numerical simulation of the Föppl-Von Kármán equations is performed and reproduces qualitatively and quantitatively the experimental results when the experimentally measured damping rate of waves $γ_\mathbf{k}= a + bk^2$ is used. This confirms that the Föppl-Von Kármán equations are a valid theoretical framework to describe our experiments. When we progressively tune the dissipation so that to localize it at the smallest scales, we observe a gradual transition between the experimental spectrum and the Kolmogorov-Zakharov prediction. Thus it is shown dissipation has a major influence on the scaling properties stationary solutions of weakly non linear wave turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Structures and dynamics of small scales in decaying magnetohydrodynamic turbulence

The topological and dynamical features of small scales are studied in the context of decaying magnetohydrodynamic turbulent flows using direct numerical simulations. Joint probability density functions (PDFs) of the invariants of gradient quantities related to the velocity and the magnetic fields demonstrate that structures and dynamics at the time of maximum dissipation depend on the large scale initial conditions. This is evident in particular from the fact that each flow has a different shape for the joint PDF of the invariants of the velocity gradient in contrast to the universal teardrop shape of hydrodynamic turbulence. The general picture that emerges from the analysis of the invariants is that regions of high vorticity are correlated with regions of high strain rate $\bm S$ also in contrast to hydrodynamic turbulent flows. Magnetic strain dominated regions are also well correlated with region of high current density $\bm j$. Viscous dissipation ($\propto \bm S^2$) as well as Ohmic dissipation ($\propto \bm j^2$) reside in regions where strain and rotation are locally almost in balance. The structures related to the velocity gradient possess different characteristics than those associated with the magnetic field gradient with the latter being locally more quasi-two dimensional.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Symmetry breaking and universality of decaying MHD Taylor-Green flows

We investigate the evolution and stability of a decaying magnetohydrodynamic (MHD) Taylor-Green flow. The chosen flow has been shown to result in a steep total energy spectrum with power law behaviour $k^{-2}$. We investigate the symmetry breaking of this flow by exciting perturbations of different amplitudes. It is shown that for any finite amplitude perturbation there is a high enough Reynolds number for which the perturbation will grow enough at the peak of dissipation rate resulting to a non-linear feedback in the flow and subsequently break the Taylor-Green symmetries. In particular, we show that symmetry breaking at large scales occurs if the amplitude of the perturbation is $ρ_{crit} \sim Re^{-1}$ and at small scales occurs if $ρ_{crit} \sim Re^{-3/2}$. This symmetry breaking modifies the scaling laws of the energy spectra at the peak of dissipation rate away from the $k^{-2}$ scaling and towards the classical $k^{-5/3}$ and $k^{-3/2}$ power laws.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The origins of $k^{-2}$ spectrum in the decaying Taylor-Green magnetohydrodynamic turbulent flows

We investigate the origins of $k^{-2}$ spectrum in a decaying Taylor-Green magnetohydrodynamic flow with zero large scale magnetic flux that was reported in Lee et al. (2010). A possible candidate for this scaling exponent has been the weak turbulence phenomenology. From our numerical simulations, we observe that current sheets in the magnetic Taylor-Green flow are formed in regions of magnetic discontinuities. Based on this observation and by studying the influence of the current sheets on the energy spectrum, using a filtering technique, we argue that the discontinuities are responsible for the $-2$ power law scaling of the energy spectra of this flow.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The signature of initial conditions on magnetohydrodynamic turbulence

We demonstrate that special correlations in the initial conditions of freely evolving, homogeneous magnetohydrodynamic (MHD) turbulence can lead to the formation of enormous current sheets. These coherent structures are observed at the peak of the energy dissipation rate and are the carriers of long-range correlations despite all the non-linear interactions during the formation of turbulence. Even though the largest scale separation has been achieved at this point, these structures are coherent with a size that spans our computational domain dominating the scaling of the energy spectrum, which follows a $E \propto k^{-2}$ power law. As Reynolds number increases curling of the current sheets, due to Kelvin-Helmhotlz type instabilities and reconnection, modifies the scaling of the energy spectrum from $k^{-2}$ towards $k^{-5/3}$. This transition occurs at the highest Reynolds numbers of direct numerical simulations with resolutions up to $2048^3$ grid points. Finite Reynolds number behaviour is observed due to the initial correlations without reaching a finite asymptote for the energy dissipation rate and with an unexpected scaling between the Taylor and the integral scale Reynolds numbers, i.e. $Re_λ\propto Re^{2/3}$. Our results, therefore, demonstrate that even state-of-the-art numerical simulations of the highest resolution can be influenced by the choice of initial conditions and consequently they are inadequate to deduce unequivocally the fate of universality in MHD turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transition from wave turbulence to dynamical crumpling in vibrated elastic plates

We study the dynamical regime of wave turbulence of a vibrated thin elastic plate based on experimental and numerical observations. We focus our study to the strongly non linear regime described in a previous letter by N. Yokoyama & M. Takaoka. At small forcing, a weakly non linear regime is compatible with the Weak Turbulence Theory when the dissipation is localized at high wavenumber. When the forcing intensity is increased, a strongly non linear regime emerges: singular structures dominate the dynamics at large scale whereas at small scales the weak turbulence is still present. A turbulence of singular structures, with folds and D-cones, develops that alters significantly the energy spectra and causes the emergence of intermittency.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Critical exponents in zero dimensions

In the vicinity of the onset of an instability, we investigate the effect of colored multiplicative noise on the scaling of the moments of the unstable mode amplitude. We introduce a family of zero dimensional models for which we can calculate the exact value of the critical exponents $β_m$ for all the moments. The results are obtained through asymptotic expansions that use the distance to onset as a small parameter. The examined family displays a variety of behaviors of the critical exponents that includes anomalous exponents: exponents that differ from the deterministic (mean-field) prediction, and multiscaling: non-linear dependence of the exponents on the order of the moment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From Weakly to Strongly Magnetized Isotropic MHD Turbulence

High Reynolds number isotropic magneto-hydro-dynamic turbulence in the presence of large scale magnetic fields is investigated as a function of the magnetic field strength. For a variety of flow configurations the energy dissipation rate ε, follows the Kolmogorov scaling ε~ U^3/L even when the large scale magnetic field energy is twenty times larger than the kinetic. Further increase of the magnetic energy showed a transition to the ε~ U^2 B / L scaling implying that magnetic shear becomes more efficient at this point at cascading the energy than the velocity fluctuations. Strongly helical configurations form helicity condensates that deviate from these scalings. Weak turbulence scaling was absent from the investigation. Finally, the magnetic energy spectra showed support for the Kolmogorov spectrum k^{-5/3} while kinetic energy spectra are closer to the Iroshnikov-Kraichnan spectrum k^{-3/2}.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bounding the scalar dissipation scale for mixing flows in the presence of sources

We investigate the mixing properties of scalars stirred by spatially smooth, divergence-free flows and maintained by a steady source-sink distribution. We focus on the spatial variation of the scalar field, described by the {\it dissipation wavenumber}, $k_d$, that we define as a function of the mean variance of the scalar and its gradient. We derive a set of upper bounds that for large Péclet number ($\Pe\gg1$) yield four distinct regimes for the scaling behavior of $k_d$, one of which corresponds to the Batchelor regime. The transition between these regimes is controlled by the value of $\Pe$ and the ratio $ρ=\ell_u/\ell_s$, where $\ell_u$ and $\ell_s$ are respectively, the characteristic lengthscales of the velocity and source fields. A fifth regime is revealed by homogenization theory. These regimes reflect the balance between different processes: scalar injection, molecular diffusion, stirring and bulk transport from the sources to the sinks. We verify the relevance of these bounds by numerical simulations for a {two-dimensional, chaotically mixing} example flow and discuss their relation to previous bounds. Finally, we note some implications for three dimensional turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the 2D behavior of 3D MHD with a strong guiding field

The Magneto-hydrodynamic (MHD) equations in the presence of a guiding magnetic field are investigated by means of direct numerical simulations. The basis of the investigation consists of 9 runs forced at the small scales. The results demonstrate that for a large enough uniform magnetic field the large scale flow behaves as a two dimensional (non-MHD) fluid exhibiting an inverse cascade of energy in the direction perpendicular to the magnetic field, while the small scales behave like a three dimensional MHD-fluid cascading the energy forwards. The amplitude of the inverse cascade is sensitive to the magnetic field amplitude, the domain size, the forcing mechanism, and the forcing scale. All these dependencies are demonstrated by the varying parameters of simulations. Furthermore, in the case that the system is forced anisotropically in the small parallel scales an inverse cascade in the parallel direction is observed that is feeding the 2D modes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Searching for the fastest dynamo: Laminar ABC flows

The growth rate of the dynamo instability as a function of the magnetic Reynolds number Rm is investigated by means of numerical simulations for the family of the ABC flows and for 2 different forcing scales. For the ABC flows that are driven at the largest available length scale it is found that as the magnetic Reynolds number is increased: (a) The flow that results first in dynamo is the 2.5D flow for which A=B and C=0 (and all permutations). (b) The second type of flow that results in dynamo is the one for which A=B=2C/5 (and permutations). (c) The most symmetric flow A=B=C is the third type of flow that results in dynamo. (d) As Rm is increased, the A=B=C flow stops being a dynamo and transitions from a local maximum to a third-order saddle point. (e) At larger Rm the A=B=C flow re-establishes its self as a dynamo but remains a saddle point. (f) At the largest examined Rm the growth rate of the 2.5D flows starts to decrease, the A=B=C flow comes close to a local maximum again and the flow A=B=2C/5 (and permutations) results in the fastest dynamo with growth rate $γ~0.12$ at the largest examined Rm. For the ABC flows that are driven at the second largest available length scale it is found that (a) the 2.5D flows A=B, C=0 (and permutations) are again the first flows that result in dynamo with a decreased onset. (b) The most symmetric flow A=B=C is the second type of flow that results in dynamo. It is and remains a local maximum. (c) At larger Rm the flow A=B=2C/5 (and permutations) appears as the third type of flow that results in dynamo. As Rm is increased it becomes the flow with the largest growth rate. The growth rates appear to have some correlation with the Lyaponov exponents but constructive re-folding of the field lines appears equally important in determining the fasted dynamo flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nonlinear dynamos at infinite magnetic Prandtl number

The dynamo instability is investigated in the limit of infinite magnetic Prandtl number. In this limit the fluid is assumed to be very viscous so that the inertial terms can be neglected and the flow is slaved to the forcing. The forcing consist of an external forcing function that drives the dynamo flow and the resulting Lorentz force caused by the back reaction of the magnetic field. The flows under investigation are the Archontis flow, and the ABC flow forced at two different scales. The investigation covers roughly three orders of magnitude of the magnetic Reynolds number above onset. All flows show a weak increase of the averaged magnetic energy as the magnetic Reynolds number is increased. Most of the magnetic energy is concentrated in flat elongated structures that produce a Lorentz force with small solenoidal projection so that the resulting magnetic field configuration was almost force-free. Although the examined system has zero kinetic Reynolds number at sufficiently large magnetic Reynolds number the structures are unstable to small scale fluctuations that result in a chaotic temporal behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stratified shear flow instabilities at large Richardson numbers

Numerical simulations of stratified shear flow instabilities are performed in two dimensions in the Boussinesq limit. The density variation length scale is chosen to be four times smaller than the velocity variation length scale so that Holmboe or Kelvin-Helmholtz unstable modes are present depending on the choice of the global Richardson number Ri. Three different values of Ri were examined Ri =0.2, 2, 20. The flows for the three examined values are all unstable due to different modes namely: the Kelvin-Helmholtz mode for Ri=0.2, the first Holmboe mode for Ri=2, and the second Holmboe mode for Ri=20 that has been discovered recently and it is the first time that it is examined in the non-linear stage. It is found that the amplitude of the velocity perturbation of the second Holmboe mode at the non-linear stage is smaller but comparable to first Holmboe mode. The increase of the potential energy however due to the second Holmboe modes is greater than that of the first mode. The Kelvin-Helmholtz mode is larger by two orders of magnitude in kinetic energy than the Holmboe modes and about ten times larger in potential energy than the Holmboe modes. The results in this paper suggest that although mixing is suppressed at large Richardson numbers it is not negligible, and turbulent mixing processes in strongly stratified environments can not be excluded.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Shell to shell energy transfer in MHD, Part I: steady state turbulence

We investigate the transfer of energy from large scales to small scales in fully developed forced three-dimensional MHD-turbulence by analyzing the results of direct numerical simulations in the absence of an externally imposed uniform magnetic field. Our results show that the transfer of kinetic energy from the large scales to kinetic energy at smaller scales, and the transfer of magnetic energy from the large scales to magnetic energy at smaller scales, are local, as is also found in the case of neutral fluids, and in a way that is compatible with Kolmogorov (1941) theory of turbulence. However, the transfer of energy from the velocity field to the magnetic field is a highly non-local process in Fourier space. Energy from the velocity field at large scales can be transfered directly into small scale magnetic fields without the participation of intermediate scales. Some implications of our results to MHD turbulence modeling are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Shell to shell energy transfer in MHD, Part II: Kinematic dynamo

We study the transfer of energy between different scales for forced three-dimensional MHD turbulent flows in the kinematic dynamo regime. Two different forces are examined: a non-helical Taylor Green flow with magnetic Prandtl number P_M=0.4, and a helical ABC flow with P_M=1. This analysis allows us to examine which scales of the velocity flow are responsible for dynamo action, and identify which scales of the magnetic field receive energy directly from the velocity field and which scales receive magnetic energy through the cascade of the magnetic field from large to small scales. Our results show that the turbulent velocity fluctuations are responsible for the magnetic field amplification in the small scales (small scale dynamo) while the large scale field is amplified mostly due to the large scale flow. A direct cascade of the magnetic field energy from large to small scales is also present and is a complementary mechanism for the increase of the magnetic field in the small scales. Input of energy from the velocity field in the small magnetic scales dominates over the energy that is cascaded down from the large scales until the large-scale peak of the magnetic energy spectrum is reached. At even smaller scales, most of the magnetic energy input is from the cascading process.

0 citations0 reviewsNo code linkedOpen access