Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boltzmann Transport 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Coupled confined phase behavior and transport of methane in slit nanopores

The proximity of the order of magnitude molecular mean free path and pore sizes in nanopores leads to remarkable interactions between molecules and walls. In such systems, the thermodynamic property and the transport behavior of fluids deviate from those at bulk conditions. Molecular dynamics (MD) simulation may be used to investigate the effects of confinement on fluid physics in nanosize pores. However, translating these subscale observations into a larger scale of interest remains a challenging task. In this work, we propose a modified extension of Peng-Robinson equation of state (PR-EOS) motivated and guided by MD simulation results. The shift of critical properties, i.e., pressure and temperature, are evaluated independently. A temperature dependent parameter is introduced to account for the existence of capillary pressure in the two-phase region. Here, this formulation is capable of capturing the shift in critical properties as well as density phase diagrams under various confinement scenarios. We incorporate the proposed EOS in the lattice Boltzmann method (LBM) to study coupled confined phase behavior and transport of methane in nano-size slit pores. Adsorptive strength is determined such that the resulting density ratio matches that from MD simulation. Results indicate that the use of the proposed EOS leads to smaller fluid velocities compared to PR-EOS. This effect is due to stronger interactions between fluid particles under confinement. Also, results show that transport characteristics are impacted by the pressure of the system: in systems with a relatively low pressure, transport seems to be dominated by Knudsen diffusion; however, at a higher pressure, the contribution from viscous flow increases as the pore widens while the influence of Knudsen and surface diffusion diminishes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

An Improved Elastic and Nonelastic Neutron Transport Algorithm for Space Radiation

A neutron transport algorithm including both elastic and nonelastic particle interaction processes for use in space radiation protection for arbitrary shield material is developed. The algorithm is based upon a multiple energy grouping and analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. The algorithm is then coupled to the Langley HZETRN code through a bidirectional neutron evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for an aluminum water shield-target configuration is then compared with MCNPX and LAHET Monte Carlo calculations for the same shield-target configuration. With the Monte Carlo calculation as a benchmark, the algorithm developed in this paper showed a great improvement in results over the unmodified HZETRN solution. In addition, a high-energy bidirectional neutron source based on a formula by Ranft showed even further improvement of the fluence results over previous results near the front of the water target where diffusion out the front surface is important. Effects of improved interaction cross sections are modest compared with the addition of the high-energy bidirectional source terms.

Clowdsley, Martha S.↗
NASA NTRS1975

Theory of photoelectron thermalization and transport in the ionosphere

A new theoretical approach for calculating the equilibrium photoelectron flux energy and pitch angle distribution in the ionosphere is presented. Photoelectron transport, secondary electron production, and energy degradation by the excitation of the discrete energy states of the neutral atmospheric constituents and by continuous energy transfer to the ambient thermal electron gas are included in a manner consistent with the Boltzmann equation which constitutes the foundation of the theory. A difference equation, suitable for numerical solution, is given, and a numerical method for the solution of this equation is discussed in detail.

Mantas, G. P.↗
NASA NTRS1973

Composition dependence of ion transport coefficients in gas mixtures

A simple momentum-transfer theory for the composition dependence of ion mobilities and diffusion coefficients in gas mixtures at arbitrary field strengths is corrected, extended, and compared with a similar theory based on momentum and energy transfer, and with results based on direct solution of the Boltzmann equation by Kihara's method. Final equations are recommended for predicting composition dependences, given only results on ion mobilities and diffusion coefficients in the pure component gases.

Whealton, J. H.↗
NASA NTRS1974

Composition dependence of ion-transport coefficients in gas mixtures

Simple momentum-transfer theory for the composition dependence of ion mobilities and diffusion coefficients in gas mixtures at arbitrary field strengths is corrected and extended, and compared with a similar theory based on momentum and energy transfer, and with results based on direct solution of the Boltzmann equation by Kihara's method. Final equations are recommended for predicting composition dependences, given only results on ion mobilities and diffusion coefficients in the pure component gases.

Whealton, J. H.↗
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 OSTI2022

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

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

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Potential quantum advantage for simulation of fluid dynamics

Numerical simulation of turbulent fluid dynamics needs to either parametrize turbulence—which introduces large uncertainties—or explicitly resolve the smallest scales—which is prohibitively expensive. Here, we provide evidence through analytic bounds and numerical studies that a potential quantum speedup can be achieved to simulate fluid dynamics using quantum computing. Specifically, we provide a lattice Boltzmann formulation of fluid dynamics for which we give evidence that low-order Carleman linearization is much more accurate than previously believed for these systems. This is achieved via a combination of reformulating the Navier-Stokes nonlinearity (u·$\triangledown$u) to lattice-Boltzmann nonlinearity (u 2 ) and accurately linearizing the dynamical equations, which effectively trades nonlinearity for additional degrees of freedom that add negligible expense in the quantum solver. Based on this, we apply a quantum algorithm for simulating the Carleman-linearized lattice Boltzmann equation and provide evidence that its cost scales logarithmically with system size compared with polynomial scaling in the best known classical algorithms. In this paper, we suggest that a quantum advantage may exist for simulating fluid dynamics, paving the way for simulating nonlinear multiscale transport phenomena in a wide range of disciplines using quantum computing.

42 ENGINEERING↗
NASA NTRS1979

Calculation of radiative properties of nonequilibrium hydrogen plasma

A computer program called NEQRAP is described that calculates the radiative properties of nonequilibrium ionized hydrogen. From the given electron temperature, electron density, and atom density values (which do not necessarily satisfy the equilibrium relationship) and intensities of incident radiation, the non-Boltzmann populations of electronic states are computed by solving the equation of quasi-steady-state population distribution. Emission and absorption coefficients are determined as functions of wavelength by invoking the principle of detailed balance between the upper and lower states of each radiative transition. Radiative transport through the medium is computed assuming a one-dimensional uniform slab. The rate of ionic reaction is also computed. When used on a sample case, the program shows that there is a large difference between the calculated intensities of radiation emitted by a bulk of equilibrium and nonequilibrium hydrogen. The accuracy of the program is estimated to be better than 10%.

Park, C.↗
NASA NTRS2002

Lattice Boltzmann Equation On a 2D Rectangular Grid

We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitras~y inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.

Bouzidi, MHamed↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2022

Intermediate polaronic charge transport in organic crystals from a many-body first-principles approach

Abstract Charge transport in organic molecular crystals (OMCs) is conventionally categorized into two limiting regimes − band transport, characterized by weak electron-phonon (e-ph) interactions, and charge hopping due to localized polarons formed by strong e-ph interactions. However, between these two limiting cases there is a less well understood intermediate regime where polarons are present but transport does not occur via hopping. Here we show a many-body first-principles approach that can accurately predict the carrier mobility in this intermediate regime and shed light on its microscopic origin. Our approach combines a finite-temperature cumulant method to describe strong e-ph interactions with Green-Kubo transport calculations. We apply this parameter-free framework to naphthalene crystal, demonstrating electron mobility predictions within a factor of 1.5−2 of experiment between 100 and 300 K. Our analysis reveals the formation of a broad polaron satellite peak in the electron spectral function and the failure of the Boltzmann equation in the intermediate regime.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Transition regime droplet growth and evaporation - An integrodifferential variational approach

A variational principle based on the integrodifferential form of the Boltzmann equation, which allows very general forms of the collision integral and of the boundary operator, has been used to compute the mass flux to a spherical particle. For simplicity the problem has been restricted to the one-speed, constant cross-section approximation, the black sphere problem of neutron transport theory. Using the Hilbert expansion of the solution far from the sphere as a trial function leads to a rational expression in the inverse Knudsen number for the extrapolation distance, which is simply related to the number flux to the sphere. A five-term expansion gives results that are no more than 6 percent in error when compared to an accurate numerical solution. The technique can be applied to other forms of the collision model and other boundary conditions.

Cipolla, J. W., Jr.↗
NASA NTRSDate not supplied

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗
DOE OSTI2022

The theory of transport in helical spin-structure crystals

Abstract We study helical structures in spin-spiral single crystals. In the continuum approach for the helicity potential energy the simple electronic band splits into two non-parabolic bands. For low exchange integrals, the lower band is described by a surface with a saddle shape in the direction of the helicity axis. Using the Boltzmann equation with the relaxation due to acoustic phonons, we discover the dependence of the current on the angle between the electric field and helicity axis leading to the both parallel and perpendicular to the electric field components in the electroconductivity. The latter can be interpreted as a planar Hall effect. In addition, we find that the transition rates depend on an electron spin allowing the transition between the bands. The electric conductivities exhibit nonlinear behaviors with respect to chemical potential µ . We explain this effect as the interference of the band anisotropy, spin conservation, and interband transitions. The proposed theory with the spherical model in the effective mass approximation for conduction electrons can elucidate nonlinear dependencies that can be identified in experiments. We find the excellent agreement between the theoretical and experimental data for parallel resistivity depending on temperature at the phase transition from helical to ferromagnetic state in a M n P single crystal. In addition, we predict that the perpendicular resistivity abruptly drops to zero in the ferromagnetic phase.

Physics↗
DOE OSTI2021

Novel Relaxation Time Approximation to the Relativistic Boltzmann Equation

In this work, we show that the widely used relaxation time approximation to the relativistic Boltzmann equation contains basic flaws, being incompatible with micro- and macroscopic conservation laws if the relaxation time depends on energy or general matching conditions are applied. We propose a new approximation that fixes such fundamental issues and maintains the basic properties of the linearized Boltzmann collision operator. We show how this correction affects transport coefficients, such as the bulk viscosity and particle diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗