Engineering Papers⌕ Search

Engineering topics

Beresnyak, Andrey

Publications and source records attributed to Beresnyak, Andrey.

Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma

The Naval Research Laboratory ‘Mag Noh problem’, described in this paper, is a self-similar magnetized implosion flow, which contains a fast magnetohydrodynamic (MHD) outward propagating shock of constant velocity. We generalize the classic Noh ( OSTI Tech. Rep. 577058, 1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$ . Our self-similar solutions are indeed realized when we solve the initial value MHD problem with the finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field and a shock wave.

Mechanics↗

Three-dimensional spontaneous magnetic reconnection

Magnetic reconnection is best known from observations of the Sun where it causes solar flares. Observations estimate the reconnection rate as a small, but non-negligible fraction of the Alfvén speed, so-called fast reconnection. Until recently, the prevailing pictures of reconnection were either of resistivity or plasma microscopic effects, which was contradictory to the observed rates. Alternative pictures were either of reconnection due to the stochasticity of magnetic field lines in turbulence or the tearing instability of the thin current sheet. In this paper we simulate long-term three-dimensional nonlinear evolution of a thin, planar current sheet subject to a fast oblique tearing instability using direct numerical simulations of resistive-viscous magnetohydrodynamics. The late-time evolution resembles generic turbulence with a −5/3 power spectrum and scale-dependent anisotropy, so we conclude that the tearing-driven reconnection becomes turbulent reconnection. The turbulence is local in scale, so microscopic diffusivity should not affect large-scale quantities. This is confirmed by convergence of the reconnection rate toward ∼0.015v{sub A} with increasing Lundquist number. In this spontaneous reconnection, with mean field and without driving, the dissipation rate per unit area also converges to ∼0.006ρv{sub A}{sup 3}, and the dimensionless constants 0.015 and 0.006 are governed only by self-driven nonlinear dynamics of the sheared magnetic field. Remarkably, this also means that a thin current sheet has a universal fluid resistance depending only on its length to width ratio and to v{sub A}/c.

79 ASTRONOMY AND ASTROPHYSICS↗

FIRST-ORDER PARTICLE ACCELERATION IN MAGNETICALLY DRIVEN FLOWS

We demonstrate that particles are regularly accelerated while experiencing curvature drift in flows driven by magnetic tension. Some examples of such flows include spontaneous turbulent reconnection and decaying magnetohydrodynamic turbulence, where a magnetic field relaxes to a lower-energy configuration and transfers part of its energy to kinetic motions of the fluid. We show that this energy transfer, which normally causes turbulent cascade and heating of the fluid, also results in a first-order acceleration of non-thermal particles. Since it is generic, this acceleration mechanism is likely to play a role in the production of non-thermal particle distribution in magnetically dominant environments such as the solar chromosphere, pulsar magnetospheres, jets from supermassive black holes, and γ-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

TURBULENT AMPLIFICATION AND STRUCTURE OF THE INTRACLUSTER MAGNETIC FIELD

We compare DNS calculations of homogeneous isotropic turbulence with the statistical properties of intracluster turbulence from the Matryoshka Run and find remarkable similarities between their inertial ranges. This allowed us to use the time-dependent statistical properties of intracluster turbulence to evaluate dynamo action in the intracluster medium, based on earlier results from a numerically resolved nonlinear magneto-hydrodynamic turbulent dynamo. We argue that this approach is necessary (a) to properly normalize dynamo action to the available intracluster turbulent energy and (b) to overcome the limitations of low Re affecting current numerical models of the intracluster medium. We find that while the properties of intracluster magnetic field are largely insensitive to the value and origin of the seed field, the resulting values for the Alfvén speed and the outer scale of the magnetic field are consistent with current observational estimates, basically confirming the idea that the magnetic field in today’s galaxy clusters is a record of its past turbulent activity.

79 ASTRONOMY AND ASTROPHYSICS↗

Acceleration during magnetic reconnection

The presentation begins with colorful depictions of solar x-ray flares and references to pulsar phenomena. Plasma reconnection is complex, could be x-point dominated or turbulent, field lines could break due to either resistivity or non-ideal effects, such as electron pressure anisotropy. Electron acceleration is sometimes observed, and sometimes not. One way to study this complex problem is to have many examples of the process (reconnection) and compare them; the other way is to simplify and come to something robust. Ideal MHD (E=0) turbulence driven by magnetic energy is assumed, and the first-order acceleration is sought. It is found that dissipation in big (length >100 ion skin depths) current sheets is universal and independent on microscopic resistivity and the mean imposed field; particles are regularly accelerated while experiencing curvature drift in flows driven by magnetic tension. One example of such flow is spontaneous reconnection. This explains hot electrons with a power-law tail in solar flares, as well as ultrashort time variability in some astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗

THE FLOW AROUND A COSMIC STRING. I. HYDRODYNAMIC SOLUTION

Cosmic strings are linear topological defects which are hypothesized to be produced during inflation. Most searches for strings have relied on the string’s lensing of background galaxies or the cosmic microwave background. In this paper, I obtained a solution for the supersonic flow of collisional gas past the cosmic string which has two planar shocks with a shock compression ratio that depends on the angle defect of the string and its speed. The shocks result in the compression and heating of the gas and, given favorable conditions, particle acceleration. Gas heating and over-density in an unusual wedge shape can be detected by observing the Hi line at high redshifts. Particle acceleration can occur in the present-day universe when the string crosses the hot gas contained in galaxy clusters and, since the consequences of such a collision persist for cosmological timescales, could be located by looking at unusual large-scale radio sources situated on a single spatial plane.

79 ASTRONOMY AND ASTROPHYSICS↗

Basic Bell-MHD turbulence

Nonresonant current instability was identified by Bell as an important mechanism for magnetic field amplification in supernova remnants. In this paper we focus on studying the nonlinear stage of this instability using the incompressible MHD formulation. We demonstrate that the evolution of magnetic turbulence driven by the Bell instability resembles turbulence driven on large scales. More importantly, we demonstrate that the energy-containing scale for magnetic fields is proportional to the square root of the magnetic energy density. Given the observational constraints of the possible field amplification, this new relation allows us to directly estimate the maximum energy of particles scattered by such fields, and this estimate is normally below the average particle energy. This implies that, without taking into account the feedback to cosmic rays (CRs), the typical scales of Bell fields, in either the linear or nonlinear regime, will be too small to affect high-energy particle acceleration. We mention several scenarios of back reaction to CRs that could be important.

79 ASTRONOMY AND ASTROPHYSICS↗

SPECTRA OF STRONG MAGNETOHYDRODYNAMIC TURBULENCE FROM HIGH-RESOLUTION SIMULATIONS

Magnetohydrodynamic (MHD) turbulence is present in a variety of solar and astrophysical environments. Solar wind fluctuations with frequencies lower than 0.1 Hz are believed to be mostly governed by Alfvénic turbulence with particle transport depending on the power spectrum and the anisotropy of such turbulence. Recently, conflicting spectral slopes for the inertial range of MHD turbulence have been reported by different groups. Spectral shapes from earlier simulations showed that MHD turbulence is less scale-local compared with hydrodynamic turbulence. This is why higher-resolution simulations, and careful and rigorous numerical analysis is especially needed for the MHD case. In this Letter, we present two groups of simulations with resolution up to 4096{sup 3}, which are numerically well-resolved and have been analyzed with an exact and well-tested method of scaling study. Our results from both simulation groups indicate that the asymptotic power spectral slope for all energy-related quantities, such as total energy and residual energy, is around –1.7, close to Kolmogorov's –5/3. This suggests that residual energy is a constant fraction of the total energy and that in the asymptotic regime of Alfvénic turbulence magnetic and kinetic spectra have the same scaling. The –1.5 slope for energy and the –2 slope for residual energy, which have been suggested earlier, are incompatible with our numerics.

79 ASTRONOMY AND ASTROPHYSICS↗