Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma simulations”

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 361 records · Page 20

Simulating Boundary Region Plasma Instabilities

The boundary regions of low-temperature plasmas are known to be susceptible to kinetic instabilities, which can affect the energies and fluxes of particles directed at the material boundary. For example, both the ion acoustic instability as well as an instability near the electron plasma frequency have been observed. Particle-in-cell (PIC) simulation is a tool that, alongside experiments, can capture the effects these instabilities have on the particle distribution functions. Ultimately, simulations can determine under what conditions these effects are significant by comparing to theoretical predictions and explore conditions unamenable to experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of energetic particle transport models for time-dependent integrated simulations of burning plasmas (Final ST Report)

The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental simulation of space plasma interactions with high voltage solar arrays

Operating high voltage solar arrays in the space environment can result in anomalously large currents being collected through small insulation defects. Tests of simulated defects have been conducted in a 45-cm vacuum chamber with plasma densities of 100,000 to 1,000,000/cu cm. Plasmas were generated using an argon hollow cathode. The solar array elements were simulated by placing a thin sheet of polyimide (Kapton) insulation with a small hole in it over a conductor. Parameters tested were: hole size, adhesive, surface roughening, sample temperature, insulator thickness, insulator area. These results are discussed along with some preliminary empirical correlations.

Stillwell, R. P.↗

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beam-plasma interactions in space experiments - A simulation study

The plasma environment in the vicinity of a spacecraft during the injection of dense electron beams is studied using a two-dimensional, isolated-system electrostatic simulation model. The dependence of the beam stagnation time on the beam width and energy is examined. It is found that the relative size of the beam stagnation time and the ambient-plasma response time determines the environment of the spacecraft. The case of cross-field injection with beam stagnation time greater than plasma response time is discussed in detail. Also, the nature of the beam properties, plasma response, and wave spectra are considered.

Pritchett, P. L.↗

Monte Carlo Simulation of Nonlinear Radiation Induced Plasmas

A Monte Carlo simulation model for radiation induced plasmas with nonlinear properties due to recombination was, employing a piecewise linearized predict-correct iterative technique. Several important variance reduction techniques were developed and incorporated into the model, including an antithetic variates technique. This approach is especially efficient for plasma systems with inhomogeneous media, multidimensions, and irregular boundaries. The Monte Carlo code developed has been applied to the determination of the electron energy distribution function and related parameters for a noble gas plasma created by alpha-particle irradiation. The characteristics of the radiation induced plasma involved are given.

Wang, B. S.↗

Observations of heteroclinic bifurcations in resistive magnetohydrodynamic simulations of the plasma response to resonant magnetic perturbations

A class of topological magnetic island bifurcations that has not previously been observed in toroidal plasmas is described. Increasing an externally applied three-dimensional magnetic field in resistive magnetohydrodynamic simulations results in the asymmetric elongation of resonant island flux surfaces followed by a sequence of heteroclinic bifurcations. Here, these bifurcations produce new sets of hyperbolic-elliptic fixed points as predicted by the Poincaré-Birkoff fixed point theorem. Field line calculations verify that the new fixed points do not connect to those of the prebifurcated islands as required for heteroclinic bifurcations on a torus with winding numbers composed of common integer factors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation studies of plasma waves in the electron foreshock - The generation of Langmuir waves by a gentle bump-on-tail electron distribution

Particle simulation experiments were used to study the basic physical ingredients needed for building a global model of foreshock wave phenomena. In particular, the generation of Langmuir waves by a gentle bump-on-tail electron distribution is analyzed. It is shown that, with appropriately designed simulations experiments, quasi-linear theory can be quantitatively verified for parameters corresponding to the electron foreshock.

Dum, C. T.↗

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)↗

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.↗

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.↗

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↗

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.↗

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.↗

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.↗

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.↗

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↗