Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMHDSTS: a new explicit numerical scheme for simulations of partially ionised solar plasma
254
0
0.0
(
0
)
Added by
Pedro A. Gonz\\'alez-Morales
Publication date
2018
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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 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.
For time-dependent problems with high-contrast multiscale coefficients, the time step size for explicit methods is affected by the magnitude of the coefficient parameter. With a suitable construction of multiscale space, one can achieve a stable temporal splitting scheme where the time step size is independent of the contrast. Consider the parabolic equation with heterogeneous diffusion parameter, the flow rates vary significantly in different regions due to the high-contrast features of the diffusivity. In this work, we aim to introduce a multirate partially explicit splitting scheme to achieve efficient simulation with the desired accuracy. We first design multiscale subspaces to handle flow with different speed. For the fast flow, we obtain a low-dimensional subspace with respect to the high-diffusive component and adopt an implicit time discretization scheme. The other multiscale subspace will take care of the slow flow, and the corresponding degrees of freedom are treated explicitly. Then a multirate time stepping is introduced for the two parts. The stability of the multirate methods is analyzed for the partially explicit scheme. Moreover, we derive local error estimators corresponding to the two components of the solutions and provide an upper bound of the errors. An adaptive local temporal refinement framework is then proposed to achieve higher computational efficiency. Several numerical tests are presented to demonstrate the performance of the proposed method.
We derive self-consistent formalism for the description of multi-component partially ionized solar plasma, by means of the coupled equations for the charged and neutral components for an arbitrary number of chemical species, and the radiation field. All approximations and assumptions are carefully considered. Generalized Ohms law is derived for the single-fluid and two-fluid formalism. Our approach is analytical with some order-of-magnitude support calculations. After general equations are developed we particularize to some frequently considered cases as for the interaction of matter and radiation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
