Search NASA⌕ Search

SEARCH · Search NASA

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

RAM‐SCB simulations of electron transport and plasma wave scattering during the October 2012 “double‐dip” storm

Abstract Mechanisms for electron injection, trapping, and loss in the near‐Earth space environment are investigated during the October 2012 “double‐dip” storm using our ring current‐atmosphere interactions model with self‐consistent magnetic field (RAM‐SCB). Pitch angle and energy scattering are included for the first time in RAM‐SCB using L and magnetic local time (MLT)‐dependent event‐specific chorus wave models inferred from NOAA Polar‐orbiting Operational Environmental Satellites (POES) and Van Allen Probes Electric and Magnetic Field Instrument Suite and Integrated Science observations. The dynamics of the source (approximately tens of keV) and seed (approximately hundreds of keV) populations of the radiation belts simulated with RAM‐SCB is compared with Van Allen Probes Magnetic Electron Ion Spectrometer observations in the morning sector and with measurements from NOAA 15 satellite in the predawn and afternoon MLT sectors. We find that although the low‐energy ( E < 100 keV) electron fluxes are in good agreement with observations, increasing significantly by magnetospheric convection during both SYM‐H dips while decreasing during the intermediate recovery phase, the injection of high‐energy electrons is underestimated by this mechanism throughout the storm. Local acceleration by chorus waves intensifies the electron fluxes at E ≥50 keV considerably, and RAM‐SCB simulations overestimate the observed trapped fluxes by more than an order of magnitude; the precipitating fluxes simulated with RAM‐SCB are weaker, and their temporal and spatial evolutions agree well with POES/Medium Energy Proton and Electron Detectors data.

58 GEOSCIENCES↗
DOE OSTI2015

Comprehensive comparisons of geodesic acoustic mode characteristics and dynamics between Tore Supra experiments and gyrokinetic simulations

In a dedicated collisionality scan in Tore Supra, the geodesic acoustic mode (GAM) is detected and identified with the Doppler backscattering technique. Observations are compared to the results of a simulation with the gyrokinetic code GYSELA. We found that the GAM frequency in experiments is lower than predicted by simulation and theory. Moreover, the disagreement is higher in the low collisionality scenario. Bursts of non harmonic GAM oscillations have been characterized with filtering techniques, such as the Hilbert-Huang transform. When comparing this dynamical behaviour between experiments and simulation, the probability density function of GAM amplitude and the burst autocorrelation time are found to be remarkably similar. In the simulation, where the radial profile of GAM frequency is continuous, we observed a phenomenon of radial phase mixing of the GAM oscillations, which could influence the burst autocorrelation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Computing the Paschen curve for argon with speed-limited particle-in-cell simulation

Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The 1D3V simulations modeled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional PIC methods computationally slow. To decrease this timescale difference and speed up computation, here we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully kinetic treatment of dynamics and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Optimisation of confinement in a fusion reactor using a nonlinear turbulence model

We report the confinement of heat in the core of a magnetic fusion reactor is optimised using a multidimensional optimisation algorithm. For the first time in such a study, the loss of heat due to turbulence is modelled at every stage using first-principles nonlinear simulations which accurately capture the turbulent cascade and large-scale zonal flows. The simulations utilise a novel approach, with gyrofluid treatment of the small-scale drift waves and gyrokinetic treatment of the large-scale zonal flows. A simple near-circular equilibrium with standard parameters is chosen as the initial condition. The figure of merit, fusion power per unit volume, is calculated, and then two control parameters, the elongation and triangularity of the outer flux surface, are varied, with the algorithm seeking to optimise the chosen figure of merit. A twofold increase in the plasma power per unit volume is achieved by moving to higher elongation and strongly negative triangularity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Temporal and spatial profiles of emission intensities in atmospheric pressure helium plasma jet driven by microsecond pulse: Experiment and simulation

A needle-circular electrode structure helium plasma jet driven by microsecond pulsed power is studied. Spatially resolved emission results show that the emission intensity of He(3{sup 3}S{sub 1}) line decreases monotonically along the axial direction, while those of N{sub 2}(C{sup 3}Π{sub u}), N{sub 2}{sup +}(B{sup 2}∑{sup +}{sub u}), and O(3p{sup 5}P) reach their maxima at 3 cm, 2.6 cm, and 1.4 cm, respectively. The plasma plume of the four species shows different characteristics: The N{sub 2} emission plume travels at a fast speed along the entire plasma jet; the N{sub 2}{sup +} emission plume is composed of a bright head and relatively weak tail and travels a shorter distance than the N{sub 2} emission plume; the He emission plume travels at a slower speed for only a very short distance; propagation of the O emission plume is not observed. Results of calculation of radiation fluxes emitted by positive streamers propagating along helium plasma jets are presented. It is shown, in agreement with the results of the present experiment and with other available experimental data, that the intensities of radiation of N{sub 2}(C{sup 3}Π{sub u}) molecules and He(3{sup 3}S{sub 1}) atoms vary with time (along the plasma jet) quite differently. The factors resulting in this difference are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Effects of eletron heating on the current driven electrostatic ion cyclotron instability and plasma transport processes along auroral field lines

Fluid simulations of the plasma along auroral field lines in the return current region have been performed. It is shown that the onset of electrostatic ion cyclotron (EIC) related anomalous resistivity and the consequent heating of electrons leads to a transverse ion temperature that is much higher than that produced by the current driven EIC instability (CDICI) alone. Two processes are presented for the enhancement of ion heating by anomalous resistivity. The anomalous resistivity associated with the turbulence is limited by electron heating, so that CDICI saturates at transverse temperature that is substantially higher than in the absence of resistivity. It is suggested that this process demonstrates a positive feedback loop in the interaction between CDICI, anomalous resistivity, and parallel large-scale dynamics in the topside ionosphere.

Ganguli, Supriya B.↗
DOE OSTI2015

Progress in preparing scenarios for operation of the International Thermonuclear Experimental Reactor

The development of operating scenarios is one of the key issues in the research for ITER which aims to achieve a fusion gain (Q) of ∼10, while producing 500 MW of fusion power for ≥300 s. The ITER Research plan proposes a success oriented schedule starting in hydrogen and helium, to be followed by a nuclear operation phase with a rapid development towards Q ∼ 10 in deuterium/tritium. The Integrated Operation Scenarios Topical Group of the International Tokamak Physics Activity initiates joint activities among worldwide institutions and experiments to prepare ITER operation. Plasma formation studies report robust plasma breakdown in devices with metal walls over a wide range of conditions, while other experiments use an inclined EC launch angle at plasma formation to mimic the conditions in ITER. Simulations of the plasma burn-through predict that at least 4 MW of Electron Cyclotron heating (EC) assist would be required in ITER. For H-modes at q{sub 95} ∼ 3, many experiments have demonstrated operation with scaled parameters for the ITER baseline scenario at n{sub e}/n{sub GW} ∼ 0.85. Most experiments, however, obtain stable discharges at H{sub 98(y,2)} ∼ 1.0 only for β{sub N} = 2.0–2.2. For the rampup in ITER, early X-point formation is recommended, allowing auxiliary heating to reduce the flux consumption. A range of plasma inductance (l{sub i}(3)) can be obtained from 0.65 to 1.0, with the lowest values obtained in H-mode operation. For the rampdown, the plasma should stay diverted maintaining H-mode together with a reduction of the elongation from 1.85 to 1.4. Simulations show that the proposed rampup and rampdown schemes developed since 2007 are compatible with the present ITER design for the poloidal field coils. At 13–15 MA and densities down to n{sub e}/n{sub GW} ∼ 0.5, long pulse operation (>1000 s) in ITER is possible at Q ∼ 5, useful to provide neutron fluence for Test Blanket Module assessments. ITER scenario preparation in hydrogen and helium requires high input power (>50 MW). H-mode operation in helium may be possible at input powers above 35 MW at a toroidal field of 2.65 T, for studying H-modes and ELM mitigation. In hydrogen, H-mode operation is expected to be marginal, even at 2.65 T with 60 MW of input power. Simulation code benchmark studies using hybrid and steady state scenario parameters have proved to be a very challenging and lengthy task of testing suites of codes, consisting of tens of sophisticated modules. Nevertheless, the general basis of the modelling appears sound, with substantial consistency among codes developed by different groups. For a hybrid scenario at 12 MA, the code simulations give a range for Q = 6.5–8.3, using 30 MW neutral beam injection and 20 MW ICRH. For non-inductive operation at 7–9 MA, the simulation results show more variation. At high edge pedestal pressure (T{sub ped} ∼ 7 keV), the codes predict Q = 3.3–3.8 using 33 MW NB, 20 MW EC, and 20 MW ion cyclotron to demonstrate the feasibility of steady-state operation with the day-1 heating systems in ITER. Simulations using a lower edge pedestal temperature (∼3 keV) but improved core confinement obtain Q = 5–6.5, when ECCD is concentrated at mid-radius and ∼20 MW off-axis current drive (ECCD or LHCD) is added. Several issues remain to be studied, including plasmas with dominant electron heating, mitigation of transient heat loads integrated in scenario demonstrations and (burn) control simulations in ITER scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Characterization of uncertainties in electron-argon collision cross sections

Abstract The predictive capability of a plasma discharge model depends on accurate representations of electron-impact collision cross sections, which determine the corresponding reaction rates and electron transport properties. The values of cross sections can be known only approximately either through experiments or simulations and are thus subject to uncertainties. Quantifying the uncertainties in plasma simulations allows us to assess the reliability of simulations and to provide a basis for interpreting discrepancies between simulations and experiments. For such uncertainty quantification of plasma simulations, it is essential to quantify the uncertainties of the underlying cross sections. Although much effort has been committed to calibrate the cross section values, their uncertainties are not well investigated. We characterize uncertainties in electron-argon atom collision cross sections using a Bayesian framework. Six collision processes—elastic momentum transfer, ionization, and four excitations—are characterized with semi-empirical models, which effectively capture the features important to the macroscopic properties of the plasma. A probability model for the uncertain parameters of these semi-empirical models is developed. Specifically, a Gaussian-process likelihood model is proposed to capture discrepancies among data sets, as well as the model-form inadequacies of the semi-empirical models. Two other likelihood models are compared with the proposed Gaussian-process model, to illustrate the importance of the choice of the likelihood model. The cross section models are calibrated using the electron-beam experiments and ab-inito quantum simulations. The resulting calibrated uncertainties capture well the scattering among the data sets. The calibrated cross section models are further validated against swarm-parameter experiments and zero-dimensional Boltzmann equation simulations of widely used cross section datasets.

Chung, Seung Whan (ORCID:0000000302501549)↗
NASA NTRS2010

Review of Leading Approaches for Mitigating Hypersonic Vehicle Communications Blackout and a Method of Ceramic Particulate Injection Via Cathode Spot Arcs for Blackout Mitigation

Vehicles flying at hypersonic velocities within the atmosphere become enveloped in a "plasma sheath" that prevents radio communication, telemetry, and most importantly, GPS signal reception for navigation. This radio "blackout" period has been a problem since the dawn of the manned space program and was an especially significant hindrance during the days of the Apollo missions. An appropriate mitigation method must allow for spacecraft to ground control and ground control to spacecraft communications through the reentry plasma sheath. Many mitigation techniques have been proposed, including but not limited to, aerodynamic shaping, magnetic windows, and liquid injection. The research performed on these mitigation techniques over the years will be reviewed and summarized, along with the advantages and obstacles that each technique will need to overcome to be practically implemented. A unique approach for mitigating the blackout communications problem is presented herein along with research results associated with this method. The novel method involves the injection of ceramic metal-oxide particulate into a simulated reentry plasma to quench the reentry plasma. Injection of the solid ceramic particulates is achieved by entrainment within induced, energetic cathode spot flows.

Gillman, Eric D.↗
DOE OSTI2017

Study of Conditions for Obtaining Quasi-Stationary Scenarios in T-15 MD Tokamak

The study of conditions for obtaining quasi-stationary scenarios in the T-15MD tokamak is performed. Results of simulation and optimization of T-15MD regimes with a fully noninductive current drive using NBICD systems with the power of 6–8 MW, ECCD system with the power ~8 MW, and RFCD system with power up to 7–8 MW are presented. It is shown that, at the somewhat reduced value as compared to the basic value of the toroidal magnetic field B{sub t} ~ 1.5 T in the T-15MD tokamak, discharges can be produced with a fully noninductive current of about 1 MA, plasma temperature of several keV, plasma density of about (3–7) × 10{sup 19} m{sup –3}, and discharge duration of about 20 s. Such discharges are of interest for further experimental study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Electromagnetics Simulations with Vector-valued Finite Elements in MOOSE

Vector-valued finite elements have become an essential part of modern electromagnetic (EM) numerical simulation. This talk will describe the development and initial validation of a MOOSE-based EM solver application, ELK (Electromagnetics Library for Kinetics and fluids), that utilizes the recently added vector-valued finite element system. The primary motivation for the development of ELK was to study EM wave propagation in radio-frequency (rf) driven systems, such as those encountered in the study of antenna structures, waveguides, and plasma physics. In simulating these scenarios with traditional scalar-valued finite elements, boundary conditions on tangential field components can sometimes be difficult to define, and field discontinuities encountered with certain geometric or structural features (e.g., abrupt corners and inhomogeneous materials) can often lead to non-physical results. Vector-valued elements can largely mitigate and even eliminate these difficulties. ELK verification and validation efforts thus far have focused on standard wave propagation benchmarks for antennas and waveguides, and the impact of vector-valued elements on the implementation and results of these simulations will be described. Progress using ELK for plasma physics simulation, where inhomogeneous simulation domains are commonplace, will also be discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Plasmoid instability in the semi-collisional regime

Here, we investigate analytically and numerically the semi-collisional regime of the plasmoid instability, defined by the inequality $δ_{SP} \gg ρ_s \gg δ_{in}$, where $δ_{SP}$ is the width of a Sweet–Parker current sheet, $ρ_s$ is the ion sound Larmor radius and $δ_{in}$ is the width of the boundary layer that arises in the plasmoid instability analysis. Theoretically, this regime is predicted to exist if the Lundquist number S and the length of the current sheet L are such that ($L/ρ_s$) 14/9 < S < ($L/ρ_s$) 2 (for a sinusoidal-like magnetic configuration; for a Harris-type sheet the lower bound is replaced with ($L/ρ_s$) 8/5 ). These bounds are validated numerically by means of simulations using a reduced gyrokinetic model conducted with the code Viriato. Importantly, this regime is conjectured to allow for plasmoid formation at relatively low, experimentally accessible, values of the Lundquist number. Our simulations obtain plasmoid instability at values of S as low as ~250. The simulations do not prescribe a Sweet–Parker sheet; rather, one is formed self-consistently during the nonlinear evolution of the initial tearing mode configuration. This proves that this regime of the plasmoid instability is realizable, at least at the relatively low values of the Lundquist number that are accessible to current dedicated experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

Assessment of Pole Erosion in a Magnetically Shielded Hall Thruster

Numerical simulations of a 6-kW laboratory Hall thruster called H6 have been performed to quantify the erosion rate at the inner pole. The assessments have been made in two versions of the thruster, namely the unshielded (H6US) and magnetically shielded (H6MS) configurations. The simulations have been performed with the 2-D axisymmetric code Hall2De which employs a new multi-fluid ion algorithm to capture the presence of low-energy ions in the vicinity of the poles. It is found that the maximum computed erosion rate at the inner pole of the H6MS exceeds the measured rate of back-sputtered deposits by ~4.5 times. This explains only part of the surface roughening that was observed after a 150-h wear test, which covered most of the pole area exposed to the plasma. For the majority of the pole surface the computed erosion rates are found to be below the back-sputter rate and comparable to those in the H6US which exhibited little to no sputtering in previous tests. Possible explanations for the discrepancy are discussed.

plasma simulations↗
DOE OSTI2024

Optimization of nonlinear turbulence in stellarators

We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native gyrokinetic code GX with the stellarator equilibrium and optimization code DESC . Since using GX allows for fast nonlinear simulations, we directly optimize for reduced nonlinear heat fluxes. To handle the noisy heat flux traces returned by these simulations, we employ the simultaneous perturbation stochastic approximation (SPSA) method that only uses two objective function evaluations for a simple estimate of the gradient. We show several examples that optimize for both reduced heat fluxes and good quasi-symmetry as a proxy for low neoclassical transport. Finally, we run full transport simulations using the T3D stellarator transport code to evaluate the changes in the macroscopic profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Power combination of a self-coherent high power microwave source

In our previous work, generating two phase-locked high power microwaves (HPMs) in a single self-coherent HPM device has been demonstrated. In this paper, after optimizing the structure of the previous self-coherent source, we design a power combiner with a folded phase-adjustment waveguide to realize power combination between its two sub-sources. Further particle-in-cell simulation of the combined source shows that when the diode voltage is 687 kV and the axial magnetic field is 0.8 T, a combined output microwave with 3.59 GW and 9.72 GHz is generated. The impedance of the combined device is 36 Ω and the total power conversion efficiency is 28%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The appearance, motion, and disappearance of three-dimensional magnetic null points

While theoretical models and simulations of magnetic reconnection often assume symmetry such that the magnetic null point when present is co-located with a flow stagnation point, the introduction of asymmetry typically leads to non-ideal flows across the null point. To understand this behavior, we present exact expressions for the motion of three-dimensional linear null points. The most general expression shows that linear null points move in the direction along which the magnetic field and its time derivative are antiparallel. Null point motion in resistive magnetohydrodynamics results from advection by the bulk plasma flow and resistive diffusion of the magnetic field, which allows non-ideal flows across topological boundaries. Null point motion is described intrinsically by parameters evaluated locally; however, global dynamics help set the local conditions at the null point. During a bifurcation of a degenerate null point into a null-null pair or the reverse, the instantaneous velocity of separation or convergence of the null-null pair will typically be infinite along the null space of the Jacobian matrix of the magnetic field, but with finite components in the directions orthogonal to the null space. Not all bifurcating null-null pairs are connected by a separator. Furthermore, except under special circumstances, there will not exist a straight line separator connecting a bifurcating null-null pair. The motion of separators cannot be described using solely local parameters because the identification of a particular field line as a separator may change as a result of non-ideal behavior elsewhere along the field line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

Preface: C/NOFS Results and Equatorial Ionospheric Dynamics

The Communication/Navigation Outage Forecasting System (C/NOFS) satellite was launched into orbit in April 2008 as part of an ongoing effort to understand and identify plasma irregularities that adversely impact the propagation of radio waves in the upper atmosphere. Combined with recent improvements in radar, airglow, and ground-based studies, as well as state-of-the-art modeling techniques, the C/NOFS mission has led to new insights into equatorial ionospheric electrodynamics. In order to document these advances, the C/NOFS Results and Equatorial Dynamics Technical Interchange Meeting was held in Albuquerque, New Mexico from 12 to 14 March 2013. The meeting was a great success with 55 talks and 22 posters, and covered topics including the numerical simulations of plasma irregularities, the effects of atmospheric tides, stratospheric phenomena, and magnetic storms on the upper atmosphere, causes and predictions of scintillation-causing ionospheric irregularities, current and future instrumentation efforts in the equatorial region. The talks were broken into the following three topical sessions: A. Ambient Ionosphere and Thermosphere B. Transient Phenomena in the Low-Latitude Ionosphere C. New Missions, New Sensors, New Science and Engineering Issues. The following special issue was planned as a follow-up to the meeting. We would like to thank Mike Pinnock, the editors and staff of Copernicus, and our reviewers for their work in bringing this special issue to the scientific community. Our thanks also go to Patricia Doherty and the meeting organizing committee for arranging the C/NOFS Technical Interchange Meeting.

CNOFS↗
NASA NTRS2014

Reconstruction of a Broadband Spectrum of Alfvenic Fluctuations

Alfvenic fluctuations in the solar wind exhibit a high degree of velocities and magnetic field correlations consistent with Alfven waves propagating away and toward the Sun. Two remarkable properties of these fluctuations are the tendencies to have either positive or negative magnetic helicity (-1 less than or equal to sigma(sub m) less than or equal to +1) associated with either left- or right- topological handedness of the fluctuations and to have a constant magnetic field magnitude. This paper provides, for the first time, a theoretical framework for reconstructing both the magnetic and velocity field fluctuations with a divergence-free magnetic field, with any specified power spectral index and normalized magnetic- and cross-helicity spectrum field fluctuations for any plasma species. The spectrum is constructed in the Fourier domain by imposing two conditions-a divergence-free magnetic field and the preservation of the sense of magnetic helicity in both spaces-as well as using Parseval's theorem for the conservation of energy between configuration and Fourier spaces. Applications to the one-dimensional spatial Alfvenic propagation are presented. The theoretical construction is in agreement with typical time series and power spectra properties observed in the solar wind. The theoretical ideas presented in this spectral reconstruction provide a foundation for more realistic simulations of plasma waves, solar wind turbulence, and the propagation of energetic particles in such fluctuating fields.

waves↗