Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA SIMULATION”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Full wave simulation of waves in ECRIS plasmas based on the finite element method

This paper describes the modeling and the full wave numerical simulation of electromagnetic waves propagation and absorption in an anisotropic magnetized plasma filling the resonant cavity of an electron cyclotron resonance ion source (ECRIS). The model assumes inhomogeneous, dispersive and tensorial constitutive relations. Maxwell's equations are solved by the finite element method (FEM), using the COMSOL Multiphysics{sup ®} suite. All the relevant details have been considered in the model, including the non uniform external magnetostatic field used for plasma confinement, the local electron density profile resulting in the full-3D non uniform magnetized plasma complex dielectric tensor. The more accurate plasma simulations clearly show the importance of cavity effect on wave propagation and the effects of a resonant surface. These studies are the pillars for an improved ECRIS plasma modeling, that is mandatory to optimize the ion source output (beam intensity distribution and charge state, especially). Any new project concerning the advanced ECRIS design will take benefit by an adequate modeling of self-consistent wave absorption simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Simulation of the Plasma Afterglow in the Discharge Gap of a Subnanosecond Switch Based on an Open Discharge in Helium

The phenomenon of subnanosecond electrical breakdown in a strong electric field observed in an open discharge in helium at pressures of 6–20 Torr can be used to create ultrafast plasma switches triggering into a conducting state for a time shorter than 1 ns. To evaluate the possible repetition rate of such a subnanosecond switch, it is interesting to study the decay dynamics of the plasma remaining in the discharge gap after ultrafast breakdown. In this paper, a kinetic model based on the particle-in-cell Monte Carlo collision method is used to study the dynamics of the plasma afterglow in the discharge gap of a subnanosecond switch operating with helium at a pressure of 6 Torr. The simulation results show that the radiative, collisional-radiative, and three-body collision recombination mechanisms significantly contribute to the afterglow decay only while the plasma density remains higher than 10{sup 12} cm{sup −3}; the main mechanism of the further plasma decay is diffusion of plasma particles onto the wall. Therefore, the effect of recombination in the plasma bulk is observed only during the first 10–20 μs of the afterglow. Over nearly the same time, plasma electrons become thermalized. The afterglow time can be substantially reduced by applying a positive voltage U{sub c} to the cathode. Since diffusive losses are limited by the ion mobility, the additional ion drift toward the wall significantly accelerates plasma decay. As U{sub c} increases from 0 to +500 V, the characteristic time of plasma decay is reduced from 35 to 10 μs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

High-resolution hybrid simulations of kinetic plasma turbulence at proton scales

We investigate properties of plasma turbulence from magnetohydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe an MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of −5/3 and −3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ∼−3. A −2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell, and the resolution employed, the perpendicular proton temperature strongly depends on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Numerical error analysis of SOLPS-ITER simulations of EAST

Abstract Plasma edge simulations with codes like SOLPS-ITER are widely employed to interpret fusion experiments. However, numerical errors appearing in such simulations are rarely investigated, despite their potential large impact on simulation results. These errors consist of the statistical error and the bias, both resulting from the finite number of employed EIRENE Monte Carlo particles and incomplete convergence, and the discretization error due to the finite resolution of the computational grids. In this contribution, the resulting numerical errors on simulations of pure deuterium and neon seeded H-mode EAST discharges are examined. The statistical error can be kept small compared to other numerical error contributions by averaging the plasma profiles. This allows investigating the bias and discretization errors using Richardson extrapolation. It is shown that grid refinement and the number of employed Monte Carlo particles have the largest influence on the result, in agreement with similar studies of an ITER deuterium case. For the first time, numerical error bars on the entire simulated target profiles are determined showing that the largest numerical error is 17.9%, mainly due to the plasma grid discretization. On top, also numerical errors on simulated neutral pressures are investigated in detail, for which the statistical error is dominant. The analysis demonstrates which setup is needed to keep numerical errors limited: the SOLPS-ITER averaging procedure should be employed including enough EIRENE particles, and the involved grids should be sufficiently refined to reduce discretization errors.

Boeyaert, Dieter (ORCID:0000000309208660)↗
NASA NTRS1983

Numerical simulations of counterstreaming plasmas and their relevance to interhemispheric flows

The collisionless expansion of counterstreaming plasmas has been studied in order to elucidate the basic physical processes that may be operating during the initial refilling of depleted flux tubes after a magnetic storm. The numerical technique applied is briefly described, and simulation results are presented. The simulation geometry consisted of two high-density H(+)-O(+) electron plasmas separated by a low-density H(+)-electron plasma. The temporal evolution of the expanding plasmas and the electostatic potential in the region between the two sources is described. The main interacting streams are found to be stable with respect to both the ion acoustic and ion cyclotron modes and only the suprathermal forerunner ions are unstable with respect to the ion cyclotron mode. The results also suggest that a localized potential hill can form at the equator and that this potential hill can play an important role in the subsequent trapping and thermalization of the ion streams.

Singh, N.↗
NASA NTRS1976

Plasma core reactor simulations using RF uranium seeded argon discharges

Experimental results are described in which pure uranium hexafluoride was injected into an argon-confined, steady-state, RF-heated plasma to investigate characteristics of plasma core nuclear reactors. The 80 kW (13.56 MHz) and 1.2 MW (5.51 MHz) rf induction heater facilities were used to determine a test chamber flow scheme which offered best uranium confinement with minimum wall coating. The cylindrical fused-silica test chamber walls were 5.7-cm-ID by 10-cm-long. Test conditions included RF powers of 2-85 kW, chamber pressures of 1-12 atm, and uranium hexafluoride mass-flow rates of 0.005-0.13 g/s. Successful techniques were developed for fluid-mechanical confinement of RF-heated plasmas with pure uranium hexafluoride injection.

Roman, W. C.↗
DOE OSTI2025

Quantum Information for Fusion Energy Sciences (Final Technical Report)

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Surface wave and linear operating mode of a plasma antenna

The relation between the propagation conditions of a surface electromagnetic wave along a finiteradius plasma cylinder and the linear operating mode of a plasma antenna is investigated. The solution to the dispersion relation for a surface wave propagating along a finite-radius plasma cylinder is analyzed for weakly and strongly collisional plasmas. Computer simulations of an asymmetrical plasma dipole antenna are performed using the KARAT code, wherein the dielectric properties of plasma are described in terms of the Drude model. The plasma parameters corresponding to the linear operating mode of a plasma antenna are determined. It is demonstrated that the characteristics of the plasma antenna in this mode are close to those of an analogous metal antenna.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Plasma Energetic Particles Simulation Center (PEPSC)

The main effort of the Texas group was to develop theoretical and simplified numerical models to understand chirping phenomena often seen for Alfven and geodesic acoustic waves in experimental plasmas such as D-III-D, NSTX and JET. Its main numerical effort was to modify the AEGIS code, which was originally developed as an eigenvalue solver. To apply to the chirping problem this code has to be able to treat the linear response to the continuum and the response of the plasma to external drive or to an internal drive that comes from the formation of phase space chirping structures. The theoretical underpinning of this investigation still needed to be more fully developed to understand how to best formulate the theoretical problem. Considerable progress was made on this front by B.N. Breizman and his collaborators and a new reduced model was developed by H. L. Berk and his PhD student, G. Wang which can be uses as simplified model to describe chirping in a large aspect ratio tokamak. This final report will concentrate on these two directions that were developed as well as results that were found in the work with the AEGIS code and in the progress in developing a novel quasi-linear formulation for a description of Alfvenic modes destabilized by energetic particles, such as alpha particles in a burning plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Three-dimensional simulation study of ionospheric plasma clouds

The results of fully three-dimensional numerical simulations of ionospheric plasma cloud evolution are presented. The evolution of the plasma cloud considered by Drake and Huba (1987) in the limit of vanishingly small ion compressibility is discussed. Simulations support the results of the analytical theory: finite plasma temperature, combined with fully three-dimensional plasma dynamics, is a stabilizing influence on plasma cloud evolution. This stability is associated with sheared azimuthal ion flows in the vicinity of the cloud surface. Cloud evolution using realistic values of ion compressibility show that the cloud rapidly diffuses to a state in which the sheared azimuthal flow is substantially reduced; subsequently, the cloud becomes unstable and structures.

Zalesak, S. T.↗
NASA NTRS1972

Probe measurements of a cesium plasma in a simulated thermionic energy converter

Cesium-filled thermionic energy converters are considered as electrical energy sources in future spacecraft requiring tens to hundreds of kilowatts of electric power. The high operating temperatures necessary for a large specific power and high efficiency impose stringent constraints on the converter fabrication. The converter physics for reducing operating temperatures and cesium plasma losses are being studied to achieve high reliability without sacrificing the power performance of the converters. Various cesium parameters which affect the converter performance are: (1) electron temperatures, (2) plasma ion densities, and (3) electric potential profiles. These were investigated using a Langmuir probe in a simulated converter. The parameters were measured in different cesium discharge modes.

Shimada, K.↗
DOE OSTI2018

NIMROD Development AND Applications for Advanced Simulations of Tokamak Plasmas

Operation in a regime of zero off-normal events is an important goal for ITER, and after that DEMO. Edge Localized Modes (ELMs) are off-normal events of particular concern due to their ubiquity in the highest performing discharges and their localized energy deposition on vessel walls. To avoid damage from these events, ITER has invested in costly coils in addition to exploring other mitigation techniques. A more attractive approach is to find an operational regime that has the requisite high performance while avoiding ELMs entirely. DIII-D experimentalists have found such a regime, wherein the ELMs are prevented from occurring by the existence of a saturated instability. Although this regime is attractive, how this saturated instability occurs and how to extrapolate to ITER is not well understood. NIMROD is well suited for studying this mode, and significant progress has already been made in modeling the equilibria that give rise to this instability. Another off-normal event of concern to ITER is the sawtooth mode in the presence of energetic particles. These modes are ubiquitous in near-steady-state tokamaks, and their behavior in the presence of the energetic particles that will be found in burning plasma experiments is poorly understood. The internal kink is notorious for requiring all extended MHD effects (two-fluid, closures, vacuum regime, energetic particles) in a production level simulation. NIMROD has matured to the point where it can now address this problem, though significant expertise in the MHD physics, numerical algorithms, and large-scale computation will be required. The expertise of Tech-X in attacking the most difficult NIMROD problems has already resulted in significant improvements to the NIMROD code that benefit the larger fusion community, and we expect further improvements as we attack these challenging problems. To maximize the community benefits, Tech-X provided community outreach and code management efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Plasma core reactor simulations using RF uranium seeded argon discharges

An experimental investigation was conducted using the United Technologies Research Center (UTRC) 80 kW and 1.2 MW RF induction heater systems to aid in developing the technology necessary for designing a self-critical fissioning uranium plasma core reactor (PCR). A nonfissioning, steady-state RF-heated argon plasma seeded with pure uranium hexafluoride (UF6) was used. An overall objective was to achieve maximum confinement of uranium vapor within the plasma while simultaneously minimizing the uranium compound wall deposition. Exploratory tests were conducted using the 80 kW RF induction heater with the test chamber at approximately atmospheric pressure and discharge power levels on the order of 10 kW. Four different test chamber flow configurations were tested to permit selection of the configuration offering the best confinement characteristics for subsequent tests at higher pressure and power in the 1.2 MW RF induction heater facility.

Roman, W. C.↗
DOE OSTI2013

Dense electron-positron plasmas and bursts of gamma-rays from laser-generated quantum electrodynamic plasmas

In simulations of a 12.5 PW laser (focussed intensity I=4×10{sup 23}Wcm{sup −2}) striking a solid aluminum target, 10% of the laser energy is converted to gamma-rays. A dense electron-positron plasma is generated with a maximum density of 10{sup 26}m{sup −3}, seven orders of magnitude denser than pure e{sup −} e{sup +} plasmas generated with 1PW lasers. When the laser power is increased to 320 PW (I=10{sup 25}Wcm{sup −2}), 40% of the laser energy is converted to gamma-ray photons and 10% to electron-positron pairs. In both cases, there is strong feedback between the QED emission processes and the plasma physics, the defining feature of the new “QED-plasma” regime reached in these interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Simulation of direct plasma injection for laser ion beam acceleration with a radio frequency quadrupole

The direct plasma injection scheme (DPIS) has been being studied at Institute of Modern Physics since several years ago. A C{sup 6+} beam with peak current of 13 mA, energy of 593 keV/u has been successfully achieved after acceleration with DPIS method. To understand the process of DPIS, some simulations have been done as follows. First, with the total current intensity and the relative yields of different charge states for carbon ions measured at the different distance from the target, the absolute current intensities and time-dependences for different charge states are scaled to the exit of the laser ion source in the DPIS. Then with these derived values as the input parameters, the extraction of carbon beam from the laser ion source to the radio frequency quadrupole with DPIS is simulated, which is well agreed with the experiment results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2005

Improved Simulations of Astrophysical Plasmas: Computation of New Atomic Data

Our research program is designed to carry out state-of-the-art atomic physics calculations crucial to advancing our understanding of fundamental astrophysical problems. We redress the present inadequacies in the atomic data base along two important areas: dielectronic recombination and inner-shell photoionization and multiple electron ejection/Auger fluorescence therefrom. All of these data are disseminated to the astrophysical community in the proper format for implementation in spectral simulation code.

Gorczyca, Thomas W.↗
NASA NTRS1993

Simulations of a molecular plasma in collisional-radiative nonequilibrium

A code for the simulation of nonequilibrium plasmas is being developed, with the capability to couple the plasma fluid-dynamics for a single fluid with a collisional-radiative model, where electronic states are treated as separate species. The model allows for non-Boltzmann distribution of the electronic states. Deviations from the Boltzmann distributions are expected to occur in the rapidly ionizing regime behind a strong shock or in the recombining regime during a fast expansion. This additional step in modeling complexity is expected to yield more accurate predictions of the nonequilibrium state and the radiation spectrum and intensity. An attempt at extending the code to molecular plasma flows is presented. The numerical techniques used, the thermochemical model, and the results of some numerical tests are described.

Cambier, Jean-Luc↗
DOE OSTI2017

Partnership for Edge Physics (EPSI), University of Texas Final Report

Simulations of tokamak plasmas require a number of inputs whose values are uncertain. The effects of these input uncertainties on the reliability of model predictions is of great importance when validating predictions by comparison to experimental observations, and when using the predictions for design and operation of devices. However, high fidelity simulation of tokamak plasmas, particular those aimed at characterization of the edge plasma physics, are computationally expensive, so lower cost surrogates are required to enable practical uncertainty estimates. Two surrogate modeling techniques have been explored in the context of tokamak plasma simulations using the XGC family of plasma simulation codes. The first is a response surface surrogate, and the second is an augmented surrogate relying on scenario extrapolation. In addition, to reduce the costs of the XGC simulations, a particle resampling algorithm was developed, which allows marker particle distributions to be adjusted to maintain optimal importance sampling. This means that the total number of particles in and therefore the cost of a simulation can be reduced while maintaining the same accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗