Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

CME Flux Rope and Shock Identifications and Locations: Comparison of White Light Data, Graduated Cylindrical Shell Model, and MHD Simulations

Coronal mass ejections (CMEs) are major transient phenomena in the solar corona that are observed with ground-based and spacecraft-based coronagraphs in white light or with in situ measurements by spacecraft. CMEs transport mass and momentum and often drive shocks. In order to derive the CME and shock trajectories with high precision, we apply the graduated cylindrical shell (GCS) model to fit a flux rope to the CME directed toward STEREO A after about 19:00 UT on 29 November 2013 and check the quality of the heliocentric distance-time evaluations by carrying out a three-dimensional magnetohydrodynamic (MHD) simulation of the same CME with the Block Adaptive Tree Solar-Wind Roe Upwind Scheme (BATS-R-US) code. Heliocentric distances of the CME and shock leading edges are determined from the simulated white light images and magnetic field strength data. We find very good agreement between the predicted and observed heliocentric distances, showing that the GCS model and the BATS-R-US simulation approach work very well and are consistent. In order to assess the validity of CME and shock identification criteria in coronagraph images, we also compute synthetic white light images of the CME and shock. We find that the outer edge of a cloud-like illuminated area in the observed and predicted images in fact coincides with the leading edge of the CME flux rope and that the outer edge of a faint illuminated band in front of the CME leading edge coincides with the CME-driven shock front.

Schmidt, J. M.↗
NASA NTRS2017

The Mars Crustal Magnetic Field Control of Plasma Boundary Locations and Atmospheric Loss: MHD Prediction and Comparison with MAVEN

We present results from a global Mars time-dependent MHD simulation under constant solar wind and solar radiation impact considering inherent magnetic field variations due to continuous planetary rotation. We calculate the 3-D shapes and locations of the bow shock (BS) and the induced magnetospheric boundary (IMB) and then examine their dynamic changes with time. We develop a physics-based, empirical algorithm to effectively summarize the multidimensional crustal field distribution. It is found that by organizing the model results using this new approach, the Mars crustal field shows a clear, significant influence on both the IMB and the BS. Specifically, quantitative relationships have been established between the field distribution and the mean boundary distances and the cross-section areas in the terminator plane for both of the boundaries. The model-predicted relationships are further verified by the observations from the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. Our analysis shows that the boundaries are collectively affected by the global crustal field distribution, which, however, cannot be simply parameterized by a local parameter like the widely used subsolar longitude. Our calculations show that the variability of the intrinsic crustal field distribution in Mars-centered Solar Orbital itself may account for approx.60% of the variation in total atmospheric loss, when external drivers are static. It is found that the crustal field has not only a shielding effect for atmospheric loss but also an escape-fostering effect by positively affecting the transterminator ion flow cross-section area.

Global↗
NASA NTRSDate not supplied

Physical Properties of the Solar Atmosphere Derived from Comparison of Spectro-Polarimetric SDO/HMI Observables with 3D Radiative MHD Simulations

In this study, we compare the SDO/HMI line-of-sight observables (magnetic field, velocity, continuum intensity, and line depth) with the related physical properties for several dynamo simulation runs performed using the "StellarBox" 3D Radiative MHD code. The modeling of the Fe I 6173 A Stokes profiles is performed using the SPINOR radiative transfer code in the LTE approximation. The reproduced SDO/HMI line-of-sight pipeline is applied to the modeled spectra, and the observables are synthesized with high (numerical) and SDO/HMI (instrumental) resolutions. Correlations between the observables and the physical properties at various heights in the atmosphere are studied for a set of view angles (0, 30, 45, 60, 70, and 80 degrees away from the solar disk center). It is found that the SDO/HMI magnetic field and velocity measurements are unambiguously correlated with physical parameters at certain heights of the solar atmosphere. These heights increase from about 100 km above the photosphere for the disk center case to 300-600 km above the photosphere for the 80-degree case. The heights are found to be slightly lower in regions where stronger magnetic fields are found. The comparison of the photospheric magnetic flux and integrated continuum intensity derived from the SDO/HMI observables and high-resolution observations and spectra is discussed. The results of our study improve physics-based interpretations of the SDO/HMI observables and provide a better understanding of the physical properties of the solar atmosphere.

SMD↗
NASA NTRS2023

Investigating Potential Causes for the Prediction of Spurious Magnetopause Crossings at Geosynchronous Orbit in MHD Simulations

During intense geomagnetic storms, the magnetopause can move in as far as geosynchronous orbit, leaving the satellites in that orbit out in the magnetosheath. Spacecraft operators turn to numerical models to predict the response of the magnetopause to solar wind conditions, but the predictions of the models are not always accurate. This study investigates four storms with a magnetopause crossing by at least one GOES satellite, using four magnetohydrodynamic models at NASA's Community Coordinated Modeling Center to simulate the events, and analyzes the results to investigate the reasons for errors in the predictions. Two main reasons can explain most of the erroneous predictions. First, the solar wind input to the simulations often contains features measured near the L1 point that did not eventually arrive at Earth; incorrect predictions during such periods are due to the solar wind input rather than to the models themselves. Second, while the models do well when the primary driver of magnetopause motion is a variation in the solar wind density, they tend to overpredict or underpredict the integrated Birkeland currents and their effects during times of strong negative interplanetary magnetic field (IMF) Bz, leading to poorer prediction capability. Coupling the MHD codes to a ring current model, when such a coupling is available, generally will improve the predictions but will not always entirely correct them. More work is needed to fully characterize the response of each code under strong southward IMF conditions as it relates to prediction of magnetopause location.

magnetopause↗
DOE OSTI2015

Particle acceleration and plasma energization in substorms: MHD and test particle studies

The author organizes his slide presentation under the following topics: background, MHD simulation, orbit integration, typical orbits, spatial and temporal features, acceleration mechanisms, source locations, and source energies. Field-­aligned energetic particle fluxes are shown for 45-keV electrons and 80-keV protons. It is concluded that the onset from local thin current sheet is electron tearing. Acceleration is mainly from field collapse, governed by E y = -v x XB z : importance of localization; betatron acceleration (similar if nonadiabatic); 1st order Fermi, type B (or A; current sheet acceleration). There are two source regions (of comparable importance in magnetotail): - flanks, inner tail - drift entry - early, higher energy - outer plasma sheet - reconnection entry - later, lower energy. Both thermal and suprathermal sources are important, with limited energy range for acceleration

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Simulating the Heliosphere with Kinetic Hydrogen and Dynamic MHD Source Terms

The interaction between the ionized plasma of the solar wind (SW) emanating from the sun and the partially ionized plasma of the local interstellar medium (LISM) creates the heliosphere. The heliospheric interface is characterized by the tangential discontinuity known as the heliopause that separates the SW and LISM plasmas, and a termination shock on the SW side along with a possible bow shock on the LISM side. Neutral Hydrogen of interstellar origin plays a critical role in shaping the heliospheric interface, since it freely traverses the heliopause. Charge-exchange between H-atoms and plasma protons couples the ions and neutrals, but the mean free paths are large, resulting in non-equilibrated energetic ion and neutral components. In our model, source terms for the MHD equations are generated using a kinetic approach for hydrogen, and the key computational challenge is to resolve these sources with sufficient statistics. For steady-state simulations, statistics can accumulate over arbitrarily long time intervals. In this paper we discuss an approach for improving the statistics in time-dependent calculations, and present results from simulations of the heliosphere where the SW conditions at the inner boundary of the computation vary according to an idealized solar cycle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Nonlinear 3D M3D-C1 simulations of tokamak plasmas crossing a MHD linear stability boundary

The goal of the present work is to better understand and develop a predictive capability for when approaching and crossing a MHD linear instability boundary leads to a thermal quench and subsequent disruption (hard limit), and when it just leads to increased transport or small amplitude oscillations(soft). Understanding the difference between hard and soft limits is crucial for effective disruption prediction and avoidance. We present several examples of both hard and soft beta limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
DOE OSTI2011

Ideal MHD Stability of ITER Steady State Scenarios with ITBs

One of ITER goals is to demonstrate feasibility of continuous operations using non-inductive current drive. Two main candidates have been identified for advanced operations: the long duration, high neutron fluency hybrid scenario and the steady state scenario, both operating at a plasma current lower than the reference ELMy scenario [1][2] to minimize the required current drive. The steady state scenario targets plasmas with current 7-10 MA in the flat-top, 50% of which will be provided by the self-generated, pressure-driven bootstrap current. It has been estimated that, in order to obtain a fusion gain Q > 5 at a current of 9 MA, it should be ΒN > 2.5 and H > 1.5 [3]. This implies the presence of an Internal Transport Barrier (ITB). This work discusses how the stability of steady state scenarios with ITBs is affected by the external heating sources and by perturbations of the equilibrium profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A BOUT++ extension for full annular tokamak edge MHD and turbulence simulations

For tokamak edge plasma simulation, a plasma simulation framework BOUT++ employs a dual coordinate system to simulate moderate-n and high-n plasma instability with reasonable computational cost, where n is the toroidal mode number. This coordinate system however limits the computational domain to the toroidal wedge (full torus divided into N parts in the toroidal direction) for computational efficiency and the use of flute-ordering approximation in the field solver calculating the flow potential from the vorticity which may not be valid for low-n modes. Improving numerical treatment of low-n modes is however indispensable to address simulations of low-n current-driven edge localized mode (ELM), ELM control by resonant magnetic perturbations (RMPs), edge turbulence with RMPs and so on. In this work, BOUT++ is extended to simulate the interplay between $n=0$, low-n and high-n plasma components in a full annular tokamak edge domain through hybrid modeling of the flow potential and the vorticity. Low-n modes of flow potential are calculated in an orthogonal flux surface coordinate and high-n modes in the dual coordinate system separately in Fourier space. Finally, the proposed scheme can capture an interplay between $n=1$ global modes and high-n turbulence during pedestal collapse in a full annular torus domain with a circular cross section.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Importance of Ambipolar Electric Field in Driving Ion Loss From Mars: Results From a Multifluid MHD Model With the Electron Pressure Equation Included

The multi-fluid (MF) magnetohydrodynamic model of Mars is improved by solving an additional electron pressure equation. Through the electron pressure equation, the electron temperature is calculated based on the effects from various electron-related heating and cooling processes (e.g., photo-electron heating, electron-neutral collision, and electron-ion collision), and thus, the improved model can calculate the electron temperature and the electron pressure force terms self-consistently. Model results of a typical case using the MF with electron pressure equation included model are compared in detail to identical cases using the MF and multi-species models to identify the effect of the improved physics. Here, we find that when the electron pressure equation is included, the general interaction patterns are similar to those with no electron pressure equation. However, the MF with electron pressure equation included model predicts that the electron temperature is much larger than the ion temperature in the ionosphere, consistent with both Viking and Mars Atmosphere and Volatile EvolutioN (MAVEN) observations. Using our numerical model, we also examined in detail the relative importance of different forces in the plasma interaction region. All three models are also applied to a MAVEN event study using identical input conditions; overall, the improved model matches best with MAVEN observations. All of the simulation cases are examined in terms of the total ion loss, and the results show that the inclusion of the electron pressure equation increases the escape rates by 50–110% in total mass, depending on solar condition and strong crustal field orientation, clearly demonstrating the importance of the ambipolar electric field in facilitating ion escape.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The role of MHD in 3D aspects of massive gas injection

Simulations of massive gas injection for disruption mitigation in DIII-D are carried out to compare the toroidal peaking of radiated power for the cases of one and two gas jets. The radiation toroidal peaking factor (TPF) results from a combination of the distribution of impurities and the distribution of heat flux associated with the n = 1 mode. When we ignored the effects of strong uni-directional neutral beam injection and rotation present in the experiment, the injected impurities were found to spread helically along field lines preferentially toward the high-field-side, which is explained in terms of a nozzle equation. Thus when considering the plasma rest frame, reversing the current direction also reverses the toroidal direction of impurity spreading. During the pre-thermal quench phase of the disruption, the toroidal peaking of radiated power is reduced in a straightforward manner by increasing from one to two gas jets. But, during the thermal quench phase, reduction in the TPF is achieved only for a particular arrangement of the two gas valves with respect to the field line pitch. Particularly, the relationship between the two valve locations and the 1/1 mode phase is critical, where gas valve spacing that is coherent with 1/1 symmetry effectively reduces TPF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Role of electric fields in the MHD evolution of the kink instability

Here, the discovery of electrostatic fields playing a crucial role in establishing plasma motion in the flux conversion and dynamo processes in reversed field pinches is revisited. In order to further elucidate the role of the electrostatic fields, a flux rope configuration susceptible to the kink instability is numerically studied with anMHDcode. Simulated nonlinear evolution of the kink instability is found to confirm the crucial role of the electrostatic fields. Anew insight is gained on the special function of the electrostatic fields: they lead the plasma towards the reconnection site at the mode resonant surface. Without this step the plasma column could not relax to its nonlinear state, since no other agent is present to perform this role. While the inductive field generated directly by the kink instability is the dominant flow driver, the electrostatic field is found to allow the motion in the vicinity of the reconnection region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

MHD simulations of cold bubble formation from 2/1 tearing mode during massive gas injection in a tokamak

Massive gas injection (MGI) experiments have been carried out in many tokamaks to study disruption dynamics and mitigation schemes. Two events often observed in those experiments are the excitation of the m = 2, n = 1 magnetohydrodynamic mode, and the formation of cold bubble structure in the temperature distribution before the thermal quench (TQ). Here m is the poloidal mode number, n the toroidal mode number. The physics mechanisms underlying those phenomena, however, have not been entirely clear. In this work, our recent NIMROD simulations of the MGI process in a tokamak have reproduced the main features of both events, which has allowed us to examine and establish the causal relation between them. In these simulations, the 3/1 and 2/1 islands are found to form successively after the arrival of impurity ion cold front at the corresponding q = 3 and q = 2 rational surfaces. At the interface between impurity and plasma, a local thin current sheet forms due to an enhanced local pressure gradient and moves inward following the gas cold front, this may contribute to the formation of a dominant 2/1 mode. Following the growth of the 2/1 tearing mode, the impurity penetration into the core region inside the q = 2 surface gives rise to the formation of the cold bubble temperature structure and initiates the final TQ. Here, a subdominant 1/1 mode developed earlier near the q = 1 surface alone does not cause such a cold bubble formation, however, the exact manner of the preceding impurity penetration depends on the nature of the 1/1 mode: kink-tearing or quasi-interchange.

2/1 tearing mode↗