Engineering Papers⌕ Search

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.

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

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1988

Computational design of an experimental laser-powered thruster

An extensive numerical experiment, using the developed computer code, was conducted to design an optimized laser-sustained hydrogen plasma thruster. The plasma was sustained using a 30 kW CO2 laser beam operated at 10.6 micrometers focused inside the thruster. The adopted physical model considers two-dimensional compressible Navier-Stokes equations coupled with the laser power absorption process, geometric ray tracing for the laser beam, and the thermodynamically equilibrium (LTE) assumption for the plasma thermophysical and optical properties. A pressure based Navier-Stokes solver using body-fitted coordinate was used to calculate the laser-supported rocket flow which consists of both recirculating and transonic flow regions. The computer code was used to study the behavior of laser-sustained plasmas within a pipe over a wide range of forced convection and optical arrangements before it was applied to the thruster design, and these theoretical calculations agree well with existing experimental results. Several different throat size thrusters operated at 150 and 300 kPa chamber pressure were evaluated in the numerical experiment. It is found that the thruster performance (vacuum specific impulse) is highly dependent on the operating conditions, and that an adequately designed laser-supported thruster can have a specific impulse around 1500 sec. The heat loading on the wall of the calculated thrusters were also estimated, and it is comparable to heat loading on the conventional chemical rocket. It was also found that the specific impulse of the calculated thrusters can be reduced by 200 secs due to the finite chemical reaction rate.

Jeng, San-Mou↗
DOE OSTI2020

Application of laser-driven capacitor-coil to target normal sheath acceleration

A laser-driven accelerator generates protons with tens of MeV in energy by a compact, strong, and transient accelerating electric field produced as a result of laser–plasma interactions at relativistic intensities. In previous studies, two- and three-dimensional particle-in-cell simulations revealed that the application of a kT-level axial magnetic field results in an enhancement of proton acceleration via the target normal sheath acceleration mechanism due to reduced lateral electron divergence and improved electron heating efficiency. An experimental investigation of this scheme on the GEKKO-XII and the LFEX facilities found that the number and maximum energy of the accelerated protons decreased with increasing the temporal delay between the pulse driving the external magnetic-field and the pulse accelerating the protons, contrary to the theoretical and numerical expectations. Finally, we identify sources responsible for the degradation of the proton beam performance and we propose an alternative experimental setup to mitigate the degradation in future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Magnetic Reconnection in Strongly-Magnetized, Weakly-Collisional Plasmas: Onset, Turbulence, and Energy-Partition in 3D, Plasmoid-Dominated Regimes

This project is concerned with theoretical and numerical studies of magnetic reconnection. Key areas of focus are the reconnection onset (how reconnection gets started), energy partition (how the magnetic energy gets distributed during a reconnection event), and the effect of background turbulence on reconnection. We aim to shed light on these issues in kinetic plasmas (i.e., plasmas where the collisional mean free path is long compared to the length-scales of interest), with emphasis on 3D geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Magnetic mapping and Birkeland currents in the Toffoletto-Hill and Tsyganenko magnetosphere models

We investigate the mapping of magnetic flux tubes between the high-latitude ionosphere and the equatorial magnetosphere, using the Toffoletto and Hill (1989) theoretical model and the Tsyganenko (1987, 1989) empirical models. In all models, strong distortion of flux tubes occurs near the magnetopause and in the distant tail. These effects are illustrated for various interplanetary magnetic field directions in the Toffoletto and Hill (1989) model and for various Kp values in the Tsyganenko (1987, 1989) model. We calculate the Birkeland current that is actually present in each model (from Ampere's law), as well as the Birkeland current that would be implied if the model were in magnetostatic equilibrium with isotropic plasma pressure (from the Vasyliunas equation). Comparison of the actual and implied Birkeland currents indicates that no physically significant Birkeland current is included in any of the models tested. Each model contains a 'noise' component of Birkeland current that is, however, much smaller in magnitude than the 'real' Birkeland currents that remain to be included.

Ding, Cheng↗
DOE OSTI2026

Overview of ST40 results and future: expanding the physics basis of high-field spherical tokamaks

The goal of the ST40 programme is to explore the physics of high-field spherical tokamaks (STs), to validate empirical and theoretical models and, hence, to build confidence in predictions required to support the design of future generations of STs. ST40 is a compact high-field ST that has achieved the following parameters: R 0 = 0.4–0.55 m, I p = 0.20–0.85 MA, B t (R= 0.4 m) = 0.7–2.1 T, κ ⩽ 1.9, and A = 1.6–1.9. Highlights of recent experimental results include (i) H-mode and confinement studies at B t ⩽ 2.1 T, (ii) observation of bifurcation of the scrape-off-layer power fall-off width, λ q , into a ‘wide’ branch that follows existing H-mode scalings and a ‘narrow’ branch that exhibits λ q values that are up to 10 times lower than the predictions of established scalings, (iii) development of high-performance scenarios with plasma current, I p , up to 0.85 MA, (iv) development of highly non-inductive scenarios with high β p , and (v) the first ST40 experiments utilising the newly commissioned impurity powder dropper. The work on all these topics has been supported by a number of advancements in ST40 hardware and software, from plasma control to data analysis and interpretation. At the end of 2025, ST40 embarked on a major upgrade to further expand its capabilities by introducing, among other improvements, all-metal plasma-facing components, 1 MW of electron cyclotron heating, a pellet injector, and a pair of lithium evaporators for wall conditioning.

confinement↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fast Particle-Wave Interactions and Alfvén Eigenmodes in JET Tokamak Plasmas

This document serves as the closeout report of DOE Grant Award No. DE-FG02-99ER54563 with project period 1 April 2015 through 31 March 2020. The project comprised the international collaboration between MIT and EU scientists on the JET facility to improve our understanding of the physics of energetic particle-wave interactions by measuring the damping rates of stable Alfvén Eigenmodes (AEs) and unstable energetic particle driven modes. More specifically, this project involved the continued participation of MIT, the Culham Centre for Fusion Energy, the Swiss Plasma Center, and theorists from various European laboratories and from UC Irvine. Past contributions from the University of São Paulo are gratefully acknowledged. The Alfvén Eigenmode Active Diagnostic (AEAD) on JET discharges was successfully upgraded in this period and has extracted physically useful information, in primarily deuterium plasma discharges. It is expected that such experiments would be continued during the DT campaign in CY 2021 to assess the damping rates of similar modes in the presence of alpha particles, a product of fusing burning plasma. Only JET would carry out such experiments in the near term in the world. It is important to note that this collaboration was continued under DOE Grant Award No. DE-SC0014264 from 1 April 2020 through the present. In this grant period, the AEAD was upgraded with a set of individual amplifiers for each of a set of six antennas in two toroidally opposite locations. These new amplifiers allowed targeted selection of the antennas’ toroidal spectrum for AEs of interest with toroidal mode number | n | ≤ 20. The resonant detection and measurement of the damping rates of AEs was obtained from magnetic probes that could be compared with theory and simulations. This is a key topic of investigation for ITER and all other next-step fusion experiments, where such modes will interact with energetic particles produced by fusion reactions (alpha particles), Neutral Beam Injection (NBI) and Ion Cyclotron Resonance Heating (ICRH). The majority of the observations during this period were of Toroidal Alfvén Eigenmodes (TAEs) as these are most commonly observed on JET; a number of measurements were made during dedicated TAE experiments. However, a new set of lower frequency band filters were procured with the goal of studying Geodesic Acoustic Modes (GAMs), Beta (Acoustic) AEs (BAE/BAAEs), and Reverse Shear AEs (RSAEs), also predicted by theory. Commissioning and optimizations were successfully completed in this period, and the diagnostic has been in operation during the more recent JET campaigns. Initial results have been obtained in dedicated TAE experiments and successfully compared to drift-kinetic theory. A wide range of theoretical studies have been undertaken in support of the upcoming JET campaigns experiments. To supplement the ongoing use of ideal MHD codes, such as MISHKA during the studies of TAEs, gyrokinetic simulations of low frequency AEs were performed in collaboration with UC Irvine. The Gyrokinetic Toroidal Code (GTC) was used to determine the structure, frequencies and stability of AEs in JET plasmas. Thus, a solid scientific foundation has been laid for future DT campaigns in JET in the 2021 operation period when the damping rates of relevant modes in the presence of alpha particles could be assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Heavy quark transport: a theoretical overview

A brief overview on heavy-quark-QGP interactions is presented in this work. Theories at different momentum scales are summarized, including the flavor hierarchy of jet energy loss at high pT, hadronization and the charmed hadron chemistry at medium to low pT, and the extraction of the heavy-quark-QGP potential at low pT. Additionally, recent theoretical developments on heavy quark theories are discussed, such as multi-scale transport models of heavy quarks, probing initial state effects and cold nuclear matter effects using heavy quarks, as well as extending studies from single heavy flavor hadrons to heavy flavor jets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Integrating novel stellarator single-stage optimization algorithms to design the Columbia stellarator experiment

Abstract The Columbia Stellarator eXperiment (CSX), currently being designed at Columbia University, aims to test theoretical predictions related to QA plasma behavior, and to pioneer the construction of an optimized stellarator using three-dimensional, non-insulated high-temperature superconducting (NI-HTS) coils. The magnetic configuration is generated by a combination of two circular planar poloidal field (PF) coils and two 3D-shaped interlinked (IL) coils, with the possibility to add windowpane coils to enhance shaping and experimental flexibility. The PF coils and vacuum vessel are repurposed from the former Columbia Non-Neutral Torus experiment, while the IL coils will be custom-wound in-house using NI-HTS tapes. To obtain a plasma shape that meets the physics objectives with a limited number of coils, novel single-stage optimization techniques are employed, optimizing both the plasma and coils concurrently, in particular targeting a tight aspect ratio QA plasma and minimized strain on the HTS tape. Despite the increased complexity due to the expanded degrees of freedom, these methods successfully identify optimized plasma geometries that can be realized by coils meeting engineering specifications. This paper discusses the derivation of the constraints and objectives specific to CSX, and describe how two recently developed single-stage optimization methodologies are applied to the design of CSX. A set of selected configurations for CSX is then described in detail.

Baillod, A. (ORCID:0000000303529180)↗
DOE OSTI2023

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Exploring theoretical uncertainties in the hydrodynamic description of relativistic heavy-ion collisions

We explore theoretical uncertainties in the hydrodynamic description of relativistic heavy-ion collisions by examining the full nonlinear causality conditions and quantifying the second-order transport coefficients' role on flow observables. The causality conditions impose physical constraints on the maximum allowed values of inverse Reynolds numbers during the hydrodynamic evolution. Including additional second-order gradient terms in the Denicol-Niemi-Molnár-Rischke (DNMR) theory significantly shrinks the casual regions compared to those in the Israel-Stewart hydrodynamics. For Au + Au collisions, we find the variations of flow observables are small with and without imposing the necessary causality conditions, suggesting a robust extraction of the quark-gluon plasma's transport coefficients in previous model-to-data comparisons. However, sizable sensitivity is present in small p + Au collisions, which poses challenges to study the small systems' collectivity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hybrid simulation of Alfvén wave parametric decay instability in a laboratory relevant plasma

Large-amplitude Alfvén waves are subject to parametric decays, which can have important consequences in space, astrophysical, and fusion plasmas. Although the Alfvén wave parametric decay instability (PDI) was predicted decades ago, observational evidence is limited, stimulating considerable interest in laboratory demonstration of the instability and associated numerical modeling. Here, we report an important step toward direct hybrid simulation of the Alfvén wave PDI in a laboratory plasma, using antenna-like wave injection of a circularly polarized wave and realistic wave-plasma parameters. Considering collisionless damping, we identify the threshold Alfvén wave amplitudes and frequencies required for triggering the instability. These threshold behaviors are corroborated by simple theoretical analysis, considering the balance between PDI growth and ion Landau damping and PDI development in a bounded plasma. Other effects not included in the present model such as finite transverse wave scale and ion–neutral collisions are briefly discussed. These hybrid simulations promise to be a useful tool for investigating laboratory Alfvén wave dynamics and may provide guidance for future laboratory demonstration of the PDI process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Theoretical justification for heat flux limiter 0.15

In this section, we describe a semi-nonlocal kinetic model for the electron heat flux. The model has no free parameters. It takes the more physically motivated approach of placing a limiter on the perturbation to the electron distribution function, rather than the usual ad-hoc limiter on the heat flux. If the kinetic model is fitted by a heat flux limiter model, we find that the higher value of f e = 0.15 is a much closer fit than other models with f e ≃ 0.02.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hot QCD White Paper

Hot QCD physics studies the nuclear strong force under extreme temperature and densities. Experimentally these conditions are achieved via high-energy collisions of heavy ions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). In the past decade, a unique and substantial suite of data was collected at RHIC and the LHC, probing hydrodynamics at the nucleon scale, the temperature dependence of the transport properties of quark-gluon plasma, the phase diagram of nuclear matter, the interaction of quarks and gluons at different scales and much more. This document, as part of the 2023 nuclear science long range planning process, was written to review the progress in hot QCD since the 2015 Long Range Plan for Nuclear Science, as well as highlight the realization of previous recommendations, and present opportunities for the next decade, building on the accomplishments and investments made in theoretical developments and the construction of new detectors. Furthermore, this document provides additional context to support the recommendations voted on at the Joint Hot and Cold QCD Town Hall Meeting, which are reported in a separate document.

FOS: Physical sciences↗
NASA NTRS1981

Physics of auroral arc formation; Proceedings of the Chapman Conference on Formation of Auroral Arcs, Fairbanks, AK, July 21-25, 1980

The results of satellite, balloon-borne, and ground-based observations of auroral phenomena are discussed in terms of the morphology of auroral arcs and the behavior and characteristics of auroral electrons and ions. Further attention is given to auroral electric fields and field aligned currents, and to the development of models of auroral potential structures and energization of auroral particles. A simulation of space plasma phenomena is presented, along with numerical simulations of auroral potential structures and related phenomena, including the V-potential double layers and auroral arc deformations, and a simulation of auroral arcs. Plasma waves observed on auroral field lines and in laboratory conditions are reported, and theoretical studies of waves and turbulence in auroral plasmas are described.

Akasofu, S.-I.↗
DOE OSTI2024

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. Here, we develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The slow wave resonance cone in the collisional regime

In the low-density edge plasma of tokamaks, ion cyclotron range of frequencies actuators may parasitically emit slow waves. If the density is sufficiently low, which may be common in large future devices such as international thermonuclear experimental reactor (ITER), these slow waves take the form of so-called resonance cones. The traditional theoretical description of this wave mode relies on formally relating an electrostatic approximation of the frequency-domain wave equation to a time-domain wave equation and relating the cone angle to the wave speed in the time-domain wave equation. In the cold plasma collisional regime, that wave speed is complex. We investigate that scenario in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗