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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Non-linear dynamics of compound sawteeth in tokamaks

Compound sawteeth is studied with the XTOR-2F code. Non-linear full 3D magnetohydrodynamic simulations show that the plasma hot core is radially displaced and rotates during the partial crash, but is not fully expelled out of the q = 1 surface. Partial crashes occur when the radius of the q = 1 surface exceeds a critical value, at fixed poloidal beta. This critical value depends on the plasma elongation. The partial crash time is larger than the collapse time of an ordinary sawtooth, likely due to a weaker diamagnetic stabilization. This suggests that partial crashes result from a competition between destabilizing effects such as the q = 1 radius and diamagnetic stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Elimination of numerical Cherenkov instability in flowing-plasma particle-in-cell simulations by using Galilean coordinates

Particle-in-cell (PIC) simulations of relativistic flowing plasmas are of key interest to several fields of physics (including, e.g., laser-wakefield acceleration, when viewed in a Lorentz-boosted frame) but remain sometimes infeasible due to the well-known numerical Cherenkov instability (NCI). In this article, we show that, for a plasma drifting at a uniform relativistic velocity, the NCI can be eliminated by simply integrating the PIC equations in Galilean coordinates that follow the plasma (also sometimes known as comoving coordinates ) within a spectral analytical framework. The elimination of the NCI is verified empirically and confirmed by a theoretical analysis of the instability. Moreover, it is shown that this method is applicable both to Cartesian geometry and to cylindrical geometry with azimuthal Fourier decomposition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Global Field Aligned Mesh and Gyrokinetic Field Solver in a Tokamak Edge Geometry (Final Technical Report)

This project was a collaboration between researchers at the California Institute of Technology and the University of California, Irvine to investigate the utility of a global field-aligned mesh and gyrokinetic field solver for simulations of the tokamak plasma edge region. Mesh generation software from UC Irvine was tested with specific tokamak edge magnetic geometry scenarios and the quality of the meshes and the solutions to the gyrokinetic Poisson equation were evaluated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Plasma viscosity elevations with simulated weightlessness

A hypothesis correlating an increase in blood viscosity during bed rest to a decrease in aerobic capacity during simulated weightlessness is tested. Eight human subjects were studied on the sixth day of bed rest during two consecutive 10-d bed rest periods separated by a 14-d recovery interval designed to simulate the flight-layover schedule of Shuttle astronauts. Plasma viscosity and volume were measured, together with maximal aerobic capacity (VO2max). An increase in hematocrit, plasma protein, and fibrinogen concentrations was found, contributing to an elevation in plasma viscosity. VO2max decreased significantly in the first, but not the second bed rest cycle, and though many individuals exhibited a decrease in plasma volume and aerobic capacity coupled with elevated plasma viscosity, correlations between these variables were lacking. It is concluded that the decrease in VO2max observed following simulated weightlessness cannot be attributed to alterations in muscle blood flow resulting from increased blood viscosity.

Martin, D. G.↗
NASA NTRS1989

Space plasma contractor research, 1988

Results of experiments conducted on hollow cathode-based plasma contractors are reported. Specific tests in which attempts were made to vary plasma conditions in the simulated ionospheric plasma are described. Experimental results showing the effects of contractor flowrate and ion collecting surface size on contactor performance and contactor plasma plume geometry are presented. In addition to this work, one-dimensional solutions to spherical and cylindircal space-charge limited double-sheath problems are developed. A technique is proposed that can be used to apply these solutions to the problem of current flow through elongated double-sheaths that separate two cold plasmas. Two conference papers which describe the essential features of the plasma contacting process and present data that should facilitate calibration of comprehensive numerical models of the plasma contacting process are also included.

Williams, John D.↗
NASA NTRS1997

An Analysis of Closure Mechanisms in the Plasma Wake of the TSS-1R Satellite

Collisionless Plasma Expansion (CPE), also known as "plasma expansion into a vacuum," in its most fundamental form, is the process by which a plasma expands into a void, or region highly depleted of particles. In CPE, the expansion is driven by the highly mobile electron constituent of the plasma as it moves across the density gradient and into the void region. A bi-polar electric field is set up between this rapidly expanding electron front and the massive ions, which have a low thermal speed and lag behind in the region of the plasma-void interface. These ions are quickly accelerated by the expansion electric field and can, theoretically, approach the thermal speed of the electrons. CPE is a fundamental and a very robust process in that it only requires that a sharp density gradient exist in a plasma. It has already been observed to exist in various plasmas ranging over five orders of magnitude in density. The range of phenomena in which CPE has been observed includes the closure of plasma wakes created in simulated space plasmas in the laboratory; it was found to be the process controlling the refilling of the wake created by the Space Shuttle in the ionospheric plasma; and recent in situ data, obtained by the WIND spacecraft as it passed through the wake of the earth's moon last year, have shown CPE to be involved in the closure of the Lunar wake. CPE is also expected to influence the morphology and physics of the solar wind interactions with Mercury and Mars, and has been suggested as a potential factor in such wide ranging phenomena as the acceleration of plasma away from the coma of comets and the energization of ions that underlies the upward expansion of the polar wind into the earth's magnetosphere. The Tethered Satellite System (TSS)-1R data provides an opportunity to test this process in situ under semi-controlled conditions. The available data shows the intensity, drift energy and angle of inclination to the wake of the converging ion streams at a down stream distance of one satellite radius (80 CM).

Stone, N. H.↗
NASA NTRS2002

Magnetized Target Fusion

Magnetized target fusion (MTF) is under consideration as a means of building a low mass, high specific impulse, and high thrust propulsion system for interplanetary travel. This unique combination is the result of the generation of a high temperature plasma by the nuclear fusion process. This plasma can then be deflected by magnetic fields to provide thrust. Fusion is initiated by a small traction of the energy generated in the magnetic coils due to the plasma's compression of the magnetic field. The power gain from a fusion reaction is such that inefficiencies due to thermal neutrons and coil losses can be overcome. Since the fusion reaction products are directly used for propulsion and the power to initiate the reaction is directly obtained from the thrust generation, no massive power supply for energy conversion is required. The result should be a low engine mass, high specific impulse and high thrust system. The key is to successfully initiate fusion as a proof-of-principle for this application. Currently MSFC is implementing MTF proof-of-principle experiments. This involves many technical details and ancillary investigations. Of these, selected pertinent issues include the properties, orientation and timing of the plasma guns and the convergence and interface development of the "pusher" plasma. Computer simulations of the target plasma's behavior under compression and the convergence and mixing of the gun plasma are under investigation. This work is to focus on the gun characterization and development as it relates to plasma initiation and repeatability.

Steven T. Griffin↗
DOE OSTI2014

Verification of gyrokinetic particle simulation of current-driven instability in fusion plasmas. I. Internal kink mode

The gyrokinetic toroidal code (GTC) capability has been extended for simulating internal kink instability with kinetic effects in toroidal geometry. The global simulation domain covers the magnetic axis, which is necessary for simulating current-driven instabilities. Here, GTC simulation in the fluid limit of the kink modes in cylindrical geometry is verified by benchmarking with a magnetohydrodynamic eigenvalue code. Gyrokinetic simulations of the kink modes in the toroidal geometry confirum that ion kinetic effects greatly reduce the growth rate even when the banana orbit width is much smaller than the radial width of the perturbed current layer at the mode rational surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Numerical simulations of output pulse extraction from a high-power microwave compressor with a plasma switch

Numerical simulations of the process of electromagnetic energy release from a high-power microwave pulse compressor comprising a gas-filled cavity and interference switch were carried out. A microwave plasma discharge in a rectangular waveguide H-plane tee was modeled with the use of the fully electromagnetic particle-in-cell code MAGIC. The gas ionization, plasma evolution, and interaction with RF fields accumulated within the compressor were simulated using different approaches provided by the MAGIC code: particle-in-cell approach accounting for electron-neutral collisions, gas conductivity model based on the concept of mobility, and hybrid modeling. The dependences of the microwave output pulse peak power and waveform on parameters that can be controlled in experiments, such as an external ionization rate, RF field amplitude, and background gas pressure, were investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Formation of v-shaped potentials

The V-shaped potential structures formed by the injection of a non-neutral electron current into a cold background plasma were simulated numerically. The injection disturbs the initial quasi-neutral plasma, leading to the excitation of strong turbulences which heat the plasma. This leads to expulsion of the plasma from the simulation region. Due to ambipolar electric fields the current injection is interrupted and the initial background plasma is extracted from the system. A particle composition with the characteristics of the two plasma reservoirs now represents the plasma in the simulation region. The interaction of the electron beam with this plasma excites turbulences of smaller amplitudes. A nearly constant time averaged potential drop with nonstationary distribution develops across the system. Single and multiple double layers may form for the duration of one ion plasma period.

Thiemann, H.↗
DOE OSTI2015

Measurements of the populations of metastable and resonance levels in the plasma of an RF capacitive discharge in argon

The behavior of the populations of two metastable and two lower resonance levels of argon atoms in the plasma of an RF capacitive discharge was studied. The populations were measured by two methods: the method of emission self-absorption and the method based on measurements of the intensity ratios of spectral lines. It is shown that the populations of resonance levels increase with increasing power deposited in the discharge, whereas the populations of metastable levels is independent of the RF power. The distribution of the populations over energy levels is not equilibrium under these conditions. The population kinetics of argon atomic levels in the discharge plasma is simulated numerically. The distribution function of plasma electrons recovered from the measured populations of atomic levels and numerical simulations is found to be non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Parameterizing the Morse potential for coarse-grained modeling of blood plasma

Multiscale simulations of fluids such as blood represent a major computational challenge of coupling the disparate spatiotemporal scales between molecular and macroscopic transport phenomena characterizing such complex fluids. In this paper, a coarse-grained (CG) particle model is developed for simulating blood flow by modifying the Morse potential, traditionally used in Molecular Dynamics for modeling vibrating structures. The modified Morse potential is parameterized with effective mass scales for reproducing blood viscous flow properties, including density, pressure, viscosity, compressibility and characteristic flow dynamics of human blood plasma fluid. The parameterization follows a standard inverse-problem approach in which the optimal micro parameters are systematically searched, by gradually decoupling loosely correlated parameter spaces, to match the macro physical quantities of viscous blood flow. The predictions of this particle based multiscale model compare favorably to classic viscous flow solutions such as Counter-Poiseuille and Couette flows. It demonstrates that such coarse grained particle model can be applied to replicate the dynamics of viscous blood flow, with the advantage of bridging the gap between macroscopic flow scales and the cellular scales characterizing blood flow that continuum based models fail to handle adequately.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Fluid modeling of radio frequency and direct currents in a biased magnetized plasma

This model aims at simulating a magnetized plasma column connected on one side to a probe and on the other side to an ICRH (Ion Cyclotron Resonant Heating) antenna. This is a double probe modeling of a plasma flux tube exchanging perpendicular currents: rf polarization current and dc currents (inertia, viscous, and friction) perpendicular to the magnetic field. A self consistent solution for the rectified potential and the density is obtained under the assumptions of flute hypothesis, inertialess electrons, and no collision in parallel direction. The main effect of rf biasing on the antenna side is to shift the IV characteristic so that the floating potential can be increased up to ln(I{sub 0}(eV{sub rf}/(k{sub B}T{sub e}))), with I{sub 0} the modified Bessel function of the first kind. On the contrary, the effect of dc currents is to decrease the plasma potential and the sheath potential which can be lower than 3k{sub B}T{sub e}/e or even be negative. Experimental characteristics are well matched by the 1D fluid code and exhibit very high negative currents (more than 30 j{sub i} the ion saturation current) for high positive biasing of the probe and for long flux tube (10 m). The non-saturation of the electron current is here due to high transverse transport enhanced by convective fluxes and dc currents able to bring an amount of density around the biased flux tube. During comparisons with experiments, the floating potential measured by a reciprocating probe is recovered by the code revealing that for a 120 V measured peak potential on the probe, the rf potential on the ICRH antenna is twice this value. Finally, the density profile can be flattened or steepened as a function of the transverse dc current direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Fundamental processes in the expansion, energization, and coupling of single- and multi-Ion plasmas in space: Laboratory simulation experiments

We have conducted a laboratory investigation into the physics of plasma expansions and their associated energization processes. We studied single- and multi-ion plasma processes in self-expansions, and included light and heavy ions and heavy/light mixtures to encompass the phenomenological regimes of the solar and polar winds and the AMPTE and CRRES chemical release programs. The laboratory experiments provided spatially-distributed time-dependent measurements of total plasma density, temperature, and density fluctuation power spectra with the data confirming the long-theorized electron energization process in an expanding cloud - a result that was impossible to determine in spaceborne experiments (as e.g., in the CRRES program). These results provided the missing link in previous laboratory and spaceborne programs. confirming important elements in our understanding of such solar-terrestrial processes as manifested in expanding plasmas in the solar wind (e.g., CMES) and in ionospheric outflow in plasmaspheric fluctuate refilling after a storm. The energization signatures were seen in an entire series of runs that varied the ion species (Ar', Xe', Kr' and Ne'), and correlative studies included spectral analyses of electrostatic waves collocated with the energized electron distributions. In all cases wave energies were most intense during the times in which the suprathermal populations were present, with wave intensity increasing with the intensity of the suprathermal electron population. This is consistent with theoretical expectations wherein the energization process is directly attributable to wave particle interactions. No resonance conditions were observed, in an overall framework in which the general wave characteristics were broadband with power decreasing with increasing frequency.

Szuszczewicz, E. P.↗
DOE OSTI2025

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗
DOE OSTI2014

Three-dimensional MHD simulation of the Caltech plasma jet experiment: first results

Magnetic fields are believed to play an essential role in astrophysical jets with observations suggesting the presence of helical magnetic fields. Here, we present three-dimensional (3D) ideal MHD simulations of the Caltech plasma jet experiment using a magnetic tower scenario as the baseline model. Magnetic fields consist of an initially localized dipole-like poloidal component and a toroidal component that is continuously being injected into the domain. This flux injection mimics the poloidal currents driven by the anode-cathode voltage drop in the experiment. The injected toroidal field stretches the poloidal fields to large distances, while forming a collimated jet along with several other key features. Detailed comparisons between 3D MHD simulations and experimental measurements provide a comprehensive description of the interplay among magnetic force, pressure, and flow effects. In particular, we delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms. With suitably chosen parameters that are derived from experiments, the jet in the simulation agrees quantitatively with the experimental jet in terms of magnetic/kinetic/inertial energy, total poloidal current, voltage, jet radius, and jet propagation velocity. Specifically, the jet velocity in the simulation is proportional to the poloidal current divided by the square root of the jet density, in agreement with both the experiment and analytical theory. This work provides a new and quantitative method for relating experiments, numerical simulations, and astrophysical observation, and demonstrates the possibility of using terrestrial laboratory experiments to study astrophysical jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE POWER OF IMAGING: CONSTRAINING THE PLASMA PROPERTIES OF GRMHD SIMULATIONS USING EHT OBSERVATIONS OF Sgr A*

Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons in the low density accreting plasma. This poses a significant challenge in predicting accretion flow images and spectra from first principles. Because of this, simplified emission models have often been used, with widely different configurations (e.g., disk- versus jet-dominated emission), and were able to account for the observed spectral properties of accreting black holes. Exploring the large parameter space introduced by such models, however, requires significant computational power that exceeds conventional computational facilities. In this paper, we use GRay, a fast graphics processing unit (GPU) based ray-tracing algorithm, on the GPU cluster El Gato, to compute images and spectra for a set of six general relativistic magnetohydrodynamic simulations with different magnetic field configurations and black-hole spins. We also employ two different parametric models for the plasma thermodynamics in each of the simulations. We show that, if only the spectral properties of Sgr A* are used, all 12 models tested here can fit the spectra equally well. However, when combined with the measurement of the image size of the emission using the Event Horizon Telescope, current observations rule out all models with strong funnel emission, because the funnels are typically very extended. Our study shows that images of accretion flows with horizon-scale resolution offer a powerful tool in understanding accretion flows around black holes and their thermodynamic properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Laser produced plasma sources for nanolithography—Recent integrated simulation and benchmarking

Photon sources for extreme ultraviolet lithography (EUVL) are still facing challenging problems to achieve high volume manufacturing in the semiconductor industry. The requirements for high EUV power, longer optical system and components lifetime, and efficient mechanisms for target delivery have narrowed investigators towards the development and optimization of dual-pulse laser sources with high repetition rate of small liquid tin droplets and the use of multi-layer mirror optical system for collecting EUV photons. We comprehensively simulated laser-produced plasma sources in full 3D configuration using 10–50 μm tin droplet targets as single droplets as well as, for the first time, distributed fragmented microdroplets with equivalent mass. The latter is to examine the effects of droplet fragmentation resulting from the first pulse and prior to the incident second main laser pulse. We studied the dependence of target mass and size, laser parameters, and dual pulse system configuration on EUV radiation output and on atomic and ionic debris generation. Our modeling and simulation included all phases of laser target evolution: from laser/droplet interaction, energy deposition, target vaporization, ionization, plasma hydrodynamic expansion, thermal and radiation energy redistribution, and EUV photons collection as well as detail mapping of photons source size and location. We also simulated and predicted the potential damage to the optical mirror collection system from plasma thermal and energetic debris and the requirements for mitigating systems to reduce debris fluence. The debris effect on mirror collection system is analyzed using our three-dimensional ITMC-DYN Monte Carlo package. Modeling results were benchmarked against our CMUXE laboratory experimental studies for the EUV photons production and for debris and ions generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗