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 235 records · Page 13

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗

Voltage-Dependent First-Principles Simulation of Insertion of Chloride Ions into Al/Al 2 O 3 Interfaces Using the Quantum Continuum Approximation

Experiments have shown that pitting corrosion can develop in aluminum surfaces at potentials > − 0.5 V relative to the standard hydrogen electrode (SHE). Until recently, the onset of pitting corrosion in aluminum has not been rigorously explored at an atomistic scale because of the difficulty of incorporating a voltage into density functional theory (DFT) calculations. We introduce the Quantum Continuum Approximation (QCA) which self-consistently couples explicit DFT calculations of the metal-insulator and insulator-solution interfaces to continuum Poisson-Boltzmann electrostatic distributions describing the bulk of the insulating region. By decreasing the number of atoms necessary to explicitly simulate with DFT by an order of magnitude, QCA makes the first-principles prediction of the voltage of realistic electrochemical interfaces feasible. After developing this technique, we apply QCA to predict the formation energy of chloride atoms inserting into oxygen vacancies in Al(111)/α-Al 2 O 3 (0001) interfaces as a function of applied voltage. We predict that chloride insertion is only favorable in systems with a grain boundary in the Al 2 O 3 for voltages > − 0.2 V (SHE). Here our results roughly agree with the experimentally demonstrated onset of corrosion, demonstrating QCA's utility in modeling realistic electrochemical systems at reasonable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗

Stochastic renormalization group and gradient flow

A non-perturbative and continuous definition of RG transformations as stochastic processes is proposed, inspired by the observation that the functional RG equations for effective Boltzmann factors may be interpreted as Fokker-Planck equations. The result implies a new approach to Monte Carlo RG that is amenable to lattice simulation. Long-distance correlations of the effective theory are shown to approach gradient-flowed correlations, which are simpler to measure. The Markov property of the stochastic RG transformation implies an RG scaling formula which allows for the measurement of anomalous dimensions when transcribed into gradient flow expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scale-invariant machine-learning model accelerates the discovery of quaternary chalcogenides with ultralow lattice thermal conductivity

We design an advanced machine-learning (ML) model based on crystal graph convolutional neural network that is insensitive to volumes (i.e., scale) of the input crystal structures to discover novel quaternary chalcogenides, AMM'Q 3 (A/M/M' = alkali, alkaline earth, post-transition metals, lanthanides, and Q = chalcogens). These compounds are shown to possess ultralow lattice thermal conductivity (κ l ), a desired requirement for thermal-barrier coatings and thermoelectrics. Upon screening the thermodynamic stability of ~1 million compounds using the ML model iteratively and performing density-functional theory (DFT) calculations for a small fraction of compounds, we discover 99 compounds that are validated to be stable in DFT. Taking several DFT-stable compounds, we calculate their κ l using Peierls–Boltzmann transport equation, which reveals ultralow κ l (<2 Wm -1 K -1 at room temperature) due to their soft elasticity and strong phonon anharmonicity. Our work demonstrates the high efficiency of scale-invariant ML model in predicting novel compounds and presents experimental-research opportunities with these new compounds.

36 MATERIALS SCIENCE↗

Simulating excited-state complex ensembles: Fluorescence and solvatochromism in amine-arene exciplexes

Exciplexes are excited-state complexes formed as a result of partial charge transfer from the donor to the acceptor species when one moiety of the donor–acceptor pair is electronically excited. The arene–amine exciplex formed between oligo-(p-phenylene) (OPP) and triethylamine (TEA) is of interest in the catalytic photoreduction of CO2 because it can compete with complete electron transfer to the OPP catalyst. Therefore, formation of the exciplex can hinder the generation of a radical anion OPP·− necessary for subsequent CO2 reduction. We report an implementation of a workflow automating quantum-chemistry calculations that generate and characterize an ensemble of structures to represent this exciplex state. We use FireWorks, Pymatgen, and Custodian Python packages for high-throughput ensemble generation. The workflow includes time-dependent density functional theory optimization, verification of excited-state minima, and exciplex characterization with natural transition orbitals, exciton analysis, excited-state Mulliken charges, and energy decomposition analysis. Fluorescence spectra computed for these ensembles using Boltzmann-weighted contributions of each structure agree better with experiment than our previous calculations based on a single representative exciplex structure [Kron et al., J. Phys. Chem. A 126, 2319–2329 (2022)]. The ensemble description of the exciplex state also reproduces an experimentally observed red shift of the emission spectrum of [OPP-4–TEA]* relative to [OPP-3–TEA]*. The workflow developed here streamlines otherwise labor-intensive calculations that would require significant user involvement and intervention.

Chemistry↗

An Improved Green's Function for Ion Beam Transport

Ion beam transport theory allows testing of material transmission properties in the laboratory environment generated by particle accelerators. This is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are the main emphasis of the present work. In consequence, an analytic solution of the linear Boltzmann equation is pursued in the form of a Green's function allowing flexibility in application to a broad range of boundary value problems. It has been established that simple solutions can be found for the high charge and energy (HZE) by ignoring nuclear energy downshifts and dispersion. Such solutions were found to be supported by experimental evidence with HZE ion beams when multiple scattering was added. Lacking from the prior solutions were range and energy straggling and energy downshift with dispersion associated with nuclear events. Recently, we have found global solutions including these effects providing a broader class of HZE ion solutions.

Tweed, J.↗

An improved Green's function for ion beam transport

Ion beam transport theory allows testing of material transmission properties in the laboratory environment generated by particle accelerators. This is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are the main emphasis of the present work. In consequence, an analytic solution of the linear Boltzmann equation is pursued in the form of a Green's function allowing flexibility in application to a broad range of boundary value problems. It has been established that simple solutions can be found for high charge and energy (HZE) ions by ignoring nuclear energy downshifts and dispersion. Such solutions were found to be supported by experimental evidence with HZE ion beams when multiple scattering was added. Lacking from the prior solutions were range and energy straggling and energy downshift with dispersion associated with nuclear events. Recently, we have found global solutions including these effects providing a broader class of HZE ion solutions. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

NASA Center LaRC↗

Blending and cavern formation within non-Newtonian fluids in stirred tanks: Application to nuclear waste fluid processing

A numerical approach for predicting the time-accurate fluid flow and mixing properties of non-Newtonian fluids is presented. This approach, which is based on the lattice Boltzmann method, is used to characterize the blending in shear-thinning fluids, as well as cavern formation in yield-stress fluids. Predictions compare favorably to measured data and expectations from first-principles theory. Importantly, across the range of fluids and systems here, the simulations require no reparameterization or retuning between scenarios. This generality is exploited to model a nuclear waste processing unit operation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Algorithm for the Transport of Anisotropic Neutrons

One major obstacle to human space exploration is the possible limitations imposed by the adverse effect of long-term exposure to the space environment. Even before human spaceflight began, the potentially brief exposure of astronauts to the very intense random solar particle events (SPE) were of great concern. A new challenge appears in deep space exploration from exposure to the low-intensity heavy-ion flux of the galactic cosmic rays (GCR) since the missions are of long duration and the accumulated GCR exposures can be high. Because cancer induction rates increase behind low to rather large thicknesses of aluminum shielding, according to available biological data on mammalian exposures to GCR like ions, the shield requirements for a Mars mission are prohibitively expensive in terms of mission launch costs. Therefore, a critical issue in the Human Exploration and Development of Space enterprise is cost effective mitigation of risk associated with ionizing radiation exposure. In order to estimate astronaut risk to GCR exposure and associated cancer risks and health hazards, it is necessary to do shield material studies. To determine an optimum radiation shield material it is necessary to understand nuclear interaction processes such as fragmentation and secondary particle production which is a function of energy dependent cross sections. This requires knowledge of material transmission characteristics either through laboratory testing or improved theoretical modeling. Here ion beam transport theory is of importance in that testing of materials in the laboratory environment generated by particle accelerators is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are a major emphasis of the present work.

Tweed, J.↗

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↗

Entropy Generation Across Earth's Bow Shock

Earth's bow shock is a transition layer that causes an irreversible change in the state of plasma that is stationary in time. Theories predict entropy increases across the bow shock but entropy has never been directly measured. Cluster and Double Star plasma experiments measure 3D plasma distributions upstream and downstream of the bow shock that allow calculation of Boltzmann's entropy function H and his famous H-theorem, dH/dt ~ O. We present the first direct measurements of entropy density changes across Earth's bow shock. We will show that this entropy generation may be part of the processes that produce the non-thermal plasma distributions is consistent with a kinetic entropy flux model derived from the collisionless Boltzmann equation, giving strong support that solar wind's total entropy across the bow shock remains unchanged. As far as we know, our results are not explained by any existing shock models and should be of interests to theorists.

Parks, George K.↗

Semiconductor Nanostructures by Scientific Design (Final Technical Report)

We developed and used quantum simulation techniques to predict structural, electronic, optical and thermal properties affecting the scientific design of semiconductor nanostructures. We focused on the study of surfaces and interfaces at the nanoscale, in particular embedding and solvation of semiconductor nanoparticles, and provided a microscopic description of nanostructured materials in realistic environments, directly comparable with experimental conditions. We addressed three major problems: (i) the design of realistic structural models of integrated nanostructures, e.g. quantum dots embedded in solid matrices and nanostructures immersed in fluids; (ii) the development of efficient and accurate calculations to describe the excited state properties of nanostructures, beyond standard Density Functional Theory, so as to understand and predict electronic properties involved in renewable energy applications, for example solar conversion devices. In particular we developed and used methods based on a many body perturbation theory. (iii) the development of calculations of heat transport properties of nanostructured materials, for thermoelectric applications, including predictive calculations of thermal and electrical conductivity and of Seebeck coefficients. Specifically, we developed and used methods based on the Boltzmann transport equation and molecular dynamics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗

Casimir effect in Very Special Relativity at finite temperature

Recently a great deal of interest in field theories formulated in a Lorentz violating framework has been developed. Here the Very Special Relativity (VSR) is considered. The main aspect in the VSR proposal is that laws of physics are invariant under the subgroups of the Poincaré group preserving the basic elements of special relativity. An important point is that the photon acquires a mass. In this context, the energy–momentum tensor for the electromagnetic field is calculated. From this, the Stefan–Boltzmann law and Casimir effect at finite temperature in VSR are obtained. The effects of temperature are introduced using the Thermo Field Dynamics (TFD) formalism. A comparative analysis with the Casimir effect for the standard electromagnetic case is developed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Fate of Nucleated Black Holes in de Sitter Quantum Gravity

The Euclidean Nariai geometry has long been proposed as the instanton describing the nucleation of maximal-mass black holes in de Sitter space. We place this interpretation on firmer footing by showing that, once an observer is included, the gravitational path integral produces the imaginary phase required for a transition rate. As a warmup, we revisit the Hawking-Moss instanton and, as a byproduct, find that scalar fields can enhance black-hole nucleation, suggesting a quantum-gravity bound on scalar potentials with de Sitter solutions. We then study the subsequent semiclassical evolution of the nucleated black hole. We show that the previously claimed "anti-evaporation" channel is unphysical, arising from a quantum state with singular horizons. In a smooth state, the black hole instead undergoes standard thermal Hawking evaporation. We verify explicit agreement with the no-boundary state and argue that this evaporation is not subject to large quantum-gravity corrections. The nucleated black hole thus evaporates completely back to the maximally-entropic empty de Sitter vacuum, making the full process a Boltzmann fluctuation.

FOS: Physical sciences↗

First-Principles Modeling of IAT-Driven Anomalous Resistivity in Hollow Cathode Discharges I: Theory

Self-consistent expressions are derived for the impact of ion acoustic turbulence (IAT) on the fluid properties of electrons and ions in the plume of a hollow cathode. This instability is believed to have a dominant role in the plasmadynamics in this region. The random phase approximation is applied to the perturbed Boltzmann equation to determine the effective drag and heating of the plasma species induced by the IAT. These driving terms are then incorporated into momentum and energy fluid equations for the ions and electrons. A fluid-like relation governing the evolution of the total IAT wave energy density is also presented. The physical significance of the IAT-based terms is discussed in the context of energy and momentum conservation, and a path forward is outlined for implementing the effects of IAT in the fluid-based model in Part II.

Jorns, Benjamin A.↗

Transport coefficients of gaseous ions in an electric field

A general theory of ion mobility formulated by Kihara (1953) is extended to ion diffusion and to mixtures of neutral gases. The theory assumes that only binary collisions between ions and neutral particles need to be taken into account and that the velocity distribution function of the neutral particles is Maxwellian. These assumptions make it possible to use a linearized Boltzmann equation. Questions of mobility are considered along with aspects of diffusion and deviations from Fick's law of diffusion.

Whealton, J. H.↗