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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Nonlinear excitation of low-n harmonics in reduced magnetohydrodynamic simulations of edge-localized modes

Nonlinear simulations of the early edge-localized mode (ELM) phase based on a typical type-I ELMy ASDEX Upgrade discharge have been carried out using the reduced MHD code JOREK. The analysis is focused on the evolution of the toroidal Fourier spectrum. It is found that during the nonlinear evolution, linearly subdominant low-n Fourier components, in particular the n = 1, grow to energies comparable with linearly dominant harmonics. A simple model is developed, based on the idea that energy is transferred among the toroidal harmonics via second order nonlinear interaction. The simple model reproduces and explains very well the early nonlinear evolution of the toroidal spectrum in the JOREK simulations. Furthermore, it is shown for the n = 1 harmonic, that its spatial structure changes significantly during the transition from linear to nonlinearly driven growth. The rigidly growing structure of the linearly barely unstable n = 1 reaches far into the plasma core. In contrast, the nonlinearly driven n= 1 has a rigidly growing structure localized at the plasma edge, where the dominant toroidal harmonics driving the n = 1 are maximal and in phase. The presented quadratic coupling model might explain the recent experimental observation of strong low-n components in magnetic measurements [Wenninger et al., “Non-linear magnetic perturbations during edge localized modes in TCV dominated by low n mode components,” Nucl. Fusion (submitted)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Laser-plasma interactions with a Fourier-Bessel particle-in-cell method

A new spectral particle-in-cell (PIC) method for plasma modeling is presented and discussed. In the proposed scheme, the Fourier-Bessel transform is used to translate the Maxwell equations to the quasi-cylindrical spectral domain. In this domain, the equations are solved analytically in time, and the spatial derivatives are approximated with high accuracy. In contrast to the finite-difference time domain (FDTD) methods, that are used commonly in PIC, the developed method does not produce numerical dispersion and does not involve grid staggering for the electric and magnetic fields. These features are especially valuable in modeling the wakefield acceleration of particles in plasmas. The proposed algorithm is implemented in the code PLARES-PIC, and the test simulations of laser plasma interactions are compared to the ones done with the quasi-cylindrical FDTD PIC code CALDER-CIRC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

A Coupled Plasma and Sheath Model for High Density Reactors

We present a coupled plasma and collisionless; sheath model for the simulation of high density plasma processing reactors. Due to inefficiencies in numerical schemes and the resulting computational burden, a coupled multidimensional plasma and sheath simulation has not been possible model for gas mixtures and high density reactors of practical interest. In this work we demonstrate that with a fully implicit algorithm and a refined computational mesh, a self-consistent plasma and sheath simulation is feasible. We discuss the details of the model equations, the importance of ion inertia, and the resulting sheath profiles for argon and chlorine plasmas. We find that at low operating pressures (10-30 mTorr), the charge separation occurs only within a 0.5 mm layer near the surface in a 300 mm inductively coupled plasma etch reactor. A unified model eliminates the use of off-line or loosely coupled sheath models with simplifying assumptions which generally lead to uncertainties in ion flux and sheath electrical properties.

Deepak, Bose↗
DOE OSTI2020

Plasma decay in hydrocarbons and hydrocarbon- and H 2 O-containing mixtures excited by high-voltage nanosecond discharge at elevated gas temperatures

Plasma decay after a high-voltage nanosecond discharge was experimentally and numerically studied in pure hydrocarbons (C 2 H 6 and C 3 H 8 ), and H 2 O:N 2 and C 3 H 8 :O 2 mixtures for pressures in the range 2–4 Torr and gas temperatures from 300 to 600 K. In a stoichiometric C 3 H 8 :O 2 mixture, plasma decay was also studied in a repetitively pulsed discharge for varing numbers of discharge pulses (varying degrees of fuel oxidation). The rate of plasma decay was determined from the temporal evolution of electron density measured using the microwave interferometer. It was observed that gas heating to 600 K leads to a decrease in the rate of plasma decay in all cases. The effect of heating on plasma decay was most profound in the H 2 O:N 2 mixture (after a single discharge pulse) and in the C 3 H 8 :O 2 mixture for high degrees of fuel oxidation. A kinetic scheme was developed to numerically simulate the plasma decay in hydrocarbons and combustible mixtures. Numerical analysis showed that, under the conditions studied, plasma decay was controlled by dissociative electron recombination with simple molecular and cluster ions. Gas heating led to a decrease in the rate of the electron-ion recombination and the rate of conversion of molecular ions to cluster ions. As a result, the gas temperature increase caused a decrease in the fraction of cluster ions for which the recombination coefficients are an order of magnitude higher than the recombination coefficients for molecular ions. The influence of gas heating on the decrease of the amount of cluster ions was more important when the ion composition was dominated by hydrated H 3 O + (H 2 O) k ions. The rates of the formation of these ions are extremely sensitive to any variations in gas temperature. A result, in agreement with our observations, gas heating led to an anomalous decrease in the rates of plasma decay in the H 2 O:N 2 mixture, as well as in the C 3 H 8 :O 2 mixture when H 2 O molecules were produced due to fuel oxidation.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2015

THIN CURRENT SHEETS AND ASSOCIATED ELECTRON HEATING IN TURBULENT SPACE PLASMA

Intermittent structures, such as thin current sheets, are abundant in turbulent plasmas. Numerical simulations indicate that such current sheets are important sites of energy dissipation and particle heating occurring at kinetic scales. However, direct evidence of dissipation and associated heating within current sheets is scarce. Here, we show a new statistical study of local electron heating within proton-scale current sheets by using high-resolution spacecraft data. Current sheets are detected using the Partial Variance of Increments (PVI) method which identifies regions of strong intermittency. We find that strong electron heating occurs in high PVI (>3) current sheets while no significant heating occurs in low PVI cases (<3), indicating that the former are dominant for energy dissipation. Current sheets corresponding to very high PVI (>5) show the strongest heating and most of the time are consistent with ongoing magnetic reconnection. This suggests that reconnection is important for electron heating and dissipation at kinetic scales in turbulent plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Effect of laser light polarization on generation of relativistic ion beams driven by an ultraintense laser

The effect of laser light polarization on properties of proton and carbon ion beams generated from a CH target irradiated by a 130 fs laser pulse of ultra-relativistic intensity ({approx}10{sup 22}-10{sup 23} W/cm{sup 2}) is investigated using particle-in-cell simulations. It is shown that only circular light polarization ensures the production of quasi-monoenergetic relativistic beams of both protons and carbon ions from such a target while using the linear one results in the generation of quasi-monoenergetic protons accompanied with carbon ions of complex and broad energy spectrum. The influence of the target thickness and laser intensity on the ion energy spectrum and the laser-ions energy conversion efficiency is examined.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Electrochemical Corrosion Behavior of Plasma-Sprayed FeCrNiMoCBSi Amorphous/Nanocrystalline Coatings in Simulated Seawater Medium

A novel FeCrNiMoCBSi coating was fabricated on F22 marine steel by plasma spraying. The coating microstructure and composition were characterized by scanning electron microscope, transmission electron microscope, x-ray diffraction, differential scanning calorimetry, and x-ray photoelectron spectroscopy. The results showed that the coating consisted of amorphous and nanocrystalline phases. Additionally, the corrosion resistance in simulated seawater was evaluated by electrochemical corrosion methods. It was found that the coating exhibited an excellent resistance to pitting corrosion. For long-time immersion tests, the corrosion resistance was gradually decreased with the increasing through-porosity due to the permeation of the corrosive solution into the coating. The coating lost its efficacy through three stages, namely the protection stage, the permeation stage, and the failure stage.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Development of electron cyclotron emission imaging system on the HL-2A tokamak

A 2D electron cyclotron emission imaging (ECEI) system has been developed for measurement of electron temperature fluctuations in the HL-2A tokamak. It is comprised of a front-end 24 channel heterodyne imaging array with a tunable RF range spanning 75–110 GHz, and a set of back-end ECEI electronics that together generate 24 × 8 = 192 channel images of the 2nd harmonic X-mode electron cyclotron emission from the HL-2A plasma. The simulated performance of the local oscillator (LO) optics and radio frequency (RF) optics is presented, together with the laboratory characterization results. The Gaussian beams from the LO optics are observed to properly cover the entire detector array. The ECE signals from the plasma are mixed with the LO signal in the array box, then delivered to the electronics system by low-loss microwave cables, and finally to the digitizers. The ECEI system can achieve temporal resolutions of ∼μs, and spatial resolutions of 1 cm (radially) and 2 cm (poloidally)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

THE PLASMA ENVIRONMENT IN COMETS OVER A WIDE RANGE OF HELIOCENTRIC DISTANCES: APPLICATION TO COMET C/2006 P1 (MCNAUGHT)

On 2007 January 12, comet C/2006 P1 (McNaught) passed its perihelion at 0.17 AU. Abundant remote observations offer plenty of information on the neutral composition and neutral velocities within 1 million kilometers of the comet nucleus. In early February, the Ulysses spacecraft made an in situ measurement of the ion composition, plasma velocity, and magnetic field when passing through the distant ion tail and the ambient solar wind. The measurement by Ulysses was made when the comet was at around 0.8 AU. With the constraints provided by remote and in situ observations, we simulated the plasma environment of Comet C/2006 P1 (McNaught) using a multi-species comet MHD model over a wide range of heliocentric distances from 0.17 to 1.75 AU. The solar wind interaction of the comet at various locations is characterized and typical subsolar standoff distances of the bow shock and contact surface are presented and compared to analytic solutions. We find the variation in the bow shock standoff distances at different heliocentric distances is smaller than the contact surface. In addition, we modified the multi-species model for the case when the comet was at 0.7 AU and achieved comparable water group ion abundances, proton densities, plasma velocities, and plasma temperatures to the Ulysses/SWICS and SWOOPS observations. We discuss the dominating chemical reactions throughout the comet-solar wind interaction region and demonstrate the link between the ion composition near the comet and in the distant tail as measured by Ulysses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

MAGNETOHYDRODYNAMIC WAVES AND CORONAL HEATING: UNIFYING EMPIRICAL AND MHD TURBULENCE MODELS

We present a new global model of the solar corona, including the low corona, the transition region, and the top of the chromosphere. The realistic three-dimensional magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating MHD Alfven wave turbulence. We assume this turbulence and its nonlinear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the nonlinear cascade rate degrades in strongly anisotropic (imbalanced) turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes, from which the fast solar wind originates. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and STEREO EUVI instruments for the Carrington rotation 2107.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Final Report for Grant DE-FG02-93ER54215

The research supported by this grant led to many important publications in Physical Review Letters, Physics of Plasmas and Nuclear Fusion. Two invited talks at the American Physical Society Division of Plasma Physics, oral presentations at the Sherwood Fusion Theory Conference and MHD workshops. The most important achievements include: MHD equilibrium code FLOW for tokamak equilibria with arbitrary poloidal and toroidal flow; Two fluid tokamak equilibria with arbitrary flows; 2D Compressible MHD code to simulate tokamak plasma evolution with poloidal flow and the MHD edge pedestal; Resistive Wall Mode fluid theory in the presence of flow and differentially rotating (liquid) walls; Kinetic theory of the Resistive Wall Mode; Generalized kinetic energy principle for Resistive Wall Mode stability and the RWM kinetic stability code MISK; Kinetic theory of the m=1 mode for ITER including trapped thermal ions and electrons and fishbone oscillations driven by tangential beam injection; Two-fluid Lawson ignition criterion for tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

A magnetospheric magnetic field model with flexible current systems driven by independent physical parameters

A tilt-dependent magnetic field model of the Earth's magnetosphere with variable magnetopause standoff distance is presented. Flexible analytic representations for the ring and cross-tail currents, each composed of the elements derived from the Tsyganenko and Usmanov (1982) model, are combined with the fully shielded vacuum dipole configurations of Voigt (1981). Although the current sheet does not warp in the y-z plane, changes in the shape and position of the neutral sheet with dipole tilt are consistent with both MHD equilibrium theory and observations. In addition, there is good agreement with observed Delta B profiles and the average equatorial contours of magnetic field magnitude. While the dipole field is rigorously shielded within the defined magnetopause, the ring and cross-tails currents are not similarly confined, consequently, the model's region of validity is limited to the inner magnetosphere. The model depends on four independent external parameters. We present a simple but limited method of simulating several substorm related magnetic field changes associated with the disrupion of the near-Earth cross-tail current sheet and collapse of the midnight magnetotail field region. This feature further facilitates the generation of magnetic field configuration time sequences useful in plasma convection simulations of real magnetospheric events.

Hilmer, Robert V.↗
NASA NTRS2005

Particle acceleration, magnetic field generation, and emission in relativistic pair jets

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Recent simulations show that the Weibel instability created by relativistic pair jets is responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet propagating through an ambient plasma with and without initial magnetic fields. The growth rates of the Weibel instability depends on the distribution of pair jets. The Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. This instability is also responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The jitter radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗
DOE OSTI2015

A comparative study of capacitively coupled HBr/He, HBr/Ar plasmas for etching applications: Numerical investigation by fluid model

Fluid model has been applied to perform a comparative study of hydrogen bromide (HBr)/He and HBr/Ar capacitively coupled plasma discharges that are being used for anisotropic etching process. This model has been used to identify the most dominant species in HBr based plasmas. Our simulation results show that the neutral species like H and Br, which are the key player in chemical etching, have bell shape distribution, while ions like HBr{sup +}, Br{sup +}, which play a dominant rule in the physical etching, have double humped distribution and show peaks near electrodes. It was found that the dilution of HBr by Ar and/or He results in an increase in electron density and electron temperature, which results in more ionization and dissociation and hence higher densities of neutral and charged species can be achieved. The ratio of positive ion flux to the neutral flux increases with an increase in additive gas fraction. Compare to HBr/He plasma, the HBr/Ar plasma shows a maximum change in the ion density and flux and hence the etching rate can be considered in the ion-assisted and in the ion-flux etch regime in HBr/Ar discharge. The densities of electron and other dominant species in HBr/Ar plasma are higher than those of HBr/He plasma. The densities and fluxes of the active neutrals and positive ions for etching and subsequently chemical etching versus physical sputtering in HBr/Ar and HBr/He plasmas discharge can be controlled by tuning gas mixture ratio and the desire etching can be achieved.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Validation of NSFsim as a Grad-Shafranov equilibrium solver at DIII-D

Plasma shape is a significant factor that must be considered for any Fusion Pilot Plant (FPP) as it has significant consequences for plasma stability and core confinement. A new simulator, NSFsim, has been developed based on a historically successful code, DINA [1], offering tools to simulate both transport and plasma shape. Specifically, NSFsim is a free boundary equilibrium and transport solver and has been configured to match the properties of the DIII-D tokamak. This paper is focused on validating the Grad-Shafranov (GS) solver of NSFsim by analyzing its ability to recreate the plasma shape, the poloidal flux distribution, and the measurements of the simulated diagnostic signals originating from flux loops and magnetic probes in DIII-D. Five different plasma shapes are simulated to show the robustness of NSFsim to different plasma conditions; these shapes are Lower Single Null (LSN), Upper Single Null (USN), Double Null (DN), Inner Wall Limited (IWL), and Negative Triangularity (NT). The NSFsim results are compared against real measured signals, magnetic profile fits from EFIT [2], and another plasma equilibrium simulator, GSevolve [3]. EFIT reconstructions of shots are readily available at DIII-D, but GSevolve was manually ran by us to provide simulation data to compare against.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Kinetic Physics of Homogeneous Turbulence in Collisionless Plasmas

We present results of large scale, fully kinetic 3D simulations of decaying plasma turbulence that span a range of scales from tens of ion gyroradii to the electron Debye length. A significant fraction of the dissipated energy is transferred to ions, in contrast to the scenario where the dominant dissipation is due to Landau damping of wave-like fluctuations. Intense current sheets are formed in the simulation with thicknesses over a range of scales approximately between electron and ion inertial scales. The rate of ion heating is well correlated with intermittency level. Geometrical analysis of the current sheets suggest that magnetic reconnection is active in the present simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗