The Wang-Sheeley-Arge (WSA)-ENLIL+Cone model is used extensively in space weather operations world-wide to model CME propagation. As such, it is important to assess its performance. We present validation results of the WSA-ENLIL+Cone model installed at the Community Coordinated Modeling Center (CCMC) and executed in real-time by the CCMC space weather team. CCMC uses the WSA-ENLIL+Cone model to predict CME arrivals at NASA missions throughout the inner heliosphere. In this work we compare model predicted CME arrival-times to in-situ ICME leading edge measurements at STEREO-A, STEREO-B, and Earth (Wind and ACE) for simulations completed between March 2010-December 2016 (over 1,800 CMEs). We report hit, miss, false alarm, and correct rejection statistics for all three locations. For all predicted CME arrivals, the hit rate is 0.5, and the false alarm rate is 0.1. For the 273 events where the CME was predicted to arrive at Earth, STEREO-A, or STEREO-B, and was actually observed (hit event), the mean absolute arrival-time prediction error was 10.4 +/- 0.9 hours, with a tendency to early prediction error of -4.0 hours. We show the dependence of the arrival-time error on CME input parameters. We also explore the impact of the multi-spacecraft observations used to initialize the model CME inputs by comparing model verification results before and after the STEREO-B communication loss (since September 2014) and STEREO-A sidelobe operations (August 2014-December 2015). There is an increase of 1.7 hours in the CME arrival time error during single, or limited two-viewpoint periods, compared to the three-spacecraft viewpoint period. This trend would apply to a future space weather mission at L5 or L4 as another coronagraph viewpoint to reduce CME arrival time errors compared to a single L1 viewpoint.

On 7 January 2014 an X1.2 flare and CME with a radial speed $\approx$2500 km s$^{-1}$ was observed from near an active region close to disk center. This led many forecasters to estimate a rapid arrival at Earth ($\approx$36 hours) and predict a strong geomagnetic storm. However, only a glancing CME arrival was observed at Earth with a transit time of $\approx$49 hours and a $K_{\rm P}$ geomagnetic index of only $3-$. We study the interplanetary propagation of this CME using the ensemble Wang-Sheeley-Arge (WSA)-ENLIL+Cone model, that allows a sampling of CME parameter uncertainties. We explore a series of simulations to isolate the effects of the background solar wind solution, CME shape, tilt, location, size, and speed, and the results are compared with observed in-situ arrivals at Venus, Earth, and Mars. Our results show that a tilted ellipsoid CME shape improves the initial real-time prediction to better reflect the observed in-situ signatures and the geomagnetic storm strength. CME parameters from the Graduated Cylindrical Shell model used as input to WSA--ENLIL+Cone, along with a tilted ellipsoid cloud shape, improve the arrival-time error by 14.5, 18.7, 23.4 hours for Venus, Earth, and Mars respectively. These results highlight that CME orientation and directionality with respect to observatories play an important role in understanding the propagation of this CME, and for forecasting other glancing CME arrivals. This study also demonstrates the importance of three-dimensional CME fitting made possible by multiple viewpoint imaging.

Ensemble modeling of CMEs provides a probabilistic forecast of CME arrival time which includes an estimation of arrival time uncertainty from the spread and distribution of predictions and forecast confidence in the likelihood of CME arrival. The real-time ensemble modeling of CME propagation uses the WSA-ENLIL+Cone model installed at the CCMC and executed in real-time. The current implementation evaluates the sensitivity of WSA-ENLIL+Cone model simulations of CME propagation to initial CME parameters. We discuss the results of real-time ensemble simulations for a total of 35 CME events between January 2013 - July 2014. For the 17 events where the CME was predicted to arrive at Earth, the mean absolute arrival time prediction error was 12.3 hours, which is comparable to the errors reported in other studies. For predictions of CME arrival at Earth the correct rejection rate is 62% and the false-alarm rate is 38%. The arrival time was within the range of the ensemble arrival predictions for 8 out of 17 events. The Brier Score for CME arrival predictions is 0.15 (where 1 is a perfect forecast), indicating that on average, the predicted likelihood of CME arrival is fairly accurate. The reliability of ensemble CME arrival predictions is heavily dependent on the initial distribution of CME input parameters, particularly the median and spread. Preliminary analysis of the probabilistic forecasts suggests undervariability, indicating that these ensembles do not sample a wide enough spread in CME input parameters. Prediction errors can also arise from ambient model parameters, the accuracy of the solar wind background derived from coronal maps, or other model limitations. Finally, predictions of the Kp geomagnetic index differ from observed values by less than one for 11 out of 17 of the ensembles and Kp prediction errors computed from the mean predicted Kp show a mean absolute error of 1.3.

The solar coronal magnetic field produces solar activity, including extremely energetic solar flares and coronal mass ejections (CMEs). Knowledge of the structure and evolution of the magnetic field of the solar corona is important for investigating and understanding the origins of space weather. Although the coronal field remains difficult to measure directly, there is considerable interest in accurate modeling of magnetic fields in and around sunspot regions on the Sun using photospheric vector magnetograms as boundary data. In this work, we investigate effects of the size of the domain chosen for coronal magnetic field modeling on resulting model solution. We apply spherical Optimization procedure to vector magnetogram data of Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO) with four Active Region observed on 09 March 2012 at 20:55UT. The results imply that quantities like magnetic flux density, electric current density and free magnetic energy density of ARs of interest are significantly different from the corresponding quantities obtained in the same region within the wider field of view. The difference is even more pronounced in the regions where there are connections to outside the domain.

Solar eruptive phenomena, like flares and coronal mass ejections(CMEs) are governed by magnetic fields. To describe the structure of these phenomena one needs information on the magnetic flux density and the electric current density vector components in three dimensions throughout the atmosphere. However, current spectro-polarimetric measurements typically limit the determination of the vector magnetic field only to the photosphere. Therefore, there is considerable interest in accurate modeling of the solar coronal magnetic field using photospheric vector magnetograms as boundary data. In this work, we model the coronal magnetic field for global solar atmosphere using a nonlinear force-free field(NLFFF) extrapolation codes implemented to a synoptic maps of photospheric vector magnetic field synthesized from Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) as boundary condition. Using the resulting three dimensional magnetic field, we calculate the three dimensional electric current density and free magnetic energy throughout the solar atmosphere for Carrington rotation 2124. We found that spatially, the low-lying, current-carrying core field demonstrates strong concentration of free energy in the AR core, from the photosphere to the lower corona; the coronal field becomes slightly more sheared in the lowest layer and it relaxes to the potential field configuration with height. The free energy density appears largely co-spatial with the electric current distribution.

The magnetic field permeating the solar atmosphere is generally thought to provide the energy for much of the activity seen in the solar corona, such as flares, coronal mass ejections (CMEs), etc. To overcome the unavailability of coronal magnetic field measurements, photospheric magnetic field vector data can be used to reconstruct the coronal field. Currently there are several modelling techniques being used to calculate three-dimension of the field lines into the solar atmosphere. For the first time, synoptic maps of photospheric vector magnetic field synthesized from Vector Spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) are used to model the coronal magnetic field and estimate free magnetic energy in the global scale. The free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. We solve the nonlinear force-free field equations using optimization principle in spherical geometry. The resulting three-dimensional magnetic fields are used to estimate the magnetic free energy content E_{free}=E_{nlfff}-E_{pot}, i.e., the difference of the magnetic energies between the nonpotential field and the potential field in the global solar corona. For comparison, we overlay the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the Atmospheric Imaging Assembly on board SDO. For a single Carrington rotation 2121, we find that the global NLFFF magnetic energy density is 10.3% higher than the potential one. Most of this free energy is located in active regions.

The magnetic fields in the solar atmosphere structure the plasma, store free magnetic energy and produce a wide variety of active solar phenomena, like flare and coronal mass ejections(CMEs). The distribution and strength of magnetic fields are routinely measured in the solar surface(photosphere). Therefore, there is considerable interest in accurately modeling the 3D structure of the coronal magnetic field using photospheric vector magnetograms. Knowledge of the 3D structure of magnetic field lines also help us to interpret other coronal observations, e.g., EUV images of the radiating coronal plasma. Nonlinear force-free field (NLFFF) models are thought to be viable tools for those task. Usually those models use Cartesian geometry. However, the spherical nature of the solar surface cannot be neglected when the field of view is large. In this work, we model the coronal magnetic field above multiple active regions using NLFFF extrapolation code using vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/ Vector Spectromagnetograph (VSM) as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with their respective coronal EUV-images from the Atmospheric Imaging Assembly (SDO/AIA) observed on October 11, 2011 and November 13, 2012. This study has found that the NLFFF model in spherical geometry reconstructs the magnetic configurations for several active regions which agrees with observations. During October 11, 2011 observation, there are substantial number of trans-equatorial loops carrying electric current.

Measurements of magnetic fields and electric currents in the pre-eruptive corona are crucial to study solar eruptive phenomena, like flare and coronal mass ejections(CMEs). However, spectro-polarimetric measurements of certain photospheric lines permit a determination of the vector magnetic field at the photosphere. Thus, substantial collection of magnetograms relate to the photospheric surface field only. Numerical modeling is carried out by applying state-of-the-art nonlinear force-free field (NLFFF) reconstruction. Cartesian nonlinear force-free field (NLFFF) codes are not well suited for larger domains, since the spherical nature of the solar surface cannot be neglected when the field of view is large. One of the most significant results of Solar Dynamic Observatory (SDO) mission to date has been repeated observations of large, almost global scale events in which large scale connection between active regions may play fundamental role. Therefore, it appears prudent to implement a NLFFF procedure in spherical geometry for use when large scale boundary data are available, such as from the Helioseismic and Magnetic Imager (HMI) on board SDO. In this work, we model the coronal magnetic field above multiple active regions with the help of a potential field and a NLFFF extrapolation codes in a full-disk using HMI data as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with full-disk coronal images from the Atmospheric Imaging Assembly (SDO/AIA) for both models. This study has found that the NLFFF model reconstructs the magnetic configuration better than the potential field model. We have concluded that much of trans-equatorial loops connecting the two solar hemispheres are current-free.

Extrapolation codes in Cartesian geometry for modelling the magnetic field in the corona do not take the curvature of the Sun's surface into account and can only be applied to relatively small areas, e.g., a single active region. We compare the analysis of the photospheric magnetic field and subsequent force-free modeling based on full-disk vector maps from Helioseismic and Magnetic Imager (HMI) on board solar dynamics observatory (SDO) and Vector Spectromagnetograph (VSM) of the Synoptic Optical Long-term Investigations of the Sun (SOLIS). We use Helioseismic and Magnetic Imager and Vector Spectromagnetograph photospheric magnetic field measurements to model the force-free coronal field above multiple solar active regions, assuming magnetic forces to dominate. We solve the nonlinear force-free field equations by minimizing a functional in spherical coordinates over a full disk excluding the poles. After searching for the optimum modeling parameters for the particular data sets, we compare the resulting nonlinear force-free model fields. We compare quantities like the total magnetic energy content and free magnetic energy, the longitudinal distribution of the magnetic pressure and surface electric current density using our spherical geometry extrapolation code. The magnetic field lines obtained from nonlinear force-free extrapolation based on Helioseismic and Magnetic Imager and Vector Spectromagnetograph data have good agreement. However, the nonlinear force-free extrapolation based on Helioseismic and Magnetic Imager data have more contents of total magnetic energy, free magnetic energy, the longitudinal distribution of the magnetic pressure and surface electric current density compared to the one from Vector Spectromagnetograph data.

We present the results of a simulation of a wind-driven non-linear gravity wave breaking on the surface of a white dwarf. The ``wind'' consists of H/He from an accreted envelope, and the simulation demonstrates that this breaking wave mechanism can produce a well-mixed layer of H/He with C/O from the white dwarf above the surface. Material from this mixed layer may then be transported throughout the accreted envelope by convection, which would enrich the C/O abundance of the envelope as is expected from observations of novae.

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.

We present direct numerical simulations of mixing at Rayleigh-Taylor unstable interfaces performed with the FLASH code, developed at the ASCI/Alliances Center for Astrophysical Thermonuclear Flashes at the University of Chicago. We present initial results of single-mode studies in two and three dimensions. Our results indicate that three-dimensional instabilities grow significantly faster than two-dimensional instabilities and that grid resolution can have a significant effect on instability growth rates. We also find that unphysical diffusive mixing occurs at the fluid interface, particularly in poorly resolved simulations.

We discuss some of the computational challenges encountered in simulating the evolution of clusters of galaxies. Eulerian adaptive mesh refinement (AMR) techniques can successfully address these challenges but are currently being used by only a few groups. We describe our publicly available AMR code, FLASH, which uses an object-oriented framework to manage its AMR library, physics modules, and automated verification. We outline the development of the FLASH framework to include collisionless particles, permitting it to be used for cluster simulation.

Adaptive techniques are crucial for successful numerical modeling of gravitational waves from astrophysical sources such as coalescing compact binaries, since the radiation typically has wavelengths much larger than the scale of the sources. We have carried out an important step toward this goal, the evolution of weak gravitational waves using adaptive mesh refinement in the Einstein equations. The 2-level adaptive simulation is compared with unigrid runs at coarse and fine resolution, and is shown to track closely the features of the fine grid run.

We present simulations of a model for the formation of a prominence condensation in a coronal loop. The key idea behind the model is that the spatial localization of loop heating near the chromosphere leads to a catastrophic cooling in the corona (Antiochos & Klimchuk 1991). Using a new adaptive grid code, we simulate the complete growth of a condensation, and find that after approx. 5,000 s it reaches a quasi-steady state. We show that the size and the growth time of the condensation are in good agreement with data, and discuss the implications of the model for coronal heating and SOHO/TRACE observations.

We present a method for integrating the cosmological hydrodynamical equations including a collisionless dark matter component. For modelling the baryonic matter component, we use the Piecewise Parabolic Method (PPM) which is a high-accuracy shock capturing technique. The dark matter component is modeled using gravitationally interacting particles whose evolution is determined using standard particle-in-cell techniques. We discuss details of the inclusion of gravity and expansion in the PPM code and give results of a number of tests of the code. This code has been developed for a massively parallel, SIMD supercomputer: the MasPar MP-2 parallel processor. We present details of the techniques we have used to implement the code for this architecture and discuss performance of the code on the MP-2. The code processes $5.0 \times 10^4$ grid zones per second and requires 53 seconds of machine time for a single timestep in a $128^3$ simulation.