Minping Wan

A public turbulence database cluster and applications to study Lagrangian evolution of velocity increments in turbulence

A public database system archiving a direct numerical simulation (DNS) data set of isotropic, forced turbulence is described in this paper. The data set consists of the DNS output on 10243 spatial points and 1024 time samples spanning about one large-scale turnover time. This complete 10244 spacetime history of turbulence is accessible to users remotely through an interface that is based on the Web-services model. Users may write and execute analysis programs on their host computers, while the programs make subroutine-like calls that request desired parts of the data over the network. The users are thus able to perform numerical experiments by accessing the 27 terabytes (TB) of DNS data using regular platforms such as laptops. The architecture of the database is explained, as are some of the locally defined functions, such as differentiation and interpolation. Test calculations are performed to illustrate the usage of the system and to verify the accuracy of the methods. The database is then used to analyze a dynamical model for small-scale intermittency in turbulence. Specifically, the dynamical effects of pressure and viscous terms on the Lagrangian evolution of velocity increments are evaluated using conditional averages calculated from the DNS data in the database. It is shown that these effects differ considerably among themselves and thus require different modeling strategies in Lagrangian models of velocity increments and intermittency.

Coherent structures, intermittent turbulence, and dissipation in high-temperature plasmas

An unsolved problem in plasma turbulence is how energy is dissipated at small scales. Particle collisions are too infrequent in hot plasmas to provide the necessary dissipation. Simulations either treat the fluid scales and impose an ad hoc form of dissipation (e.g., resistivity) or consider dissipation arising from resonant damping of small amplitude disturbances where damping rates are found to be comparable to that predicted from linear theory. Here, we report kinetic simulations that span the macroscopic fluid scales down to the motion of electrons. We find that turbulent cascade leads to generation of coherent structures in the form of current sheets that steepen to electron scales, triggering strong localized heating of the plasma. The dominant heating mechanism is due to parallel electric fields associated with the current sheets, leading to anisotropic electron and ion distributions which can be measured with NASAs upcoming Magnetospheric Multiscale mission. The motion of coherent structures also generates waves that are emitted into the ambient plasma in form of highly oblique compressional and shear Alfven modes. In 3D, modes propagating at other angles can also be generated. This indicates that intermittent plasma turbulence will in general consist of both coherent structures and waves. However, the current sheet heating is found to be locally several orders of magnitude more efficient than wave damping and is sufficient to explain the observed heating rates in the solar wind.

The link between shocks, turbulence, and magnetic reconnection in collisionless plasmas

Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence, and magnetic reconnection. In particular, collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These “wavefronts” maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, “stir” the downstream region, creating large scale disturbances ...

INTERMITTENT HEATING IN SOLAR WIND AND KINETIC SIMULATIONS

Low-density astrophysical plasmas may be described by magnetohydrodynamics at large scales, but require kinetic description at ion scales in order to include dissipative processes that terminate the cascade. Here kinetic plasma simulations and high-resolution spacecraft observations are compared to facilitate the interpretation of signatures of various dissipation mechanisms. Kurtosis of increments indicates that kinetic scale coherent structures are present, with some suggestion of incoherent activity near ion scales. Conditioned proton temperature distributions suggest heating associated with coherent structures. The results reinforce the association of intermittent turbulence, coherent structures, and plasma dissipation.

Intermittency and local heating in the solar wind

Evidence for nonuniform heating in the solar wind plasma near current sheets dynamically generated by magnetohydrodynamic (MHD) turbulence is obtained using measurements from the ACE spacecraft. These coherent structures only constitute 19% of the data, but contribute 50% of the total plasma internal energy. Intermittent heating manifests as elevations in proton temperature near current sheets, resulting in regional heating and temperature enhancements extending over several hours. The number density of non-Gaussian structures is found to be proportional to the mean proton temperature and solar wind speed. These results suggest magnetofluid turbulence drives intermittent dissipation through a hierarchy of coherent structures, which collectively could be a significant source of coronal and solar wind heating.

Intermittency, nonlinear dynamics and dissipation in the solar wind and astrophysical plasmas

An overview is given of important properties of spatial and temporal intermittency, including evidence of its appearance in fluids, magnetofluids and plasmas, and its implications for understanding of heliospheric plasmas. Spatial intermittency is generally associated with formation of sharp gradients and coherent structures. The basic physics of structure generation is ideal, but when dissipation is present it is usually concentrated in regions of strong gradients. This essential feature of spatial intermittency in fluids has been shown recently to carry over to the realm of kinetic plasma, where the dissipation function is not known from first principles. Spatial structures produced in intermittent plasma influence dissipation, heating, and transport and acceleration of charged particles. Temporal intermittency can give rise to very long time correlations or a delayed approach to steady-state conditions, and has been associated with inverse cascade or quasi-inverse cascade systems, with possible implications for heliospheric prediction.

LOCAL ANISOTROPY, HIGHER ORDER STATISTICS, AND TURBULENCE SPECTRA

Correlation anisotropy emerges dynamically in magnetohydrodynamics (MHD), producing stronger gradients across the large-scale mean magnetic field than along it. This occurs both globally and locally, and has significant implications in space and astrophysical plasmas, including particle scattering and transport, and theories of turbulence. Properties of local correlation anisotropy are further documented here by showing through numerical experiments that the effect is intensified in more localized estimates of the mean field. The mathematical formulation of this property shows that local anisotropy mixes second-order with higher order correlations. Sensitivity of local statistical estimates to higher order correlations can be understood in connection with the stochastic coordinate system inherent in such formulations. We demonstrate this in specific cases, and illustrate the connection to higher order statistics by showing the sensitivity of local anisotropy to phase randomization, after which the global measure of anisotropy is recovered at all scales of averaging. This establishes that anisotropy of the local structure function is not a measure of anisotropy of the energy spectrum. Evidently, the local enhancement of correlation anisotropy is of substantial fundamental interest and must be understood in terms of higher order correlations, specifically fourth-order and above.

On the accuracy of simulations of turbulence

The widely recognized issue of adequate spatial resolution in numerical simulations of turbulence is studied in the context of two-dimensional magnetohydrodynamics. The familiar criterion that the dissipation scale should be resolved enables accurate computation of the spectrum, but fails for precise determination of higher-order statistical quantities. Examination of two straightforward diagnostics, the maximum of the kurtosis and the scale-dependent kurtosis, enables the development of simple tests for assessing adequacy of spatial resolution. The efficacy of the tests is confirmed by examining a sample problem, the distribution of magnetic reconnection rates in turbulence. Oversampling the Kolmogorov dissipation scale by a factor of 3 allows accurate computation of the kurtosis, the scale-dependent kurtosis, and the reconnection rates. These tests may provide useful guidance for resolution requirements in many plasma computations involving turbulence and reconnection.

Physical mechanism of the inverse energy cascade of two-dimensional turbulence: a numerical investigation

We report an investigation of inverse energy cascade in steady-state two-dimensional turbulence by direct numerical simulation (DNS) of the two-dimensional Navier–Stokes equation, with small-scale forcing and large-scale damping. We employed several types of damping and dissipation mechanisms in simulations up to 2048 2 resolution. For all these simulations we obtained a wavenumber range for which the mean spectral energy flux is a negative constant and the energy spectrum scales as k −5/3 , consistent with the predictions of Kraichnan ( Phys. Fluids , vol. 439, 1967, p. 1417). To gain further insight, we investigated the energy cascade in physical space, employing a local energy flux defined by smooth filtering. We found that the inverse energy cascade is scale local, but that the strongly local contribution vanishes identically, as argued by Kraichnan ( J. Fluid Mech ., vol. 47, 1971, p. 525). The mean flux across a length scale l was shown to be due mainly to interactions with modes two to eight times smaller. A major part of our investigation was devoted to identifying the physical mechanism of the two-dimensional inverse energy cascade. One popular idea is that inverse energy cascade proceeds via merger of like-sign vortices. We made a quantitative study employing a precise topological criterion of merger events. Our statistical analysis showed that vortex mergers play a negligible direct role in producing mean inverse energy flux in our simulations. Instead, we obtained with the help of other works considerable evidence in favour of a ‘vortex thinning’ mechanism, according to which the large-scale strains do negative work against turbulent stress as they stretch out the isolines of small-scale vorticity. In particular, we studied a multi-scale gradient (MSG) expansion developed by Eyink ( J. Fluid Mech ., vol. 549, 2006 a , p. 159) for the turbulent stress, whose contributions to inverse cascade can all be explained by ‘thinning’. The MSG expansion up to second order in space gradients was found to predict well the magnitude, spatial structure and scale distribution of the local energy flux. The majority of mean flux was found to be due to the relative rotation of strain matrices at different length scales, a first-order effect of ‘thinning’. The remainder arose from two second-order effects, differential strain rotation and vorticity gradient stretching. Our findings give strong support to vortex thinning as the fundamental mechanism of two-dimensional inverse energy cascade.

ASSOCIATION OF SUPRATHERMAL PARTICLES WITH COHERENT STRUCTURES AND SHOCKS

Various mechanisms have been proposed to explain observed suprathermal particle populations in the solar wind, including direct acceleration at flares, stochastic acceleration, shock acceleration, and acceleration by random compression or reconnection sites. Using magnetic field and suprathermal particle data from the Advanced Composition Explorer (ACE), we identify coherent structures and interplanetary shocks, and analyze the temporal association of energetic particle fluxes with these coherent structures. Coherent structures having a range of intensities are identified using the magnetic Partial Variance of Increments statistic, essentially a normalized vector increment. A stronger association of energetic particle flux in the 0.047-4.75 MeV range is found with intense magnetic discontinuities than is found with shocks. Nevertheless, the average profile of suprathermals near shocks is quite consistent with standard models of diffusive shock acceleration, while a significant amount of the energetic particles measured and strong discontinuities are found by ACE within six hours of a shock. This evidence supports the view that multiple mechanisms contribute to the acceleration and transport of interplanetary suprathermal particles.