Energy dynamics in linear MHD with ion parallel viscosity

Analytic results for the time dependences of the kinetic and magnetic energies of an incompressible magneto fluid threaded by a strong uniform magnetic field B₀ are obtained. The governing equations are the linearised magnetohydrodynamic (MHD) ones, but with the conventional Laplacian dissipation replaced by ion parallel viscous effects. The behaviour is shown to depend on the relative sizes of the Alfvén frequency and the viscous and resistive dissipation rates. For many cases equipartition of the kinetic and magnetic energy holds at the (Fourier) modal level. An important exception to this behaviour occurs for two-dimensional fluctuations, that is when the wavevectors are perpendicular to B₀

Kinetic helicity and MHD turbulence

The issue of dynamical anisotropy in helical three-dimensional magnetohydrodynamic turbulence with a mean magnetic field B₀ is investigated. Using high-resolution direct numerical simulations, we follow the evolution of various isotropic initial states characterized by their different values of the kinetic helicity. The cross helicity and magnetic helicity of the initial conditions are also varied. In agreement with earlier work, we find that such initial states become anisotropic in of order an eddy-turnover time, with correlation lengths parallel to B₀ remaining largely unchanged while finer scales are excited in the perpendicular directions. Moreover, it is found that the development of both the anisotropy and the energy are essentially independent of the initial level of kinetic helicity. The physics associated with this latter feature is discussed

Nonlocality and the critical Reynolds numbers of the minimum state magnetohydrodynamic turbulence

Magnetohydrodynamic (MHD) systems can be strongly nonlinear (turbulent) when their kinetic and magnetic Reynolds numbers are high, as is the case in many astrophysical and space plasma flows. Unfortunately these high Reynolds numbers are typically much greater than those currently attainable in numerical simulations of MHD turbulence. A natural question to ask is how can researchers be sure that their simulations have reproduced all of the most influential physics of the flows and magnetic fields? In this paper, a metric is defined to indicate whether the necessary physics of interest has been captured. It is found that current computing resources will typically not be sufficient to achieve this minimum state metric

Parallel and perpendicular cascades in solar wind turbulence

MHD-scale fluctuations in the velocity, magnetic, and density fields of the solar wind are routinely observed. The evolution of these fluctuations, as they are transported radially outwards by the solar wind, is believed to involve both wave and turbulence processes. The presence of an average magnetic field has important implications for the anisotropy of the fluctuations and the nature of the turbulent wavenumber cascades in the directions parallel and perpendicular to this field. In particular, if the ratio of the rms magnetic fluctuation strength to the mean field is small, then the parallel wavenumber cascade is expected to be weak and there are difficulties in obtaining a cascade in frequency. The latter has been invoked in order to explain the heating of solar wind fluctuations (above adiabatic levels) via energy transfer to scales where ion-cyclotron damping can occur. Following a brief review of classical hydrodynamic and magnetohydrodynamic (MHD) cascade theories, we discuss the distinct nature of parallel and perpendicular cascades and their roles in the evolution of solar wind fluctuations

Properties of mass-loading shocks: 1. Hydrodynamic considerations

The one-dimensional hydrodynamics of flows subjected to mass loading are considered anew, with particular emphasis placed on determining the properties of mass-loading shocks. This work has been motivated by recent observations of the outbound Halley bow shock (Neubauer et al., 1990), which cannot be understood in terms of simple hydrodynamical or magnetohydrodynamical descriptions. By including mass injection at the shock, we have investigated the properties of the Rankine-Hugoniot conditions on the basis of a geometric formulation of the entropy condition. Such a condition, which is more powerful than the usual thermodynamical formulation, serves to determine those solutions to the Rankine-Hugoniot conditions which correspond to a physically realizable downstream state. On this basis a concise theoretical description of hydrodynamic mass-loading shocks is obtained. We show that mass-loading shocks have more in common with combustion shocks than with ordinary nonreacting gas dynamical shocks. It is shown that for decelerated solutions to the Rankine-Hugoniot conditions to exist, the upstream flow speed u0 must satisfy u0 > ucrit > cs, where cs is the sound speed. Besides the usual supersonic-subsonic transition, mass-loading fronts can also admit a decelerating supersonic-supersonic transition, the structure of which consists of a sharp decrease in the flow velocity preceding a recovery and an increase in the final downstream flow speed. We suggest the possibility that such structures may describe the inbound Halley bow shock (Coates et al., 1987a). Both parallel and oblique shocks are considered, the primary difference being that oblique shocks are subjected to a shearing stress due to mass loading. It is conjectured that such a shearing may destabilize the shock

Propagation and dissipation of Alfvén waves in stellar atmospheres permeated by isothermal winds

We investigate the nonlinear evolution of Alfvén waves in a radially stratified isothermal atmosphere with wind, from the atmospheric base out to the Alfvénic point. Nonlinear interactions, triggered by wave reflection due to the atmospheric gradients, are assumed to occur mainly in directions perpendicular to the mean radial magnetic field. The nonlinear coupling between waves propagating in opposite directions is modeled by a phenomenological term, containing an integral turbulent length scale, which acts as a dissipative coefficient for waves of a given frequency. Although the wind acceleration profile is not determined self-consistently one may estimate the dissipation rate inside the layer and follow the evolution of an initial frequency spectrum. Reflection of low frequency waves drives dissipation across the whole spectrum, and steeper gradients, i.e. lower coronal temperatures, enhance the dissipation rate. Moreover, when reasonable wave amplitudes are considered, waves of all frequencies damp at the same rate and the spectrum is not modified substantially during propagation. Therefore the sub-Alfvénic coronal layer acts differently when waves interact nonlinearly, no longer behaving as a frequency dependent filter once reflection-generated nonlinear interactions are included, at least within the classes of models discussed here

Current singularities at finitely compressible three-dimensional magnetic null points

The formation of current singularities at line-tied two- and three-dimensional (2D and 3D, respectively) magnetic null points in a nonresistive magnetohydrodynamic environment is explored. It is shown that, despite the different separatrix structures of 2D and 3D null points, current singularities may be initiated in a formally equivalent manner. This is true no matter whether the collapse is triggered by flux imbalance within closed, line-tied null points or driven by externally imposed velocity fields in open, incompressible geometries. A Lagrangian numerical code is used to investigate the finite amplitude perturbations that lead to singular current sheets in collapsing 2D and 3D null points. The form of the singular current distribution is analyzed as a function of the spatial anisotropy of the null point, and the effects of finite gas pressure are quantified. It is pointed out that the pressure force, while never stopping the formation of the singularity, significantly alters the morphology of the current distribution as well as dramatically weakening its strength. The impact of these findings on 2D and 3D magnetic reconnection models is discussed

Reduced magnetohydrodynamics and parallel spectral transfer

The self-consistency of the reduced magnetohydrodynamics (RMHD) model is explored by examining whether (parallel) spectral transfer might invalidate the assumptions employed in deriving it. Using direct numerical simulations we find that transfer of energy to structures with high parallel wavenumber is in fact limited by ongoing perpendicular transfer. Thus, the dynamics associated with RMHD models remains consistent with the underlying assumptions of RMHD. In particular, in well-resolved simulations it is neither necessary nor correct to introduce additional dissipation terms that (artificially) damp spectral transfer parallel to the mean magnetic field B0

A two-component phenomenology for homogeneous magnetohydrodynamic turbulence

A one-point closure model for energy decay in three-dimensional magnetohydrodynamic (MHD) turbulence is developed. The model allows for influence of a large-scale magnetic field that may be of strength sufficient to induce Alfvén wave propagation effects, and takes into account components of turbulence in which either the wave-like character is negligible or is dominant. This two-component model evolves energy and characteristic length scales, and may be useful as a simple description of homogeneous MHD turbulent decay. In concert with spatial transport models, it can form the basis for approximate treatment of low-frequency plasma turbulence in a variety of solar, space, and astrophysical contexts

Direct comparisons of compressible magnetohydrodynamics and reduced magnetohydrodynamics turbulence

Direct numerical simulations of low Mach number compressible three-dimensional magnetohydrodynamic (CMHD3D) turbulence in the presence of a strong mean magnetic field are compared with simulations of reduced magnetohydrodynamics (RMHD). Periodic boundary conditions in the three spatial coordinates are considered. Different sets of initial conditions are chosen to explore the applicability of RMHD and to study how close the solution remains to the full compressible MHD solution as both freely evolve in time. In a first set, the initial state is prepared to satisfy the conditions assumed in the derivation of RMHD, namely, a strong mean magnetic field and plane-polarized fluctuations, varying weakly along the mean magnetic field. In those circumstances, simulations show that RMHD and CMHD3D evolve almost indistinguishably from one another. When some of the conditions are relaxed the agreement worsens but RMHD remains fairly close to CMHD3D, especially when the mean magnetic field is large enough. Moreover, the well-known spectral anisotropy effect promotes the dynamical attainment of the conditions for RMHD applicability. Global quantities (mean energies, mean-square current, and vorticity) and energy spectra from the two solutions are compared and point-to-point separation estimations are computed. The specific results shown here give support to the use of RMHD as a valid approximation of compressible MHD with a mean magnetic field under certain but quite practical conditions