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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2024

High Field Transport in (Ultra) Wide Bandgap Semiconductors: Diamond Versus Cubic GaN

We provide an analysis of nonlinear transport in diamond and cubic GaN (c-GaN) with emphasis on the different types of optical phonon scattering [i.e., optical deformation potential (ODP) scattering versus polar optic phonon (POP) scattering] limiting the carrier velocity. Both types of carrier mobilities and carrier saturation (peak) velocities in diamond and c-GaN as functions of different doping types and concentrations are obtained by directly solving the Boltzmann equation. Our model indicates that the nonrandomizing nature of POP scattering causes carrier temperature cooling, resulting in higher carrier drift velocity than with ODP scattering. This effect, in addition to the small carrier effective masses and large optical phonon energies, is responsible for the higher peak velocities in c-GaN, compared to carrier drift velocity in diamond.

42 ENGINEERING↗
DOE OSTI2024

Direct and indirect spin current generation and spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers

Spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers are studied using a combination of symmetry analysis, circuit theory, semiclassical simulations, and first-principles calculations using the nonequilibrium Green's function method with supercell disorder averaging. Here, we focus on unconventional processes involving the interplay between the two ferromagnetic layers, which are classified into direct and indirect mechanisms. The direct mechanism involves spin current generation by one ferromagnetic layer and its subsequent absorption by the other. In the indirect mechanism, the in-plane spin-polarized current from one ferromagnetic layer “leaks” into the other layer, where it is converted into an out-of-plane spin current and reabsorbed by the original layer. The direct mechanism results in a predominantly dampinglike torque, which damps the magnetization towards a certain direction 𝐬 𝑑 . The indirect mechanism results in a predominantly fieldlike torque with respect to a generally different direction 𝐬 𝑓 . Similarly to the current-in-plane giant magnetoresistance, the indirect mechanism is only active if the thickness of the nonmagnetic spacer is smaller than or comparable to the mean free path. Numerical calculations for a semiclassical model based on the Boltzmann equation confirm the presence of both direct and indirect mechanisms of spin current generation. First-principles calculations reveal sizable unconventional spin-orbit torques in Co/Cu/Co, Py/Cu/Py, and Co/Pt/Co trilayers and provide strong evidence of indirect spin current generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Thermalization of Interacting Quasi-One-Dimensional Systems

Many experimentally relevant systems are quasi-one-dimensional, consisting of nearly decoupled chains. In these systems, there is a natural separation of scales between the strong intrachain interactions and the weak interchain coupling. When the intrachain interactions are integrable, weak interchain couplings play a crucial part in thermalizing the system. Here, we develop a Boltzmann-equation formalism involving a collision integral that is asymptotically exact for any interacting integrable system, and apply it to develop a quantitative theory of relaxation in coupled Bose gases in the experimentally relevant Newton's cradle setup. Here we find that relaxation involves a broad spectrum of timescales. We provide evidence that the Markov process governing relaxation at late times is gapless; thus, the approach to equilibrium is generally nonexponential, even for spatially uniform perturbations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2023

Theory of the ion–electron temperature relaxation rate in strongly magnetized plasmas

Recent works have shown that strongly magnetized plasmas characterized by having a gyrofrequency greater than the plasma frequency exhibit novel transport properties. One example is that the friction force on a test charge shifts, obtaining components perpendicular to its velocity in addition to the typical stopping power component antiparallel to its velocity. Here, we apply a recent generalization of the Boltzmann equation for strongly magnetized plasmas to calculate the ion–electron temperature relaxation rate. Strong magnetization is generally found to increase the temperature relaxation rate perpendicular to the magnetic field and to cause the temperatures parallel and perpendicular to the magnetic field to not relax at equal rates. This, in turn, causes a temperature anisotropy to develop during the equilibration. Strong magnetization also breaks the symmetry of independence of the sign of the charges of the interacting particles on the collision rate, commonly known as the “Barkas effect.” Furthermore, it is found that the combination of oppositely charged interaction and strong magnetization causes the ion–electron parallel temperature relaxation rate to be significantly suppressed, scaling inversely proportional to the magnetic field strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anderson acceleration with approximate calculations: Applications to scientific computing

Here we provide rigorous theoretical bounds for Anderson acceleration (AA) that allow for approximate calculations when applied to solve linear problems. We show that, when the approximate calculations satisfy the provided error bounds, the convergence of AA is maintained while the computational time could be reduced. We also provide computable heuristic quantities, guided by the theoretical error bounds, which can be used to automate the tuning of accuracy while performing approximate calculations. For linear problems, the use of heuristics to monitor the error introduced by approximate calculations, combined with the check on monotonicity of the residual, ensures the convergence of the numerical scheme within a prescribed residual tolerance. Motivated by the theoretical studies, we propose a reduced variant of AA, which consists in projecting the least-squares used to compute the Anderson mixing onto a subspace of reduced dimension. The dimensionality of this subspace adapts dynamically at each iteration as prescribed by the computable heuristic quantities. We numerically show and assess the performance of AA with approximate calculations on: (i) linear deterministic fixed-point iterations arising from the Richardson's scheme to solve linear systems with open-source benchmark matrices with various preconditioners and (ii) non-linear deterministic fixed-point iterations arising from non-linear time-dependent Boltzmann equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Multi-component Dark Matter and small scale structure formation

We consider the evolution of non-thermal dark matter perturbations in models which contain both Weakly Interacting Massive Particles (WIMPs) and axions. Using constraints from existing observations we examine the percentage of WIMPs and axions that may comprise the cosmological dark matter budget in models with an Early Matter Dominated Epoch (EMDE) — where entropy production is important. After carefully tracking the thermal evolution of the various species by solving the Boltzmann equations, we consider the enhancement of perturbations that may have led to early structure formation for axions and WIMPs. We investigate the impact of enhanced perturbations on the parameter space of both species, after imposing existing constraints from indirect detection experiments. Given these constraints we establish the feasibility of axions to form miniclusters in the early universe in EMDEs for a given percentage of allowed WIMPs. We find that EMDEs with low reheat temperatures near the BBN limit are preferred for axion minicluster formation. When the EMDE is caused by string moduli, the WIMP contribution to the relic density is set by the moduli branching to dark matter at the level of ≲ $\mathcal{O}$ (10 −4 ).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dark matter and dark radiation from the early universe with a modulus coupled to the PQMSSM

The supersymmetrized DFSZ axion model is especially compelling in that it contains 1. the SUSY solution to the gauge hierarchy problem, 2. the Peccei-Quinn (PQ) solution to the strong CP problem and 3. the Kim-Nilles solution to the SUSY μ problem. In a string setting, where a discrete R-symmetry ($Z$$^{R}_{24}$ for example) may emerge from the compactification process, a high-quality accidental axion (accion) can emerge from the accidental, approximate remnant global U(1) PQ symmetry where the decay constant f a is linked to the SUSY breaking scale, and is within the cosmological sweet zone. In this setup, one also expects the presence of stringy remnant moduli fields Φ i . Here, we consider the situation of a single light modulus Φ coupled to the PQMSSM in the early universe, with mixed axion plus higgsino-like WIMP dark matter. We evaluate dark matter and dark radiation production via nine coupled Boltzmann equations and assess the severity of the cosmological moduli problem (CMP) along with dark matter and dark radiation production rates. We find that typically the light modulus mass should be m Φ ≳ 10 4 TeV to avoid the moduli-induced dark matter overproduction problem. If one is able to (anthropically) tune the modulus field amplitude, we find a value of Φ 0 ≲ 10 –7 m P would be required to solve the overall CMP.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Global flow structure and exact formal transseries of the Gubser flow in kinetic theory

In this work we introduce the generic conditions for the existence of a non-equilibrium attractor that is an invariant manifold determined by the long-wavelength modes of the physical system. We investigate the topological properties of the global flow structure of the Gubser flow for the Israel-Stewart theory and a kinetic model for the Boltzmann equation by employing Morse-Smale theory. We present a complete classification of the invariant submanifolds of the flow and determine all the possible flow lines connecting any pair of UV/IR fixed points. The formal transseries solutions to the Gubser dynamical system around the early-time (UV) and late-time (IR) fixed points are constructed and analyzed. It is proven that these solutions are purely perturbative (or power-law asymptotic) series with a finite radius of convergence. Based on these analyses, we find that Gubser-like expanding kinetic systems do not hydrodynamize owing to the failure of the hydrodynamization process which heavily relies on the classification of (non)hydrodynamic modes in the IR regime. This is in contrast to longitudinal boost-invariant plasmas where the asymptotic dynamics is described by a few terms of the hydrodynamic gradient expansion. We finally compare our results for both Bjorken and Gubser conformal kinetic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

FIMP dark matter from flavon portals

We investigate the phenomenology of a non-thermal dark matter (DM) candidate in the context of flavor models that explain the hierarchy in the masses and mixings of quarks and leptons via the Froggatt-Nielsen (FN) mechanism. A flavor-dependent U(1) FN symmetry explains the fermion mass and mixing hierarchy, and also provides a mechanism for suppressed interactions of the DM, assumed to be a Majorana fermion, with the Standard Model (SM) particles, resulting in its FIMP (feebly interacting massive particle) character. Such feeble interactions are mediated by a flavon field through higher dimensional operators governed by the U(1) FN charges. We point out a natural stabilizing mechanism for the DM within this framework with the choice of half-integer U(1) FN charge n for the DM fermion, along with integer charges for the SM fermions and the flavon field. In this flavon portal scenario, the DM is non-thermally produced from the decay of the flavon in the early universe which becomes a relic through the freeze-in mechanism. We explore the allowed parameter space for this DM candidate from relic abundance by solving the relevant Boltzmann equations. We find that reproducing the correct relic density requires the DM mass to be in the range (100 – 300) keV for n = 7.5 and (3 – 10) MeV for n = 8.5 where n is the U(1) FN charge of the DM fermion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Origin of nontopological soliton dark matter: solitosynthesis or phase transition

This work demonstrates that nontopological solitons with large global charges and masses, even above the Planck scale, can form in the early universe and dominate the dark matter abundance. In solitosynthesis, solitons prefer to grow as large as possible under equilibrium dynamics when an initial global charge asymmetry is present. Their abundance is set by when soliton formation via particle fusion freezes out, and their charges are set by the time it takes to accumulate free particles. This work improves the estimation of both quantities, and in particular shows that much larger-charged solitons form than previously thought. The results are estimated analytically and validated numerically by solving the coupled Boltzmann equations. Without solitosynthesis, phase transitions can still form solitons from particles left inside false-vacuum pockets and determine their present-day abundance and properties. Even with zero charge asymmetry, solitons formed in this way can have very large charges on account of statistical fluctuations in the numbers of (anti)particles inside each pocket.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗
DOE OSTI2020

Molecular Statics Analyses of Thermodynamics and Kinetics of Hydrogen Cottrell Atmosphere Formation Around Edge Dislocations in Aluminum

Aluminum alloys are being explored as lightweight structural materials for use in hydrogen-containing environments. In this paper, to understand hydrogen effects on deformation, we perform molecular statics studies of the hydrogen Cottrell atmosphere around edge dislocations in aluminum. First, we calculate the hydrogen binding energies at all interstitial sites in a periodic aluminum crystal containing an edge dislocation dipole. This allows us to use the Boltzmann equation to quantify the hydrogen Cottrell atmosphere. Based on these binding energies, we then construct a continuum model to study the kinetics of the hydrogen Cottrell atmosphere formation. Finally, we compare our results with existing theories and discuss the effects of hydrogen on deformation of aluminum-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Current-conserving relativistic linear response for collisional plasmas

Here we investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi–Dirac equilibrium distribution in terms of the average relaxation rate k. We study the resulting covariant, gauge invariant, and current conserving form of the polarization tensor in the ultrarelativistic and non-relativistic limits. We evaluate the susceptibility in the ultrarelativistic limit and explore their dependence on k. Finally, we study the dispersion relations for the longitudinal and transverse poles of the propagator. We show that for k >2ω p , where ω p is the plasma frequency, the plasma wave modes are overdamped. In the opposite case, k$\ll$ ω p , the propagating plasma modes are weakly damped.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Spin-dependent Seebeck and Nernst effects in an ideal skyrmion gas

In this work, we theoretically and numerically study spin-dependent Seebeck and Nernst effects in 2D ferromagnetic materials with the topological spin texture (skyrmion and vortex) ideal gas. From the numerical solution of the matrix Boltzmann equation for a nonequilibrium distribution function and the Lippmann–Schwinger equation for a T-matrix we find the strong nonlinear behaviors in the thermoelectric coefficients depending on skyrmion/vortex diameters and electron concentrations. In particular, the dramatic dependences in the Seebeck and Nernst coefficients take place at larger magnetic texture sizes where the abrupt sign flip in the vortex Seebeck and Nernst coefficients occurs in the narrow region of electron concentrations. In this case the normalized Nernst coefficient changes from +5 to -7. The spin-dependent thermoelectric coefficients are proportional to T at low temperatures for all skyrmion/vortex sizes.

36 MATERIALS SCIENCE↗