No Arabic abstract
Observations of coronal jets increasingly suggest that local fragmentation and intermittency play an important role in the dynamics of these events. In this work we investigate this fragmentation in high-resolution simulations of jets in the closed-field corona. We study two realizations of the embedded-bipole model, whereby impulsive helical outflows are driven by reconnection between twisted and untwisted field across the domed fan plane of a magnetic null. We find that the reconnection region fragments following the onset of a tearing-like instability, producing multiple magnetic null points and flux-rope structures within the current layer. The flux ropes formed within the weak-field region in the center of the current layer are associated with blobs of density enhancement that become filamentary threads as the flux ropes are ejected from the layer, whereupon new flux ropes form behind them. This repeated formation and ejection of flux ropes provides a natural explanation for the intermittent outflows, bright blobs of emission, and filamentary structure observed in some jets. Additional observational signatures of this process are discussed. Essentially all jet models invoke reconnection between regions of locally closed and locally open field as the jet-generation mechanism. Therefore, we suggest that this repeated tearing process should occur at the separatrix surface between the two flux systems in all jets. A schematic picture of tearing-mediated jet reconnection in three dimensions is outlined.
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.
We numerically investigate the excitation and temporal evolution of oscillations in a two-dimensional coronal arcade by including the three-dimensional propagation of perturbations. The time evolution of impulsively generated perturbations is studied by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-beta approximation. As we neglect gas pressure the slow mode is absent and therefore only coupled MHD fast and Alfven modes remain. Two types of numerical experiments are performed. First, the resonant wave energy transfer between a fast normal mode of the system and local Alfven waves is analyzed. It is seen how, because of resonant coupling, the fast wave with global character transfers its energy to Alfvenic oscillations localized around a particular magnetic surface within the arcade, thus producing the damping of the initial fast MHD mode. Second, the time evolution of a localized impulsive excitation, trying to mimic a nearby coronal disturbance, is considered. In this case, the generated fast wavefront leaves its energy on several magnetic surfaces within the arcade. The system is therefore able to trap energy in the form of Alfvenic oscillations, even in the absence of a density enhancement such as that of a coronal loop. These local oscillations are subsequently phase-mixed to smaller spatial scales. The amount of wave energy trapped by the system via wave energy conversion strongly depends on the wavelength of perturbations in the perpendicular direction, but is almost independent from the ratio of the magnetic to density scale heights.
We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth which are independent of $S$ provided this is high enough, as it should be natural having the ideal case as a limit for $Stoinfty$. Our simulations confirm that the tearing instability growth rate can be as fast as $gammaapprox 0.6,{tau_A}^{-1}$, where $tau_A$ is the ideal Alfvenic time set by the macroscopic scales, for our least diffusive case with $S=10^7$. The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of $S$ explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with $S$, here calculated on the emph{local}, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfvenic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.
Coronal jets represent important manifestations of ubiquitous solar transients, which may be the source of significant mass and energy input to the upper solar atmosphere and the solar wind. While the energy involved in a jet-like event is smaller than that of nominal solar flares and coronal mass ejections (CMEs), jets share many common properties with these phenomena, in particular, the explosive magnetically driven dynamics. Studies of jets could, therefore, provide critical insight for understanding the larger, more complex drivers of the solar activity. On the other side of the size-spectrum, the study of jets could also supply important clues on the physics of transients close or at the limit of the current spatial resolution such as spicules. Furthermore, jet phenomena may hint to basic process for heating the corona and accelerating the solar wind; consequently their study gives us the opportunity to attack a broad range of solar-heliospheric problems.
We present the first three-dimensional (3D), hydrodynamic simulations of the core convection zone (CZ) and extended radiative zone spanning from 1% to 90% of the stellar radius of an intermediate mass (3 $mathrm{M}_odot$) star. This allows us to self-consistently follow the generation of internal gravity waves (IGWs) at the convective boundary and their propagation to the surface. We find that convection in the core is dominated by plumes. The frequency spectrum in the CZ and that of IGW generation is a double power law as seen in previous two-dimensional (2D) simulations. The spectrum is significantly flatter than theoretical predictions using excitation through Reynolds stresses induced by convective eddies alone. It is compatible with excitation through plume penetration. An empirically determined distribution of plume frequencies generally matches the one necessary to explain a large part of the observed spectrum. We observe waves propagating in the radiation zone and excited standing modes, which can be identified as gravity and fundamental modes. They show similar frequencies and node patterns to those predicted by the stellar oscillation code GYRE. The continuous part of the spectrum fulfills the IGW dispersion relation. A spectrum of tangential velocity and temperature fluctuations close to the surface is extracted, which are directly related to observable brightness variations in stars. Unlike 2D simulations we do not see the high frequencies associated with wave breaking, likely because these 3D simulations are more heavily damped.