Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plasma kinetic 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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Perturbed Ion Temperature and Toroidal Flow Profile Measurements in Rotating Neoclassical Tearing Mode Magnetic Islands

The perturbed ion temperature and toroidal flow were measured in rotating Neoclassical Tearing Modes (NTM) in a tokamak for the first time. These toroidally and radially resolved profiles were obtained by impurity ion spectroscopy in a 2,1 NTM in DIII-D. In agreement with drift-kinetic simulations, the electron temperature profile is flat, while the ion temperature gradient is restored across the magnetic island O-point in the presence of fast ions; the perturbed flow has minima in the O-points and maxima at the X-points. Furthermore, these measurements provide the first confirmation of the theoretically expected ion temperature and flow response to a magnetic island needed to predict the NTM on set threshold scaling for ITER and other future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Comparing XRISM Cluster Velocity Dispersions with Predictions from Cosmological Simulations: Are Feedback Models Too Ejective?

The dynamics of the intracluster medium (ICM), the hot plasma that fills galaxy clusters, are shaped by gravity-driven cluster mergers and feedback from supermassive black holes (SMBHs) in the cluster cores. XRISM measurements of ICM velocities in several clusters offer insights into these processes. We compare XRISM measurements for nine galaxy clusters (Virgo, Perseus, Centaurus, Hydra-A, PKS 0745–19, A2029, Coma, A2319, and Ophiuchus) with predictions from three state-of-the-art cosmological simulation suites, TNG-Cluster, the Three Hundred Project GADGET-X, and GIZMO-SIMBA, that employ different models of feedback. In cool cores, XRISM reveals systematically lower velocity dispersions than the simulations predict, with all 10 measurements below the median simulated values by a factor of 1.5–1.7 on average and all falling within the bottom 10% of the predicted distributions. The observed kinetic-to-total pressure ratio is also lower, with a median value of 2.2%, compared to the predicted 5.0%–6.5% for the three simulations. Outside the cool cores and in non-cool-core (NCC) clusters, simulations show better agreement with XRISM measurements, except for the outskirts of the relaxed, cool-core cluster A2029, which exhibits an exceptionally low kinetic pressure support (<1%), with none of the simulated systems in either of the three suites reaching such low levels. The NCC Coma and A2319 exhibit dispersions at the lower end but within the simulated spread. Our comparison suggests that the three numerical models may overestimate the kinetic effects of SMBH feedback in cluster cores. Additional XRISM observations of NCC clusters will clarify if there is a systematic tension in the gravity-dominated regime as well.

galaxy clusters↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Collaborative Research: Vlasov–Maxwell Simulations to resolve electron heating and dissipation, in quasiperpendicular shocks (Final Technical Report)

The project is to study electron heating at quasi-perpendicular shocks. We have analyzed measurements from the Magnetosphere Multiscale (MMS) spacecraft and particle-in-cell simulations to address the question. One key discovery is that waves and are excited in the shock transition region, which further develop into complicated structures and thin current layers, which contribute to accelerate and heat electrons. We further investigated the excitation and evolution of these waves and structures, and their association with ion and electron dynamics. In addition, the analysis technique of Vlasov equation terms is further applied to study electron dynamics at the dayside magnetopause, to understand the effect of density, velocity and temperature variations on electron distributions and kinetics.

58 GEOSCIENCES↗
DOE OSTI2022

The density distribution and physical origins of intermittency in supersonic, highly magnetized turbulence with diverse modes of driving

The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Second Year Report/ Final Report)

This research will focus on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and ow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal ow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aim to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

A Local Macroscopic Conservative (LoMaC) Low Rank Tensor Method for the Vlasov Dynamics

Abstract In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum and energy at the discrete level. This is a follow-up work of our previous development of a conservative low rank tensor approach for Vlasov dynamics ( arXiv:2201.10397 ). In that work, we applied a low rank tensor method with a conservative singular value decomposition to the high dimensional VP system to mitigate the curse of dimensionality, while maintaining the local conservation of mass and momentum. However, energy conservation is not guaranteed, which is a critical property to avoid unphysical plasma self-heating or cooling. The new ingredient in the LoMaC low rank tensor algorithm is that we simultaneously evolve the macroscopic conservation laws of mass, momentum and energy using a flux-difference form with kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables by a conservative orthogonal projection. The algorithm is extended to the high dimensional problems by hierarchical Tuck decomposition of solution tensors and a corresponding conservative projection algorithm. Extensive numerical tests on the VP system are showcased for the algorithm’s efficacy.

Guo, Wei↗
DOE OSTI2022

Measurement and simulation of small cryogenic neon pellet Ne-I 640 nm photon efficiency during ablation in DIII-D plasma

Small (~1 mm) neon pellet fragments are fired into DIII-D H-mode plasmas, and resulting trajectory-averaged photon efficiencies (neutral neon ionization events for every photon emitted) of S/XB ≈ 85 are estimated for Ne-I 640 nm by dividing the estimated initial pellet fragment mass by the measured number of emitted Ne-I photons. The experiments are modeled by running the Lagrangian particle (LP) fluid/magneto-hydrodynamic pellet code to estimate axial ablation plume neon density profiles and temperature profiles at each pellet position. These solutions are then fed into the PrismSPECT collisional-radiative code, which calculates resulting neon charge states and photon emission rates, giving a profile-average of S/XB ≈ 109. The burnthrough plasma minor radius predicted by LP (ρ ≈ 0.63) is reasonably close to the experimental observation ρ ≈ 0.6. Furthermore, the modeling indicates that local S/XB is not constant along the pellet trajectory but tends to increase with increasing ablation rate. Non-equilibrium kinetics are predicted to be very important, while line trapping is predicted to be relatively unimportant (for Ne-I 640 nm S/XB).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Study of the Transition to Turbulence in a Bed of 67 Spherical Pebbles

Packed beds are commonly found in many engineering systems and have been widely studied for decades. A relatively new packed bed system is the Pebble Bed Reactor, a type of generation-IV nuclear reactor. Unlike many of the packed beds encountered in chemical and process engineering applications, Pebble Bed Reactors are larger and operate at significantly higher Reynolds numbers. As a result of these differences, there is a very limited amount of information on the detailed flow physics that exist in these complex geometries. This work seeks to contribute to a growing database of flow data for Pebble Bed Reactor systems by performing Direct Numerical Simulations of the flow in an experimental bed of 67 pebbles for a range of conditions. Simulations are performed at a Prandtl number of 0.66 and Reynolds numbers from 300–600. These Reynolds numbers are chosen to gain additional knowledge on the spatial development of turbulence in these systems. Analysis of the Turbulent Kinetic Energy, turbulence anisotropy, and Turbulent Heat Flux is performed. Results demonstrate significant development of the TKE across the tested range of Reynolds numbers. Examination of both the TKE and THF reveal that development first occurs near the center of the bed and propagates radially as the flow moves further into the bed. Notable regions of negative production of turbulent kinetic energy are observed in regions where flow accelerates around pebble contact points. Furthermore, these regions are found to coincide with regions of 1-component turbulence.Kindly check and confirm, all authors email id is correctly identified.These are correct

Direct numberical simulation↗
DOE OSTI2025

Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code

The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Broadening of the drumhead-mode spectrum due to in-plane thermal fluctuations of two-dimensional trapped ion crystals in a Penning trap

Two-dimensional crystals of ions stored in Penning traps are a leading platform for quantum simulation and sensing experiments. For small amplitudes, the out-of-plane motion of such crystals can be described by a discrete set of normal modes called the drumhead modes, which can be used to implement a range of quantum information protocols. However, experimental observations of crystals with Doppler-cooled and even nearground- state-cooled drumhead modes reveal an unresolved drumhead-mode spectrum. In this work, we establish in-plane thermal fluctuations in ion positions as a major contributor to the broadening of the drumhead-mode spectrum. In the process, we demonstrate how the confining magnetic field leads to unconventional in-plane normal modes, whose average potential and kinetic energies are not equal. This property, in turn, has implications for the sampling procedure required to choose the in-plane initial conditions for molecular-dynamics simulations. Further, for current operating conditions of the NIST Penning trap, our study suggests that the two-dimensional crystals produced in this trap undergo in-plane potential-energy fluctuations of the order of 10mK. Our study therefore motivates the need for designing improved techniques to cool the in-plane degrees of freedom.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Flow Channels and the Generation of Alfvénic Turbulence Along Storm‐Time Inner Magnetospheric Field‐Lines

Abstract A feature of Earth's storm‐time magnetosphere outside the plasmapause is the occurrence of broad‐spectrum Alfvénic fluctuations. In this letter observations from the Van Allen Probes are compared with 3‐D fluid‐kinetic simulations of an evolving convective flow channel to investigate the mechanisms generating the observed spectrum. It is shown how narrow channels of fast convection are unstable to the Kelvin‐Helmholtz instability which on closed field‐lines initiates a cascade to small scales. Sustained driving of the flow combined with reflection from the topside ionosphere leads to the generation of an intensified spectrum of electromagnetic structures having similar spectral and morphological characteristics to those observed. This process couples enhanced magnetospheric convection to kinetic scale electromagnetic fluctuations that drive particle transport, scattering and energization through the outer radiation belt and ring current during geomagnetic storms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗