Direct molecular gas dynamics simulations of re-entry vehicles via the Boltzmann equation
Explore the source record for details and available documents.
Engineering - Aeronautics, Mathematics and Computi↗
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.
Pure quantum extension of the semiclassical Boltzmann-Uehling-Uhlenbeck equation
The Boltzmann equation is the traditional framework in which one extends the time-dependent mean field classical description of a many-body system to include the effect of particle-particle collisions in an approximate manner. A semiclassical extension of this approach to quantum many-body systems was suggested by Uehling and Uhlenbeck in 1933 for both Fermi and Bose statistics, and many further developments of this approach are known as the Boltzmann-Uehling-Uhlenbeck (BUU) equations. In this paper I introduce a pure quantum version of the BUU type of equations, which is mathematically equivalent to a generalized time-dependent density functional theory extended to superfluid systems. As expected, during nonequilibrium processes the quantum Boltzmann one-body entropy increases during evolution.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Fusion and fission events regulate endosome maturation and viral escape
Endosomes are intracellular vesicles that mediate the communication of the cell with its extracellular environment. They are an essential part of the cell’s machinery regulating intracellular trafficking via the endocytic pathway. Many viruses, which in order to replicate require a host cell, attach themselves to the cellular membrane; an event which usually initiates uptake of a viral particle through the endocytic pathway. In this way viruses hijack endosomes for their journey towards intracellular sites of replication and avoid degradation without host detection by escaping the endosomal compartment. Recent experimental techniques have defined the role of endosomal maturation in the ability of enveloped viruses to release their genetic material into the cytoplasm. Endosome maturation depends on a family of small hydrolase enzymes (or GTPases) called Rab proteins, arranged on the cytoplasmic surface of its membrane. Here, we model endosomes as intracellular compartments described by two variables (its levels of active Rab5 and Rab7 proteins) and which can undergo coagulation (or fusion) and fragmentation (or fission). The key element in our approach is the “per-cell endosomal distribution” and its dynamical (Boltzmann) equation. The Boltzmann equation allows us to derive the dynamics of the total number of endosomes in a cell, as well as the mean and the standard deviation of its active Rab5 and Rab7 levels. We compare our mathematical results with experiments of Dengue viral escape from endosomes. The relationship between endosomal active Rab levels and pH suggests a mechanism that can account for the observed variability in viral escape times, which in turn regulate the viability of a viral intracellular infection.
59 BASIC BIOLOGICAL SCIENCES↗
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↗
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↗
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↗
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↗
Scaling and adiabaticity in a rapidly expanding gluon plasma
In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Post-Newtonian kinetic theory
A kinetic theory for relativistic gases in the presence of gravitational fields is developed in the second post-Newtonian approximation. The corresponding Boltzmann equation is determined from the evolution of the one-particle distribution function with respect to the proper time along the world line of the particle. From the knowledge of the equilibrium Maxwell–Jüttner distribution function in the second post-Newtonian approximation the components of the particle four-flow and energy–momentum tensor are obtained. The Eulerian hydrodynamic equations for the mass density, mass-energy density and momentum density in the second post-Newtonian approximation are determined from the Boltzmann equation. It is shown that the combination of the hydrodynamic equations of mass and mass-energy densities leads to the hydrodynamic equation for the internal energy density in the first post-Newtonian approximation.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation
We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Distribution function of nuclei from $\mathcal{e}^±$ scattering in the presence of a strong primordial magnetic field
The amplitude of the primordial magnetic field (PMF) is constrained from observational limits on primordial nuclear abundances. Within this constraint, it is possible that nuclear motion is regulated by Coulomb scattering with electrons and positrons (e ± s), while e ± s are affected by a PMF rather than collisions. For example, at a temperature of 10 9 K, thermal nuclei typically experience ~ 10 21 scatterings per second that are dominated by very small angle scattering leading to minuscule changes in the nuclear kinetic energy of order $\mathcal{O}$(1) eV. In this paper the upper limit on the effects of a possible discretization of the e ± momenta by the PMF on the nuclear momentum distribution is estimated under the extreme assumptions that the momentum of the e ± is relaxed before and after Coulomb scattering to Landau levels, and that during Coulomb scattering the PMF is neglected. This assumption explicitly breaks the time reversal invariance of Coulomb scattering, and the Maxwell-Boltzmann distribution is not a trivial steady state solution of the Boltzmann equation under these assumptions. Here we numerically evaluate the collision terms in the Boltzmann equation, and show that the introduction of a special direction in the e ± distribution by the PMF generates no directional dependence of the collisional destruction term of nuclei. Large anisotropies in the nuclear distribution function are then constrained from big bang nucleosynthesis. Ultimately, we conclude that a PMF does not significantly affect the isotropy or big bang nucleosynthesis.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Semi-Lagrangian nodal discontinuous Galerkin method for the BGK model
In this work, we propose a semi-Lagrangian (SL) nodal discontinuous Galerkin (DG) solver for the BGK equation. The BGK model was introduced by Bhatnagar, Gross, and Krook [1] as a relaxation model for the fundamental Boltzmann equation [5], which describes the kinetic dynamic of rarefied gases with a probability distribution function. The challenges of designing efficient numerical schemes for the Boltzmann equation mainly come from its high dimensionality and complicated nonlinear collision operator. The BGK model gains interests since it has much lower computational cost, due to the relatively simple structure of the relaxation operator in replacement of the collision operator, while simultaneously preserving several important physical properties, such as macroscopic quantities and dissipation of entropy.
97 MATHEMATICS AND COMPUTING↗
ULYSSES: Universal LeptogeneSiS Equation Solver
ULYSSES is a python package that calculates the baryon asymmetry produced from leptogenesis in the context of a type-I seesaw mechanism. The code solves the semi-classical Boltzmann equations for points in the model parameter space as specified by the user. We provide a selection of predefined Boltzmann equations as well as a plugin mechanism for externally provided models of leptogenesis. Furthermore, the ULYSSES code provides tools for multi-dimensional parameter space exploration. The emphasis of the code is on user flexibility and rapid evaluation. It is publicly available at https://github.com/earlyuniverse/ulysses
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A quick introduction to linear transport phenomena
The goal of this report is to provide an introductory-level overview of the linear Boltzmann equation in the context of neutron transport. After deriving the transport equation, we discuss some of its basic applications in reactor theory. Finally, we review several simplifying approximations of the linearized Boltzmann equation that are essential to solving it in many applications. Although the context of neutron transport is called upon to add concreteness to our discussion, many of the concepts and approximations that we discuss remain relevant in other many other phenomena.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Hydrodynamic generators in relativistic kinetic theory
Here, we resum the nonequilibrium gradient corrections to a single-particle distribution function evolved by the Boltzmann equation in the relaxation time approximation (RTA). We first study a system undergoing Bjorken expansion and show that, for a constant relaxation time, the exact solution of the RTA Boltzmann equation at late times (i.e., after the decay of nonhydrodynamic modes) generates the Borel resummed Chapman-Enskog series. Extending this correspondence to systems without Bjorken symmetry, we construct a ( 3 + 1 )-dimensional hydrodynamic generator for RTA kinetic theory, which is an integral representation of the Chapman-Enskog series in the limit of vanishing nonhydrodynamic modes. Relaxing this limit, we discover at earlier times a set of nonhydrodynamic modes coupled to the RTA Chapman-Enskog expansion. Including the dynamics of these nonhydrodynamic modes is shown to control the emergence of hydrodynamics as an effective field theory description of nonequilibrium fluids, which works well even for far-off-equilibrium situations where the Knudsen number is initially large.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
A Plasma Modeling Hierarchy and Verification Approach
This report reviews a hierarchy of formal mathematical models for describing plasma phenomena. Starting with the Boltzmann equation, a sequence of approximations and modeling assumptions can be made that progressively reduce to the equations for magnetohydrodynamics. Understanding the assumptions behind each of these models and their mathematical form is essential to appropriate use of each level of the hierarchy. A sequence of moment models of the Boltzmann equation are presented, then focused into a generalized three-fluid model for neutral species, electrons, and ions. This model is then further reduced to a two-fluid model, for which Braginskii described a useful closure. Further reduction of the two-fluid model yields a Generalized Ohm's Law model, which provides a connection to magnetohydrodynamic approaches. A verification approach based on linear plasma waves is presented alongside the model hierarchy, which is intended as an initial and necessary but not sufficient step for verification of plasma models within this hierarchy.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Low-temperature T -linear resistivity due to umklapp scattering from a critical mode
In this study, we consider the transport properties of a model of fermions scattered by a critical bosonic mode. The mode is overdamped and scattering is mainly in the forward direction. Such a mode appears at the quantum critical point for an electronic nematic phase transition and in gauge theories for a U(1) spin liquid. It leads to a short fermion lifetime, violating Landau's criterion for a Fermi liquid. In spite of this, transport can be described by a Boltzmann equation. We include momentum relaxation by umklapp scattering, supplemented by weak impurity scattering. We find that above a very low temperature which scales with Δ q 3 , where Δ q is the minimum umklapp scattering vector, the resistivity is linear in T with a coefficient which is independent of the amount of disorder. This picture holds until the temperature dependent part of the resistivity exceeds that due to impurity scattering. We compare the relaxation time approximation with an exact numerical solution of the Boltzmann equation. Surprisingly we find that, unlike the resistivity, the Hall coefficient strongly deviates from the relaxation time approximation and shows a strong reduction with increasing temperature. We comment on possible comparisons with experiments on high T c cuprates.
36 MATERIALS SCIENCE↗
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↗