Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOLTZMANN EQUATION”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

High-order diamond differencing schemes for the Boltzmann Fokker-Planck equation in 3D Cartesian geometries

The Boltzmann Fokker-Planck, an approximate form of the linear Boltzmann equation is commonly used to treat efficiently the transport of charged particles in matter. This paper introduces the application of high-order diamond differencing schemes (HODD), specifically the DD1 and DD2 schemes which are 4- and 6-order accurate respectively, to handle the spatial discretization of that equation in 3D Cartesian geometries. The energy deposition solutions for the coupled transport of electrons and photons presented in this work shows that HODD, compared to classical DD scheme, provides correction to the oscillations and a reduced propensity to yield negative fluxes. They are useful tools to minimize local error, notably in regions with abrupt variations of the flux solution. They also can be used to reduced execution time by decreasing the needed number of voxels to obtain a fixed accuracy. On the tested benchmarks, the DD1 scheme is 87%- 92%-91% more accurate than the classical DD scheme for total, mean per-voxels and maximum deviation of energy deposition values respectively. For comparison, a calculation with 8 times more voxels, requiring roughly 2.5 times more time to execute, is 92%-90%-77% more accurate. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Revisiting shear stress tensor evolution: Nonresistive magnetohydrodynamics with momentum-dependent relaxation time

This study aims to develop second-order relativistic viscous magnetohydrodynamics (MHD) derived from kinetic theory within an extended relaxation time approximation (momentum/energy dependent) for the collision kernel. The investigation involves a detailed examination of shear stress tensor evolution equations and associated transport coefficients. The Boltzmann equation is solved using a Chapman-Enskog-like gradient expansion for a charge-conserved conformal system, incorporating a momentum-dependent relaxation time. The derived relativistic nonresistive, viscous second-order MHD equations for the shear stress tensor reveal significant modifications in the coupling with dissipative charge current and magnetic field due to the momentum dependence of the relaxation time. By utilizing a power law parametrization to quantify the momentum dependence of the relaxation time, the anisotropic magnetic field-dependent shear coefficients in the Navier-Stokes limit have been investigated. The resulting viscous coefficients are seen to be sensitive to the momentum dependence of the relaxation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Semi-classical Kinetic Theory for Massive Spin-half Fermions with Leading-order Spin Effects

We consider the quantum kinetic-theory description for interacting massive spin-half fermions using the Wigner function formalism. We derive a general kinetic theory description assuming that the spin effects appear at the classical and quantum level. To track the effect of such different contributions we use the semi-classical expansion method to obtain the generalized dynamical equations including spin, analogous to classical Boltzmann equation. This approach can be used to obtain a collision kernel involving local as well as non-local collisions among the microscopic constituent of the system and eventually, a framework of spin hydrodynamics ensuring the conservation of the energy-momentum tensor and total angular momentum tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

MARBLES (Multi-scale Adaptively Refined Boltzmann LatticE Solver) [SWR-23-37]

MARBLES (Multi-scale Adaptively Refined Boltzmann LatticE Solver) is an open-source computational fluid dynamics package powered by the lattice Boltzmann equations and built on AMReX. In the lattice Boltzmann method, local collisions between meso-scale fictitious particles drive the governing equations which enables MARBLES to easily simulate flow around complex and/or moving geometry without the generation of a body-conforming mesh. Using AMReX data structures and operations ensures a high level of computational performance and parallel scaling on heterogenous architectures while also naturally supporting locally enhanced grid resolution and fidelity through automatic mesh refinement. New domains and problem definitions are easily specified through an input file with examples and guidance on all options and variables provided in the MARBLES documentation.

Henry de Frahan, Marc↗
DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Spectroscopic Determination of Stellar Parameters and Oxygen Abundances for Hyades/Field G–K Dwarfs

It has been suggested that Fe abundances of K dwarfs derived from Fe i and Fe ii lines show considerable discrepancies, and oxygen abundances determined from high-excitation O i 7771–5 triplet lines are appreciably overestimated (the problem becoming more serious toward lower T {sub eff}), which, however, has not yet been widely confirmed. With the aim of clarifying this issue, we spectroscopically determined the atmospheric parameters of 148 G–K dwarfs (Hyades cluster stars and field stars) by assuming the classical Fe i/Fe ii ionization equilibrium as usual, and determined their oxygen abundances by applying the non-local thermal equilibrium spectrum fitting analysis to O i 7771–5 lines. It turned out that the resulting parameters did not show any significant inconsistency with those determined by other methods (for example, the mean differences in T {sub eff} and logg from the well-determined solutions of Hyades dwarfs are mostly ≲100 K and ≲0.1 dex). Likewise, the oxygen abundances of Hyades stars are around [O/H] ∼ +0.2 dex (consistent with the metallicity of this cluster) without exhibiting any systematic T {sub eff}-dependence. Accordingly, we conclude that parameters can be spectroscopically evaluated to a sufficient precision in the conventional manner (based on the Saha–Boltzmann equation for Fe i/Fe ii) and oxygen abundances can be reliably determined from the O i 7771–5 triplet for K dwarfs as far as stars of T {sub eff} ≳ 4500 K are concerned. We suspect that previously reported strongly T {sub eff}-dependent discrepancies may have stemmed mainly from overestimation of weak-line strengths and/or improper T {sub eff} scale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Quantum kinetics of anomalous and nonlinear Hall effects in topological semimetals

Highlights: • Linear and photogalvanic anomalous Hall responses are systematically derived. • Extrinsic mechanisms of AHE include Gaussian, diffractive, hybrid skew scattering. • Diagrammatic calculations are matched to semiclassical picture of AHE. • The Pancharatnam phase of multifold fermions determines the skew scattering amplitude. • Photon-induced interband scattering are accompanied by coordinate shifts. We present a systematic microscopic derivation of the semiclassical Boltzmann equation for band structures with the finite Berry curvature based on Keldysh technique of nonequilibrium systems. In the analysis, an AC electrical driving field is kept up to quadratic order, and both cases of small and large frequencies corresponding to intra- and interband transitions are considered. In particular, this formulation is suitable for the study of nonlinear Hall effect and photogalvanic phenomena. The role of impurity scattering is carefully addressed. Specifically, in addition to previously studied side-jump and skew-scattering processes, quantum interference diffractive contributions are now explicitly incorporated within the developed framework. This theory is applied to multifold fermions in topological semimetals, for which the generic formula for the skew scattering rate from the Pancharatnam phase is obtained along with the corresponding anomalous Hall conductivity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nonlinear Poisson–Boltzmann solutions for charged parallel plates: When opposite charges repel

I present an exact solution of the Poisson–Boltzmann equation for two parallel plates and discuss the solution properties. I discuss in more detail plates with opposite charges: In this case, there are two critical separations, L c,1 < L c,2 . For separations less than L c,1 , the force between plates is repulsive. It switches to attractive at L c,1 , but with the electric potential having the same sign on both plates. For L > L c,2 , the force remains attractive, and the potential at the plates has the same sign as the charge on each plate. I also describe charge regulation, determined by pK a , and provide formulas for both the critical distance where oppositely charged plates repel and their charging process. Finally, the implications of these results for the nanoparticle assembly, as driven by electrostatic interactions, are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2026

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2024

A collision-based hybrid method for the BGK equation

In this article, we apply the collision-based hybrid method introduced by Hauck and McClarren to the Boltzmann equation with the BGK operator and a hyperbolic scaling. An implicit treatment of the source term is used to handle stiffness associated with the BGK operator. Although it helps the numerical scheme become stable with a large time step size, it is still not obvious to achieve the desired order of accuracy due to the relationship between the size of the spatial cell and the mean free path. Without asymptotic preserving property, a very restricted grid size is required to resolve the mean free path, which is not practical. Our approaches are based on the noncollision-collision decomposition of the BGK equation. We introduce the arbitrary order of nodal discontinuous Galerkin (DG) discretization in space with a semi-implicit time-stepping method; we employ the backward Euler time integration for the uncollided equation and the 2nd order predictor-corrector scheme for the collided equation, i.e., both source terms in uncollided and collided equations are treated implicitly and only streaming term in the collided equation is solved explicitly. This improves the computational efficiency without the complexity of the numerical implementation. Numerical results are presented for various Knudsen numbers to present the effectiveness and accuracy of our hybrid method. Also, we compare the solutions of the hybrid and non-hybrid schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Finite scale theory: Predicting nature’s shocks

In this paper, we describe a continuum model that accurately reproduces the experimentally measured structure of physical shocks in a perfect gas. In this work, we begin by presenting a history of shock structure research, theoretical, experimental and numerical, to quantify the significant discrepancies between Navier–Stokes predictions and laboratory measurements. In our first main result, we discuss modifications that generalize the Chapman–Enskog approximation and lead to a new continuum model termed the Finite Scale Equations. In our second main result, we use this continuum model to calculate shock structure with results comparable in accuracy to numerical simulations based on the Boltzmann equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Evolution of Bose–Einstein condensate systems beyond the Gross–Pitaevskii equation

While many phenomena in cold atoms and other Bose–Einstein condensate (BEC) systems are often described using the mean-field approaches, understanding the kinetics of BECs requires the inclusion of particle scattering via the collision integral of the quantum Boltzmann equation. A rigorous approach for many problems in the dynamics of the BEC, such as the nucleation of the condensate or the decay of the persistent current, requires, in the presence of factors making a symmetry breaking possible, considering collisions with thermal atoms via the collision integral. These collisions permit the emergence of vorticity or other signatures of long-range order in the nucleation of the BEC or the transfer of angular momentum to thermal atoms in the decay of persistent current, due to corresponding terms in system Hamiltonians. Here, we also discuss the kinetics of spin–orbit-coupled BEC. The kinetic equation for the particle spin density matrix is derived. Numerical simulations demonstrate significant effects of the collision integral on the dynamics of the spin–orbit-coupled BEC upon quenching of the Raman coupling that generates synthetic electric and magnetic fields.

Bose–Einstein condensation↗
DOE OSTI2022

Nonextensive hydrodynamics of boost-invariant plasmas

We use quasiparticle anisotropic hydrodynamics to study the non-conformal and non-extensive dynamics of a system undergoing boost-invariant Bjorken expansion. To introduce nonextensivity, we use an underlying Tsallis distribution with a time-dependent nonextensivity parameter q. By taking moments of the quasiparticle Boltzmann equation in the relaxation-time approximation, we obtain dynamical equations which allow us to determine the time evolution of all microscopic parameters including q. We compare numerical solutions for bulk observables obtained using the nonextensive evolution with results obtained using quasiparticle anisotropic hydrodynamics with a Boltzmann distribution function (q → 1). We show that the evolution of the temperature, pressure ratio, and scaled energy density, are quite insensitive to which distribution function is assumed. However, we find significant differences in the early-time evolution of the bulk pressure which are observed for even small deviations from the Boltzmann distribution function. Finally, we discuss the existence of non-conformal hydrodynamic attractors for the longitudinal and transverse pressures, the bulk and shear viscous corrections, and the nonextensivity parameter q.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗