Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma global models”

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 127 records · Page 7

Magnetic field and current structures in the magnetosphere of Uranus

Voyager 2 magnetic field and plasma data are compared with theoretical model calculations for the magnetosphere of Uranus to derive a global picture from the limited set of measurements. The results suggest that Voyager 2 entered the Uranian magnetosphere during a relaxation phase which followed a compression of the entire magnetosphere. The plasma beta values in the Uranian neutral sheet were found to be smallar by a factor of 3-8 compared to corresponding beta values in the average terrestrial neutral sheet. The excellent agreement found between observed and calculated magnetic tail lobe field strengths indicates that the Uranian magnetosphere reaches the state of quasi-static (i.e., slowly time-dependent MHD equilibrium).

Voigt, G.-H.↗

Strong IMF By-Related Plasma Convection in the Ionosphere and Cusp Field-Aligned Currents Under Northward IMF Conditions

We present in this paper an investigation of IMF-By related plasma convection and cusp field-aligned currents using FAST data and AMIE model during a prolonged interval with large positive IMF By and northward Bz conditions (By/Bz much greater than 1). Using the FAST single trajectory observations to validate the global convection patterns at key times and key locations, we have demonstrated that the AMIE procedure provides a reasonably good description of plasma circulations in the ionosphere during this interval. Our results show that the plasma convection in the ionosphere is consistent with the anti-parallel merging model. When the IMF has a strongly positive By component under northward conditions, we find that the global plasma convection forms two cells oriented nearly along the Sun-earth line in the ionosphere. In the northern hemisphere, the dayside cell has clockwise convection mainly circulating within the polar cap on open field lines. A second cell with counterclockwise convection is located in the nightside circulating across the polar cap boundary, The observed two-cell convection pattern appears to be driven by the reconnection along the anti-parallel merging lines poleward of the cusp extending toward the dusk side when IMF By/Bz much greater than 1. The magnetic tension force on the newly reconnected field lines drives the plasma to move from dusk to dawn in the polar cusp region near the polar cap boundary. The field-aligned currents in the cusp region flow downward into the ionosphere. The return field-aligned currents extend into the polar cap in the center of the dayside convection cell. The field-aligned currents are closed through the Peterson currents in the ionosphere, which flow poleward from the polar cap boundary along the electric field direction.

Le, G.↗

Dynamical Evolution of the Inner Heliosphere Approaching Solar Activity Maximum: Interpreting Ulysses Observations Using a Global MHD Model

In this study we describe a series of MHD simulations covering the time period from 12 January 1999 to 19 September 2001 (Carrington Rotation 1945 to 1980). This interval coincided with: (1) the Sun s approach toward solar maximum; and (2) Ulysses second descent to the southern polar regions, rapid latitude scan, and arrival into the northern polar regions. We focus on the evolution of several key parameters during this time, including the photospheric magnetic field, the computed coronal hole boundaries, the computed velocity profile near the Sun, and the plasma and magnetic field parameters at the location of Ulysses. The model results provide a global context for interpreting the often complex in situ measurements. We also present a heuristic explanation of stream dynamics to describe the morphology of interaction regions at solar maximum and contrast it with the picture that resulted from Ulysses first orbit, which occurred during more quiescent solar conditions. The simulation results described here are available at: http://sun.saic.com.

Riley, Pete↗

Electric field and turbulence in global Braginskii simulations across the ASDEX Upgrade edge and scrape-off layer

Turbulence simulations in diverted geometry across the edge and scrape-off layer (SOL) of ASDEX Upgrade are performed with the GRILLIX code (Stegmeir et al 2019 Phys. Plasmas 26 052517). The underlying global (full-f) drift-reduced Braginskii model allows to concurrently study the self-consistent dynamics of the turbulence and the background as well as the evolution of toroidal and zonal flows. Different contributions to the radial electric field are identified. The dominant contribution on closed flux surfaces comes from the ion pressure gradient, due to the diamagnetic drift in the curved magnetic field. Large deviations can be induced, in particular, by the polarization particle flux, leading to zonal flows. The latter are driven by small-scale eddies, but do not exhibit much impact on the overall transport which is driven by ballooning modes at larger scales. Ion viscosity is found to be important in damping poloidal rotation through adjusting of the parallel velocity profile, but not via direct vorticity damping. The zonal flow drive peaks at the separatrix, where a strong shear layer forms due to the sheath-induced counter-propagating SOL flow, allowing for the formation of a transport barrier. The temperature profile across the separatrix is determined by the competition between cross-field transport and outflow in the SOL, the latter being largely controlled by the parallel heat conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Role of Self-Organized Criticality in the Substorm Phenomenon and its Relation to Localized Reconnection in the Plasma Sheet

Recent observations of the magnetotail plasma sheet have shown it to be a dynamic and turbulent region. Research has found strong turbulence in the plasma sheet at approximately 20 Earth's Radius tailward of Earth; the turbulence is observed at all activity levels. The existence of strong turbulence in the plasma sheet in the region associated with substorm onset might be thought difficult to reconcile with the coherence and repeatability of the substorm cycle. We review a variety of evidence that strongly suggests the magnetotail is driven, through magnetic flux transfer, into a state of "self-organized criticality" (SOC). It is an important property of physical systems that evolve into SOC that they self-organize into a unique, global dynamic state. This global state is inevitable, and repeatable. In this state, however, small-spatiotemporal-scale system fluctuations are unpredictable and can be only described statistically. This is the basis, we think, for the global coherence and repeatability of the substorm phenomenon in the turbulent plasma sheet. At, or near, substorm onset the plasma sheet can be described by a global SOC state containing significant small scale turbulence. In several recent studies, "sandpile" models were driven into SOC and then shown to reproduce various measures of substorm activity. We discuss the plasma physical foundation of these sandpile models. The evolution of simple continuum plasma sheet models into SOC-like states of many small reconnection events in the turbulent plasma sheet under the will be demonstrated. We view the substorm phenomenon as an avalanche assumption that the plasma sheet is in a SOC state.

Klimas, A. J.↗

Ion-neutral momentum coupling near discrete high-latitude ionospheric features

A two-dimensional numerical model is developed to study the momentum coupling between the ionosphere and neutral atmosphere in the vicinity of discrete high-latitude features, such as convection channels and plasma density troughs. Based on generalized magnetohydrodynamic equations the model takes account of global pressure gradients, viscous dissipation, ion drag, the Coriolis force, and electrodynamic drifts. Among the findings of an initial steady state investigation are the following: (1) in convection channels, significant shears and rotations of the thermospheric flow can occur below 200 km if a minimum in the electron density profile is present between the E and F regions; (2) in convection channels, the thermospheric wind decreases with height in the F region owing to the effects of horizontal viscosity; and (3) at low altitudes, the boundaries of convection channels may produce Ekman spirals.

St-Maurice, J.-P.↗

A three-dimensional time-dependent model of the polar wind

A time-dependent three-dimensional multiion model of the polar wind was developed, which covers the altitude range of from 120 to 9000 km and takes into account supersonic ion outflow, shock formation, and ion energization during plasma expansion events. The model was used to study the temporal response of global polar wind to changing magnetospheric conditions, for the winter solstice and for solar-minimum conditions in the northern polar region. Graphs illustrating temporal changes with changes in T(e), T(i), and T(n) along the dawn, the trough, and the dusk convection trajectories and in the O(+), O, and H densities along the same convection trajectories are presented together with conntours of the H(+) and the O(+) densities along the three convection trajectories.

Schunk, R. W.↗

Carbon-Nitrogen-Oxygen Line Radiation and the X-ray Bowen Fluorescence Mechanism in Optically Thick, Highly Ionized Media

Radiative transfer effects due to overlapping X-ray lines in a high-temperature, optically thick, highly ionized medium are investigated. One particular example, in which the O VIII Lyalpha doublet (2(sup 2) P(sub 1/2,3/2)-1(sup 2) S(sub 1/2) coincides in frequency with the N VII Lyzeta lines (7(sup 2) P(sub 1/2,3/2)-1(sup 2) S(sub 1/2) is studied in detail to illustrate the effects on the properties of the emergent line spectrum. We solve the radiative transfer equation to study the energy transport of resonance-line radiation in a static, infinite, plane-parallel geometry, which is used to compute the destruction/escape probabilities for each of the lines for various total optical thicknesses of the medium, as well as destruction probabilities by sources of underlying photoelectric opacity. It is found that a large fraction of the O vIII Lyalpha line radiation can be destroyed by N VII, which can result in a reversal of the O VIII Lyalpha/N VII Lyalpha line intensity ratio similar to what may be seen under nonsolar abundances. Photoelectric absorption by ionized carbon and nitrogen can also subsequently increase the emission-line intensities of these ions. We show that line ratios, which are directly proportional to the abundance ratios in optically thin plasmas, are not good indicators of the true CNO abundances. Conversely, global spectral modeling that assumes optically thin conditions may yield incorrect abundance estimates when compared with observations, especially if the optical depth is large. Other potentially important overlapping lines and continua in the X-ray band are also identified, and their possible relevance to recent high-resolution spectroscopic observations with Chandra and XMM-Newton are briefly discussed.

Sako, Masao↗

Science Enhancements by the MAVEN Participating Scientists

NASA implemented a Participating Scientist Program and released a solicitation for the Mars Atmosphere and Volatile EvolutioN mission (MAVEN) proposals on February 14, 2013. After a NASA peer review panel evaluated the proposals, NASA Headquarters selected nine on June 12, 2013. The program's intent is to enhance the science return from the mission by including new investigations that broaden and/or complement the baseline investigations, while still addressing key science goals. The selections cover a broad range of science investigations. Included are: a patching of a 3D exosphere model to an improved global ionosphere-thermosphere model to study the generation of the exosphere and calculate the escape rates; the addition of a focused study of upper atmosphere variability and waves; improvement of a multi-fluid magnetohydrodynamic model that will be adjusted according to MAVEN observations to enhance the understanding of the solar-wind plasma interaction; a global study of the state of the ionosphere; folding MAVEN measurements into the Mars International Reference Ionosphere under development; quantification of atmospheric loss by pick-up using ion cyclotron wave observations; the reconciliation of remote and in situ observations of the upper atmosphere; the application of precise orbit determination of the spacecraft to measure upper atmospheric density and in conjunction with other Mars missions improve the static gravity field model of Mars; and an integrated ion/neutral study of ionospheric flows and resultant heavy ion escape. Descriptions of each of these investigations are given showing how each adds to and fits seamlessly into MAVEN mission science design.

Thermosphere↗

Science Enhancements by the MAVEN Participating Scientists

NASA implemented a Participating Scientist Program and released a solicitation for the Mars Atmosphere and Volatile EvolutioN mission (MAVEN) proposals on February 14, 2013. After a NASA peer review panel evaluated the proposals, NASA Headquarters selected nine on June 12, 2013. The program’s intent is to enhance the science return from the mission by including new investigations that broaden and/or complement the base-line investigations, while still addressing key science goals. The selections cover a broad range of science investigations. Included are: a patching of a 3D exosphere model to an improved global ionosphere-thermosphere model to study the generation of the exosphere and calculate the escape rates; the addition of a focused study of upper atmosphere variability and waves; improvement of a multi-fluid magnetohydrodynamic model that will be adjusted according to MAVEN observations to enhance the understanding of the solar-wind plasma interaction; a global study of the state of the ionosphere; folding MAVEN measurements into the Mars International Reference Ionosphere under development; quantification of atmospheric loss by pick-up using ion cyclotron wave observations; the reconciliation of remote and in situ observations of the upper atmosphere; the application of precise orbit determination of the spacecraft to measure upper atmospheric density and in conjunction with other Mars missions improve the static gravity field model of Mars; and an integrated ion/neutral study of ionospheric flows and resultant heavy ion escape. Descriptions of each of these investigations are given showing how each adds to and fits seamlessly into MAVEN mission science design.

Mars↗

Validation of prediction capability of operating space for plasma initiation in MAST-U

DYON is a plasma initiation modelling code that solves the differential equation system of the full circuit equations (plasma current, active coil currents and eddy currents in full passive structures) and 0D global energy and particle balance equations (Kim 2022 Nucl. Fusion 62 126012). In order to test the capability of the full electromagnetic plasma initiation model to predict individual discharges in experiments and thus the operating space in the device, a dedicated experimental database was built in MAST-U by scanning the prefilled gas pressure p 0 and the induced loop voltage V loop . In the experimental operating space of p 0 and V loop the lower and the upper limits of p 0 are determined by the plasma breakdown failure and the plasma burn-through failure, respectively. The lower limit of V loop is determined by the plasma burn-through failure. By directly reading the control room data used in each discharge (i.e. currents in the solenoid, poloidal field coils, and toroidal field coils, p 0 , and gas puffing rate), the full electromagnetic DYON consistently predicted the failed breakdown, failed burn-through, and successful plasma initiation discharges in the experimental database, demonstrating its capability to predict the operating space for inductive plasma initiation. The Paschen curve calculated with the effective connection length in MAST-U indicates a much higher p 0 required for plasma breakdown than the experimental data, indicating that individual field line evaluation is necessary to calculate the quantitative requirements for Townsend breakdown. The demonstration in this paper shows that the full electromagnetic DYON could be a useful simulation tool to assess the feasibility of inductive plasma initiation and to optimise operating scenarios in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analytical Ion Thruster Discharge Performance Model

A particle and energy balance model of the plasma discharge in magnetic ring-cusp ion thrusters has been developed. The model follows the original work of Brophy in the development of global 0-D discharge models that utilize conservation of particles into and out of the thruster and conservation of energy into the discharge and out of the plasma in the form of charged particles to the walls and beam and plasma radiation. The present model is significantly expanded over Brophy's original work by including self-consistent calculations of the internal neutral pressure, electron temperature, primary electron density, electrostatic ion confinement (due to the ring-cusp fields), plasma potential, discharge stability, and time dependent behavior during recycling. The model only requires information on the thruster geometry, ion optics performance and electrical inputs such as discharge voltage and currents, etc. to produce accurate performance curves of discharge loss versus mass utilization efficiency. The model has been benchmarked against the NEXIS Laboratory Model (LM) and Development Model (DM) thrusters, and successfully predicts the thruster discharge loss as a function of mass utilization efficiency for a variety of thrusters. The discharge performance model will be presented and results showing ion thruster performance and stability given.

discharge chambers↗

Reply to Comments by Tsurutani et al. on "Modeling Extreme 'Carrington-Type' Space Weather Events Using Three-Dimensional Global MHD Simulations"

In this response, we address the three main comments by Tsurutani et al. (2018, http://doi.org/10.1002/2017JA024779) namely, unusually high plasma density, interplanetary magnetic field intensity, and fast storm recovery phase. The authors agree that there is room to improve the modeling by taking into account these comments and other aspects that were not fully explored during our initial work. We are already in the process of undertaking a more comprehensive modeling project.

Ngwira, Chigomezyo M.↗

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗

Development of a global 3-D magnetohydrodynamic computational model for solar wind/cometary and planetary studies

Progress is reported in the development of a global, 3-D magnetohydrodynamic computational model to quantitatively describe the detailed continuum field and plasma interaction process of the solar wind with cometary and planetary magneto/ionopause shapes. Activities described include enhancements to the present Level 1 model, development of the 3-D gas dynamic flow field solvers for the Level 2 model, and continued development of the 3-D frozen magnetic field solver for the Level 2 model.

Stahara, S. S.↗

Artificial Intelligence/Deep Learning FRNN Software for Prediction & Real-Time Control of DIII-D Plasma Control System (PCS)

This collaborative project integrated an improved version of the Artificial Intelligence/Deep Learning FRNN prediction and control software into the real-time DIII-D PCS (plasma control system). A key AI/DL software challenge is to build a modern high-performance computing (HPC) enabled “synthetic plasma simulator” capable of carrying out HPC-driven real-time plasma control applications. This involves development of a deep learning framework to train the surrogate model for a first-principles-based instability analysis simulator (“SGTC”) derived from the global gyrokinetic code GTC. The role of SGTC is to provide accurate and detailed plasma instability information from a real-time AI-based simulator capability to complement the deep learning prediction and control from experimentally-measured signals, such as ECE Imaging, supplemented by synthetic SGTC-ECEI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AMPPED LDRD: Advanced Models of the Physics & Phenomena of Electrical Discharge.

This Laboratory Directed Research and Development (LDRD) project aimed at addressing an unmet need for physics-based models that enable prediction of electrical breakdown in the presence of high voltage gradients and real-world environments. This project focused on understanding coupled plasma-surface feedback phenomena critical to electrical discharge processes and design of breakdown voltages, breakdown times, and breakdown paths. New secondary electron emission and coupled plasma-surface measurement methods and models were developed. Highly reproducible discharge experiments were developed to provide current-voltage data as a function of time, and methods were developed to condition electrodes for predictable discharge. To assess plasma-surface feedback mechanisms, new surface physics capabilities, including measurement systems for local and global photoemission yield and ion-induced secondary electron emission yield were developed. Discharge models were developed to incorporate these yields for predicting discharge behavior to understand influence of surface physics on discharge current and time evolution, and methods to accelerate these simulations were demonstrated. A key finding was that the initiation process of a non-vacuum discharge appears highly dependent on photoemission at the cathode. Advanced diagnostics were developed, including spatial electrical field imaging, spectroscopy of plasmas, and high-speed imaging of discharge event s in metal and dielectric environments. These diagnostics and improved electrical models are a step toward enabling quantitative design margins, assurance of safety architectures, and implementation of novel materials and technologies for wide ranging energy and electronics applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The impact of kinetic and global effects on ideal ballooning 2nd stable pedestals of conventional and low aspect-ratio tokamaks

The EPED model (Snyder et al 2011 Nucl. Fusion 51 103016) has had success in describing the pedestal structure of Type-I ELM and QH-mode plasmas in tokamaks, by combining kinetic ballooning mode (KBM) and peeling–ballooning (PB) constraints. Within EPED, the KBM constraint is often approximated by a technique using calculated ideal ballooning mode (IBM) thresholds, fit to a functional form designed to capture non-local effects. It has been noted that quantitative differences between local ideal MHD and gyro-kinetic (GK) ballooning stability can be larger at low aspect ratio. KBM critical pedestals are consistent with observations in initial studies on conventional and spherical tokamaks. In this work, the application of a reduced model for the calculation of the kinetic ballooning stability boundary is presented based on a novel and newly developed gyro-fluid system (GFS) code (Staebler et al 2023 Phys. Plasmas 30 102501). GFS is observed to capture KBMs in DIII-D as well as the NSTX pedestals, opening a route to integrating this model into EPED. Finally, high but finite $n$ global ballooning modes are observed to limit the access to the local 2nd stability and thus provide a transport mechanism that constrains the pedestal evolution with $β_{p,\textrm{ped}}$. The high $n$ global ballooning stability is approximated by its ideal MHD analog using ELITE. It is shown that nearly-local high $n$ modes with $k_yρ_s$ ~ $0.25 - 0.5$ can provide a proxy for the critical $β_{p,\textrm{ped}}$ when a 2nd stable access exists on DIII-D plasmas. The use of GFS and ELITE scaling in EPED provides improved agreement to EPED1 in an initial comparison with DIII-D pedestal data.

Anastopoulos Tzanis, Michail Savvas [Oak Ridge Nat↗