SEARCH · Engineering Papers

Results for “Theoretical plasma physics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Physical properties of comets

Until they reach positions a few AU from the sun, comets may be icy bodies a few km or less in diameter. At smaller distances from the sun, the nucleus of a comet will begin to vaporize, thereby generating a large neutral atmosphere, or coma. Two tails are also produced, one of which consists of dust particles entrained in the coma gas during vaporization and then separated from the coma by radiation pressure, while the other consists of ionized molecules swept back by the solar wind. Attention is presently given to the features noted by studies undertaken since 1980 in the nucleus, coma, dust tail and plasma tail regions of observable comets, as well as to the results of theoretical discussions concerning the origin and evolution of comets. Also evaluated are current preparations for the observation of Comet Halley's 1986 arrival at perihelion.

Butterworth, P. S.↗

DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗

DOE OSTI2026

Role of Neutrals Versus Transport in Determining the Pedestal Density Structure: Final Technical Report

In fusion devices the plasma density plays a crucial role in determining the fusion reaction rate and has a direct impact on the fusion gain of a given device. This density is in general regulated by the particle sources and transport near the plasma edge, which give rise to an edge density pedestal. When predicting the performance of future devices, this density pedestal is often prescribed, rather than predicted, due to a lack of models which allow confident extrapolation. This project aims to advance these models through the focused validation of theoretical models related to the transport of fueling neutral particles, and through interpretive transport modeling in present day fusion plasmas, in which the penetration of neutrals is altered to better simulate future reactor-like conditions. Achievements in theory and model validation under this project have advanced our understanding of how much of the edge density profile is set by transport versus direct ionization, enabling interesting projections to future burning plasma devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Transport barrier and spinning blob dynamics in the tokamak edge

Abstract In this work, we investigate the dynamics of plasma blobs in the edge of magnetic confinement devices using a full-f gyrokinetic particle-in-cell code with X-point geometry. In simulations, the evolution of a seeded blob is followed as it approaches a naturally-forming zonal shear layer near the separatrix, where the blob is stabilized by a large spin induced by the self-consistent adiabatic electron response, and blob bifurcation and trapping are observed during the cross-field propagation of blobs. A new theoretical explanation in both the zonal free and zonal shear layer is constructed, where the dominant E × B spin motion is included. A theoretical condition for a transport barrier induced by the interaction between spinning blobs and the zonal shear layer is obtained, and its scaling is verified with simulations. The new theoretical framework, especially the transport barrier, is applicable to explain and predict various experimental phenomena. In particular, the transport barrier condition calculated with experimental parameters demonstrates that the blob radial transport for H mode is smaller than L mode in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A Low Fuel Convergence Path to Inertial Confinement Fusion on the National Ignition Facility (Final Report)

After a decade of trying, the grand challenge problem of achieving hot spot ignition on NIF has not been achieved using indirect-drive of a high-convergence single shell capsule inside a hohlraum. Here, the alternate concept of obtaining volume-like ignition of liquid DT fuel via multi-shell implosion using laser polar direct drive (PDD) on NIF was assessed through theoretical analysis, computer simulation, fabrication development and implosion experiments on both the Laboratory for Laser Energetics’ (LLE) Omega laser and the National Ignition Facility (NIF). This research has shown that multi-shell ignition has merit and should be further advanced by the National ICF Program to determine its ultimate limitations for achieving ignition on NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

The dynamics of planetary rings

The physical processes that occur in planetary rings are discussed. The theoretical arguments leading to the conclusion that Saturn's rings are solid particles in nonuniform rotation are summarized, and the optical depth, thickness, and particle size of the rings are discussed. The influence of nearby satellites on the rings is analyzed, and asymmetries in the rings are briefly discussed. What is known of the rings of Uranus and of Jupiter's ring is summarized. Some of the dynamical processes and influences that are expected to be incorporated in more advanced theories of planetary rings are reviewed in detail, including radiation drag, plasma drag, interparticle collisions and viscosity, resonances with external satellites, shepherd satellites and moonlets. Finally, the orbital evolution of the shepherd satellites caused by the rings is estimated.

Goldreich, P.↗

NASA NTRS1991

Extension of convection modeling into the high-latitude ionosphere - Some theoretical difficulties

The Rice Convection Model (RCM) is extended and merged with empirical models so as to cover the entire high-latitude ionosphere with the aim of providing precipitation and electric field inputs for ionosphere and thermosphere modelers and producing a model in which the boundaries of the precipitation and electric field patterns maintain physically consistent relationships to each other. The computed auroral electron energy flux, plotted as a function of latitude, exhibited an exaggerated two-peak structure. When no floor was placed under the precipitation rate, the minimum between the two peaks was much too deep to be consistent with typical observations. The regions of excessively weak precipitation map to equatorial distances of 15-35 RE and thus to the regions of the plasma sheet that were not included in previous self-consistent convection calculations.

Wolf, R. A.↗

NASA NTRS2012

Space Weather Impacts on Spacecraft Design and Operations in Auroral Charging Environments

Spacecraft in low altitude, high inclination (including sun-synchronous) orbits are widely used for remote sensing of the Earth s land surface and oceans, monitoring weather and climate, communications, scientific studies of the upper atmosphere and ionosphere, and a variety of other scientific, commercial, and military applications. These systems are episodically exposed to environments characterized by a high flux of energetic (approx.1 to 10 s kilovolt) electrons in regions of very low background plasma density which is similar in some ways to the space weather conditions in geostationary orbit responsible for spacecraft charging to kilovolt levels. While it is well established that charging conditions in geostationary orbit are responsible for many anomalies and even spacecraft failures, to date there have been relatively few such reports due to charging in auroral environments. This presentation first reviews the physics of the space environment and its interactions with spacecraft materials that control auroral charging rates and the anticipated maximum potentials that should be observed on spacecraft surfaces during disturbed space weather conditions. We then describe how the theoretical values compare to the observational history of extreme charging in auroral environments and discuss how space weather impacts both spacecraft design and operations for vehicles on orbital trajectories that traverse auroral charging environments.

Minow, Joseph I.↗

DOE OSTI2026

New approaches to Bayesian uncertainty quantification for Nuclear Science (Final Technical Report)

Inverse problems play a central role in experimentation and theory/data comparisons for many areas of modern Nuclear Physics (NP) and High-Energy Physics (HEP). Bayes’s Theorem is a powerful tool for solving Inverse Problems, providing conceptually transparent and unbiased constraints on theoretical parameters and their uncertainties (“Bayesian Inference”) and enabling the quantification of agreement or tension between models and data. However, analyses based on Bayesian Inference are often challenging for NP and HEP applications, either because of the large number of parameters in the problem, the high computational cost, or both. We propose a multi-institutional collaboration to develop and deploy novel Bayesian analysis tools that advance the scientific scope of a broad range of current and future NP experiments. This project brings together NP domain scientists working on several high-profile NP projects for which new, high-performance Bayesian Uncertainty Quantification (“Bayesian UQ”) methods are essential to carry out the science, and data scientists who are developing state-of-the-art methods applicable to these problems. The NP projects in this proposal comprise measurements of the mass and fundamental nature of the neutrino; study of the Quark-Gluon Plasma that filled the early universe; and mapping of natural and anthropogenic radiation environments. While these NP projects have very different scientific goals, with datasets and analysis approaches that differ significantly, they share common requirements for improving computationally intensive Bayesian analyses using advanced Machine Learning algorithms and will benefit strongly from a coherent effort to develop general solutions. This proposal brings together these projects and forefront ML-based data science algorithms to develop such general solutions. The methods developed in this project will also be more widely applicable, thereby advancing science in the larger Nuclear Physics portfolio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

3D Equilibrium Reconstruction and MHD Studies in Stellarators (Final Report)

This is the final technical report for the grant Three Dimensional Equilibrium Reconstruction and MHD Studies in Stellarators. The major goals of this project were to: further develop V3FIT code capabilities, continue to support experiments using the V3FIT code, continue theoretical support for the CTH experiment at Auburn, particularly in the area of NIMROD simulations of CTH discharges. Dr. Lopez Ortiz has made significant progress on his analysis of NIMROD simulations of CTH MHD phenomena. In particular, his study has focused on the fixed point dynamics and energy flow amongst toroidal Fourier modes in sawtooth oscillations simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

What is the Quark-Gluon Plasma made of?

This article surveys our present understanding of the internal structure of the fully developed quark-gluon plasma at temperatures outside the crossover region. The theoretical part of the review covers perturbative and nonperturbative approaches to quark-gluon plasma structure, in particular, hard-thermal loop effective theory, lattice QCD and the functional renormalization group. The phenomenological part of the review scrutinizes the information that has been derived from bulk observables and hard probes in relativistic heavy ion collisions in terms of how it informs our knowledge about the structure of the quark-gluon plasma. The final section lists possible avenues for future progress.

FOS: Physical sciences↗

NASA NTRS1999

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗

NASA NTRS2016

Fast Ionized X-Ray Absorbers in AGNs

We investigate the physics of the X-ray ionized absorbers often identified as warm absorbers (WAs) and ultra-fast outflows (UFOs) in Seyfert AGNs from spectroscopic studies in the context of magnetically-driven accretion-disk wind scenario. Launched and accelerated by the action of a global magnetic field anchored to an underlying accretion disk around a black hole, outflowing plasma is irradiated and ionized by an AGN radiation field characterized by its spectral energy density (SED). By numerically solving the Grad-Shafranov equation in the magnetohydrodynamic (MHD) framework, the physical property of the magnetized disk-wind is determined by a wind parameter set, which is then incorporated into radiative transfer calculations with xstar photoionization code under heating-cooling equilibrium state to compute the absorber's properties such as column density N(sub H), line-of-sight (LoS) velocity v, ionization parameter xi, among others. Assuming that the wind density scales as n varies as r(exp. -1), we calculate theoretical absorption measure distribution (AMD) for various ions seen in AGNs as well as line spectra especially for the Fe K alpha absorption feature by focusing on a bright quasar PG 1211+143 as a case study and show the model's plausibility. In this note we demonstrate that the proposed MHD-driven disk-wind scenario is not only consistent with the observed X-ray data, but also help better constrain the underlying nature of the AGN environment in a close proximity to a central engine.

accretion↗

NASA NTRS1980

Plasma physics analysis of SERT-2 operation

An analysis of the major plasma processes involved in the SERT 2 spacecraft experiments was conducted to aid in the interpretation of recent data. A plume penetration model was developed for neutralization electron conduction to the ion beam and showed qualitative agreement with flight data. In the SERT 2 configuration conduction of neutralization electrons between thrusters was experimentally demonstrated in space. The analysis of this configuration suggests that the relative orientation of the two magnetic fields was an important factor in the observed results. Specifically, the opposed field orientation appeared to provide a high conductivity channel between thrusters and a barrier to the ambient low energy electrons in space. The SERT 2 neutralizer currents with negative neutralizer biases were up to about twice the theoretical prediction for electron collection by the ground screen. An explanation for the higher experimental values was a possible conductive path from the neutralizer plume to a nearby part of the ground screen. Plasma probe measurements of SERT 2 gave the clearest indication of plasma electron temperature, with normal operation being near 5 eV and discharge only operation near 2 eV.

Kaufman, H. R.↗

DOE OSTI2025

Bayesian inference analysis of jet quenching using inclusive jet and hadron suppression measurements

The JETSCAPE Collaboration reports a new determination of the jet transport parameter $\hat{q}$ in the quark-gluon plasma (QGP) using Bayesian inference, incorporating all available inclusive hadron and jet yield suppression data measured in heavy-ion collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC). This multi-observable analysis extends the previously published JETSCAPE Bayesian inference determination of $\hat{q}$, which was based solely on a selection of inclusive hadron suppression data. jetscape is a modular framework incorporating detailed dynamical models of QGP formation and evolution, and jet propagation and interaction in the QGP. Virtuality-dependent partonic energy loss in the QGP is modeled as a thermalized weakly coupled plasma, with parameters determined from Bayesian calibration using soft-sector observables. This Bayesian calibration of $\hat{q}$ utilizes active learning, a machine-learning approach, for efficient exploitation of computing resources. The experimental data included in this analysis span a broad range in collision energy and centrality, and in transverse momentum. In order to explore the systematic dependence of the extracted parameter posterior distributions, several different calibrations are reported, based on combined jet and hadron data; on jet or hadron data separately; and on restricted kinematic or centrality ranges of the jet and hadron data. Tension is observed in comparison of these variations, providing new insights into the physics of jet transport in the QGP and its theoretical formulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗

NASA NTRS1982

Highlights of theoretical progress related to the International Magnetospheric Study

U.S. theoretical research efforts have addressed three areas within the International Magnetospheric Study. The first, solar wind/magnetosphere interaction, is presently concerned with the suggestion that magnetic merging may predominantly occur near the polar cusps rather than near the subsolar point. Mechanisms have been proposed for noncollisional diffusion of solar wind plasma across the closed magnetopause entailed by such a phenomenon. The second area considers the importance to magnetotail dynamics of a continuous source of solar wind plasma, and of sporadic plasma loss associated with an unsteady convection cycle. In the third area, the electrodynamic magnetosphere/ionosphere interaction, an advanced state has been reached in the understanding of the relevant physics, with respect both to coupling in the subauroral region and the large scale structure of auroral zone electric fields parallel, and perpendicular to, the magnetic field.

Hill, T. W.↗