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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Origin of Pulsed Radio Emission from Magnetars

Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10 30 erg s −1 , and a broad range of frequencies extending up to ∼100 GHz.

Astronomical simulations↗
DOE OSTI2021

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

The Flow of Plasma in the Solar-Terrestrial Environment

The overall goal of our NASA theory research is to trace the flow of mass, momentum, and energy through the magnetosphere-ionosphere-atmosphere system taking into account the coupling, time delays, and feedback mechanisms that are characteristic of the system. Our approach is to model the magnetosphere-ionosphere-atmosphere (M-I-A) system in a self-consistent quantitative manner using unique global models that allow us to study the coupling between the different regions on a range of spatial and temporal scales. The uniqueness of our global models stems from their high spatial and temporal resolutions, the physical processes included, and the numerical techniques employed. Currently, we have time-dependent global models of the ionosphere, thermosphere, polar wind, plasmasphere, and electrodynamics. It is now becoming clear that a significant fraction of the flow of mass, momentum, and energy in the M-I-A system occurs on relatively small spatial scales. Therefore, an important aspect of our NASA Theory program concerns the effect that mesoscale (100-l000 km) density structures have on the macroscopic flows in the ionosphere, thermosphere, and polar wind. The structures can be created either by structured magnetospheric inputs (i.e., structured electric field, precipitation, or Birkeland current patterns) or by time variations of these inputs due to geomagnetic storms and substorms. Some of the mesoscale structures of interest include sun-aligned polar cap arcs, propagating plasma patches, traveling convection vortices, subauroral ion drift (SAID) channels, gravity waves, and the polar hole.

Schunk, Robert W.↗
NASA NTRS1998

Streamer Evaporation

Evaporation is the consequence of slow plasma heating near the tops of streamers where the plasma is only weakly contained by the magnetic field. The form it takes is the slow opening of field lines at the top of the streamer and transient formation of new solar wind. It was discovered in polytropic model calculations, where due to the absence of other energy loss mechanisms in magnetostatic streamers, its ultimate endpoint is the complete evaporation of the streamer. This takes, for plausible heating rates, weeks to months in these models. Of course streamers do not behave this way, for more than one reason. One is that there are losses due to thermal conduction to the base of the streamer and radiation from the transition region. Another is that streamer heating must have a characteristic time constant and depend on the ambient physical conditions. We use our global Magnetohydrodynamics (MHD) model with thermal conduction to examine a few examples of the effect of changing the heating scale height and of making ad hoc choices for how the heating depends on ambient conditions. At the same time, we apply and extend the analytic model of streamers, which showed that streamers will be unable to contain plasma for temperatures near the cusp greater than about 2xl0(exp 6) K. Slow solar wind is observed to come from streamers through transient releases. A scenario for this that is consistent with the above physical process is that heating increases the near-cusp temperature until field lines there are forced open. The subsequent evacuation of the flux tubes by the newly forming slow wind decreases the temperature and heating until the flux tubes are able to reclose. Then, over a longer time scale, heating begins to again refill the flux tubes with plasma and increase the temperature until the cycle repeats itself. The calculations we report here are first steps towards quantitative evaluation of this scenario.

Suess, Steven T.↗
NASA NTRS1983

Global MHD model of the earth's magnetosphere

A global MHD model of the earth's magnetosphere is defined. An introduction to numerical methods for solving the MHD equations is given with emphasis on the shock-capturing technique. Finally, results concerning the shape of the magnetosphere and the plasma flows inside the magnetosphere are presented.

Wu, C. C.↗
NASA NTRS1985

Space plasma simulations; Proceedings of the Second International School for Space Simulations, Kapaa, HI, February 4-15, 1985. Parts 1 & 2

Space plasma simulations, observations, and theories are discussed. Papers are presented on the capabilities of various types of simulation codes and simulation models. Consideration is given to plasma waves in the earth's magnetotail, outer planet magnetosphere, geospace, and the auroral and polar cap regions. Topics discussed include space plasma turbulent dissipation, the kinetics of plasma waves, wave-particle interactions, whistler mode propagation, global energy regulation, and auroral arc formation.

Ashour-Abdalla, M.↗
NASA NTRS1988

The role of magnetohydrodynamics in heliospheric space plasma physics research

Magnetohydrodynamics (MHD) is a fairly recent extension of the field of fluid mechanics. While much remains to be done, it has successfully been applied to the contemporary field of heliospheric space plasma research to evaluate the 'macroscopic picture' of some vital topics via the use of conducting fluid equations and numerical modeling and simulations. Some representative examples from solar and interplanetary physics are described to demonstrate that the continuum approach to global problems (while keeping in mind the assumptions and limitations therein) can be very successful in providing insight and large scale interpretations of otherwise intractable problems in space physics.

Dryer, Murray↗
NASA NTRS2010

Martian Dust Devil Electron Avalanche Process and Associated Electrochemistry

Mars' dynamic atmosphere displays localized dust devils and larger, global dust storms. Based on terrestrial analog studies, electrostatic modeling, and laboratory work these features will contain large electrostatic fields formed via triboelectric processes. In the low-pressure Martian atmosphere, these fields may create an electron avalanche and collisional plasma due to an increase in electron density driven by the internal electrical forces. To test the hypothesis that an electron avalanche is sustained under these conditions, a self-consistent atmospheric process model is created including electron impact ionization sources and electron losses via dust absorption, electron dissociation attachment, and electron/ion recombination. This new model is called the Dust Devil Electron Avalanche Model (DDEAM). This model solves simultaneously nine continuity equations describing the evolution of the primary gaseous chemical species involved in the electrochemistry. DDEAM monitors the evolution of the electrons and primary gas constituents, including electron/water interactions. We especially focus on electron dynamics and follow the electrons as they evolve in the E field driven collisional gas. When sources and losses are self-consistently included in the electron continuity equation, the electron density grows exponentially with increasing electric field, reaching an equilibrium that forms a sustained time-stable collisional plasma. However, the character of this plasma differs depending upon the assumed growth rate saturation process (chemical saturation versus space charge). DDEAM also shows the possibility of the loss of atmospheric methane as a function of electric field due to electron dissociative attachment of the hydrocarbon. The methane destruction rates are presented and can be included in other larger atmospheric models.

Jackson, Telana L.↗
NASA NTRS1995

Heliospheric Termination Shock Motion Due to Fluctuations in the Solar Wind Upstream Conditions: Spherically Symmetric Model

Large-scale fluctuations in the solar wind plasma upstream of the heliospheric termination shock (TS) will cause inward and outward motions of the shock. Using numerical techniques, we extend an earlier strictly one-dimensional (planar) analytic gas dynamic model to spherical symmetry to investigate the features of global behavior of shock motion. Our starting point is to establish a steady numerical solution of the gasdynamic equations describing the interaction between the solar wind and the interstellar medium. We then introduce disturbances of the solar wind dynamic pressure at an inner boundary, and follow the subsequent evolution of the system, especially the motion of the termination shock. Our model solves spherically symmetric gasdynamic equations as an initial-boundary value problem. The equations in conservative form are solved using a fully implicit Total Variation Diminishing (TVD) upwind scheme with Roe-type Riemann solver. Boundary conditions are given by the solar wind parameters on an inner spherical boundary, where they are allowed to vary with time for unsteady calculations, and by a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. We find that immediately after the interaction, the shock moves with speeds given by the earlier analogous analytic models. However, as the termination shock propagates it begins to slow down, seeking a new equilibrium position. In addition, the disturbance transmitted through the TS, either a shock or rarefaction wave, will encounter the heliopause boundary and be reflected back. The reflected signal will encounter the TS, causing it to oscillate. The phenomenon may be repeated for a number of reflections, resulting in a "ringing" of the outer heliosphere.

Ratkiewicz, R.↗
NASA NTRS2022

Direct Observations of a Complex Coronal Web Driving Highly Structured Slow Solar Wind

The solar wind is formed by continuous streams of charged particles that escape into the heliosphere from the Sun. Fast wind with speeds more than 500 km s -1 emerges from the interiors of coronal holes. Near the ecliptic plane, fast wind from low-latitude coronal holes is interspersed with a highly structured slow solar wind. The source regions and drivers of this slow wind are poorly understood. Here we report observations from an off-point campaign of the GOES Solar Ultraviolet Imager that reveal a spatially complex web of magnetized plasma structures, over a low-latitude coronal hole active region system, that persistently interact and reconnect in the middle corona. SOHO coronagraphic images show concurrent emergence of slow wind streams over these coronal web structures. With advanced global MHD coronal models, we demonstrate that the observed coronal web is a direct imprint of the magnetic separatrix web (S-web). By revealing a highly dynamic portion of S-web, our observations open a window into important middle coronal processes that appear to play key role in driving structured slow solar wind.

Solar physics↗
DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Increasing the density in Wendelstein 7-X: benefits and limitations

In stellarators, increasing the density is beneficial for the energy confinement. While there is no single reason for this observation, it is still very robust across different devices and this is reflected in the empirical energy confinement time scaling for stellarators, ISS04. In order to study whether this is also true for Wendelstein 7-X, here the density scaling of the energy confinement time is analyzed and compared to ISS04 for the first divertor experiments. When the density is increased beyond a critical density, however, radiative collapses are frequently observed. Existing analytical models for the critical density are revisited to assess whether they can predict the accessible density range. Furthermore, since close to the collapse the radiation losses increase substantially, the impact on the global energy confinement is investigated. It is found that in plasmas with high radiation the density scaling of the energy confinement time becomes weaker, the reason for this observation is not yet clear. In the second half of the first divertor campaign, boronization was applied to W7-X for the first time. This broadened the operational window, allowing for operation at higher density and, hence, higher stored energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Real-time confinement regime detection in fusion plasmas with convolutional neural networks and high-bandwidth edge fluctuation measurements

Abstract A real-time detection of the plasma confinement regime can enable new advanced plasma control capabilities for both the access to and sustainment of enhanced confinement regimes in fusion devices. For example, a real-time indication of the confinement regime can facilitate transition to the high-performing wide-pedestal (WP) quiescent H-mode, or avoid unwanted transitions to lower confinement regimes that may induce plasma termination. To demonstrate real-time confinement regime detection, we use the 2D beam emission spectroscopy (BES) diagnostic system to capture localized density fluctuations of long wavelength turbulent modes in the edge region at a 1 MHz sampling rate. BES data from 330 discharges in either L-mode, H-mode, quiescent H (QH)-mode, or WP QH-mode were collected from the DIII-D tokamak and curated to develop a high-quality database to train a deep-learning classification model for real-time confinement detection. We utilize the 6×8 spatial configuration with a time window of 1024 µ s and recast the input to obtain spectral-like features via fast Fourier transform preprocessing. We employ a shallow 3D convolutional neural network for the multivariate time-series classification task and utilize a softmax in the final dense layer to retrieve a probability distribution over the different confinement regimes. Our model classifies the global confinement state on 44 unseen test discharges with an average F 1 score of 0.94, using only ∼1 ms snippets of BES data at a time. This activity demonstrates the feasibility for real-time data analysis of fluctuation diagnostics in future devices such as ITER, where the need for reliable and advanced plasma control is urgent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Toroidal modeling of Alfvén eigenmodes excited by runaway electrons in DIII-D and ITER

The non-perturbative MHD-kinetic hybrid code MARS-K (Liu et al 2014 Phys. Plasmas 21 056105) is updated to include relativistic effects for kinetic fast particles, enabling the code to model excitation of Alfvén eigenmodes (AEs) by runaway electrons (REs) in post-disruption tokamak plasmas. Applying the updated code to RE beams in both DIII-D and ITER, a zoo of AE modes triggered by trapped REs due to precessional drift-kinetic resonances is computed while scanning the RE energy. At fixed RE energy, multiple unstable roots are also excited. These AE modes possess radially different eigenmode structures, ranging from global modes to core-localized ones. The computed mode frequency is in the Alfvén frequency range, increasing with the assumed RE energy in a staircase fashion and quantitatively matching the experimental measurement (in DIII-D). At the (more relevant) high-frequency range (above 1 MHz), the modeled eigenmodes are identified as compressional AEs (CAEs) in DIII-D and a mixture of CAE and shear Alfvén waves in ITER.

Alfvén eigenmodes↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗