Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Simulations of saturated MHD activity in the HBT-EP tokamak

The NIMROD code is used to perform simulations of MHD activity in the HBT-EP tokamak, including the effect of a resistive wall. Linear simulations are conducted to scan current and pressure profiles to establish self-consistent baseline equilibria that result in MHD instabilities with low error between experimental diagnostics and synthetic equivalents. A family of such equilibria is found with relatively peaked current and pressure profiles. Nonlinear simulations are performed using the optimal equilibrium, which is linearly unstable to a Resistive Wall Mode—stable in the ideal wall limit and with an Alfvénic growth rate in the no-wall limit. Using resistive wall boundary conditions, nonlinear Hall MHD simulations yield a rotating, saturated n=1 mode similar to the type observed experimentally in HBT-EP discharges. An n=1 mode around the q=2 surface mediates a localized axisymmetric perturbation that flattens the plasma current from a linearly unstable state to a 3D, rotating, stable state. During the saturated period, periodic changes in inter-diagnostic phase lag allow for qualitative estimations of effective, macroscopic transport coefficients used in MHD closure models and could be extended to model post-disruption MHD structure rotation after thermal quenches seen in HBT-EP.

Arnold, David A. (ORCID:0009000059535761)↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Generation and optimization of cross-sections for electron-C 4 F 7 N collisions

Abstract A complete and consistent set of electron-neutral collision cross-sections for the novel insulating gas C 4 F 7 N is reported. The set is composed of a combination of cross-sections previously reported in literature, optimized via a genetic algorithm in conjunction with a multi-term Boltzmann equation solver, and calculatedab initiousing the R-matrix code Quantemol-EC. The finalized set accurately reproduces reported macroscopic rate and transport coefficients as well as Townsend coefficients and critical electric field strengths in C 4 F 7 N and its mixtures with nitrogen, carbon dioxide, and argon.

Physics↗
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↗