Source author record

Gregory L. Eyink

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

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

23 published item(s)

preprint2025arXiv

Beyond chaos: fluctuations, anomalies and spontaneous stochasticity in fluid turbulence

In this perspective, we consider the development of statistical hydrodynamics, focusing on the way in which the intrinsic stochasticity of turbulent phenomena was identified and is being explored. A major purpose of our discussion is to bring out the role of anomalies in turbulent phenomena, in ways that are not usually done, and to emphasize how the description of turbulent phenomena requires delicate considerations of asymptotic limits. The scope of our narrative includes selected historical aspects that are not usually emphasized, primarily due to G.I. Taylor, as well as discussions of certain aspects of the laminar-turbulent transition, the behaviour of turbulent drag at intermediate Reynolds numbers, and the statistics of fully-developed turbulence that exhibit spontaneous stochasticity.

preprint2022arXiv

Onsager Theory of Turbulence, the Josephson-Anderson Relation, and the D'Alembert Paradox

The Josephson-Anderson relation, valid for the incompressible Navier-Stokes solutions which describe flow around a solid body, instantaneously equates the power dissipated by drag to the flux of vorticity across the flow lines of the potential Euler solution considered by d'Alembert. Its derivation involves a decomposition of the velocity field into this background potential-flow field and a solenoidal field corresponding to the rotational wake behind the body, with the flux term describing transfer from the interaction energy between the two fields and into kinetic energy of the rotational flow. We establish the validity of the Josephson-Anderson relation for the weak solutions of the Euler equations obtained in the zero-viscosity limit, with one transfer term due to inviscid vorticity flux and the other due to a viscous skin-friction anomaly. Furthermore, we establish weak forms of the local balance equations for both interaction and rotational energies. We define nonlinear spatial fluxes of these energies and show that the asymptotic flux of interaction energy to the wall equals the anomalous skin-friction term in the Josephson-Anderson relation. However, when the Euler solution satisfies suitably the no-flow-through condition at the wall, then the anomalous term vanishes. In this case, we can show also that the asymptotic flux of rotational energy to the wall must vanish and we obtain in the rotational wake the Onsager-Duchon-Robert relation between viscous dissipation anomaly and inertial dissipation due to scale-cascade. In this way we establish a precise connection between the Josephson-Anderson relation and the Onsager theory of turbulence, and we provide a novel resolution of the d'Alembert paradox.

preprint2022arXiv

Origin of Enhanced Skin Friction at the Onset of Boundary-Layer Transition

Boundary-layer transition is accompanied by a significant increase in skin friction whose origin is rigorously explained using the stochastic Lagrangian formulation of the Navier-Stokes equations. This formulation permits the exact analysis of vorticity dynamics in individual realizations of a viscous incompressible fluid flow. The Lagrangian reconstruction formula for vorticity is here extended for the first time to Neumann boundary conditions (Lighthill source). We can thus express the wall vorticity, and therefore the wall stress, as the expectation of a stochastic Cauchy invariant in backward time, with contributions from (a) wall-vorticity flux (Lighthill source) and (b) interior vorticity that has been evolved by nonlinear advection, viscous diffusion, vortex stretching and tilting. We consider the origin of stress maxima in the transitional region, examining a sufficient number of events to represent the increased skin friction. The stochastic Cauchy analysis is applied to each event to trace the origin of the wall vorticity. We find that the Lighthill source, vortex tilting, diffusion and advection of outer vorticity make minor contributions. They are less important than spanwise stretching of near-wall spanwise vorticity, which is the dominant source of skin-friction increase during laminar-to-turbulent transition. Our analysis should assist more generally in understanding drag generation and reduction strategies and flow separation in terms of near-wall vorticity dynamics.

preprint2022arXiv

Path large deviations for the kinetic theory of weak turbulence

We consider a generic Hamiltonian system of nonlinear interacting waves with 3-wave interactions. In the kinetic regime of wave turbulence, which assumes weak nonlinearity and large system size, the relevant observable associated with the wave amplitude is the empirical spectral density that appears as the natural precursor of the spectral density, or spectrum, for finite system size. Following classical derivations of the Peierls equation for the moment generating function of the wave amplitudes in the kinetic regime, we propose a large deviation estimate for the dynamics of the empirical spectral density, where the number of admissible wavenumbers, which is proportional to the volume of the system, appears as the natural large deviation parameter. The large deviation stochastic Hamiltonian that quantifies the minus of the log probability of a trajectory is computed within the kinetic regime which assumes the Random Phase approximation for weak nonlinearity. We compare this Hamiltonian with the one for a system of modes interacting in a mean-field way with the empirical spectrum. Its relationship with the Random Phase and Amplitude approximation is discussed. Moreover, for the specific case when no forces and dissipation are present, a few fundamental properties of the large deviation dynamics are investigated. We show that the latter conserves total energy and momentum, as expected for a 3-wave interacting systems. In addition, we compute the equilibrium quasipotential and check that global detailed balance is satisfied at the large deviation level. Finally, we discuss briefly some physical applications of the theory.

preprint2020arXiv

3D Turbulent Reconnection: Theory, Tests and Astrophysical Implications

Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws.

preprint2015arXiv

Inertial-Range Reconnection in Magnetohydrodynamic Turbulence and in the Solar Wind

In situ spacecraft data on the solar wind show events identified as magnetic reconnection with outflows and apparent "`$X$-lines" $10^{3-4}$ times ion scales. To understand the role of turbulence at these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic (MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard magnetic flux-freezing due to Richardson dispersion, and a broadened reconnection zone containing many current sheets. The coarse-grain magnetic geometry is like large-scale reconnection in the solar wind, however, with a hyperbolic flux-tube or "$X$-line" extending over integral length-scales.

preprint2015arXiv

Quantum Spontaneous Stochasticity

The quantum wave-function of a massive particle with small initial uncertainties (consistent with the uncertainty relation) is believed to spread very slowly, so that the dynamics is deterministic. This assumes that the classical motions for given initial data are unique. In fluid turbulence non-uniqueness due to "roughness" of the advecting velocity field is known to lead to stochastic motion of classical particles. Vanishingly small random perturbations are magnified by Richardson diffusion in a "nearly rough" velocity field so that motion remains stochastic as the noise disappears, or classical spontaneous stochasticity, . Analogies between stochastic particle motion in turbulence and quantum evolution suggest that there should be quantum spontaneous stochasticity (QSS). We show this for 1D models of a particle in a repulsive potential that is "nearly rough" with $V(x) \sim C|x|^{1+α}$ at distances $|x|\gg \ell$ , for some UV cut-off $\ell$, and for initial Gaussian wave-packet centered at 0. We consider the WKB limit with $\hbar/m \to 0$, then position-spread $σ\to 0$. The limit of the position density is non-deterministic, with equal probabilities of the particle following either of two, non-unique classical solutions and with "Richardson-like" spreading of the wave-packet. The splitting of the wave-packet occurs in a very short time $t\sim(\ell^{1+α}/C)^{1/2}$ for small $\ell$. QSS also occurs in the scattering of a wave-packet off the rough potential, for careful fine-tuning. We also consider other semi-classical limits by a numerical solution of the Schrödinger equation, observing QSS both in position- and momentum-space. Although the wave-function remains split into two widely separated branches in the classical limit, rapid phase oscillations within each branch prevent any coherent superposition. These effects should be observable in laboratory experiments.

preprint2015arXiv

Turbulent Reconnection and Its Implications

Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes in magnetized plasmas. In most astrophysical environments the Reynolds numbers are large and therefore the transition to turbulence is inevitable. This turbulence must be taken into account for any theory of magnetic reconnection, since the initially laminar configurations can transit to the turbulence state, what is demonstrated by 3D high resolution numerical simulations. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (1999) reconnection model and present numerical evidence supporting the model and demonstrate that it is closely connected to the concept of Richardson diffusion and compatible with the Lagrangian dynamics of magnetized fluids. We point out that the Generalized Ohm's Law, that accounts for turbulent motion, predicts the subdominance of the microphysical plasma effects for a realistically turbulent media. We show that on of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a corner stone of most theories dealing with magnetized plasmas and therefore its change induces fundamental shifts in accepted paradigms like turbulent reconnection entailing the diffusion process that is essential for understanding star formation. We argue, that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares predicted by turbulent reconnection and relate them to solar flares and gamma ray bursts. We analyze solar observations, measurements in the solar wind or heliospheric current sheet, and show their correspondence with turbulent reconnection predictions. Finally, we discuss 1st Order Fermi acceleration as a natural consequence of the turbulent reconnection.

preprint2014arXiv

Kinematic Morphology of Large-scale Structure: Evolution from Potential to Rotational Flow

As an alternative way of describing the cosmological velocity field, we discuss the evolution of rotational invariants constructed from the velocity gradient tensor. Compared with the traditional divergence-vorticity decomposition, these invariants, defined as coefficients of characteristic equation of the velocity gradient tensor, enable a complete classification of all possible flow patterns in the dark-matter comoving frame, including both potential and vortical flows. We show that this tool, first introduced in turbulence two decades ago, proves to be very useful in understanding the evolution of the cosmic web structure, and in classifying its morphology. Before shell-crossing, different categories of potential flow are highly associated with cosmic web structure, because of the coherent evolution of density and velocity. This correspondence is even preserved at some level when vorticity is generated after shell-crossing. The evolution from the potential to vortical flow can be traced continuously by these invariants. With the help of this tool, we show that the vorticity is generated in a particular way that is highly correlated with the large-scale structure. This includes a distinct spatial distribution and different types of alignment between cosmic web and vorticity direction for various vortical flows. Incorporating shell-crossing into closed dynamical systems is highly non-trivial, but we propose a possible statistical explanation for some of these phenomena relating to the internal structure of the three-dimensional invariants space.

preprint2014arXiv

Turbulent General Magnetic Reconnection

Plasma flows with an MHD-like turbulent inertial range, such as the solar wind, require a generalization of General Magnetic Reconnection (GMR) theory. We introduce the slip-velocity source vector, which gives the rate of development of slip velocity per unit arc length of field line. The slip source vector is the ratio of the curl of the non ideal electric field in the Generalized Ohm's Law and the magnetic field strength. It diverges at magnetic nulls, unifying GMR with magnetic null-point reconnection. Only under restrictive assumptions is the slip velocity related to the gradient of the quasi potential (integral of parallel electric field along field lines). In a turbulent inertial range the curl becomes extremely large while the parallel component is tiny, so that line slippage occurs even while ideal MHD becomes accurate. The resolution of this paradox is that ideal MHD is valid for a turbulent inertial-range only in a weak sense which does not imply magnetic line freezing. The notion of weak solution is explained in terms of spatial coarse-graining and renormalization group (RG) theory. We give an argument for the weak validity of the ideal Ohm's law in the inertial range, via rigorous estimates of the terms in the Generalized Ohm's Law for an electron-ion plasma. All of the nonideal terms (from collisional resistivity, Hall field, electron pressure anisotropy, and electron inertia) are shown to be irrelevant in the RG sense and large-scale reconnection is thus governed solely by ideal dynamics. We briefly discuss some implications for heliospheric reconnection, in particular for deviations from the Parker spiral model of interplanetary magnetic field. Solar wind observations show that reconnection in a turbulence broadened heliospheric current sheet, consistent with the Lazarian-Vishniac theory, leads to slip velocities that cause field lines to lag relative to the spiral model.

preprint2013arXiv

The Diffusion Approximation in Turbulent Two-Particle Dispersion

We solve an inverse problem for fluid particle pair-statistics: we show that a time sequence of probability density functions (PDF's) of separations can be exactly reproduced by solving the diffusion equation with a suitable time-dependent diffusivity. The diffusivity tensor is given by a time-integral of a conditional Lagrangian velocity structure-function, weighted by a ratio of PDF's. Physical hypotheses for hydrodynamic turbulence (sweeping, short memory, mean-field) yield simpler integral formulas, including one of Kraichnan and Lundgren. We evaluate the latter using a spacetime database from a numerical Navier-Stokes solution for driven turbulence. This diffusion theory reproduces PDF's well at rms separations, but growth rate of mean-square dispersion is overpredicted due to neglect of memory effects. More general applications of our approach are sketched.

preprint2012arXiv

Do Sweeping Effects Suppress Particle Dispersion in Synthetic Turbulence?

Synthetic models of Eulerian turbulence like so called "Kinematic Simulations" (KS) have been criticized by Thomson & Devenish (TD) (2005), who argued that sweeping decorrelation effects suppress pair dispersion in such models. We derive analytical results for Eulerian turbulence modeled by Gaussian random fields for the case with zero mean velocity. Our starting point is an exact integrodifferential equation for the particle pair separation distribution. When memory times of particle locations are short, a Markovian approximation leads to a Richardson-type diffusion model. We obtain a diffusivity tensor of the form $K_{ij}(r,t)=S_{ij}(r)τ(r,t)$ where $S_{ij}(r)$ is the structure-function tensor and $τ(r,t)$ is an effective correlation time of velocity increments. This is found to be the minimum value of three times: the intrinsic turnover time $τ_{eddy}(r)$ at separation $r$, the overall evolution time $t,$ and the sweeping time $r/v_0$ with $v_0$ the rms velocity. We study the diffusion model numerically by a Monte Carlo method. With moderate inertial-ranges like those achieved in current KS, our model is found to reproduce the $t^{9/2}$ power-law for pair dispersion predicted by TD and observed in the KS. However, for much longer ranges, our model exhibits three distinct pair-dispersion laws in the inertial-range: a Batchelor $t^2$-regime, followed by a $t^1$ diffusive regime, and then a $t^6$ regime. Finally, outside the inertial-range, there is another $t^1$ regime with particles undergoing independent Taylor diffusion. These scalings are exactly the same as those predicted by TD for KS with large mean velocities, which we argue hold also for KS with zero mean velocity. Our results support the basic conclusion that sweeping effects make Lagrangian properties of KS fundamentally different from hydrodynamic turbulence for very extended inertial-ranges.

preprint2012arXiv

Kinetic Wave Turbulence

We consider a general model of Hamiltonian wave systems with triple resonances, in the standard kinetic limit of a continuum of weakly interacting dispersive waves with random phases. In this asymptotic limit we show that the correct dynamical equation for multimode amplitude distributions is not the well-known equation of Peierls but is instead a reduced equation with only a subset of the terms in that equation. The equations that we derive are the direct analogue of the Boltzmann hierarchy obtained from the BBGKY hierarchy in the low-density limit for gases. We show that the asymptotic multimode equations possess factorized solutions for factorized initial data, which correspond to preservation in time of the property of "random phases & amplitudes". The factors satisfy the equations for the 1-mode probability density functions previously derived by Jakobsen & Newell and Choi et al. We show that the factorization of the hierarchy equations implies that these quantities are self-averaging: they satisfy the wave-kinetic closure equations of the spectrum and 1-mode PDF for almost any selection of phases and amplitudes from the initial ensemble. We show that both of these closure equations satisfy an H-theorem for an entropy defined by Boltzmann's prescription. We also characterize the general solutions of our multimode distribution equations, for initial conditions with random phases but with no assumptions on the amplitudes. Analogous to a result of Spohn for the Boltzmann hierarchy, these are "super-statistical solutions" that correspond to ensembles of solutions of the wave-kinetic closure equations with random initial conditions or random forces. On the basis of our results, we discuss possible kinetic explanations of intermittency and non-Gaussian statistics in wave turbulence.

preprint2011arXiv

Fast Magnetic Reconnection and Spontaneous Stochasticity

Magnetic field-lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion gyroradius. The rapid reconnection of magnetic flux structures with dimensions vastly larger than the gyroradius requires a breakdown in the standard Alfvén flux-freezing law. We attribute this breakdown to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic energy spectra. Lagrangian particle trajectories in such environments become "spontaneously stochastic", so that infinitely-many magnetic field-lines are advected to each point and must be averaged to obtain the resultant magnetic field. The relative distance between initial magnetic field lines which arrive to the same final point depends upon the properties of two-particle turbulent dispersion. We develop predictions based on the phenomenological Goldreich & Sridhar theory of strong MHD turbulence and on weak MHD turbulence theory. We recover the predictions of the Lazarian & Vishniac theory for the reconnection rate of large-scale magnetic structures. Lazarian & Vishniac also invoked "spontaneous stochasticity", but of the field-lines rather than of the Lagrangian trajectories. More recent theories of fast magnetic reconnection appeal to microscopic plasma processes that lead to additional terms in the generalized Ohm's law, such as the collisionless Hall term. We estimate quantitatively the effect of such processes on the inertial-range turbulence dynamics and find them to be negligible in most astrophysical environments. For example, the predictions of the Lazarian-Vishniac theory are unchanged in Hall MHD turbulence with an extended inertial range, whenever the ion skin depth $δ_i$ is much smaller than the turbulent integral length or injection-scale $L_i.$

preprint2010arXiv

Fluctuation Dynamo and Turbulent Induction at Small Prandtl Number

We study the Lagrangian mechanism of the fluctuation dynamo at zero Prandtl number and infinite magnetic Reynolds number, in the Kazantsev-Kraichnan model of white-noise advection. With a rough velocity field corresponding to a turbulent inertial-range, flux-freezing holds only in a stochastic sense. We show that field-lines arriving to the same point which were initially separated by many resistive lengths are important to the dynamo. Magnetic vectors of the seed field that point parallel to the initial separation vector arrive anti-correlated and produce an "anti-dynamo" effect. We also study the problem of "magnetic induction" of a spatially uniform seed field. We find no essential distinction between this process and fluctuation dynamo, both producing the same growth-rates and small-scale magnetic correlations. In the regime of very rough velocity fields where fluctuation dynamo fails, we obtain the induced magnetic energy spectra. We use these results to evaluate theories proposed for magnetic spectra in laboratory experiments of turbulent induction

preprint2010arXiv

Stochastic Flux-Freezing and Magnetic Dynamo

We argue that magnetic flux-conservation in turbulent plasmas at high magnetic Reynolds numbers neither holds in the conventional sense nor is entirely broken, but instead is valid in a novel statistical sense associated to the "spontaneous stochasticity" of Lagrangian particle tra jectories. The latter phenomenon is due to the explosive separation of particles undergoing turbulent Richardson diffusion, which leads to a breakdown of Laplacian determinism for classical dynamics. We discuss empirical evidence for spontaneous stochasticity, including our own new numerical results. We then use a Lagrangian path-integral approach to establish stochastic flux-freezing for resistive hydromagnetic equations and to argue, based on the properties of Richardson diffusion, that flux-conservation must remain stochastic at infinite magnetic Reynolds number. As an important application of these results we consider the kinematic, fluctuation dynamo in non-helical, incompressible turbulence at unit magnetic Prandtl number. We present results on the Lagrangian dynamo mechanisms by a stochastic particle method which demonstrate a strong similarity between the Pr = 1 and Pr = 0 dynamos. Stochasticity of field-line motion is an essential ingredient of both. We finally consider briefly some consequences for nonlinear MHD turbulence, dynamo and reconnection

preprint2009arXiv

Localness of energy cascade in hydrodynamic turbulence, I. Smooth coarse-graining

We introduce a novel approach to scale-decomposition of the fluid kinetic energy (or other quadratic integrals) into band-pass contributions from a series of length-scales. Our decomposition is based on a multiscale generalization of the ``Germano identity'' for smooth, graded filter kernels. We employ this method to derive a budget equation that describes the transfers of turbulent kinetic energy both in space and in scale. It is shown that the inter-scale energy transfer is dominated by local triadic interactions, assuming only the scaling properties expected in a turbulent inertial-range. We derive rigorous upper bounds on the contributions of non-local triads, extending the work of Eyink (2005) for low-pass filtering. We also propose a physical explanation of the differing exponents for our rigorous upper bounds and for the scaling predictions of Kraichnan (1966,1971). The faster decay predicted by Kraichnan is argued to be the consequence of additional cancellations in the signed contributions to transfer from non-local triads, after averaging over space. This picture is supported by data from a $512^3$ pseudospectral simulation of Navier-Stokes turbulence with phase-shift dealiasing.

preprint2009arXiv

Localness of energy cascade in hydrodynamic turbulence, II. Sharp spectral filter

We investigate the scale-locality of subgrid-scale (SGS) energy flux and inter-band energy transfers defined by the sharp spectral filter. We show by rigorous bounds, physical arguments and numerical simulations that the spectral SGS flux is dominated by local triadic interactions in an extended turbulent inertial-range. Inter-band energy transfers are also shown to be dominated by local triads if the spectral bands have constant width on a logarithmic scale. We disprove in particular an alternative picture of ``local transfer by nonlocal triads,'' with the advecting wavenumber mode at the energy peak. Although such triads have the largest transfer rates of all {\it individual} wavenumber triads, we show rigorously that, due to their restricted number, they make an asymptotically negligible contribution to energy flux and log-banded energy transfers at high wavenumbers in the inertial-range. We show that it is only the aggregate effect of a geometrically increasing number of local wavenumber triads which can sustain an energy cascade to small scales. Furthermore, non-local triads are argued to contribute even less to the space-average energy flux than is implied by our rigorous bounds, because of additional cancellations from scale-decorrelation effects. We can thus recover the -4/3 scaling of nonlocal contributions to spectral energy flux predicted by Kraichnan's ALHDIA and TFM closures. We support our results with numerical data from a $512^3$ pseudospectral simulation of isotropic turbulence with phase-shift dealiasing. We conclude that the sharp spectral filter has a firm theoretical basis for use in large-eddy simulation (LES) modeling of turbulent flows.

preprint2009arXiv

Scale-Locality of Magnetohydrodynamic Turbulence

We investigate the scale-locality of cascades of conserved invariants at high kinetic and magnetic Reynolds numbers in the ``inertial-inductive range'' of magnetohydrodynamic (MHD) turbulence, where velocity and magnetic field increments exhibit suitable power-law scaling. We prove that fluxes of total energy and cross-helicity--or, equivalently, fluxes of Elsässer energies--are dominated by the contributions of local triads. Corresponding spectral transfers are also scale-local when defined using octave wavenumber bands. Flux and transfer of magnetic helicity may be dominated by non-local triads. The magnetic stretching term also may be dominated by non-local triads but we prove that it can convert energy only between velocity and magnetic modes at comparable scales. We explain the disagreement with numerical studies that have claimed conversion nonlocally between disparate scales. We present supporting data from a $1024^3$ simulation of forced MHD turbulence.

preprint2009arXiv

Small-Scale Kinematic Dynamo and Non-Dynamo in Inertial-Range Turbulence

We investigate the Lagrangian mechanism of the kinematic ``fluctuation'' magnetic dynamo in turbulent plasma flow at small magnetic Prandtl numbers. The combined effect of turbulent advection and plasma resistivity is to carry infinitely many field lines to each space point, with the resultant magnetic field at that point given by the average over all the individual line vectors. As a consequence of the roughness of the advecting velocity, this remains true even in the limit of zero resistivity. We show that the presence of dynamo effect requires sufficient angular correlation of the passive line-vectors that arrive simultaneously at the SAME space point. We demonstrate this in detail for the Kazantsev-Kraichnan model of kinematic dynamo with a Gaussian advecting velocity that is spatially rough and white-noise in time. In the regime where dynamo action fails, we also obtain the precise rate of decay of the magnetic energy. These exact results for the model are obtained by a generalization of the ``slow-mode expansion'' of Bernard, Gawedzki and Kupiainen to non-Hermitian evolution. Much of our analysis applies also to magnetohydrodynamic turbulence.

preprint2008arXiv

Stochastic Line-Motion and Stochastic Conservation Laws for Non-Ideal Hydromagnetic Models. I. Incompressible Fluids and Isotropic Transport Coefficients

We prove that smooth solutions of non-ideal (viscous and resistive) incompressible magnetohydrodynamic equations satisfy a stochastic law of flux conservation. This property involves an ensemble of surfaces obtained from a given, fixed surface by advecting it backward in time under the plasma velocity perturbed with a random white-noise. It is shown that the magnetic flux through the fixed surface is equal to the average of the magnetic fluxes through the ensemble of surfaces at earlier times. This result is an analogue of the well-known Alfven theorem of ideal MHD and is valid for any value of the magnetic Prandtl number. A second stochastic conservation law is shown to hold at unit Prandtl number, a random version of the generalized Kelvin theorem derived by Bekenstein-Oron for ideal MHD. These stochastic conservation laws are not only shown to be consequences of the non-ideal MHD equations, but are proved in fact to be equivalent to those equations. We derive similar results for two more refined hydromagnetic models, Hall magnetohydrodynamics and the two-fluid plasma model, still assuming incompressible velocities and isotropic transport coefficients. Finally, we use these results to discuss briefly the infinite-Reynolds-number limit of hydromagnetic turbulence and to support the conjecture that flux-conservation remains stochastic in that limit.

preprint2006arXiv

The Breakdown of Alfven's Theorem in Ideal Plasma Flows

This paper presents both rigorous results and physical theory on the breakdown of magnetic flux conservation for ideal plasmas, by nonlinear effects. Our analysis is based upon an effective equation for magnetohydrodynamic (MHD) modes at length-scales $>\ell,$ with smaller scales eliminated, as in renormalization-group methodology. We prove that flux-conservation can be violated for an arbitrarily small length-scale $\ell,$ and in the absence of any non-ideality, but only if singular current sheets and vortex sheets both exist and intersect in sets of large enough dimension. This result gives analytical support to and rigorous constraints on theories of fast turbulent reconnection. Mathematically, our theorem is analogous to Onsager's result on energy dissipation anomaly in hydrodynamic turbulence. As a physical phenomenon, the breakdown of magnetic-flux conservation in ideal MHD is similar to the decay of magnetic flux through a narrow superconducting ring, by phase-slip of quantized flux lines. The effect should be observable both in numerical MHD simulations and in laboratory plasma experiments at moderately high magnetic Reynolds numbers.

preprint2003arXiv

Recovering Isotropic Statistics in Turbulence Simulations: The Kolmogorov 4/5th-Law

One of the main benchmarks in direct numerical simulations of three-dimensional turbulence is the Kolmogorov 1941 prediction for third-order structure functions with homogeneous and isotropic statistics in the infinite-Reynolds number limit. Previous DNS techniques to obtain isotropic statistics have relied on time-averaging structure functions in a few directions over many eddy turnover times, using forcing schemes carefully constructed to generate isotropic data. Motivated by recent theoretical work which removes isotropy requirements by spherically averaging structure functions over all directions, we will present results which supplement long-time averaging by angle-averaging over up to 73 directions from a single flow snapshot. The directions are among those natural to a square computational grid, and are weighted to approximate the spherical average. The averaging process is cheap, and for the Kolmogorov 1941 4/5ths law, reasonable results can be obtained from a single snapshot of data. This procedure may be used to investigate the isotropic statistics of any quantity of interest.