In a recent numerical study [Ng et al., Astrophys. J. {bf 747}, 109, 2012], with a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on
a series of runs in which random photospheric motions are imposed for hundreds to thousands of al time in order to obtain converged statistical values. The heating rate found in these simulations saturate to a level that is independent of the Lundquist number. This scaling result was also supported by an analysis with the assumption of the Sweet-Parker scaling of the current sheets, as well as how the width, length and number of current sheets scale with Lundquist number. In order to test these assumptions, we have implemented an automated routine to analyze thousands of current sheets in these simulations and return statistical scalings for these quantities. It is found that the Sweet-Parker scaling is justified. However, some discrepancies are also found and require further study.
New developments in the theory and numerical simulation of a recently proposed one-dimensional nonlinear time-dependent fluid model [K. Avinash, A. Bhattacharjee, and S. Hu, Phys. Rev. Lett. 90, 075001 (2003)] for void formation in dusty plasmas are
presented. The model describes an initial instability caused by the ion drag, rapid nonlinear growth, and a nonlinear saturation mechanism that realizes a quasi-steady state containing a void. The earlier one-dimensional model has been extended to two and three dimensions (the latter, assuming spherical symmetry), using a more complete set of dynamical equations than was used in the earlier one-dimensional formulation. The present set of equations includes an ion continuity equation and a nonlinear ion drag operator. Qualitative features of void formation are shown to be robust with respect to different functional forms of the ion drag operator.
An analytical and numerical treatment is given of a constrained version of the tectonics model developed by Priest, Heyvaerts, & Title [2002]. We begin with an initial uniform magnetic field ${bf B} = B_0 hat{bf z}$ that is line-tied at the surfaces
$z = 0$ and $z = L$. This initial configuration is twisted by photospheric footpoint motion that is assumed to depend on only one coordinate ($x$) transverse to the initial magnetic field. The geometric constraints imposed by our assumption precludes the occurrence of reconnection and secondary instabilities, but enables us to follow for long times the dissipation of energy due to the effects of resistivity and viscosity. In this limit, we demonstrate that when the coherence time of random photospheric footpoint motion is much smaller by several orders of magnitude compared with the resistive diffusion time, the heating due to Ohmic and viscous dissipation becomes independent of the resistivity of the plasma. Furthermore, we obtain scaling relations that suggest that even if reconnection and/or secondary instabilities were to limit the build-up of magnetic energy in such a model, the overall heating rate will still be independent of the resistivity.
Whether the phenomenology governing MHD turbulence is Kolmogorov or Iroshnikov-Kraichnan (IK) remains an open question, theoretically as well as observationally. The ion heating profile observed in the solar wind provides a quantitative, if indirect,
observational constraint on the relevant phenomenology. Recently, a solar wind heating model based on Kolmogorov spectral scaling has produced reasonably good agreement with observations, provided the effect of turbulence generation due to pickup ions is included in the model. Without including the pickup ion contributions, the Kolmogorov scaling predicts a proton temperature profile that decays too rapidly beyond a radial distance of 15 AU. In the present study, we alter the heating model by applying an energy cascade rate based on IK scaling, and show that the model yields higher proton temperatures, within the range of observations, with or without the inclusion of the effect due to pickup ions. Furthermore, the turbulence correlation length based on IK scaling seems to follow the trend of observations better.
In recent years, a phenomenological solar wind heating model based on a turbulent energy cascade prescribed by the Kolmogorov theory has produced reasonably good agreement with observations on proton temperatures out to distances around 70 AU, provid
ed the effect of turbulence generation due to pickup ions is included in the model. In a recent study [Ng et al., J. Geophys. Res., 115, A02101 (2010)], we have incorporated in the heating model the energy cascade rate based on Iroshnikov-Kraichnan (IK) scaling. We showed that the IK cascade rate can also produce good agreement with observations, with or without the inclusion of pickup ions. This effect was confirmed both by integrating the model using average boundary conditions at 1 AU, and by applying a method [Smith et al., Astrophys. J., 638, 508 (2006)] that uses directly observed values as boundary conditions. The effects due to pickup ions is found to be less important for the IK spectrum, which is shallower than the Kolmogorov spectrum. In this paper, we will present calculations of the pickup ions effect in more details, and discuss the physical reason why a shallower spectrum generates less waves and turbulence.
A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability
(MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.
