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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Simulating Sources of Superstorm Plasmas

We evaluated the contributions to magnetospheric pressure (ring current) of the solar wind, polar wind, auroral wind, and plasmaspheric wind, with the surprising result that the main phase pressure is dominated by plasmaspheric protons. We used global simulation fields from the LFM single fluid ideal MHD model. We embedded the Comprehensive Ring Current Model within it, driven by the LFM transpolar potential, and supplied with plasmas at its boundary including solar wind protons, polar wind protons, auroral wind O+, and plasmaspheric protons. We included auroral outflows and acceleration driven by the LFM ionospheric boundary condition, including parallel ion acceleration driven by upward currents. Our plasmasphere model runs within the CRCM and is driven by it. Ionospheric sources were treated using our Global Ion Kinetics code based on full equations of motion. This treatment neglects inertial loading and pressure exerted by the ionospheric plasmas, and will be superceded by multifluid simulations that include those effects. However, these simulations provide new insights into the respective role of ionospheric sources in storm-time magnetospheric dynamics.

Fok, Mei-Ching↗
DOE OSTI2019

Experiments and Simulations of Hypervelocity Impact Plasmas: Final Report

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, 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 the first 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. 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. Furthermore, we combined the collected data with results from numerical simulations using particle-in-cell (PIC) and smoothed particle hydrodynamics (SPH) techniques. Using results from both experimental measurements and simulations, we addressed the following outstanding science questions: 1) What are the properties of hypervelocity impact plasmas? 2) What conditions cause impact plasmas to create RF emission and what is the power dependence on frequency? 3) How does the ambient environment around the spacecraft, including the background magnetic and electric fields and plasma density, influence the hypervelocity impact plasmas?

42 ENGINEERING↗
DOE OSTI2021

Analysis of Design Alternatives of Actively Cooled RF Window for MPEX

The Materials Plasma Exposure eXperiment (MPEX) has been designed as a linear plasma divertor simulator in order to address plasma material interaction (PMI) science for next-generation fusion devices. It will have the capability to test neutron irradiated samples with plasma fluxes of greater than 10 24 m -2 s -1 . It is expected to operate steady state for up to 10 6 s to consider PMI affects through reactor end of life. The conceptual design of MPEX was completed in 2019, with preliminary design having begun in 2020. The plasma source for MPEX is a helicon antenna, where the energized helical antenna sits outside of the vacuum in order to minimize impurities in the plasma. It is expected to receive up to 200 kW of continuous power, and so the antenna and the window must be actively cooled. The water-cooled copper antenna has been operated at full power on the Proto-MPEX device (which is a test facility to demonstrate the plasma source and heating systems). The water-cooled window, however, is a novel component that must meet numerous competing requirements. It requires a low dielectric loss to allow the Radio Frequency (RF) power to create the plasma within the vacuum boundary. It must be structurally robust to handle the significant heat flux from the plasma and any heat from dielectric coupling. It must be compatible with the coolant (preferably water). It requires a vacuum seal that minimizes impurities into the plasma and does not compromise the structural integrity of the window. Two window designs have been tested. Results from these tests, where temperatures are measured and heat fluxes inferred from infrared camera data, have been correlated with thermal-structural simulations. When these simulations are extrapolated to the full power steady-state heat fluxes that are expected in MPEX, the designs do not appear to have the necessary structural robustness. This study explores design alternatives for the MPEX helicon antenna window, presents analysis results for several of the alternatives, and shows a viable solution that satisfies the requirements for MPEX operation.

Tipton Jr, Joseph↗
NASA NTRS2010

Integration Assessment of Visiting Vehicle Induced Electrical Charging of the International Space Station Structure

The International Space Station (ISS) undergoes electrical charging in low Earth orbit (LEO) due to positively biased, exposed conductors on solar arrays that collect electrical charges from the space plasma. Exposed solar array conductors predominately collect negatively charged electrons and thus drive the metal ISS structure electrical ground to a negative floating potential (FP) relative to plasma. This FP is variable in location and time as a result of local ionospheric conditions. ISS motion through Earth s magnetic field creates an addition inductive voltage up to 20 positive and negative volts across ISS structure depending on its attitude and location in orbit. ISS Visiting Vehicles (VVs), such as the planned Orion crew exploration vehicle, contribute to the ISS plasma charging processes. Upon physical contact with ISS, the current collection properties of VVs combine with ISS. This is an ISS integration concern as FP must be controlled to minimize arcing of ISS surfaces and ensure proper management of extra vehicular activity crewman shock hazards. This report is an assessment of ISS induced charging from docked Orion vehicles employing negatively grounded, 130 volt class, UltraFlex (ATK Space Systems) solar arrays. To assess plasma electron current collection characteristics, Orion solar cell test coupons were constructed and subjected to plasma chamber current collection measurements. During these tests, coupon solar cells were biased between 0 and 120 V while immersed in a simulated LEO plasma. Tests were performed using several different simulated LEO plasma densities and temperatures. These data and associated theoretical scaling of plasma properties, were combined in a numerical model which was integrated into the Boeing Plasma Interaction Model. It was found that the solar array design for Orion will not affect the ISS FP by more than about 2 V during worst case charging conditions. This assessment also motivated a trade study to determine acceptable plasma electron current levels that can be collected by a single or combined fleet of ISS-docked VVs.

Kramer, Leonard↗
DOE OSTI2022

Eulerian simulations of electrostatic waves in plasmas with a single sign of charge

An Eulerian, numerical simulation is used to model the launching of plasma waves in a non-neutral plasma that is confined in a Penning–Malmberg trap. The waves are launched by applying an oscillating potential to an electrically isolated sector at one end of the conducting cylinder that bounds the confinement region and are received by another electrically isolated sector at the other end of the cylinder. The launching of both Trivelpiece–Gould waves and electron acoustic waves is investigated. Adopting a stratagem, the simulation captures essential features of the finite length plasma, while retaining the numerical advantages of a simulation employing periodic spatial boundary conditions. As a benchmark test of the simulation, the results for launched Trivelpiece–Gould waves of small amplitude are successfully compared to a linearized analytic solution for these fluctuations.

Physics↗
DOE OSTI2025

Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

A First-Principles Model Based on Saturation of the Electron Cyclotron Drift Instability for Electron Transport in Hydrodynamics Simulations of Hall Thruster Plasmas

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Lopez Ortega, Alejandro↗
DOE OSTI2025

SciDAC ISEP: Integrated Simulation of Energetic Particles in Burning Plasmas

The objective of the SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP) is to improve physics understanding of energetic particle (EP) confinement and EP interactions with burning thermal plasmas through large-scale simulations. The ISEP center will develop a multiscale and multiphysics ISEP framework for a predictive capability of EP physics and deliver an EP module incorporating both first-principles simulations and high fidelity reduced transport models to the fusion whole device modeling (WDM) project. The ISEP framework will enable us to perform long time, global kinetic simulations of EP physics in burning plasmas, by utilizing the full power of the next generation supercomputers. Our research and development activities will build on fruitful collaborations with computer scientists and applied mathematicians to offer enabling technologies for performance scalability, portability, solvers, coupling for integration with the fusion WDM project, and long-term preservation of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Magnetogasdynamic compression of a coaxial plasma accelerator flow for micrometeoroid simulation.

A new configuration of a coaxial plasma accelerator with self-energized magnetic compressor coil attached is described. It is shown that the circuit may be treated theoretically by analyzing an equivalent circuit mesh. The results obtained from the theoretical analysis are shown to compare favorably with the results measured experimentally. Using this accelerator configuration, glass beads of 125-micron diameter were accelerated to velocities as high as 11 km/sec, while 700-micron-diam glass beads were accelerated to velocities as high as 5 km/sec. The velocities are within the hypervelocity regime of meteroids.

Igenbergs, E. B.↗
NASA NTRS1974

Magnetogasdynamic compression of a coaxial plasma accelerator flow for micrometeoroid simulation

A new configuration of a coaxial plasma accelerator with self-energized magnetic compressor coil attached is described. It is shown that the circuit may be treated theoretically by analyzing an equivalent circuit mesh. The results obtained from the theoretical analysis compare favorably with the results measured experimentally. Using this accelerator configuration, glass beads of 125 micron diameter were accelerated to velocities as high as 11 kilometers per second, while 700 micron diameter glass beads were accelerated to velocities as high as 5 kilometers per second. The velocities are within the hypervelocity regime of meteoroids.

Igenbergs, E. B.↗
NASA NTRS1992

Titan's organic chemistry: Results of simulation experiments

Recent low pressure continuous low plasma discharge simulations of the auroral electron driven organic chemistry in Titan's mesosphere are reviewed. These simulations yielded results in good accord with Voyager observations of gas phase organic species. Optical constants of the brownish solid tholins produced in similar experiments are in good accord with Voyager observations of the Titan haze. Titan tholins are rich in prebiotic organic constituents; the Huygens entry probe may shed light on some of the processes that led to the origin of life on Earth.

Sagan, Carl↗
DOE OSTI2020

Electron-Only Reconnection in Plasma Turbulence

Hybrid-Vlasov–Maxwell simulations of magnetized plasma turbulence including non-linear electron-inertia effects in a generalized Ohm's law are presented. When fluctuation energy is injected on scales sufficiently close to ion-kinetic scales, the ions efficiently become de-magnetized and electron-scale current sheets largely dominate the distribution of the emerging current structures, in contrast to the usual picture, where a full hierarchy of structure sizes is generally observed. These current sheets are shown to be the sites of electron-only reconnection (e-rec), in which the usual electron exhausts are unaccompanied by ion outflows and which are in qualitative agreement with those recently observed by MMS in the Earth's turbulent magnetosheath, downstream of the bow shock. Some features of the e-rec phenomenology are shown to be consistent with an electron magnetohydrodynamic description. Simulations suggest that this regime of collisionless reconnection may be found in turbulent systems where plasma processes, such as micro-instabilities and/or shocks, overpower the more customary turbulent cascade by directly injecting energy close to the ion-kinetic scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Core plasma physics basis and its impacts on the FNSF

The FNSF core plasma physics is reportede on with detailed analysis to establish the basis for the reference configuration, R = 4.8 m, a = 1.2 m, Ip = 7.87 MA, and B T = 7.5 T, q 95 = 6, f BS = 0.52, β N < 2.7 established in Ref. [1]. Central solenoid (CS) and poloidal field (PF) coils are far from the plasma as in a power plant, and acceptable coil currents are determined for the rampup and flattop fiducial states. Time-dependent free-boundary plasma evolution simulations demonstrate that the FNSF plasma can be established, ramped up, and relaxed into flattop, including vertical stabilizers, internal feedback coils and feedback control on plasma current, position, and shape. A range of density (n o / = 1.3–1.5) and temperature (T o / = 2.2–2.7) profiles are examined, indicating that energy confinement of H 98 = 1.1–1.2 is required to provide 100% non-inductive plasma current in the FNSF. GLF23 theory based transport model predicted lower energy confinement of H 98 ~0.6–0.85. The EPED analysis shows that the pedestal temperature ranges from 4.0–4.7 keV for pedestal densities of 1.7–1.0 × 10 20 /m 3 . The n = 1 kink stability shows no-wall beta limits, using the pressure and current profiles associated with the transport and current drive sources, ranging from β N ~2.25–2.55 depending on l i . A conducting wall can extend these limits by 10–40% depending on li and wall location. At the lower beta’s of the reference plasma, a combination of 50 MW of NB, 30 MW of LH, 20 MW of ICRF, 20 MW of EC, and bootstrap current, are found to provide 100% of the plasma current with a stable current profile. Impacts on the FNSF of plasma physics are discussed and R&D challenges are recorded.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Power-transfer and fixed-point analysis of sawtooth simulations of a current-carrying stellarator

Power-transfer and fixed-point analysis of previous NIMROD simulations (Roberds et al. , Phys. Plasmas , vol. 23, issue 9, 2016, 092513) improved the understanding of the effect of 3D (non-axisymmetric equilibrium) magnetic fields on sawtooth oscillations in the Compact Toroidal Hybrid (CTH) experiment. Computing the locations of order-1 fixed points, their Greene's residues, and local values for the rotational transform results in a description of CTH sawteeth consistent with Kadomtsev's model. A power-transfer analysis quantifies the distribution of energy among toroidal Fourier modes and their nonlinear interactions. The Lorentz power transfer drives sawtooth growth, and it is unambiguously interpreted as the flow of energy from toroidal mode $n'$ to mode $n$ , catalysed by $\boldsymbol {B}_{n-n'}$ . It has been reported previously that the CTH sawtooth frequency increases with the 3D field strength. This is attributed to an increased growth rate of the internal kink that drives sawtooth oscillations. Here, 3D fields remove energy from the kink, eliminating the possibility that these fields are an additional energy source that drives growth. Instead, 3D fields catalyse energy transfer from large-to-small scales, where magnetic reconnection is stronger. It is proposed that this energy transfer increases the reconnection rate at small scales, which is consistent with the increased growth rate observed at higher 3D field strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Edge equilibrium code for tokamaks

The edge equilibrium code (EEC) described in this paper is developed for simulations of the near edge plasma using the finite element method. It solves the Grad-Shafranov equation in toroidal coordinate and uses adaptive grids aligned with magnetic field lines. Hermite finite elements are chosen for the numerical scheme. A fast Newton scheme which is the same as implemented in the equilibrium and stability code (ESC) is applied here to adjust the grids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Student Research Awards - Sherwood Fusion Theory Conference

Student research awards were given for outstanding student presentations at the Annual International Sherwood Fusion Theory Conference, held May 1-3rd 2017 in Annapolis, MD. The goal of controlled thermonuclear fusion research is to develop a safe, environmentally benign, and essentially inexhaustible source of energy for the future. The Annual International Sherwood Fusion Theory Conference is organized for educational and scientific purposes in support of fusion energy and plasma physics research. Specifically, the objective of the Annual Conference shall be the communication of recent results in the physics of controlled fusion energy. This conference is conducted in accordance with the International Council of Scientific Unions guidelines on free circulation of scientists. Papers are solicited on the following topics for the International Sherwood Fusion Theory Conference: Plasma Properties, Equilibrium, Stability, and Transport; Physics of Plasma Edge and Divertor Region; Plasma Production and Heating; Computer Simulation of Plasmas. Student research awards at Sherwood serve several purposes. First, student research awards generally help recruit and retain top students in fusion energy sciences, because they highlight excellence in research and help to develop the CVs of young scientists. Second, the Sherwood Fusion Theory Conference has a long history of student research awards, and continuing this tradition helps to maintain the integrity and impact of past awards. Third, research awards in theory are especially important because they highlight not only individual students, but the entire theory field in fusion energy sciences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Numerical simulations of a nonequilibrium argon plasma in a shock-tube experiment

A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.

Cambier, Jean-Luc↗
DOE OSTI2018

Path Integral Monte Carlo Simulations of Warm Dense Matter and Plasmas

New path integral Monte Carlo simulation (PIMC) techniques will be developed and applied to derive the equation of state (EOS) for the regime of warm dense matter and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 10,000 K, this method is not applicable at much higher temperature where electronic excitations become important because the number of partially occupied electronic orbitals reaches intractably large numbers and, more importantly, the use of zero-temperature exchange-correlation functionals introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it has commonly been assumed such methods can only be applied to elements without core electrons like hydrogen and helium, we recently showed how to extend PIMC to heavier elements by performing the first PIMC simulations of carbon and water plasmas [Driver, Militzer, Phys. Rev. Lett. 108 (2012) 115502]. Here we propose to continue this important development to extend the reach of PIMC simulations to yet heavier elements and also lower temperatures. The goal is to provide a robust first-principles simulation method that can accurately and efficiently study materials with excited electrons at solid-state densities in order to access parts of the phase diagram such the regime of warm dense matter and plasmas where so far only more approximate, semi-analytical methods could be applied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗