Researcher profile

S. A. Gilchrist

S. A. Gilchrist contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.SR

Trust snapshot

Quick read

Trust 21 - Emerging

6works

0followers

1topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2016arXiv

A Fixed-point Scheme for the Numerical Construction of Magnetohydrostatic Atmospheres in Three Dimensions

Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin (Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular.

preprint2015arXiv

Critical decay index at the onset of solar eruptions

Magnetic flux ropes are topological structures consisting of twisted magnetic field lines that globally wrap around an axis. The torus instability model predicts that a magnetic flux rope of major radius $R$ undergoes an eruption when its axis reaches a location where the decay index $-d(\ln B_{ex})/d(\ln R)$ of the ambient magnetic field $B_{ex}$ is larger than a critical value. In the current-wire model, the critical value depends on the thickness and time-evolution of the current channel. We use magneto-hydrodynamic (MHD) simulations to investigate if the critical value of the decay index at the onset of the eruption is affected by the magnetic flux rope's internal current profile and/or by the particular pre-eruptive photospheric dynamics. The evolution of an asymmetric, bipolar active region is driven by applying different classes of photospheric motions. We find that the critical value of the decay index at the onset of the eruption is not significantly affected by either the pre-eruptive photospheric evolution of the active region or by the resulting different magnetic flux ropes. As in the case of the current-wire model, we find that there is a `critical range' $ [1.3-1.5]$, rather than a `critical value' for the onset of the torus instability. This range is in good agreement with the predictions of the current-wire model, despite the inclusion of line-tying effects and the occurrence of tether-cutting magnetic reconnection.

preprint2015arXiv

The Influence of Spatial Resolution on Nonlinear Force-Free Modeling

The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting.

preprint2013arXiv

Nonlinear force-free modeling of the corona in spherical coordinates

We present a code for solving the nonlinear force-free equations in spherical polar geometry, with the motivation of modeling the magnetic field in the corona. The code is an implementation of the Grad-Rubin method. Our method is applicable to a spherical domain of arbitrary angular size. The implementation is based on a global spectral representation for the magnetic field which makes no explicit assumptions about the form of the magnetic field at the transverse boundaries of the domain. We apply the code to a bipolar test case with analytic boundary conditions, and we demonstrate the convergence of the Grad-Rubin method, and the self-consistency of the resulting numerical solution.

preprint2012arXiv

A Magnetostatic Grad-Rubin Code for Coronal Magnetic Field Extrapolations

The coronal magnetic field cannot be directly observed, but in principle it can be reconstructed from the comparatively well observed photospheric magnetic field. A popular approach uses a nonlinear force-free model. Non-magnetic forces at the photosphere are significant meaning the photospheric data are inconsistent with the force-free model, and this causes problems with the modeling (De Rosa et al., Astrophys. J. 696, 1780, 2009). In this paper we present a numerical implementation of the Grad-Rubin method for reconstructing the coronal magnetic field using a magnetostatic model. This model includes a pressure force and a non-zero magnetic Lorentz force. We demonstrate our implementation on a simple analytic test case and obtain the speed and numerical error scaling as a function of the grid size.

preprint2011arXiv

The free energy of NOAA active region AR 11029

The NOAA active region AR 11029 was a small but highly active sunspot region which produced 73 GOES soft X-ray flares. The flares appear to show a departure from the well known power-law frequency-size distribution. Specifically, too few GOES C-class and no M-class flares were observed by comparison with a power-law distribution (Wheatland in Astrophys. J. 710, 1324, 2010). This was conjectured to be due to the region having insufficient magnetic energy to power large events. We construct nonlinear force-free extrapolations of the coronal magnetic field of active region AR 11029 using data taken on 24 October by the SOLIS Vector-SpectroMagnetograph (SOLIS/VSM), and data taken on 27 October by the Hinode Solar Optical Telescope SpectroPolarimeter (Hinode/SP). Force-free modeling with photospheric magnetogram data encounters problems because the magnetogram data are inconsistent with a force-free model, and we employ a recently developed `self-consistency' procedure which addresses this and accommodates uncertainties in the boundary data (Wheatland and Regnier in Astrophys. J. 700, L88, 2009). We calculate the total energy and free energy of the self-consistent solution and find that the free energy was 4x10^29 erg on 24 October, and 7x10^31 erg on 27 October. An order of magnitude scaling between RHESSI non-thermal energy and GOES peak X-ray flux is established from a sample of flares from the literature and is used to estimate flare energies from observed GOES peak X-ray flux. Based on the scaling, we conclude that the estimated free energy of AR 11029 on 27 October when the flaring rate peaked is sufficient to power M-class or X-class flares, and hence the modeling does not appear to support the hypothesis that the absence of large flares is due to the region having limited energy.