Search NASA⌕ Search

SEARCH · Search NASA

Results for “Anisotropic plasmas”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Novel Approach to Fast Calculation of High-Order Q-Cumulants

The method of Q-cumulants has been shown as a powerful tool to study the fine details of the azimuthal anisotropies in high-energy nucleus-nucleus collisions. A new method for the fast calculation of arbitrary order Q-cumulant $v_{n}\{2k\}$ values, based on the partition of a non-negative integer l $\le$ m for calculation of the 2m-particle azimuthal correlations is presented in this paper. Unlike the standard Q-cumulants method in which the calculation of high-order multi-particle calculations is impractical, the newly proposed method enables easy calculation. The validity of the method is proven via a toy model that uses the elliptic power distribution to simulate anisotropic emission of particles. The method enables the study of fine details of the $v_{2}$ distribution, such as higher-order central moments of the $v_{2}$ distribution, as well as the hydrodynamic behavior of the Quark-Gluon Plasma.

FOS: Physical sciences↗
DOE OSTI2021

Diagnosing plasma magnetization in inertial confinement fusion implosions using secondary deuterium-tritium reactions

Diagnosing plasma magnetization in inertial confinement fusion implosions is important for understanding how magnetic fields affect implosion dynamics and to assess plasma conditions in magnetized implosion experiments. Secondary deuterium–tritium (DT) reactions provide two diagnostic signatures to infer neutron-averaged magnetization. Magnetically confining fusion tritons from deuterium–deuterium (DD) reactions in the hot spot increases their path lengths and energy loss, leading to an increase in the secondary DT reaction yield. In addition, the distribution of magnetically confined DD-triton is anisotropic, and this drives anisotropy in the secondary DT neutron spectra along different lines of sight. Implosion parameter space as well as sensitivity to the applied B-field, fuel ρR, temperature, and hot-spot shape will be examined using Monte Carlo and 2D radiation-magnetohydrodynamic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Progress with the 5D full-F continuum gyrokinetic code COGENT

COGENT is an Eulerian gyrokinetic code being developed for edge plasma modelling. The code is distinguished by the use of a high-order finite-volume (conservative) discretization combined with mapped multi-block grid technology. Our recent work is focused on the development of a 5D full-F COGENT version.Anumerical algorithm utilizing locally a field-aligned multi-block coordinate system is implemented to facilitate simulations of highly anisotropic microturbulence in the presence of a strong magnetic shear. In this approach, the toroidal direction is divided into blocks such that,within each block, the cells are field-aligned and a non-matching (non-conformal) grid interface is allowed at the block boundaries. In this paper we report on details of the numerical implementation and present preliminary results of verification studies performed for the case of the ion temperature gradient (ITG) instability in a sheared toroidal annulus geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP

Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of electrode porosity architecture on electrochemical performances of 1 mm-thick LiFePO4 binder-free Li-ion electrodes fabricated by Spark Plasma Sintering

Thick electrodes with high active material loadings have been intensively studied over the last couple of decades in pursuit of achieving high energy density systems. To optimize and enhance the electrochemical performance of such electrodes, one has to control the pore morphology by, for example, varying the pore size and shape, and the level of porosity. In the present work, the fabrication of thick binder-free LiFePO4 (LFP) electrodes with two different pore sizes (12 and 20 mu m) and porosities (21 vol% and 44 vol%) using Spark Plasma Sintering (SPS) and templating approach is reported. The well-controlled porous architecture inside the thick electrodes is realized by fine-tuning experimental parameters. The impact of porosity architecture on electrochemical performance is quantified and correlated with the 3D tortuosity values determined from both micro-computed tomography and electrochemical impedance-based experimental methods. Based on the micro-computed tomography data analysis, estimated tortuosity values along X, Y, and Z axes reveal an anisotropic effect perpendicularly to the SPS compression axis (Z-direction). This is particularly profoundly observed in the samples with larger pores (20 mu m). The correlation between morphological properties and the rate capability performance is established indicating that the capacity loss happens mainly due to the Li-ion diffusion limitations.

Cathode material↗
DOE OSTI2020

Machine Learning for Memory Reduction in the Implicit Monte Carlo Simulations of Thermal Radiative Transfer [Slides]

Project Goal: Use parametric Machine Learning methods in order to reduce memory requirements at checkpointing & restarting in the IMC simulations of Thermal Radiative Transfer using: Expectation Maximization and Weighted Gaussian Mixture Model-based approach for `particle-data compression', introduced in Plasma Physics to model Maxwellian particle distributions by Luis Chacon and Guangye Chen; Expectation Maximization with Weighted Hyper-Erlang Model in order to compress isotropic IMC particle data in the frequency domain; and Expectation Maximization and von Mises-Fisher Mixture Model for compression of anisotropic IMC particle data in the angular domain (work-in-progress).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Constraints on ion velocity distributions from fusion product spectroscopy

Abstract Recent inertial confinement fusion experiments have shown primary fusion spectral moments which are incompatible with a Maxwellian velocity distribution description. These results show that an ion kinetic description of the reacting ions is necessary. We develop a theoretical classification of non-Maxwellian ion velocity distributions using the spectral moments. At the mesoscopic level, a monoenergetic decomposition of the velocity distribution reveals there are constraints on the space of spectral moments accessible by isotropic distributions. General expressions for the directionally dependent spectral moments of anisotropic distributions are derived. At the macroscopic level, a distribution of fluid element velocities modifies the spectral moments in a constrained manner. Experimental observations can be compared to these constraints to identify the character and isotropy of the underlying reactant ion velocity distribution and determine if the plasma is hydrodynamic or kinetic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Impact of QCD energy evolution on observables in heavy-ion collisions

Here, we study how the inclusion of energy dependence as dictated by quantum chromodynamic (QCD) small-x evolution equations affects key observables in ultra-relativistic heavy-ion collisions. Specifically, we incorporate JIMWLK evolution into the IP-Glasma framework, which serves as the initial condition for a simulation pipeline that includes viscous relativistic hydrodynamics and a hadronic afterburner. This approach enables a consistent modeling of highly energetic nuclei across varying Bjorken-x values, which are relevant for different collision energies and rapidity regions. In comparison to the standard IP-Glasma setup without small-x evolution, we observe pronounced changes in particle multiplicities and spectral distributions, especially in smaller systems and at the highest available energies. We further explore effects on anisotropic flow observables and correlations between mean transverse momentum and elliptic flow. Our findings underscore the critical role of nonlinear QCD evolution in accurately modeling the early stages of heavy-ion collisions, as well as its implications for extracting transport properties of the quark-gluon plasma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Observation of Hydrodynamic Flows in Imploding Fusion Plasmas on the National Ignition Facility [plus supplemental information]

Inertial confinement fusion implosions designed to have minimal fluid motion at peak compression often show significant linear flows in the laboratory, attributable per simulations to percent-level imbalances in the laser drive illumination symmetry. In this work, we present experimental results which intentionally varied the mode 1 drive imbalance by up to 4% to test hydrodynamic predictions of flows and the resultant imploded core asymmetries and performance, as measured by a combination of DT neutron spectroscopy and high-resolution x-ray core imaging. Neutron yields decrease by up to 50%, and anisotropic neutron Doppler broadening increases by 20%, in agreement with simulations. Furthermore, a tracer jet from the capsule fill-tube perturbation that is entrained by the hot-spot flow confirms the average flow speeds deduced from neutron spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Surface Josephson plasma waves in a high-temperature superconductor

Abstract Electron density oscillations with acoustic dispersions and sustained at boundaries between different media provide information about surface and interface properties of heterostructures. In ultrathin metallic films these plasmonic excitations are heavily damped. Superconductivity is predicted to reduce dissipation allowing detection of these resonances. Emerging low-loss interface Cooper-pair waves have been studied before, however, the observation of surface-confined Josephson plasmons in highly anisotropic superconductors has remained elusive. Here, we report on generation and coupling to these excitations in an ultrathin single-crystal film of high-temperature superconductor La 1.85 Sr 0.15 CuO 4 . The film becomes brighter than Au below the critical temperature when probed with sub-gap THz photons. We show that the enhanced signal in the superconducting state, which can be visualized with a spatial resolution better than λ/3000, originates from near-field coupling of light to surface Josephson plasmons. Our results open a path towards non-invasive investigation of enhanced superconductivity in artificial multilayers, buried interface states in topological heterostructures, and non-linear phenomena in Josephson devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A Performance Portable, Fully Implicit Landau Collision Operator with Batched Linear Solvers

Modern accelerators use hierarchical parallel programming models that enable massive multithreading within a processing element (PE), with multiple PEs per device driven by traditional processes. Batching is a technique for exposing PE-level parallelism in algorithms that have traditionally run on MPI processes or multiple threads within a single process. Opportunities for batching arise in, for example, kinetic discretizations of magnetized plasmas where collisions are advanced in velocity space at each spatial point independently. This paper builds on previous work on a high-performance, fully nonlinear, Landau collision operator by batching the linear solver, as well as batching the spatial point problems and adding new support for multiple grids for multiscale, multispecies problems. An anisotropic relaxation verification test that agrees well with previously published results and analytical models is presented. The performance results from NVIDIA A100 and AMD MI250X nodes are presented with hardware utilization analysis for each architecture. Finally, the entire implicit Landau operator time advance is implemented in Kokkos for performance portability, running entirely on the device and is available in the PETSc numerical library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Axion-photon conversion in strongly magnetised plasmas

Axion dark matter can resonantly convert to photons in the magnetosphere of neutron stars, possibly giving rise to radio signals observable on Earth. This method for the indirect detection of axion dark matter has recently received significant attention in the literature. The calculation of the radio signal is complicated by a number of effects; most importantly, the gravitational infall of the axions onto the neutron star accelerates them to semi-relativistic speed, and the neutron star magnetosphere is highly anisotropic. Both of these factors complicate the calculation of the conversion of axions to photons. In this work, we present the first fully three-dimensional calculation of the axion-photon conversion in highly magnetised anisotropic media. Depending on the axion trajectory, this calculation leads to orders-of-magnitude differences in the conversion compared to the simplified one-dimensional calculation used so far in the literature, altering the directionality of the produced photons. Our results will have important implications for the radio signal one would observe in a telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Heat conduction in an irregular magnetic field. Part 2. Heat transport as a measure of the effective non-integrable volume

Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyse critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $\kappa _{\perp } \rightarrow 0$ , where $\kappa_{\perp}$ is the perpendicular diffusion coefficient. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Stability of beta-induced Alfvén eigenmodes (BAE) in DIII-D

Although the stability of ellipticity, toroidal and reversed-shear Alfvén eigenmodes (EAE, TAE, RSAE) are relatively well understood, less is known about the stability of lower-frequency modes such as the beta-induced Alfvén eigenmode (BAE) but, because they are often unstable in present devices and are implicated in fast-ion transport, understanding their stability is vital. BAE stability is studied in primarily weak or reversed shear DIII-D plasmas with sub-Alfvénic deuterium beams. Modes are classified based on electron cyclotron emission, beam emission spectroscopy, magnetics, and interferometer data. The study is limited to the initial two seconds of the discharge, where the evolving q profile provides an effective scan of the dependence of stability upon q. In a dedicated experiment, BAEs are unstable at times in the discharge when the minimum of the safety factor q min is close to a rational number. The observed mode frequencies are usually close to analytic estimates of the BAE accumulation point and the eigenfunction peaks in the vicinity of q min . Unstable BAEs usually occur in bursts that chirp rapidly in frequency. To isolate the importance of thermal and beam gradients in driving the modes, the beam and electron cyclotron heating power is altered for 50–100 ms durations in reproducible discharges. As expected from the resonance condition, BAEs depend sensitively on the beam power and injection geometry. Modes only persist for ~25 ms because the anisotropic beam population only interacts strongly with the modes over a relatively narrow range of q. Lastly, a database of over 1000 beam-heated discharges shows that BAEs are more likely to be unstable when the poloidal beta exceeds 0.5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Weibel instability beyond bi-Maxwellian anisotropy

The shape of the anisotropic velocity distribution function, beyond the realm of strict Maxwellians can play a significant role in determining the evolution of the Weibel instability dictating the dynamics of self-generated magnetic fields. For non-Maxwellian distribution functions, we show that the direction of the maximum growth rate wave vector changes with shape. In this work, we investigate different laser-plasma interaction model distributions which show that their Weibel generated magnetic fields may require closer scrutiny beyond the second moment (temperature) anisotropy ratio characterization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Metal-Insulation REBCO Pancake Coil With Solder Surface Shunt: Testing and Modeling

Rare-earth barium copper oxide (REBCO) high-temperature superconducting (HTS) magnets can generate very high magnetic fields with superior engineering current density and robust mechanical integrity. The no-insulation (NI) winding technique, including its metal-insulation (MI) variant, enables “self-protection” against localized overheating through inter-turn current sharing, improving defect tolerance and allowing temporary current overloads. In REBCO pancake coils, the current-sharing path is predominantly confined to the edges of the REBCO tapes due to insulating buffer layers. Applying a surface shunt using solder can effectively enhance current sharing in NI/MI coils via the REBCO tape edges. Here, in this paper, we present both experimental and numerical research on the current distribution in surface-shunt metal-insulation (SSMI) REBCO pancake coils. We tested a stainless steel co-wound MI REBCO double-pancake coil, before and after applying a soldered surface shunt, in a saturated liquid nitrogen bath. The SSMI coil was modeled using the H-formulation and rotated anisotropic resistivity to study its dynamic behavior, including inter-turn current sharing and intra-turn screening currents. Benchmarked by experimental data, the numerical model is able to study the defect tolerance of the SSMI coil during energization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗