-
Mixing, heating and ion-neutral decoupling induced by Rayleigh-Taylor instability in prominence-corona transition regions
Authors:
V. S. Lukin,
E. Khomenko,
B. Popescu Braileanu
Abstract:
This study explores non-linear development of the magnetized Rayleigh-Taylor instability (RTI) in a prominence-corona transition region. Using a two-fluid model of a partially ionized plasma, we compare RTI simulations for several different magnetic field configurations. We follow prior descriptions of the numerical prominence model [Popescu Braileanu et al., 2021a,b, 2023] and explore the charged…
▽ More
This study explores non-linear development of the magnetized Rayleigh-Taylor instability (RTI) in a prominence-corona transition region. Using a two-fluid model of a partially ionized plasma, we compare RTI simulations for several different magnetic field configurations. We follow prior descriptions of the numerical prominence model [Popescu Braileanu et al., 2021a,b, 2023] and explore the charged-neutral fluid coupling and plasma heating in a two-dimensional mixing layer for different magnetic field configurations. We also investigate how the shear in magnetic field surrounding a prominence may impact the release of the gravitational potential energy of the prominence material. We show that the flow decoupling is strongest in the plane normal to the direction of the magnetic field, where neutral pressure gradients drive ion-neutral drifts and frictional heating is balanced by adiabatic cooling of the expanding prominence material. We also show that magnetic field within the mixing plane can lead to faster non-linear release of the gravitational energy driving the RTI, while more efficiently heating the plasma via viscous dissipation of associated plasma flows. We relate the computational results to potential observables by highlighting how integrating over under-resolved two-fluid sub-structure may lead to misinterpretation of observational data.
△ Less
Submitted 27 February, 2024; v1 submitted 6 January, 2024;
originally announced January 2024.
-
Magnetic field amplification and structure formation by the Rayleigh-Taylor instability
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko
Abstract:
We report on results of high resolution two fluid non-linear simulations of the Rayleigh Taylor Instability (RTI) at the interface between a solar prominence and the corona. These follow results reported earlier by Popescu Braileanu et al. (2021a,b) on linear and early non-linear RTI dynamics in this environment. The simulations use a two fluid model that includes collisions between neutrals and c…
▽ More
We report on results of high resolution two fluid non-linear simulations of the Rayleigh Taylor Instability (RTI) at the interface between a solar prominence and the corona. These follow results reported earlier by Popescu Braileanu et al. (2021a,b) on linear and early non-linear RTI dynamics in this environment. The simulations use a two fluid model that includes collisions between neutrals and charges, including ionization/recombination, energy and momentum transfer, and frictional heating. High resolution 2.5D magnetized RTI simulations with the magnetic field dominantly normal to and slightly sheared with respect to the prominence plane demonstrate that in a fully developed state of RTI a large fraction of the gravitational energy of a prominence thread can be converted into quasi-turbulent energy of the magnetic field. RTI magnetic energy generation is further accompanied by magnetic and plasma density structure formation, including dynamic formation, break-up, and merging of current sheets and plasmoid sub-structures. The simulations show the role of flow decoupling and ionization/recombination reactions between the neutrals and charges on the structure formation in magnetized RTI. We provide a careful examination of sources and form of numerical dissipation of the evolving magnetic field structures.
△ Less
Submitted 22 December, 2022; v1 submitted 24 December, 2021;
originally announced December 2021.
-
Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread II. Effects of collisionality
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid m…
▽ More
In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account neutral viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper II, the sensitivity of the RTI dynamics to collisional effects for different magnetic field configurations supporting the prominence thread is explored. This is done by artificially varying, or eliminating, effects of both elastic and inelastic collisions by modifying the model equations. We find that ionization and recombination reactions between ionized and neutral fluids, if in equilibrium prior to the onset of the instability, do not substantially impact the development of the primary RTI. However, such reactions can impact development of secondary structures during mixing of the cold prominence and hotter surrounding coronal material. We find that collisionality within and between ionized and neutral particle populations play an important role in both linear and non-linear development of RTI, with ion-neutral collision frequency as the primary determining factor in development or damping of small scale structures. We also observe that degree and signatures of flow decoupling between ion and neutral fluids can depend both on the inter-particle collisionality and the magnetic field configuration of the prominence thread.
△ Less
Submitted 8 April, 2021; v1 submitted 29 January, 2021;
originally announced January 2021.
-
Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread. I. Effects of prominence magnetization and mass loading
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) sim…
▽ More
Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) simulations of the Rayleigh-Taylor instability (RTI) at a smoothly changing interface between a solar prominence thread and the corona. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper I, the sensitivity of the RTI dynamics to the prominence equilibrium configuration, including the impact of the magnetic field strength and shear supporting the prominence thread, and the amount of prominence mass-loading is explored. We show that, at small scales, a realistically smooth prominence-corona interface leads to qualitatively different linear RTI evolution than that expected for a discontinuous interface, while magnetic field shear has the stabilizing effect of reducing the growth rate or eliminating the instability. In the non-linear phase, we observe that in the presence of field shear the development of the instability leads to formation of coherent and interacting 2.5D magnetic structures, which, in turn, can lead to substantial plasma flow across magnetic field lines and associated decoupling of the fluid velocities of charges and neutrals.
△ Less
Submitted 15 December, 2020; v1 submitted 31 July, 2020;
originally announced July 2020.
-
Two-fluid simulations of waves in the solar chromosphere II. Propagation and damping of fast magneto-acoustic waves and shocks
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with propertie…
▽ More
Waves and shocks traveling through the solar chromospheric plasma are influenced by its partial ionization and weak collisional coupling, and may become susceptible to multi-fluid effects, similar to interstellar shock waves. In this study, we consider fast magneto-acoustic shock wave formation and propagation in a stratified medium, that is permeated by a horizontal magnetic field, with properties similar to that of the solar chromosphere. The evolution of plasma and neutrals is modeled using a two-fluid code that evolves a set of coupled equations for two separate fluids. We observed that waves in neutrals and plasma, initially coupled at the upper photosphere, become uncoupled at higher heights in the chromosphere. This decoupling can be a consequence of either the characteristic spatial scale at the shock front, that becomes similar to the collisional scale, or the change in the relation between the wave frequency, ion cyclotron frequency, and the collisional frequency with height. The decoupling height is a sensitive function of the wave frequency, wave amplitude, and the magnetic field strength. We observed that decoupling causes damping of waves and an increase in the background temperature due to the frictional heating. The comparison between analytical and numerical results allows us to separate the role of the nonlinear effects from the linear ones on the decoupling and damping of waves.
△ Less
Submitted 14 August, 2019;
originally announced August 2019.
-
Two-fluid simulations of waves in the solar chromosphere I: numerical code verification
Authors:
B. Popescu Braileanu,
V. S. Lukin,
E. Khomenko,
A. de Vicente
Abstract:
Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the pu…
▽ More
Solar chromosphere consists of a partially ionized plasma, which makes modeling the solar chromosphere a particularly challenging numerical task. Here we numerically model chromospheric waves using a two-fluid approach with a newly developed numerical code. The code solves two-fluid equations of conservation of mass, momentum and energy, together with the induction equation, for the case of the purely hydrogen plasma with collisional coupling between the charged and neutral fluid components. The implementation of a semi-implicit algorithm allows us to overcome the numerical stability constraints due to the stiff collisional terms. We test the code against analytical solutions of acoustic and Alfvén wave propagation in uniform medium in several regimes of collisional coupling.The results of our simulations are consistent with the analytical estimates, and with other results described in the literature. In the limit of a large collisional frequency, the waves propagate with a common speed of a single fluid. In the other limit of a vanishingly small collisional frequency, the Alfvén waves propagate with an Alfvén speed of the charged fluid only, while the perturbation in neutral fluid is very small. The acoustic waves in these limits propagate with the sound speed corresponding to either the charges or the neutrals, while the perturbation in the other fluid component is very small. Otherwise, when the collision frequency is similar to the real part of the wave frequency, the interaction between charges and neutrals through momentum transfer collisions cause alterations of the waves frequencies and damping of the wave amplitudes.
△ Less
Submitted 9 May, 2019;
originally announced May 2019.
-
Onset of secondary instabilities and plasma heating during magnetic reconnection in strongly magnetized regions of the low solar atmosphere
Authors:
Lei Ni,
Vyacheslav S. Lukin
Abstract:
We numerically study magnetic reconnection on different spatial scales and at different heights in the weakly ionized plasma of the low solar atmosphere (around $300-800$~km above the solar surface) within a reactive 2.5 D multi-fluid plasma-neutral model. We consider a strongly magnetized plasma ($β\sim 6\% $) evolving from a force-free magnetic configuration and perturbed to initialize formation…
▽ More
We numerically study magnetic reconnection on different spatial scales and at different heights in the weakly ionized plasma of the low solar atmosphere (around $300-800$~km above the solar surface) within a reactive 2.5 D multi-fluid plasma-neutral model. We consider a strongly magnetized plasma ($β\sim 6\% $) evolving from a force-free magnetic configuration and perturbed to initialize formation of a reconnection current sheet. On large scales, the resulting current sheets are observed to undergo a secondary 'plasmoid' instability. A series of simulations at different scales demonstrate a cascading current sheet formation process that terminates for current sheets with width of ~2m and length of $\sim100$~m, corresponding to the critical current sheet aspect ratio of $\sim50$. We also observe that the plasmoid instability is the primary physical mechanism accelerating the magnetic reconnection in this plasma parameter regime. After plasmoid instabilities appear, the reconnection rate sharply increases to a value of $\sim$ 0.035, observed to be independent of the Lundquist number. These characteristics are very similar to magnetic reconnection in fully ionized plasmas. In this low $β$ guide field reconnection regime, both the recombination and collisionless effects are observed to have a small contribution to the reconnection rate. The simulations show that it is difficult to heat the dense weakly ionized photospheric plasmas to above $2\times10^4$~K during the magnetic reconnection process. However, the plasmas in the low solar chromosphere can be heated above $3\times10^4$~K with reconnection magnetic fields of $500$~G or stronger.
△ Less
Submitted 25 October, 2018; v1 submitted 23 October, 2018;
originally announced October 2018.
-
Magnetic reconnection in the low solar chromosphere with a more realistic radiative cooling model
Authors:
Lei Ni,
Vyacheslav S. Lukin,
Nicholas A. Murphy,
Jun Lin
Abstract:
Magnetic reconnection is the most likely mechanism responsible for the high temperature events that are observed in strongly magnetized locations around the temperature minimum in the low solar chromosphere. This work improves upon our previous work ["Magnetic Reconnection in Strongly Magnetized Regions of the Low Solar Chromosphere", The Astrophysical Journal 852, 95 (2018) ] by using a more real…
▽ More
Magnetic reconnection is the most likely mechanism responsible for the high temperature events that are observed in strongly magnetized locations around the temperature minimum in the low solar chromosphere. This work improves upon our previous work ["Magnetic Reconnection in Strongly Magnetized Regions of the Low Solar Chromosphere", The Astrophysical Journal 852, 95 (2018) ] by using a more realistic radiative cooling model computed from the OPACITY project and the CHIANTI database. We find that the rate of ionization of the neutral component of the plasma is still faster than recombination within the current sheet region. For low $β$ plasmas, the ionized and neutral fluid flows are well-coupled throughout the reconnection region resembling the single-fluid Sweet-Parker model dynamics. Decoupling of the ion and neutral inflows appears in the higher $β$ case with $β_0=1.46$, which leads to a reconnection rate about three times faster than the rate predicted by the Sweet-Parker model. In all cases, the plasma temperature increases with time inside the current sheet, and the maximum value is above $2\times10^4$~K when the reconnection magnetic field strength is greater than $500$~G. While the more realistic radiative cooling model does not result in qualitative changes of the characteristics of magnetic reconnection, it is necessary for studying the variations of the plasma temperature and ionization fraction inside current sheets in strongly magnetized regions of the low solar atmosphere. It is also important for studying energy conversion during the magnetic reconnection process when the hydrogen-dominated plasma approaches full ionization.
△ Less
Submitted 25 October, 2018; v1 submitted 16 April, 2018;
originally announced April 2018.
-
Reflection of fast magnetosonic waves near magnetic reconnection region
Authors:
Elena Provornikova,
J. Martin Laming,
Vyacheslav S. Lukin
Abstract:
Magnetic reconnection in the solar corona is thought to be unstable to the formation of multiple interacting plasmoids, and previous studies have shown that plasmoid dynamics can trigger MHD waves of different modes propagating outward from the reconnection site. However, variations in plasma parameters and magnetic field strength in the vicinity of a coronal reconnection site may lead to wave ref…
▽ More
Magnetic reconnection in the solar corona is thought to be unstable to the formation of multiple interacting plasmoids, and previous studies have shown that plasmoid dynamics can trigger MHD waves of different modes propagating outward from the reconnection site. However, variations in plasma parameters and magnetic field strength in the vicinity of a coronal reconnection site may lead to wave reflection and mode conversion. In this paper we investigate the reflection and refraction of fast magnetoacoustic waves near a reconnection site. Under a justified assumption of an analytically specified Alfvén speed profile, we derive and solve analytically the full wave equation governing propagation of fast mode waves in a non-uniform background plasma without recourse to the small-wavelength approximation. We show that the waves undergo reflection near the reconnection current sheet due to the Alfvén speed gradient and that the reflection efficiently depends on the plasma-$β$ parameter as well as on the wave frequency. In particular, we find that waves are reflected more efficiently near reconnection sites in a low-$β$ plasma which is typical for the solar coronal conditions. Also, the reflection is larger for lower frequency waves while high frequency waves propagate outward from the reconnection region almost without the reflection. We discuss the implications of efficient wave reflection near magnetic reconnection sites in strongly magnetized coronal plasma for particle acceleration, and also the effect this might have on First Ionization Potential (FIP) fractionation by the ponderomotive force of these waves in the chromosphere.
△ Less
Submitted 29 March, 2018;
originally announced March 2018.
-
Magnetic reconnection in strongly magnetized regions of the low solar chromosphere
Authors:
Lei Ni,
Vyacheslav S. Lukin,
Nicholas A. Murphy,
Jun Lin
Abstract:
Magnetic reconnection in strongly magnetized regions around the temperature minimum region of the low solar atmosphere is studied by employing MHD-based simulations of a partially ionized plasma within a reactive 2.5D multi-fluid model. It is shown that in the absence of magnetic nulls in a low $β$ plasma the ionized and neutral fluid flows are well-coupled throughout the reconnection region. Howe…
▽ More
Magnetic reconnection in strongly magnetized regions around the temperature minimum region of the low solar atmosphere is studied by employing MHD-based simulations of a partially ionized plasma within a reactive 2.5D multi-fluid model. It is shown that in the absence of magnetic nulls in a low $β$ plasma the ionized and neutral fluid flows are well-coupled throughout the reconnection region. However, non-equilibrium ionization-recombination dynamics play a critical role in determining the structure of the reconnection region, lead to much lower temperature increases and a faster magnetic reconnection rate as compared to simulations that assume plasma to be in ionization-recombination equilibrium. The rate of ionization of the neutral component of the plasma is always faster than recombination within the current sheet region even when the initial plasma $β$ is as high as $β_0=1.46$. When the reconnecting magnetic field is in excess of a kilogauss and the plasma $β$ is lower than 0.0145, the initially weakly ionized plasmas can become fully ionized within the reconnection region and the current sheet can be strongly heated to above $2.5\times10^4$~K, even as most of the collisionally dissipated magnetic energy is radiated away. The Hall effect increases the reconnection rate slightly, but in the absence of magnetic nulls it does not result in significant asymmetries or change the characteristics of the reconnection current sheet down to meter scales.
△ Less
Submitted 22 December, 2017; v1 submitted 2 December, 2017;
originally announced December 2017.
-
Plasma compression in magnetic reconnection regions in the solar corona
Authors:
Elena Provornikova,
John Martin Laming,
Vyacheslav S. Lukin
Abstract:
It has been proposed that particles bouncing between magnetized flows converging in a reconnection region can be accelerated by the first order Fermi mechanism. Analytical considerations of this mechanism have shown that the spectral index of accelerated particles is related to the total plasma compression within the reconnection region similarly to the case of diffusive shock acceleration mechani…
▽ More
It has been proposed that particles bouncing between magnetized flows converging in a reconnection region can be accelerated by the first order Fermi mechanism. Analytical considerations of this mechanism have shown that the spectral index of accelerated particles is related to the total plasma compression within the reconnection region similarly to the case of diffusive shock acceleration mechanism. As a first step to investigate the efficiency of Fermi acceleration in reconnection regions in producing hard energy spectra of particles in the solar corona, we explore the degree of plasma compression that can be achieved at reconnection sites. In particular, we aim to determine the conditions for the strong compressions to form. Using a two-dimensional resistive MHD numerical model we consider a set of magnetic field configurations where magnetic reconnection can occur including a Harris current sheet, a force-free current sheet, and two merging flux ropes. Plasma parameters are taken to be characteristic of the solar corona. Numerical simulations show that strong plasma compressions ($\geq 4$) in the reconnection regions can form when the plasma heating due to reconnection is efficiently removed by fast thermal conduction or radiative cooling process. The radiative cooling process which is negligible in the typical 1 MK corona can play an important role in the low corona/transition region. It is found that plasma compression is expected to be strongest in low-beta plasma $β\sim 0.01-0.07$ at reconnection magnetic nulls.
△ Less
Submitted 28 April, 2016; v1 submitted 25 April, 2016;
originally announced April 2016.
-
Asymmetric Magnetic Reconnection in Weakly Ionized Chromospheric Plasmas
Authors:
Nicholas A. Murphy,
Vyacheslav S. Lukin
Abstract:
Realistic models of magnetic reconnection in the solar chromosphere must take into account that the plasma is partially ionized and that plasma conditions within any two magnetic flux bundles undergoing reconnection may not be the same. Asymmetric reconnection in the chromosphere may occur when newly emerged flux interacts with pre-existing, overlying flux. We present 2.5D simulations of asymmetri…
▽ More
Realistic models of magnetic reconnection in the solar chromosphere must take into account that the plasma is partially ionized and that plasma conditions within any two magnetic flux bundles undergoing reconnection may not be the same. Asymmetric reconnection in the chromosphere may occur when newly emerged flux interacts with pre-existing, overlying flux. We present 2.5D simulations of asymmetric reconnection in weakly ionized, reacting plasmas where the magnetic field strengths, ion and neutral densities, and temperatures are different in each upstream region. The plasma and neutral components are evolved separately to allow non-equilibrium ionization. As in previous simulations of chromospheric reconnection, the current sheet thins to the scale of the neutral-ion mean free path and the ion and neutral outflows are strongly coupled. However, the ion and neutral inflows are asymmetrically decoupled. In cases with magnetic asymmetry, a net flow of neutrals through the current sheet from the weak field (high density) upstream region into the strong field upstream region results from a neutral pressure gradient. Consequently, neutrals dragged along with the outflow are more likely to originate from the weak field region. The Hall effect leads to the development of a characteristic quadrupole magnetic field modified by asymmetry, but the X-point geometry expected during Hall reconnection does not occur. All simulations show the development of plasmoids after an initial laminar phase.
△ Less
Submitted 6 April, 2015;
originally announced April 2015.
-
On Flux Rope Stability and Atmospheric Stratification in Models of Coronal Mass Ejections Triggered by Flux Emergence
Authors:
E. Lee,
V. S. Lukin,
M. G. Linton
Abstract:
Flux emergence is widely recognized to play an important role in the initiation of coronal mass ejections. The Chen-Shibata (2000) model, which addresses the connection between emerging flux and flux rope eruptions, can be implemented numerically to study how emerging flux through the photosphere can impact the eruption of a pre-existing coronal flux rope. The model's sensitivity to the initial co…
▽ More
Flux emergence is widely recognized to play an important role in the initiation of coronal mass ejections. The Chen-Shibata (2000) model, which addresses the connection between emerging flux and flux rope eruptions, can be implemented numerically to study how emerging flux through the photosphere can impact the eruption of a pre-existing coronal flux rope. The model's sensitivity to the initial conditions and reconnection micro-physics is investigated with a parameter study. In particular, we aim to understand the stability of the coronal flux rope in the context of X-point collapse and the effects of boundary driving in both unstratified and stratified atmospheres. In the absence of driving, we assess the behavior of waves in the vicinity of the X-point. With boundary driving applied, we study the effects of reconnection micro-physics and atmospheric stratification on the eruption. We find that the Chen-Shibata equilibrium can be unstable to an X-point collapse even in the absence of driving due to wave accumulation at the X-point. However, the equilibrium can be stabilized by reducing the compressibility of the plasma, which allows small-amplitude waves to pass through the X-point without accumulation. Simulations with the photospheric boundary driving evaluate the impact of reconnection micro-physics and atmospheric stratification on the resulting dynamics: we show the evolution of the system to be determined primarily by the structure of the global magnetic fields with little sensitivity to the micro-physics of magnetic reconnection; and in a stratified atmosphere, we identify a novel mechanism for producing quasi-periodic behavior at the reconnection site behind a rising flux rope as a possible explanation of similar phenomena observed in solar and stellar flares.
△ Less
Submitted 14 August, 2014; v1 submitted 2 March, 2014;
originally announced March 2014.
-
Ionized Plasma and Neutral Gas Coupling in the Sun's Chromosphere and Earth's Ionosphere/Thermosphere
Authors:
J. E. Leake,
C. R. DeVore,
J. P. Thayer,
A. G. Burns,
G. Crowley,
H. R. Gilbert,
J. D Huba,
J. Krall,
M. G. Linton,
V. S. Lukin,
W. Wang
Abstract:
We review our understanding of ionized plasma and neutral gas coupling in the weakly ionized, stratified, electromagnetically-permeated regions of the Sun's chromosphere and Earth's ionosphere/thermosphere. Using representative models for each environment we derive fundamental descriptions of the coupling of the constituent parts to each other and to the electric and magnetic fields, and we examin…
▽ More
We review our understanding of ionized plasma and neutral gas coupling in the weakly ionized, stratified, electromagnetically-permeated regions of the Sun's chromosphere and Earth's ionosphere/thermosphere. Using representative models for each environment we derive fundamental descriptions of the coupling of the constituent parts to each other and to the electric and magnetic fields, and we examine the variation in magnetization of the ionized component. Using these descriptions we compare related phenomena in the two environments, and discuss electric currents, energy transfer and dissipation. We present a coupled theoretical and numerical study of plasma instabilities in the two environments that serves as an example of how the chromospheric and ionospheric communities can further collaborate. We also suggest future collaborative studies that will help improve our understanding of these two different atmospheres which share many similarities, but have large disparities in key quantities.
△ Less
Submitted 29 August, 2016; v1 submitted 1 October, 2013;
originally announced October 2013.
-
Two-fluid simulations of driven reconnection in the Mega-Ampere Spherical Tokamak
Authors:
A. Stanier,
P. Browning,
M. Gordovskyy,
K. G. McClements,
M. P. Gryaznevich,
V. S. Lukin
Abstract:
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment.
In this paper, 2D fluid simulatio…
▽ More
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment.
In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities.
In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.
△ Less
Submitted 13 August, 2013;
originally announced August 2013.
-
Magnetic Reconnection in a Weakly Ionized Plasma
Authors:
James E Leake,
Vyacheslav S Lukin,
Mark G Linton
Abstract:
Magnetic reconnection in partially ionized plasmas is a ubiquitous phenomenon spanning the range from laboratory to intergalactic scales, yet it remains poorly understood and relatively little studied. Here, we present results from a self-consistent multi-fluid simulation of magnetic reconnection in a weakly ionized reacting plasma with a particular focus on the parameter regime of the solar chrom…
▽ More
Magnetic reconnection in partially ionized plasmas is a ubiquitous phenomenon spanning the range from laboratory to intergalactic scales, yet it remains poorly understood and relatively little studied. Here, we present results from a self-consistent multi-fluid simulation of magnetic reconnection in a weakly ionized reacting plasma with a particular focus on the parameter regime of the solar chromosphere. The numerical model includes collisional transport, interaction and reactions between the species, and optically thin radiative losses. This model improves upon our previous work in Leake et al. 2012 \cite{Leake2012} by considering realistic chromospheric transport coefficients, and by solving a generalized Ohm's law that accounts for finite ion-inertia and electron-neutral drag.
We find that during the two dimensional reconnection of a Harris current sheet with an initial width larger than the neutral-ion collisional coupling scale, the current sheet thins until its width becomes less than this coupling scale, and the neutral and ion fluids decouple upstream from the reconnection site. During this process of decoupling, we observe reconnection faster than the single-fluid Sweet-Parker prediction, with recombination and plasma outflow both playing a role in determining the reconnection rate.
As the current sheet thins further and elongates it becomes unstable to the secondary tearing instability, and plasmoids are seen. The reconnection rate, outflows and plasmoids observed in this simulation provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets.
△ Less
Submitted 13 February, 2013;
originally announced February 2013.
-
Multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma
Authors:
James E. Leake,
Vyacheslav S. Lukin,
Mark G. Linton,
Eric T. Meier
Abstract:
We present results from the first self-consistent multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma. We simulate two dimensional magnetic reconnection in a Harris current sheet with a numerical model which includes ion-neutral scattering collisions, ionization, recombination, optically thin radiative loss, collisional heating, and thermal conduction…
▽ More
We present results from the first self-consistent multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma. We simulate two dimensional magnetic reconnection in a Harris current sheet with a numerical model which includes ion-neutral scattering collisions, ionization, recombination, optically thin radiative loss, collisional heating, and thermal conduction. In the resulting tearing mode reconnection the neutral and ion fluids become decoupled upstream from the reconnection site, creating an excess of ions in the reconnection region and therefore an ionization imbalance. Ion recombination in the reconnection region, combined with Alfvénic outflows, quickly removes ions from the reconnection site, leading to a fast reconnection rate independent of Lundquist number. In addition to allowing fast reconnection, we find that these non-equilibria partial ionization effects lead to the onset of the nonlinear secondary tearing instability at lower values of the Lundquist number than has been found in fully ionized plasmas.These simulations provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets.
△ Less
Submitted 5 October, 2012;
originally announced October 2012.
