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

Conservative 3 + 1 general relativistic Boltzmann equation

We present a new derivation of the conservative form of the general relativistic Boltzmann equation and specialize it to the 3+1 metric. The resulting transport equation is intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are lab frame coordinate basis spacetime position components and comoving frame curvilinear momentum space coordinates. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we examine the relationship between lepton number and four-momentum exchange between neutrinos and the fluid.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multidisciplinary benchmarks of a conservative spectral solver for the nonlinear Boltzmann equation

The Boltzmann equation describes the evolution of the phase-space probability distribution of classical particles under binary collisions. Approximations to it underlie the basis for several scholarly fields, including aerodynamics and plasma physics. While these approximations are appropriate in their respective domains, they can be violated in niche but diverse applications which require direct numerical solution of the original nonlinear Boltzmann equation. An expanded implementation of the Galerkin–Petrov conservative spectral algorithm is employed to study a wide variety of physical problems. Enabled by distributed precomputation, solutions of the spatially homogeneous Boltzmann equation can be achieved in seconds on modern personal hardware, while spatially-inhomogeneous problems are solvable in minutes. Here, several benchmarks are presented focusing on accuracy compared to both analytic theoretical predictions and other Boltzmann solvers. These benchmarks span several physical domains including weakly ionized plasma, gaseous fluids, and atomic-plasma interaction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Boltzmann equation solver adapted to emergent chemical non-equilibrium

We present a novel method to solve the spatially homogeneous and isotropic relativistic Boltzmann equation. We employ a basis set of orthogonal polynomials dynamically adapted to allow for emergence of chemical non-equilibrium. Two time dependent parameters characterize the set of orthogonal polynomials, the effective temperature T(t) and phase space occupation factor ϒ(t). In this first paper we address (effectively) massless fermions and derive dynamical equations for T(t) and ϒ(t) such that the zeroth order term of the basis alone captures the particle number density and energy density of each particle distribution. We validate our method and illustrate the reduced computational cost and the ability to easily represent final state chemical non-equilibrium by studying a model problem that is motivated by the physics of the neutrino freeze-out processes in the early Universe, where the essential physical characteristics include reheating from another disappearing particle component (e{sup ±}-annihilation)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Nonlinear dynamics from the relativistic Boltzmann equation in the Friedmann-Lemaître-Robertson-Walker spacetime

The dissipative dynamics of an expanding massless gas with constant cross section in a spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) universe is researched. The mathematical problem of solving the full nonlinear relativistic Boltzmann equation is recast into an infinite set of nonlinear ordinary differential equations for the moments of the one-particle distribution function. Momentum-space resolution is determined by the number of nonhydrodynamic modes included in the moment hierarchy, i.e., by the truncation order. We show that in the FLRW spacetime the nonhydrodynamic modes decouple completely from the hydrodynamic degrees of freedom. This results in the system flowing as an ideal fluid while at the same time producing entropy. The solutions to the nonlinear Boltzmann equation exhibit transient tails of the distribution function with nontrivial momentum dependence. The evolution of this tail is not correctly captured by the relaxation time approximation nor by the linearized Boltzmann equation. Yet, the latter probes additional high-momentum details unresolved by the relaxation time approximation. While the expansion of the FLRW spacetime is slow enough for the system to move towards (and not away from) local thermal equilibrium, it is not sufficiently slow for the system to actually ever reach complete local equilibrium. Equilibration is fastest in the relaxation time approximation, followed, in turn, by kinetic evolution with a linearized and a fully nonlinear Boltzmann collision term.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Chapman–Enskog Analyses on the Gray Lattice Boltzmann Equation Method for Fluid Flow in Porous Media

The gray lattice Boltzmann equation (GLBE) method has recently been used to simulate fluid flow in porous media. It employs a partial bounce-back of populations (through a fractional coefficient θ, which represents the fraction of populations being reflected by the solid phase) in the evolution equation to account for the linear drag of the medium. Several particular GLBE schemes have been proposed in the literature and these schemes are very easy to implement; but there exists uncertainty about the need for redefining the macroscopic velocity as there has been no systematic analysis to recover the Brinkman equation from the various GLBE schemes. Rigorous Chapman–Enskog analyses are carried out to show that the momentum equation recovered from these schemes can satisfy Brinkman equation to second order in ϵ only if θ=O(ϵ) in which ϵ is the ratio of the lattice spacing to the characteristic length of physical dimension. The need for redefining macroscopic velocity is shown to be scheme-dependent. When a body force is encountered such as the gravitational force or that caused by a pressure gradient, different forms of forcing redefinitions are required for each GLBE scheme.

42 ENGINEERING↗
DOE OSTI2015

Bianchi type-I magnetized cosmological models for the Einstein-Boltzmann equation with the cosmological constant

Global solutions regular for the Einstein-Boltzmann equation on a magnetized Bianchi type-I cosmological model with the cosmological constant are investigated. We suppose that the metric is locally rotationally symmetric. The Einstein-Boltzmann equation has been already considered by some authors. But, in general Bancel and Choquet-Bruhat [Ann. Henri Poincaré XVIII(3), 263 (1973); Commun. Math. Phys. 33, 83 (1973)], they proved only the local existence, and in the case of the nonrelativistic Boltzmann equation. Mucha [Global existence of solutions of the Einstein-Boltzmann equation in the spatially homogeneous case. Evolution equation, existence, regularity and singularities (Banach Center Publications, Institute of Mathematics, Polish Academy of Science, 2000), Vol. 52] obtained a global existence result, for the relativistic Boltzmann equation coupled with the Einstein equations and using the Yosida operator, but confusing unfortunately with the nonrelativistic case. Noutchegueme and Dongho [Classical Quantum Gravity 23, 2979 (2006)] and Noutchegueme, Dongho, and Takou [Gen. Relativ. Gravitation 37, 2047 (2005)], have obtained a global solution in time, but still using the Yosida operator and considering only the uncharged case. Noutchegueme and Ayissi [Adv. Stud. Theor. Phys. 4, 855 (2010)] also proved a global existence of solutions to the Maxwell-Boltzmann system using the characteristic method. In this paper, we obtain using a method totally different from those used in the works of Noutchegueme and Dongho [Classical Quantum Gravity 23, 2979 (2006)], Noutchegueme, Dongho, and Takou [Gen. Relativ. Gravitation 37, 2047 (2005)], Noutchegueme and Ayissi [Adv. Stud. Theor. Phys. 4, 855 (2010)], and Mucha [Global existence of solutions of the Einstein-Boltzmann equation in the spatially homogeneous case. Evolution equation, existence, regularity and singularities (Banach Center Publications, Institute of Mathematics, Polish Academy of Science, 2000), Vol. 52] the global in time existence and uniqueness of a regular solution to the Einstein-Maxwell-Boltzmann system with the cosmological constant. We define and we use the weighted Sobolev separable spaces for the Boltzmann equation; some special spaces for the Einstein equations, then we clearly display all the proofs leading to the global existence theorems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Properties of the Boltzmann equation in the classical approximation

We examine the Boltzmann equation with elastic point-like scalar interactions in two different versions of the the classical approximation. Although solving numerically the Boltzmann equation with the unapproximated collision term poses no problem, this allows one to study the effect of the ultraviolet cutoff in these approximations. This cutoff dependence in the classical approximations of the Boltzmann equation is closely related to the non-renormalizability of the classical statistical approximation of the underlying quantum field theory. The kinetic theory setup that we consider here allows one to study in a much simpler way the dependence on the ultraviolet cutoff, since one has also access to the non-approximated result for comparison.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

CMB spectral distortions as solutions to the Boltzmann equations

We propose to re-interpret the cosmic microwave background spectral distortions as solutions to the Boltzmann equation. This approach makes it possible to solve the second order Boltzmann equation explicitly, with the spectral y distortion and the momentum independent second order temperature perturbation, while generation of μ distortion cannot be explained even at second order in this framework. We also extend our method to higher order Boltzmann equations systematically and find new type spectral distortions, assuming that the collision term is linear in the photon distribution functions, namely, in the Thomson scattering limit. As an example, we concretely construct solutions to the cubic order Boltzmann equation and show that the equations are closed with additional three parameters composed of a cubic order temperature perturbation and two cubic order spectral distortions. The linear Sunyaev-Zel'dovich effect whose momentum dependence is different from the usual y distortion is also discussed in the presence of the next leading order Kompaneets terms, and we show that higher order spectral distortions are also generated as a result of the diffusion process in a framework of higher order Boltzmann equations. The method may be applicable to a wider class of problems and has potential to give a general prescription to non-equilibrium physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A bi-fidelity method for the multiscale Boltzmann equation with random parameters

In this paper, we study the multiscale Boltzmann equation with multi-dimensional random parameters by a bi-fidelity stochastic collocation (SC) method developed in [52], [70], [71]. By choosing the compressible Euler system as the low-fidelity model, we adapt the bi-fidelity SC method to combine computational efficiency of the low-fidelity model with high accuracy of the high-fidelity (Boltzmann) model. With only a small number of high-fidelity asymptotic-preserving solver runs for the Boltzmann equation, the bi-fidelity approximation can capture well the macroscopic quantities of the solution to the Boltzmann equation in the random space. Additionally, a priori estimate on the accuracy between the high- and bi-fidelity solutions together with a convergence analysis is established. Finally, we present extensive numerical experiments to verify the efficiency and accuracy of our proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Transient oscillations in a macroscopic effective theory of the Boltzmann equation

A new transient effective theory of the relativistic Boltzmann equation is derived for locally momentum-anisotropic systems. In the expansion of the distribution function around a local “quasi-equilibrium” state, a nonhydrodynamic dynamical degree of freedom is introduced at leading order that breaks local momentum isotropy. By replacing the deviation of the distribution function from this quasi-equilibrium state in terms of moments of the leading-order distribution and applying a systematic power-counting scheme that orders the nonhydrodynamic modes by their microscopic time scales, a closed set of equations for the dynamical degrees of freedom is obtained. Truncating this set at the level of the slowest nonhydroynamic mode, we find that it exhibits transient oscillatory behavior—a phenomenon previously found only in strongly coupled theories, where it appears to be generic. In weakly coupled systems described by the Boltzmann equation, these transient oscillations depend on the breaking of local momentum isotropy being treated nonperturbatively at leading order in the expansion of the distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Classical non-Markovian Boltzmann equation

The modeling of particle transport involves anomalous diffusion, (x²(t) ) ∝ t{sup α} with α ≠ 1, with subdiffusive transport corresponding to 0 < α < 1 and superdiffusive transport to α > 1. These anomalies give rise to fractional advection-dispersion equations with memory in space and time. The usual Boltzmann equation, with only isolated binary collisions, is Markovian and, in particular, the contributions of the three-particle distribution function are neglected. We show that the inclusion of higher-order distribution functions give rise to an exact, non-Markovian Boltzmann equation with resulting transport equations for mass, momentum, and kinetic energy with memory in both time and space. The two- and the three-particle distribution functions are considered under the assumption that the two- and the three-particle correlation functions are translationally invariant that allows us to obtain advection-dispersion equations for modeling transport in terms of spatial and temporal fractional derivatives.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Fixed points and flow analysis on off-equilibrium dynamics in the boson Boltzmann equation

We consider fixed points of steady solutions and flow directions using the boson Boltzmann equation that is a one-dimensionally reduced kinetic equation after the angular integration. With an elastic collision integral of the two-to-two scattering process, in the dense (dilute) regime where the distribution function is large (small), the boson Boltzmann equation has approximate fixed points with a power-law spectrum in addition to the thermal distribution function. We argue that the power-law fixed point can be exact in special cases. We elaborate a graphical presentation to display evolving flow directions similarly to the renormalization group flow, which explicitly exhibits how fixed points are connected and parameter space is separated by critical lines. We discuss that such a flow diagram contains useful information on thermalization processes out of equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Analysis of the multicomponent relativistic Boltzmann equation for electron scattering in big bang nucleosynthesis

Big-bang nucleosynthesis (BBN) is valuable as a means to constrain the physics of the early universe and it is the only probe of the radiation-dominated epoch. A fundamental assumption in BBN is that the nuclear velocity distributions obey Maxwell-Boltzmann (MB) statistics as they do in stars. Specifically, the BBN epoch is characterized by a dilute baryon plasma for which the velocity distribution of nuclei is mainly determined by the dominant Coulomb elastic scattering with mildly relativistic electrons. One must therefore deduce the momentum distribution for reacting nuclei from the multicomponent relativistic Boltzmann equation. However, the full multicomponent relativistic Boltzmann equation has only recently been analyzed and its solution has only been worked out in special cases. Moreover, a variety of schemes have been proposed that introduce nonthermal components into the BBN environment which can alter the thermal distribution of reacting nuclei. Here, we construct the relativistic Boltzmann equation in the context of BBN. We also derive a Langevin model and perform relativistic Monte-Carlo simulations which clarify the baryon distribution during BBN and can be used to analyze any relaxation from a nonthermal injection. Further, we show by these analyses that the thermalization process leads to a nuclear distribution function that remains very close to MB statistics even during the most relativistic environment relevant to BBN. Hence, the predictions of standard BBN remain unchanged.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Solution of the nonlinear Poisson-Boltzmann equation: Application to ionic diffusion in cementitious materials

A robust numerical solution of the nonlinear Poisson-Boltzmann equation for asymmetric polyelectrolyte solutions in discrete pore geometries is presented. Comparisons to the linearized approximation of the Poisson-Boltzmann equation reveal that the assumptions leading to linearization may not be appropriate for the electrochemical regime in many cementitious materials. Implications of the electric double layer on both partitioning of species and on diffusive release are discussed. The influence of the electric double layer on anion diffusion relative to cation diffusion is examined.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The halo Boltzmann equation

Dark matter halos are the building blocks of the universe as they host galaxies and clusters. The knowledge of the clustering properties of halos is therefore essential for the understanding of the galaxy statistical properties. We derive an effective halo Boltzmann equation which can be used to describe the halo clustering statistics. In particular, we show how the halo Boltzmann equation encodes a statistically biased gravitational force which generates a bias in the peculiar velocities of virialized halos with respect to the underlying dark matter, as recently observed in N-body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Perfect fluidity of a dissipative system: Analytical solution for the Boltzmann equation in AdS 2 Ⓧ S 2

In this paper we obtain an analytical solution of the relativistic Boltzmann equation under the relaxation time approximation that describes the out-of-equilibrium dynamics of a radially expanding massless gas. This solution is found by mapping this expanding system in flat spacetime to a static flow in the curved spacetime AdS 2 Ⓧ S 2 . We further derive explicit analytic expressions for the momentum dependence of the single-particle distribution function as well as for the spatial dependence of its moments. We find that this dissipative system has the ability to flow as a perfect fluid even though its entropy density does not match the equilibrium form. The nonequilibrium contribution to the entropy density is shown to be due to higher-order scalar moments (which possess no hydrodynamical interpretation) of the Boltzmann equation that can remain out of equilibrium but do not couple to the energy-momentum tensor of the system. Furthermore, in this system the slowly moving hydrodynamic degrees of freedom can exhibit true perfect fluidity while being totally decoupled from the fast moving, nonhydrodynamical microscopic degrees of freedom that lead to entropy production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗