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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗
NASA NTRS1998

Transport of Light Ions in Matter

A recent set of light ion experiments are analyzed using the Green's function method of solving the Boltzmann equation for ions of high charge and energy (the GRNTRN transport code) and the NUCFRG2 fragmentation database generator code. Although the NUCFRG2 code reasonably represents the fragmentation of heavy ions, the effects of light ion fragmentation requires a more detailed nuclear model including shell structure and short range correlations appearing as tightly bound clusters in the light ion nucleus. The most recent NTJCFRG2 code is augmented with a quasielastic alpha knockout model and semiempirical adjustments (up to 30 percent in charge removal) in the fragmentation process allowing reasonable agreement with the experiments to be obtained. A final resolution of the appropriate cross sections must await the full development of a coupled channel reaction model in which shell structure and clustering can be accurately evaluated.

Wilson, J. W.↗
NASA NTRS1994

A multi-ion generalized transport model of the polar wind

The higher-order generalizations of the equations of standard hydrodynamics, known collectively as generalized transport theories, have been used since the early 1980s to describe the terrestrial polar wind. Inherent in the structure of generalized transport theories is the ability to describe not only interparticle collisions but also certain non-Maxwellian processes, such as heat flow and viscous stress, that are characteristic of any plasma flow that is not collision dominated. Because the polar wind exhibits a transition from collision-dominated to collisionless flow, generalized transport theories possess advantages for polar wind modeling not shared by either collision-dominated models (such as standard hydrodynamics) or collisionless models (such as those based on solving the collisionless Boltzmann equation). In general, previous polar wind models have used generalized transport equations to describe electrons and only one species of ion (H(+)). If other ion species were included in the models at all, it was in a simplified or semiempirical manner. The model described in this paper is the first polar wind model that uses a generalized transport theory (bi-Maxwellian-based 16-moment theory) to describe all of the species, both major and minor, in the polar wind plasma. In the model, electrons and three ion species (H(+), He(+), O(+)) are assumed to be major and several ion species are assumed to be minor (NO(+), Fe(+), O(++)). For all species, a complete 16-moment transport formulation is used, so that profiles of density, drift velocity, parallel and perpendicular temperatures, and the field-aligned parallel and perpendicular energy flows are obtained. In the results presented here, emphasis is placed on describing those constituents of the polar wind that have received little attention in past studies. In particular, characteristic solutions are presented for supersonic H(+) outflow and for both supersonic and subsonic outflows of the major ion He(+). Solutions are also presented for various minor ions, both atomic and molecular and both singly and multiply charged.

Demars, H. G.↗
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 OSTI2022

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

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

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Kinetic effects of self-organization in low-temperature magnetized plasmas (Final Report for DE-SC0019045: 9/1/2018-8/31/2019)

The research goal of this Early Career Research Program project is to advance the science of low temperature magnetized plasmas using first-principles kinetic, e.g., solving the Boltzmann equation, models. A better understanding of the plasma transport properties, such as electrical resistivity and thermal conductivity must be obtained to better control and design magnetized plasma sources. A key physical phenomenon in low-temperature magnetized plasmas is the self-organization of plasma flows, e.g., coherent plasma structures and transition between different discharge modes. Understanding the physics of coherent plasma structures has been identified as one of the four important frontier topics in a recent workshop by the Office of Fusion Energy Science. These self-emerging plasma patterns results from the nonlinear coupling between plasma constituents (ions, electrons, and neutral atoms), electromagnetic fields, and materials. The physics of low-temperature magnetized plasmas is particularly complex because (i) the collisionless instabilities affect collisional phenomena (ionization, transport, etc.) and vice versa and (ii) the presence of materials (sheath, electron injection, etc.) influences the bulk plasma properties in low-temperature plasma (LTP) devices. Computational models of these plasma flows remain challenging because the plasma density can vary several orders of magnitudes, a wide range of temporal and spatial scales must be resolved, the phenomena are inherently three-dimensional, and various physical processes, such as wall interaction, instabilities, inelastic collisions, etc., must be simultaneously and self-consistently taken into account. Understanding the fundamental transport mechanisms in cross-field configurations enables additional controllability of plasma properties and electron energy distribution functions (EEDFs).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Validation of the SCEPTRE Boltzmann-CSD Solver

A new Boltzmann-CSD solver has been developed within the SCEPTRE radiation-transport code, based on the 1st-order form of the transport equation, using discontinuous finite elements in space and energy and discrete ordinates in angle. The Boltzmann-CSD solver has been validated against experimental data for electron energy deposition distributions and for electron emission spectra. Comparison of the calculated results with experimental data shows excellent agreement for many of the test configurations and reasonable agreement for other test configurations. The tests have also been modeled with the ITS Monte Carlo code, which also shows excellent to reasonable agreement with the SCEPTRE results and experimental data. The SCEPTRE Boltzmann-CSD solver relies on electron cross sections generated by the legacy CEPXS code, which currently is limited to electron-only Boltzmann-CSD cross sections. Performing full electron-photon radiation transport with the Boltzmann-CSD solver will require further development in the cross section generating code. For the energy- deposition calculations, neglecting photon transport results in at most about 5% overprediction of the energy deposition for high-energy electrons on high-Z targets, and relatively insignificant difference for the other test configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2000

An Improved Neutron Transport Algorithm for Space Radiation

A low-energy neutron transport algorithm for use in space radiation protection is developed. The algorithm is based upon a multigroup analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. This analysis is accomplished by solving a realistic but simplified neutron transport test problem. The test problem is analyzed by using numerical and analytical procedures to obtain an accurate solution within specified error bounds. Results from the test problem are then used for determining mean values associated with rescattering terms that are associated with a multigroup solution of the straight-ahead Boltzmann equation. The algorithm is then coupled to the Langley HZETRN code through the evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for a water and an aluminum-water shield-target configuration is then compared with LAHET and MCNPX Monte Carlo code calculations for the same shield-target configuration. The algorithm developed showed a great improvement in results over the unmodified HZETRN solution. In addition, a two-directional solution of the evaporation source showed even further improvement of the fluence near the front of the water target where diffusion from the front surface is important.

Heinbockel, John H.↗
NASA NTRS1986

Consistent transport coefficients in astrophysics

A consistent theory for dealing with transport phenomena in stellar atmospheres starting with the kinetic equations and introducing three cases (LTE, partial LTE, and non-LTE) was developed. The consistent hydrodynamical equations were presented for partial-LTE, the transport coefficients defined, and a method shown to calculate them. The method is based on the numerical solution of kinetic equations considering Landau, Boltzmann, and Focker-Planck collision terms. Finally a set of results for the transport coefficients derived for a partially ionized hydrogen gas with radiation was shown, considering ionization and recombination as well as elastic collisions. The results obtained imply major changes is some types of theoretical model calculations and can resolve some important current problems concerning energy and mass balance in the solar atmosphere. It is shown that energy balance in the lower solar transition region can be fully explained by means of radiation losses and conductive flux.

Fontenla, Juan M.↗
DOE OSTI2022

Investigating the AGN-201M research reactor's unique dominance ratio

The AGN-201M reactor is a small, thermal spectrum reactor located at the University of New Mexico (UNM). It is moderated by polyethylene, reflected by graphite, with fuel comprised of uranium microspheres embedded in polyethylene plates that are separated by an aluminum baffle. The dominance ratio of a multiplying system is the ratio of the first higher-mode eigenvalue to the fundamental, k{sub 1}/k{sub 0}. It is a key parameter in characterizing the rate of convergence, analyzing system stability, and determining how tightly coupled the system's neutron fission is. The purpose of this paper is to investigate the system's unusually low dominance ratio of 0.632. This investigation varied the reactor's physical parameters, such as the fuel geometry, fuel density, and reflector thickness, and examined their impact on the reactor's dominance ratio. Additionally, this study examined neutronically similar systems in hopes of identifying common causes for systems with low dominance ratios. The reason why the dominance ratio of the reactor is small compared to larger thermal reactors was determined to be because of its significantly smaller geometry and the composition of the fuel plates. The reflector's effect on the dominance ratio is small in comparison to the other factors, but was found to have a non-zero effect. Furthermore, the AGN-201M was found to have a significantly lower dominance ratio than systems in which it shares a very high (c{sub k} >95%) degree of neutronic similarity.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1976

Kinetic theory of turbulent compressible flows and comparison with classical theory

A complete set of turbulent correlations is given in terms of expansion coefficients in a double series of Hermite polynomials of a two-particle correlation in phase space. Only two of these coefficients, corresponding to Reynolds stress and turbulent heat flux, are shown to appear in gasdynamic equations of turbulence. This is accomplished by incorporating the phase-space correlation in the collision integral of Boltzmann's equation and by deriving generalized Navier-Stokes and Fourier transport relations. In contrast to the classical formalism, the simple expression remains invariant in form whether the flow is compressible or not. Nevertheless, a bilinear transformation of fluctuating quantities shows that the two formalisms are identical with regard to terms of double correlations. The kinetic theory justifies nonexistence of the higher-order correlations characteristic of compressible-turbulence equations in the classical regime.

Tsuge, S.↗
NASA NTRS1981

Analytical description of charged particle transport along arbitrary guiding field configurations

Aspects of charged particle transport theory have been studied by Earl (1973, 1974, 1976) in a number of investigations. The investigations were conducted under the assumption that particles propagate along a DC guiding field with negligible perpendicular transport and that pitch-angle scattering by random fields can be described by the Fokker-Planck formalism. The present study provides a generalized and simplified treatment of focused transport. The first step is taken with respect to a perturbation analysis by which the Boltzmann equation is solved in terms of successively improved approximations to its exact solution. The first approximation is the sum of isotropic pseudodiffusive and anisotropic supercoherent components whose coefficients must satisfy two coupled partial differential equations. The normalization constant K and the pseudodiffusive velocity, which appear in the first-order transport equations, are given by simple integrals over the pitch angle.

Earl, J. A.↗
DOE OSTI2021

Analysis of a New Implicit Solver for a Semiconductor Model

Here, we present and analyze a new iterative solver for implicit discretizations of a simplified Boltzmann--Poisson system. The algorithm builds on recent work that incorporated a sweeping algorithm for the Vlasov--Poisson equations as part of nested inner-outer iterative solvers for the Boltzmann--Poisson equations. The new method eliminates the need for nesting and requires only one transport sweep per iteration. It arises as a new fixed-point formulation of the discretized system which we prove to be contractive for a given electric potential. We also derive an accelerator to improve the convergence rate for systems in the drift-diffusion regime. We numerically compare the efficiency of the new solver, with and without acceleration, against a recently developed nested iterative solver.

97 MATHEMATICS AND COMPUTING↗
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↗
NASA NTRS1993

Comprehensive investigation of HgCdTe metalorganic chemical vapor deposition

The principal objective of this experimental and theoretical research program was to explore the possibility of depositing high quality epitaxial CdTe and HgCdTe at very low pressures through metalorganic chemical vapor deposition (MOCVD). We explored two important aspects of this potential process: (1) the interaction of molecular flow transport and deposition in an MOCVD reactor with a commercial configuration, and (2) the kinetics of metal alkyl source gas adsorption, decomposition and desorption from the growing film surface using ultra high vacuum surface science reaction techniques. To explore the transport-reaction issue, we have developed a reaction engineering analysis of a multiple wafer-in-tube ultrahigh vacuum chemical vapor deposition (UHV/CVD) reactor which allows an estimate of wafer or substrate throughput for a reactor of fixed geometry and a given deposition chemistry with specified film thickness uniformity constraints. The model employs a description of ballistic transport and reaction based on the pseudo-steady approximation to the Boltzmann equation in the limit of pure molecular flow. The model representation takes the form of an integral equation for the flux of each reactant or intermediate species to the wafer surfaces. Expressions for the reactive sticking coefficients (RSC) for each species must be incorporated in the term which represents reemission from a wafer surface. The interactions of MOCVD precursors with Si and CdTe were investigated using temperature programmed desorption (TPD) in ultra high vacuum combined with Auger electron spectroscopy (AES). These studies revealed that diethyltellurium (DETe) and dimethylcadmium (DMCd) adsorb weakly on clean Si(100) and desorb upon heating without decomposing. These precursors adsorb both weakly and strongly on CdTe(111)A, with DMCd exhibiting the stronger interaction with the surface than DETe.

Raupp, Gregory B.↗
NASA NTRS1988

Charged-particle transport in one dimension

A numerical solution to high-energy charged-particle transport is found by evaluation of the integral equations obtained by inverting the Boltzmann differential operator. An algorithm has been written for either continuous or discrete spectra at the boundary, which allows efficient and accurate evaluation of this rather complicated problem. Present results are compared with analytic solutions based on the perturbation theory.

NASA Discipline Radiation Health↗