Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boltzmann transport”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulating transient heat transfer in graphene at finite Knudsen number via the Boltzmann transport equation and molecular dynamics

The phonon Boltzmann transport equation (BTE) with the relaxation time approximation (RTA) scattering model is used to calculate transient temperature profiles in graphene, and the results are compared to analogous molecular dynamics (MD) simulations. For the BTE calculations, the phonon dispersion relation and frequency-dependent scattering rates are obtained from a combination of MD data and semi-empirical power-law expressions for the normal and Umklapp phonon lifetimes. The dimensions and initial temperature conditions of graphene are varied to study the size and temperature dependence of thermal transport physics at the mesoscopic scale. Good quantitative agreement to within 5% is found between the BTE and MD results, over a wide range of temperatures and lengthscales of the temperature variation in the graphene sheet. Small differences are attributed to the inaccuracy of the RTA as applied to graphene, and to neglecting four-phonon scattering in the BTE simulations. The present results may further understanding in applications such as the transient heating of nanoelectronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Thermodynamic and electron transport properties of Ca 3 Ru 2 O 7 from first-principles phonon calculations and Boltzmann transport theory

Here, this work demonstrates a first-principles-based approach to obtaining finite temperature thermal and electronic transport properties which can be employed to model and understand mesoscale structural evolution during electronic, magnetic, and structural phase transitions. A computationally tractable model was introduced to estimate electron relaxation time and its temperature dependence. The model is applied to Ca 3 Ru 2 O 7 with a focus on understanding its electrical resistivity across the electronic phase transition at 48 K. A quasiharmonic phonon approach to the lattice vibrations was employed to account for thermal expansion while the Boltzmann transport theory including spin-orbit coupling was used to calculate the electron-transport properties, including the temperature dependence of electrical conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

The elphbolt ab initio solver for the coupled electron-phonon Boltzmann transport equations

elphbolt is a modern Fortran (2018 standard) code for efficiently solving the coupled electron–phonon Boltzmann transport equations from first principles. Using results from density functional and density functional perturbation theory as inputs, it can calculate the effect of the non-equilibrium phonons on the electronic transport (phonon drag) and non-equilibrium electrons on the phononic transport (electron drag) in a fully self-consistent manner and obeying the constraints mandated by thermodynamics. It can calculate the lattice, charge, and thermoelectric transport coefficients for the temperature gradient and electric fields, and the effect of the mutual electron–phonon drag on these transport properties. The code fully exploits the symmetries of the crystal and the transport-active window to allow the sampling of extremely fine electron and phonon wave vector meshes required for accurately capturing the drag phenomena. The corray feature of modern Fortran, which offers native and convenient support for parallelization, is utilized. The code is compact, readable, well-documented, and extensible by design.

electron drag↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Linear Boltzmann transport for jet propagation in the quark-gluon plasma: Inelastic processes and jet modification

A linear Boltzmann transport (LBT) Monte Carlo model has been developed to describe jet propagation and interaction with the quark-gluon plasma (QGP) in relativistic heavy-ion collisions. A complete set of elastic-scattering processes and medium-induced gluon emissions based on the higher-twist formalism are incorporated for both jet shower and medium recoil partons. It has been employed to describe experimental data on large transverse momentum hadron and jet spectra, correlation and jet substructures in high-energy heavy-ion collisions. Here, we document in detail the structure of the model and validation of the Monte Carlo implementations of the physics processes in LBT, in particular, the inelastic process of medium-induced gluon radiation. We carry out a comprehensive examination of the jet-medium interaction as implemented in LBT through energy loss and momentum broadening of a single hard parton, the energy and transverse momentum transfer from leading partons to medium-induced gluons and jet-induced medium excitation, and medium modification of reconstructed jets in a static and uniform medium. With realistic and event-by-event hydrodynamic medium in heavy-ion collisions, we compute and compare with experimental data on the jet cone-size dependence of the single inclusive jet suppression at both the BNL Relativistic Heavy-Ion Collider (RHIC) and the CERN Large Hadron Collider (LHC), the dijet asymmetry at the LHC and γ-jet correlation at RHIC. Effects of medium-induced gluon emissions and jet-induced medium excitation on jet observables are systematically examined. Rescatterings of the radiated gluons and recoil partons with the QGP are found essential to account for the enhancement of soft particle yield toward the edge of the jet cone.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

QLBT: a linear Boltzmann transport model for heavy quarks in a quark-gluon plasma of quasi-particles

Abstract We develop a new heavy quark transport model, QLBT, to simulate the dynamical propagation of heavy quarks inside the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Our QLBT model is based on the linear Boltzmann transport (LBT) model with the ideal QGP replaced by a collection of quasi-particles to account for the non-perturbative interactions among quarks and gluons of the hot QGP. The thermal masses of quasi-particles are fitted to the equation of state from lattice QCD simulations using the Bayesian statistical analysis method. Combining QLBT with our advanced hybrid fragmentation-coalescence hadronization approach, we calculate the nuclear modification factor $$R_\mathrm {AA}$$ R AA and the elliptic flow $$v_2$$ v 2 of D mesons at the Relativistic Heavy-Ion Collider and the Large Hadron Collider. By comparing our QLBT calculation to the experimental data on the D meson $$R_\mathrm {AA}$$ R AA and $$v_2$$ v 2 , we extract the heavy quark transport parameter $$\hat{q}$$ q ^ and diffusion coefficient $$D_\mathrm {s}$$ D s in the temperature range of $$1-4~T_\mathrm {c}$$ 1 - 4 T c , and compare them with the lattice QCD results and other phenomenological studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A Stochastic Calculus Approach to Boltzmann Transport

Traditional Monte Carlo methods for particle transport utilize source iteration to express the solution, the flux density, of the transport equation as a Neumann series. Our contribution is to show that the particle paths simulated within source iteration are associated with the adjoint flux density and the adjoint particle paths are associated with the flux density. Here, we make our assertion rigorous through the use of stochastic calculus by representing the particle path used in source iteration as a solution to a stochastic differential equation (SDE). The solution to the adjoint Boltzmann equation is then expressed in terms of the same SDE, and the solution to the Boltzmann equation is expressed in terms of the SDE associated with the adjoint particle process. An important consequence is that the particle paths used within source iteration simultaneously provide Monte Carlo samples of the flux density and adjoint flux density in the detector and source regions, respectively. The significant practical implication is that particle trajectories can be reused to obtain both forward and adjoint quantities of interest. To the best our knowledge, the reuse of entire particles paths has not appeared in the literature. Monte Carlo simulations are presented to support the reuse of the particle paths.

Boltzmann transport↗
DOE OSTI2021

Coupled Boltzmann transport equations of heavy quarks and quarkonia in quark-gluon plasma

We develop a framework of coupled transport equations for open heavy flavor and quarkonium states, in order to describe their transport inside the quark-gluon plasma. Our framework is capable of studying simultaneously both open and hidden heavy flavor observables in heavy-ion collision experiments and can account for both, uncorrelated and correlated recombination. Our recombination implementation depends on real-time open heavy quark and antiquark distributions. We carry out consistency tests to show how the interplay among open heavy flavor transport, quarkonium dissociation and recombination drives the system to equilibrium. We then apply our framework to study bottomonium production in heavy-ion collisions. We include Υ(1S), Υ(2S), Υ(3S), χb (1P) and χb (2P) in the framework and take feed-down contributions during the hadronic gas stage into account. Cold nuclear matter effects are included by using nuclear parton distribution functions for the initial primordial heavy flavor production. A calibrated 2 + 1 dimensional viscous hydrodynamics is used to describe the bulk QCD medium. We calculate both the nuclear modification factor R AA of all bottomonia states and the azimuthal angular anisotropy coefficient v 2 of the Υ(1S) state and find that our results agree reasonably with experimental measurements. Our calculations indicate that correlated cross-talk recombination is an important production mechanism of bottomonium in current heavy-ion experiments. The importance of correlated recombination can be tested experimentally by measuring the ratio of R AA ( χb (1P)) and R AA (Υ(2S)).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part IV: PFNS and $\bar{ν}$ Uncertainty Propagation

An analytic benchmark with continuous-energy cross sections was previously derived to validate criticality calculations. Here, to extend the utility of the analytic benchmark to verify the implementation of $\bar{ν}$ and prompt fission neutron spectrum (PFNS) uncertainty propagation methods, new simplified forms that are dependent on the incident (fission-causing) neutron energy, as well as the outgoing neutron energy for the PFNS, are introduced in this work. The analytical forms for the flux and adjoint flux are derived for the extended benchmark and used to determine the 𝑘-eigenvalue sensitivity to $\bar{ν}$ and PFNS. The 𝑘-eigenvalue uncertainty due to $\bar{ν}$ and PFNS is calculated for the analytic benchmark using simplified$\bar{ν}$ and PFNS representations based on the ENDF-B/VIII.0 239 Pu evaluation. Because of the low sensitivity of the analytic benchmark to the physical PFNS, a nonphysical high-sensitivity PFNS is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Phoebe: a high-performance framework for solving phonon and electron Boltzmann transport equations

Understanding the electrical and thermal transport properties of materials is critical to the design of electronics, sensors, and energy conversion devices. Computational modeling can accurately predict material properties but, in order to be reliable, requires accurate descriptions of electron and phonon states and their interactions. While first-principles methods are capable of describing the energy spectrum of each carrier, using them to compute transport properties is still a formidable task, both computationally demanding and memory intensive, requiring integration of fine microscopic scattering details for estimation of macroscopic transport properties. To address this challenge, we present Phoebe—a newly developed software package that includes the effects of electron–phonon, phonon–phonon, boundary, and isotope scattering in computations of electrical and thermal transport properties of materials with a variety of available methods and approximations. This open source C++ code combines MPI-OpenMP hybrid parallelization with GPU acceleration and distributed memory structures to manage computational cost, allowing Phoebe to effectively take advantage of contemporary computing infrastructures. We demonstrate that Phoebe accurately and efficiently predicts a wide range of transport properties, opening avenues for accelerated computational analysis of complex crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamics of nonequilibrium magnons in gapped Heisenberg antiferromagnets

Nonequilibrium dynamics in spin systems is a topic currently under intense investigation as it provides fundamental insights into thermalization, universality, and exotic transport phenomena. While most of the studies have been focused on ideal closed quantum many-body systems such as ultracold atomic quantum gases and one-dimensional spin chains, driven-dissipative Bose gases in steady states away from equilibrium in classical systems also lead to intriguing nonequilibrium physics. Here, in this work, we theoretically investigate out-of-equilibrium dynamics of magnons in a gapped Heisenberg quantum antiferromagnet based on Boltzmann transport theory. We show that, by treating scattering terms beyond the relaxation-time approximation in the Boltzmann transport equation, energy and particle number conservation mandate that nonequilibrium magnons cannot relax to equilibrium, but decay to other nonequilibrium stationary states. The only decay channel for these stationary states back to equilibrium is through the nonconserving interactions (i.e., changing particle number and/or energy within the magnon system) such as boundary or magnon-phonon scattering. At low temperatures, these nonconserving interactions are much slower processes than intrinsic magnon-magnon interaction in a gapped spin system. Using magnon-phonon interaction as a quintessential type of nonconserving interaction, we then propose that nonequilibrium steady states of magnons can be maintained and tailored using periodic driving at frequencies faster than relaxation due to phonon interactions. These findings reveal a class of classical material systems that are suitable platforms to study nonequilibrium statistical physics and macroscopic phenomena such as classical Bose-Einstein condensation of quasiparticles and magnon supercurrents that are relevant for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Efficient Space–Time Reduced Order Model for Linear Dynamical Systems in Python Using Less than 120 Lines of Code

A classical reduced order model (ROM) for dynamical problems typically involves only the spatial reduction of a given problem. Recently, a novel space–time ROM for linear dynamical problems has been developed [Choi et al., Space–tume reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems, Journal of Computational Physics, 2020], which further reduces the problem size by introducing a temporal reduction in addition to a spatial reduction without much loss in accuracy. The authors show an order of a thousand speed-up with a relative error of less than 10−5 for a large-scale Boltzmann transport problem. In this work, we present for the first time the derivation of the space–time least-squares Petrov–Galerkin (LSPG) projection for linear dynamical systems and its corresponding block structures. Utilizing these block structures, we demonstrate the ease of construction of the space–time ROM method with two model problems: 2D diffusion and 2D convection diffusion, with and without a linear source term. For each problem, we demonstrate the entire process of generating the full order model (FOM) data, constructing the space–time ROM, and predicting the reduced-order solutions, all in less than 120 lines of Python code. We compare our LSPG method with the traditional Galerkin method and show that the space–time ROMs can achieve O(10−3) to O(10−4) relative errors for these problems. Depending on parameter–separability, online speed-ups may or may not be achieved. For the FOMs with parameter–separability, the space–time ROMs can achieve O(10) online speed-ups. Finally, we present an error analysis for the space–time LSPG projection and derive an error bound, which shows an improvement compared to traditional spatial Galerkin ROM methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Solution Irregularity Remediation for Spatial Discretization Error Estimation for S N Transport Solutions

The discrete ordinates linear Boltzmann transport equation is typically solved in its spatially discretized form, incurring spatial discretization error. Quantification of this error for purposes such as adaptive mesh refinement or error analysis requires an a posteriori estimator, which utilizes the numerical solution to the spatially discretized equation to compute an estimate. Because the quality of the numerical solution informs the error estimate, irregularities, present in the true solution for any realistic problem configuration, tend to cause the largest deviation in the error estimate vis-a-vis the true error. In this paper, an analytical partial singular characteristic tracking (pSCT) procedure for reducing the estimator’s error is implemented within our novel residual source estimator for a zeroth-order discontinuous Galerkin scheme, at the additional cost of a single inner iteration. Here, a metric-based evaluation of the pSCT scheme versus the standard residual source estimator is performed over the parameter range of a Method of Manufactured Solutions test suite. The pSCT scheme generates near-ideal accuracy in the estimate in problems where the dominant source of the estimator’s error is the solution irregularity, namely, problems where the true solution is discontinuous and problems where the true solution’s first derivative is discontinuous and the scattering ratio is low. In problems where the scattering ratio is high and the true solution is discontinuous in the first derivative, the error in the scattering source, which is not converged by the pSCT scheme, is greater than the error incurred due to the irregularity. Ultimately, a pSCT scheme is judged to be useful for error estimation in problems where the computational cost of the scheme is justified. In the presence of many irregularities, such a scheme may be intractable for general use, but in benchmarks, as an analytical tool, or in problems that have nondissipative discontinuities, the scheme may prove invaluable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗