Search

Intermittency as it occurs in fast dynamos in the MHD framework is evaluated through the examination of relations between normalized moments at third order (skewness S) and fourth order (kurtosis K) for both the velocity and magnetic field, and for their local dissipations. As investigated by several authors in various physical contexts such as fusion plasmas (Krommes (2008)), climate evolution (Sura & Sardeshmukh (2008)), fluid turbulence or rotating stratified flows (Pouquet et al. (2023)), approximate parabolic K(S) ~ alpha laws emerge whose origin may be related to the applicability of intermittency models to their dynamics. The results analyzed herein are obtained through direct numerical simulations of MHD flows for both Taylor-Green and Beltrami ABC forcing at moderate Reynolds numbers, and for up to 314000 turn-over times. We observe for the dissipation 0.2 < alpha <3.0, an evaluation that varies with the field, the forcing, and when filtering for high-skewness intermittent structures. When using the She & Leveque (1994) intermittency model, one can compute alpha analytically; we then find alpha = 2.5, clearly differing from a (strict) parabolic scaling, a result consistent with the numerical data.

Nonlinear phenomena and turbulence are central to our understanding and modeling the dynamics of fluids and plasmas, and yet they still resist analytical resolutions in many instances. However, progress has been made recently, displaying a richness of phenomena which was somewhat unexpected a few years back, such as the double constant-flux cascades of a same invariant to both the large and to the small scales, or the presence of non-Gaussian wings in the large-scale fields, for fluids and plasmas. Here, I will concentrate on the direct measurement of the magnitude of dissipation and an evaluation of intermittency in a turbulent plasma using exact laws stemming from invariance principles and involving cross-correlation tensors with both the velocity and the magnetic fields. I will illustrate these points through scaling laws, together with data analysis from existing experiments, observations and numerical simulations. Finally, I will also briefly explore the possible implications for validity and use of several modeling strategies. To appear, Plasma Physics and Controlled Fusion, 2023.

Space plasmas are three-dimensional dynamic entities. Except under very special circumstances, their structure in space and their behavior in time are not related in any simple way. Therefore, single spacecraft in situ measurements cannot unambiguously unravel the full space-time structure of the heliospheric plasmas of interest in the inner heliosphere, in the Geospace environment, or the outer heliosphere. This shortcoming leaves numerous central questions incompletely answered. Deficiencies remain in at least two important subjects, Space Weather and fundamental plasma turbulence theory, due to a lack of a more complete understanding of the space-time structure of dynamic plasmas. Only with multispacecraft measurements over suitable spans of spatial separation and temporal duration can these ambiguities be resolved. We note that these characterizations apply to turbulence across a wide range of scales, and also equally well to shocks, flux ropes, magnetic clouds, current sheets, stream interactions, etc. In the following, we will describe the basic requirements for resolving space-time structure in general, using turbulence' as both an example and a principal target or study. Several types of missions are suggested to resolve space-time structure throughout the Heliosphere., Comment: White Paper submitted to: Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033. arXiv admin note: substantial text overlap with arXiv:1903.06890

We observe the emergence of strong vertical drafts in direct numerical simulations of the Boussinesq equations in a range of parameters of geophysical interest. These structures, which appear intermittently in space and time, generate turbulence and enhance kinetic and potential energy dissipation, providing a possible explanation for the observed variability of the local energy dissipation in the bulk of oceanic flows, and the modulation of its probability distribution function. We show how, due to the extreme drafts, in runs with Froude numbers observable in geophysical scenarios, roughly 10% of the domain flow can account for up to 50% of the global volume dissipation, reminiscent of estimates based on oceanic models., Comment: 16 pages, 6 figures

Non-Gaussian statistics of large-scale fields are routinely observed in data from atmospheric and oceanic campaigns and global models. Recent direct numerical simulations (DNSs) showed that large-scale intermittency in stably stratified flows is due to the emergence of sporadic, extreme events in the form of bursts in the vertical velocity and the temperature. This phenomenon results from the interplay between waves and turbulent motions, affecting mixing. We provide evidence of the enhancement of the classical small-scale (or internal) intermittency due to the emergence of large-scale drafts, connecting large- and small-scale bursts. To this aim we analyze a large set of DNSs of the stably stratified Boussinesq equations over a wide range of values of the Froude number ($Fr\approx 0.01-1$). The variation of the buoyancy field kurtosis with $Fr$ is similar to (though with smaller values than) the kurtosis of the vertical velocity, both showing a non-monotonic trend. We present a mechanism for the generation of extreme vertical drafts and vorticity enhancements which follows from the exact equations for field gradients., Comment: 7 pages, 4 figures

Helicity, a measure of the breakage of reflectional symmetry representing the topology of turbulent flows, contributes in a crucial way to their dynamics and to their fundamental statistical properties. We review several of their main features, both new and old, as the discovery of bi-directional cascades or the role of helical vortices in the enhancement of large-scale magnetic fields in the dynamo problem. The dynamical contribution in magnetohydrodynamic (MHD) of the cross-correlation between velocity and induction is discussed as well. We consider next how turbulent transport is affected by helical constraints, in particular in the context of magnetic reconnection and fusion plasmas {under one- and two-fluid approximations}. Central issues on how to construct turbulence models for non-reflectionally symmetric helical flows are reviewed, including in the presence of shear, and we finally briefly mention the possible role of helicity in the development of strongly localized quasi-singular structures {at small scale)., Comment: 14 pages, 0 figure 0 table

In turbulence, for neutral or conducting fluids, a large ratio of scales is excited because of the possible occurrence of inverse cascades to large, global scales together with direct cascades to small, dissipative scales, as observed in the atmosphere and oceans, or in the solar environment. In this context, using direct numerical simulations with forcing, we analyze scale dynamics in the presence of magnetic fields with a generalized Ohm's law including a Hall current. The ion inertial length epsilon_H serves as the control parameter at fixed Reynolds number. Both the magnetic and generalized helicity -- invariants in the ideal case -- grow linearly with time, as expected from classical arguments. The cross-correlation between the velocity and magnetic field grows as well, more so in relative terms for a stronger Hall current. We find that the helical growth rates vary exponentially with epsilon_H, provided the ion inertial scale resides within the inverse cascade range. These exponential variations are recovered phenomenologically using simple scaling arguments. They are directly linked to the wavenumber power-law dependence of generalized and magnetic helicity, k^(-2), in their inverse ranges. This illustrates and confirms the important role of the interplay between large and small scales in the dynamics of turbulent flows.

We study the spatio-temporal behavior of the Els\"asser variables describing magnetic and velocity field fluctuations, using direct numerical simulations of three-dimensional magnetohydrodynamic turbulence. We consider cases with relatively small, intermediate, and large values of a mean background magnetic field, and with null, small, and high cross-helicity (correlations between the velocity and the magnetic field). Wavenumber-dependent time correlation functions are computed for the different simulations. From these correlation functions, the decorrelation time is computed and compared with different theoretical characteristic times: the local non-linear time, the random-sweeping time, and the Alfv\'enic time. It is found that decorrelation times are dominated by sweeping effects for low values of the mean magnetic field and for low values of the cross-helicity, while for large values of the background field or of the cross-helicity and for wave vectors sufficiently aligned with the guide field, decorrelation times are controlled by Alfv\'enic effects. Finally, we observe counter-propagation of Alfv\'enic fluctuations due to reflections produced by inhomogeneities in the total magnetic field. This effect becomes more prominent in flows with large cross-helicity, strongly modifying the propagation of waves in turbulent magnetohydrodynamic flows., Comment: 16 pages, 40 figures

Geophysical fluid flows are predominantly turbulent and often strongly affected by the Earth's rotation, as well as by stable density stratification. Using direct numerical simulations of forced Boussinesq equations, we study the influence of these effects on the motion of fluid particles, focusing on cases where the frequencies associated with rotation and stratification (RaS), $N$ and $f$ respectively, are held at a fixed ratio $N/f=5$. As the intensity of RaS increases, a sharp transition is observed between a regime dominated by eddies to a regime dominated by waves, which can also be seemingly described by simply comparing the time scale $1/N$ and $\tau_\eta$ (the Kolmogorov time scale). We perform a detailed study of Lagrangian statistics of acceleration, velocity and related quantities in the two regimes. The flow anisotropy induces a clear difference between particle motion in the horizontal and vertical directions. In the regime $N\tau_\eta <1$, acceleration statistics in both horizontal and vertical directions, exhibit well known characteristics of isotropic turbulence. In contrast for $N\tau_\eta >1$, they are directly influenced by imposed RaS. The Lagrangian velocity statistics exhibit visible anisotropy for all runs; nevertheless the degree of anisotropy becomes very strong in the regime $N\tau_\eta >1$. We find that in the regime $N\tau_\eta <1$, rotation enhances the mean displacement of particles in horizontal planes at short times, but inhibits them at longer times. This inhibition of horizontal displacement becomes stronger for $N\tau_\eta >1$, with no clear diffusive behavior. Displacements in the vertical direction are always inhibited. The inhibition becomes extremely strong when $N\tau_\eta >1$, with the particles almost being trapped horizontally., Comment: 24 pages, 10 figures

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence., Comment: White paper submitted to the PLASMA 2020 Decadal Survey Committee

An existing hybrid MPI-OpenMP scheme is augmented with a CUDA-based fine grain parallelization approach for multidimensional distributed Fourier transforms, in a well-characterized pseudospectral fluid turbulence code. Basics of the hybrid scheme are reviewed, and heuristics provided to show a potential benefit of the CUDA implementation. The method draws heavily on the CUDA runtime library to handle memory management, and on the cuFFT library for computing local FFTs. The manner in which the interfaces are constructed to these libraries, and ISO bindings utilized to facilitate platform portability, are discussed. CUDA streams are implemented to overlap data transfer with cuFFT computation. Testing with a baseline solver demonstrates significant aggregate speed-up over the hybrid MPI-OpenMP solver by offloading to GPUs on an NVLink-based test system. While the batch streamed approach provides little benefit with NVLink, we see a performance gain of 30% when tuned for the optimal number of streams on a PCIe-based system. It is found that strong GPU scaling is ideal, or slightly better than ideal, in all cases. In addition to speed-up measurements for the fiducial solver, we also consider several other solvers with different numbers of transform operations and find that aggregate speed-ups are nearly constant for all solvers., Comment: 17 pages, 8 figures; submitted to Parallel Computing

In this paper, we present an overview of concepts and data concerning inverse cascades of excitation towards scales larger than the forcing scale in a variety of contexts, from two-dimensional fluids and wave turbulence, to geophysical flows in the presence of rotation and stratification. We briefly discuss the role of anisotropy in the occurrence and properties of such cascades. We then show that the cascade of some invariant, for example the total energy, may be transferred through nonlinear interactions to both the small scales and the large scales, with in each case a constant flux. This is in contrast to the classical picture, and we illustrate such a dual cascade in the context of atmospheric and oceanic observations, direct numerical simulations and modeling. We also show that this dual cascade of total energy can in fact be decomposed in some cases into separate cascades of the kinetic and potential energies, provided the Froude and Rossby numbers are small enough. In all cases, the potential energy flux remains small, of the order of 10% or less relative to the kinetic energy flux. Finally, we demonstrate that, in the small-scale inertial range, approximate equipartition between potential and kinetic modes is obtained, leading to an energy ratio close to one, with strong departure at large scales due to the dominant kinetic energy inverse cascade and piling-up at the lowest spatial frequency, and at small scales due to unbalanced dissipation processes, even though the Prandtl number is equal to one.

The effects of changing the orientation of the rotation axis on homogeneous turbulence is considered. We perform direct numerical simulations on a periodic box of $1024^3$ grid points, where the orientation of the rotation axis is changed (a) at a fixed time instant (b) regularly at time intervals commensurate with the rotation time scale. The former is characterized by a dominant inverse energy cascade whereas in the latter, the inverse cascade is stymied due to the recurrent changes in the rotation axis resulting in a strong forward energy transfer and large scale structures that resemble those of isotropic turbulence., Comment: 7 pages, 8 figures, The European Physical Journal E (EPJ E)

The interplay between waves and eddies in stably stratified rotating flows is investigated by means of world-class direct numerical simulations using up to $3072^3$ grid points. Strikingly, we find that the shift from vortex to wave dominated dynamics occurs at a wavenumber $k_R$ which does not depend on Reynolds number, suggesting that partition of energy between wave and vortical modes is not sensitive to the development of turbulence at the smaller scales. We also show that $k_R$ is comparable to the wavenumber at which exchanges between kinetic and potential modes stabilize at close to equipartition, emphasizing the role of potential energy, as conjectured in the atmosphere and the oceans. Moreover, $k_R$ varies as the inverse of the Froude number as explained by the scaling prediction proposed, consistent with recent observations and modeling of the Mesosphere-Lower Thermosphere and of the ocean.

We study the partition of energy between waves and vortices in stratified turbulence, with or without rotation, for a variety of parameters, focusing on the behavior of the waves and vortices in the inverse cascade of energy towards the large scales. To this end, we use direct numerical simulations in a cubic box at a Reynolds number Re=1000, with the ratio between the Brunt-V\"ais\"al\"a frequency N and the inertial frequency f varying from 1/4 to 20, together with a purely stratified run. The Froude number, measuring the strength of the stratification, varies within the range 0.02 < Fr < 0.32. We find that the inverse cascade is dominated by the slow quasi-geostrophic modes. Their energy spectra and fluxes exhibit characteristics of an inverse cascade, even though their energy is not conserved. Surprisingly, the slow vortices still dominate when the ratio N/f increases, also in the stratified case, although less and less so. However, when N/f increases, the inverse cascade of the slow modes becomes weaker and weaker, and it vanishes in the purely stratified case. We discuss how the disappearance of the inverse cascade of energy with increasing N/f can be interpreted in terms of the waves and vortices, and identify three major effects that can explain this transition based on inviscid invariants arguments.

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions., Comment: 51 pages, 6 figures (Figs 2, 3, and 4 color), accepted to Space Science Reviews (in press). The paper is a product of a workshop on "LES Modeling in MHD Turbulence" held in Boulder, CO in May, 2013, sponsored by the Geophysical Turbulence Program at the National Center for Atmospheric Research

Most of the turbulent flows appearing in nature (e.g. geophysical and astrophysical flows) are subjected to strong rotation and stratification. These effects break the symmetries of classical, homogenous isotropic turbulence. In doing so, they introduce a natural decomposition of phase space in terms of wave modes and potential vorticity modes. The appearance of a new time scale associated to the propagation of waves, in addition to the eddy turnover time, increases the complexity of the energy transfers between the various scales; nonlinearly interacting waves may dominate at some scales while balanced motion may prevail at others. In the end, it is difficult to predict \emph{a priori} if the energy cascades downscale as in homogeneous isotropic turbulence, upscale as expected from balanced dynamics, or follows yet another phenomenology. In this paper, we suggest a theoretical approach based on equilibrium statistical mechanics for the ideal system, inspired from the restricted partition function formalism introduced in metastability studies. In this framework, we show analytically that in the presence of rotation, when the dynamics is restricted to the slow modes, the equilibrium energy spectrum features an infrared divergence characteristic of an inverse cascade regime, whereas this is not the case for purely stratified flows.

We perform numerical simulations of decaying rotating stratified turbulence and show, in the Boussinesq framework, that helicity (velocity-vorticity correlation), as observed in super-cell storms and hurricanes, is spontaneously created due to an interplay between buoyancy and rotation common to large-scale atmospheric and oceanic flows. Helicity emerges from the joint action of eddies and of inertia-gravity waves (with inertia and gravity with respective associated frequencies $f$ and $N$), and it occurs when the waves are sufficiently strong. For $N/f < 3$ the amount of helicity produced is correctly predicted by a quasi-linear balance equation. Outside this regime, and up to the highest Reynolds number obtained in this study, namely $Re\approx 10000$, helicity production is found to be persistent for $N/f$ as large as $\approx 17$, and for $ReFr^2$ and $ReRo^2 $ respectively as large as $\approx 100$ and $\approx 24000$., Comment: 10 pages, 5 figures

A review of some of the issues that have arisen over the years concerning the energy distribution among scales for magnetohydrodynamics (MHD) turbulence is given here. A variety of tools are employed to that effect, and a central role is played by taking into consideration the ideal (non-dissipative) invariants, namely the total energy, the magnetic helicity and the cross-correlations between the velocity and the magnetic field (concentrating on the three-dimensional case). These concepts, based mostly on theory, models and direct numerical simulations, are briefly put in the context of observations, in particular the solar wind, and some of the remaining open questions are delineated as well. New results on ideal MHD dynamics in three dimensions on equivalent grids of up to $6144^3$ points using the Taylor-Green flow generalized to MHD are also mentioned., Comment: 17 pages, 4 figures, review

We investigate long-time properties of three-dimensional MHD turbulence in the absence of forcing and examine in particular the role played by the quadratic invariants of the system and by the symmetries of the initial configurations. We observe that, when sufficient accuracy is used, initial conditions with a high degree of symmetries, as in the absence of helicity, do not travel through parameter space over time whereas by perturbing these solutions either explicitly or implicitly using for example single precision for long times, the flows depart from their original behavior and can become either strongly helical, or have a strong alignment between the velocity and the magnetic field. When the symmetries are broken, the flows evolve towards different end states, as predicted by statistical arguments for non-dissipative systems with the addition of an energy minimization principle, as already analyzed in \cite{stribling_90} for random initial conditions using a moderate number of Fourier modes. Furthermore, the alignment properties of these flows, between velocity, vorticity, magnetic potential, induction and current, correspond to the dominance of two main regimes, one helically dominated and one in quasi-equipartition of kinetic and magnetic energy. We also contrast the scaling of the ratio of magnetic energy to kinetic energy as a function of wavenumber to the ratio of eddy turn-over time to Alfv\'en time as a function of wavenumber. We find that the former ratio is constant with an approximate equipartition for scales smaller than the largest scale of the flow whereas the ratio of time scales increases with increasing wavenumber., Comment: 14 pages, 6 figures

Rapidly rotating turbulent flow is characterized by the emergence of columnar structures that are representative of quasi-two dimensional behavior of the flow. It is known that when energy is injected into the fluid at an intermediate scale $L_f$, it cascades towards smaller as well as larger scales. In this paper we analyze the flow in the \textit{inverse cascade} range at a small but fixed Rossby number, {$\mathcal{R}o_f \approx 0.05$}. Several {numerical simulations with} helical and non-helical forcing functions are considered in periodic boxes with unit aspect ratio. In order to resolve the inverse cascade range with {reasonably} large Reynolds number, the analysis is based on large eddy simulations which include the effect of helicity on eddy viscosity and eddy noise. Thus, we model the small scales and resolve explicitly the large scales. We show that the large-scale energy spectrum has at least two solutions: one that is consistent with Kolmogorov-Kraichnan-Batchelor-Leith phenomenology for the inverse cascade of energy in two-dimensional (2D) turbulence with a {$\sim k_{\perp}^{-5/3}$} scaling, and the other that corresponds to a steeper {$\sim k_{\perp}^{-3}$} spectrum in which the three-dimensional (3D) modes release a substantial fraction of their energy per unit time to 2D modes. {The spectrum that} emerges {depends on} the anisotropy of the forcing function{,} the former solution prevailing for forcings in which more energy is injected into 2D modes while the latter prevails for isotropic forcing. {In the case of anisotropic forcing, whence the energy} goes from the 2D to the 3D modes at low wavenumbers, large-scale shear is created resulting in another time scale $\tau_{sh}$, associated with shear, {thereby producing} a $\sim k^{-1}$ spectrum for the {total energy} with the 2D modes still following a {$\sim k_{\perp}^{-5/3}$} scaling.

Understanding the in situ amplification of large scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f_h,C has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that f_h,C decreases as the ratio of forcing to large scale wave numbers k_F/k_min increases. From the condition that a large scale helical dynamo must overcome the backreaction from any non-helical field on the large scales, we develop a theory that can explain the simulations. For k_F/k_min>8 we find f_h,C< 3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate scale separation., Comment: 5 pages, 4 figures

We investigate numerically the dynamics of two-dimensional Euler and ideal magnetohydrodynamics (MHD) flows in systems with a finite number of modes, up to $4096^2$, for which several quadratic invariants are preserved by the truncation and the statistical equilibria are known. Initial conditions are the Orszag-Tang vortex with a neutral X-point centered on a stagnation point of the velocity field in the large scales. In MHD, we observe that the total energy spectra at intermediate times and intermediate scales correspond to the interactions of eddies and waves, $E_T(k)\sim k^{-3/2}$. Moreover, no dissipative range is visible neither for Euler nor for MHD in two dimensions; in the former case, this may be linked to the existence of a vanishing turbulent viscosity whereas in MHD, the numerical resolution employed may be insufficient. When imposing a uniform magnetic field to the flow, we observe a lack of saturation of the formation of small scales together with a significant slowing-down of their equilibration, with however a cut-off independent partial thermalization being reached at intermediate scales.

We present a systematic analysis of statistical properties of turbulent current and vorticity structures at a given time using cluster analysis. The data stems from numerical simulations of decaying three-dimensional (3D) magnetohydrodynamic turbulence in the absence of an imposed uniform magnetic field; the magnetic Prandtl number is taken equal to unity, and we use a periodic box with grids of up to 1536^3 points, and with Taylor Reynolds numbers up to 1100. The initial conditions are either an X-point configuration embedded in 3D, the so-called Orszag-Tang vortex, or an Arn'old-Beltrami-Childress configuration with a fully helical velocity and magnetic field. In each case two snapshots are analyzed, separated by one turn-over time, starting just after the peak of dissipation. We show that the algorithm is able to select a large number of structures (in excess of 8,000) for each snapshot and that the statistical properties of these clusters are remarkably similar for the two snapshots as well as for the two flows under study in terms of scaling laws for the cluster characteristics, with the structures in the vorticity and in the current behaving in the same way. We also study the effect of Reynolds number on cluster statistics, and we finally analyze the properties of these clusters in terms of their velocity-magnetic field correlation. Self-organized criticality features have been identified in the dissipative range of scales. A different scaling arises in the inertial range, which cannot be identified for the moment with a known self-organized criticality class consistent with MHD. We suggest that this range can be governed by turbulence dynamics as opposed to criticality, and propose an interpretation of intermittency in terms of propagation of local instabilities., Comment: 17 pages, 9 figures, 5 tables

A hybrid scheme that utilizes MPI for distributed memory parallelism and OpenMP for shared memory parallelism is presented. The work is motivated by the desire to achieve exceptionally high Reynolds numbers in pseudospectral computations of fluid turbulence on emerging petascale, high core-count, massively parallel processing systems. The hybrid implementation derives from and augments a well-tested scalable MPI-parallelized pseudospectral code. The hybrid paradigm leads to a new picture for the domain decomposition of the pseudospectral grids, which is helpful in understanding, among other things, the 3D transpose of the global data that is necessary for the parallel fast Fourier transforms that are the central component of the numerical discretizations. Details of the hybrid implementation are provided, and performance tests illustrate the utility of the method. It is shown that the hybrid scheme achieves near ideal scalability up to ~20000 compute cores with a maximum mean efficiency of 83%. Data are presented that demonstrate how to choose the optimal number of MPI processes and OpenMP threads in order to optimize code performance on two different platforms., Comment: Submitted to Parallel Computing

The subfilter-scale (SFS) physics of regularization models are investigated to understand the regularizations' performance as SFS models. The strong suppression of spectrally local SFS interactions and the conservation of small-scale circulation in the Lagrangian-averaged Navier-Stokes alpha-model (LANS-alpha) is found to lead to the formation of rigid bodies. These contaminate the superfilter-scale energy spectrum with a scaling that approaches k^1 as the SFS spectra is resolved. The Clark-alpha and Leray-alpha models, truncations of LANS-alpha, do not conserve small-scale circulation and do not develop rigid bodies. LANS-alpha, however, is closest to Navier-Stokes in intermittency properties. All three models are found to be stable at high Reynolds number. Differences between L^2 and H^1 norm models are clarified. For magnetohydrodynamics (MHD), the presence of the Lorentz force as a source (or sink) for circulation and as a facilitator of both spectrally nonlocal large to small scale interactions as well as local SFS interactions prevents the formation of rigid bodies in Lagrangian-averaged MHD (LAMHD-alpha). We find LAMHD-alpha performs well as a predictor of superfilter-scale energy spectra and of intermittent current sheets at high Reynolds numbers. We expect it may prove to be a generally applicable MHD-LES., Comment: Proceedings of the Quality and Reliability of Large-Eddy Simulations II 2009 Workshop, Springer ERCOFTAC volume 16 (2011); submitted to Journal of Scientific Computing; 22 pages, 8 figures

We show in this paper that third- and fourth-order low storage Runge-Kutta algorithms can be built specifically for quadratic nonlinear operators, at the expense of roughly doubling the time needed for evaluating the temporal derivatives. The resulting algorithms are especially well suited for computational fluid dynamics. Examples are given for the H\'enon-Heiles Hamiltonian system and, in one and two space dimensions, for the Burgers equation using both a pseudo-spectral code and a spectral element code, respectively. The scheme is also shown to be practical in three space solving the incompressible Euler equation using a fully parallelized pseudo-spectral code.

We demonstrate that, for the case of quasi-equipartition between the velocity and the magnetic field, the Lagrangian-averaged magnetohydrodynamics alpha-model (LAMHD) reproduces well both the large-scale and small-scale properties of turbulent flows; in particular, it displays no increased (super-filter) bottleneck effect with its ensuing enhanced energy spectrum at the onset of the sub-filter-scales. This is in contrast to the case of the neutral fluid in which the Lagrangian-averaged Navier-Stokes alpha-model is somewhat limited in its applications because of the formation of spatial regions with no internal degrees of freedom and subsequent contamination of super-filter-scale spectral properties. No such regions are found in LAMHD, making this method capable of large reductions in required numerical degrees of freedom; specifically, we find a reduction factor of 200 when compared to a direct numerical simulation on a large grid of 1536^3 points at the same Reynolds number., Comment: 22 pages, 9 figures; accepted Phys.Rev.E

New generalized equations of motion for the Weber-Clebsch potentials that describe both the Navier-Stokes and MHD dynamics are derived. These depend on a new parameter, which has dimensions of time for Navier-Stokes and inverse velocity for MHD. Direct numerical simulations are performed. For Navier-Stokes, the generalized formalism captures the intense reconnection of vortices of the Boratav, Pelz and Zabusky flow, in agreement with the previous study by Ohkitani and Constantin. For MHD, the new formalism is used to detect magnetic reconnection in several flows: the 3D Arnold, Beltrami and Childress (ABC) flow and the (2D and 3D) Orszag-Tang vortex. It is concluded that periods of intense activity in the magnetic enstrophy are correlated with periods of increasingly frequent resettings. Finally, the positive correlation between the sharpness of the increase in resetting frequency and the spatial localization of the reconnection region is discussed.

We determine how the differences in the treatment of the subfilter-scale physics affect the properties of the flow for three closely related regularizations of Navier-Stokes. The consequences on the applicability of the regularizations as SGS models are also shown by examining their effects on superfilter-scale properties. Numerical solutions of the Clark-alpha model are compared to two previously employed regularizations, LANS-alpha and Leray-alpha (at Re ~ 3300, Taylor Re ~ 790) and to a DNS. We derive the Karman-Howarth equation for both the Clark-alpha and Leray-alpha models. We confirm one of two possible scalings resulting from this equation for Clark as well as its associated k^(-1) energy spectrum. At sub-filter scales, Clark-alpha possesses similar total dissipation and characteristic time to reach a statistical turbulent steady-state as Navier-Stokes, but exhibits greater intermittency. As a SGS model, Clark reproduces the energy spectrum and intermittency properties of the DNS. For the Leray model, increasing the filter width decreases the nonlinearity and the effective Re is substantially decreased. Even for the smallest value of alpha studied, Leray-alpha was inadequate as a SGS model. The LANS energy spectrum k^1, consistent with its so-called "rigid bodies," precludes a reproduction of the large-scale energy spectrum of the DNS at high Re while achieving a large reduction in resolution. However, that this same feature reduces its intermittency compared to Clark-alpha (which shares a similar Karman-Howarth equation). Clark is found to be the best approximation for reproducing the total dissipation rate and the energy spectrum at scales larger than alpha, whereas high-order intermittency properties for larger values of alpha are best reproduced by LANS-alpha., Comment: 21 pages, 8 figures

We present a new version of a dynamical spectral model for Large Eddy Simulation based on the Eddy Damped Quasi Normal Markovian approximation \cite{sao,chollet_lesieur}. Three distinct modifications are implemented and tested. On the one hand, whereas in current approaches, a Kolmogorov-like energy spectrum is usually assumed in order to evaluate the nonlocal transfer, in our method the energy spectrum of the subgrid scales adapts itself dynamically to the large-scale resolved spectrum; this first modification allows in particular for a better treatment of transient phases and instabilities, as shown on one specific example. Moreover, the model takes into account the phase relationships of the small-scales, embodied for example in strong localized structures such as vortex filaments. To that effect, phase information is implemented in the treatment of the so-called eddy noise in the closure model. Finally, we also consider the role that helical small scales may play in the evaluation of the transfer of energy and helicity, the two invariants of the primitive equations in the inviscid case; this leads as well to intrinsic variations in the development of helicity spectra. Therefore, our model allows for simulations of flows for a variety of circumstances and a priori at any given Reynolds number. Comparisons with Direct Numerical Simulations of the three-dimensional Navier-Stokes equation are performed on fluids driven by an ABC (Beltrami) flow which is a prototype of fully helical flows. Good agreements are obtained for physical and spectral behavior of the large scales.

We compute solutions of the Lagrangian-Averaged Navier-Stokes alpha-model (LANS) for significantly higher Reynolds numbers (up to Re 8300) than have previously been accomplished. This allows sufficient separation of scales to observe a Navier-Stokes (NS) inertial range followed by a 2nd LANS inertial range. The analysis of the third-order structure function scaling supports the predicted l^3 scaling; it corresponds to a k^(-1) scaling of the energy spectrum. The energy spectrum itself shows a different scaling which goes as k^1. This latter spectrum is consistent with the absence of stretching in the sub-filter scales due to the Taylor frozen-in hypothesis employed as a closure in the derivation of LANS. These two scalings are conjectured to coexist in different spatial portions of the flow. The l^3 (E(k) k^(-1)) scaling is subdominant to k^1 in the energy spectrum, but the l^3 scaling is responsible for the direct energy cascade, as no cascade can result from motions with no internal degrees of freedom. We verify the prediction for the size of the LANS attractor resulting from this scaling. From this, we give a methodology either for arriving at grid-independent solutions for LANS, or for obtaining a formulation of a LES optimal in the context of the alpha models. The fully converged grid-independent LANS may not be the best approximation to a direct numerical simulation of the NS equations since the minimum error is a balance between truncation errors and the approximation error due to using LANS instead of the primitive equations. Furthermore, the small-scale behavior of LANS contributes to a reduction of flux at constant energy, leading to a shallower energy spectrum for large alpha. These small-scale features, do not preclude LANS to reproduce correctly the intermittency properties of high Re flow., Comment: 37 pages, 17 figures

We compute numerically the threshold for dynamo action in Taylor-Green swirling flows. Kinematic calculations, for which the flow field is fixed to its time averaged profile, are compared to dynamical runs for which both the Navier-Stokes and the induction equations are jointly solved. The kinematic instability is found to have two branches, for all explored Reynolds numbers. The dynamical dynamo threshold follows these branches: at low Reynolds number it lies within the low branch while at high kinetic Reynolds number it is close to the high branch., Comment: 4 pages, 4 figures

We study the inverse cascade of magnetic helicity in conducting fluids by investigating the detailed transfer of helicity between different spherical shells in Fourier space in direct numerical simulations of three-dimensional magnetohydrodynamics (MHD). Two different numerical simulations are used, one where the system is forced with an electromotive force in the induction equation, and one in which the system is forced mechanically with an ABC flow and the magnetic field is solely sustained by a dynamo action. The magnetic helicity cascade at the initial stages of both simulations is observed to be inverse and local (in scale space) in the large scales, and direct and local in the small scales. When saturation is approached most of the helicity is concentrated in the large scales and the cascade is non-local. Helicity is transfered directly from the forced scales to the largest scales. At the same time, a smaller in amplitude direct cascade is observed from the largest scale to small scales., Comment: Submitted to PRE

We present an object-oriented geophysical and astrophysical spectral-element adaptive refinement (GASpAR) code for application to turbulent flows. Like most spectral-element codes, GASpAR combines finite-element efficiency with spectral-method accuracy. It is also designed to be flexible enough for a range of geophysics and astrophysics applications where turbulence or other complex multiscale problems arise. For extensibility and flexibilty the code is designed in an object-oriented manner. The computational core is based on spectral-element operators, which are represented as objects. The formalism accommodates both conforming and nonconforming elements and their associated data structures for handling interelement communications in a parallel environment. Many aspects of this code are a synthesis of existing methods; however, we focus on a new formulation of dynamic adaptive refinement (DARe) of nonconforming h-type. This paper presents the code and its algorithms; we do not consider parallel efficiency metrics or performance. As a demonstration of the code we offer several two-dimensional test cases that we propose as standard test problems for comparable DARe codes. The suitability of these test problems for turbulent flow simulation is considered., Comment: 51 pages, 12 figures, submitted to J. Comp. Phys., 2005/5/11

We present direct numerical simulations and Lagrangian averaged (also known as alpha-model) simulations of forced and free decaying magnetohydrodynamic turbulence in two dimensions. The statistics of sign cancellations of the current at small scales is studied using both the cancellation exponent and the fractal dimension of the structures. The alpha-model is found to have the same scaling behavior between positive and negative contributions as the direct numerical simulations. The alpha-model is also able to reproduce the time evolution of these quantities in free decaying turbulence. At large Reynolds numbers, an independence of the cancellation exponent with the Reynolds numbers is observed., Comment: Finite size box effects have been taken into account in the definition of the partition function. This has resulted in a more clear scaling in all figures. Several points are clarified in the text

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., Comment: Submitted to PRE

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., Comment: Submitted to PRE

A heuristic model is given for anisotropic magnetohydrodynamics (MHD) turbulence in the presence of a uniform external magnetic field $B_0 {\bf {\hat e}_{\pa}}$. The model is valid for both moderate and strong $B_0$ and is able to describe both the strong and weak wave turbulence regimes as well as the transition between them. The main ingredient of the model is the assumption of constant ratio at all scales between \add{the} linear wave period and \add{the} nonlinear turnover timescale. Contrary to the model of critical balance introduced by Goldreich and Sridhar [P. Goldreich and S. Sridhar, ApJ {\bf 438}, 763 (1995)], it is not assumed in addition that this ratio be equal to unity at all scales which allows us to use the Iroshnikov-Kraichnan phenomenology. It is then possible to recover the widely observed anisotropic scaling law $\kpa \propto \kpe^{2/3}$ between parallel and perpendicular wavenumbers (with reference to $B_0 {\bf {\hat e}_{\pa}}$) and to obtain the universal prediction, $3\alpha + 2\beta = 7$, for the total energy spectrum $E(\kpe,\kpa) \sim \kpe^{-\alpha} \kpa^{-\beta}$. In particular, with such a prediction the weak Alfv\'en wave turbulence constant-flux solution is recovered and, for the first time, a possible explanation to its precursor found numerically by Galtier et al [S. Galtier et al., J. Plasma Phys. {\bf 63}, 447 (2000)] is given.

By re-interpreting a recent successful closure procedure in terms of local spatial averaging and the neglect of fluctuations about that average, it is shown how the results of that closure scheme (the ''alpha model'') may be more simply derived., Comment: To appear in Physics of Fluids (2002)