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,243 records · Page 69

An X-band high-impedance relativistic klystron amplifier with an annular explosive cathode

The feasibility of employing an annular beam instead of a solid one in the X-band high-impedance relativistic klystron amplifier (RKA) is investigated in theory and simulation. Small-signal theory analysis indicates that the optimum bunching distance, fundamental current modulation depth, beam-coupling coefficient, and beam-loaded quality factor of annular beams are all larger than the corresponding parameters of solid beams at the same beam voltage and current. An annular beam RKA and a solid beam RKA with almost the same geometric parameters are compared in particle-in-cell simulation. Output microwave power of 100 MW, gain of 50 dB, and power conversion efficiency of 42% are obtained in an annular beam RKA. The annular beam needs a 15% lower uniform guiding magnetic field than the solid beam. Our investigations demonstrate that we are able to use a simple annular explosive cathode immersed in a lower uniform magnetic field instead of a solid thermionic cathode in a complicated partially shielding magnetic field for designing high-impedance RKA, which avoids high temperature requirement, complicated electron-optical system, large area convergence, high current density, and emission uniformity for the solid beam. An equivalent method for the annular beam and the solid beam on bunching features is proposed and agrees with the simulation. The annular beam has the primary advantages over the solid beam that it can employ the immersing uniform magnetic field avoiding the complicated shielding magnetic field system and needs a lower optimum guiding field due to the smaller space charge effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimation of China Fusion Engineering Test Reactor performance and burning fraction in different pellet fueling scenarios by a multi-species radial transport model

Tritium self-sufficiency in future deuterium–tritium fusion reactors is a crucial challenge. As an engineering test reactor, the China Fusion Engineering Test Reactor requires a burning fraction of 3% for the goal to test the accessibility to the future fusion plant. To self-consistently simulate burning plasmas with profile changes in pellet injection scenarios and to estimate the corresponding burning fraction, a one-dimensional multi-species radial transport model is developed in the BOUT++ framework. Here several pellet-fueling scenarios are then tested in the model. The results show that the increased fueling depth improves the burning fraction by particle confinement improvement and fusion power increase. Nevertheless, by increasing the depth, the pellet cooling-down may significantly lower the temperature in the core region. Taking the density perturbation into consideration, the reasonable parameters of the fueling scenario in these simulations are estimated as pellet radius ${r}_{{\rm{p}}}=3\,{\rm{mm}},$ injection rate $=\,4\,\mathrm{Hz},$ and pellet injection velocity $=\,1000\mbox{-}2000\,{\rm{m}}\,{{\rm{s}}}^{-1}$ without drift or $450\,{\rm{m}}\,{{\rm{s}}}^{-1}$ with high-field-side drift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-consistent evolution of plasma discharge and electromagnetic fields in a microwave pulse compressor

Nanosecond-scale evolution of plasma and RF electromagnetic fields during the release of energy from a microwave pulse compressor with a plasma interference switch was investigated numerically using the code MAGIC. The plasma was simulated in the scope of the gas conductivity model in MAGIC. The compressor embodied an S-band cavity and H-plane waveguide tee with a shorted side arm filled with pressurized gas. In a simplified approach, the gas discharge was initiated by setting an external ionization rate in a layer crossing the side arm waveguide in the location of the electric field antinode. It was found that with increasing ionization rate, the microwave energy absorbed by the plasma in the first few nanoseconds increases, but the absorption for the whole duration of energy release, on the contrary, decreases. In a hybrid approach modeling laser ignition of the discharge, seed electrons were set around the electric field antinode. In this case, the plasma extends along the field forming a filament and the plasma density increases up to the level at which the electric field within the plasma decreases due to the skin effect. Then, the avalanche rate decreases but the density still rises until the microwave energy release begins and the electric field becomes insufficient to support the avalanche process. The extraction of the microwave pulse limits its own power by terminating the rise of the plasma density and filament length. For efficient extraction, a sufficiently long filament of dense plasma must have sufficient time to be formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of radiation flux in certain band via the preheat of aluminum sample

Quantitative evaluation of the fractions of high energy x-rays in a hohlraum is crucial to the indirect driven-drive scheme of inertial confinement fusion and many other applications in high energy density physics. Preheat of a sample due to x-rays sensitively depends on optical thin photons. Analyzing the motion of a sample due to preheat can thus provide valuable information of those x-rays. In this article, we propose a method to infer the temporal evolution of the x-ray fluxes in the bands of our interest. By matching the simulation results to the motions of an aluminum sample, we can infer the time-resolved x-ray fluxes around the aluminum K-edge and the gold M-band inside the hohlraum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Verification of electromagnetic fluid-kinetic hybrid electron model in global gyrokinetic particle simulation

The fluid-kinetic hybrid electron model is verified in global gyrokinetic particle simulation of linear electromagnetic drift-Alfvenic instabilities in tokamak. In particular, we have recovered the {beta}-stabilization of the ion temperature gradient mode, transition to collisionless trapped electron mode, and the onset of kinetic ballooning mode as {beta}{sub e} (ratio of electron kinetic pressure to magnetic pressure) increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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↗

An improved model to determine the inception of positive upward leader-streamer system considering the leader propagation during dark period

Stem-leader transition and front-streamer inception are two essential conditions for the inception of positive upward leader-streamer system (LSS). Previous models have not considered the initial-leader propagation during dark period and have not been verified systematically. In this paper, a series of positive upward discharge simulation experiments was designed and carried out. Characteristic parameters of the discharge process related to the inception of positive upward LSS, namely, the first-corona inception voltage, the first-corona charge, the dark period, and the LSS inception voltage, were obtained. By comparing these experiment results with simulation results calculated using previous models, it was found that it is improper to assume that the length of the initial leader is a fixed value. Finally, an improved inception model of positive upward LSS considering the leader propagation during dark period was developed and verified with experiment results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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.↗

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↗

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)↗

Initial conditions of radiative shock experiments

We performed experiments at the Omega Laser Facility to characterize the initial, laser-driven state of a radiative shock experiment. These experiments aimed to measure the shock breakout time from a thin, laser-irradiated Be disk. The data are then used to inform a range of valid model parameters, such as electron flux limiter and polytropic γ, used when simulating radiative shock experiments using radiation hydrodynamics codes. The characterization experiment and the radiative shock experiment use a laser irradiance of ∼7 × 10{sup 14} W cm{sup −2} to launch a shock in the Be disk. A velocity interferometer and a streaked optical pyrometer were used to infer the amount of time for the shock to move through the Be disk. The experimental results were compared with simulation results from the Hyades code, which can be used to model the initial conditions of a radiative shock system using the CRASH code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics simulation of corona discharge induced ionic wind

Ionic wind devices or electrostatic fluid accelerators are becoming of increasing interest as tools for thermal management, in particular for semiconductor devices. In this work, we present a numerical model for predicting the performance of such devices; its main benefit is the ability to accurately predict the amount of charge injected from the corona electrode. Our multiphysics numerical model consists of a highly nonlinear, strongly coupled set of partial differential equations including the Navier-Stokes equations for fluid flow, Poisson's equation for electrostatic potential, charge continuity, and heat transfer equations. To solve this system we employ a staggered solution algorithm that generalizes Gummel's algorithm for charge transport in semiconductors. Predictions of our simulations are verified and validated by comparison with experimental measurements of integral physical quantities, which are shown to closely match.

36 MATERIALS SCIENCE↗

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.↗