NASA's Parker Solar Probe (Parker) spacecraft reached its first perihelion of 35.7 solar radii on November 5th, 2018. To aid in mission planning, and in anticipation of the unprecedented measurements to be returned, in late October, we developed a three-dimensional magnetohydrodynamic (MHD) solution for the solar corona and inner heliosphere, driven by the then available observations of the Sun's photospheric magnetic field. Our model incorporates a wave-turbulence-driven (WTD) model to heat the corona. Here, we present our predictions for the structure of the solar corona and the likely {\it in situ} measurements that Parker will be returning over the next few months. We infer that, in the days prior to 1st Encounter, Parker was immersed in wind emanating from a positive-polarity equatorial coronal hole. During the encounter, however, field lines from the spacecraft mapped to another, negative-polarity equatorial coronal hole. Following the encounter, Parker was magnetically connected to the large, positive-polarity northern polar coronal hole. When the Parker data become available, these model results can be used to assist in their calibration and interpretation, and, additionally, provide a global context for interpreting the localized {\it in situ} measurements. In particular, we can identify what types of solar wind Parker encountered, what the underlying magnetic structure was, and how complexities in the orbital trajectory can be interpreted within a global, inertial frame. Ultimately, the measurements returned by Parker can be used to constrain current theories for heating the solar corona and accelerating the solar wind.

We describe the initial version of the Solar Particle Event (SPE) Threat Assessment Tool or STAT. STAT relies on elements of Corona-Heliosphere (CORHEL) and the Earth-Moon-Mars Radiation Environment Module (EMMREM), and allows users to investigate coronal mass ejection (CME) driven SPEs using coupled magnetohydrodynamic (MHD) and focused transport solutions. At the present time STAT focuses on modeling solar energetic particle (SEP) acceleration in and transport from the low corona, where the highest energy SEP events are generated. We illustrate STAT's capabilities with a model of the July 14, 2000 "Bastille Day" event, including innovative diagnostics for understanding the three-dimensional distribution of particle fluxes and their relation to the structure of the underlying CME driver. A preliminary comparison with NOAA GOES measurements is shown.

Slitless spectrometers can provide simultaneous imaging and spectral data over an extended field of view, thereby allowing rapid data acquisition for extended sources. In some instances, when the object is greatly extended or the spectral dispersion is too small, there may be locations in the focal plane where contributions from emission lines at different wavelengths contribute. It is then desirable to unfold the overlapped regions in order to isolate the contributions from the individual wavelengths. In this paper, we describe a method for such an unfolding, using an inversion technique developed for an extreme ultraviolet imaging spectrometer and coronagraph named the COronal Spectroscopic Imager in the EUV (COSIE). The COSIE spectrometer wavelength range (18.6 - 20.5 nm) contains a number of strong coronal emission lines and several density sensitive lines. We focus on optimizing the unfolding process to retrieve emission measure maps at constant temperature, maps of spectrally pure intensity in the Fe XII and Fe XIII lines and density maps based on both Fe XII and Fe XIII diagnostics.

Ion fractional charge states, measured in situ in the heliosphere, depend on the properties of the plasma in the inner corona. As the ions travel outward in the solar wind and the electron density drops, the charge states remain essentially unaltered or "frozen in". Thus they can provide a powerful constraint on heating models of the corona and acceleration of the solar wind. We have implemented non-equilibrium ionization calculations into a 1D wave-turbulence-driven (WTD) hydrodynamic solar wind model and compared modeled charge states with the Ulysses 1994-5 in situ measurements. We have found that modeled charge state ratios of $C^{6+}/C^{5+}$ and $O^{7+}/O^{6+}$, among others, were too low compared with Ulysses measurements. However, a heuristic reduction of the plasma flow speed has been able to bring the modeled results in line with observations, though other ideas have been proposed to address this discrepancy. We discuss implications of our results and the prospect of including ion charge state calculations into our 3D MHD model of the inner heliosphere.

GPU accelerators have had a notable impact on high-performance computing across many disciplines. They provide high performance with low cost/power, and therefore have become a primary compute resource on many of the largest supercomputers. Here, we implement multi-GPU acceleration into our Solar MHD code (MAS) using OpenACC in a fully portable, single-source manner. Our preliminary implementation is focused on MAS running in a reduced physics "zero-beta" mode. While valuable on its own, our main goal is to pave the way for a full physics, thermodynamic MHD implementation. We describe the OpenACC implementation methodology and challenges. "Time-to-solution" performance results of a production-level flux rope eruption simulation on multi-CPU and multi-GPU systems are shown. We find that the GPU-accelerated MAS code has the ability to run "zero-beta" simulations on a single multi-GPU server at speeds previously requiring multiple CPU server-nodes of a supercomputer.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

Observations of coronal loops have identified several common loop characteristics, including that loops appear to cool and have higher than expected densities. Two potential heating scenarios have been suggested to explain these observations. One scenario is that the loops are heated by a series of small-scale impulsive heating events, or nanoflares. Another hypothesis is that the heating is quasi-steady and highly-stratified, i.e., ``footpoint heating'. The goal of this paper is to identify observables that can be used to differentiate between these two heating scenarios. For footpoint heating, we vary the heating magnitude and stratification, for impulsive heating, we vary the heating magnitude. We use one-dimensional hydrodynamic codes to calculate the resulting temperature and density evolution and expected lightcurves in four channels of AIA and one channel of XRT. We consider two principal diagnostics: the time lag between the appearance of the loop in two different channels, and the ratio of the peak intensities of the loop in the two channels. We find that 1) footpoint heating can predict longer time lags than impulsive heating, 2) footpoint heating can predict zero or negative time lags, 3) the intensity ratio expected from impulsive heating is confined to a narrow range, while footpoint heating predicts a wider range of intensity ratios, and 4) the range of temperatures expected in impulsive heating is broader than the range of temperatures expected in footpoint heating. This preliminary study identifies observables that may be useful in discriminating between heating models in future work.

We report SDO/AIA observations of an extraordinary global extreme ultraviolet (EUV) wave triggered by the X8.2+ flare-CME eruption on 2017 September 10. This was one of the best EUV waves ever observed with modern instruments, yet it was likely the last one of such magnitudes of Solar Cycle 24 as the Sun heads toward the minimum. Its remarkable characteristics include the following. (1) The wave was observed, for the first time, to traverse the full-Sun corona over the entire visible solar disk and off-limb circumference, manifesting a truly global nature, owing to its exceptionally large amplitude, e.g., with EUV enhancements by up to 300% at 1.1 Rsun from the eruption. (2) This leads to strong transmissions (in addition to commonly observed reflections) in and out of both polar coronal holes, which are usually devoid of EUV waves. It has elevated wave speeds >2000 km/s within them, consistent with the expected higher fast-mode magnetosonic wave speeds. The coronal holes essentially serve as new "radiation centers" for the waves being refracted out of them, which then travel toward the equator and collide head-on, causing additional EUV enhancements. (3) The wave produces significant compressional heating to local plasma upon its impact, indicated by long-lasting EUV intensity changes and differential emission measure increases at higher temperatures (e.g., log T=6.2) accompanied by decreases at lower temperatures (e.g., log T=6.0). These characteristics signify the potential of such EUV waves for novel magnetic and thermal diagnostics of the solar corona {\it on global scales}.

Seven different models are applied to the same problem of simulating the Sun's coronal magnetic field during the solar eclipse on 2015 March 20. All of the models are non-potential, allowing for free magnetic energy, but the associated electric currents are developed in significantly different ways. This is not a direct comparison of the coronal modelling techniques, in that the different models also use different photospheric boundary conditions, reflecting the range of approaches currently used in the community. Despite the significant differences, the results show broad agreement in the overall magnetic topology. Among those models with significant volume currents in much of the corona, there is general agreement that the ratio of total to potential magnetic energy should be approximately 1.4. However, there are significant differences in the electric current distributions; while static extrapolations are best able to reproduce active regions, they are unable to recover sheared magnetic fields in filament channels using currently available vector magnetogram data. By contrast, time-evolving simulations can recover the filament channel fields at the expense of not matching the observed vector magnetic fields within active regions. We suggest that, at present, the best approach may be a hybrid model using static extrapolations but with additional energization informed by simplified evolution models. This is demonstrated by one of the models.

We describe and analyze observations by the Solar Dynamics Observatory of the emergence of a small, bipolar active region within an area of unipolar magnetic flux that was surrounded by a circular, quiescent filament. Within only eight hours of the start of the emergence, a partial splitting of the filament and two consecutive coronal mass ejections took place. We argue that all three dynamic events occurred as a result of particular magnetic-reconnection episodes between the emerging bipole and the pre-existing coronal magnetic field. In order to substantiate our interpretation, we consider three-dimensional magnetohydrodynamic simulations that model the emergence of magnetic flux in the vicinity of a large-scale coronal flux rope. The simulations qualitatively reproduce most of the reconnection episodes suggested by the observations; as well as the filament-splitting, the first eruption, and the formation of sheared/twisted fields that may have played a role in the second eruption. Our results suggest that the position of emerging flux with respect to the background magnetic configuration is a crucial factor for the resulting evolution, while previous results suggest that parameters such as the orientation or the amount of emerging flux are important as well. This poses a challenge for predicting the onset of eruptions that are triggered by flux emergence, and it calls for a detailed survey of the relevant parameter space by means of numerical simulations.

Solar eruptions are the main driver of space-weather disturbances at the Earth. Extreme events are of particular interest, not only because of the scientific challenges they pose, but also because of their possible societal consequences. Here we present a magnetohydrodynamic (MHD) simulation of the 14 July 2000 Bastille Day eruption, which produced a very strong geomagnetic storm. After constructing a thermodynamic MHD model of the corona and solar wind, we insert a magnetically stable flux rope along the polarity inversion line of the eruption's source region and initiate the eruption by boundary flows. More than 10^33 ergs of magnetic energy are released in the eruption within a few minutes, driving a flare, an EUV wave, and a coronal mass ejection (CME) that travels in the outer corona at about 1500 km/s, close to the observed speed. We then propagate the CME to Earth, using a heliospheric MHD code. Our simulation thus provides the opportunity to test how well in situ observations of extreme events are matched if the eruption is initiated from a stable magnetic-equilibrium state. We find that the flux-rope center is very similar in character to the observed magnetic cloud, but arrives about 8.5 hours later and about 15 degrees too far to the North, with field strengths that are too weak by a factor of about 1.6. The front of the flux rope is highly distorted, exhibiting localized magnetic-field concentrations as it passes 1 AU. We discuss these properties with regard to the development of space-weather predictions based on MHD simulations of solar eruptions.

Many existing models assume that magnetic flux ropes play a key role in solar flares and coronal mass ejections (CMEs). It is therefore important to develop efficient methods for constructing flux-rope configurations constrained by observed magnetic data and the morphology of the pre-eruptive source region. For this purpose, we have derived and implemented a compact analytical form that represents the magnetic field of a thin flux rope with an axis of arbitrary shape and circular cross-sections. This form implies that the flux rope carries axial current $I$ and axial flux $F$, so that the respective magnetic field is the curl of the sum of axial and azimuthal vector potentials proportional to $I$ and $F$, respectively. We expressed the vector potentials in terms of modified Biot-Savart laws whose kernels are regularized at the axis in such a way that, when the axis is straight, these laws define a cylindrical force-free flux rope with a parabolic profile for the axial current density. For the cases we have studied so far, we determined the shape of the rope axis by following the polarity inversion line of the eruptions' source region, using observed magnetograms. The height variation along the axis and other flux-rope parameters are estimated by means of potential field extrapolations. Using this heuristic approach, we were able to construct pre-eruption configurations for the 2009 February 13 and 2011 October 1 CME events. These applications demonstrate the flexibility and efficiency of our new method for energizing pre-eruptive configurations in simulations of CMEs.

The Miniature X-ray Solar Spectrometer (MinXSS) is a 3 Unit (3U) CubeSat designed for a 3-month mission to study solar soft X-ray spectral irradiance. The first of the two flight models was deployed from the International Space Station in 2016 May and operated for one year before its natural deorbiting. This was the first flight of the Blue Canyon Technologies XACT 3-axis attitude determination and control system -- a commercially available, high-precision pointing system. We characterized the performance of the pointing system on orbit including performance at low altitudes where drag torque builds up. We found that the pointing accuracy was 0.0042\degree\ - 0.0117\degree\ (15$''$ - 42$''$, 3$\sigma$, axis dependent) consistently from 190 km - 410 km, slightly better than the specification sheet states. Peak-to-peak jitter was estimated to be 0.0073\degree\ (10 s$^{-1}$) - 0.0183\degree\ (10 s$^{-1}$) (26$''$ (10 s$^{-1}$) - 66$''$ (10 s$^{-1}$), 3$\sigma$). The system was capable of dumping momentum until an altitude of 185 km. We found small amounts of sensor degradation in the star tracker and coarse sun sensor. Our mission profile did not require high-agility maneuvers so we are unable to characterize this metric. Without a GPS receiver, it was necessary to periodically upload ephemeris information to update the orbit propagation model and maintain pointing. At 400 km, these uploads were required once every other week. At $\sim$270 km, they were required every day. We also characterized the power performance of our electric power system, which includes a novel pseudo-peak power tracker -- a resistor that limited the current draw from the battery on the solar panels. With 19 30\% efficient solar cells and an 8 W system load, the power balance had 65\% of margin on orbit. We present several recommendations to other CubeSat programs throughout.

The heliospheric magnetic field is of pivotal importance in solar and space physics. The field is rooted in the Sun's photosphere, where it has been observed for many years. Global maps of the solar magnetic field based on full disk magnetograms are commonly used as boundary conditions for coronal and solar wind models. Two primary observational constraints on the models are (1) the open field regions in the model should approximately correspond to coronal holes observed in emission, and (2) the magnitude of the open magnetic flux in the model should match that inferred from in situ spacecraft measurements. In this study, we calculate both MHD and PFSS solutions using fourteen different magnetic maps produced from five different types of observatory magnetograms, for the time period surrounding July, 2010. We have found that for all of the model/map combinations, models that have coronal hole areas close to observations underestimate the interplanetary magnetic flux, or, conversely, for models to match the interplanetary flux, the modeled open field regions are larger than coronal holes observed in EUV emission. In an alternative approach, we estimate the open magnetic flux entirely from solar observations by combining automatically detected coronal holes for Carrington rotation 2098 with observatory synoptic magnetic maps. This approach also underestimates the interplanetary magnetic flux. Our results imply that either typical observatory maps underestimate the Sun's magnetic flux, or a significant portion of the open magnetic flux is not rooted in regions that are obviously dark in EUV and X-ray emission.

For almost 20 years the physical nature of globally propagating waves in the solar corona (commonly called "EIT waves") has been controversial and subject to debate. Additional theories have been proposed over the years to explain observations that did not fit with the originally proposed fast-mode wave interpretation. However, the incompatibility of observations made using the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory with the fast-mode wave interpretation was challenged by differing viewpoints from the twin Solar Terrestrial Relations Observatory spacecraft and higher spatial/temporal resolution data from the Solar Dynamics Observatory. In this article, we reexamine the theories proposed to explain "EIT waves" to identify measurable properties and behaviours that can be compared to current and future observations. Most of us conclude that "EIT waves" are best described as fast-mode large-amplitude waves/shocks that are initially driven by the impulsive expansion of an erupting coronal mass ejection in the low corona.

To simulate the energy balance of coronal plasmas on macroscopic scales, we often require the specification of the coronal heating mechanism in some functional form. To go beyond empirical formulations and to build a more physically motivated heating function, we investigate the wave-turbulence-driven (WTD) phenomenology for the heating of closed coronal loops. Our implementation is designed to capture the large-scale propagation, reflection, and dissipation of wave turbulence along a loop. The parameter space of this model is explored by solving the coupled WTD and hydrodynamic evolution in 1D for an idealized loop. The relevance to a range of solar conditions is also established by computing solutions for over one hundred loops extracted from a realistic 3D coronal field. Due to the implicit dependence of the WTD heating model on loop geometry and plasma properties along the loop and at the footpoints, we find that this model can significantly reduce the number of free parameters when compared to traditional empirical heating models, and still robustly describe a broad range of quiet-sun and active region conditions. The importance of the self-reflection term in producing relatively short heating scale heights and thermal nonequilibrium cycles is also discussed.

We describe a method for the automatic mapping of coronal holes (CH) using simultaneous multi-instrument EUV imaging data. Synchronized EUV images from STEREO/EUVI A&B 195A and SDO/AIA 193A are preprocessed, including PSF deconvolution and the application of data-derived intensity corrections that account for center-to-limb variations (limb brightening) and inter-instrument intensity normalization. We systematically derive a robust limb-brightening correction that takes advantage of unbiased long-term averages of data and respects the physical nature of the problem. The new preprocessing greatly assists in CH detection, allowing for the use of a simplified variable-connectivity two-threshold region growing image segmentation algorithm to obtain consistent detection results. We generate synchronic EUV and CH maps, and show a preliminary analysis of CH evolution. Several data and code products are made available to the community (www.predsci.com/chd): For the period of this study (06/10/2010 to 08/18/14) we provide synchronic EUV and coronal hole map data at 6-hour cadence, data-derived limb-brightening corrections for STEREO/EUVI A&B 195A and SDO/AIA 193A, and inter-instrument correction factors to equate their intensities. We also provide the coronal hole image segmentation code module (ezseg) implemented in both FORTRAN OpenMP and GPU-accelerated C-CUDA. A complete implementation of our coronal hole detection pipeline in the form of a ready-to-use MATLAB driver script euv2chm utilizing ezseg is also made available.

Quasi-periodic, fast-mode, propagating wave trains (QFPs) are a new observational phenomenon recently discovered in the solar corona by the Solar Dynamics Observatory with extreme ultraviolet (EUV) imaging observations. They originate from flares and propagate at speeds up to ~2000 km/s within funnel-shaped waveguides in the wakes of coronal mass ejections (CMEs). QFPs can carry sufficient energy fluxes required for coronal heating during their occurrences. They can provide new diagnostics for the solar corona and their associated flares. We present recent observations of QFPs focusing on their spatio-temporal properties, temperature dependence, and statistical correlation with flares and CMEs. Of particular interest is the 2010-Aug-01 C3.2 flare with correlated QFPs and drifting zebra and fiber radio bursts, which might be different manifestations of the same fast-mode wave trains. We also discuss the potential roles of QFPs in accelerating and/or modulating the solar wind.

Although it is widely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process, the detailed mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still not fully understood. Cranmer et al. (2007) developed a sophisticated, 1D, time-steady model of the solar wind with turbulence dissipation. By varying the coronal magnetic field, they obtain, for a single choice of wave properties, a realistic range of slow and fast wind conditions with a sharp latitudinal transition between the two streams. Using a 1D, time-dependent model of the solar wind of Lionello et al. (2014), which incorporates turbulent dissipation of Alfv\'en waves to provide heating and acceleration of the plasma, we have explored a similar configuration, obtaining qualitatively equivalent results. However, our calculations suggest that the rapid transition between slow and fast wind suggested by this 1D model may be disrupted in multidimensional MHD simulations by the requirement of transverse force balance.

EUV images of Comet Lovejoy (C/2011 W3) from the AIA show striations related to the magnetic field structure in both open and closed magnetic regions. The brightness contrast implies coronal density contrasts of at least a factor of 6 between neighboring flux tubes over scales of a few thousand km. These density structures imply variations in the Alfven speed on a similar scale. They will drastically affect the propagation and dissipation of Alfven waves, and that should be taken into account in models of coronal heating and solar wind acceleration. In each striation, the cometary emission moves along the magnetic field and broadens with time. The speed and the rate of broadening are related to the parallel and perpendicular components of the velocities of the cometary neutrals when they become ionized. We use an MHD model of the coronal magnetic field and the theory of pickup ions to compare the measurements with theoretical predictions, in particular with the energy lost to Alfven waves as the cometary ions isotropize.

Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. \citet{2010ApJ...708L.116V} presented a model for heating and accelerating the solar wind based on the turbulent dissipation of Alfv\'en waves. We first use a time-dependent model of the solar wind to reproduce one of \citeauthor{2010ApJ...708L.116V}'s solutions; then we extend its application to the case when the energy equation includes thermal conduction and radiation losses, and the upper chromosphere is part of the computational domain. Using this model, we explore parameter space and describe the characteristics of a fast-solar-wind solution. We discuss how this formulation may be applied to a 3D MHD model of the coron a and solar wind \citep{2009ApJ...690..902L}.

We study the influence of the spatial resolution on scales of $5\deg$ and smaller of solar surface magnetic field maps on global magnetohydrodynamic solar wind models, and on a model of coronal heating and X-ray emission. We compare the solutions driven by a low-resolution Wilcox Solar Observatory magnetic map, the same map with spatial resolution artificially increased by a refinement algorithm, and a high-resolution Solar and Heliospheric Observatory Michelson Doppler Imager map. We find that both the wind structure and the X-ray morphology are affected by the fine-scale surface magnetic structure. Moreover, the X-ray morphology is dominated by the closed loop structure between mixed polarities on smaller scales and shows significant changes between high and low resolution maps. We conclude that three-dimensional modeling of coronal X-ray emission has greater surface magnetic field spatial resolution requirements than wind modeling, and can be unreliable unless the dominant mixed polarity magnetic flux is properly resolved.

We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107.

We present a numerical investigation of the coronal evolution of a coronal mass ejection (CME) on 2005 August 22 using a 3-D thermodynamics magnetohydrodynamic model, the SWMF. The source region of the eruption was anemone active region (AR) 10798, which emerged inside a coronal hole. We validate our modeled corona by producing synthetic extreme ultraviolet (EUV) images, which we compare to EIT images. We initiate the CME with an out-of-equilibrium flux rope with an orientation and chirality chosen in agreement with observations of a H-alpha filament. During the eruption, one footpoint of the flux rope reconnects with streamer magnetic field lines and with open field lines from the adjacent coronal hole. It yields an eruption which has a mix of closed and open twisted field lines due to interchange reconnection and only one footpoint line-tied to the source region. Even with the large-scale reconnection, we find no evidence of strong rotation of the CME as it propagates. We study the CME deflection and find that the effect of the Lorentz force is a deflection of the CME by about 3 deg/Rsun towards the East during the first 30 minutes of the propagation. We also produce coronagraphic and EUV images of the CME, which we compare with real images, identifying a dimming region associated with the reconnection process. We discuss the implication of our results for the arrival at Earth of CMEs originating from the limb and for models to explain the presence of open field lines in magnetic clouds.

We present an inquiry lab activity on Circuit Design that was conducted in Fall 2009 with first-year community college students majoring in Electrical Engineering Technology. This inquiry emphasized the use of engineering process skills, including circuit assembly and problem solving, while learning technical content. Content goals of the inquiry emphasized understanding voltage dividers (Kirchoff's voltage law) and analysis and optimization of resistive networks (Thevenin equivalence). We assumed prior exposure to series and parallel circuits and Ohm's law (the relationship between voltage, current, and resistance) and designed the inquiry to develop these skills. The inquiry utilized selection of engineering challenges on a specific circuit (the Wheatstone Bridge) to realize these learning goals. Students generated questions and observations during the starters, which were categorized into four engineering challenges or design goals. The students formed teams and chose one challenge to focus on during the inquiry. We created a rubric for summative assessment which helped to clarify and solidify project goals while designing the inquiry and aided in formative assessment during the activity. After describing implementation, we compare and contrast engineering-oriented inquiry design as opposed to activities geared toward science learning.

We perform a numerical study of the evolution of a Coronal Mass Ejection (CME) and its interaction with the coronal magnetic field based on the May 12, 1997, CME event using a global MagnetoHydroDynamic (MHD) model for the solar corona. The ambient solar wind steady-state solution is driven by photospheric magnetic field data, while the solar eruption is obtained by superimposing an unstable flux rope onto the steady-state solution. During the initial stage of CME expansion, the core flux rope reconnects with the neighboring field, which facilitates lateral expansion of the CME footprint in the low corona. The flux rope field also reconnects with the oppositely orientated overlying magnetic field in the manner of the breakout model. During this stage of the eruption, the simulated CME rotates counter-clockwise to achieve an orientation that is in agreement with the interplanetary flux rope observed at 1 AU. A significant component of the CME that expands into interplanetary space comprises one of the side lobes created mainly as a result of reconnection with the overlying field. Within 3 hours, reconnection effectively modifies the CME connectivity from the initial condition where both footpoints are rooted in the active region to a situation where one footpoint is displaced into the quiet Sun, at a significant distance ($\approx 1R_\odot$) from the original source region. The expansion and rotation due to interaction with the overlying magnetic field stops when the CME reaches the outer edge of the helmet streamer belt, where the field is organized on a global scale. The simulation thus offers a new view of the role reconnection plays in rotating a CME flux rope and transporting its footpoints while preserving its core structure.

In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model; one set in the upper transition region (the Radiative Energy Balance model), as well as a uniform chromospheric condition where the transition region can be modeled in its entirety. Via observation synthesis from model results and the subsequent comparison to full sun extreme ultraviolet (EUV) and soft X-Ray observations of Carrington Rotation (CR) 1913 centered on Aug 27, 1996, we demonstrate the need for these additional considerations when using global MHD models to describe the unique conditions in the low corona. Through multiple simulations we examine ability of the LC model to asses and discriminate between coronal heating models, and find that a relative simple empirical heating model is adequate in reproducing structures observed in the low corona. We show that the interplay between coronal heating and electron heat conduction provides significant feedback onto the 3D magnetic topology in the low corona as compared to a potential field extrapolation, and that this feedback is largely dependent on the amount of mechanical energy introduced into the corona.