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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Numerical studies of electron dynamics in oblique quasi-perpendicular collisionless shock waves

Linear and nonlinear electron damping of the whistler precursor wave train to low Mach number quasi-perpendicular oblique shocks is studied using a one-dimensional electromagnetic plasma simulation code with particle electrons and ions. In some parameter regimes, electrons are observed to trap along the magnetic field lines in the potential of the whistler precursor wave train. This trapping can lead to significant electron heating in front of the shock for low beta(e). Use of a 64-processor hypercube concurrent computer has enabled long runs using realistic mass ratios in the full particle in-cell code and thus simulate shock parameter regimes and phenomena not previously studied numerically.

Liewer, P. C.↗
NASA NTRS1992

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗
NASA NTRS1992

Perspectives of future observations on particle acceleration in the heliosphere

The development of the quantitative study of particle acceleration processes in the outer heliosphere into the foreseeable future is outlined. In the near future, the combination of spacecraft, new improved instrumentation, and powerful plasma simulations are expected to enhance the understanding of acceleration processes in the near-earth environment, on the sun, and in the outer heliosphere. The multispacecraft mission CLUSTER will allow the most detailed view of the bow shock, while spacecraft passing through various regions of the inner and outer heliosphere are to provide the large-scale view needed to understand the heliospheric boundary. Improvements in ionic charge resolution and collecting power are to elucidate fractionation processes.

Moebius, Eberhard↗
NASA NTRS2005

Particle acceleration magnetic field generation, and emission in Relativistic pair jets

Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) are responsible for particle acceleration in relativistic pair jets. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic pair jet propagating through a pair plasma. Simulations show that the Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. Simulation results show that this instability generates and amplifies highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The "jitter' I radiation from deflected electrons can have different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. The growth rate of the Weibel instability and the resulting particle acceleration depend on the magnetic field strength and orientation, and on the initial particle distribution function. In this presentation we explore some of the dependencies of the Weibel instability and resulting particle acceleration on the magnetic field strength and orientation, and the particle distribution function.

Nishikawa, K.-I.↗
NASA NTRS2005

Ion thruster life models

In this paper we review the most recent 2-D and 3-D plasma simulation models for ion thruster hollow cathodes, discharge chambers, and ion optics.

grid life↗
NASA NTRS2009

Space Flight Plasma Data Analysis

This slide presentation reviews a method to analyze the plasma data that is reported on board the International Space station (ISS). The Floating Potential Measurement Unit (FPMU), the role of which is to obtain floating potential and ionosphere plasma measurements for validation of the ISS charging model, assess photo voltaic array variability and interpreting IRI predictions, is composed of four probes: Floating Potential Probe (FPP), Wide-sweep Langmuir Probe (WLP), Narrow-sweep Langmuir Probe (NLP) and the Plasma Impedance Probe (PIP). This gives redundant measurements of each parameter. There are also many 'boxes' that the data must pass through before being captured by the ground station, which leads to telemetry noise. Methods of analysis for the various signals from the different sets are reviewed. There is also a brief discussion of LP analysis of Low Earth Orbit plasma simulation source.

Wright, Kenneth H.↗
NASA NTRS2009

Splanchnic Compression Improves the Efficacy of Compression Stockings to Prevent Orthostatic Intolerance

Purpose: Post-spaceflight orthostatic intolerance (OI) is observed in 20-30% of astronauts. Previous data from our laboratory suggests that this is largely a result of decreased venous return. Currently, NASA astronauts wear an anti-gravity suit (AGS) which consists of inflatable air bladders over the calves, thighs and abdomen, typically pressurized from 26 to 78 mmHg. We recently determined that, thigh-high graded compression stockings (JOBST , 55 mmHg at ankle, 6 mmHg at top of thigh) were effective, though to a lesser degree than the AGS. The purpose of this study was to evaluate the addition of splanchnic compression to prevent orthostatic intolerance. Methods: Ten healthy volunteers (6M, 4F) participated in three 80 head-up tilts on separate days while (1) normovolemic (2) hypovolemic w/ breast-high compression stockings (BS)(JOBST(R), 55 mmHg at the ankle, 6 mmHg at top of thigh, 12 mmHg over abdomen) (3) hypovolemic w/o stockings. Hypovolemia was induced by IV infusion of furosemide (0.5 mg/kg) and 48 hrs of a low salt diet to simulate plasma volume loss following space flight. Hypovolemic testing occurred 24 and 48 hrs after furosemide. One-way repeated measures ANOVA, with Bonferroni corrections, was used to test for differences in blood pressure and heart rate responses to head-up tilt, stand times were compared using a Kaplan-Meyer survival analysis. Results: BS were effective in preventing OI and presyncope in hypovolemic test subjects ( p = 0.015). BS prevented the decrease in systolic blood pressure seen during tilt in normovolemia (p < 0.001) and hypovolemia w/o countermeasure (p = 0.005). BS also prevented the decrease in diastolic blood pressure seen during tilt in normovolemia (p = 0.006) and hypovolemia w/o countermeasure (p = 0.041). Hypovolemia w/o countermeasure showed a higher tilt-induced heart rate increase (p = 0.022) than seen in normovolemia; heart rate while wearing BS was not different than normovolemia (p = 0.353). Conclusion: BS may be an effective countermeasure to post-space flight OI. The addition of splanchnic compression is more effective than the previous thigh-high garments. These stockings are readily available, inexpensive, and can be worn for days following landing as astronauts re-adapt to Earth gravity.

Platts, Steven H.↗
NASA NTRS2010

Space Weather Models and Their Validation and Verification at the CCMC

The Community Coordinated l\lodeling Center (CCMC) is a US multi-agency activity with a dual mission. With equal emphasis, CCMC strives to provide science support to the international space research community through the execution of advanced space plasma simulations, and it endeavors to support the space weather needs of the CS and partners. Space weather support involves a broad spectrum, from designing robust forecasting systems and transitioning them to forecasters, to providing space weather updates and forecasts to NASA's robotic mission operators. All of these activities have to rely on validation and verification of models and their products, so users and forecasters have the means to assign confidence levels to the space weather information. In this presentation, we provide an overview of space weather models resident at CCMC, as well as of validation and verification activities undertaken at CCMC or through the use of CCMC services.

Hesse, Michael↗
NASA NTRS2012

Comparison of Computed and Measured Performance of a Pulsed Inductive Thruster Operating on Argon Propellant

Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic field. A recent review of the developmental history of planar-geometry pulsed inductive thrusters, where the coil take the shape of a flat spiral, can be found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT)[2, 3] and the Faraday Accelerator with Radio-frequency Assisted Discharge (FARAD)[4]. There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et al.[5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 1 provides a description of the electrical coupling between the current flowing in the thruster I1 and the plasma current I2. Recently, the model was upgraded to include an equation governing the deposition of energy into various modes present in a pulsed inductive thruster system (acceleration, magnetic flux generation, resistive heating, etc.)[7]. An MHD description of the plasma energy density evolution was tailored to the thruster geometry by assuming only one-dimensional motion and averaging the plasma properties over the spatial dimensions of the current sheet to obtain an equation for the time-evolution of the total energy. The equation set governing the dynamics of the coupled electrodynamic-current sheet system is composed of first-order, coupled ordinary differential equations that can be easily solved numerically without having to resort to much more complex 2-D finite element plasma simulations.

Polzin, Kurt A.↗
NASA NTRS2021

Long-Duration Wear Testing of the ASTRAEUS Thruster, Phase I: 50 kg Xe Throughput

The thruster component of the Ascendant Sub-kW Transcelestial Electric Propulsion System successfully completed a propellant throughput demonstration of 50 kg Xe. The thruster was operated for a total of 3,287 h, of which 97.9% was spent at 300 V – 1000 W, and generated 790 kN-s of total impulse. 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. High-fidelity plasma simulations of the erosion of the thruster’s downstream faces 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 50 kg Xe throughput to take these erosion measurements before the test was resumed at a different operating condition.

Goebel, Dan↗
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 OSTI2014

A Screen Space GPGPU Surface LIC Algorithm for Distributed Memory Data Parallel Sort Last Rendering Infrastructures

The surface line integral convolution(LIC) visualization technique produces dense visualization of vector fields on arbitrary surfaces. We present a screen space surface LIC algorithm for use in distributed memory data parallel sort last rendering infrastructures. The motivations for our work are to support analysis of datasets that are too large to fit in the main memory of a single computer and compatibility with prevalent parallel scientific visualization tools such as ParaView and VisIt. By working in screen space using OpenGL we can leverage the computational power of GPUs when they are available and run without them when they are not. We address efficiency and performance issues that arise from the transformation of data from physical to screen space by selecting an alternate screen space domain decomposition. We analyze the algorithm's scaling behavior with and without GPUs on two high performance computing systems using data from turbulent plasma simulations.

Loring, Burlen↗
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 OSTI2016

Dust on a surface in a plasma: A charge simulation

An electrically isolated particle on a substrate surface will be electrically charged when a plasma is applied above it. We show that the magnitude of this charge is determined by a balance of the impeding ion and electron fluxes that are strongly influenced by the nearby substrate. By simulating this process, we find that the charge density of the particle can be much higher than that of the substrate. This is due to the height of the particle, which causes additional electron collection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗