SEARCH · Search NASA
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.
Thin film and semiconductor microelectronics, radar scattering, radome thermal stress, boundary layer phenomena, guided missile parts, turbulent mixing, antenna systems, and plasma dynamics
An experimental plasma accelerator for magnetic target fusion (MTF) applications under development at the NASA Marshall Space Flight Center is described. The accelerator is a pulsed plasma thruster and has been tested experimentally and plasma jet velocities of approximately 50 km/sec have been obtained. The plasma jet structure has been photographed with 10 ns exposure times to reveal a stable and repeatable plasma structure. Data for velocity profile information has been obtained using light pipes embedded in the gun walls to record the plasma transit at various barrel locations. Preliminary spatially resolved spectral data and magnetic field probe data are also presented. A high speed triggering system has been developed and tested as a means of reducing the gun "jitter". This jitter is being characterized and future work for second generation "ultra-low jitter" gun development is being identified.
An experimental plasma accelerator for magnetic target fusion (MTF) applications under development at the NASA Marshall Space Flight Center is described. The accelerator is a coaxial pulsed plasma thruster (Figure 1). It has been tested experimentally and plasma jet velocities of approx.50 km/sec have been obtained. The plasma jet has been photographed with 10-ns exposure times to reveal a stable and repeatable plasma structure (Figure 2). Data for velocity profile information has been obtained using light pipes and magnetic probes embedded in the gun walls to record the plasma and current transit respectively at various barrel locations. Preliminary spatially resolved spectral data and magnetic field probe data are also presented. A high speed triggering system has been developed and tested as a means of reducing the gun "jitter". This jitter is being characterized and future work for second generation "ultra-low jitter" gun development is being identified.
A conceptual model of the formation of the plasma sheet and of its dynamical behavior in association with magnetospheric substorms is proposed. The two essential assumptions of the model are that (1) the plasma mantle is the primary source of plasma-sheet particles and (2) the momentum tends to inhibit the rate of magnetic field annihilation (merging) in the tail current sheet. It is found that there exists no steady state configuration for the plasma sheet, which must instead shrink continuously in thickness until the near-earth portion of the current sheet is disrupted by the formation of a magnetic neutral line. The thinning of the plasma sheet and the resulting current-sheet disruption are proposed to be the direct causes of the growth phase and expansion phase, respectively, of substorms.
Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.
We have imaged the evolution of plasma following nanosecond pulsed laser ablation of cerium in air with sub-10 ns resolution. We find the ablated plasma plume expansion and shape to be strongly dependent on the time after mechanical polishing of the sample surface. After an initial shock front common to both new and aged samples (velocities up to 16 km s ₋1 ), we observe a second front that is relatively localized to the surface vicinity of newly polished samples, but in the case of older samples it moves outward at velocities of up to 4 km s ₋1 . We attribute this behavior to the formation and growth of a native oxide layer on the order of hours after polishing. These results demonstrate that plasma imaging can be used to diagnose the presence of, and evaluate the extent of, thin surface oxides or other compounds that form soon after exposure to some atmosphere. These findings also highlight the need to explicitly consider the presence of such overlayers when modeling laser/metal interaction, when making measurements of vapor phase chemistry, and when analyzing the composition, phase, and morphology of solid particulates formed after the ablation of reactive metals.
Due to the non-linear nature of relativistic laser induced plasma processes, the development of laser-plasma accelerators requires precise numerical modeling. Especially high intensity laser-solid interactions are sensitive to the temporal laser rising edge and the predictive capability of simulations suffers from incomplete information on the plasma state at the onset of the relativistic interaction. Experimental diagnostics utilizing ultra-fast optical backlighters can help to ease this challenge by providing temporally resolved inside into the plasma density evolution. We present the successful implementation of an off-harmonic optical probe laser setup to investigate the interaction of a high-intensity laser at 5.4 x 10 21 W/cm 2 peak intensity with a solid-density cylindrical cryogenic hydrogen jet target of 5μm diameter as a target test bed. The temporal synchronization of pump and probe laser, spectral filtering and spectrally resolved data of the parasitic plasma self-emission are discussed. The probing technique mitigates detector saturation by self-emission and allowed to record a temporal scan of shadowgraphy data revealing details of the target ionization and expansion dynamics that were so far not accessible for the given laser intensity. Plasma expansion speeds of up to (2.3 ± 0.4)x 10 7 m/s followed by full target transparency at 1.4ps after the high intensity laser peak are observed. A three dimensional particle-in-cell simulation initiated with the diagnosed target pre-expansion at -0.2ps and post processed by ray tracing simulations supports the experimental observations and demonstrates the capability of time resolved optical diagnostics to provide quantitative input and feedback to the numerical treatment within the time frame of the relativistic laser-plasma interaction.
Plasma wakefield dynamics over timescales up to 800 ps, approximately 100 plasma periods, are studied experimentally at the Advanced Wakefield Experiment (AWAKE). The development of the longitudinal wakefield amplitude driven by a self-modulated proton bunch is measured using the external injection of witness electrons that sample the fields. In simulation, resonant excitation of the wakefield causes plasma electron trajectory crossing, resulting in the development of a potential outside the plasma boundary as electrons are transversely ejected. Trends consistent with the presence of this potential are experimentally measured and their dependence on wakefield amplitude are studied via seed laser timing scans and electron injection delay scans.
The stability of the weakly ionized plasma and the transition from a stable homogeneous discharge to unstable filaments play an important role in gas laser physics, plasma-assisted combustion, chemical reforming, and material synthesis. Here, theoretical stability analysis and thermal-chemical mode analysis were performed to understand the mechanism of plasma thermal-chemical instability by using a zero-dimensional plasma system with both simplified and detailed chemical kinetics of H 2 /O 2 /N 2 mixtures. The plasma dynamic and kinetic models accounted for multiple physical mechanisms in the chemically-reactive weakly ionized plasma, including ionization, attachment/detachment, recombination, vibrational and electronic energy relaxation, convective and diffusive species/heat removal, Joule heating, and detailed chemical kinetics. An analytical criterion and the explosive mode species/temperature pointers were formulated while the representative active species were identified for different thermal-chemical modes. The results showed that in addition to the classical thermal-ionization mechanism, various chemical modes from chemical heat imbalance and elementary kinetics significantly modified the time dynamics and the stability of the weakly ionized plasma. The present analysis provides insights and guidance to control plasma instability using chemical kinetics.
Theories which describe currents collected by conducting and non-conducting bodies immersed in plasmas have many of their concepts based upon the fundamentals of sheath-potential distributions and charged-particle behavior in superimposed electric and magnetic fields. Those current-collecting bodies (or electrodes) may be Langmuir probes, electric field detectors, aperture plates on ion mass spectrometers and retarding potential analyzers, or spacecraft and their rigid and tethered appendages. Often the models are incomplete in representing the conditions under which the current-voltage characteristics of the electrode and its system are to be measured. In such cases, the experimenter must carefully take into account magnetic field effects and particle anisotropies, perturbations caused by the current collection process itself and contamination on electrode surfaces, the complexities of non-Maxwellian plasma distributions, and the temporal variability of the local plasma density, temperature, composition and fields. This set of variables is by no means all-inclusive, but it represents a collection of circumstances guaranteed to accompany experiments involving energetic particle beams, plasma discharges, chemical releases, wave injection and various events of controlled and uncontrolled spacecraft charging. Here, an attempt is made to synopsize these diagnostic challenges and frame them within a perspective that focuses on the physics under investigation and the requirements on the parameters to be measured. Examples include laboratory and spaceborne applications, with specific interest in dynamic and unstable plasma environments.
The composition of exhausted gas is a key parameter in long-pulse plasma fusion experiments, and its evolution shall be monitored at timescales relevant to plasma dynamics and plasma-wall interactions. A diagnostic residual gas analyzer (DRGA) is a multisensor instrument particularly suited to these studies, and ITER will adopt DRGAs in the equatorial and in the divertor tokamak regions. In this work, we have revisited the design of the ITER divertor DRGA through simple vacuum analytical considerations supported by simulations conducted with Molflow+, a test particle Monte Carlo (TPMC) simulation code commonly used in the particle accelerator community. Here, starting with recommendations on the manufacturing of the vacuum piping of the DRGA, this work is followed by a complete vacuum characterization of the diagnostic vacuum setup (pressure profiles at base pressure and during sampling, orifice diameter, and length optimization), and finally, the in-vessel residence time of the most important gas species is simulated. These studies have allowed us to give insights into some experimental results recently found on the prototype DRGA installed in the Wendelstein W7-X stellarator.
Here, a neural network version of the Multi-Mode Anomalous Transport Module, known as MMMnet, has been developed to calculate plasma turbulent diffusivities in DIII-D with a calculation time suitable for control applications. MMMnet uses a simple artificial neural network structure to predict the ion thermal, electron thermal, and toroidal momentum diffusivities while reproducing Multi-Mode Model (MMM) data with good accuracy and keeping the calculation time as a fraction of that associated with MMM. Model-based control techniques require models with fast calculation times, making many existing physics-oriented predictive codes unsuitable. The control-oriented predictive code COTSIM (Control Oriented Transport SIMulator) calculates the most significant plasma dynamics in response to the different actuators while running at a speed useful for control design. In order to achieve this calculation speed, COTSIM often relies on scaling laws and control-level models. Replacing some of these scaling laws and control-level models with neural network versions of more complex physics-level models has the potential of increasing the range of validity and the level of accuracy of COTSIM without compromising its computational speed. In this work, MMMnet is integrated into COTSIM to improve the turbulent diffusivity predictions, which will in turn improve the prediction accuracy associated with the dynamics of many plasma properties.
With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.
Electron density images and plasma dynamics in the middle and high latitude dayside magnetosphere are studied with the IMAGE/RPI remote measurements of the electron densities along magnetic field lines, measured before and during a magnetic storm when the solar wind and interplanetary magnetic field (IMF) impinging on the magnetopause varied considerably. Several regions of different density distribution characteristics, including plasmasphere, plasma trough, sub-auroral/auroral density depletion, density enhancements in the aurora/cusp, and polar cap, are identified in 'two dimensional images', i.e., along the satellite orbit and field lines. The plasma dynamics, such as the plasma refilling in the outer plasmasphere and the plasma acceleration in the aurora/cusp region are inferred from density gradients along the field lines. It is shown that the densities and locations of the plasma regions vary in accordance with the solar wind, particularly with the IMF variations for the case examined. The changes in these regions reflect the way the magnetosphere reconfigures in response to changes in the solar wind.
The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.
Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) MK-1 pulsed plasma thruster. Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.
Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) pulsed electromagnetic plasma accelerator (PEPA-0). Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.
Two phases of substorm-associated magnetospheric dynamics are discussed in terms of particles and fields at synchronous orbit. The first phase corresponds to the 'decreases' of energetic particle flux. The second phase is the conventional expansion phase that begins with the 'onset', characterized in the study by (1) a sudden decrease in the tail current and a return of the inflated magnetosphere to a dipolelike configuration, (2) a sudden shift of trapped high-energy particles toward the tail again following contours of constant B, and at the same time (3) a surge of tail plasma toward the earth as the induced electric field now increases the total convection field.