Three-dimensional unstructured grid generation via incremental insertion and local optimization

Algorithms for the generation of 3D unstructured surface and volume grids are discussed. These algorithms are based on incremental insertion and local optimization. The present algorithms are very general and permit local grid optimization based on various measures of grid quality. This is very important; unlike the 2D Delaunay triangulation, the 3D Delaunay triangulation appears not to have a lexicographic characterization of angularity. (The Delaunay triangulation is known to minimize that maximum containment sphere, but unfortunately this is not true lexicographically). Consequently, Delaunay triangulations in three-space can result in poorly shaped tetrahedral elements. Using the present algorithms, 3D meshes can be constructed which optimize a certain angle measure, albeit locally. We also discuss the combinatorial aspects of the algorithm as well as implementational details

Comparison of Birkeland current observations during two magnetic cloud events with MHD simulations

Low altitude field-aligned current densities ob- tained from global magnetospheric simulations are compared with two-dimensional distributions of Birkeland currents at the topside ionosphere derived from magnetic field observa- tions by the constellation of Iridium satellites. We present the analysis of two magnetic cloud events, 17–19 August 2003 and 19–21 March 2001, where the interplanetary magnetic field (IMF) rotates slowly (∼10◦/h) to avoid time-aliasing in the magnetic perturbations used to calculate the Birkeland currents. In the August 2003 event the IMF rotates from southward to northward while maintaining a negative IMF By during much of the interval. During the March 2001 event the IMF direction varies from dawnward to southward to duskward. We find that the distributions of the Birkeland current densities in the simulations agree qualitatively with the observations for northward IMF. For southward IMF, the dayside Region-1 currents are reproduced in the simu- ◦ the ionospheric grids in the simulations and the observations is shown to have only secondary effect on the magnitudes of the Birkeland currents. The electric potentials in the simu- lation for southward IMF periods are twice as large as those obtained from measurements of the plasma drift velocities by DMSP, implying that the reconnection rates in the simulation are too large. Keywords. Ionosphere (Electric fields and currents; Ionosphere-magnetosphere interactions; Modeling and forecasting) 1 Introduction Global magnetohydrodynamic (MHD) models are the most comprehensive numerical tool for studying the coupling of energy and momentum of the solar wind into the Earth’s magnetosphere and ionosphere. A particular advantage of global MHD simulations is the ability to provide continu- ous temporal and spatial coverage of key physical parame- ters over the entire simulation volume. For this reason, MHD simulations have become one of the principal tools for study- ing space weather events such as the interaction of the Earth’s magnetosphere with coronal mass ejections (CMEs) (Ridley et al., 2002) as well as magnetic storms (Slinker et al., 1998; Goodrich et al., 1998) and substorms (Lyon et al., 1998; Lopez et al., 1998; Wiltberger et al., 2000). Since the simula- tion results are frequently used to interpret physical processes in the magnetosphere–ionosphere system, assessing their ac- curacy by comparison with observations is an important task. A number of such studies have been carried out in the past us- ing space-based (Frank et al., 1995; Raeder et al., 1997) and ground-based observations (Ridley et al., 2001), or a com- bination thereof (Fedder et al., 1998; Slinker et al., 1999). However, interpreting the discrepancies between model and observations is not straightforward because the observational lation, but appear on average 5 served location, while the nightside Region-1 currents and the Region-2 currents are largely under-represented. Com- parison of the observed and simulated Birkeland current dis- tributions, which are intimately related to the plasma drifts at the ionosphere, shows that the ionospheric convection pat- tern in the MHD model and its dependence on the IMF ori- entation is essentially correct. The Birkeland total currents in the simulations are about a factor of 2 larger than observed during southward IMF. For Bz\u3e0 the disparity in the total current is reduced and the simulations for purely northward IMF agree with the observations to within 10%. The dispar- ities in the magnitudes of the Birkeland currents between the observations and the simulation results are a combined effect of the simulation overestimating the ionospheric electric field and of the Iridium fits underestimating the magnetic pertur- bations

A novel metric for coronal MHD models

[1] In the interest of quantitatively assessing the capabilities of coronal MHD models, we have developed a metric that compares the structures of the white light corona observed with SOHO LASCO C2 to model predictions. The MAS model is compared to C2 observations from two Carrington rotations during solar cycle 23, CR1913 and CR1984, which were near the minimum and maximum of solar activity, respectively, for three radial heights, 2.5 R⊙, 3.0 R⊙, and 4.5 R⊙. In addition to simulated polarization brightness images, we create a synthetic image based on the field topology along the line of sight in the model. This open-closed brightness is also compared to LASCO C2 after renormalization. In general, the model\u27s magnetic structure is a closer match to observed coronal structures than the model\u27s density structure. This is expected from the simplified energy equations used in current global corona MHD models

Magnetohydrodynamic Modeling of Three Van Allen Probes Storms in 2012 and 2013

Coronal mass ejection (CME)-shock compression of the dayside magnetopause has been observed to cause both prompt enhancement of radiation belt electron flux due to inward radial transport of electrons conserving their first adiabatic invariant and prompt losses which at times entirely eliminate the outer zone. Recent numerical studies suggest that enhanced ultra-low frequency (ULF) wave activity is necessary to explain electron losses deeper inside the magnetosphere than magnetopause incursion following CME-shock arrival. A combination of radial transport and magnetopause shadowing can account for losses observed at radial distances into L=4.5, well within the computed magnetopause location. We compare ULF wave power from the Electric Field and Waves (EFW) electric field instrument on the Van Allen Probes for the 8 October 2013 storm with ULF wave power simulated using the Lyon–Fedder–Mobarry (LFM) global magnetohydrodynamic (MHD) magnetospheric simulation code coupled to the Rice Convection Model (RCM). Two other storms with strong magnetopause compression, 8–9 October 2012 and 17–18 March 2013, are also examined. We show that the global MHD model captures the azimuthal magnetosonic impulse propagation speed and amplitude observed by the Van Allen Probes which is responsible for prompt acceleration at MeV energies reported for the 8 October 2013 storm. The simulation also captures the ULF wave power in the azimuthal component of the electric field, responsible for acceleration and radial transport of electrons, at frequencies comparable to the electron drift period. This electric field impulse has been shown to explain observations in related studies (Foster et al., 2015) of electron acceleration and drift phase bunching by the Energetic Particle, Composition, and Thermal Plasma Suite (ECT) instrument on the Van Allen Probes

Effects of electrojet turbulence on a magnetosphere-ionosphere simulation of a geomagnetic storm

Ionospheric conductance plays an important role in regulating the response of the magnetosphere‐ionosphere system to solar wind driving. Typically, models of magnetosphere‐ionosphere coupling include changes to ionospheric conductance driven by extreme ultraviolet ionization and electron precipitation. This paper shows that effects driven by the Farley‐Buneman instability can also create significant enhancements in the ionospheric conductance, with substantial impacts on geospace. We have implemented a method of including electrojet turbulence (ET) effects into the ionospheric conductance model utilized within geospace simulations. Our particular implementation is tested with simulations of the Lyon‐Fedder‐Mobarry global magnetosphere model coupled with the Rice Convection Model of the inner magnetosphere. We examine the impact of including ET‐modified conductances in a case study of the geomagnetic storm of 17 March 2013. Simulations with ET show a 13% reduction in the cross polar cap potential at the beginning of the storm and up to 20% increases in the Pedersen and Hall conductance. These simulation results show better agreement with Defense Meteorological Satellite Program observations, including capturing features of subauroral polarization streams. The field‐aligned current (FAC) patterns show little differences during the peak of storm and agree well with Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE) reconstructions. Typically, the simulated FAC densities are stronger and at slightly higher latitudes than shown by AMPERE. The inner magnetospheric pressures derived from Tsyganenko‐Sitnov empirical magnetic field model show that the inclusion of the ET effects increases the peak pressure and brings the results into better agreement with the empirical model.This material is based upon work supported by NASA grants NNX14AI13G, NNX13AF92G, and NNX16AB80G. The National Center for Atmospheric Research is sponsored by the National Science Foundation. This work used the XSEDE and TACC computational facilities, supported by National Science Foundation grant ACI-1053575. We would like to acknowledge high-performance computing support from Yellowstone (ark:/85065/d7wd3xhc) provided by NCAR's Computational and Information Systems Laboratory, sponsored by the National Science Foundation. We thank the AMPERE team and the AMPERE Science Center for providing the Iridium derived data products. All model output, simulation codes, and analysis routines are being preserved on the NCAR High-Performance Storage System and will be made available upon written request to the lead author of this publication. (NNX14AI13G - NASA; NNX13AF92G - NASA; NNX16AB80G - NASA; National Science Foundation; ACI-1053575 - National Science Foundation

An event study to provide validation of TING and CMIT geomagnetic middle-latitude electron densities at the F2 peak

[1] The coupled thermosphere-ionosphere magnetosphere (CMIT) model and the Thermosphere Ionosphere Nested Grid (TING) model have been run to simulate the 15 May 1997 interplanetary coronal mass ejection\u27s (ICME) effects on the Earth\u27s ionosphere and thermosphere. Comparisons were made between these model runs, the IRI-2007 model, and geomagnetic middle-latitude ionosonde data (NmF2) from the World Data Center to determine how well the models simulated the event and to understand the causes of model-data disagreement. The following conclusions were drawn from this study: (1) skill scores were more often negative than positive on average; (2) the best and the worst skill scores occurred on the recovery day; (3) the line plots comparing models to data look better than the skill scores might suggest; (4) skill scores are significantly affected by timing issues and large, short-duration variability; (5) skill scores give an indication of the relative ability of one model relative to another, rather than an absolute statement of model accuracy; (6) the models capture negative storm effects better than they capture positive storm effects; (7) the TING model captured many short duration features seen in the data at high middle latitude stations that result from changes in the size of the auroral oval; (8) CMIT overestimates the energy driving changes in NmF2, whereas TING provides approximately the correct energy input as a result of the saturation effects on potential that are included in TING; and (9) both TING and CMIT electron densities decreased too rapidly after sunset

Predicting magnetopause crossings at geosynchronous orbit during the Halloween storms

[1] In late October and early November of 2003, the Sun unleashed a powerful series of events known as the Halloween storms. The coronal mass ejections launched by the Sun produced several severe compressions of the magnetosphere that moved the magnetopause inside of geosynchronous orbit. Such events are of interest to satellite operators, and the ability to predict magnetopause crossings along a given orbit is an important space weather capability. In this paper we compare geosynchronous observations of magnetopause crossings during the Halloween storms to crossings determined from the Lyon-Fedder-Mobarry global magnetohydrodynamic simulation of the magnetosphere as well to predictions of several empirical models of the magnetopause position. We calculate basic statistical information about the predictions as well as several standard skill scores. We find that the current Lyon-Fedder-Mobarry simulation of the storm provides a slightly better prediction of the magnetopause position than the empirical models we examined for the extreme conditions present in this study. While this is not surprising, given that conditions during the Halloween storms were well outside the parameter space of the empirical models, it does point out the need for physics-based models that can predict the effects of the most extreme events that are of significant interest to users of space weather forecasts

Posteruptive phenomena in coronal mass ejections and substorms: Indicators of a universal process?

[1] We examine phenomena associated with eruptions in the two different regimes of the solar corona and the terrestrial magnetosphere. We find striking similarities between the speeds of shrinking magnetic field lines in the corona and dipolarization fronts traversing the magnetosphere. We also examine the similarities between supra-arcade downflows observed during solar flares and bursty bulk flows seen in the magnetotail and find that these phenomena have remarkably similar speeds, velocity profiles, and size scales. Thus we show manifest similarities in the magnetic reconfiguration in response to the ejection of coronal mass ejections in the corona and the ejection of plasmoids in the magnetotail. The subsequent return of loops to a quasi-potential state in the corona and field dipolarization in the magnetotail are physical analogs and trigger similar phenomena such as downflows, which provides key insights into the underlying drivers of the plasma dynamics

Magnetic Flux Circulation During Dawn-Dusk Oriented Interplanetary Magnetic Field

Magnetic flux circulation is a primary mode of energy transfer from the solar wind into the ionosphere and inner magnetosphere. For southward interplanetary magnetic field (IMF), magnetic flux circulation is described by the Dungey cycle (dayside merging, night side reconnection, and magnetospheric convection), and both the ionosphere and inner magnetosphere receive energy. For dawn-dusk oriented IMF, magnetic flux circulation is not well understood, and the inner magnetosphere does not receive energy. Several models have been suggested for possible reconnection patterns; the general pattern is: dayside merging; reconnection on the dayside or along the dawn/dusk regions; and, return flow on dayside only. These models are consistent with the lack of energy in the inner magnetosphere. We will present evidence that the Dungey cycle does not explain the energy transfer during dawn-dusk oriented IMF. We will also present evidence of how magnetic flux does circulate during dawn-dusk oriented IMF, specifically how the magnetic flux reconnects and circulates back