Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boltzmann theory”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Flow of rarefied gases over two-dimensional bodies

A kinetic-theory analysis is made of the flow of rarefied gases over two-dimensional bodies of arbitrary curvature. The Boltzmann equation simplified by a model collision integral is written in an arbitrary orthogonal curvilinear coordinate system, and solved by means of finite-difference approximation with the discrete ordinate method. A numerical code is developed which can be applied to any two-dimensional submerged body of arbitrary curvature for the flow regimes from free-molecular to slip at transonic Mach numbers. Predictions are made for the case of a right circular cylinder.

Jeng, Duen-Ren↗
NASA NTRS1999

A Scale-Invariant Model of Statistical Mechanics and Modified Forms of the First and the Second Laws of Thermodynamics

A scale-invariant statistical theory of fields is presented that leads to invariant definition of density, velocity, temperature, and pressure, The definition of Boltzmann constant is introduced as k(sub k) = m(sub k)v(sub k)c = 1.381 x 10(exp -23) J x K(exp -1), suggesting that the Kelvin absolute temperature scale is equivalent to a length scale. Two new state variables called the reversible heat Q(sub rev) = TS and the reversible work W(sub rev) = PV are introduced. The modified forms of the first and second law of thermodynamics are presented. The microscopic definition of heat (work) is presented as the kinetic energy due to the random (peculiar) translational, rotational, and pulsational motions. The Gibbs free energy of an element at scale Beta is identified as the total system energy at scale (Beta-1), thus leading to an invariant form of the first law of thermodynamics U(sub Beta) = Q(sub Beta) - W(sub Beta) +N(e3)U(sub Beta-1).

Sohrab, Siavash H.↗
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

(Photo)Electrocatalytic CO 2 Reduction at the Defective Anatase TiO 2 (101) Surface

Excessive carbon dioxide (CO 2 ) emissions by combustion of fossil fuels are linked to global warming and rapid climate change. One promising route to lowering the concentration of CO 2 in the atmosphere is to reduce it to useful small molecules via photoelectrocatalytic hydrogenation, which would enable solar energy storage with a zero-carbon emission cycle and perform a more efficient separation of the photogenerated electron and hole pair than pure photocatalysis. Indeed, photoelectrocatalytic CO 2 reduction has been an intense focus of research. Using the density functional theory (DFT), we studied the CO 2 reduction reaction on the defective anatase TiO 2 (101) surface, at both the solvent/catalyst and the electrolyte/catalyst interfaces. The analysis of the electronic structure of the surface shows a contrast between the solvent/catalyst and the electrolyte/catalyst interfaces, which results in the two corresponding catalytic cycles being distinct. Furthermore, our study explains at the electronic and mechanistic levels why methanol is the main product in the presence of the electrolyte and why the overpotential is not only controlled by the reaction process but also by the diffusion process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

QCD axion: Some like it hot

We compare the quantum chromodynamics (QCD) axion phase-space distribution from unitarized next-to-leading order chiral perturbation theory with the one extracted from pion-scattering data. We derive a robust bound by confronting momentum-dependent Boltzmann equations against up-to-date observations of the cosmic microwave background, of the baryonic acoustic oscillations and of primordial abundances. These datasets imply 𝑚 𝑎 ≤ 0.16 eV for the 95% credible interval, i.e., ∼ 30% stronger bound than what previously found. We present forecasts using dedicated likelihoods for future cosmological surveys and the sphaleron rate from unquenched lattice QCD.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A new approach for resonance treatment of doubly heterogeneous fuel using the RSE method

A new resonance calculation method for the doubly-heterogeneous (DH) fuels such as high- temperature gas-cooled reactor fuel is proposed based on the resonance calculation based on Spectral Expansion (RSE) method. The concept of pointwise disadvantage factor for fuel grain is taken into account to treat the DH fuels. The verification calculation is carried out for simplified single fuel cell and fuel compact consisting of five fuel cells and graphite moderator. The calculation results indicate that the present method can appropriately handle the space- dependent self-shielding effect for DH fuels. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Performance enhancement of planar MOC calculation in hexagonal geometries by assembly-wise domain decomposition

The parallelization of planar Method of Characteristics (MOC) calculation with the higher order scattering in the hexagonal version of the nTRACER direct whole core calculation code is enhanced by adopting assembly-wise domain decomposition and other techniques. In order to augment the hexagonal domain decomposition, the three-color scheme, the angular decomposition, and the angular flux storage scheme are examined. The parallel performance assessed for the VVER benchmark problems demonstrates that the hexagonal assembly-wise domain decomposition significantly reduces the computing time in ray tracing by up to 64 %. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Variance reduction techniques for Monte Carlo neutron noise simulations

The small fluctuations of the neutron flux caused by small perturbations of the macroscopic cross-sections take the name of neutron noise. Advanced Monte Carlo methods have been recently proposed in order to solve the neutron noise equations in the frequency domain, which allows establishing reference solutions to validate faster but approximate deterministic solvers. Due to the presence of particles carrying two statistical weights (for the real and imaginary components of the noise field), both of which may be positive or negative, Monte-Carlo simulations of neutron noise pose distinct challenges in terms of variance reduction. In this work we investigate two variance-reduction techniques, namely branchless collisions and weight cancellation, and probe their effectiveness for a benchmark con- figuration concerning the noise field induced by a pin with oscillating cross sections in a fuel assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

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↗
NASA NTRS1973

The formation of ion acoustic shocks

Recent experiments performed in the double plasma (DP) device have verified the existence of electrostatic ion acoustic laminar shocks. The influence of the piston on the shock structure is investigated by modeling the DP device and by numerically solving the temporal and spatial evolution of the shock. In order to isolate piston effects, as opposed to kinetic theory effects such as reflected ions and trapped electrons, the DP plasma is modeled as a cold ion fluid with isothermal Boltzmann electrons. It is shown that laminar shock transitions with structure agreeing with DP shock experiments can be excited.

White, R. B.↗
NASA NTRS1981

Information theory lateral density distribution for Earth inferred from global gravity field

Information Theory Inference, better known as the Maximum Entropy Method, was used to infer the lateral density distribution inside the Earth. The approach assumed that the Earth consists of indistinguishable Maxwell-Boltzmann particles populating infinitesimal volume elements, and followed the standard methods of statistical mechanics (maximizing the entropy function). The GEM 10B spherical harmonic gravity field coefficients, complete to degree and order 36, were used as constraints on the lateral density distribution. The spherically symmetric part of the density distribution was assumed to be known. The lateral density variation was assumed to be small compared to the spherically symmetric part. The resulting information theory density distribution for the cases of no crust removed, 30 km of compensated crust removed, and 30 km of uncompensated crust removed all gave broad density anomalies extending deep into the mantle, but with the density contrasts being the greatest towards the surface (typically + or 0.004 g cm 3 in the first two cases and + or - 0.04 g cm 3 in the third). None of the density distributions resemble classical organized convection cells. The information theory approach may have use in choosing Standard Earth Models, but, the inclusion of seismic data into the approach appears difficult.

Rubincam, D. P.↗
DOE OSTI2022

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↗
DOE OSTI2022

Reduced-order modeling of neutron transport separated in energy by Minimax Proper Generalized Decomposition

In this article, we demonstrate a Petrov-Galerkin Proper Generalized Decomposition (PGD) known as Minimax PGD for modeling neutron transport separated in energy. To compare the Minimax with the classical Galerkin PGD, we assess both on a model problem of UO{sub 2} or Mixed Oxide (MOX) fuel pins in an infinite lattice with 3 industry-standard energy meshes. We find the Minimax PGD achieves a superior decomposition to Galerkin PGD, both with and without update of the energy modes. This suggests Minimax PGD may be more computationally efficient, provided this reduction in modes (to achieve a given accuracy) outweighs the cost of solving the necessary adjoint problems. In either case, we note that PGD offers an a priori Reduced-Order Model (ROM) which may be dramatically cheaper to solve than the full-order model, especially in problems with fine to ultrafine energy meshes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1977

Quantum statistical theory of semiconductor junctions in thermal equilibrium

Free carrier and electric field distributions of one-dimensional semiconductor junctions are evaluated using a quantum mechanical phase-space distribution and its corresponding Boltzmann equation. Attention is given to quantum and exchange corrections in cases of high doping concentrations when carrier densities become degenerate. Quantitative differences between degenerate and classical junction characteristics, e.g., maximum electric field and built-in voltage and carrier concentration within the transition region, are evaluated numerically.

Von Roos, O.↗
DOE OSTI2025

Exploring the energy landscape of RBMs: reciprocal space insights into bosons, hierarchical learning and symmetry breaking

Deep generative models have become ubiquitous due to their ability to learn and sample from complex distributions. Despite the proliferation of various frameworks, the relationships among these models remain largely unexplored, a gap that hinders the development of a unified theory of AI learning. In this work, we address two central challenges: clarifying the connections between different deep generative models and deepening our understanding of their learning mechanisms. We focus on Restricted Boltzmann Machines (RBMs), a class of generative models known for their universal approximation capabilities for discrete distributions. By introducing a reciprocal space formulation for RBMs, we reveal a connection between these models, diffusion processes, and systems of coupled bosons. Our analysis shows that at initialization, the RBM operates at a saddle point, where the local curvature is determined by the singular values of the weight matrix, whose distribution follows the Marc̆enko-Pastur law and exhibits rotational symmetry. During training, this rotational symmetry is broken due to hierarchical learning, where different degrees of freedom progressively capture features at multiple levels of abstraction. This leads to a symmetry breaking in the energy landscape, reminiscent of Landau’s theory. This symmetry breaking in the energy landscape is characterized by the singular values and the weight matrix eigenvector matrix. We derive the corresponding free energy in a mean-field approximation. We show that in the limit of infinite size RBM, the reciprocal variables are Gaussian distributed. Our findings indicate that in this regime, there will be some modes for which the diffusion process will not converge to the Boltzmann distribution. To illustrate our results, we trained replicas of RBMs with different hidden layer sizes using the MNIST dataset. Our findings not only bridge the gap between disparate generative frameworks but also shed light on the fundamental processes underpinning learning in deep generative models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering

Angularly resolved Thomson scattering is a novel extension of Thomson scattering, enabling the measurement of the electron velocity distribution function over many orders of magnitude. Here, details of the theoretical basis of the technique and the instrument designed for this measurement are described. Angularly resolved Thomson-scattering data from several experiments are shown with descriptions of the corresponding distribution functions. A reduced model describing the distribution function is given and used to perform a Monte Carlo analysis of the uncertainty in the measurements. The electron density and temperature were determined to a precision of 12% and 21%, respectively, on average, while all other parameters defining the distribution function were generally determined to better than 20%. It was found that these uncertainties were primarily due to limited signal to noise and instrumental effects. Furthermore, measurements with this level of precision were sufficient to distinguish between Maxwellian and non-Maxwellian distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

On the Cauchy problem for Boltzmann equation modeling a polyatomic gas

In the present article, we consider the Boltzmann equation that models a polyatomic gas by introducing one additional continuous variable, referred to as microscopic internal energy. We establish existence and uniqueness theory in the space homogeneous setting for the full non-linear case, under an extended Grad-type assumption on transition probability rates, which comprises hard potentials for both the relative speed and internal energy with the rate in the interval 0,2, multiplied by an integrable angular part and integrable partition functions. The Cauchy problem is resolved by means of an abstract ordinary differential equation (ODE) theory in Banach spaces for the initial data with finite and strictly positive gas mass and energy, finite momentum, and additionally finite K* polynomial moment, with K* depending on the rate of the transition probability and the structure of a polyatomic molecule or its internal degrees of freedom. Moreover, we prove that polynomially and exponentially weighted Banach space norms associated with the solution are both generated and propagated uniformly in time.

Physics↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗