Source author record

Andrew Hillier

Andrew Hillier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2023arXiv

Connecting theory of plasmoid-modulated reconnection to observations of solar flares

The short timescale of the solar flare reconnection process has long proved to be a puzzle. Recent studies suggest the importance of the formation of plasmoids in the reconnecting current sheet, with quantifying the aspect ratio of the width to length of the current sheet in terms of a negative power $α$ of the Lundquist number, i.e. $S^{-α}$, being key to understanding the onset of plasmoids formation. In this paper we make the first application of theoretical scalings for this aspect ratio to observed flares to evaluate how plasmoid formation may connect with observations. We find that for three different flares showing plasmoids a range of $α$ values of $α= 0.27$ to $0.31$. The values in this small range implies that plasmoids may be forming before the theoretically predicted critical aspect ratio ($α=1/3$) has been reached, potentially presenting a challenge for the theoretical models.

preprint2020arXiv

Estimating the energy dissipation {from Kelvin-Helmholtz instability induced} turbulence in oscillating coronal loops}

Kelvin-Helmholtz {instability induced} turbulence is one promising mechanism by which loops in the solar corona can be heated by MHD waves. In this paper we present an analytical model of the dissipation rate of {Kelvin-Helmholtz instability induced} turbulence $\varepsilon_{\rm D}$, finding it scales as the wave amplitude ($d$) to the third power ($\varepsilon_{\rm D}\propto d^3$). Based on the concept of steady-state turbulence, we expect the turbulence heating throughout the volume of {the} loop to match the total energy injected through its footpoints. In situations where this holds, the wave amplitude has to vary as the cube-root of the injected energy. Comparing the analytic results with those of simulations shows that our analytic formulation captures the key aspects of the turbulent dissipation from the numerical work. Applying this model to the observed characteristics of decayless kink waves we predict that the amplitudes of these observed waves is insufficient to turbulently heat the solar corona.

preprint2020arXiv

Mode conversion of two-fluid shocks in a partially-ionised, isothermal, stratified atmosphere

The plasma of the lower solar atmosphere consists of mostly neutral particles, whereas the upper solar atmosphere is mostly ionised particles and electrons. A shock that propagates upwards in the solar atmosphere therefore undergoes a transition where the dominant fluid is either neutral or ionised. An upwards propagating shock also passes a point where the sound and Alfvén speed are equal. At this point the energy of the acoustic shock can separated into fast and slow components. How the energy is distributed between the two modes depends on the angle of magnetic field. The separation of neutral and ionised species in a gravitationally stratified atmosphere is investigated. The role of two-fluid effects on the structure of the shocks post-mode-conversion and the frictional heating is quantified for different levels of collisional coupling. Two-fluid numerical simulations are performed using the (P\underline{I}P) code of a wave steepening into a shock in an isothermal, partially-ionised atmosphere. The collisional coefficient is varied to investigate the regimes where the plasma and neutral species are weakly, strongly and finitely coupled. The propagation speeds of the compressional waves hosted by neutral and ionised species vary, therefore velocity drift between the two species is produced as the plasma attempts to propagate faster than the neutrals. This is most extreme for a fast-mode shock. We find that the collisional coefficient drastically changes the features present in the system, specifically the mode conversion height, type of shocks present, and the finite shock widths created by the two-fluid effects. In the finitely-coupled regime fast-mode shock widths can exceed the pressure scale height leading to a new potential observable of two-fluid effects in the lower solar atmosphere.

preprint2019arXiv

Coronal cooling as a result of mixing by the nonlinear Kelvin--Helmholtz instability

Recent observations show cool, oscillating prominence threads fading when observed in cool spectral lines and appearing in warm spectral lines. A proposed mechanism to explain this evolution is that the threads were heated by turbulence driven by the Kelvin--Helmholtz instability that developed as a result of wave-driven shear flows on the surface of the thread. As the Kelvin--Helmholtz instability is an instability that works to mix the fluids, in the solar corona it can be expected to work by mixing the cool prominence material with that of the hot corona to form a warm boundary layer. In this paper we develop a simple phenomenological model of nonlinear Kelvin--Helmholtz mixing, using it to determine the characteristic density and temperature of the mixing layer, which for the case under study with constant pressure across the two fluids are $ρ_{\rm mixed}=\sqrt{ρ_1ρ_2}$ and $T_{\rm mixed}=\sqrt{T_1T_2}$. One result from the model is that it provides an accurate, as determined by comparison with simulation results, determination of the kinetic energy in the mean velocity field. A consequence of this is that the magnitude of turbulence, and with it the energy that can be dissipated on fast time-scales, as driven by this instability can be determined. For the prominence-corona system, the mean temperature rise possible from turbulent heating is estimated to be less than 1\% of the characteristic temperature (which is found to be $10^5$\,K). These results highlight that mixing, and not heating, are likely to be the cause of the observed transition between cool to warm material in Okamoto et. al (2015). One consequence of this result is that the mixing creates a region with higher radiative loss rates on average than either of the original fluids, meaning that this instability could contribute a net loss of thermal energy from the corona, i.e. coronal cooling.

preprint2019arXiv

Ion-neutral decoupling in the nonlinear Kelvin--Helmholtz instability: Case of field-aligned flow

The nonlinear magnetic Kelvin-Helmholtz instability (KHi), and the turbulence it creates, appears in many astrophysical systems. This includes those systems where the local plasma conditions are such that the plasma is not fully ionised, for example in the lower solar atmosphere and molecular clouds. In a partially ionised system, the fluids couple via collisions which occur at characteristic frequencies, therefore neutral and plasma species become decoupled for sufficiently high-frequency dynamics. Here we present high-resolution 2D two-fluid simulations of the nonlinear KHi for a system that traverses the dynamic scales between decoupled fluids and coupled dynamics. We discover some interesting phenomena, including the presence of a density coupling that is independent of the velocity coupling. Using these simulations we analyse the heating rate, and two regimes appear. The first is a regime where the neutral flow is decoupled from the magnetic field that is characterised with a constant heating rate, then at larger scales the strong coupling approximation holds and the heating rate. At large scales with the KHi layer width to the $-2$ power. There is an energy cascade in the simulation, but the nature of the frictional heating means the heating rate is determined by the largest scale of the turbulent motions, a fact that has consequences for understanding turbulent dissipation in multi-fluid systems.

preprint2016arXiv

Nonlinear Instability and Intermittent Nature of Magnetic Reconnection in Solar Chromosphere

The recent observations of Singh et al. (2012) have shown multiple plasma ejections and the intermittent nature of magnetic reconnection in the solar chromosphere, highlighting the need for fast reconnection to occur in highly collisional plasma. However, the physical process through which fast magnetic reconnection occurs in partially ionized plasma, like the solar chromosphere, is still poorly understood. It has been shown that for sufficiently high magnetic Reynolds numbers, Sweet-Parker current sheets can become unstable leading to tearing mode instability and plasmoid formation, but when dealing with a partially ionized plasma the strength of coupling between the ions and neutrals plays a fundamental role in determining the dynamics of the system. We propose that as the reconnecting current sheet thins and the tearing instability develops, plasmoid formation passes through strongly, intermediately, and weakly coupled (or decoupled) regimes, with the time scale for the tearing mode instability depending on the frictional coupling between ions and neutrals. We present calculations for the relevant time scales for fractal tearing in all three regimes. We show that as a result of the tearing mode instability and the subsequent non-linear instability due to the plasmoid-dominated reconnection, the Sweet-Parker current sheet tends to have a fractal-like structure, and when the chromospheric magnetic field is sufficiently strong the tearing instability can reach down to kinetic scales, which are hypothesized to be necessary for fast reconnection.

preprint2016arXiv

On the nature of the magnetic Rayleigh-Taylor instability in Astrophysical Plasma: The case of uniform magnetic field strength

The magnetic Rayleigh-Taylor instability has been shown to play a key role in many astrophysical systems. The equation for the growth rate of this instability in the incompressible limit, and the most-unstable mode that can be derived from it, are often used to estimate the strength of the magnetic field that is associated with the observed dynamics. However, there are some issues with the interpretations given. Here we show that the class of most unstable modes $k_u$ for a given $θ$, the class of modes often used to estimate the strength of the magnetic field from observations, for the system leads to the instability growing as $σ^2=1/2 A g k_u$, a growth rate which is independent of the strength of the magnetic field and which highlights that small scales are preferred by the system, but not does not give the fastest growing mode for that given $k$. We also highlight that outside of the interchange ($\mathbf{k}\cdot\mathbf{B}=0$) and undular ($\mathbf{k}$ parallel to $\mathbf{B}$) modes, all the other modes have a perturbation pair of the same wavenumber and growth rate that when excited in the linear regime can result in an interference pattern that gives field aligned filamentary structure often seen in 3D simulations. The analysis was extended to a sheared magnetic field, where it was found that it was possible to extend the results for a non-sheared field to this case. We suggest that without magnetic shear it is too simplistic to be used to infer magnetic field strengths in astrophysical systems.

preprint2016arXiv

The formation and evolution of reconnection-driven slow-mode shocks in a partially ionised plasma

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.

preprint2014arXiv

Investigating the Dynamics and Density Evolution of Returning Plasma Blobs from the 2011 June 7 Eruption

This work examines infalling matter following an enormous Coronal Mass Ejection (CME) on 2011 June 7. The material formed discrete concentrations, or blobs, in the corona and fell back to the surface, appearing as dark clouds against the bright corona. In this work we examined the density and dynamic evolution of these blobs in order to formally assess the intriguing morphology displayed throughout their descent. The blobs were studied in five wavelengths (94, 131, 171, 193 and 211 Å) using the Solar Dynamics Observatory Atmospheric Imaging Assembly (SDO/AIA), comparing background emission to attenuated emission as a function of wavelength to calculate column densities across the descent of four separate blobs. We found the material to have a column density of hydrogen of approximately 2 $\times$ 10$^{19}$ cm$^{-2}$, which is comparable with typical pre-eruption filament column densities. Repeated splitting of the returning material is seen in a manner consistent with the Rayleigh-Taylor instability. Furthermore, the observed distribution of density and its evolution are also a signature of this instability. By approximating the three-dimensional geometry (with data from STEREO-A), volumetric densities were found to be approximately 2 $\times$ 10$^{-14}$ g cm$^{-3}$, and this, along with observed dominant length-scales of the instability, was used to infer a magnetic field of the order 1 G associated with the descending blobs.

preprint2014arXiv

Superflare occurrence and energies on G-, K- and M-type dwarfs

Kepler data from G-, K- and M-type stars are used to study conditions that lead to superflares with energies above $10^{34} {\rm erg}$. From the 117,661 stars included, 380 show superflares with a total of 1690 such events. We study whether parameters, like effective temperature or the rotation rate, have any effect on the superflare occurrence rate or energy. With increasing effective temperature we observe a decrease in the superflare rate, which is analogous to the previous findings of a decrease in dynamo activity with increasing effective temperature. For slowly rotating stars, we find a quadratic increase of the mean occurrence rate with the rotation rate up to a critical point, after which the rate decreases linearly. Motivated by standard dynamo theory, we study the behavior of the relative starspot coverage, approximated as the relative brightness variation. For faster rotating stars, an increased fraction of stars shows higher spot coverage, which leads to higher superflare rates. A turbulent dynamo is used to study the dependence of the Ohmic dissipation as a proxy of the flare energy on the differential rotation or shear rate. The resulting statistics of the dissipation energy as a function of dynamo number is similar to the observed flare statistics as a function of the inverse Rossby number and shows similarly strong fluctuations. This supports the idea that superflares might well be possible for solar-type G stars.

preprint2013arXiv

On the Support of Solar Prominence Material by the Dips of a Coronal Flux Tube

The dense prominence material is believed to be supported against gravity through the magnetic tension of dipped coronal magnetic field. For quiescent prominences, which exhibit many gravity-driven flows, hydrodynamic forces are likely to play an important role in the determination of both the large and small scale magnetic field distributions. In this study, we present the first steps toward creating three-dimensional magneto-hydrostatic prominence model where the prominence is formed in the dips of a coronal flux tube. Here 2.5D equilibria are created by adding mass to an initially force-free magnetic field, then performing a secondary magnetohydrodynamic relaxation. Two inverse polarity magnetic field configurations are studied in detail, a simple o-point configuration with a ratio of the horizontal field (B_x) to the axial field (B_y) of 1:2 and a more complex model that also has an x-point with a ratio of 1:11. The models show that support against gravity is either by total pressure or tension, with only tension support resembling observed quiescent prominences. The o-point of the coronal flux tube was pulled down by the prominence material, leading to compression of the magnetic field at the base of the prominence. Therefore tension support comes from the small curvature of the compressed magnetic field at the bottom and the larger curvature of the stretched magnetic field at the top of the prominence. It was found that this method does not guarantee convergence to a prominence-like equilibrium in the case where an x-point exists below the prominence flux tube. The results imply that a plasma beta of ~ 0.1 is necessary to support prominences through magnetic tension.

preprint2012arXiv

Determination of Prominence Plasma Beta from the Dynamics of Rising Plumes

The launch of Hinode satellite led to the discovery of rising plumes, dark in chromospheric lines, in quiescent prominences that propagate from large (~10 Mm) bubbles that form at the base of the prominences. These plumes present a very interesting opportunity to study Magnetohydrodynamic (MHD) phenomena in quiescent prominences, but obstacles still remain. One of the biggest issues is that of the magnetic field strength, which is not easily measurable in prominences. In this paper we present a method that may be used to determine a prominence's plasma βwhen rising plumes are observed. Using the classic fluid dynamic solution for flow around a circular cylinder with an MHD correction, the compression of the prominence material can be estimated. This has been successfully confirmed through simulations; application to a prominence gave an estimate of the plasma beta as β=0.47 to 1.13 with an error of 0.080 for the range γ=1.4 to 1.7. Using this method it may be possible to estimate the plasma beta of observed prominences, therefore helping our understanding of a prominence's dynamics in terms of MHD phenomena.

preprint2011arXiv

Numerical Simulations of the Magnetic Rayleigh-Taylor Instability in the Kippenhahn-Schlüter Prominence Model

The launch of the Hinode satellite has allowed unprecedented high resolution, stable images of solar quiescent prominences to be taken over extended periods of time. These new images led to the discovery of dark upflows that propagated from the base of prominences developing highly turbulent profiles. As yet, how these flows are driven is not fully understood. To study the physics behind this phenomena we use 3-D magnetohydrodynamic (MHD) simulations to investigate the nonlinear stability of the Kippenhahn-Shlüter prominence model to the magnetic Rayleigh-Taylor instability. The model simulates the rise of a buoyant tube inside a quiescent prominence, where the upper boundary between the tube and prominence model is perturbed to excite the interchange of magnetic field lines. We found upflows of constant velocity (maximum found 6\,km\,s$^{-1}$) and a maximum plume width $\approx 1500$\,km which propagate through a height of approximately 6\,Mm in the no guide field case. The case with the strong guide field (initially $B_y=2B_x$) results in a large plume that rises though the prominence model at $\sim 5$\,km\,s$^{-1}$ with width $\sim 900$\,km (resulting in width of 2400\,km when viewed along the axis of the prominence) reaching a height of $\sim 3.1$\,Mm. In both cases nonlinear processes were important for determining plume dynamics.

preprint2011arXiv

Observations of Plasma Blob Ejection from a Quiescent Prominence by Hinode SOT

We report findings from 0.2" resolution observations of the 2007 October 03 quiescent prominence observed with the Solar Optical Telescope on the Hinode satellite. The observations show clear ejections from the top of the quiescent prominence of plasma blobs. The ejections, originating from the top of rising prominence threads, are impulsively accelerated to Alfvenic velocities and then undergo ballistic motion. The ejections have a characteristic size between ~ 1000 - 2000 km. These characteristics are similar to downwardly propagating knots (typical size ~ 700 km) that have been observed in prominence threads, we suggest that the plasma blob ejections could be the upward moving counterpart to the downwardly propagating knots. We discuss the tearing instability as a possible mechanism to explain the ejections.

preprint2010arXiv

Spicule Dynamics over Plage Region

We studied spicular jets over a plage area and derived their dynamic characteristics using Hinode Solar Optical Telescope (SOT) high-resolution images. The target plage region was near the west limb of the solar disk. This location permitted us to study the dynamics of spicular jets without the overlapping effect of spicular structures along the line of sight. In this work, to increase the ease with which we can identify spicules on the disk, we applied the image processing method `MadMax' developed by Koutchmy et al. (1989). It enhances fine, slender structures (like jets), over a diffuse background. We identified 169 spicules over the target plage. This sample permits us to derive statistically reliable results regarding spicular dynamics. The properties of plage spicules can be summarized as follows: (1) In a plage area, we clearly identified spicular jet features. (2) They were shorter in length than the quiet region limb spicules, and followed ballistic motion under constant deceleration. (3) The majority (80%) of the plage spicules showed the cycle of rise and retreat, while 10% of them faded out without a complete retreat phase. (4) The deceleration of the spicule was proportional to the velocity of ejection (i.e. the initial velocity).