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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Plasma Exhaust in a Sheared-Flow-Stabilized Z Pinch

The fusion Z-pinch experiment device produces a 0.3 cm radius by 50 cm-long Z pinch between the end of the inner electrode of a coaxial plasma gun and an end wall 50 cm away. The plasma column is stabilized for thousands of instability growth times by an embedded radially sheared axial plasma flow. To investigate the effect of end wall design on Z-pinch plasma behavior, the solid end wall with a central exhaust port is replaced with a spoked end wall with a solid hub. While the Z-pinch plasma behavior was largely unaffected, the plasma exhausted from the Z-pinch provides an experimental platform to study magnetic confinement and detachment. Here, current and density signals are compared between two cases: a low energy and a high energy case. Plasma is found to be frozen-in flux in the Z-pinch assembly region. The exhaust of plasma from the flux conserving region is found to be dependent on the ratio of plasma ram and thermal pressures to the local magnetic pressure, calculated from an equilibrium model at the end wall. Plasma exhaust is observed to increase with the spoked end wall only for the high energy conditions as its ratio of ram and thermal pressures to the magnetic pressure increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Verification and benchmarking relativistic electron beam transport through a background gas

It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. Additionally, we also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas

A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.

Solmaz, Evrim (ORCID:0000000211830881)↗
DOE OSTI2026

Life cycle analysis of hydrogen production via methane pyrolysis using plasma arc

Steam methane reforming of natural gas is the primary method of producing hydrogen in the United States, accounting for 95% of all hydrogen produced there. Methane pyrolysis, an alternative production pathway that decomposes natural gas into solid carbon and hydrogen, both eliminates CO 2 emissions associated with methane reforming and allows for additional income from carbon black. A life-cycle inventory of this process has been developed using ASPEN Plus to model the methane pyrolysis (plasma arc) process. From well to gate, hydrogen production via methane pyrolysis produces 2.78 kg CO 2 e/kg H 2 of greenhouse gas emissions using mass allocation of emissions between hydrogen and carbon black coproducts. The well-to-gate emissions are mainly driven by electricity consumption (∼38 kW h/kg H 2 ), which accounts for 81% of the emissions; if renewable electricity is used, well-to-gate emissions can be reduced to −0.448 kg CO 2 e/kg H 2 .

08 HYDROGEN↗
DOE OSTI2022

Data-driven surrogate modeling of hPIC ion energy-angle distributions for high-dimensional sensitivity analysis of plasma parameters' uncertainty

In this work, we present a data-driven strategy for effective construction of a surrogate model in high-dimensional parameter space for the ion energy-angle distribution (IEAD) output of hPIC simulations of plasma-surface interactions. The methodology is based on a bin-by-bin least-squares fitting of the IEAD in the parameter space. The fitting is performed in a transformed coordinate system to normalize the IEAD, and it employs sparse grids for sampling the parameter space to overcome sampling challenges in high dimensions. The surrogate model is significantly cheaper computationally than direct hPIC simulations yet maintains high fidelity to them, providing a fast emulator for hPIC simulations. Sensitivity analysis based on the surrogate model is utilized to characterize the dependence of the ion impact angle and energy moments on the physical parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toward fusion plasma scenario planning for NSTX-U using machine-learning-accelerated models

One of the most promising devices for realizing power production through nuclear fusion is the tokamak. To maximize performance, it is preferable that tokamak reactors achieve advanced operating scenarios characterized by good plasma confinement, improved magnetohydrodynamic (MHD) stability, and a largely non-inductively driven plasma current. Such scenarios could enable steady-state reactor operation with high \emph{fusion gain} --- the ratio of produced fusion power to the external power provided through the plasma boundary. Precise and robust control of the evolution of the plasma boundary shape as well as the spatial distribution of the plasma current, density, temperature, and rotation will be essential to achieving and maintaining such scenarios. The complexity of the evolution of tokamak plasmas, arising due to nonlinearities and coupling between various parameters, motivates the use of model-based control algorithms that can account for the system dynamics. In this work, a learning-based accelerated model trained on data from the National Spherical Torus Experiment Upgrade (NSTX-U) is employed to develop planning and control strategies for regulating the density and temperature profile evolution around desired trajectories. The proposed model combines empirical scaling laws developed across multiple devices with neural networks trained on empirical data from NSTX-U and a database of first-principles-based computationally intensive simulations. The reduced execution time of the accelerated model will enable practical application of optimization algorithms and reinforcement learning approaches for scenario planning and control development. An initial demonstration of applying optimization approaches to the learning-based model is presented, including a strategy for mitigating the effect of leaving the finite validity range of the accelerated model. The approach shows promise for actuator planning between experiments and in real-time.

machine learning↗
DOE OSTI2022

Investigation of resistive magnetic field generation by intense proton beams in dense plasmas

Current and future applications of intense proton sources abound, including radiography, cancer therapy, warm dense matter generation, and inertial confinement fusion. With increasingly efficient acceleration and focusing mechanisms, proton current densities may soon approach and exceed 10 10 A/cm 2 , e.g., via intense laser drivers. Simulations have previously shown that in this current density regime, beam-induced field generation plays a significant role in beam transport through dense plasmas. Here, we present a theoretical model for the generation of resistive magnetic fields by intense proton beam transport through solid density plasmas. The theoretical evolution of the magnetic field profile is calculated using an analytic model for aluminum resistivity, heat capacity, and stopping power, applicable from cold matter to hot plasma. The effects of various beam and material parameters on the field are investigated and explained for both monoenergetic and Maxwellian proton beams. For a proton beam with Maxwellian temperature 5 MeV and total energy 10 J, the model calculates resistive magnetic fields up to 150 T in aluminum. In conclusion, the calculated field profiles from several beam cases are compared with 2D hybrid particle-in-cell simulations, with good agreement found in magnitude and time scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experiment and model for a Stokes layer in a strongly coupled dusty plasma

A Stokes layer, which is a flow pattern that arises in a viscous fluid adjacent to an oscillatory boundary, was observed in an experiment using a two-dimensional strongly coupled dusty plasma. Liquid conditions were maintained using laser heating, while a separate laser manipulation applied an oscillatory shear that was localized and sinusoidal. The evolution of the resulting flow was analyzed using space-time diagrams. These figures provide an intuitive visualization of a Stokes layer, including features such as the depth of penetration and wavelength. Another feature, the characteristic speed for the penetration of the oscillatory flow, also appears prominently in space-time diagrams. To model the experiment, the Maxwell-fluid model of a Stokes layer was generalized to describe a two-phase liquid. In our experiment, the phases were gas and dust, where the dust cloud was viscoelastic due to strong Coulomb coupling. Furthermore, the model is found to agree with the experiment, in the appearance of the space-time diagrams, and in the values of the characteristic speed, depth of penetration, and wavelength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The influence of surface impurities on photoelectric currents driven by intense soft x rays

In an x-ray driven cavity experiment, an intense flux of soft x rays on the emitting surface produces significant emission of photoelectrons having several kiloelectronvolts of kinetic energy. At the same time, rapid heating of the emitting surface occurs, resulting in the release of adsorbed surface impurities and subsequent formation of an impurity plasma. This numerical study explores a simple model for the photoelectric currents and the impurity plasma. In this work, attention is given to the effect of varying the composition of the impurity plasma. The presence of protons or hydrogen molecular ions leads to a substantially enhanced cavity current, while heavier plasma ions are seen to have a limited effect on the cavity current due to their lower mobility. Additionally, it is demonstrated that an additional peak in the current waveform can appear due to the impurity plasma. A correlation between the impurity plasma composition and the timing of this peak is elucidated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of ionization waves on dust chain formation in a DC discharge

An interesting aspect of complex plasma is its ability to self-organize into a variety of structural configurations and undergo transitions between these states. A striking phenomenon is the isotropic-to-string transition observed in electrorheological complex plasma under the influence of a symmetric ion wake field. Such transitions have been investigated using the Plasma Kristall-4 (PK-4) microgravity laboratory on the International Space Station. Recent experiments and numerical simulations have shown that, under PK-4-relevant discharge conditions, the seemingly homogeneous direct current discharge column is highly inhomogeneous, with large axial electric field oscillations associated with ionization waves occurring on microsecond time scales. here, a multi-scale numerical model of the dust–plasma interactions is employed to investigate the role of the electric field in the charge of individual dust grains, the ion wake field and the order of string-like structures. Results are compared with those for dust strings formed in similar conditions in the PK-4 experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Two Mechanisms Limiting the Emitted Electron Current from a Cathode to an Anode

It is known that the current of emitted electrons flowing through a plasma can saturate upon formation of a potential well adjacent to the cathode (the “space charge effect”). Here, we demonstrate another saturation mechanism that will often set a more restrictive limit on the global current. When “backflow saturation” occurs, the cathode sheath weakens to allow emitted electrons that already entered the plasma to backflow to the cathode. This effect could not be captured by studies modelling the cathode sheath by itself because its origin is coupled to processes in the interior plasma and anode sheath. By modeling a full plasma diode, we show that depending on conditions the global current can be limited in four ways; by backflow alone, by space charge alone, by both mechanisms in a stable cooperative form, or by both in a competing oscillatory form. Published by the American Physical Society 2025

Laboratory plasma↗
DOE OSTI2022

Faraday-effect polarimetry for current profile measurement in the tokamak plasma edge

Toroidal current profile measurements in the tokamak plasma edge are critical for fusion plasma physics research and model validation. A three-wave Faraday-effect polarimeter-interferometer with a sub-centimeter spatial resolution is proposed on the DIII-D tokamak to determine the edge current profile via Abel inversion. By using probe beams with 316 µm wavelength, a low-field-side, vertical-view, single-pass optical layout covering the plasma edge region (R = 2.15–2.27 m) is assessed. Measurements with no greater than 0.1° polarimetric systematic uncertainty, no greater than 0.01° polarimetric root-mean-square noise (1 kHz bandwidth), and a 0.8 cm radial chord spacing are considered feasible based on the achieved performance of existing systems using similar wavelengths on fusion devices. Furthermore, synthetic diagnostic calculations taking various factors into account, such as diagnostic uncertainty and quality of magnetic flux surfaces, find that the edge current profile can be determined with up to 0.12 MA/m 2 uncertainty, or about 10% of the peak current density in the pedestal of an investigated high-confinement plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Surrogate Modeling of Landau Damping with Deep Operator Networks

Kinetic simulations excel at capturing microscale plasma physics phenomena with high accuracy, but their computational demands make them impractical for modeling large-scale space and astrophysical systems. In this context, we build a surrogate model, using Deep Operator Networks (DeepONets), based upon the Vlasov–Poisson simulation data to model the dynamical evolution of plasmas, focusing on the Landau damping process—a fundamental kinetic phenomenon in space and astrophysical plasmas. The trained DeepONets are able to capture the evolution of electric field energy in both linear and nonlinear regimes under various conditions. Extensive validation highlights DeepONets’ robust performance in reproducing complex plasma behaviors with high accuracy, paving the way for large-scale modeling of space and astrophysical plasmas.

plasma astrophysics↗
DOE OSTI2021

Modeling of thermal and kinetic processes in non-equilibrium plasma ignition applied to a lean combustion engine

In recent years novel ignition systems have been developed to enable stable and efficient engine operations with lean mixtures. Among them, radio-frequency corona ignition systems create discharges that involve a much wider region compared to traditional spark, and produce non equilibrium plasma with high levels of active radicals and excited species. These devices considerably increase the early flame growth speed and extend stable operating limits. With the aim of expanding the knowledge on high efficiency lean-burn SI engines, this paper investigates and compares the combustion development generated by spark and corona ignitions through computational fluid dynamics, within the Reynolds-Averaged Navier-Stokes framework for turbulence modeling. In order to simultaneously take thermal and chemical effects into account, the Perfectly Stirred Reactor combustion model is used. Experimental data are also collected for validation in an optical access engine, for different mixture levels, from stoichiometric to very lean. Furthermore, the faster burn rate generated by the corona system in the initial stage of the combustion is well predicted by the simulations, in all the relative air-fuel ratio conditions. Remarkably, as the mixture becomes lean, simulations are able to capture the non-linear transition from fast to slow kernel growth, before a self-sustainable flame propagation is established. This correlates very well with the measured engine cyclic variability and the corresponding steep change in the duration of the flame kernel formation. Ultimately, this study highlights the important role of the atomic oxygen, as active radical, in promoting and enhancing the combustion initiated by a corona discharge, in addition to the volumetric ignition effect. By contrast, the validated simulations allow to explain that the high-temperature thermal plasma generated in a traditional spark discharge is insensitive to kinetic aspects.

42 ENGINEERING↗
DOE OSTI2020

The effect of enrofloxacin on enteric Escherichia coli : Fitting a mathematical model to in vivo data

Antimicrobial drugs administered systemically may cause the emergence and dissemination of antimicrobial resistance among enteric bacteria. To develop logical, research-based recommendations for food animal veterinarians, we must understand how to maximize antimicrobial drug efficacy while minimizing risk of antimicrobial resistance. Our objective is to evaluate the effect of two approved dosing regimens of enrofloxacin (a single high dose or three low doses) on Escherichia coli in cattle. We look specifically at bacteria above and below the epidemiological cutoff (ECOFF), above which the bacteria are likely to have an acquired or mutational resistance to enrofloxacin. We developed a differential equation model for the antimicrobial drug concentrations in plasma and colon, and bacteria populations in the feces. The model was fit to animal data of drug concentrations in the plasma and colon obtained using ultrafiltration probes. Fecal E. coli counts and minimum inhibitory concentrations were measured for the week after receiving the antimicrobial drug. We predict that the antimicrobial susceptibility of the bacteria above the ECOFF pre-treatment strongly affects the composition of the bacteria following treatment. Faster removal of the antimicrobial drugs from the colon throughout the study leads to improved clearance of bacteria above the ECOFF in the low dose regimen. If we assume a fitness cost is associated with bacteria above the ECOFF, the increased fitness costs leads to reduction of bacteria above the ECOFF in the low dose study. These results suggest the initial E. coli susceptibility is a strong indicator of how steers respond to antimicrobial drug treatment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Self-consistent simulation of resistive kink instabilities with runaway electrons

A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Here, given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evaluating the Impact of Tritium Permeation Membrane Performance and Direct Internal Recycling on Fusion Fuel Cycle Efficiency Using TMAP8

An efficient fuel cycle is vital to sustainable and cost-effective energy generation in fusion systems. Since tritium is not widely available, fusion systems must breed their own tritium for sustainable fusion deuterium-tritium reactions. An inefficient fuel cycle increases the tritium inventory needed for operations, which increases costs, constraints on tritium management systems, and safety concerns. A fuel cycle model is a powerful tool for understanding tritium inventories and flow rates across all systems in the fuel cycle. By simplifying the technical details into time-dependent tritium flow rates and inventories, the model can simulate the entire fuel cycle with high computational efficiency, even for technologies that are still under development. It can therefore quantify the impact of new tritium management technologies on fuel cycle efficiency. To evaluate the impact of key components on reducing tritium inventory, we are using and expanding an existing fuel cycle models based on latest advancements in fuel cycle research. The new model integrates Tritium Permeation Membrane (TPM) and Direct Internal Recycling (DIR) to enhance tritium transport from blanket breeders and plasma exhaust. These fuel cycle models are implemented in TMAP8 (Tritium Migration Analysis Program, version 8), a MOOSE-based open-source application designed to provide cutting-edge capabilities for tritium transport and fuel cycle modeling. The study aims to demonstrate the extensibility of existing fuel cycle modeling capability in TMAP8 and to offer a proof-of-principle design for future fusion plant systems. The presentation will cover the performance of fuel cycle modeling capabilities available in TMAP8, highlight advancements in fuel cycle research, and present a sensitivity analysis of these models. The results underline potential approaches and technology solutions to lower tritium inventory requirements, highlighting their role in shaping the future of fusion energy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Understanding ICF hohlraums using NIF gated laser-entrance-hole images

The newly available ns-gated laser-entrance-hole (LEH) imager on the National Ignition Facility provides routine, non-perturbative measurements of the x-ray emission from laser-heated plasmas inside the hohlraum as viewed at 19° to the hohlraum axis through one of its LEHs. Multiple images are acquired for a series of times and filter-selected x-ray energy bands within a single shot. The images provide time dependent data on phenomena including the effective radius of the LEH, the length of the gold-plasma “bubble” evolving off the interior wall surface heated by the outer beams, the evolving radius of the x-ray heated hohlraum wall, and the radius of the ablation front of the fusion capsule. Here, these measurements are explained and illustrated with sample data. These techniques are then applied to understand hohlraum behavior as a function of gas fill. For hohlraums with helium gas fill densities of 0.15 to 0.30 mg/cm 3 , synthetic images computed from simulations agree well with experimental gated LEH images when an inhibited heat transport model [Jones et al., Phys. Plasmas 24, 056312 (2017)] is used. This model can be adjusted to reproduce the expansion rate of the laser-heated plasma bubble in such a way as to improve agreement with the images. At the higher 0.6 mg/cc gas fill, the experimental images show more pronounced 3D features, resulting in slightly less good agreement with the 2D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗