Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Asymptotic preserving methods for fluid electron-fluid models in the large magnetic field limit with mathematically guaranteed properties (Final Report)

The current manuscript is a final report on the activities carried out under the Project LDRD-CIS #226834. In scientific terms, the work reported in this manuscript is a continuation of the efforts started with Project LDRD-express #223796 with final report of activities SAND2021-11481, see [83]. In this section we briefly explain what pre-existing developments motivated the current body of work and provide an overview of the activities developed with the funds provided. The overarching goal of the current project LDRD-CIS #226834 and the previous project LDRD-express #223796 is the development of numerical methods with mathematically guaranteed properties in order to solve the Euler-Maxwell system of plasma physics and generalizations thereof. Even though Project #223796 laid out general foundations of space and time discretization of Euler-Maxwell system, overall, it was focused on the development of numerical schemes for purely electrostatic fluid-plasma models. In particular, the project developed a family of schemes with mathematically guaranteed robustness in order to solve the Euler-Poisson model. This model is an asymptotic limit where only electrostatic response of the plasma is considered. Its primary feature is the presence of a non-local force, the electrostatic force, which introduces effects with infinite speed propagation into the problem. Even though instantaneous propagation of perturbations may be considered nonphysical, there are plenty of physical regimes of technical interest where such an approximation is perfectly valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Hot ion plasma heating experiments in SUMMA

Initial results are presented for the hot-ion plasma heating experiments conducted in the new SUMMA (superconducting magnetic mirror apparatus) at NASA Lewis Research Center. A discharge is formed by applying a radially inward dc electric field between cylindrical anodes and hallow cathodes located at the peak of the mirrors. Data were obtained at midplane magnetic field strengths from 1.0 to 3.5 tesla. Charge-exchange neutral particle energy analyzer data were reduced to ion temperatures using a plasma model that included a Maxwellian energy distribution superimposed on an azimuthal drift, finite ion orbits, and radial variations in density and electric field. The best ion temperatures in a helium plasma were 5 keV and in hydrogen the H2(+) and H(+) ions were 1.2 keV and 1 keV respectively. Optical spectroscopy line broadening measurements yielded ion temperatures about 50 percent higher than the charge-exchange neutral particle analyzer results. Spectroscopically obtained electron temperature ranged from 3 to 30 eV. Ion temperature was found to scale roughly linearly with the ratio of power input-to-magnetic field strength, P/B.

Reinmann, J. J.↗
DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Framework for assessment of magnetic equilibrium controller performance on the MAST upgrade spherical tokamak

Here, in this work we present the assessment framework for magnetic equilibrium controllers on MAST Upgrade spherical tokamak (MAST-U) spherical tokamak. Such controllers are essential for the MAST-U since exhaust physics and core-edge integration studies require advanced divertor plasma configurations. The developed framework is based on the TokSys suite of plasma control codes, which was adapted and upgraded for MAST-U. However, extra capabilities were added on top of TokSys to support the development of new control algorithms, deployment of controllers to the plasma control system (PCS) and evaluation of their performance. The controller assessment was realized via closed-loop integrated control simulations with the actual MAST-U PCS and different physics-based plasma models. Since all components of the assessment chain were experimentally validated, these simulations provide qualified controllers applicable for direct use in the experiment. This resulted in the successful experimental demonstration of advanced plasma shape control on MAST-U with minimal on-machine development time. A similar methodology would be beneficial to other tokamaks, both existing and future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Vacuum ultraviolet line radiation measurements of a shock-heated nitrogen plasma

Line radiation, in the wavelength region from 1040 to 2500 A from nitrogen plasmas, was measured at conditions typical of those produced in the shock layer in front of vehicles entering the earth's atmosphere at superorbital velocities. The radiation was also predicted with a typical radiation transport computer program to determine whether such calculations adequately model plasmas for the conditions tested. The results of the comparison show that the radiant intensities of the lines between 1040 and 1700 A are actually lower than are predicted by such computer models.

Mcclenahan, J. O.↗
NASA NTRS1980

An extended soft X-ray source in Delphinus - H2027+19

A new extended soft X-ray source has been observed with the HEAO 1 A-2 experiment. The source, H2027+19, emits primarily in the 0.16-0.4 keV band with a total flux in this band of 2 x 10 to the -11th erg/sq cm s. It is found that both simple continuum and coronal plasma models provide good fits to the observed pulse-height spectrum. The most likely physical models are either that the source is an old supernova remnant or that it is a region of enhanced soft X-ray emission surrounding an H I cloud imbedded in a coronal plasma, as suggested by Hayakawa et al. (1979) for the Lupus Loop.

Stern, R. A.↗
NASA NTRS1986

Ionospheric effects to antenna impedance

The reciprocity between high power satellite antennas and the surrounding plasma are examined. The relevant plasma states for antenna impedance calculations are presented and plasma models, and hydrodynamic and kinetic theory, are discussed. A theory from which a variation in antenna impedance with regard to the radiated power can be calculated for a frequency range well above the plasma resonance frequency is give. The theory can include photo and secondary emission effects in antenna impedance calculations.

Bethke, K. H.↗
NASA NTRS2020

X-ray Spectra of the Fe-L Complex II: Atomic Data Constraints from the EBIT Experiment and X-Ray Grating Observations of Capella

The Hitomi results for the Perseus cluster have shown that accurate atomic models are essential to the success of X-ray spectroscopic missions and just as important as the store of knowledge on instrumental calibration and astrophysical modeling. Preparing the models requires a multifaceted approach, including theoretical calculations, laboratory measurements, and calibration using real observations. In a previous paper, we presented a calculation of the electron impact cross sections on the transitions forming the Fe-L complex. In the present work, we systematically tested the calculation against cross-sections of ions measured in an electron beam ion trap experiment. A two-dimensional analysis in the electron beam energies and X-ray photon energies was utilized to disentangle radiative channels following dielectronic recombination, direct electron-impact excitation, and resonant excitation processes in the experimental data. The data calibrated through laboratory measurements were further fed into a global modeling of the Chandra grating spectrum of Capella. We investigated and compared the fit quality, as well as the sensitivity of the derived physical parameters to the underlying atomic data and the astrophysical plasma modeling. We further list the potential areas of disagreement between the observations and the present calculations, which, in turn, calls for renewed efforts with regard to theoretical calculations and targeted laboratory measurements.

Liyi Gu↗
NASA NTRS1996

Modelling the Soft X-Ray and EUV Emission in Classical Novae: EUVE and ROSAT Observations of V1974 Cygni

We have conducted an extensive analysis of the observability of Classical Novae with the EUVE Lex/B and Al/Ti/C detectors. Predicted count rates have been computed using optically thin, isothermal plasma models for solar and metal-rich compositions, and hot ONeMg white dwarf model atmospheres. We find EUVE to be quite sensitive to both the EUV and soft X-ray emission emitted by the underlying hot white dwarf during novae outbursts, except for the coolest temperatures with very high intervening hydrogen column density. These results are used to interpret the emission detected during the EUVE all-sky survey of Nova Cygni 1992 (approx. = V1974 Cyg), 279-290 days after visual maximum. We find the best fit to the observed emission from V1974 Cyg arises from a hot ONeMg white dwarf with surface temperature approx. 4 x l0(exp 5) K and a mas of approx. 1.2 solar mass, and derive an interstellar hydrogen column density of approx. 3 x 10(exp 21)/sq cm. Virtually all this emission arises from supersoft X-rays rather than the EUV. We also report the detection of V1974 Cyg with the EUVE Deep Survey detector at 549 days after visual maximum. This observation is compatible with the above properties, indicating that the mechanism responsible for the soft X-ray emission, connected with the underlying white dwarf, had not yet entirely turned off. We also present analysis of a ROSAT PSPC observation which is contemporaneous with the EUVE survey observations; this independently confirms the high column density we derived from the EUVE survey observation. Light curves for the EUVE and ROSAT observations are presented. Statistical tests for variability show that all of these observations are indeed highly variable over various time scales. The EUVE survey data shows one day variations, the EUVE DS data show approx. 30 minute fluctuations, while the ROSAT data vary rapidly on time scales of seconds. The EUVE data shows no periodic variability on any time scale. The implications of the rapid variability are briefly discussed.

Stringfellow, Guy S.↗
DOE OSTI2024

Impurity transport in PISCES-RF

Linear plasma devices (LPD) utilizing a helicon plasma source, a high density light ion source, can generate impurities due to progressive erosion of the radio frequency (RF) transmission window caused by rectified sheath voltage. These source-born impurities can entrain and be transported by the plasma toward a target, affecting plasma-material interaction studies. Earlier work on material testing in Prototype-Materials Plasma Exposure eXperiment at ORNL revealed significant source impurity deposition on downstream targets. However, using a similar RF source, no target impurity deposition is observed in Plasma Interaction Surface Component Experimental Station (PISCES)-RF despite evidence of RF window erosion in the source region, thereby motivating the present work. Experimentally, using various magnetic field configurations upstream of the PISCES-RF plasma source and seeding titanium (Ti) impurities at various axial locations, impurity transport and deposition along the machine axis were investigated. It was found that Ti deposition was localized to the side of the plasma source where the Ti impurity was seeded. In contrast, aluminum (Al) deposition, originating from the sputtering of the helicon window, occurred predominantly upstream of the plasma source, suggesting an asymmetry in the axial transport of eroded RF window material. These observations suggest a stagnation of the parallel plasma flow immediately downstream of the plasma source, with impurity ions remaining unmagnetized near the source upstream. Al deposition in magnetic field-free regions in PISCES-RF indicates that sputtered Al impurities likely remained neutral due to their large ionization mean-free path under PISCES-RF conditions. Plasma modeling and simulation supported this, indicating that Al-neutrals transport toward the helicon source upstream for low electron density cases. It was found that the Larmor radius of the Al ions was greater than the plasma radius towards the source upstream and remained weakly magnetized in PISCES-RF, meaning that plasma source-born impurities are not efficiently entrained in the plasma flow. These findings provide critical insights into impurity transport in helicon plasma-based LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modeling of the cold electron plasma density for radiation belt physics

This review focusses strictly on existing plasma density models, including ionospheric source models, empirical density models, physics-based and machine-learning density models. This review is framed in the context of radiation belt physics and space weather codes. The review is limited to the most commonly used models or to models recently developed and promising. A great variety of conditions is considered such as the magnetic local time variation, geomagnetic conditions, ionospheric source regions, radial and latitudinal dependence, and collisional vs. collisionless conditions. These models can serve to complement satellite observations of the electron plasma density when data are lacking, are for most of them commonly used in radiation belt physics simulations and can improve our understanding of the plasmasphere dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

Systematic Comparison of Photoionized Plasma Codes with Application to Spectroscopic Studies of AGN in X-Rays

Atomic data and plasma models play a crucial role in the diagnosis and interpretation of astrophysical spectra, thus influencing our understanding of the Universe. In this investigation we present a systematic comparison of the leading photoionization codes to determine how much their intrinsic differences impact X-ray spectroscopic studies of hot plasmas in photoionization equilibrium. We carry out our computations using the Cloudy, SPEX, and XSTAR photoionization codes, and compare their derived thermal and ionization states for various ionizing spectral energy distributions. We examine the resulting absorption-line spectra from these codes for the case of ionized outflows in active galactic nuclei. By comparing the ionic abundances as a function of ionization parameter, we find that on average there is about 30 deviation between the codes in where ionic abundances peak. For H-like to B-like sequence ions alone, this deviation in is smaller at about 10 on average. The comparison of the absorption-line spectra in the X-ray band shows that there is on average about 30 deviation between the codes in the optical depth of the lines produced at log 1 to 2, reducing to about 20 deviation at log 3. We also simulate spectra of the ionized outflows with the current and upcoming high-resolution X-ray spectrometers, on board XMM-Newton, Chandra, Hitomi, and Athena. From these simulations we obtain the deviation on the best-fit model parameters, arising from the use of different photoionization codes, which is about 10 to40. We compare the modeling uncertainties with the observational uncertainties from the simulations. The results highlight the importance of continuous development and enhancement of photoionization codes for the upcoming era of X-ray astronomy with Athena.

plasmas – atomic processes – atomic data – t↗
DOE OSTI2024

Excited states in warm and hot dense matter

Accurate modeling of warm and hot dense matter is challenging in part due to the multitude of excited states that must be considered. Here, we present a variational framework that models these excited states. In this framework an excited state is defined by a set of effective one-electron occupation factors, and the corresponding energy is defined by the effective one-body energy with an exchange and correlation term. Here, the variational framework is applied to an atom-in-plasma model (a generalization of the so-called average atom model). Comparisons with a density functional theory based average atom model generally reveal good agreement in the calculated pressure, but our model also gives access to the excitation energies and charge state distributions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

A mirror instability associated with newly created ions in a moving plasma

Using a plasma model in which the ambient magnetic field is parallel to the z axis and the wave vector, an instability that resembles the usual hydromagnetic mirror instability is demonstrated. The instability is caused by freshly created ions in a moving plasma; the source of its free energy is associated with the ring distribution of the newly created ions, resulting in a high ion kinetic temperature in the direction transverse to the ambient magnetic field. The results of a stability analysis indicate that this instability can lead to the amplification of magnetosonic waves.

Wu, C. S.↗
NASA NTRS1991

Particle flows to shape and voltage surface discontinuities in the electron sheath surrounding a high voltage solar array in LEO

This paper discusses the numerical modeling of electron flows from the sheath surrounding high positively biased objects in LEO (Low Earth Orbit) to regions of voltage or shape discontinuity on the biased surfaces. The sheath equations are derived from the Two-fluid, Warm Plasma Model. An equipotential corner and a plane containing strips of alternating voltage bias are treated in two dimensions. A self-consistent field solution of the sheath equations is outlined and is pursued through one cycle. The electron density field is determined by numerical solution of Poisson's equation for the electrostatic potential in the sheath using the NASCAP-LEO relation between electrostatic potential and charge density. Electron flows are calculated numerically from the electron continuity equation. Magnetic field effects are not treated.

Metz, Roger N.↗
NASA NTRS1998

Physics of the Inner Heliosphere 1-10 Rs: Plasma Diagnostics and Models

While the mechanisms responsible for the solar corona and the high-speed solar wind streams are still unknown, model computations offer means of predicting the properties of such mechanisms in light of the empirical constraints currently available. Modeling and data analysis efforts were aimed at understanding the plasma properties of the acceleration of the solar wind, its filamentary nature, and the conditions needed to account for a rapidly accelerating solar wind, reaching its terminal speed within 10 R(sub s). A sequence of models ranging from steady one-fluid descriptions of the solar wind to multi-fluid time-dependent models were developed. Plasma diagnostics evolved from the analysis of data acquired from Skylab to SOHO, and complemented by ground-based observations.

Habbal, Shadia R.↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗