Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal plasma”

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 109 records · Page 6

Energetic Particles - microturbulence interaction thrust (EP SIWG white paper)

This white paper outlines a strategic approach to Energetic Particle (EP) research, aimed at applying the results to the future burning devices such as ITER. Our goal is to deepen the understanding of the interplay between plasma background microturbulence and diverse species of energetic ions in burning plasma devices. This strategy is poised to advance our knowledge in this crucial area of fusion science, paving the way for significant breakthroughs in plasma physics and fusion energy research. We consider the confinement of energetic ions such as auxiliary heating beams, minority ions from ICRH, and fusion product alpha particles. Recent studies have underscored the potential of EP/microturbulence interaction to enhance fusion plasma performance by stabilizing microturbulence. However, it also poses risks, such as exacerbating Alfvénic eigenmode (AE) instabilities, which could compromise the sustainability of plasma discharge by causing EP losses. The interplay between EPs and background microturbulence is evident in the phenomenon of effective pitch angle scattering, a crucial aspect of the quasilinear (QL) theory. This theory is integral to developing numerically-efficient yet comprehensive and self-consistent approaches, recently employed to investigate the relaxation of energetic particle populations in the holistic modeling of fusion-grade plasmas. The primary objectives of this white paper thrust encompass a dual focus. First is that we need to delve into the mechanisms through which microturbulence engenders effective pitch angle scattering. This investigation will entail the examination of the presence of known or self-consistently predicted spectra of modes accountable for microturbulence. In these explorations, electrostatic microturbulence serves as a logical initial stride towards achieving the trust's objectives. This endeavor is anticipated to yield formulations expressing the parametric dependencies of the effective pitch angle scattering frequency on variables such as thermal plasma electron and/or ion temperatures, as well as their respective thermal conductivities. In our second objective, we aim to delve into the intricate formation of zonal flow (ZF) structures amidst the complex interplay of microturbulence and Alfvenic eigenmodes (AEs). This endeavor poses greater challenges as we seek to unravel the macroscopic manifestations influenced by microturbulence, termed as zonal structures (ZS), stemming from microturbulence-induced ZF. Drawing from this understanding, we anticipate employing a QL approach to yield comprehensive insights into the distribution function of energetic particles (EP) within phase space. This method entails resolving the dominant multidimensional phase space diffusion processes while effectively averaging over the rapid ballistic responses. Nevertheless, substantial strides remain imperative to realize a comprehensive whole-device modeling framework. This entails meticulous verification and validation exercises against experimental observations, as well as rigorous benchmarking against theoretical frameworks and numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of physics results from MAST upgrade towards core-pedestal-exhaust integration

Recent results from MAST Upgrade are presented, emphasising understanding the capabilities of this new device and deepening understanding of key physics issues for the operation of ITER and the design of future fusion power plants. The impact of MHD instabilities on fast ion confinement have been studied, including the first observation of fast ion losses correlated with Compressional and Global Alfvén Eigenmodes. High-performance plasma scenarios have been developed by tailoring the early plasma current ramp phase to avoid internal reconnection events, resulting in a more monotonic q profile with low central shear. The impact of m/n = 3/2, 2/1 and 1/1 modes on thermal plasma confinement and rotation profiles has been quantified, and scenarios optimised to avoid them have transiently reached values of normalised beta approaching 4.2. In pedestal and ELM physics, a maximum pedestal top temperature of ~350 eV has been achieved, exceeding the value achieved on MAST at similar heating power. Mitigation of type-I ELMs with n = 1 RMPs has been observed. Studies of plasma exhaust have concentrated on comparing conventional and Super-X divertor configurations, while X-point target, X-divertor and snowflake configurations have been developed and studied in parallel. In L-mode discharges, the separatrix density required to detach the outer divertors is approximately a factor 2 lower in the Super-X than the conventional configuration, in agreement with simulations. Detailed analysis of spectroscopy data from studies of the Super-X configuration reveal the importance of including plasma-molecule interactions and D 2 Fulcher band emission to properly quantify the rates of ionisation, plasma-molecule interactions and volumetric recombination processes governing divertor detachment. In H-mode with conventional and Super-X configurations, the outer divertors are attached in the former and detached in the latter with no impact on core or pedestal confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bespoke, High-Purity Coal Tar with Microwave Plasma Pyrolysis

H Quest Vanguard, Inc is developing microwave plasma driven ultra-rapid (100-500 ms residence times) coal pyrolysis. In this innovative process, entrained in an oxygen-free gas stream, the coal particles enter a non-thermal plasma and are rapidly heated (~10 5 K/sec), resulting in a carbonized exterior and electrically conductive coal particles. Microwave radiation can then directly couple to entrained particles creating localized high energy regions. Aromatic and polyaromatic tar components are stripped from the macromolecular structure of coal and ejected into the gas stream leaving behind porous char particles

Viswanathan, Vignesh↗

Generalized gluon distribution for quarkonium dynamics in strongly coupled N = 4 Yang-Mills theory

We study the generalized gluon distribution that governs the dynamics of quarkonium inside a non-Abelian thermal plasma characterizing its dissociation and recombination rates. This gluon distribution can be written in terms of a correlation function of two chromoelectric fields connected by an adjoint Wilson line. We formulate and calculate this object in N = 4 supersymmetric Yang-Mills theory at strong coupling using the AdS/CFT correspondence, allowing for a nonzero center-of-mass velocity v of the heavy quark pair relative to the medium. The effect of a moving medium on the dynamics of the heavy quark pair is described by the simple substitution T → γ T in agreement with previous calculations of other observables at strong coupling, where T is the temperature of the plasma in its rest frame, and γ = ( 1 − v 2 ) − 1 / 2 is the Lorentz boost factor. Such a velocity dependence can be important when the quarkonium momentum is larger than its mass. Contrary to general expectations for open quantum systems weakly coupled with large thermal environments, the contributions to the transition rates that are usually thought of as the leading ones in Markovian descriptions vanish in this strongly coupled plasma. This calls for new theoretical developments to assess the effects of strongly coupled non-Abelian plasmas on in-medium quarkonium dynamics. Finally, we compare our results with those from weakly coupled QCD, and find that the QCD result moves toward the N = 4 strongly coupled result as the coupling constant is increased within the regime of applicability of perturbation theory. This behavior makes it even more pressing to develop a non-Markovian description of quarkonium in-medium dynamics. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collaborative Research Proposal: Time Resolved Optical Emissions Spectroscopy and Laser Induced Fluorescence Spectroscopy of Nanosecond Pulsed Discharges in a Gas-Liquid Water Film Reactor

Electrical discharge plasma formed in contact with liquid water is of interest for a wide range of applications in chemical, biomedical, agricultural, electrical, and materials science and engineering. Such plasma reactors are of very timely importance since they also have significant disinfection capability by inactivating bacteria, viruses and other pathogens. Many types of plasma sources including those driven by AC, DC, RF, microwave, and pulsed electrical power supplies coupled to a wide range of different electrode configurations and reactor designs have been developed and explored which contact the plasma with liquid water. Recent roadmaps have recommended that further work is needed to develop our understanding of the fundamental chemical and physical process which occur at the interface of non-thermal plasma with liquid water solutions in order to advance this large diversity of applications which ultimately depend upon efficient production of key reactive chemical species. There is a wide range of interacting factors that affect the chemical reactions that occur in the plasma, in the liquid phase, and at the interface. These factors ultimately govern the key reactive chemical species formed and used in the various applications and they include a) the reactor design and input parameters, b) the discharge and transport processes, c) the plasma properties, and d) the resulting chemical reactions. For pulsed discharges, the power supply design and output parameters control the applied voltage, frequency, rise time, and width (duration) of the applied pulses. The reactor design involves specification of the gas-liquid contacting methods, electrode gap distance, reactor volume and shape, and gas and liquid flow rates achievable. The gas and liquid compositions as well as the liquid properties such as pH and conductivity are also of key importance in determination of the resulting chemical reactions. In addition, the important plasma properties include plasma gas temperature, electron density, electron energy (distribution), and size of the plasma channels which are all affected by the discharge and transport properties. Many studies have focused on specific aspects of these various processes. Of particular important and relevance to the proposed work is the utilization of nanosecond pulses to generate plasma in gas-liquid systems, liquid bubbles, underwater, and in gases. Recent advances in nanosecond pulses provide significant advantages in utilization with liquid water. For example, fast rise time and short pulses are less sensitive to water conductivity, such short pulses may provide advantages in fast temporal quenching of the plasma, and fundamental analysis of pulse properties, including pulse shape and width, may be facilitated by investigation of single filamentary fast pulses. Many studies have also dealt with the role of the gas composition on formation of reactive oxygen species (ROS) (i.e., hydroxyl radicals – ·OH, hydrogen peroxide – H 2 O 2 , various atomic oxygen species, ozone-O 3 , hydroperoxyl radicals HO 2 ·) and reactive nitrogen species (RNS) (i.e., nitrogen oxides – NO, NO 2 , N 2 O - collectively termed NO X , nitrite-NO 2 - , nitrate-NO 3 - , peroxynitrite-ONOO - ). The hydroxyl radical is the critical species in many chemical oxidation reactions for chemical degradation of toxic compounds in water and gases and for synthesis of some compounds. The mixture of various nitrogen oxide species is important for many biochemical and biological processes involved in biomedical and agricultural applications including disinfection and fertilizer production. The present proposal focuses on determination of the effects of time resolved electron density and hydroxyl radicals on plasma chemical reactions through collaboration with the Princeton Collaborative Low Temperature Plasma Research Facility (PCRF) at the Princeton Plasma Physics Laboratory (PPPL). In order to further investigate the role of the plasma generated electrons and hydroxyl radicals on the overall formation of the key species including hydrogen peroxide, hydroxyl radicals, and nitrogen oxides, the proposed work, thus seeks to determine high resolution time resolved electron density by optical emissions spectroscopy and time resolved hydroxyl radicals using laser induced fluorescence in the nanosecond discharge reactor. The combination of data on electron density and hydroxyl radical concentration will be utilized in the present work to more fully characterize the chemical reaction processes in this system and to advance the design, development, and operation of such chemical reactors for a wide range of applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Oxygenate Production from Plasma-Activated Reaction of CO 2 and Ethane

Upgrading ethane with CO 2 as a soft oxidant represents a desirable means of obtaining oxygenated hydrocarbons. This reaction is not thermodynamically feasible under mild conditions and has not been previously achieved as a one-step process. Non-thermal plasma was implemented as an alternative means of supplying energy to overcome activation barriers, leading to the production of alcohols, aldehydes, and acids as well as C 1 -C 5+ hydrocarbons under ambient pressure, with a maximum total oxygenate selectivity of 12%. A plasma chemical kinetic computational model was developed and found to be in good agreement with the experimental trends. Results from this study illustrate the potential to use plasma for the direct synthesis of value-added alcohols, acids, and aldehydes from ethane and CO 2 under mild conditions.

54 ENVIRONMENTAL SCIENCES↗

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗

From thermal catalysis to plasma catalysis: a review of surface processes and their characterizations

The use of atmospheric pressure plasma to enhance catalytic chemical reactions involves complex surface processes induced by the interactions of plasma-generated fluxes with catalyst surfaces. Industrial implementation of plasma catalysis necessitates optimizing the design and realization of plasma catalytic reactors that enable chemical reactions that are superior to conventional thermal catalysis approaches. This requires the fundamental understanding of essential plasma-surface interaction mechanisms of plasma catalysis from the aspect of experimental investigation and theoretical analysis or computational modeling. In addition, experimental results are essential to validate the relative theoretical models and hypotheses of plasma catalysis that was rarely understood so far, compared to conventional thermal catalysis. This overview focuses on two important application areas, nitrogen fixation and methane reforming, and presents a comparison of important aspects of the state of knowledge of these applications when performed using either plasma-catalysis or conventional thermal catalysis. We discuss the potential advantage of plasma catalysis over thermal catalysis from the aspects of plasma induced synergistic effect and in situ catalyst regeneration. In-situ/operando surface characterization of catalysts in plasma catalytic reactors is a significant challenge since the high pressure of realistic plasma catalysis systems preclude the application of many standard surface characterization techniques that operate in a low-pressure environment. Here, we present a review of the status of experimental approaches to probe gas-surface interaction mechanisms of plasma catalysis, including an appraisal of demonstrated approaches for integrating surface diagnostic tools into plasma catalytic reactors.

54 ENVIRONMENTAL SCIENCES↗

X-Ray Spectroscopic Studies of X-Pinch Plasmas with 3-5 Picosecond Resolution: A Quest for Clear Experimental Evidence for Radiative Collapse in the X-ray Spectra (Final Report)

Dense Z-pinch plasmas produced from current-carrying exploding wires frequently produce very brief (less than 0.1 nanosecond), intense, bursts of soft X-rays from tiny (about 1 micrometer), very hot (10-30 million K) regions of the plasma. This project addresses the physical processes that lead to those tiny, high-energy density plasmas and develops means to measure the conditions just before and during the X-ray bursts. One possible relevant process is radiative collapse, which starts if the power radiated by a dense Z-pinch, proportional to the density-squared, exceeds the resistive heating rate from the current flow, thereby cooling the plasma and enabling the magnetic pinching force to exceed the outward plasma thermal pressure. As the resulting radius reduction further increases the radiation rate, this process produces a runaway magnetic implosion. A competing mechanism is a sausage-like instability that can develop in a Z-pinch plasma at a lower current than is needed to induce radiative collapse. This mechanism can also produce tiny hot spots in the neck regions of the sausage-shaped plasma column. It is also possible that these two mechanisms work together to produce the X-ray emitting hot spots. If radiative collapse does contribute to hot spot formation, it would be terminated when the plasma becomes so dense that the radiation cannot escape, enabling the plasma kinetic pressure to build up and stop the current-driven implosion. This project aims to study the development of the hot spots by means of high temporal and spatial resolution X-ray spectroscopy to determine if there is evidence for the presence of the radiative collapse mechanism. Progress toward this goal is presented. To summarize, during the four-year period 09/15/2017 - 9/14/2021, we have collected necessary spectroscopic data to determine the plasma conditions and continue to do so. We are also now in the process of calibrating diagnostic systems (especially the x-ray streak camera system) and analyzing data to determine plasma conditions as a function of time before, during and after the X-ray burst.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards the first plasma-electron screening experiment

The enhancement of fusion reaction rates in a thermonuclear plasma by electron screening of the Coulomb barrier is an important plasma-nuclear effect that is present in stellar models but has not been experimentally observed. Experiments using inertial confinement fusion (ICF) implosions may provide a unique opportunity to observe this important plasma-nuclear effect. Herein, we show that experiments at the National Ignition Facility (NIF) have reached the relevant physical regime, with respect to the density and temperature conditions, but the estimated impacts of plasma screening on nuclear reaction rates are currently too small and need to be increased to lower the expected measurement uncertainty. Detailed radiation hydrodynamics simulations show that practical target changes, like adding readily available high-Z gases, and significantly slowing the inflight implosion velocity, while maintaining inflight kinetic energy, might be able to push these conditions to those where plasma screening effects may be measurable. We also perform synthetic data exercises to help understand where the anticipated experimental uncertainties will become important. But challenges remain, such as the detectability of the reaction products, non-thermal plasma effects, species separation, and impacts of spatial and temporal gradients. This work lays the foundation for future efforts to develop an important platform capable of the first plasma electron screening observation.

inertial confinement fusion↗

Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model

The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energetic particle physics: Chapter 7 of the special issue: on the path to tokamak burning plasma operation

We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A drift kinetic model for the expander region of a magnetic mirror

Here, we present a drift kinetic model for the free expansion of a thermal plasma out of a magnetic nozzle. This problem relates to plasma space propulsion systems, natural environments such as the solar wind, and end losses from the expander region of mirror magnetically confined fusion concepts such as the gas dynamic trap. The model incorporates trapped and passing orbit types encountered in the mirror expander geometry and maps to an upstream thermal distribution. This boundary condition and quasineutrality require the generation of an ambipolar potential drop of ~5T e /e, forming a thermal barrier for the electrons. The model for the electron and ion velocity distributions and fluid moments is confirmed with data from a fully kinetic simulation. Finally, the model is extended to account for a population of fast sloshing ions arising from neutral beam heating within a magnetic mirror, again resulting in good agreement with a corresponding kinetic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spheroidization of Nd–Fe–B particles

Spherical powders are required for many advanced manufacturing techniques due to their inherent requirement of flowability, either within feed tubes or during powder spreading. As advanced manufacturing of magnets continues to develop, new production methods for feedstocks are also sought. Plasma spheroidization is a high-yield method to produce spherical Nd–Fe–B powders from irregularly shaped particles, with advantages including high throughput and a well-controlled size distribution. Highly spherical Nd–Fe–B powders with large scale production (i.e., kg) have been demonstrated using an inductively coupled thermal plasma system; however, the magnetic properties of the output powder display significant degradation. The coercivity was decreased from the initial 8 kOe (636 kA/m) of the as-received to 0.7 kOe (55 kA/m) for spheroidized powders. Microstructural investigation reveals 6% Nd depletion caused by the extreme temperatures of the plasma, leading to the formation of low-coercivity α-Fe and a subsequent decrease in energy product. In conclusion, post-spheroidization heat treatments with Nd can partially mitigate the coercivity degradation, increasing to 1.7 kOe (135 kA/m), potentially offering a pathway toward spherical powders for a range of applications.

36 MATERIALS SCIENCE↗

Atmospheric Pressure Plasma Deposition of TiO2: A Review

Atmospheric pressure plasma (APP) deposition techniques are useful today because of their simplicity and their time and cost savings, particularly for growth of oxide films. Among the oxide materials, titanium dioxide (TiO2) has a wide range of applications in electronics, solar cells, and photocatalysis, which has made it an extremely popular research topic for decades. Here, we provide an overview of non-thermal APP deposition techniques for TiO2 thin film, some historical background, and some very recent findings and developments. First, we define non-thermal plasma, and then we describe the advantages of APP deposition. In addition, we explain the importance of TiO2 and then describe briefly the three deposition techniques used to date. We also compare the structural, electronic, and optical properties of TiO2 films deposited by different APP methods. Lastly, we examine the status of current research related to the effects of such deposition parameters as plasma power, feed gas, bias voltage, gas flow rate, and substrate temperature on the deposition rate, crystal phase, and other film properties. The examples given cover the most common APP deposition techniques for TiO2 growth to understand their advantages for specific applications. In addition, we discuss the important challenges that APP deposition is facing in this rapidly growing field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validation of energetic particle transport models for time-dependent integrated simulations of burning plasmas (Final ST Report)

The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Acceleration of the particle-in-cell code Osiris with graphics processing units

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗