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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

FreeDam – A webtool for free-electron laser-induced damage in femtosecond X-ray crystallography

Over the last decade X-ray free-electron laser (XFEL) sources have been made available to the scientific community. One of the most successful uses of these new machines has been protein crystallography. When samples are exposed to the intense short X-ray pulses provided by the XFELs, the sample quickly becomes highly ionized and the atomic structure is affected. Here we present a webtool dubbed FreeDam based on non-thermal plasma simulations, for estimation of radiation damage in free-electron laser experiments in terms of ionization, temperatures and atomic displacements. The aim is to make this tool easily accessible to scientists who are planning and performing experiments at XFELs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Diffraction of radio frequency waves by spatially modulated interfaces in the plasma edge in tokamaks

We report the use of radio frequency (RF) waves in fusion plasmas for heating, for non-inductive current generation, for profile control and for diagnostics has been well established. The RF waves, excited by antenna structures placed near the wall of a fusion device, have to propagate through density fluctuations at the plasma edge. These fluctuations can modify the properties of the RF waves that propagate towards the core of the plasma. A full-wave electromagnetic computational code ScaRF based on the finite difference frequency domain (FDFD) method has been developed to study the effect of density turbulence on RF waves. The anisotropic plasma permittivity used in the scattering studies is that for a magnetized, cold plasma. The code is used to study the propagation of an RF plane wave through a modulated, spatially periodic density interface. Such an interface could arise in the edge region due to magnetohydrodynamic instability or drift waves. The frequency of the plane wave is taken to be in the range of the electron cyclotron frequency. The scattering analysis is applicable to ITER-like plasmas, as well as to plasmas in medium sized tokamaks such as TCV, ASDEX-U and DIII-D. The effect of different density contrasts across the interface and of different spatial modulations are discussed. While ScaRF is used to study a periodic density fluctuation, the code is general enough to include different varieties of density fluctuations in the edge region – such as blobs and filaments, and spatially random fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A hybrid full-wave Markov chain approach to calculating radio-frequency wave scattering from scrape-off layer filaments

The interaction of radio-frequency (RF) waves with edge turbulence modifies the incident wave spectrum, and can significantly affect RF heating and current drive in tokamaks. Previous lower hybrid (LH) scattering models have either used the weak-turbulence approximation, or treated more realistic, filamentary turbulence in the ray tracing limit. In this study, a new model is introduced which retains full-wave effects of RF scattering in filamentary turbulence. First, a Mie-scattering technique models the interaction of an incident wave with a single Gaussian filament. Next, an effective differential scattering width is derived for a statistical ensemble of filaments. Lastly, a Markov chain solves for the transmitted wave spectrum in slab geometry. This model is applied to LH launching for current drive. The resulting wave spectrum is asymmetrically broadened in angular wavenumber space. This asymmetry is not accounted for in previous LH scattering models. The modified wave spectrum is coupled to a ray tracing/Fokker–Planck solver (GENRAY/CQL3D) to study its impact on current drive. The resulting current profile is greatly altered, and there is significant increase in the on-axis current and decrease in the off-axis peaks. This is attributed to a portion of the modified wave spectrum that is strongly dampened on-axis during the first pass.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Response to 'Comment on 'Undamped electrostatic plasma waves''[Phys. Plasmas 20, 034701 (2013)]

Numerical and experimental evidence is given for the occurrence of the plateau states and concomitant corner modes proposed in Valentini et al.[Phys. Plasmas 19, 092103 (2012)]. It is argued that these states provide a better description of reality for small amplitude off-dispersion disturbances than the conventional Bernstein-Greene-Kruskal or cnoidal states such as those proposed in Schamel [Phys. Plasmas 20, 034701 (2013)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Evanescent wave decomposition in a novel resonator comprising unmagnetized and magnetized plasma layers

A 4 Multiplication-Sign 4 transfer matrix method has been applied to study the decomposition of any elliptically polarized wave in a magnetized resonator. When the incident elliptically polarized wave passes through the structure, it is orthogonally decomposed into two circular polarizations at two resonance frequencies. Without changing the structure of the resonator, the positions of the resonant frequencies of the right- and left-handed circularly polarized waves can be modulated by changing the external magnetized field. The results show that the proposed magnetized structure can be used to design a novel resonator, which can be applied in the decomposition of polarized electromagnetic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

On the multistream approach of relativistic Weibel instability. II. Bernstein-Greene-Kruskal-type waves in magnetic trapping

The stationary state with magnetically trapped particles is investigated at the saturation of the relativistic Weibel instability, within the “multiring” model in a Hamiltonian framework. The multistream model and its multiring extension have been developed in Paper I, under the assumption that the generalized canonical momentum is conserved in the perpendicular direction. One dimensional relativistic Bernstein-Greene-Kruskal waves with deeply trapped particles are addressed using similar mathematical formalism developed by Lontano et al.[Phys. Plasmas 9, 2562 (2002); Phys. Plasmas 10, 639 (2003)] using several streams and in the presence of both electrostatic and magnetic trapping mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Nonlinear scattering term in the gyrokinetic Vlasov equation

Nonlinear scattering term is found from the nonlinear gyrokinetic equation by decoupling the perturbed gyrocenter motion from the unperturbed motion. The gyro-center distribution function is determined by the well-understood unperturbed motion, with the effects of fields perturbation included in the nonlinear scattering term, which explicitly reveals the nonlinear stochastic dissipation on the time scale longer than the wave correlation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

On applicability of the “thermalized potential” solver in simulations of the plasma flow in Hall thrusters

In Hall thrusters, the potential distribution plays an important role in discharge processes and ion acceleration. This paper presents a 2D potential solver in the Hall thruster instead of the “thermalized potential”, and compares equipotential contours solved by these two methods for different magnetic field conditions. The comparison results reveal that the expected “thermalized potential” works very well when the magnetic field is nearly uniform and electron temperature is constant along the magnetic field lines. However for the case with a highly non-uniform magnetic field or variable electron temperature along the magnetic field lines, the “thermalized potential” is not accurate. In some case with magnetic separatrix inside the thruster channel, the “thermalized potential” model cannot be applied at all. In those cases, a full 2D potential solver must be applied. Overall, this paper shows the limit of applicability of the “thermalized potential” model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Narrow gap electronegative capacitive discharges

Narrow gap electronegative (EN) capacitive discharges are widely used in industry and have unique features not found in conventional discharges. In this paper, plasma parameters are determined over a range of decreasing gap length L from values for which an electropositive (EP) edge exists (2-region case) to smaller L-values for which the EN region connects directly to the sheath (1-region case). Parametric studies are performed at applied voltage V{sub rf}=500 V for pressures of 10, 25, 50, and 100 mTorr, and additionally at 50 mTorr for 1000 and 2000 V. Numerical results are given for a parallel plate oxygen discharge using a planar 1D3v (1 spatial dimension, 3 velocity components) particle-in-cell (PIC) code. New interesting phenomena are found for the case in which an EP edge does not exist. This 1-region case has not previously been investigated in detail, either numerically or analytically. In particular, attachment in the sheaths is important, and the central electron density n{sub e0} is depressed below the density n{sub esh} at the sheath edge. The sheath oscillations also extend into the EN core, creating an edge region lying within the sheath and not characterized by the standard diffusion in an EN plasma. An analytical model is developed using minimal inputs from the PIC results, and compared to the PIC results for a base case at V{sub rf}=500 V and 50 mTorr, showing good agreement. Selected comparisons are made at the other voltages and pressures. A self-consistent model is also developed and compared to the PIC results, giving reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Calibration improvements expand filterscope diagnostic use

The filterscope diagnostic on DIII-D utilizes photomultiplier tubes to measure visible light emission from the plasma. The system has undergone a substantial upgrade since previous attempts to cross-calibrate the filterscope with other spectroscopic diagnostics were unsuccessful. The optics now utilize a dichroic mirror to initially split the light at nearly 99% transmission or reflectance for light below or above 550 nm. This allows the system to measure D α emission without degrading visible light emission from the plasma for wavelengths below 550 nm (to measure D β , D γ , W–I, C-III, etc.). Additional optimization of the optical components and calibration techniques reduce the error in the signal up to 10% in some channels compared to previous methods. Cross-calibration measurements with two other high resolution spectroscopic diagnostics now show excellent agreement for the first time. This expands the capabilities of the filterscope system allowing measurement of divertor detachment, emission profiles, edge-localized mode behavior, and plasma–wall interactions. Additionally, it enables direct comparisons against calculations from boundary plasma simulations. These were not possible before.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗
DOE OSTI2021

Asynchronous and Load-Balanced Union-Find for Distributed and Parallel Scientific Data Visualization and Analysis

We present a novel distributed union-find algorithm that features asynchronous parallelism and k-d tree based load balancing for scalable visualization and analysis of scientific data. Applications of union-find include level set extraction and critical point tracking, but distributed union-find can suffer from high synchronization costs and imbalanced workloads across parallel processes. In this study, we prove that global synchronizations in existing distributed union-find can be eliminated without changing final results, allowing overlapped communications and computations for scalable processing. We also use a k-d tree decomposition to redistribute inputs, in order to improve workload balancing. We benchmark the scalability of our algorithm with up to 1,024 processes using both synthetic and application data. Here, we demonstrate the use of our algorithm in critical point tracking and super-level set extraction with high-speed imaging experiments and fusion plasma simulations, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Highly selective Si 3 N 4 /SiO 2 etching using an NF 3 /N 2 /O 2 /H 2 remote plasma. II. Surface reaction mechanism

Developing processes for highly selective etching of silicon nitride (Si 3 N 4 ) with respect to silicon dioxide (SiO 2 ) is a major priority for semiconductor fabrication processing. In this paper and in Paper I [Volynets et al., J. Vac. Sci. Technol. A 38, 023007 (2020)], mechanisms are discussed for highly selective Si3N4 etching in a remote plasma based on experimental and theoretical investigations. The Si 3 N 4 /SiO 2 etch selectivity of up to 380 was experimentally produced using a remote plasma sustained in NF 3 /N 2 /O 2 /H 2 mixtures. A selectivity strongly depends on the flow rate of H 2 , an effect attributed to the formation of HF molecules in vibrationally excited states that accelerate etching reactions. Based on experimental measurements and zero-dimensional plasma simulations, an analytical etching model was developed for etch rates as a function of process parameters. Reaction rates and sticking coefficients were provided by quantum chemistry models and also fitted to the experimental results. Etch rates from the analytical model show good agreement with the experimental results and demonstrate why certain etchants accelerate or inhibit the etch process. In particular, the modeling shows the important role of HF molecules in the first vibrationally excited state [HF( v = 1)] in achieving high Si 3 N 4 /SiO 2 selectivity.

42 ENGINEERING↗
DOE OSTI2025

Minimizing CGYRO HPC Communication Costs in Ensembles with XGYRO by Sharing the Collisional Constant Tensor Structure

First-principles fusion plasma simulations are both compute and memory intensive, and CGYRO is no exception. The use of many HPC nodes to fit the problem in the available memory thus results in significant communication overhead, which is hard to avoid for any single simulation. That said, most fusion studies are composed of ensembles of simulations, so we developed a new tool, named XGYRO, that executes a whole ensemble of CGYRO simulations as a single HPC job. By treating the ensemble as a unit, XGYRO can alter the global buffer distribution logic and apply optimizations that are not feasible on any single simulation, but only on the ensemble as a whole. The main saving comes from the sharing of the collisional constant tensor structure, since its values are typically identical between parameter-sweep simulations. This data structure dominates the memory consumption of CGYRO simulations, so distributing it among the whole ensemble results in drastic memory savings for each simulation, which in turn results in overall lower communication overhead.

CGYRO↗