Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic simulations”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The role of electron-scale turbulence in the JET tokamak: experiments and modelling

Dedicated electron heat transport experiments have been carried out in L- and H-mode Deuterium plasmas of the JET-ILW tokamak to identify the amount of electron heat carried by electron-scale electron temperature gradient (ETG) modes. Ion cyclotron resonance heating at different positions has been used to probe the response of the electron temperature inverse gradient length R/L Te to changes in electron heat flux qe, while different amounts of neutral beam heating allowed to scan the ratio of ion to electron temperature T e /T i , which is a key parameter for the onset of ETGs. Results indicate a steepening of the normalized qe vs R/L Te curve above R/LT e ~ 8 for T e /T i ≤ 1, suggestive of the ETG onset. Ion-scale gyro-kinetic (GK) simulations match the ion heat flux and the low-R/L Te part of the qe curve, but do not reproduce such steepening at high R/L Te . Multi-scale GK simulations covering both ion and electron scales and including one impurity bundling light and heavy species indicate an ETG contribution only for R/L Te values larger than the experimental ones. Sensitivity studies of such result are difficult to achieve due to limitation in numerical resources. The quasi-linear TGLF model has been used for sensitivity studies. With the same bundled impurity as the GK multi-scale, TGLF shows the qe steepening at much larger R/L Te values than in experiment, but when using the real mix of light impurities neglecting the heavy impurities, TGLF gets closer to the experimental results. Profile simulations with TGLF including both light and heavy impurities show over-prediction of T e profiles and in some cases also of density, but good T i predictions, confirming issues with the model electron stiffness for these plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A novel measurement of marginal Alfvén eigenmode stability during high power auxiliary heating in JET

The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Saturation of fishbone instability through zonal flows driven by energetic particle transport in tokamak plasmas

Gyrokinetic and kinetic-MHD simulations are performed for the fishbone instability in the DIII-D discharge #178631, chosen for validation of first-principles simulations to predict the energetic particle (EP) transport in an ITER prefusion baseline scenario. Fishbone modes are found to generate zonal flows, which dominate the fishbone saturation. The underlying mechanisms of the two-way fishbone-zonal flows nonlinear interplay are discussed in details. Numerical and analytical analyses identify the fishbone-induced EP redistribution as the dominant generation mechanism for zonal flows. The zonal flows modify the nonlinear dynamics of phase space zonal structures, which reduces the amount of EPs able to resonate with the mode, leading to an early fishbone saturation. Simulation results including zonal flows agree quantitatively with DIII-D experimental measurements of the fishbone saturation amplitude and EP transport, supporting this novel saturation mechanism by self- generated zonal flows. Moreover, the wave-particle mode-locking mechanism is shown to determine quantitatively the fishbone frequency down-chirping, as evident in GTC simulation results in agreement with predictions from analytical theory. Finally, the&#xD;fishbone-induced zonal flows are possibly responsible for the formation of an ion-ITB in the DIII-D discharge. Based on the low EP transport and the large zonal flow shearing rates associated with the fishbone instability in gyrokinetic simulations of the ITER scenario, it is conjectured that high performance scenarios could be designed in ITER burning plasmas through fishbone-induced ITBs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

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↗
DOE OSTI2021

Visualization of Fast Ion Phase-Space Flow Driven by Alfvén Instabilities

Fast ion phase-space flow, driven by Alfven eigenmodes (AEs), is measured by an imaging neutral particle analyzer in the DIII-D tokamak. The flow firstly appears near the minimum safety factor at the injection energy of neutral beams, and then moves radially inward and outward by gaining and losing energy, respectively. The flow trajectories in phase space align well with the intersection lines of the constant magnetic moment surfaces and constant E – (ω/n)P ζ surfaces, where E, P ζ are energy and toroidal canonical momentum of ions; ω and n are angular frequencies and toroidal mode numbers of AEs. It is found that the flow is so destructive that the thermalization of fast ions is no longer observed in regions of strong interaction. Here, the measured phase-space flow is consistent with nonlinear hybrid kinetic-magnetohydrodynamics simulation. Calculations of the relatively narrow phase-space islands reveal that fast ions must transition between different flow trajectories to experience large-scale phase-space transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Active Control of Alfvén Eigenmodes by Externally Applied 3D Magnetic Perturbations

The suppression and excitation of Alfvén eigenmodes have been experimentally obtained, for the first time, by means of externally applied 3D perturbative fields with different spatial spectra in a tokamak plasma. The applied perturbation causes an internal fast-ion redistribution that modifies the phase-space gradients responsible for driving the modes, determining, ultimately their existence. Hybrid kinetic-magnetohydrodynamic simulations reveal an edge resonant transport layer activated by the 3D perturbative field as the responsible mechanism for the fast-ion redistribution. Finally, the results presented here may help to control fast-ion driven Alfvénic instabilities in future burning plasmas with a significant fusion born alpha particle population.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Study on Particle Trajectory Error in Finite-Element Particle-in-Cell Algorithms

Particle-in-cell (PIC) algorithms are widely used for the simulation of kinetic plasmas. PIC algorithms account for the interaction between charged particles in a plasma and the electromagnetic field in ambient space, including self-field interactions. The objective of this article is to study the error in charged particle trajectories present in finite-element (FE)-based PIC algorithms on unstructured meshes. We study how the trajectory error behaves according to the FE mesh resolution and the matrix solver employed in the FE algorithm. The study is performed by considering a trajectory established by a parabolic electric potential and an axial magnetic force acting on the charged particle. Under a proper choice combination of electric and magnetic field parameters, the 2-D particle trajectories comprise closed orbits. Numerical errors cause small perturbations on the trajectories, with cumulative effects. As a result, the resulting orbital trajectories exhibit distortions including spurious apsidal precession. These distortions provide a clear imprint of the numerical error. Here, we also study the numerical error in a quantitative fashion by computing the distance norm between the trajectories obtained by the exact fields and numerical fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Numerical Analysis of Fuel Effects on Advanced Compression Ignition Using a Cooperative Fuel Research Engine Computational Fluid Dynamics Model

Growing environmental concerns and demand for a better fuel economy are driving forces that motivate the research for more advanced engines. Multi-mode combustion strategies have gained attention for their potential to provide high thermal efficiency and low emissions for light-duty applications. These strategies target optimizing the engine performance by correlating different combustion modes to load operating conditions. The extension from boosted spark ignition (SI) mode at high loads to advanced compression ignition (ACI) mode at low loads can be achieved by increasing the compression ratio and utilizing intake air heating. Further, in order to enable an accurate control of intake charge condition for ACI mode and rapid mode-switches, it is essential to gain fundamental insights into the autoignition process. Within the scope of ACI, homogeneous charge compression ignition (HCCI) mode is of significant interest. It is known for its potential benefits, operation at low fuel consumption, low NOx, and particulate matter (PM) emissions. In the present work, a virtual Cooperative Fuel Research (CFR) engine model is used to analyze fuel effects on ACI combustion. In particular, the effect of fuel octane sensitivity (S) (at constant Research Octane Number (RON)) on autoignition propensity is assessed under beyond-RON (BRON) and beyond-MON (BMON) ACI conditions. The three-dimensional CFR engine computational fluid dynamics (CFD) model employs a finite-rate chemistry approach with a multi-zone binning strategy to capture autoignition. Two binary blends with Research Octane Number (RON) of 90 are chosen for this study: primary reference fuel (PRF) with S = 0 and toluene–heptane (TH) blend with S = 10.8, representing paraffinic and aromatic gasoline surrogates. Reduced mechanisms for these blends are generated from a detailed gasoline surrogate kinetic mechanism. Simulation results with the reduced mechanisms are validated against experimental data from an in-house CFR engine, with respect to in-cylinder pressure, heat release rate, and combustion phasing. Thereafter, the sensitivity of combustion behavior to ACI operating condition (BRON versus BMON), air-fuel ratio (λ = 2 and 3), and engine speed (600 and 900 rpm) is analyzed for both fuels. It is shown that the sensitivity of a fuel’s autoignition characteristics to λ and engine speed significantly differs at BRON and BMON conditions. Moreover, this sensitivity is found to vary among fuels, despite the same RON. It is also observed that the presence of low-temperature heat release (LTHR) under BRON condition leads to more sequential autoignition and longer combustion duration than BMON condition. Finally, the study indicates that the octane index (OI) fails to capture the trend in the variation of autoignition propensity with S under the BMON condition.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Modeling Electrokinetics of Oxygen Electrodes in Solid Oxide Electrolyzer Cells

A microscale model is presented in this study to simulate electrode kinetics of the oxygen electrode in a solid oxide electrolyzer cell (SOEC). Two mixed ionic/electronic conducting structures are examined for the oxygen producing electrode in this work: single layer porous lanthanum strontium cobalt ferrite (LSCF), and bilayer LSCF/SCT (strontium cobalt tantalum oxide) structures. A yttrium-stabilized zirconia (YSZ) electrolyte separates the hydrogen and oxygen electrodes, as well as a gadolinium doped-ceria (GDC) buffer layer on the oxygen electrode side. Electrochemical reactions occurring at the two-phase boundaries (2PBs) and three-phase boundaries (3PBs) of single-layer LSCF and bilayer LSCF/SCT oxygen electrodes are modeled under various SOEC voltages with lattice oxygen stoichiometry as the key output. The results reveal that there exists a competition in electrode kinetics between 2PBs and 3PBs, but 3PBs are the primary reactive sites for single-layer LSCF oxygen electrode under high voltages. These locations experience the greatest oxygen stoichiometry variations and are therefore the most likely locations for dimensional changes. By applying an active SCT layer over LSCF, the 2PBs become activated to compete with the 3PBs, thus alleviating oxygen stoichiometry variations and reducing the likelihood of dimensional change. This strategy could reduce lattice structural expansion, proving to be valuable for electrode-electrolyte delamination prevention and will be the focus of future work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Formation and Diagnostics of Collisionless Shocks in Laser-Produced Plasmas (Final Report)

Collisions of supersonic flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large; therefore, most shock waves in astrophysics are “collisionless” and form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at DOE high energy density (HED) facilities can achieve the conditions necessary for the formation of collisionless shocks and will provide a unique avenue for studying the nonlinear physics of shock waves. This project investigates the parameters and experimental configurations needed for the creation of collisionless shocks in the laboratory experiments. This study is being performed with state of the art 3D kinetic plasma simulations. Recent advances in ab-initio particle-in-cell (PIC) simulations of astrophysical shocks enable us to simulate shock formation on timescales and in geometries relevant to the laboratory experiments. The main goal of the proposal is to ascertain which signatures of shock physics can be extracted from the realistic experiments, accounting for limitations imposed by finite spatial and time scales in the lab. To aid the interpretation of current and future experimental data, PIC simulation data have been analyzed using the same diagnostic techniques as used in the laboratory, such as proton radiography and neutronics, to unambiguously predict the experimental signatures of relevant collisionless shock phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

w22_laser_radiography - Laser-Based MeV Radiography

This project seeks to develop novel x-ray sources based on high-power short-pulse lasers with unique advantages, including small spot size, flexible configuration, and tailorable spectra. Kinetic plasma simulations using the VPIC particle in cell code will be performed to support the R&D goals of this project. If successful, this effort will allow LANL to deploy flexible, compact x-ray sources for static tomographic radiography and dynamic multi-axis radiography at our firing sites and facilities such as pRad and U1a. The higher resolution, shorter pulse duration, and tunability of our proposed sources will enable new applications in the weapons program and broader scientific community. This effort advances the LANL FY22 Strategic Investment Plan priority of “[r]adiographic tools for dynamic and static imaging with improved spatial and temporal resolution, sensitivity, and dose.”

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Creation Synthetic Data to Train a Digital Twin to Predict Reactor Operations

Understanding techniques to strengthen the nuclear safeguards regime is crucial in preventing nuclear proliferation due to recent advancements in the nuclear energy industry such as Generation IV reactors and microreactors. Prior to the construction of a nuclear power plant, it is necessary to understand the proliferation potential of the plant's reactor. Digital twins serve as a unique solution to recognizing reactor behavior indicative of nuclear proliferation. A digital twin is defined as a virtual model that works in unison to represent a physical asset, with a transference of data between the virtual and physical assets [1]. This work serves as validation for training a digital twin on synthetic data fabricated via means of Serpent reactor physics and point kinetics equations simulations. In this case this work is based on parameters of Idaho State University's AGN-201 reactor. The synthetic data can then be utilized to train machine learning models in the future to further investigate the utility of these methods. The accuracy of the predicted data is measured against real operational data to verify the reliability of the synthetic data creation methods and decide whether these methods should be used in the future to inform inspectors of a reactor's proliferation potential.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Exploring the Acceleration Mechanisms for Particle Injection and Power-law Formation during Transrelativistic Magnetic Reconnection

Magnetic reconnection in the relativistic and transrelativistic regimes is able to accelerate particles to hard power-law energy spectra f ∝ γ -p (approaching p = 1). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during magnetic reconnection in the transrelativistic regime of a proton–electron plasma. While earlier work in this parameter regime has focused on the effects of electric field parallel to the local magnetic field on the particle injection (from thermal energy to the lower-energy bound of the power-law spectrum), here we examine the roles of both parallel and perpendicular electric fields to gain a more complete understanding on the injection process and further development of a power-law spectrum. We show that the parallel electric field does contribute significantly to particle injection, and is more important in the initial phase of magnetic reconnection. However, as the simulation proceeds, the acceleration by the perpendicular electric field becomes more important for particle injection and completely dominates the acceleration responsible for the high-energy power-law spectrum. This holds robustly, in particular for longer reconnection times and larger systems, i.e., in simulations that are more indicative of the processes in astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Fast Dissipation of Colliding Alfvén Waves in a Magnetically Dominated Plasma

Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object's luminosity. We describe a dissipation mechanism that works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields B 1 and B 2 that propagate in opposite directions along the background magnetic field B 0 and collide. The collision forms a thin current sheet perpendicular to B 0 , which absorbs the incoming waves. The current sheet is sustained by an electric field E breaking the magnetohydrodynamic condition E < B and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude A = B 1 / B 0 = B 2 /B 0 propagating in a uniform B 0 . The mechanism is activated when A > 1/2. It dissipates a large fraction of the wave energy, † = (2A - 1)/A 2 , reaching 100% when A = 1. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate A of main interest, the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at A $\gg$ 1. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A Structure-Preserving Decorated Particle Method for the Vlasov-Poisson System

We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.

65M75, 70H05, 70G65↗
DOE OSTI2021

A novel measurement of marginal Alfven Eigenmode stability during high power auxiliary heating in JET

The interaction of Alfven Eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25-250 kHz and toroidal mode numbers |n| < 20. During the 2019-2020 deuterium campaign, almost 7500 resonances and their frequencies f0, net damping rates \gamma < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient's effect on accessibility and ELM-related noise's impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized Ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ICRH and NBI) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f0, \gamma, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stout Smearing on a Quantum Computer

Smearing of gauge-field configurations in lattice field theory improves the results of lattice simulations by suppressing high energy modes from correlation functions. In quantum simulations, high kinetic energy eigenstates are introduced when the time evolution operator is approximated such as Trotterization. While improved Trotter product formulae exist to reduce the errors, they have diminishing accuracy returns with respect to resource costs. Therefore having an algorithm that has fewer resources than an improved Trotter formula is desirable. In this work I develop a representation agnostic method for quantum smearing and show that it reduces the coupling to high energy modes in the discrete nonabelian gauge theory $D_4$

Gustafson, Erik J.↗