Papers › Stability of two-fluid partially-ionised slow-mode shock fronts
Stability of two-fluid partially-ionised slow-mode shock fronts
Ben Snow, Andrew Hillier
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
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.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections