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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Magnetosheath-ionspheric plasma interactions in the cusp/cleft. 2: Mesoscale particle simulations

Ionospheric plasma flowing out from the cusp can be an important source of plasma to the magnetosphere. One source of free energy that can drive this outflow is the injection of magnetosheath plasma into the cusp. Two-dimensional (three velocity) mesoscale particle simulations are used to investigate the particle dynamics in the cusp during southward interplanetary magnetic field. This mesoscale model self-consistently incorporates (1) global influences such as the convection of plasma across the cusp, the action of the mirror force, and the injection of the magnetosheath plasma, and (2) wave-particle interactions which produce the actual coupling between the magnetosheath and ionospheric plasmas. It is shown that, because the thermal speed of the electrons is higher than the bulk motion of the magnetosheath plasma, an upward current is formed on the equatorward edge of the injection region with return currents on either side. However, the poleward return currents are the stronger due to the convection and mirroring of many of the magnetosheath electrons. The electron distribution in this latter region evolves from upward directed streams to single-sided loss cones or possibly electron conics. The ion distribution also shows a variety of distinct features that are produced by spatial and/or temporal effects associated with varying convection patterns and wave-particle interactions. On the equatorward edge the distribution has a downflowing magnetosheath component and an upflowing cold ionospheric component due to continuous convection of ionospheric plasma into the region. In the center of the magnetosheath region, heating from the development of an ion-ion streaming instability causes the suppression of the cold ionospheric component and the formation of downward ionospheric streams. Further poleward there is velocity filtering of ions with low pitch angles, so that the magnetosheath ions develop a ring-beam distribution and the ensuing wave instabilities generate downward ionospheric conics. These downward ionospheric components are eventually turned by the mirror force, leading to the production of upward conics at elevated energies throughout the region.

Winglee, R. M.↗
NASA NTRS1992

Relativistic, perpendicular shocks in electron-positron plasmas

One-dimensional particle-in-cell plasma simulations are used to examine the mechanical structure and thermalization properties of collisionless relativistic shock waves in electron-positron plasmas. Shocks propagating perpendicularly to the magnetic field direction are considered. It is shown that these shock waves exist, and that they are completely parameterized by the ratio of the upstream Poynting flux to the upstream kinetic energy flux. The way in which the Rankine-Hugoniot shock jump conditions are modified by the presence of wave fluctuations is shown, and they are used to provide a macroscopic description of these collisionless shock flows. The results of a 2D simulation that demonstrates the generality of these results beyond the assumption of the 1D case are discussed. It is suggested that the thermalization mechanism is the formation of a synchrotron maser by the coherently reflected particles in the shock front. Because the downstream medium is thermalized, it is argued that perpendicular shocks in pure electron-positron plasmas are not candidates as nonthermal particle accelerators.

Gallant, Yves A.↗
DOE OSTI2018

Complexity reduction effects on transient, atomic plasmas

Plasma simulations involving collisional-radiative processes present many computational challenges because of the wide range in scales spanned by the governing set of plasma physics equations. These numerical challenges are primarily caused by limitations in resources such as memory and computational time. Collisional-radiative models deliver a part of the challenge because of the large number of levels and transitions that are associated with atomic plasmas. Multiple complexity reduction techniques have been introduced in the past to overcome this obstacle, including the quasi-steady-state solution, uniform grouping, and the recently-devised Boltzmann grouping. These reduction schemes were tested against full numerical simulations of a simplified, time-dependent argon model, with and without an external Planckian radiation field. Lastly, results will show the advantages and disadvantages of using complexity reduction techniques when employed in various plasma regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Monte Carlo simulation of reentry plasmas

Attention is given to the treatment of ionization and plasma effects in the direct simulation Monte Carlo method. The requirements for accurate modeling of reentry plasmas are discussed along with the difficulties these requirements present. The current method for modeling such plasmas is reviewed and an alternative method is presented. Both methods are applied to the flow of a 10 km/s shock wave in air at 0.1 torr; a flowfield directly relevant to the projected aeroassisted orbital transfer vehicle. The results are compared and the differences between the methods are discussed.

Carlson, Ann B.↗
DOE OSTI2018

Continuously differentiable PIC shape functions for triangular meshes

In this study, a new class of continuously-differentiable shape functions is developed and applied to two-dimensional electrostatic PIC simulation on an unstructured simplex (triangle) mesh. It is shown that troublesome aliasing instabilities are avoided for cold plasma simulation in which the Debye length is as small as 0.01 cell sizes. These new shape functions satisfy all requirements for PIC particle shape. They are non-negative, have compact support, and partition unity. They are given explicitly by cubic expressions in the usual triangle logical (areal) coordinates. The shape functions are not finite elements because their structure depends on the topology of the mesh, in particular, the number of triangles neighboring each mesh vertex. Nevertheless, they may be useful as approximations to solution of other problems in which continuity of derivatives is required or desired.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Numerical and Experimental Investigation of Turbulent Transport Control via Shaping of Radial Plasma Flow Profiles (Final Technical Report)

The strong coupling between the different physical variables involved in the plasma transport phenomenon and the high complexity of its dynamics call for a model-based, multivariable approach to profile control where those predictive models could be exploited. The overall objective of this project has been to extend the existing body of work by investigating numerically and experimentally active control of unstable fluctuations, including fully developed turbulence and the associated cross-field particle transport, via manipulation of flow profiles in a magnetized laboratory plasma device. Fluctuations and particle transport can be monitored by an array of electrostatic probes, and E x B flow profiles can be controlled via a set of biased concentric ring electrodes that terminate the plasma column. The goals of the proposed research have been threefold: i- to develop a predictive code to simulate plasma transport in the linear HELCAT (HELicon-CAThode) plasma device at the University of New Mexico (UNM), where the experimental component of the proposed research has been carried out; ii- to establish the feasibility of using advanced model-based control algorithms to control cross-field turbulence-driven particle transport through appropriate manipulation of radial plasma flow profiles, iii- to investigate the fundamental nonlinear dynamics of turbulence and transport physics. Lehigh University (LU), including Prof. Eugenio Schuster and one full-time graduate student, has been primarily responsible for control-oriented modeling and model-based control design. Undergraduate students have also participated in this project through the National Science Foundation Research Experience for Undergraduate (REU) program. The main goal of the LU Plasma Control Group has been to study the feasibility of controlling turbulence-driven transport by shaping the radial poloidal flow profile (i.e., by controlling flow shear) via biased concentric ring electrodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Perpendicular heating of electrons by upper hybrid waves generated by a ring distribution

Satellite observations of electron conical distributions with enhanced fluxes just outside the loss cone suggest that telectrons have been heated perpendicularly to the magnetic field in the mid-altitude polar magnetosphere. To understand electron conical distributions, plasma simulations are conducted to examine an upper hybrid wave instability of a ring electron distribution perpendicular to the magnetic field in a cold electron background. The simulations indicate that both the cold and ring distributions are heated perpendicularly during the saturation stage. From the plasma data, a ring distribution can be identified as a trapped distribution function with an enhancement near 90-deg pitch angle in the phase space density plot. It is suggested that the ring distribution might provide an additional free energy source for generating upper hybrid waves associated with electron conical events.

Lin, C. S.↗
DOE OSTI2026

Collaborative Research: Enabling multi-scale studies of magnetic reconnection with interpretable data-driven models

The development of accurate reduced descriptions and improved closures for magnetic reconnection is an important and a long‐standing challenge in plasma physics. The four‐fluid approach, and associated closures, that were investigated have the potential to improve the accuracy of plasma fluid models, capturing physical effects which would otherwise require a kinetic description. If successful, this approach could have an important impact for the modeling of laboratory and space plasmas. The major goals of this project were to develop new machine learning (ML) tools based on sparse and symbolic regression techniques, and to extract interpretable and generalizable reduced models (e.g., in the form of partial differential equations - PDEs) from data generated by first principles plasma simulations. Preserving interpretability of such data‐driven models is key to addressing the long‐standing theoretical and numerical challenges. Prior proof‐of‐principle studies have demonstrated the enormous potential of this approach, by recovering the well‐established hierarchy of plasma equations (from Vlasov to MHD) from data produced by particle‐in‐cell (PIC) simulations. Our goal in this project was to extend and apply these new tools to construct better kinetic closures for magnetic reconnection; to derive better models of particle injection and acceleration by this fundamental plasma process; and to use this understanding to accelerate the development of multi‐scale plasma algorithms. While our immediate focus was on the problem of magnetic reconnection, the tools that were will developed are general and applicable to other areas of plasma physics, and more broadly to many‐body phenomena. We anticipate that the development of these multi‐scale models will have a significant impact across different areas of plasma science, from fusion to space and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Multi-Scale Investigation of Sheared Flows In Magnetized Plasmas

Flows parallel and perpendicular to magnetic fields in a plasma are important phenomena in many areas of plasma science research. The presence of these spatially inhomogeneous flows is often associated with the stability of the plasma. In fusion plasmas, these sheared flows can be stabilizing while in space plasmas, these sheared flows can be destabilizing. Because of this, there is broad interest in understanding the coupling between plasma stability and plasma flows. This research project has engaged in a study of the plasma response to spatially inhomogeneous plasma flows using three different experimental devices: the Auburn Linear Experiment for Instability Studies (ALEXIS) and the Compact Toroidal Hybrid (CTH) stellarator devices at Auburn University, and the Space Plasma Simulation Chamber (SPSC) at the Naval Research Laboratory. This work has shown that there is a commonality of the plasma response to sheared flows across a wide range of plasma parameters and magnetic field geometries. The goal of this multi-device, multi-scale project is to understand how sheared flows established by the same underlying physical mechanisms lead to different plasma responses in fusion, laboratory, and space plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Exploring the Acceleration Mechanisms for Particle Injection and Power-law Formation during Transrelativistic Magnetic Reconnection

Magnetic reconnection in the relativistic and transrelativistic regimes is able to accelerate particles to hard power-law energy spectra f ∝ γ -p (approaching p = 1). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during magnetic reconnection in the transrelativistic regime of a proton–electron plasma. While earlier work in this parameter regime has focused on the effects of electric field parallel to the local magnetic field on the particle injection (from thermal energy to the lower-energy bound of the power-law spectrum), here we examine the roles of both parallel and perpendicular electric fields to gain a more complete understanding on the injection process and further development of a power-law spectrum. We show that the parallel electric field does contribute significantly to particle injection, and is more important in the initial phase of magnetic reconnection. However, as the simulation proceeds, the acceleration by the perpendicular electric field becomes more important for particle injection and completely dominates the acceleration responsible for the high-energy power-law spectrum. This holds robustly, in particular for longer reconnection times and larger systems, i.e., in simulations that are more indicative of the processes in astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Scaled-down particle-in-cell simulation of cathode plasma expansion in magnetically insulated coaxial diode

The expansion of cathode plasma in magnetically insulated coaxial diode (MICD) is investigated in theory and particle-in-cell (PIC) simulation. The temperature and density of the cathode plasma are about several eV and 10{sup 13}–10{sup 16 }cm{sup −3}, respectively, and its expansion velocity is of the level of few cm/μs. Through hydrodynamic theory analysis, expressions of expansion velocities in axial and radial directions are obtained. The characteristics of cathode plasma expansion have been simulated through scaled-down PIC models. Simulation results indicate that the expansion velocity is dominated by the ratio of plasma density other than the static electric field. The electric field counteracts the plasma expansion reverse of it. The axial guiding magnetic field only reduces the radial transport coefficients by a correction factor, but not the axial ones. Both the outward and inward radial expansions of a MICD are suppressed by the much stronger guiding magnetic field and even cease.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Self-organized confinement in whole-device modeling of laboratory magnetospheres

Turbulent self-organization driven by global ideal interchange modes in a dipole-confined plasma is explored with self-consistent, whole-plasma simulations using a flux-tube averaged magnetohydrodynamic model in dipole magnetic geometry. We show the existence of robust particle pinch driven by ideal interchange-mode fluctuations, in which the particles are transported up the density gradient. It is found that the plasma profiles in a dipole field spontaneously relax to a marginally stable state as centrally peaked pressure and density are created by global interchange-mode transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

High-average-power, 100-Hz-repetition-rate, tabletop soft-x-ray lasers at sub-15-nm wavelengths

Efficient excitation of dense plasma columns at 100-Hz repetition rate using a tailored pump pulse profile produced a tabletop soft-x-ray laser average power of 0.1 mW at = 13.9 nm and 20 W at = 11.9 nm from transitions of Ni-like Ag and Ni-like Sn, respectively. Lasing on several other transitions with wavelengths between 10.9 and 14.7 nm was also obtained using 0.9-J pump pulses of 5-ps duration from a compact diode-pumped chirped pulse amplification Yb:YAG laser. Hydrodynamic and atomic plasma simulations show that the pump pulse profile, consisting of a nanosecond ramp followed by two peaks of picosecond duration, creates a plasma with an increased density of Ni-like ions at the time of peak temperature that results in a larger gain coefficient over a temporally and spatially enlarged space leading to a threefold increase in the soft-x-ray laser output pulse energy. The high average power of these compact soft-x-ray lasers will enable applications requiring high photon flux. These results open the path to milliwatt-average-power tabletop soft-x-ray lasers.

Reagan, Brendon↗
NASA NTRS1992

The dynamics of low-beta plasma clouds as simulated by a three-dimensional, electromagnetic particle code

A 3D electromagnetic and relativistic particle simulation code is used to investigate the dynamics of low-beta plasma clouds moving perpendicular to an ambient magnetic field in vacuum and in a background plasma. The simulations show the formation of the space charge sheaths at the sides of the cloud with the associated polarization electric field which facilitate the cross-field propagation, as well as the sheaths at the front and rear end of the cloud caused by the larger ion Larmor radius, which allows ions to move ahead and lag behind the electrons as they gyrate. It is found that, in a background plasma, electron and ion sheaths expand along the magnetic field at the same rate, whereas in vacuum the electron sheath expands much faster than the ion sheath. Sheath clouds are accelerated up to relativistic energies, which indicates that artificial plasma clouds released in the ionosphere or magnetosphere may generate optical emissions as energetic sheath electrons scatter in the upper atmosphere.

Neubert, T.↗
DOE OSTI2024

Spatial decay of electromagnetic waves from hypervelocity impact plasmas

Due to the high kinetic energy of impactors in space, hypervelocity impacts on spacecraft can ablate spacecraft material and produce plasmas. These plasmas have been observed in ground-based experiments to produce electromagnetic pulses (EMPs) and are suspected to be a cause of spacecraft anomalies. The initial EMP that is produced after plasma formation is replicated using a 2D electromagnetic plasma simulation, consisting of a discontinuous Galerkin, particle-in-cell (DG-PIC) framework. The spontaneous formation of EMPs from an initial quasineutral plasma is consistent with previous works and show behavior similar to that of waves emitted by antennas. The EMP emits from the plasma with a frequency less than the plasma frequency associated with the point of maximum density, showing the importance of the density gradient, and the far-field region of the “plasma antenna” is estimated. From the proportionality between plasma frequency and EMP frequency, there is a direct correlation between plasma density and radiation. Finally, as the plasma density decreases, the EMP frequency decreases, which suggests that several observable phenomena, such as optical flashes, radio frequency emission, and microwave emission, all result from the same plasma-related mechanism within a single hypervelocity impact.

42 ENGINEERING↗
DOE OSTI2014

Radio imaging of a type IVM radio burst on the 14th of August 2010

Propagating coronal mass ejections (CMEs) are often accompanied by burst signatures in radio spectrogram data. We present Nançay Radioheliograph observations of a moving source of broadband radio emission, commonly referred to as a type IV radio burst (type IVM), which occurred in association with a CME on the 14th of August 2010. The event was well observed at extreme ultraviolet (EUV) wavelengths by SDO/AIA and PROBA2/SWAP, and by the STEREO SECCHI and SOHO LASCO white light (WL) coronagraphs. The EUV and WL observations show the type IVM source to be cospatial with the CME core. The observed spectra is well fitted by a power law with a negative slope, which is consistent with optically thin gyrosynchrotron emission. The spectrum shows no turn over at the lowest Nançay frequencies. By comparing simulated gyrosynchrotron spectra with Nançay Radioheliograph observations, and performing a rigorous parameter search we are able to constrain several key parameters of the underlying plasma. Simulated spectra found to fit the data suggest a nonthermal electron distribution with a low energy cutoff of several tens to 100 keV, with a nonthermal electron density in the range 10{sup 0}-10{sup 2} cm{sup –3}, in a magnetic field of a few Gauss. The nonthermal energy content of the source is found to contain 0.001%-0.1% of the sources thermal energy content. Furthermore, the energy loss timescale for this distribution equates to several hours, suggesting that the electrons could be accelerated during the CME initiation or early propagation phase and become trapped in the magnetic structure of the CME core without the need to be replenished.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Numerical simulation of a plasma actuator based on ion transport

Two-dimensional numerical simulation of ion transport and flow around a single dielectric barrier discharge plasma actuator (PA) is performed. Spatial distributions of ions and electrons as well as their time evolution are obtained by solving the transport equations of monovalent positive ions, monovalent negative ions, and electrons. Voltage and frequency of the driving alternating-current signal are assumed to be 8 kV and 5 kHz, respectively. Special focus is laid upon the effect of voltage gradient dV/dt on the magnitude of the body force. The validity of steady force models often used in flow simulation is also examined. The simulation results show that the magnitude of the body force induced by the PA increases as the voltage gradient dV/dt increases and its increase rate becomes milder at higher voltage. The mechanism of body force generation is explained from the time evolution of number density fields of ions and electrons. A comparison between flow simulations using a time-resolved body force and its time-averaged counterpart demonstrates that the time-averaged model gives sufficiently accurate results when the time scale of the flow is more than 30 times greater than that of the PA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗