Search NASA⌕ Search

SEARCH · Search NASA

Results for “transport coefficients”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transport-driven toroidal rotation with general viscosity profile

Abstract Using the assumption of a weak normalized turbulent viscosity, usually valid in practice, the modulated-transport model (Stoltzfus-Dueck 2012 Phys. Plasmas 19 055908) is generalized to allow the turbulent transport coefficient to vary in an arbitrary way on radial and poloidal position. The new approach clarifies the physical interpretation of the earlier results and significantly simplifies the calculation, via a boundary-layer asymptotic method. Rigorous detailed appendices verify the result of the simple boundary-layer calculation, also demonstrating that it achieves the claimed order of accuracy and providing a concrete prediction for the strong plasma flows in the immediate vicinity of the last closed flux surface. The new formulas are used to predict plasma rotation at the core-edge boundary, in cases with and without externally applied torque. Dimensional formulas and extensive discussion are provided, to support experimental application of the new model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Density pedestal prediction model for tokamak plasmas

Abstract A model for the pedestal density prediction based on neutral penetration combined with pedestal transport is presented. The model is tested against a pedestal database of JET-ILW Type I ELMy H-modes showing good agreement over a wide range of parameters both in standalone modelling (using the experimental temperature profile) and in full Europed modelling that predicts both density and temperature pedestals simultaneously. The model is further tested for ASDEX Upgrade and MAST-U Type I ELMy H-modes and both are found to agree with the same model parameters as for JET-ILW. The JET-ILW experiment where the isotope of the main ion is varied in a D/T scan at constant gas rate and constant β N is successfully modelled as long as the separatrix density ( n e,sep ) and pedestal transport coefficient ratio ( D / χ ) are varied in accordance with the experimentally observed variation of n e,sep and the isotope dependence of D / χ found in gyrokinetic simulations. The predictions are found to be sensitive to n e,sep which is why the model is combined with an n e,sep model to predict the pedestal for the STEP fusion reactor.

Physics↗
DOE OSTI2024

Roles of non-axisymmetric perturbations in free drift vertical displacement events on EAST

Abstract The safe operation of most tokamaks, especially the large ones, relies on the feedback control of vertical displacement events (VDEs). However, most of these feedback control systems are based on axisymmetric VDE models. In this study, we use NIMROD simulations to study the role of non-axisymmetric perturbations in free drift vertical displacement events on EAST. The high- n modes in the non-axisymmetric VDE grow first, which drives the formation of high- n magnetic island chains. Subsequently, the magnetic island chains grow and overlap with each other, leading to the destruction of the magnetic flux surface, which induces a minor disruption and accelerates the start of the following major disruption. The magnetic island and the stochastic magnetic field allow the toroidally non-axisymmetric poloidal plasma current to jet towards the hoop force direction, forming finger-like and filamentary structures. Such a plasma current non-axisymmetry strongly depends on the anisotropy in the thermal transport coefficients.

Physics↗
DOE OSTI2025

Power balance and divertor asymmetries in the Super-X divertors of MAST-U using SOLPS-ITER

Spherical tokamaks (STs) present unique challenges and opportunities in the area of particle and power exhaust, intensified due to their more compact sizes. Substantial efforts are underway in STs to determine the limits in dissipative operational regimes and advanced divertor solutions, including at MAST-U which provides access to the Super-X divertor configuration. Power balance, and upper/lower divertor asymmetries have been studied using SOLPS-ITER simulations of the MAST-U Super-X divertor. A set of simulations with experimentally inferred transport coefficients with E x B and diamagnetic drifts activated, consisting of density and power scans, and high field side vs low field side gas puff locations, have been used for code experimentation to uncover trends beyond the current experimental parameter space. The upper biased asymmetry (U:L > 1) of the ratio of the peaks of the plasma energy flux densities at the outer targets increases with heating power and decreases with gas puff strength, going from symmetric to up to a factor of 15. The upper target electron temperature has been found to be a good ordering quantity for the magnitude of this asymmetry for all heating powers, gas puff strength, and gas puff locations. The lower divertor biased asymmetry (U:L < 1) of the radiation patterns processed through SOLPS-based bolometry synthetic diagnostics is in qualitative agreement with resistive bolometry experimental results, and it is in quantitative agreement with the trend of the total volume radiation within the divertors of SOLPS. However, radiation measurements alone are not sufficient to infer the magnitude of the asymmetry of the peaks of the power loads at the targets.

MAST-U↗
DOE OSTI2026

Physics insights from a large-scale 2D UEDGE simulation database for detachment control in KSTAR

A large-scale database of two-dimensional UEDGE simulations has been developed to study detachment physics in KSTAR and to support surrogate models for control applications. Nearly 70 000 steady-state solutions were generated, systematically scanning upstream density, input power, plasma current, impurity fraction, and anomalous transport coefficients, with magnetic and electric drifts across the magnetic field included. The database identifies robust detachment indicators, with strike-point electron temperature at detachment onset consistently T e,target ~ 3-4 eV, largely insensitive to upstream conditions. Scaling relations reveal weaker impurity sensitivity than one-dimensional models and show that heat flux widths follow Eich’s scaling only for uniform, low D and χ. Distinctive in–out divertor asymmetries are observed in KSTAR, differing qualitatively from DIII-D. Complementary time-dependent simulations quantify plasma response to gas puffing, with delays of 5 - 15 ms at the outer strike point and ∼ 40 ms for the low-magnetic-field-side radiation front. These dynamics are well captured by first-order-plus-dead-time models and are consistent with experimentally observed detachment-control behavior in KSTAR (Gupta et al 2025 Plasma Phys. Control. Fusion (submitted)).

Physics - Plasma physics↗
DOE OSTI2022

Integrability breaking in the Rule 54 cellular automaton

Cellular automata have recently attracted a lot of attention as testbeds to explore the emergence of many-body quantum chaos and hydrodynamics. We consider the Rule 54 model, one of the simplest interacting integrable models featuring two species of quasiparticles (solitons), in the presence of an integrability-breaking perturbation that allows solitons to backscatter. Here, we study the onset of thermalization and diffusive hydrodynamics in this model, compute perturbatively the diffusion constant of tracer particles, and comment on its relation to transport coefficients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Self-consistent and detailed opacities from a non-equilibrium average-atom model

Modern density functional theory (DFT) is a powerful tool for accurately predicting self-consistent material properties such as equations of state, transport coefficients and opacities in high energy density plasmas, but it is generally restricted to conditions of local thermodynamic equilibrium (LTE) and produces only averaged electronic states instead of detailed configurations. We propose a simple modification to the bound-state occupation factor of a DFT-based average-atom model that captures essential non-LTE effects in plasmas—including autoionization and dielectronic recombination—thus extending DFT-based models to new regimes. Finally, we then expand the self-consistent electronic orbitals of the non-LTE DFT-AA model to generate multi-configuration electronic structure and detailed opacity spectra. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗
DOE OSTI2025

Modified coherence and the transverse extent of jets

We present a study of the transverse size of parton showers and their diminishing interaction with the medium in the high virtuality stage of jet evolution. We consider the process of a hard quark produced in deep inelastic scattering off a large nucleus. Single gluon radiation from this quark, in the absence of scattering, is rederived using wave packets. This allows for a derivation of the quantum uncertainty size of the hard quark, at the point of splitting. This uncertainty size is then incorporated within a Monte Carlo shower routine yielding transverse shower sizes noticeably larger than the classical antenna size of the shower. No clear relation is found between the full uncertainty size of the shower and the virtuality of the originating parton. The single gluon emission from the hard quark is then reanalyzed for the case of single rescattering off the remainder of the nucleus. A relation is derived between the jet transport coefficient $\hat{𝑞}$ and the gluon Transverse momentum-dependent parton distribution function (gTMDPDF). Solving this relation, for a simple case, clearly demonstrates the weakening of $\hat{𝑞}$ with the virtuality of the hard splitting parton.

jet quenching↗
DOE OSTI2025

Collision-Energy Dependence in Heavy-Ion Collisions from Nonlinear QCD Evolution

We explore the effects of including the energy dependence determined from evolution equations within the color glass condensate framework on observables in ultrarelativistic heavy-ion collisions. This amounts to integrating the JIMWLK evolution equations into the impact parameter-dependent glasma model, which is then coupled to viscous relativistic hydrodynamics. This methodology allows for a systematic representation of nuclei at specific Bjorken-x values, which are probed at different center-of-mass energies of the collision and rapidities of final state particles. Comparing the methodology to the conventional impact parameter-dependent glasma model, we find significant effects on multiplicity distributions and particle spectra, especially in smaller collision systems at the highest center-of-mass energies. Our results highlight the importance of incorporating nonlinear QCD evolution in the description of heavy-ion collisions at varying center-of-mass energies, as the precise extraction of transport coefficients will be affected. This Letter establishes a robust framework for understanding the quark gluon plasma and nuclear structure at high energy, integrating small-x physics into the initial conditions of heavy-ion collisions.

Mäntysaari, Heikki↗
DOE OSTI2020

Relaxation of the degenerate one-dimensional Fermi gas

In this work, we study how a system of one-dimensional spin-1/2 fermions at temperatures well below the Fermi energy approaches thermal equilibrium. The interactions between fermions are assumed to be weak and are accounted for within the perturbation theory. In the absence of an external magnetic field, spin degeneracy strongly affects relaxation of the Fermi gas. For sufficiently short-range interactions, the rate of relaxation scales linearly with temperature. Focusing on the case of the system near equilibrium, we linearize the collision integral and find the exact solution of the resulting relaxation problem. Furthermore, we discuss the application of our results to the evaluation of the transport coefficients of the one-dimensional Fermi gas.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Unified ab initio description of Fröhlich electron-phonon interactions in two-dimensional and three-dimensional materials

Ab initio calculations of electron-phonon interactions including the polar Fröhlich coupling have advanced considerably in recent years. The Fröhlich electron-phonon matrix element is by now well understood in the case of bulk three-dimensional (3D) materials. In the case of two-dimensional (2D) materials, the standard procedure to include Fröhlich coupling is to employ Coulomb truncation, so as to eliminate artificial interactions between periodic images of the 2D layer. While these techniques are well established, the transition of the Fröhlich coupling from three to two dimensions has not been investigated. Furthermore, it remains unclear what error one makes when describing 2D systems using the standard bulk formalism in a periodic supercell geometry. In this work, we generalize previous work on the ab initio Fröhlich electron-phonon matrix element in bulk materials by investigating the electrostatic potential of atomic dipoles in a periodic supercell consisting of a 2D material and a continuum dielectric slab. We obtain a unified expression for the matrix element, which reduces to the existing formulas for three-dimensional and 2D systems when the interlayer separation tends to zero or infinity, respectively. This expression enables an accurate description of the Fröhlich matrix element in 2D systems without resorting to Coulomb truncation. We validate our approach by direct ab initio density-functional perturbation theory calculations for monolayer BN and MoS 2 , and we provide a simple expression for the 2D Fröhlich matrix element that can be used in model Hamiltonian approaches. The formalism outlined in this work may find applications in calculations of polarons, quasiparticle renormalization, transport coefficients, and superconductivity, in 2D and quasi-2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hydrodynamic electron transport in graphene Hall-bar devices

Here we consider hydrodynamic electron transport in the Hall-bar geometry. The theory is developed for systems with non-Galilean-invariant electron liquids. We show that inhomogeneity of the electron density induced by long-range disorder and gating leads to mixing between the hydrodynamic transport mode and transport relative to the electron liquid. For graphene systems near charge neutrality, these effects lead to strong coupling of the hydrodynamic flow to charge transport. As a result, the effective electrical conductivity of the system may significantly exceed the intrinsic conductivity of the electron liquid. We obtain analytic expressions for the thermoelectric transport coefficients of the system as a function of density in the full crossover region between clean and disorder-dominated regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Bayesian calibration of a hybrid nuclear collision model using $\textit{p}$₋Pb and Pb-Pb data at energies available at the CERN Large Hadron Collider

Here, we posit a unified hydrodynamic and microscopic description of the quark-gluon plasma (QGP) produced in ultrarelativistic p₋Pb and Pb-Pb collisions at $\sqrt{s}_\text{NN}$ = 5.02 TeV and evaluate our assertion using Bayesian inference. Specifically, we model the dynamics of both collision systems using initial conditions with parametric nucleon substructure, a preequilibrium free streaming stage, event-by-event viscous hydrodynamics, and a microscopic hadronic afterburner. Free parameters of the model, which describe the initial state and QGP medium are then simultaneously calibrated to fit charged-particle yields, mean $p_T$, and flow cumulants. We argue that the global agreement of the calibrated model with the experimental data strongly supports the existence of hydrodynamic flow in small collision systems at ultrarelativistic energies, and that the flow produced develops at length scales smaller than a single proton. Posterior estimates for the model's input parameters are obtained, and new insights into the temperature dependence of the QGP transport coefficients and event-by-event structure of the proton are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Energy and scale dependence of $\hat{q}$ and the “JET puzzle”

Here, we present an attempt to probe the underlying structure of the quark-gluon plasma (QGP) at high resolution, based on the extracted jet transport coefficient $\hat{q}$. We argue that the exchanged momentum $\textit{k}$ between the hard parton and the medium varies over a range of scales, and for $\textit{k}$ ≥ 1 GeV, $\hat{q}$ can be expressed in terms of a parton distribution function (PDF). Because the mass of a QGP constituent is unknown, we define a scaling variable $x_N$ to represent the ratio of the parton momentum to the momentum of a self-contained section of the plasma which has a mass of 1 GeV. This scaling variable is used to parametrize the QGP-PDF. Calculations based on this reconstructed $\hat{q}$ are compared to data sensitive to the hardcore of jets, i.e., the single hadron suppression in terms of the nuclear modification factor $R_{ A A}$ and the azimuthal anisotropy parameter $v_2$ as a function of transverse momentum $p_T$, centrality, and energy of the collision. It is demonstrated that the scale evolution of the QGP-PDF is responsible for the reduction in the normalization of $\hat{q}$ between fits to the Relativistic Heavy-Ion Collider and the Large Hadron Collider data; a puzzle, first discovered by the JET Collaboration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Minijet quenching in a concurrent jet + hydro evolution and the nonequilibrium quark-gluon plasma

Minijets, created by perturbative hard QCD collisions at moderate energies, can represent a significant portion of the total multiplicity of a heavy-ion collision event. Since their transverse momenta are initially larger than the typical saturation scale describing the bulk of the equilibrating quark-gluon plasma (QGP), they ought to be described through the physics of parton energy loss. Indeed, their typical stopping distances are larger than the usual hydrodynamization time, so they do not in general hydrodynamize at the same pace than the bulk of the collision. Therefore, in general minijets cannot be described solely by a unique preequilibrium stage that bridges the initial, overoccupied glasma state, with the hydrodynamical evolution. In this work we make use of a new concurrent minijet+hydrodynamic framework in which the properties of the hydrodynamically evolving QGP are modified due to the injection of energy and momentum from the minijets. We study the system for different choices of the minimum transverse momentum associated with minijet production. In order to achieve a realistic description of charged particle multiplicity, the amount of entropy associated to the low- x initial state needs to be reduced. Moreover, the fact that the injected momentum from the randomly oriented minijets is not correlated with the spatial gradients of the system reduces overall flow, and the value of the QGP transport coefficients needs to be reduced accordingly in order to describe the measured flow coefficients in experiments. They are, in effect, an important new source of fluctuations, resulting in a spikier, notably modified hydrodynamical evolution when compared to the scenario in which the presence of minijets is ignored. We avow that their abundance makes it necessary to include their physics in holistic descriptions of heavy-ion collisions. We discuss the impact of the minijets on a number of observables, such as p T spectra and p T -differential flow v n for a wide range of centrality classes. In contrast to elliptic, triangular or quadrangular flow, here we find that directed flow, v 1 , has the strongest potential to discriminate between different minijet production rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Comprehensive simulation of heavy-ion collisions at nonzero baryon chemical potential

Here, we present results of hydrodynamic modeling of Au-Au collisions from $\sqrt{s_{NN}}$ = 7.7 to 200 GeV. Our simulations have three novel components. First, we use a Linear EXtrapolation of Ultratrelativistic nucleon-nucleon Scattering to nucleus-nucleus collisions (LEXUS) inspired Monte Carlo initial-state model. Second, we use a crossover equation of state at finite baryon densities without a critical point. Finally, we use departure functions derived from the quasiparticle theory of transport coefficients for hadronic matter at nonzero baryon densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗