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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Long-Duration Wear Testing of the ASTRAEUS Hall Thruster, Phases II & III: 72 kg Xe Throughput & 1.2 MN-s Total Impulse

The thruster component of the Ascendant Sub-kW Transcelestial Electric Propulsion System successfully completed a propellant throughput demonstration of 72 kg Xe in a long-duration wear test planned to eventually exceed 100 kg. The thruster has been operated for a total of 5450 h over three operating conditions: 3287 h at 300 V – 1000 W, 1745 h at 500 V – 1000 W, and 417 h at 400 V – 600 W. A total of 1.2 MN-s of impulse has been demonstrated to date. Thrust, telemetry, and temperature measurements periodically recorded at four operating conditions spanning 200 – 500 V and 200 – 1000 W revealed constant performance across the full test duration. The erosion of the downstream faces of the thruster was observed to be higher than anticipated after the first 5032 h of thruster operation. This issue was traced to the application of sub-nominal magnetic field strengths throughout the LDWT, caused by an electrical-related ground-support equipment issue, which resulted in a downstream shift of the discharge plasma and enhanced erosion of the thruster’s pole covers. Nevertheless, the thruster performance and operational behaviour remained unaffected. High-fidelity plasma simulations of the erosion of the thruster’s downstream faces, supported by prior short-duration wear testing, suggest a theoretical lifetime capability of >30 kh, corresponding to propellant throughput of >450 kg Xe and a total impulse of >7 MNs, at the 300 V – 1000 W operating condition. The test was voluntarily paused after 68 kg Xe and 72 kg Xe throughput, respectively, to take these erosion measurements at the 500 V and 400 V operating conditions before the test was resumed at a different operating condition.

Goebel, Dan↗
DOE OSTI2026

Plasma processing of SRF cavities at Jefferson Lab: Experiment results and simulation insight

Plasma processing of superconducting radio frequency (SRF) cavities has been an active research effort at Jefferson Lab (JLab) since 2019, aimed at enhancing cavity performance by removing hydrocarbon contaminants and reducing field emission. In this experiment, processing using argon-oxygen and helium-oxygen gas mixtures to find minimum ignition power at different cavity pressure was investigated. Ongoing simulations are contributing to a better understanding of the plasma surface interactions and the fundamental physics behind the process. These simulations, combined with experimental studies, guide the optimization of key parameters such as gas type, RF power, and pressure to ignite plasma using selected higher-order mode (HOM) frequencies. This paper presents experimental data from argon-oxygen and helium-oxygen gas mixture C75 and C100 cavity plasma ignition studies, as well as simulation results for the C100-type cavity based on the COMSOL model previously applied to the C75 cavity.

Accelerator Physics↗
DOE OSTI2023

High-Fidelity Simulations of Plasma-Assisted Oxidation of Hydrocarbon Fuels Using Nanosecond Pulsed Discharges

Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.

adaptive mesh refinement↗
NASA NTRS1987

Numerical simulations of magnetospheric plasmas - Review and quadrennial report to the IUGG

This article provides a broad survey of U.S. progress during the quadrennium 1983-1986 in the category 'numerical simulations of magnetospheric plasmas'. There has been a substantial increase of activity in this area during this period. Simulations have been instrumental in providing valuable insights into large scale dynamic phenomena, nonlinear effects, and complex kinetic phenomena in a wide variety of subject areas, including shocks and double layers, ionosphere-magnetosphere coupling phenomena, and important microphysical processes such as broadband electrostatic noise. The methodology of computer simulation has also been advanced during this quadrennium. Vlasov algorithms have been improved; hybrid codes in 2 and 3D have been developed and applied to magnetospheric problems; and complex problems have been subjected with increasing frequency to a multipronged attack in which several types of simulation models, each designed to accurately model phenomena within a particular range of temporal or spatial scales, are employed synergistically.

Palmadesso, Peter J.↗
DOE OSTI2023

Mesh refinement for anisotropic diffusion in magnetized plasmas

Highly accurate simulation of plasma transport is needed to drive the successful design and operation of magnetically confined fusion reactors. Unfortunately, the extreme anisotropy present in magnetized plasmas results in thin boundary layers that are expensive to resolve. Here, this work investigates how various mesh refinement strategies might reduce that expense to allow for more efficient simulation by comparing standard variable refinement approaches that use a field quantity to an adaptive approach that uses an error estimator. It is first verified that higher order discretization only realizes the proper rate of convergence once the mesh resolves the thin boundary layer, therefore motivating the focusing of refinement on the boundary layer. For three two-dimensional test cases that contain characteristic features of tokamak magnetic fields, an exponential refinement strategy based on the magnetic flux function, which is the standard refinement approach in the field, is compared to an adaptive strategy utilizing the established Zienwiekicz and Zhu error estimator. The adaptive mesh refinement strategy consistently achieves the same accuracy using orders of magnitude less degrees of freedom than either exponential or uniform refinement. This result makes the adaptive refinement strategy more efficient than the exponential refinement strategy while also being more generalizable to problems with complex magnetic geometries. Scaling laws are derived that quantify the improvement in cost of the adaptive refinement strategy over other refinement approaches in 2D and 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Full-torus impurity transport simulation for optimizing plasma discharge operation using a multi-species impurity powder dropper in the large helical device

The transport of impurities supplied by a multi-species impurity powder dropper (IPD) in the large helical device (LHD) is investigated using a three-dimensional peripheral plasma fluid code (EMC3-EIRENE) coupled with a dust transport simulation code (DUSTT). The trajectories of impurity powder particles (Boron, Carbon, Iron, and Tungsten) dropped from the IPD and the impurity transport in the peripheral plasma are studied in a full-torus geometry. The simulation reveals an appropriate size of the impurity powder particles and an optimum operational range of the dust drop rates for investigating the impurity transport without inducing radiation collapse. The simulation also predicts a favourable plasma discharge condition for wall conditioning (boronization) using the IPD in order to deposit boron to high plasma flux and neutral particle density areas in the divertor region in the inboard side of the torus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Numerical simulations of auroral plasma processes. I - Ion beams and conics. II - Electric fields

The characteristics of ion beams and conics as seen in two-dimensional numerical simulations in which the plasmas are driven by current sheets of a finite thickness are described. It is shown that the most energetic ions in the simulations have pitch angles near 90 deg, implying a large perpendicular acceleration of the ions. Large perpendicular electric fields, similar to those measured in electrostatic shocks, are mainly confined near the edges of the current sheets, where E-perpendicular is much larger than E-parallel. Double layers with E-perpendicular about equal to, and spatially separated from, E-parallel form in the interior of wide current sheets. Such double layers have upward electric fields in the region of the upward current. Downward electric fields develop outside the current sheet.

Singh, Nagendra↗
DOE OSTI2026

MAPS: the MFEM Anisotropic Plasma Solver

Simulating magnetically confined fusion plasmas presents a uniquely challenging problem due to the nonlinear anisotropic heat conduction. We introduce the MAPS (MFEM Anisotropic Plasma Solver) tool, which uses a high-order finite element method to compute transport solutions on unstructured meshes. We show results for a set of three 2-D verification tests, two of which demonstrate the expected convergence properties for various mesh resolutions and polynomial degrees. We then discuss the convergence rate for the third test.

Barnett, Rhea [ORNL] (ORCID:0000000317527979)↗
DOE OSTI2021

Overcoming timestep limitations in boosted-frame particle-in-cell simulations of plasma-based acceleration

Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

High voltage solar array experiments

The interaction between the components of a high voltage solar array and a simulated space plasma is studied to obtain data for the design of a high voltage solar array capable of 15kW at 2 to 16kV. Testing was conducted in a vacuum chamber 1.5-m long by 1.5-m diameter having a plasma source which simulated the plasma conditions existing in earth orbit between 400 nautical miles and synchronous altitude. Test samples included solar array segments pinholes in insulation covering high voltage electrodes, and plain dielectric samples. Quantitative data are presented in the areas of plasma power losses, plasma and high voltage induced damage, and dielectric properties. Limitations of the investigation are described.

Kennerud, K. L.↗
NASA NTRS2006

Towards the Identification of the Keeper Erosion Cause(s): Numerical Simulations of the Plasma and Neutral Gas Using the Global Cathode Model OrCa2D-II

Numerical simulations with the time-dependent Orificed Cathode (OrCa2D-II) computer code show that classical enhancements of the plasma resistivity can not account for the elevated electron temperatures and steep plasma potential gradients measured in the plume of a 25-27.5 A discharge hollow cathode. The cathode, which employs a 0.11-in diameter orifice, was operated at 5.5 sccm without an applied magnetic field using two different anode geometries. It is found that anomalous resistivity based on electron-driven instabilities improves the comparison between theory and experiment. It is also estimated that other effects such as the Hall-effect from the self-induced magnetic field, not presently included in OrCa2D-II, may contribute to the constriction of the current density streamlines thus explaining the higher plasma densities observed along the centerline.

plasma modeling↗
DOE OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

FY2018 Fusion Theory and Simulation Milestone on Plasma Materials Interaction (PMI) (Final Report)

The interaction of the boundary plasma with the material surfaces in magnetically confined plasmas is among the most critical problems in fusion energy science. In FY2018, perform high performance computing simulations with coupled boundary plasma physics and materials surface models to predict the fuel recycling and tritium retention of the ITER divertor for D-T burning plasma conditions, accounting for erosion, re-deposition and impurity transport in the plasma boundary, and an initial evaluation of the influence of material deposition on the recycling and retention.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Time-Dependent Boundary Modeling to Inform Design of SPARC Diagnostic and Actuators

The purpose of this work was to inform the design of diagnostic and actuator systems for a new experimental facility for Commonwealth Fusion Systems (CFS), called SPARC [Creely2020], using time-dependent plasma boundary simulations with the Scrape Off Layer Plasma Simulator (SOLPS) code [Wiesen2015]. The SPARC tokamak is a device that aims to demonstrate fusion energy production and high-field plasma operating scenarios. The developed technologies will be incorporated into a net-electricity pilot that will demonstrate approximately 200MW of electric power. A major difference from other private industry fusion companies is the use of a standard aspect ratio, but at a high magnetic field strength, which provides a well-established physics basis and a modest extrapolation to reactor operation. Time-dependent simulations were used to inform the design of main-chamber plasma facing components, which are protected by feedback control of the heat and particle loads during high power operation. The simulations provide the timescales and magnitudes of plasma and neutral particle response useful for determining the position and operation of diagnostics and actuators that will enable control. The plasma response and phase-space diagrams also provide input into advanced model-based control schemes, which are envisioned for later SPARC operation. Model-based control can reduce the risk of component failure, which would require expensive in-situ repairs and the associated delays, by predicting the system response to actuators over a short time-horizon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Performance of coated columbium and tantalum alloys in plasma arc reentry simulation tests

The evaluation of coated refractory metals screened in stagnation model plasma arc tests is reported. Columbium alloys FS-85, C-129Y, and Cb-752 coated with Si-20Cr-20Fe (R512E) were tested at 1390 C. Three silicide coatings on Ta-10W were tested at 1470 C. Half-hour cycles and a 6500 N/sqm stagnation pressure were used. The best R512E coated columbium alloy was FS-85 with first local coating breakdowns occurring in 12 to 50 cycles. At coating defects, low metal recession rates (0.005 mm/min) were generally observed on coated columbium alloys while high rates (0.15 mm/min) were observed on coated Ta-10W. Coated columbium suffered large emittance losses (to below 0.7) due to surface refractory metal pentoxide formation.

Levine, S. R.↗
DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

General Relativistic Simulations of Magnetized Plasmas Around Merging Supermassive Black Holes

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this paper we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe, total amplification of the magnetic field of approx 2 orders of magnitude which is driven by the accretion onto the binary and that leads to stronger electromagnetic signals than in the force-free regime where such amplifications are not possible.

Giacomazzo, Bruno↗
NASA NTRS1990

Simulation studies of plasma waves in the electron foreshock - The generation of downshifted oscillations

The generation of waves with frequencies downshifted from the plasma frequency, as observed in the electron foreshock, is analyzed by particle simulation. Wave excitation differs fundamentally from the familiar excitation of the plasma eigenmodes by a gentle bump-on-tail electron distribution. Beam modes are destabilized by resonant interaction with bulk electrons, provided the beam velocity spread is very small. These modes are stabilized, starting with the higher frequencies, as the beam is broadened and slowed down by the interaction with the wave spectrum. Initially a very cold beam is also capable of exciting frequencies considerably above the plasma frequency, but such oscillations are quickly stabilized. Low-frequency modes persist for a long time, until the bump in the electron distribution is completely 'ironed' out. This diffusion process also is quite different from the familiar case of well-separated beam and bulk electrons. A quantitative analysis of these processes is carried out.

Dum, C. T.↗