SEARCH · Engineering Papers
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.
Together with Realta Fusion we developed advanced simulations using WarpX for magnetic confinement fusion devices and advanced neutral beam systems for plasma heating.
Any process which may be modeled statistically for a simple plasma particle may be modeled by the particle in cell technique. The success of the calculation is then dependent on having a large enough number of particles that the statistical treatment is meaningful. Thus it is possible to include the effects of secondary emission, backscattering, charge sticking and possibly dielectric breakdown, photoemission, and spallation. The first plasma dielectric interaction included in the computer code for simulating plasma insulator interactions is secondary electron emission. A calculated current density vs. voltage curve is presented and compared to an experimental curve.
Together with Realta Fusion, we set up several simulation scenarios and ran simulations using the WarpX code which is developed by LBNL and especially by members of the LBNL team. We investigated beam transport and beam neutralization and ran optimization and benchmark simulations. To achieve this, several new features were added to the warpX code. These code additions will also be beneficial to other users of warpX including several fusion companies and researchers at National Laboratories. Furthermore, we supported Realta Fusion in setting up their own instances of WarpX. We presented our results at the American Physical Society Division of Plasma Physics meeting and are in the process of writing a publication that will be submitted to a peer-reviewed journal to share our results with the broader community.
Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.
In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.
Within a permanently shadowed lunar crater the horizontal flow of solar wind is obstructed by upstream topography, forming a plasma wake that electrostatically diverts ions toward the crater floor and generates a surface potential that can reach kilovolts. In the present work kinetic plasma simulations are employed to investigate the morphology of a lunar crater wake during passage of a solar storm. Results are cast in terms of leading dimensionless ratios including the ion Mach number, ratio of crater depth to plasma Debye length, peak secondary electron yield, and electron temperature vs. electron impact energy at peak secondary yield. This small set of ratios allows generalization to a much wider range of scenarios. The kinetic simulation results are fed forward into an equivalent-circuit model of a roving astronaut. In very low-plasma-current environments triboelectric charging of the astronaut suit becomes effectively perpetual, representing a critical engineering concern for roving within shadowed lunar regions. Finally, simulated ion fluxes are used to explore sputtering and implantation processes within an idealized crater. It is suggested that the physics of plasma mini-wakes formed in the vicinity of permanently shadowed topography may play a critical role in modulating the enigmatic spatial distribution of volatiles at the lunar poles.
The particle simulations in a Variable Specific Impulse Magnetoplasma Rocket (VASIMR) currently include self-consistent calculation of. 1) stationary magnetic field in plasma, 2) ion density and velocity, 3) ion-cyclotron radio-frequency heating, 4) ambipolar electric field. The assumptions of quasineutral and collissionless plasma are based on the range of operating VASIMR parameters. The main motivation for the particle simulation in VASIMR is plasma detachment from the magnetic field in the exhaust area. The plasma detachment is caused mainly by the Larmor radius increase. The plasma beta effect on detachment is observed and investigated as well. The results of particle simulations are compared with those from MHD simulations.
The Rice convection model (RCM) has been modified to a form suitable for Jupiter (RCM-J) to study plasma interchange motion in and near the Io plasma torus. The net result of the interchange is that flux tubes, heavily loaded with torus plasma, are transported outward, to be replaced by tubes containing little low-energy (less than 1 keV) plasma. The process is numerically simulated in terms of time evolution from an initial torus that is longitudinally asymmetric and with gradually decreasing density outward from Io's orbit. In the simulations, the nonlinear stage of the instability characteristically exhibits outreaching fingers of heavily-loaded flux tubes that lengthen at an accelerating rate. The principal finding is that the primary geometrical form of outward transport of torus plasma in Jupiter's magnetosphere is through long, outward-moving fingers of plasma. In the simulations, the fingers mainly form in the active sector of the Io torus (the heavier side of the asymmetric torus), and they are spaced longitudinally roughly 20 deg apart.
The Radio Plasma Imager (RPI) will be the first-of-its kind instrument designed to use radio wave sounding techniques to perform repetitive remote sensing measurements of electron number density (N(sub e)) structures and the dynamics of the magnetosphere and plasmasphere. RPI will fly on the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission to be launched early in the year 2000. The design of the RPI is based on recent advances in radio transmitter and receiver design and modern digital processing techniques perfected for ground-based ionospheric sounding over the last two decades. Free-space electromagnetic waves transmitted by the RPI located in the low density magnetospheric cavity will be reflected at distant plasma cutoffs. The location and characteristics of the plasma at those remote reflection points can then be derived from measurements of the echo amplitude, phase, delay time, frequency, polarization, Doppler shift, and echo direction. The 500 m tip-to-tip X and Y (spin plane) antennas and 20 m boom Z axis antenna on RPI will be used to measures echoes coming from distances of several R(sub E).
Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.
A plasma simulation code is presented that models the time-dependent plasma properties in the vicinity of a spherical, charged spacecraft. After showing agreement with analytic, steady-state theories and ATS-6 satellite data, the following three problems are treated: (1) transient pulses from photoemission at various emission temperatures and ambient plasma conditions, (2) spacecharge limited emission, and (3) simulated plasma oscillations in the long wavelength limit.
The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.
Gas-core nuclear rocket engine simulation by induction coupled plasma torch, determining gas flow patterns
Electrostatic plasma double layers are numerically simulated by means of a magnetized 2 1/2 dimensional particle in cell method. The investigation of planar double layers indicates that these one dimensional potential structures are susceptible to periodic disruption by instabilities in the low potential plasmas. Only a slight increase in the double layer thickness with an increase in its obliqueness to the magnetic field is observed. Weak magnetization results in the double layer electric field alignment of accelerated particles and strong magnetization results in their magnetic field alignment. The numerical simulations of spatially periodic two dimensional double layers also exhibit cyclical instability. A morphological invariance in two dimensional double layers with respect to the degree of magnetization implies that the potential structures scale with Debye lengths rather than with gyroradii. Electron beam excited electrostatic electron cyclotron waves and (ion beam driven) solitary waves are present in the plasmas adjacent to the double layers.
We present our experience using containers to scale up a massive ensemble of coupled I/O bound workloads on the NERSC Cori supercomputer. We describe the design of a hierarchical simulation structure using the Integrated Plasma Simulator (IPS) that enables the flexible execution of coupled simulations at the system, node, and core level using the same coupling abstraction and API. The hierarchical design allows for the node-level execution to be efficiently executed using containers while not impacting the structure of the simulation at the system level. We demonstrate the viability of the approach by presenting experimental results from applications in coupled fusion plasma simulations that illustrate the performance impact of using containers to deploy the node-level workloads, in conjunction with the user mountable XFS file systems to ameliorate the load on the Lustre parallel file system. We also present results from production runs showing the ability of the ensemble simulations to scale to hundreds of Cori Haswell nodes, with little or no overhead.
The paper describes the recent development in the simulation of the ion-cyclotron acceleration of the plasma in the VASIMR experiment. The modeling is done using an improved EMIR code for RF field calculation together with particle trajectory code for plasma transport calculat ion. The simulation results correlate with experimental data on the p lasma loading and predict higher ICRH performance for a higher density plasma target. These simulations assist in optimizing the ICRF anten na so as to achieve higher VASIMR efficiency.
A mathematical model was developed to simulate flow control applications using plasma actuators. The effects of the plasma actuators on the external flow are incorporated into Navier Stokes computations as a body force vector. In order to compute this body force vector, the model solves two additional equations: one for the electric field due to the applied AC voltage at the electrodes and the other for the charge density representing the ionized air. The model is calibrated against an experiment having plasma-driven flow in a quiescent environment and is then applied to simulate a low pressure turbine flow with large flow separation. The effects of the plasma actuator on control of flow separation are demonstrated numerically.
Explore the source record for details and available documents.