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

Ion motion in the wake driven by long particle bunches in plasmas

In this paper, we explore the role of the background plasma ion motion in self-modulated plasma wakefield accelerators. We employ Dawson's plasma sheet model to derive expressions for the transverse plasma electric field and ponderomotive force in the narrow bunch limit. We use these results to determine the on-set of the ion dynamics and demonstrate that the ion motion could occur in self-modulated plasma wakefield accelerators. Simulations show the motion of the plasma ions can lead to the early suppression of the self-modulation instability and of the accelerating fields. The background plasma ion motion can nevertheless be fully mitigated by using plasmas with heavier plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas

Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Application of the parareal algorithm to simulations of ELMs in ITER plasma

This paper explores the application of the parareal algorithm to simulations of ELMs in ITER plasma. The primary focus of this research is identifying the parameters that lead to optimum performance. Since the plasma dynamics vary extremely fast during an ELM cycle, a straightforward application of the algorithm is not possible and a modification to the standard parareal correction is implemented. The size of the time chunks also have an impact on the performance and needs to be optimized. in conclusion, a computational gain of 7.8 is obtained with 48 processors to illustrate that the parareal algorithm can be successfully applied to ELM plasma.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Ground-Based Simulation of Low-Earth Orbit Plasma Conditions: Plasma Generation and Characterization

A 16-cm diameter plasma source operated on argon is described that is capable of producing a plasma environment that closely simulates the low Earth orbit (LEO) conditions experienced by satellites in the altitude range between 300 to 500 km. The plasma source uses a transverse-field magnetic filter, and has been successful in producing low electron temperature plasmas that contain streaming ion populations. Both of these characteristics are important because the plasma in LEO is relatively cold (e.g., Te approx. 0.1 eV) and the ram energy of the ions due to the motion of the satellite relative to the LEO plasma is high (e.g., 7,800 m/s which corresponds to approx. 5 eV for O+ ions). Plasma source operational conditions of flow rate and discharge power are presented that allow the electron temperature to be adjusted over a range from 0.14 to 0.4 eV. The expanding plasma flow field downstream of the source contains both low-energy, charge-exchange ions and streaming ions with energies that are adjustable over a range from 4 eV to 6 eV. At low flow rates and low facility pressures, the streaming ion component of the ion population comprises over 90% of the total plasma density. In the work described herein, a large area retarding potential analyzer was used to measure both electron and ion energy distribution functions in the low density, expanding plasma produced downstream of the plasma source. The benefits of using this type of plasma diagnostic tool in easily perturbed, low-density plasma are identified, and techniques are also discussed that can be used to perform real-time measurements of electron temperature. Finally, recommendations are made that may enable lower electron temperatures to be produced while simultaneously decreasing the plasma source flow rate below 1 to 2 sccm.

Williams, John D.↗
DOE OSTI2013

Simulation study of wave phenomena from the sheath region in single frequency capacitively coupled plasma discharges; field reversals and ion reflection

Capacitively coupled radio-frequency (RF) discharges have great significance for industrial applications. Collisionless electron heating in such discharges is important, and sometimes is the dominant mechanism. This heating is usually understood to originate in a stochastic interaction between electrons and the electric fields. However, other mechanisms may also be important. There is evidence of wave emission with a frequency near the electron plasma frequency, i.e., ω{sub pe}, from the sheath region in collisionless capacitive RF discharges. This is the result of a progressive breakdown of quasi-neutrality close to the electron sheath edge. These waves are damped in a few centimeters during their propagation from the sheath towards the bulk plasma. The damping occurs because of the Landau damping or some related mechanism. This research work reports that the emission of waves is associated with a field reversal during the expanding phase of the sheath. Trapping of electrons near to this field reversal region is observed. The amplitude of the wave increases with increasing RF current density amplitude J(tilde sign){sub 0} until some maximum is reached, beyond which the wave diminishes and a new regime appears. In this new regime, the density of the bulk plasma suddenly increases because of ion reflection, which occurs due to the presence of strong field reversal near sheath region. Our calculation shows that these waves are electron plasma waves. These phenomena occur under extreme conditions (i.e., higher J(tilde sign){sub 0} than in typical experiments) for sinusoidal current waveforms, but similar effects may occur with non-sinusoidal pulsed waveforms for conditions of experimental interest, because the rate of change of current is a relevant parameter. The effect of electron elastic collisions on plasma waves is also investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulations of Hypervelocity Impact Plasmas

The space environment is fraught with complex plasmas spanning a wide range of densities and temperatures. Much of space plasma research has focused on large-scale changes in the ambient plasma and fields from the Sun to the Earth’s atmosphere, driven primarily by the solar wind. Yet one particular type of space plasma, which forms from a hypervelocity impact (HVI), remains poorly understood. Hypervelocity impactors include both meteoroids and space debris. Meteoroids are naturally occurring objects in space that travel between 11 and 72.8 km/s and originate primarily from comets and asteroids. In contrast, space debris are human-made objects with speeds typically < 11 km/s. Hypervelocity impactors routinely hit spacecraft, yet the physics behind the formation of the plasma and the dynamics of its expansion remain largely unknown. The complexity of this phenomenon necessitates a research approach that includes both experimental studies and numerical simulation in order to understand the underlying physical processes that occur upon formation and expansion of the impact plasma. Our research has focused on providing a comprehensive understanding of plasma generated by hypervelocity impacts by meteoroids and space debris on spacecraft in order to characterize the behavior of the expanding plasma and its interactions with the ambient environment. Previously. we conducted experimental campaigns at a dust accelerator facility that can accelerate particles up to 60 km/s, which is representative of meteoroid speeds, and at a light gas gun facility that can accelerate larger projectiles up to 7 km/s, which is representative of orbital debris. The experiments included plasma, optical and radio frequency (RF) sensors in order to understand the dynamics and associated RF emission resulting from hypervelocity plasmas. For this research, we developed and applied machine learning algorithms to identify which type of impactor would produce RF and physics-based models to determine the source of the RF. This was a one-year research program that resulted in 2 refereed journal publications (uploaded separately).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Experiment and simulation on one-dimensional plasma photonic crystals

The transmission characteristics of microwaves passing through one-dimensional plasma photonic crystals (PPCs) have been investigated by experiment and simulation. The PPCs were formed by a series of discharge tubes filled with argon at 5 Torr that the plasma density in tubes can be varied by adjusting the discharge current. The transmittance of X-band microwaves through the crystal structure was measured under different discharge currents and geometrical parameters. The finite-different time-domain method was employed to analyze the detailed properties of the microwaves propagation. The results show that there exist bandgaps when the plasma is turned on. The properties of bandgaps depend on the plasma density and the geometrical parameters of the PPCs structure. The PPCs can perform as dynamical band-stop filter to control the transmission of microwaves within a wide frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗
DOE OSTI2020

Test for Reconditioning RA Waste with Simulated Bitumen and Concrete in a 1,2 MW Plasma Test Facility - 20092

The operation and maintenance of nuclear power plants, the non nuclear fuel cycle, etc generate low-level radioactive waste which, along with the historical radioactive waste from past nuclear activities, needs to be treated and stored, awaiting final disposal. Plasma technology offers a very effective way of treating this waste with a high volume reduction factor (VRF), free from organics, liquids and moisture, and meets without a doubt the acceptance criteria for safe storage and disposal. By means of a plasma beam of approximately 5000 deg. C, the inorganic materials are melted into a glassy slag, containing most of the radioactive isotopes while the organic material is gasified, oxidized and purified in an off-gas cleaning system. First the paper describes the new full-scale Plasma Melting Facility (PMF) at the Kozloduy Nuclear Power Plant in Bulgaria which was taken in nuclear operation In May 2018. The plant has a capacity of 250 tons per year and the maximum contact dose rates of the incoming waste is 2 mSv/h. Different mixtures of radioactive waste packed in 200 l drums were successfully treated resulting in a glassy slag free from liquids and organic material with an important volume reduction factor (VRF). The Project was co-financed through a grant by Kozloduy International Decommissioning fund (KIDSF) administrated by the EBRD through Bulgarian national funding. Plasma is a suitable technology for treatment of problematic waste or even reconditioning waste so Belgoprocess was contracted to do plasma tests with simulated conditioned waste types. One can do tests on a laboratory scale on smaller samples and torch capacities of e.g. 50 kW but Belgoprocess wanted to do more realistic and reliable tests. So Belgoprocess contracted Phoenix Solutions Co who has a full-scope test facility equipped with a 1200 kW plasma torch for full-scope treatment of simulated conditioned waste. For a first confidential contract simulated 200 l (55 gallon) bitumen drums were treated. The drums contained different pucks of compacted waste such as rags, used filters, granulates, etc. The pucks were stacked in the 200 l drums and subsequently embedded with bitumen. A total of 6 drums were treated in the plasma facility. For a second contract simulated homogeneous 200 l (55 gallon) concrete drums with on the one hand concentrates and on the other hand spent resins were selected. A total of 6 drums with concrete and spent resins were treated and melted in the plasma testing facility. The paper describe the test facility, volume reduction factor (VRF) of different waste streams and most important parameters. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1995

Bell Helicopter Advanced Rotocraft Transmission (ART) program

Future rotorcraft transmissions require key emerging material and component technologies using advanced and innovative design practices in order to meet the requirements for a reduced weight to power ratio, a decreased noise level, and a substantially increased reliability. The specific goals for the future rotorcraft transmission when compared with a current state-of-the-art transmission (SOAT) are: (1) a 25 percent weight reduction; (2) a 10 dB reduction in the transmitted noise level; and (3) a system reliability of 5000 hours mean-time-between-removal (MTBR) for the transmission. This report summarizes the work conducted by Bell Helicopter Textron, Inc. to achieve these goals under the Advanced Rotorcraft Transmission (ART) program from 1988 to 1995. The reference aircraft selected by BHTI for the ART program was the Tactical Tiltrotor which is a 17,000 lb gross weight aircraft. A tradeoff study was conducted comparing the ART with a Selected SOAT. The results showed the ART to be 29 percent lighter and up to 13 dB quieter with a calculated MTBR in excess of 5000 hours. The results of the following high risk component and material tests are also presented: (1) sequential meshing high contact ratio planetary with cantilevered support posts; (2) thin dense chrome plated M50 NiL double row spherical roller planetary bearings; (3) reduced kinematic error and increased bending strength spiral bevel gears; (4) high temperature WE43 magnesium housing evaluation and coupon corrosion tests; (5) flexure fatigue tests of precision forged coupons simulating precision forged gear teeth; and (6) flexure fatigue tests of plasma carburized coupons simulating plasma carburized gear teeth.

Henry, Zachary S.↗
DOE OSTI2019

Hybrid particle-in-cell simulations of laser-driven plasma interpenetration, heating, and entrainment

Kinetic-ion, quasineutral, fluid-electron particle-in-cell simulations of interpenetrating carbon–carbon plasma flows in 2D RZ cylindrical geometry are presented. The simulations are initialized with solid density targets that are subsequently irradiated by 10 14 W/cm 2 intensity lasers using a raytracing package. The ablation, interpenetration, heating, slowing, entrainment, and stagnation of the plasma flows evolve self-consistently within the code. The particle density, velocity phase space, and fits to the velocity distribution functions are used, along with analytical collisional stopping rates, to interpret the dynamics of the flow evolution. Comparisons to multifluid simulations are described and used to highlight ion-kinetic effects in the setup. Synthetic Thomson scattering diagnostic signals are generated using detailed knowledge of the plasma distribution functions. The large scale of the system, 1 × 1 mm for 2 ns, and the detailed dynamics extracted demonstrate that such hybrid codes are powerful tools for the design and evaluation of laboratory-scale high-energy-density plasma physics experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Simulation Study of Landau Damping Near the Persisting to Arrested Transition

A one-dimensional electrostatic filtered Vlasov-Poisson simulation study is discussed. The transition from persisting to arrested Landau damping that is produced by increasing the strength of a sinusoidal perturbation on a background Vlasov-Poisson equilibrium is explored. Emphasis is placed on observed features of the electron phase-space distribution when the perturbation strength is near the transition value. A single ubiquitous waveform is found perturbing the space-averaged phase-space distribution at almost any time in all of the simulations; the sole exception is the saturation stage that can occur at the end of the arrested damping scenario. This waveform contains relatively strong, very narrow structures in velocity bracketing (plus or minus velocity (sub res)) - the velocities at which electrons must move to traverse the dominant field mode wavelength in one of its oscillation periods - and propagating with (plus or minus velocity (sub res)) respectively. Local streams of electrons are found in these structures crossing the resonant velocities from low speed to high speed during Landau damping and from high speed to low speed during Landau growth. At the arrest time, when the field strength is briefly constant, these streams vanish. It is conjectured that the expected transfer of energy between electrons and field during Landau growth or damping has been visualized for the first time. No evidence is found in the phase-space distribution to support recent well-established discoveries of a second-order phase transition in the electric field evolution. While trapping is known to play a role for larger perturbation strengths, it is shown that trapping plays no role at any time in any of the simulations near the transition perturbation strength.

plasma simulation↗
DOE OSTI2015

Numerical Uncertainty Estimation for Stochastic Particle-in-Cell Simulations Applied to Verification and Validation.

Systematic veri cation and validation (V&V) is necessary to establish the credibility for high consequence simulations. In this paper, we focus on a radiation-induced plasma experimental validation exercise for simulations which uses both numerical error estimation and input parameter uncertainty quanti cation to provide a direct comparison between Particle-In-Cell (PIC) plasma simulations and experiments. This approach demonstrates how careful validation can uncover missing physics in the simulation. Three di erent validation examples are shown; a vacuum space charge limited cavity, a gas lled space charge limited cavity, and a vacuum non space charge limited cavity. Two of the example are picked to show the importance of error estimation in uncovering inaccuracy/incomplete simulation models. We also report on a newly-developed numerical error estimation approach, StREEQ, which is a notable improvement to past approaches. In the StREEQ approach, a multi- tting scheme based on L 1 , L 2 , and L$\infty$ error norms and alternate weightings is used to propagate uncertainties in the relative importance of outliers and coarse/re ned discretization levels. Bootstrap sampling is used to represent the stochasticity in the response data. The resulting method appears to robustly and conservatively predict the fully-converged response within estimated numerical error bounds for stochastic simulations. The StREEQ approach is demonstrated on two related prototype electron diode problems, and preliminary results are reported for a radiation-induced plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of an implicit electromagnetic capability for a hybrid gyrokinetic ion-fluid electron model

Here we report on the development and implementation of a hybrid kinetic ion–fluid electron model for electromagnetic COGENT simulations of edge plasmas. COGENT is a finite-volume gyrokinetic code that employs a locally field-aligned coordinate system combined with a mapped multi-block grid technology to handle strongly anisotropic edge plasma turbulence. The simulation model involves the long-wavelength limit of the ion gyrokinetic equation coupled to the vorticity and Ohm's law equations for the electromagnetic field perturbations. In order to handle the fast Alfvén wave time scales, an implicit-explicit time integration approach with a physics-based preconditioner is used. The model is successfully applied to the simulations of ion-scale resistive-drift ballooning turbulence in a toroidal annulus geometry. Substantial speed-up over a fully explicit time integration approach is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Resistive magnetohydrodynamic simulations of helicity-injected startup plasmas in National Spherical Torus eXperiment

The generation of helicity-injected startup plasmas in National Spherical Torus eXperiment (NSTX), including flux surface closure, is studied using resistive-magnetohydrodynamic simulations with plasma flows, currents, ohmic heating and anisotropic thermal conduction. An injection-voltage pulse shape is used that separates the injection and closure phases allowing elucidation of the physics. The formation of an X-point near the helicity-injection gap is triggered as the injector voltage drops to zero. Near the forming X-point, magnetic pressure due to toroidal field entrained in the E × B plasma flow from the helicity-injection gap drops, allowing resistive magnetic reconnection even though the total injected current is almost constant. Where appropriate, the simulations are compared with Transient Coaxial Helicity Injection experiments in the NSTX spherical tokamak, which have demonstrated the formation of a promising candidate for non-inductive startup plasmas [Raman et al., Phys. Rev. Lett. 90, 075005 (2003)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗