Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Using Global Simulations to Relate the Three-Part Structure of Coronal Mass Ejections to in Situ Signatures

White-light observations of coronal mass ejections (CMEs) often show the classic "three-part" structure consisting of (1) a bright front; (2) a dark cavity; and (3) a bright, compact core. It has proven difficult to unambiguously associate these features with in situ measurements of interplanetary CMEs (ICMEs), in all but a few cases. In this study we use a global MHD model to simulate the eruption and evolution of a CME out to 0.25 AU, allowing us to continuously track these features from the Sun and through the solar wind. Our results support the generally held view that the interplanetary flux rope corresponds to the dark cavity. We find that the bright front merges with solar wind material swept up by the ICME. Thus, the sheath material found ahead of fast ejecta is in fact composed from both ambient solar wind material, as well the bright front. We also note that, in this simulation, the bright front is formed from the overlying streamer configuration from within which the CME erupted and is not itself coronal material swept up during the early phase of the eruption. The conclusions reached in this study are undoubtedly sensitive to the initial configuration and mechanism used to initiate the CME, and thus care should be taken when using them to interpret specific observations. On the other hand, they provide a unique, unbroken connection between remote solar and interplanetary observations. Ultimately, detailed comparisons between observations and simulation results may be able to constrain or even rule out some mechanisms of CME initiation.

Riley, Pete↗
NASA NTRS2008

Using Global Simulations to Relate the Three-Part Structure of Coronal Mass Ejections to In Situ Signatures

White-light observations of coronal mass ejections (CMEs) often show the classic "three-part" structure consisting of (1) a bright front; (2) a dark cavity; and (3) a bright, compact core. It has proven difficult to unambiguously associate these features with in situ measurements of interplanetary CMEs (ICMEs), in all but a few cases. In this study we use a global MHD model to simulate the eruption and evolution of a CME out to 0.25 AU, allowing us to continuously track these features from the Sun and through the solar wind. Our results support the generally held view that the interplanetary flux rope corresponds to the dark cavity. We find that the bright front merges with solar wind material swept up by the ICME. Thus, the sheath material found ahead of fast ejecta is in fact composed from both ambient solar wind material, as well the bright front. We also note that, in this simulation, the bright front is formed from the overlying streamer configuration from within which the CME erupted and is not itself coronal material swept up during the early phase of the eruption. The conclusions reached in this study are undoubtedly sensitive to the initial configuration and mechanism used to initiate the CME, and thus care should be taken when using them to interpret specific observations. On the other hand, they provide a unique, unbroken connection between remote solar and interplanetary observations. Ultimately, detailed comparisons between observations and simulation results may be able to constrain or even rule out some mechanisms of CME initiation

Riley, Pete↗
NASA NTRS2008

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗
NASA NTRS2011

SDO/AIA Observation of Kelvin-Helmholtz Instability in the Solar Corona

We present observations of the formation, propagation and decay of vortex-shaped features in coronal images from the Solar Dynamics Observatory (SDO) associated with an eruption starting at about 2:30UT on Apr 8, 2010. The series of vortices formed along the interface between an erupting (dimming) region and the surrounding corona. They ranged in size from several to ten arcseconds, and traveled along the interface at 6-14 km s-1. The features were clearly visible in six out of the seven different EUV wavebands of the Atmospheric Imaging Assembly (AIA). Based on the structure, formation, propagation and decay of these features, we identified these features as the first observations of the Kelvin- Helmholtz (KH) instability in the corona in EUV. The interpretation is supported by linear analysis and by MHD model of KH instability. We conclude that the instability is driven by the velocity shear between the erupting and closed magnetic field of the Coronal Mass Ejection (CME).

Ofman, L.↗
NASA NTRS2012

Comparison of Ejection Events in the Jet and Accretion Disc Outflows in 3C 111

We present a comparison of the parameters of accretion disc outflows and the jet of the broad-line radio galaxy 3C 111 on sub-pc scales. We make use of published X-ray observations of ultra-fast outflows (UFOs) and new 43 GHz VLBA images to track the jet knots ejection. We find that the superluminal jet coexists with the mildly relativistic outflows on sub-pc scales, possibly indicating a transverse stratification of a global flow. The two are roughly in pressure equilibrium, with the UFOs potentially providing additional support for the initial jet collimation. The UFOs are much more massive than the jet, but their kinetic power is probably about an order of magnitude lower, at least for the observations considered here. However, their momentum flux is equivalent and both of them are powerful enough to exert a concurrent feedback impact on the surrounding environment. A link between these components is naturally predicted in the context of MHD models for jet/outflow formation. However, given the high radiation throughput of AGNs, radiation pressure should also be taken into account. From the comparison with the long-term 2-10 keV RXTE light curve we find that the UFOs are preferentially detected during periods of increasing flux. We also find the possibility to place the UFOs within the known X-ray dips-jet ejection cycles, which has been shown to be a strong proof of the disc-jet connection, in analogue with stellar-mass black holes. However, given the limited number of observations presently available, these relations are only tentative and additional spectral monitoring is needed to test them conclusively.

Tombesi, F.↗
NASA NTRS2016

Density Variations in the Earth's Magnetospheric Cusps

Seven years of measurements from the Polar spacecraft are surveyed to monitor the variations of plasma density within the magnetospheric cusps. The spacecraft's orbital precession from 1998 through 2005 allows for coverage of both the northern and southern cusps from low altitude out to the magnetopause. In the mid- and high- altitude cusps, plasma density scales well with the solar wind density (n(sub cusp)∕n(sub sw) approximately 0.8). This trend is fairly steady for radial distances greater then 4 R(sub E). At low altitudes (r less than 4R(sub E)) the density increases with decreasing altitude and even exceeds the solar wind density due to contributions from the ionosphere. The density of high charge state oxygen (O(greater +2) also displays a positive trend with solar wind density within the cusp. A multifluid simulation with the Block-Adaptive-Tree Solar Wind Roe-Type Upwind Scheme MHD model was run to monitor the relative contributions of the ionosphere and solar wind plasma within the cusp. The simulation provides similar results to the statistical measurements from Polar and confirms the presence of ionospheric plasma at low altitudes.

Walsh, B. M.↗
NASA NTRS2017

A Method to Predict Magnetopause Expansion in Radial IMF Events by MHD Simulations

This paper presents a method for taking into account changes of solar wind parameters in the foreshock using global MHD simulations. We simulate four events with very distant subsolar magnetopause crossings that occurred during quasi-radial interplanetary magnetic field (IMF) intervals lasting from one to several hours. Using previous statistical results, we suggest that the density and velocity in the foreshock cavity decrease to approx. 60% and approx. 94% of the ambient solar wind values when the IMF cone angle falls below 50 deg. This diminishes the solar wind dynamic pressure to 53% and causes a corresponding magnetospheric expansion. We change the upstream solar wind parameters in a global MHD model to take these foreshock effects into account. We demonstrate that the modified model predicts magnetopause distances during radial IMF intervals close to those observed by THEMIS. The strong total pressure decrease in the data seems to be a local, rather than a global, phenomenon. Although the simulations with decreased solar wind pressure generally reproduce the observed total pressure in the magnetosheath well, the total pressure in the magnetosphere often agrees better with results for nonmodified boundary conditions. The last result reveals a limitation of our method: we changed the boundary conditions along the whole inflow boundary, although a more correct approach would be to vary parameters only in the foreshock. A model with the suggested global modification of the boundary conditions better predicts the location of part of the magnetopause behind the foreshock but may fail in predicting the rest of the magnetopause.

Samsonov, A. A.↗
NASA NTRS2018

Levitation Experiments on ISS: Science and Applications

Levitation experiments on the International Space Station (ISS) are ongoing. The European Space Agency (ESA) Materials Science Laboratory Electromagnetic Levitator (MSL-EML) has been in operation since 2015. US investigators are on several of the European Topical Teams and have been heavily involved with the experiments since 2009. The experiments include magneto-hydrodynamic (MHD) modeling of macro-convection, the effects of convection on microstructure, investigation of metallic glass formation, and much more. Recently NASA selected four proposals to the MaterialsLab NASA Research Announcement (NRA) for experiments on the Japan Aerospace eXploration Agency (JAXA) Electrostatic Levitation Furnace (ELF), which was launched to the ISS in 2015. It is estimated that the four US investigators will have flight experiments on ELF within the next few years. The experiments range from a novel method to measure interfacial tension to investigations of thermophysical properties of metals, metal oxides, and non-linear optical materials. Using the MSL-EML and ELF, US investigators are studying a wide range of NASA exploration-relevant materials. Levitation experiments enable exploration in many ways. For example, high quality thermophysical properties of high-temperature materials are critical to develop accurate models of casting, welding, and metal additive manufacturing, which could lead to more efficient and more reliable production of metallic parts for exploration, commercial, and industrial applications. High-quality thermophysical properties could also lead to the development of functional oxide glass and optical materials. In many cases, the accuracy of available property data is the limiting factor in the predictive capabilities of the models. Many thermophysical properties can be measured in a levitator on Earth, but with convective contamination. This contamination plays a significant role in the formation of the intermediate phases. In particular, nucleation and viscosity measurements demand quiescent conditions that can only be attained in microgravity-based levitation systems. A brief overview of the ongoing and planned levitation experiments on the ISS will be presented, followed by the exploration-relevance.

SanSoucie, Michael P.↗
NASA NTRS2022

Preferential Acceleration of Heavy Ions in a Spontaneously Fragmenting Flare Current Sheet

We study the ion acceleration in a mesoscale, spontaneously fragmenting flare current sheet (SFCS) characterized by the presence of a plasmoid cascade. The main subject of our investigation is to determine whether and how plasmoid cascades at intermediate scales in a fragmented current sheet of a solar flare can impact the (preferential) acceleration of specific ions. The time evolution of the SFCS is obtained from high-resolution 2.5D MHD simulations. The ion trajectories (in the background fields resulting from the MHD model), energies, and pitch angles are calculated using a relativistic test-particle code based on the half-acceleration–rotation–half-acceleration method. For light ions, the main acceleration effects of electromagnetic fields within the SFCS are analyzed using the guiding center approximation. We identify regions with the most-efficient ion acceleration within the SFCS, the accelerator efficiency, and spectra of the accelerated ions. The influence of the charge-to-mass ratio on ion behavior is also studied and resulting ion abundances are compared with observational data. The main ion acceleration takes place in the regions with a strong polarization term, which is part of the first-order Fermi acceleration. Because the term is mass dependent, heavier ions undergo preferential acceleration. The ion energy spectra, abundance-enhancement factors, and differential fluxes, obtained from the model, exhibit power-law profiles, in agreement with observed solar energetic particle events. Nonetheless, the obtained slopes for the abundance-enhancement factor do not exactly match the observed data. The computed slopes and profiles are not sensitive to changes in the initial plasma temperature.

David Kramoliš↗
NASA NTRS2024

A Feasibility Study of 4-D Tomography of Soft X-ray Magnetosheath Emissivities Using Multi-Spacecraft Measurements

Upcoming heliophysics missions utilize state-of-the-art wide field-of-view (FOV) imaging technology to measure and investigate the space plasma environment on a global scale. At Earth, remote sensing of soft X-ray emissions, which are generated via the charge exchange interaction between heavy solar wind ions and exospheric neutral atoms, is a promising means to investigate the global magnetosheath structure, its response to varying solar wind conditions, and the spatiotemporal properties of the dayside magnetic reconnection. Data analysis techniques such as optical tomography can provide additional structural and time-varying information from the observed target and thus enhance the mission’s scientific return. In this work, we simulate multiple and simultaneous observations of the dayside magnetosphere using soft X-ray imagers located at long-distance vantage points to reconstruct the time-dependent, three-dimensional (3-D) structure of the magnetosheath using a dynamic tomographic approach. The OpenGCCM MHD model is used to simulate the time-varying response of the magnetosheath to solar wind conditions and, subsequently, generate synthetic soft X-ray images from multiple spacecraft vantage points separated along a common orbit. A detailed analysis is then performed to identify the nominal set of spacecraft that produces the highest fidelity tomographic reconstruction of the magnetopause. This work aims to (i) demonstrate, for the first time, the use of dynamic tomography to retrieve the time-varying magnetosheath structure and (ii) identify a nominal mission design for multi-spacecraft configurations aiming for optical tomography.

3-D/4_d tomography↗
NASA NTRS2005

A Mechanism for the Loading-Unloading Substorm Cycle Missing in MHD Global Magnetospheric Simulation Models

Loading and consequent unloading of magnetic flux is an essential element of the substorm cycle in Earth's magnetotail. We are unaware of an available global MHD magnetospheric simulation model that includes a loading- unloading cycle in its behavior. Given the central role that MHD models presently play in the development of our understanding of magnetospheric dynamics, and given the present plans for the central role that these models will play in ongoing space weather prediction programs, it is clear that this failure must be corrected. A 2-dimensional numerical driven current-sheet model has been developed that incorporates an idealized current- driven instability with a resistive MHD system. Under steady loading, the model exhibits a global loading- unloading cycle. The specific mechanism for producing the loading-unloading cycle will be discussed. It will be shown that scale-free avalanching of electromagnetic energy through the model, from loading to unloading, is carried by repetitive bursts of localized reconnection. Each burst leads, somewhat later, to a field configuration that is capable of exciting a reconnection burst again. This process repeats itself in an intermittent manner while the total field energy in the system falls. At the end of an unloading interval, the total field energy is reduced to well below that necessary to initiate the next unloading event and, thus, a loading-unloading cycle results. It will be shown that, in this model, it is the topology of bursty localized reconnection that is responsible for the appearance of the loading-unloading cycle.

Klimas, A. J.↗
NASA NTRS1988

Shock-associated MHD waves - A model for interstellar density fluctuations

The possibility that the density fluctuations responsible for radio scintillations could be due to ion-beam-generated MHD waves near interstellar shock waves is discussed. This suggestion is inspired by spacecraft observations which reveal these phenomena near shocks in the solar system. The model quite naturally accounts for the scale on which these fluctuations occur; it is dictated by the wavelength of the unstable waves.

Spangler, Steven R.↗
NASA NTRS1990

A helically distorted MHD flux rope model

A flux rope model is proposed which has a variable degree of helical distortion from axisymmetry. The basis for this suggestion is a series of numerical and analytical investigations of magnetohydrodynamic states which result when an axial electric current is directed down on dc magnetic field. The helically distorted states involve a flow velocity and seem to be favored because of their lower rate of energy dissipation. Emphasis is on the magnetometer and particle energy analyzer traces that might be characteristic of such flux ropes. It is shown that even a fractionally small helical distortion may considerably alter the traces in minimum-variance coordinates. In short, what may be fairly common MHD processes can render a flux rope almost unrecognizable under standard diagnostics, even if the departures from axisymmetry are not great.

Theobald, Michael L.↗
NASA NTRS2024

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

In this study, we simulate the Solar Wind Charge Exchange (SWCX) soft X-ray emissions at dayside magnetosheath and cusps by using magnetohydrodynamic (MHD) and LAtmos TEst Particle (LaTeP) models. MHD models are unable to resolve the particle kinetic effects, such as the different behaviors of ions with different q/m, or distinguish the magnetospheric plasma from the solar wind plasma. We investigate these effects with the LaTeP model. As the LaTeP model does not self-compute magnetic and electric field, the magnetic and electric field data obtained from Open Geospace General Circulation Model (OpenGGCM) and Lagrangian version of the piecewise parabolic method (PPMLR) MHD model are used as the input to LaTeP model. The soft X-ray emissivity maps simulated from pure OpenGGCM and PPMLR MHD approaches and from LaTeP-OpenGGCM and LaTeP-PPMLR approaches are presented and compared. The results indicate that the LaTeP model can well resolve the kinetic effects and can be used to investigate the individual spectral characteristics. Therefore, the LaTeP model is a complementary approach for simulating the X-ray emissions near the dayside magnetopause. We also calculate the ratio of integrated OVII/OVIII line intensities, produced by charge exchange of O7+ ions and O8+ ions, respectively. We find a relatively higher ratio at the bow shock compared to the surrounding areas, suggesting that this ratio can be an effective parameter to identify the bow shock location.

Qiuyu Xu↗
NASA NTRS1993

A two-dimensional MHD global coronal model - Steady-state streamers

We describe a two-dimensional time-dependent, numerical, magnetohydrodynamic model for the determination of the physical properties of coronal streamers from the top of the transition zone (solar radius = 1) to 15 solar radii. Four examples are given: for dipole, quadrupole, and hexapole initial field topologies. The computed parameters are density, temperature, velocity, and magnetic field. In addition to the properties of the solutions, their accuracy is discussed. We use the model as the basis for a general discussion of the way boundary conditions are specified in this and similar simulations.

Wang, A.-H.↗
NASA NTRS2012

Pickup Protons: Comparisons using the Three-Dimensional MHD HHMS-PI model and Ulysses SWICS Measurements

We report the first comparisons of pickup proton simulation results with in situ measurements of pickup protons obtained by the SWICS instrument on Ulysses. Simulations were run using the three dimensional (3D) time-dependent Hybrid Heliospheric Modeling System with Pickup Protons (HHMS-PI). HHMS-PI is an MHD solar wind model, expanded to include the basic physics of pickup protons from neutral hydrogen that drifts into the heliosphere from the local interstellar medium. We use the same model and input data developed by Detman et al. (2011) to now investigate the pickup protons. The simulated interval of 82 days in 2003 2004, includes both quiet solar wind (SW) and also the October November 2003 solar events (the Halloween 2003 solar storms). The HHMS-PI pickup proton simulations generally agree with the SWICS measurements and the HHMS-PI simulated solar wind generally agrees with SWOOPS (also on Ulysses) measurements. Many specific features in the observations are well represented by the model. We simulated twenty specific solar events associated with the Halloween 2003 storm. We give the specific values of the solar input parameters for the HHMS-PI simulations that provide the best combined agreement in the times of arrival of the solar-generated shocks at both ACE and Ulysses. We show graphical comparisons of simulated and observed parameters, and we give quantitative measures of the agreement of simulated with observed parameters. We suggest that some of the variations in the pickup proton density during the Halloween 2003 solar events may be attributed to depletion of the inflowing local interstellar medium (LISM) neutral hydrogen (H) caused by its increased conversion to pickup protons in the immediately preceding shock.

Intriligator, Devrie S.↗