Cristian Constantin Lalescu

Flux-freezing breakdown in high-conductivity magnetohydrodynamic turbulence

The idea of ‘frozen-in’ magnetic field lines for ideal plasmas is useful to explain diverse astrophysical phenomena, for example the shedding of excess angular momentum from protostars by twisting of field lines frozen into the interstellar medium. Frozen-in field lines, however, preclude the rapid changes in magnetic topology observed at high conductivities, as in solar flares. Microphysical plasma processes are a proposed explanation of the observed high rates, but it is an open question whether such processes can rapidly reconnect astrophysical flux structures much greater in extent than several thousand ion gyroradii. An alternative explanation is that turbulent Richardson advection brings field lines implosively together from distances far apart to separations of the order of gyroradii. Here we report an analysis of a simulation of magnetohydrodynamic turbulence at high conductivity that exhibits Richardson dispersion. This effect of advection in rough velocity fields, which appear non-differentiable in space, leads to line motions that are completely indeterministic or ‘spontaneously stochastic’, as predicted in analytical studies. The turbulent breakdown of standard flux freezing at scales greater than the ion gyroradius can explain fast reconnection of very large-scale flux structures, both observed (solar flares and coronal mass ejections) and predicted (the inner heliosheath, accretion disks, γ-ray bursts and so on). For laminar plasma flows with smooth velocity fields or for low turbulence intensity, stochastic flux freezing reduces to the usual frozen-in condition.

A Web services accessible database of turbulent channel flow and its use for testing a new integral wall model for LES

abstract The output from a direct numerical simulation (DNS) of turbulent channel flow at Reτ ≈ 1000 is used to construct a publicly and Web services accessible, spatio-temporal database for this flow. The simulated channel has a size of 8πh × 2h × 3πh, where h is the channel half-height. Data are stored at 2048 × 512 × 1536 spatial grid points for a total of 4000 time samples every 5 time steps of the DNS. These cover an entire channel flow-through time, i.e. the time it takes to traverse the entire channel length 8πh at the mean velocity of the bulk flow. Users can access the database through an interface that is based on the Web services model and perform numerical experiments on the slightly over 100 terabytes (TB) DNS data on their remote platforms, such as laptops or local desktops. Additional technical details about the pressure calculation, database interpolation, and differentiation tools are provided in several appendices. As a sample application of the channel flow database, we use it to conduct an a-priori test of a recently introduced integral wall model for large eddy simulation of wall-bounded turbulent flow. The results are compared with those of the equilibrium wall model, showing the strengths of the integral wall model as compared to the equilibrium model.

Synchronization of chaos in fully developed turbulence.

We investigate chaos synchronization of small-scale motions in the three-dimensional turbulent energy cascade, via pseudospectral simulations of the incompressible Navier-Stokes equations. The modes of the turbulent velocity field below about 20 Kolmogorov dissipation lengths are found to be slaved to the chaotic dynamics of larger-scale modes. The dynamics of all dissipation-range modes can be recovered to full numerical precision by solving small-scale dynamical equations with the given large-scale solution as an input, regardless of initial condition. The synchronization rate exponent scales with the Kolmogorov dissipation time scale, with possible weak corrections due to intermittency. Our results suggest that all sub-Kolmogorov length modes should be fully recoverable from numerical simulations with standard, Kolmogorov-length grid resolutions.

Implementation of high order spline interpolations for tracking test particles in discretized fields

A systematic approach for constructing high order spline interpolation methods is proposed for fields known on regular, rectangular grids. These interpolation methods are tested in tracking trajectories of particles submitted to a force that derives from a potential known on a grid. The interplay between the time advancement scheme and the spatial interpolation is studied in detail and it is shown how the order of the trajectory solver is directly affected by the order of the spline interpolation. It is also shown how an interpolation method that preserves topological properties of physical fields can be better exploited with these higher order spline approximations.

How tracer particles sample the complexity of turbulence

On their roller coaster ride through turbulence, tracer particles sample the fluctuations of the underlying fields in space and time. Quantitatively relating particle and field statistics remains a fundamental challenge in a large variety of turbulent flows. We quantify how tracer particles sample turbulence by expressing their temporal velocity fluctuations in terms of an effective probabilistic sampling of spatial velocity field fluctuations. To corroborate our theory, we investigate an extensive suite of direct numerical simulations of hydrodynamic turbulence covering a Taylor-scale Reynolds number range from 150 to 430. Our approach allows the assessment of particle statistics from the knowledge of flow field statistics only, therefore opening avenues to a new generation of models for transport in complex flows.

Influence of numerical schemes on statistical properties of computed charged particle trajectories in turbulent electromagnetic fields

A class of numerical schemes is developed for the study of charged particle transport in complex stationary electromagnetic fields and tested for fields obtained from a numerical solution of the magneto-hydrodynamic equation. The performances of these schemes are evaluated by analyzing the conservation of energy and the statistical properties of the trajectories. Energy conservation is affected by the interpolation technique used to estimate the field value at the particle position. However, the particle transport properties are more robust, except in the limit of low energy when a significant fraction of the particles are trapped.

Particle diffusion in prescribed electrostatic turbulence and sheared space dependent magnetic field

This study is devoted to the calculation of diffusion coefficients for a particle moving in fluctuating electrostatic field superposed to a space dependent and sheared magnetic field, using the numerical simulation.

Particle transport in incompressible MHD Kolmogorov flow

In previous work [1], an investigation of Kolmogorov flow for incompressible magnetohydrodynamics (MHD) was performed. It consists of a three-dimensional periodic flow driven by a unidirectional forcing varying on a transverse direction. In practice, the forcing is chosen to be fx = A sin(kfy), and fy = fz = 0. It was found that vorticity structures with long lifetimes can be formed in the turbulent regime. The presence of such structures may affect the transport of particles interacting with the flow as shown in this preliminary study.