Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDissipative instabilities in a partially ionised prominence plasma slab: II. The effect of compressibility
63
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This present study deals with the dissipative instability that appears in a compressible partially ionised plasma slab embedded in a uniform magnetic field, modelling the state of the plasma in solar prominences. In the partially ionised plasma, the dominant dissipative effect is the Cowling resistivity. The regions outside the slab (modelling the solar corona) are fully ionised, and the dominant mechanism of dissipation is viscosity. Analytical solutions to the extended magnetohydrodynamic (MHD) equations are found inside and outside of the slab and solutions are matched at the boundaries of the slab. The dispersion relation is derived and solutions are found analytically in the slender slab limit, while the conditions necessary for the appearance of the instability is investigated numerically for the entire parameter space. Our study is focussed on the effect of the compressibility on the generation and evolution of instabilities. We find that compressibility reduces the threshold of the equilibrium flow, where waves can be unstable, to a level that is comparable to the internal cusp speed, which is of the same order of flow speeds that are currently observed in solar prominences. Our study addresses only the slow waves, as these are the most likely perturbations to become unstable, however the time-scales of the instability are found to be rather large ranging from $10^5$-$10^7$ seconds. It is determined that the instability threshold is further influenced by the concentration of neutrals and the strength of the viscosity of the corona. Interestingly, these two latter aspects have opposite effects. Our numerical analysis shows that the interplay between the equilibrium flow, neutrals and dispersion can change considerably the nature of waves.
rate research
Read More
We investigate the nature of dissipative instability appearing in a prominence planar thread filled with partially ionised plasma in the incompressible limit. The importance of partial ionisation is investigated in terms of the ionisation factor and wavelength of waves propagating in the slab. To highlight the role of partial ionisation, we have constructed models describing various situations we can meet in solar prominence fine structure. Matching the solutions for the transversal component of the velocity and total pressure at the interfaces between the prominence slab and surrounding plasmas, we derived a dispersion relation whose imaginary part describes the evolution of the instability. Results are obtained in the limit of weak dissipation. We have investigated the appearance of instabilities in prominence dark plumes using single and two-fluid approximations. We show that dissipative instabilities appear for flow speeds that are less than the Kelvin-Helmholtz instability threshold. The onset of instability is determined by the equilibrium flow strength, the ionisation factor of the plasma, the wavelength of waves and the ion-neutral collisional rate. For a given wavelength and for ionisation degrees closer to a neutral gas, the propagating waves become unstable for a narrow band of flow speeds, meaning that neutrals have a stabilising effect. Our results show that the partially ionised plasma describing prominence dark plumes becomes unstable only in a two-fluid (charged particles-neutrals) model, that is for periods that are smaller than the ion-neutral collision time. The present study improves our understanding of stability of solar prominences and the role of partial ionisation in destabilising the plasma. We show the necessity of two-fluid approximation when discussing the nature of instabilities: waves in a single fluid approximation show a great deal of stability.
We investigate the nature of dissipative instability at the boundary (seen here as tangential discontinuity) between the viscous corona and the partially ionised prominence plasma in the incompressible limit. The importance of the partial ionisation is investigated in terms of the ionisation fraction. Matching the solutions for the transversal component of the velocity and total pressure at the interface between the prominence and coronal plasmas, we derive a dispersion relation whose imaginary part describes the evolution of the instability. Results are obtained in the limit of weak dissipation. Using simple analytical methods, we show that dissipative instabilities appear for flow speeds that are lower than the Kelvin-Helmholtz instability threshold. While viscosity tends to destabilise the plasma, the effect of partial ionisation (through the Cowling resistivity) will act towards stabilising the interface. For ionisation degrees closer to a neutral gas the interface will be unstable for larger values of equilibrium flow. The same principle is assumed when studying the appearance of instability at the interface between prominences and dark plumes. The unstable mode appearing in this case has a very small growth rate and dissipative instability cannot explain the appearance of flows in plumes. The present study improves our understanding of the complexity of dynamical processes at the interface of solar prominences and solar corona, and the role partial ionisation can have on the stability of the plasma. Our results clearly show that the problem of partial ionisation introduces new aspects of plasma stability with consequences on the evolution of solar prominences.
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.
Slow-mode shocks are important in understanding fast magnetic reconnection, jet formation and heating in the solar atmosphere, and other astrophysical systems. The atmospheric conditions in the solar chromosphere allow both ionised and neutral particles to exist and interact. Under such conditions, fine substructures exist within slow-mode shocks due to the decoupling and recoupling of the plasma and neutral species. We study numerically the fine substructure within slow-mode shocks in a partially ionised plasma, in particular, analysing the formation of an intermediate transition within the slow-mode shock. High-resolution 1D numerical simulations are performed using the (Punderline{I}P) code using a two-fluid approach. We discover that long-lived intermediate (Alfven) shocks can form within the slow-mode shock, where there is a shock transition from above to below the Alfven speed and a reversal of the magnetic field across the shock front. The collisional coupling provides frictional heating to the neutral fluid, resulting in a Sedov-Taylor-like expansion with overshoots in the neutral velocity and neutral density. The increase in density results in a decrease of the Alfven speed and with this the plasma inflow is accelerated to above the Alfven speed within the finite width of the shock leading to the intermediate transition. This process occurs for a wide range of physical parameters and an intermediate shock is present for all investigated values of plasma-$beta$, neutral fraction, and magnetic angle. As time advances the magnitude of the magnetic field reversal decreases since the neutral pressure cannot balance the Lorentz force. The intermediate shock is long-lived enough to be considered a physical structure, independent of the initial conditions.
The interaction of plasma with magnetic field in the partially ionised solar atmosphere is frequently modelled via a single-fluid approximation, which is valid for the case of a strongly coupled collisional media, such as solar photosphere and low chromosphere. Under the single-fluid formalism the main non-ideal effects are described by a series of extra terms in the generalised induction equation and in the energy conservation equation. These effects are: Ohmic diffusion, ambipolar diffusion, the Hall effect, and the Biermann battery effect. From the point of view of the numerical solution of the single-fluid equations, when ambipolar diffusion or Hall effects dominate can introduce severe restrictions on the integration time step and can compromise the stability of the numerical scheme. In this paper we introduce two numerical schemes to overcome those limitations. The first of them is known as Super Time-Stepping (STS) and it is designed to overcome the limitations imposed when the ambipolar diffusion term is dominant. The second scheme is called the Hall Diffusion Scheme (HDS) and it is used when the Hall term becomes dominant. These two numerical techniques can be used together by applying Strang operator splitting. This paper describes the implementation of the STS and HDS schemes in the single-fluid code Mancha3D. The validation for each of these schemes is provided by comparing the analytical solution with the numerical one for a suite of numerical tests.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
