Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSlow shocks and conduction fronts from Petschek reconnection of skewed magnetic fields: two-fluid effects
182
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In models of fast magnetic reconnection, flux transfer occurs within a small portion of a current sheet triggering stored magnetic energy to be thermalized by shocks. When the initial current sheet separates magnetic fields which are not perfectly anti-parallel, i.e. they are skewed, magnetic energy is first converted to bulk kinetic energy and then thermalized in slow magnetosonic shocks. We show that the latter resemble parallel shocks or hydrodynamic shocks for all skew angles except those very near the anti-parallel limit. As for parallel shocks, the structures of reconnection-driven slow shocks are best studied using two-fluid equations in which ions and electrons have independent temperature. Time-dependent solutions of these equations can be used to predict and understand the shocks from reconnection of skewed magnetic fields. The results differ from those found using a single-fluid model such as magnetohydrodynamics. In the two-fluid model electrons are heated indirectly and thus carry a heat flux always well below the free-streaming limit. The viscous stress of the ions is, however, typically near the fluid-treatable limit. We find that for a wide range of skew angles and small plasma beta an electron conduction front extends ahead of the slow shock but remains within the outflow jet. In such cases conduction will play a more limited role in driving chromospheric evaporation than has been predicted based on single-fluid, anti-parallel models.
rate research
Read More
A magnetohydrodynamic (MHD) shock front can be unstable to the corrugation instability, which causes a perturbed shock front to become increasingly corrugated with time. An ideal MHD parallel shock (where the velocity and magnetic fields are aligned) is unconditionally unstable to the corrugation instability, whereas the ideal hydrodynamic (HD) counterpart is unconditionally stable. For a partially ionised medium (for example the solar chromosphere), both hydrodynamic and magnetohydrodynamic species coexist and the stability of the system has not been studied. In this paper, we perform numerical simulations of the corrugation instability in two-fluid partially-ionised shock fronts to investigate the stability conditions, and compare the results to HD and MHD simulations. Our simulations consist of an initially steady 2D parallel shock encountering a localised upstream density perturbation. In MHD, this perturbation results in an unstable shock front and the corrugation grows with time. We find that for the two-fluid simulation, the neutral species can act to stabilise the shock front. A parameter study is performed to analyse the conditions under which the shock front is stable and unstable. We find that for very weakly coupled or very strongly coupled partially-ionised system the shock front is unstable, as the system tends towards MHD. However, for a finite coupling, we find that the neutrals can stabilise the shock front, and produce new features including shock channels in the neutral species. We derive an equation that relates the stable wavelength range to the ion-neutral and neutral-ion coupling frequencies and the Mach number. Applying this relation to umbral flashes give an estimated range of stable wavelengths between 0.6 and 56 km.
We investigate physical scaling laws for magnetic energy dissipation in solar flares, in the framework of the Sweet-Parker model and the Petschek model. We find that the total dissipated magnetic energy $E_{diss}$ in a flare depends on the mean magnetic field component $B_f$ associated with the free energy $E_f$, the length scale $L$ of the magnetic area, the hydrostatic density scale height $lambda$ of the solar corona, the Alfven Mach number $M_A=v_1/v_A$ (the ratio of the inflow speed $v_1$ to the Alfvenic outflow speed $v_A$), and the flare duration $tau_f$, i.e., $E_{diss} = (1/4pi) B_f^2 L lambda v_A M_A tau_f$, where the Alfven speed depends on the nonpotential field strength $B_{np}$ and the mean electron density $n_e$ in the reconnection outflow. Using MDI/SDO and AIA/SDO observations and 3-D magnetic field solutions obtained with the vertical-current approximation nonlinear force-free field code (VCA-NLFFF) we measure all physical parameters necessary to test scaling laws, which represents a new method to measure Alfven Mach numbers $M_A$, the reconnection rate, and the total free energy dissipated in solar flares.
The role of slow-mode MHD shocks in magnetic reconnection is one of great importance for energy conversion and transport, but in many astrophysical plasmas the plasma is not fully ionised. In this paper, we investigate, using numerical simulations, the role of collisional coupling between a proton-electron charge-neutral fluid and a neutral hydrogen fluid for the 1D Riemann problem initiated in a constant pressure and density background state by a discontinuity in the magnetic field. This system, in the MHD limit, is characterised by two waves: a fast-mode rarefaction wave that drives a flow towards a slow-mode MHD shock. The system evolves through four stage: initiation, weak coupling, intermediate coupling and a quasi steady state. The initial stages are characterised by an over-pressured neutral region that expands with characteristics of a blast wave. In the later stages, the system tends towards a self-similar solution where the main drift velocity is concentrated in the thin region of the shock front. Due to the nature of the system, the neutral fluid is overpressured by the shock when compared to a purely hydrodynamic shock which results in the neutral fluid expanding to form the shock precursor. The thickness of the shockfront once it has formed proportional to the ionisation fraction to the power -1.2, which is a smaller exponent than would be naively expected from simple scaling arguments. One interesting result is that the shock front is a continuous transition of the physical variables for sub-sonic velocity upstream of the shock front (a c-shock) to a sharp jump in the physical variables followed by a relaxation to the downstream values for supersonic upstream velocity (a j-shock). The frictional heating that results from the velocity drift across the shock front can amount to approximately two per cent of the reference magnetic energy.
Aspects of solar flare dynamics, such as chromospheric evaporation and flare light-curves, have long been studied using one-dimensional models of plasma dynamics inside a static flare loop, subjected to some energy input. While extremely successful at explaining the observed characteristics of flares, all such models so far have specified energy input ad hoc, rather than deriving it self-consistently. There is broad consensus that flares are powered by magnetic energy released through reconnection. Recent work has generalized Petscheks basic reconnection scenario, topological change followed by field line retraction and shock heating, to permit its inclusion into a one-dimensional flare loop model. Here we compare the gas dynamics driven by retraction and shocking to those from more conventional static loop models energized by ad hoc source terms. We find significant differences during the first minute, when retraction leads to larger kinetic energies and produces higher densities at the loop top, while ad hoc heating tends to rarify the loop top. The loop-top density concentration is related to the slow magnetosonic shock, characteristic of Petscheks model, but persists beyond the retraction phase occurring in the outflow jet. This offers an explanation for observed loop-top sources of X-ray and EUV emission, with advantages over that provided by ad hoc heating scenarios. The cooling phases of the two models are, however, notably similar to one another, suggesting observations at that stage will yield little information on the nature of energy input.
Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system. We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system generates a jitter/synchrotron spectrum.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
