Source author record

James E. Leake

James E. Leake 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

8works
3topics
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

8 published item(s)

preprint2026arXiv

Simulating the photospheric to coronal plasma using magnetohydrodynamic characteristics III: validation including gravity, flux emergence, and an eruption

Solar eruptions arise from instabilities or loss of equilibria in the solar atmosphere, but routinely inferring the precise magnetic and plasma properties that lead to eruptions is not currently practical using synoptic solar observations. Data driven simulations offer an appealing alternative. We test our boundary data-driven magnetohydrodynamic (MHD) approach, based on the method of characteristics, on a simulation that includes full MHD, a stratified atmosphere, and the emergence of a model solar magnetic active region, from the photosphere upwards. The driven simulation is tested against a larger, ab initio ``Ground Truth'' simulation that extends downward into the convection zone. Our driven simulation accurately reproduces the dynamic emergence of the active region above the photosphere, the formation of key topological features throughout the corona, and the subsequent eruption of mass and magnetic field. The total emerged energy matches to better than one percent, the ratio of emerged to eruptive energy is $\approx2\%$, and the actual values of each energy term agree to within $10\%$ between the two cases. Crucially, the data injection cadence, when properly scaled, matches the cadence of synoptic observations of the Sun's surface magnetic field, and is three to four orders of magnitude longer than the inherent CFL time step of the simulations. The stability of the code and fidelity of the results over an entire active region lifetime, from emergence to eruption, strongly suggests that our method will produce reliable results when driven using solar synoptic observations from existing and anticipated ground and spaced based observatories.

preprint2022arXiv

Impact of 3D Structure on Magnetic Reconnection

Results from 2.5D and 3D studies of the onset and development of the tearing instability are presented, using high fidelity resistive MHD simulations. A limited parameter study of the strength of the reconnecting field (or shear angle) was performed. An initially simple 1D equilibrium was used, consisting of a modified force-free current sheet, with periodic boundary conditions in all directions. In all cases, the linear and non-linear evolution led to a primary current sheet between two large flux ropes. The global reconnection rate during this later stage was analyzed in all simulations. It was found that in 2.5D the primary current sheet fragmented due to plasmoids, and as expected, the global reconnection rate, calculated using multiple methods, increases with the strength of the reconnecting field due to a stronger Alfvén speed. In 3D, the presence of interacting oblique modes of the tearing instability complicates the simple 2.5D picture, entangling the magnetic field of the inflow and introducing a negative effect on the reconnection rate. The two competing effects of stronger Alfvén speed and entangling, which both increase with the strength of the reconnecting field, resulted in a decrease in the reconnection rate with increasing reconnecting field. For all simulations, the 3D rates were less than in 2.5D, but suggest that as one goes to weak reconnecting field (or strong guide field), the system becomes more 2.5D like and the 2.5D and 3D rates converge. These results have relevance to situations like nano-flare heating and flare current sheets in the corona.

preprint2021arXiv

The Rise and Emergence of Untwisted Toroidal Flux Ropes on the Sun

Magnetic flux ropes (MFRs) rising buoyantly through the Sun's convection zone are thought to be subject to viscous forces preventing them from rising coherently. Numerous studies have suggested that MFRs require a minimum twist in order to remain coherent during their rise. Furthermore, even MFRs that get to the photosphere may be unable to successfully emerge into the corona unless they are at least moderately twisted, since the magnetic pressure gradient needs to overcome the weight of the photospheric plasma. To date, however, no lower limit has been placed on the critical minimum twist required for an MFR to rise coherently through the convection zone or emerge through the photosphere. In this paper, we simulate an untwisted toroidal MFR which is able to rise from the convection zone and emerge through the photosphere as an active region that resembles those observed on the Sun. We show that untwisted MFRs can remain coherent during their rise and then pile-up near the photosphere, triggering the undular instability, allowing the MFR to emerge through the photosphere. We propose that the toroidal geometry of our MFR is critical for its coherent rise. Upon emerging, a pair of lobes rises into the corona which interact and reconnect, resulting in a localized high speed jet. The resulting photospheric magnetogram displays the characteristic salt-and-pepper structure often seen in observations. Our major result is that MFRs need not be twisted to rise coherently through the convection zone and emerge through the photosphere.

preprint2014arXiv

Distribution of Electric Currents in Solar Active Regions

There has been a long-lasting debate on the question of whether or not electric currents in solar active regions are neutralized. That is, whether or not the main (or direct) coronal currents connecting the active region polarities are surrounded by shielding (or return) currents of equal total value and opposite direction. Both theory and observations are not yet fully conclusive regarding this question, and numerical simulations have, surprisingly, barely been used to address it. Here we quantify the evolution of electric currents during the formation of a bipolar active region by considering a three-dimensional magnetohydrodynamic simulation of the emergence of a sub-photospheric, current-neutralized magnetic flux rope into the solar atmosphere. We find that a strong deviation from current neutralization develops simultaneously with the onset of significant flux emergence into the corona, accompanied by the development of substantial magnetic shear along the active region's polarity inversion line. After the region has formed and flux emergence has ceased, the strong magnetic fields in the region's center are connected solely by direct currents, and the total direct current is several times larger than the total return current. These results suggest that active regions, the main sources of coronal mass ejections and flares, are born with substantial net currents, in agreement with recent observations. Furthermore, they support eruption models that employ pre-eruption magnetic fields containing such currents.

preprint2013arXiv

Simulations of Emerging Magnetic Flux. I: The Formation of Stable Coronal Flux Ropes

We present results from 3D visco-resistive magnetohydrodynamic (MHD) simulations of the emergence of a convection zone magnetic flux tube into a solar atmosphere containing a pre-existing dipole coronal field, which is orientated to minimize reconnection with the emerging field. We observe that the emergence process is capable of producing a coronal flux rope by the transfer of twist from the convection zone as found in previous simulations. We find that this flux rope is stable, with no evidence of a fast rise, and that its ultimate height in the corona is determined by the strength of the pre-existing dipole field. We also find that although the electric currents in the initial convection zone flux tube are almost perfectly neutralized, the resultant coronal flux rope carries a significant net current. These results suggest that flux tube emergence is capable of creating non-current-neutralized stable flux ropes in the corona, tethered by overlying potential fields, a magnetic configuration that is believed to be the source of coronal mass ejections.

preprint2012arXiv

Effect of ion-neutral collisions in simulations of emerging active regions

We present results of 2.5D numerical simulations of the emergence of sub-surface magnetic flux into the solar atmosphere, with emerging flux regions ranging from $10^{18}$ to $10^{21}$ Mx, representing both ephemeral and active regions. We include the presence of neutral Hydrogen in the governing equations, improve upon previous models by including the ionization in the equation of state, and use a more realistic convection zone model. We find that ionization and recombination of plasma during the rise of a convection zone flux tube reduces the rise speed of the tube's axis. The presence of neutral Hydrogen allows the effective flow of mass across fieldlines, by the addition of a Pedersen resistivity to the generalized Ohm's law, which dissipates current perpendicular to the magnetic field. This causes an increase of up to 10% in the amount of magnetic in-plane flux supplied to the corona and a reduction of up to 89% in the amount of sub-surface plasma brought up into the corona. However, it also reduces the amount of free magnetic energy supplied to the corona, and thus does not positively affect the likelihood of creating unstable coronal structures.

preprint2012arXiv

Multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma

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.

preprint2010arXiv

Tests of Dynamical Flux Emergence as a Mechanism for CME Initiation

Current coronal mass ejection (CME) models set their lower boundary to be in the lower corona. They do not calculate accurately the transfer of free magnetic energy from the convection zone to the magnetically dominated corona because they model the effects of flux emergence using kinematic boundary conditions or simply assume the appearance of flux at these heights. We test the importance of including dynamical flux emergence in CME modeling by simulating, in 2.5D, the emergence of sub-surface flux tubes into different coronal magnetic field configurations. We investigate how much free magnetic energy, in the form of shear magnetic field, is transported from the convection zone to the corona, and whether dynamical flux emergence can drive CMEs. We find that multiple coronal flux ropes can be formed during flux emergence, and although they carry some shear field into the corona, the majority of shear field is confined to the lower atmosphere. Less than 10% of the magnetic energy in the corona is in the shear field, and this, combined with the fact that the coronal flux ropes bring up significant dense material, means that they do not erupt. Our results have significant implications for all CME models which rely on the transfer of free magnetic energy from the lower atmosphere into the corona but which do not explicitly model this transfer. Such studies of flux emergence and CMEs are timely, as we have new capabilities to observe this with Hinode and SDO, and therefore to test the models against observations.