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

Hard magnetohydrodynamic limit in 1/3 sawtooth like activity in LHD

The optimization of LHD discharges in inward-shifted configurations with 1/3 sawtooth like activity is an open issue. These relaxation events limit the LHD performance driving a periodic plasma deconfinement. The aim of this study is to analyze the 1/3 sawtooth like activity in plasmas with different stability properties to foreseen the best operation conditions and minimize its undesired effects. We summarize the results of several MHD simulations for plasmas with Lundquist numbers between 10{sup 5} and 10{sup 6} in the slow reconnection regime, studying the equilibria properties during the onset of a chain of 1/3 sawtooth like events. The research conclusions point out that the hard MHD limit can be reached in the inner plasma region after the onset of a strong 1/3 resonant sawtooth like event and trigger a plasma collapse. The collapse can be avoided if the system remains in the soft MHD limit, namely, in a regime with a pressure gradient and a magnetic turbulence below the critical values to drive the soft-hard MHD transition. In the soft MHD limit the system relaxations are the non resonant 1/3 sawtooth like events or a weak version of the 1/3 resonant sawtooth like events. A system relaxation in the soft MHD regime drives a minor plasma deconfinement that does not limit the LHD performance if the event periodicity is not very high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Novel data interpretation method for DIII-D divertor retarding field energy analyzer with 3-D particle-in-cell simulations

A novel data interpretation process that utilizes comprehensive particle-in-cell (PIC) simulations is developed for the new retarding field energy analyzer (RFEA) currently being constructed at DIII-D for the lower divertor using the Divertor Material Evaluation System. Furthermore, this probe is expected to survive a heat load of up to 100 MW/m 2 for up to 5 s and reliably measure the main ion temperature (T i ) on the divertor target ranging from 10 to 200 eV. These extreme conditions posed significant engineering limitations on the probe geometry, thus extensive validation work has been performed. The conventional fitting method for the RFEA I–V characteristics is based on a simplified 1-D model without considering the ion space charge inside the probe cavity and may not be sufficient for probes designed for the DIII-D divertor environment. In this article, a more realistic description of the particle propagation process within the RFEA cavity is achieved by including both 3-D geometric effects and ion space charge in the PIC simulations, and the capability to reconstruct the ion energy distribution functions is demonstrated with reasonable consistency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Plasma ignition and steady state simulations of the Linac4 H{sup −} ion source

The RF heating of the plasma in the Linac4 H{sup −} ion source has been simulated using a particle-in-cell Monte Carlo collision method. This model is applied to investigate the plasma formation starting from an initial low electron density of 10{sup 12} m{sup −3} and its stabilization at 10{sup 18} m{sup −3}. The plasma discharge at low electron density is driven by the capacitive coupling with the electric field generated by the antenna, and as the electron density increases the capacitive electric field is shielded by the plasma and induction drives the plasma heating process. Plasma properties such as e{sup −}/ion densities and energies, sheath formation, and shielding effect are presented and provide insight to the plasma properties of the hydrogen plasma.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Chemical kinetics and relaxation of non-equilibrium air plasma generated by energetic photon and electron beams

The comprehension of electromagnetic perturbations of electronic devices, due to air plasma-induced electromagnetic field, requires a thorough study on air plasma. In the aim to understand the phenomena at the origin of the formation of non-equilibrium air plasma, we simulate, using a volume average chemical kinetics model (0D model), the time evolution of a non-equilibrium air plasma generated by an energetic X-ray flash. The simulation is undertaken in synthetic air (80% N{sub 2} and 20% O{sub 2}) at ambient temperature and atmospheric pressure. When the X-ray flash crosses the gas, non-relativistic Compton electrons (low energy) and a relativistic Compton electron beam (high energy) are simultaneously generated and interact with the gas. The considered chemical kinetics scheme involves 26 influent species (electrons, positive ions, negative ions, and neutral atoms and molecules in their ground or metastable excited states) reacting following 164 selected reactions. The kinetics model describing the plasma chemistry was coupled to the conservation equation of the electron mean energy, in order to calculate at each time step of the non-equilibrium plasma evolution, the coefficients of reactions involving electrons while the energy of the heavy species (positive and negative ions and neutral atoms and molecules) is assumed remaining close to ambient temperature. It has been shown that it is the relativistic Compton electron beam directly created by the X-ray flash which is mainly responsible for the non-equilibrium plasma formation. Indeed, the low energy electrons (i.e., the non-relativistic ones) directly ejected from molecules by Compton collisions contribute to less than 1% on the creation of electrons in the plasma. In our simulation conditions, a non-equilibrium plasma with a low electron mean energy close to 1 eV and a concentration of charged species close to 10{sup 13 }cm{sup −3} is formed a few nanoseconds after the peak of X-ray flash intensity. 200 ns after the flash application, the electrons are thermalized and their concentration has decreased from about 10{sup 13 }cm{sup −3} down to about 10{sup 12 }cm{sup −3} leaving positive and negative ionised species and atomic radicals whose recombination characteristic times are much longer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

An asymptotic preserving automatic domain decomposition method for the Vlasov–Poisson–BGK system with applications to plasmas

In this work we present an efficient strategy to deal with plasma physics simulations in which localized departures from thermodynamical equilibrium are present. The method is based on the introduction of intermediate regions which allows smooth transitions between kinetic and fluid zones. In this paper we extend Domain Decomposition techniques, obtained through dynamic coupling and buffer zones, to the study of plasmas and, moreover, we combine them with Asymptotic Preserving and Asymptotically Accurate strategies for the time integration. We use a hybrid scheme in which both kinetic and fluid descriptions are considered and coupled together while the kinetic model is solved by asymptotic preserving and accurate methods, in order to guarantee high efficiency and accuracy in all regimes. The numerical scheme is validated and its performances are analyzed by numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Thermal analysis of Test Divertor Unit Scraper Element for Wendelstein 7-X

The Wendelstein-7X stellarator will go through three operational phases, with successively more energy as more complex divertors are installed. The first operational phase operated without a divertor, instead using five limiters (one in each of the five field periods), the second phase uses ten uncooled test divertor units (TDU), and the third will include an actively cooled divertor. Simulations of plasma operation have shown that the evolution of a toroidal current can lead to significant heat fluxes sweeping across sensitive regions in the divertor, possibly exceeding their qualified limit. A “Scraper Element” (SE) has been proposed to protect these sensitive divertor regions. To examine the effectiveness of this high heat-flux SE, a TDU Scraper Element (TDU-SE) is scheduled to operate with the TDU. A thermal analysis has been performed to correlate to measurements to be obtained during operation to evaluate the location and intensity of the convective plasma flux. The analysis will give an understanding of what is occurring at the plasma edge. As a result, it will also allow for a quick evaluation of the safety margin during initial, low power pulses so that the experiment can nimbly move to higher power operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Relativistic nonlinear dynamics of an intense laser beam propagating in a hot electron-positron magnetoactive plasma

The present study is devoted to investigation of the nonlinear dynamics of an intense laser beam interacting with a hot magnetized electron-positron plasma. Propagation of the intense circularly polarized laser beam along an external magnetic field is studied using a relativistic two-fluid model. A modified nonlinear Schrödinger equation is derived based on the quasi-neutral approximation, which is valid for hot plasma. Light envelope solitary waves and modulation instability are studied, for one-dimensional case. Using a three-dimensional model, spatial-temporal development of laser pulse is investigated. Occurrence of some nonlinear phenomena such as self-focusing, self-modulation, light trapping, and filamentation of laser pulse is discussed. Also the effect of external magnetic field and plasma temperature on the nonlinear evolution of these phenomena is studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Laser-driven 1 GeV carbon ions from preheated diamond targets in the break-out afterburner regime

Experimental data are presented for laser-driven carbon C{sup 6+} ion-acceleration, verifying 2D-PIC studies for multi-species targets in the Break-Out Afterburner regime. With Trident's ultra-high contrast at relativistic intensities of 5 × 10{sup 20} W/cm{sup 2} and nm-scale diamond targets, acceleration of carbon ions has been optimized by using target laser-preheating for removal of surface proton contaminants. Using a high-resolution wide angle spectrometer, carbon C{sup 6+} ion energies exceeding 1 GeV or 83 MeV/amu have been measured, which is a 40% increase in maximum ion energy over uncleaned targets. These results are consistent with kinetic plasma modeling and analytic theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonlinear structures: Cnoidal, soliton, and periodical waves in quantum semiconductor plasma

Properties and emerging conditions of various nonlinear acoustic waves in a three dimensional quantum semiconductor plasma are explored. A plasma fluid model characterized by degenerate pressures, exchange correlation, and quantum recoil forces is established and solved. Our analysis approach is based on the reductive perturbation theory for deriving the Kadomtsev-Petviashvili equation from the fluid model and solving it by using Painlevé analysis to come up with different nonlinear solutions that describe different pulse profiles such as cnoidal, soliton, and periodical pulses. The model is then employed to recognize the possible perturbations in GaN semiconductor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Detailed modeling of electron emission for transpiration cooling of hypersonic vehicles

Electron transpiration cooling (ETC) is a recently proposed approach to manage the high heating loads experienced at the sharp leading edges of hypersonic vehicles. Computational fluid dynamics (CFD) can be used to investigate the feasibility of ETC in a hypersonic environment. In this paper, a modeling approach is presented for ETC, which includes developing the boundary conditions for electron emission from the surface, accounting for the space-charge limit effects of the near-wall plasma sheath. The space-charge limit models are assessed using 1D direct-kinetic plasma sheath simulations, taking into account the thermionically emitted electrons from the surface. The simulations agree well with the space-charge limit theory proposed by Takamura et al. for emitted electrons with a finite temperature, especially at low values of wall bias, which validates the use of the theoretical model for the hypersonic CFD code. The CFD code with the analytical sheath models is then used for a test case typical of a leading edge radius in a hypersonic flight environment. The CFD results show that ETC can lower the surface temperature of sharp leading edges of hypersonic vehicles, especially at higher velocities, due to the increase in ionized species enabling higher electron heat extraction from the surface. Finally, the CFD results also show that space-charge limit effects can limit the ETC reduction of surface temperatures, in comparison to thermionic emission assuming no effects of the electric field within the sheath.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Kinetic effects on the currents determining the stability of a magnetic island in tokamaks

The role of the bootstrap and polarization currents for the stability of neoclassical tearing modes is investigated employing both a drift kinetic and a gyrokinetic approach. The adiabatic response of the ions around the island separatrix implies, for island widths below or around the ion thermal banana width, density flattening for islands rotating at the ion diamagnetic frequency, while for islands rotating at the electron diamagnetic frequency the density is unperturbed and the only contribution to the neoclassical drive arises from electron temperature flattening. As for the polarization current, the full inclusion of finite orbit width effects in the calculation of the potential developing in a rotating island leads to a smoothing of the discontinuous derivatives exhibited by the analytic potential on which the polarization term used in the modeling is based. This leads to a reduction of the polarization-current contribution with respect to the analytic estimate, in line with other studies. Other contributions to the perpendicular ion current, related to the response of the particles around the island separatrix, are found to compete or even dominate the polarization-current term for realistic island rotation frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Investigation of SF6 Alternatives in Spark Gap Switches for GWP Reduction

This primary purpose of this project was to evaluate alternative gas mixtures to sulfur hexafluoride (SF6) developed for high voltage power delivery applications for use in high voltage spark gap switches. These SF6 alternatives lower global warming potential emissions and enable improvements to the pressure-voltage design space. A combined experimental, computational, and theoretical study was used to quantify the impact of persistent breakdown products on the breakdown distribution of SF6-replacement gas mixtures. Viable SF6 replacements suitable for use in spark gap switches were studied to enable performance and agility improvements for next-generation pulsed power research relevant to national security missions. Experimental campaign included establishing parameters of switch gases as function of concentration. Various concentrations and pressures were tested for trends in breakdown voltage, repeatability, and durability, and breakdown constituents. A zero-dimensional plasma global model was used to simulate the plasma arc decay and recombination process in spark-gap switches relevant to the Z machine. Finally, a complete and consistent set of electron-neutral collision cross-sections for the novel insulating gas C4F7N is reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Fast Plasma Investigation for MMS: Simulation of the Burst Triggering System

The Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. To meet these requirements, the Fast Plasma Instrument (FPI) consists of eight (8) identical half top-hat electron sensors and eight (8) identical ion sensors and an Instrument Data Processing Unit (IDPU). The sensors (electron or ion) are grouped into pairs whose 6 degree x 180 degree fields-of-view (FOV) are set 90 degrees apart. Each sensor is equipped with electrostatic aperture steering to allow the sensor to scan a 45 degree x 180 degree fan about the its nominal viewing (0 deflection) direction. Each pair of sensors, known as the Dual Electron Spectrometer (DES) and the Dual Ion Spectrometer (DIS), occupies a quadrant on the MMS spacecraft and the combination of the eight electron/ion sensors, employing aperture steering, image the full-sky every 30-ms (electrons) and 150-ms (ions), respectively. To probe the diffusion regions of reconnection, the highest temporal/spatial resolution mode of FPI results in the DES complement of a given spacecraft generating 6.5-Mb (raised dot) per second of electron data while the DIS generates 1.1-Mb (raised dot) per second of ion data yielding an FPI total data rate of 6.6-Mb (raised dot) per second. The FPI electron/ion data is collected by the IDPU then transmitted to the Central Data Instrument Processor (CIDP) on the spacecraft for science interest ranking. Only data sequences that contain the greatest amount of temporal/spatial structure will be intelligently down-linked by the spacecraft. This requires a data ranking process known as the burst trigger system. The burst trigger system uses pseudo physical quantities to approximate the local plasma environments. As each pseudo quantity will have a different value, a set of two scaling factors is employed for each pseudo term. These pseudo quantities are then combined at the instrument, spacecraft, and observatory level leading to a final ranking of data based on expected scientific interest. Here, we present simulations of the fixed point burst trigger system for the FPI. A variety of data sets based on previous mission data as well as analytical formulations are tested. Comparisons of floating point calculations versus the fixed point hardware simulation are shown. Analysis of the potential sources of error from overflows, quantization, etc. are examined and mitigation methods are presented. Finally a series of calibration curves are presented, showing the expected error in pseudo quantities based solely on the scale parameters chosen and the expected data range. We conclude with a presentation of the current base-lined FPI burst trigger approach.

Barrie, A. C.↗
NASA NTRS2008

Fast Plasma Instrument for MMS: Simulation Results

Magnetospheric Multiscale (MMS) mission will study small-scale reconnection structures and their rapid motions from closely spaced platforms using instruments capable of high angular, energy, and time resolution measurements. The Dual Electron Spectrometer (DES) of the Fast Plasma Instrument (FPI) for MMS meets these demanding requirements by acquiring the electron velocity distribution functions (VDFs) for the full sky with high-resolution angular measurements every 30 ms. This will provide unprecedented access to electron scale dynamics within the reconnection diffusion region. The DES consists of eight half-top-hat energy analyzers. Each analyzer has a 6 deg. x 11.25 deg. Full-sky coverage is achieved by electrostatically stepping the FOV of each of the eight sensors through four discrete deflection look directions. Data compression and burst memory management will provide approximately 30 minutes of high time resolution data during each orbit of the four MMS spacecraft. Each spacecraft will intelligently downlink the data sequences that contain the greatest amount of temporal structure. Here we present the results of a simulation of the DES analyzer measurements, data compression and decompression, as well as ground-based analysis using as a seed re-processed Cluster/PEACE electron measurements. The Cluster/PEACE electron measurements have been reprocessed through virtual DES analyzers with their proper geometrical, energy, and timing scale factors and re-mapped via interpolation to the DES angular and energy phase-space sampling measurements. The results of the simulated DES measurements are analyzed and the full moments of the simulated VDFs are compared with those obtained from the Cluster/PEACE spectrometer using a standard quadrature moment, a newly implemented spectral spherical harmonic method, and a singular value decomposition method. Our preliminary moment calculations show a remarkable agreement within the uncertainties of the measurements, with the results obtained by the Cluster/PEACE electron spectrometers. The data analyzed was selected because it represented a potential reconnection event as currently published.

Figueroa-Vinas, Adolfo↗
DOE OSTI2023

Electromagnetic turbulence simulation of tokamak edge plasma dynamics and divertor heat load during thermal quench

Abstract The edge plasma turbulence and transport dynamics, as well as the divertor power loads during the thermal-quench phase of tokamak disruptions, are numerically investigated with BOUT++’s flux-driven six-field electromagnetic turbulence model. Here, transient yet intense particle and energy sources are applied at the pedestal top to mimic the plasma power drive at the edge induced by a core thermal collapse, which flattens the core temperature profile. Interesting features, such as surging of divertor heat load (up to 50 times) and broadening of heat-flux width (up to four times) on the outer-divertor target plate, are observed in the simulation, in qualitative agreement with experimental observations. The dramatic changes in divertor heat load and width are due to the enhanced plasma turbulence activities inside the separatrix. Two cross-field transport mechanisms, namely, the E × B turbulent convection and the stochastic parallel advection/conduction, are identified to play important roles in this process. First, an elevated edge pressure gradient drives instabilities and subsequent turbulence in the entire pedestal region. The enhanced turbulence not only transports particles and energy radially across the separatrix via the E × B convection, which causes the initial divertor heat-load burst, but it also induces amplified magnetic fluctuation B ˜ . Once themagnetic fluctuation is large enough to break the magnetic flux surface, magnetic flutter effect provides an additional radial transport channel. In the late stage of our simulation, | B ˜ r / B 0 | reaches to 10 −4 level that completely breaks magnetic flux surfaces such that stochastic field lines are directly connecting pedestal top plasma to the divertor target plates or first wall, further contributing to the divertor heat-flux width broadening.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

The study of the Poincare-Lighthill-Kuo method by using the particle-in-cell simulation method in a dusty plasma

The application scope of the Poincare-Lighthill-Kuo (PLK) method is suggested by using the Particle-in-cell (PIC) numerical method to study head-on collision of two solitary waves. Comparisons between the numerical results from PIC simulations and the analytical ones from the PLK method indicate that the both are in good agreement with each other. The dependence of the phase shifts after the head-on collision on both amplitudes of two solitary waves is given from our PIC method. It is found that the phase shifts depended on the amplitude of both waves. The maximum amplitude during the colliding process is approximately equal to the sum of both amplitudes for the small amplitude solitary waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1971

Measurements of plasma parameters in a simulated thermionic converter

Cesium-filled thermionic energy converters are examined as candidate electrical energy sources in spacecraft requiring tens to hundreds of kilowatts of electric power. The high operating temperatures necessary for a large specific power and high efficiency inevitably impose stringent constraints on the converter fabrication to achieve the desired reliability of the power system. The converter physics for reducing operating temperatures and cesium plasma losses are studied to achieve high reliability without sacrificing the power performance of the converters. Various cesium parameters which affect the converter performance are: (1) electron temperatures, (2) plasma ion densities, and (3) electric potential profiles. These were investigated using a Langmuir probe in a simulated converter. The parameters were measured in different cesium discharge modes.

Shimada, K.↗
DOE OSTI2015

Physics of the interaction of ultra intense laser pulses with cold collisional plasma using large scale kinetic simulations

We present a set of 2D collisional particle-in-cell simulations of the interaction of ultra-intense laser pulses with over-dense cold collisional plasmas. The size of these simulations is about 100 times as large as those previously published. This allows studying the transport of energetic particles on time scale of the order of 400 fs without perturbations due to the influence of boundary effects and performing a very detailed analysis of the physics of the transport. We confirm the existence of a threshold in intensity close to the relativistic threshold above which the beam of energetic particles diverges when it penetrates the cold plasma. We also study the applicability of Ohm's law to compute the electric field, which is the method commonly used in hybrid codes. The heating of the cold plasma is then studied and we show that half of the heating is anomalous, i.e., not given by standard Joule effect. We discuss the previously published results in the light of these new simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗