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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Holographic open quantum systems: toy models and analytic properties of thermal correlators

We present a unified picture of open quantum systems, the theory of a system probing a noisy thermal environment, distilling lessons learnt from previous holographic analyses. Our treatment is applicable both when the system is coupled to short-lived (Markovian), and long-lived (non-Markovian) environmental degrees of freedom. The thermal environment is modeled using an asymptotically AdS black hole, and the systems of interest are simple probe field theories. The effective stochastic dynamics of the system is governed by real-time thermal correlators, which we compute using the gravitational Schwinger-Keldysh (grSK) geometry. We describe the structure of arbitrary tree-level contact and exchange Witten diagrams in the grSK geometry. In particular, we argue, that all such diagrams reduce to integrals supported on a single copy of the exterior of the black hole. The integrand is obtained as a multiple discontinuity of a function comprising ingoing boundary-bulk propagators, monodromy functions which appear as radial Boltzmann weights, and vertex factors. These results allow us to deduce the analytic structure of real-time thermal n-point functions in holographic CFTs. We illustrate the general statements by a two-dimensional toy model, dual to fields in the BTZ background, which we argue captures many of the essential features of generic open holographic QFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1977

Charge exchange in a planetary corona - Its effect on the distribution and escape of hydrogen

The theory for a spherical collisionless planetary corona is extended to include charge-exchange collisions between H(+) and H, which are assumed to constitute intermingled gases with different kinetic temperatures. The treatment is based on the conventional concept of a critical level (or exobase) above which the only collisions considered in the Boltzmann equation are those that resonantly exchange charge. Although the geometry treated is an oversimplification for a real planet, numerical examples are given for an idealized earth and Venus. For earth, an ion temperature of 4 times the neutral temperature, an ion density at the exobase of 14,000 per cu cm, and a plasmapause at 1.5 earth radii will raise the escape flux of H by a factor of 6. The total H above the exobase is changed by less than 1%. For Venus, conditions are examined that would account for the peculiar H distribution observed from Mariner 5. The plasma conditions required are not obviously outrageous by terrestrial standards, but the Mariner 5 ionosphere measurements did not show a high plasmapause at, say, 1.25 or 1.5 planetary radii, a fact that might argue against a charge-exchange model.

Chamberlain, J. W.↗
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.↗
DOE OSTI2021

Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics

Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

Liu, Wen‐Hao↗
DOE OSTI2025

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗
NASA NTRS1985

Nonlinear spiral density waves - Viscous damping

The formalism of Borderies, Goldreich, and Tremaine (1984), as simplified by Shu and Stewart (1985), is used to develop a theory for the viscous damping of nonlinear density waves in particulate disks of moderate collision frequency. The specific application is to Saturn's rings, but the development is general enough to allow application to a wider context (e.g., to gas clouds in a spiral galaxy). A Krook formulation is used rather than a Boltzmann formulation to treat the statistical effects of inelastic collisions. Issues that have arisen as a result of the study include a self-induced Q barrier in the first wavelength or two of the Mimas 5:3 density wave train and the surprising discovery that Saturn's B ring may behave almost as a superfluid, with hardly any viscous losses.

Shu, F. H.↗
NASA NTRS1974

Formation of ion-acoustic shocks

The formation of an ion-acoustic shock in a numerically modeled two-chamber double plasma device is investigated for a plasma of cold ions and isothermal Boltzmann electrons. An initial potential ramp applied to the driver chamber launches an ion-acoustic pulse into the target chamber which steepens into a shock. The quasi-steady shock structure agrees with observed double plasma shocks. An upper limit of Mach 1.6 is observed independent of the potential ramp magnitude, in agreement with theory.

White, R. B.↗
NASA NTRS1984

Acoustic-radiation stress in solids. I - Theory

The general case of acoustic-radiation stress associated with quasi-compressional and quasi-shear waves propagating in infinite and semiinfinite lossless solids of arbitrary crystalline symmetry is studied. The Boussinesq radiation stress is defined and found to depend directly on an acoustic nonlinearity parameter which characterizes the radiation-induced static strain, a stress-generalized nonlinearity parameter which characterizes the stress nonlinearity, and the energy density of the propagating wave. Application of the Boltzmann-Ehrenfest principle of adiabatic invariance to a self-constrained system described by the nonlinear equations of motion allows the acoustic-radiation-induced static strain to be identified with a self-constrained variation in the time-averaged product of the internal energy density and displacement gradient. The time-averaged product is scaled by the acoustic nonlinearity parameter and represents the first-order nonlinearity in the virial theorem. Finally, the relationship between the Boussinesq and the Cauchy radiation stress is obtained in a closed three-dimensional form.

Cantrell, J. H., Jr.↗
DOE OSTI2024

Elucidating the Impact of Cis – Trans Organic Structure Directing Agent Isomer Ratios on the Aluminum Distribution Within SSZ-39

Despite their widespread use, the mechanisms governing the synthesis of zeolite catalysts are still poorly understood. A notable example of this problem is the uncertainty surrounding the influence of synthesis conditions on the placement of Al atoms in the zeolite framework which determines the active sites available for catalytic species. In this work, the role of the cis to trans isomer ratio of the OSDA N,N-dimethyl-3-5-dimethylpiperidinium on the energetics of 26 distinct Al pair distributions in SSZ-39 is examined both in the presence and absence of Na using density functional theory calculations. The initial orientation of the OSDA was found to have a significant impact on the final energies present, necessitating the screening of a large number of initial orientations with force field calculations and single point DFT calculations. Ground state energies were found to vary significantly with the ratio of cis to trans OSDAs with a Boltzmann distribution revealing the most likely Al pair distributions shift from sharing the same 8 membered rings to sharing the same double six membered rings to having no shared subunits as one increases the amount of cis OSDA present within the framework. The presence of Na was found to favor Al pair distributions where both Als occupied the same 6-membered ring. When an implicit solvent model was used to evaluate ground state energies the ideal Na sites shifted from 6-membered rings to empty SSZ-39 cages while OSDA positions and orientations remained largely the same. To provide insight on how kinetic factors may influence Al distributions, formation energies we calculated for connected double six membered rings. Further, these formation energies revealed a preference for Al pairs to occupy the same 4-membered ring which indicates kinetic and thermodynamic control may lead to different Al distributions in SSZ-39.

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

Improving initialization and evolution accuracy of cosmological neutrino simulations

Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations — the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. Here, a recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart Script $\mathcal{O}$(1%) numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the 15% level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of m ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Polarized QED cascades

By taking the spin and polarization of the electrons, positrons and photons into account in the strong-field QED processes of nonlinear Compton emission and pair production, we find that the growth rate of QED cascades in ultra-intense laser fields can be substantially reduced. While this means that fewer particles are produced, we also found them to be highly polarized. We further find that the high-energy tail of the particle spectra is polarized opposite to that expected from Sokolov–Ternov theory, which cannot be explained by just taking into account spin-asymmetries in the pair production process, but results significantly from 'spin-straggling'. We employ a kinetic equation approach for the electron, positron and photon distributions, each of them spin/polarization-resolved, with the QED effects of photon emission and pair production modelled by a spin/polarization dependent Boltzmann-type collision operator. For photon-seeded cascades, depending on the photon polarization, we find an excess or a shortage of particle production in the early stages of cascade development, which provides a path towards a controlled experiment. Throughout this paper we focus on rotating electric field configuration, which represent an idealized model and allows for a straightforward interpretation of the observed effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

A Simulation to Study Speed Distributions in a Solar Plasma

We have carried out a numerical simulation of a plasma with characteristics similar to those found in the core of the Sun. Particular emphasis is placed on the Coulomb interaction between the ions and electrons, which could result in a relative velocity distribution different from the Maxwell-Boltzmann (MB) distribution generally assumed for a plasma. The fact that the distribution may not exactly follow the MB distribution could have very important consequences for a variety of problems in solar physics, especially the neutrino problem. Very briefly, the neutrino problem is that the observed neutrino detections from the Sun are smaller than what the standard solar theory predicts. In Section I we introduce the problem and in section II we discuss the approach to try to solve the problem: i.e., a molecular dynamics approach. In section III we provide details about the integration method, and any simplifications that can be applied to the problem. In section IV (the core of this report) we state our results. First for the specific case of 1000 particles and then for other cases with different number of particles. In section V we summarize our findings and state our conclusions. Sections VI VII and VIII provide the list of figures, reference material and acknowledgements respectively.

Cheeseman, Peter↗
NASA NTRS1999

A Simulation to Study Speed Distributions in a Solar Plasma

We (Peter Cheeseman of NASA Ames/Caelum Research) & Jose Luis Alvarellos of the SJSU Physics Department/SJSU Foundation.) have carried out a numerical simulation of a plasma with characteristics similar to those found in the core of the Sun. Particular emphasis is placed on the Coulomb interaction between the ions and electrons, which could result in a relative velocity distribution different from the Maxwell-Boltzmann (MB) distribution generally assumed for a plasma. The fact that the distribution may not exactly follow the MB distribution could have very important consequences for a variety of problems in solar physics, especially the neutrino problem. Very briefly. the neutrino problem is that the observed neutrino detections from the Sun are smaller than what the standard solar theory predicts. In Section 1 we introduce the problem and in section 2 we discuss the approach to try to solve the problem: i.e., a molecular dynamics approach. In section 3 we provide details about the integration method, and any simplifications that can be applied to the problem. In section 4 (the core of this report) we state our results, first for the specific case of 1000 particles and then for other cases with different number of particles. In section 5 we summarize our findings and state our conclusions. Sections 6 and 7 provide the list of figures, reference material and acknowledgments respectively.

Cheeseman, Peter↗
NASA NTRS1979

Ion velocity distributions in the high-latitude ionosphere

The theory and observations relating to ion velocity distributions in the high-latitude F region are reviewed. The review covers three basic aspects, including methods for calculating ion velocity distributions, experimental evidence for the existence of non-Maxwellian ion velocity distributions, and the ionospheric consequences of non-Maxwellian distributions. Specific topics covered include closed form solutions to Boltzmann's equation for a simple collision model, series solutions to the Boltzmann equation for more general collision models, retarding potential analyzer measurements, excitation of VLF electrostatic emissions, ion-neutral reaction rates, and transport effects.

St-Maurice, J.-P.↗
DOE OSTI2023

Superconductivity in a strange metal

In a gas of charged particles with the density n, mass m, and charge e, the electrical conductivity σ is given by the Drude formula: σ = ne 2 /mΓ, where Γ is the scattering rate. Standard metals are well described by Landau's Fermi Liquid (FL) theory, in which the electric current is carried by quasi-particles, low-energy excitations of the FL that resemble electrons with some effective mass m*. The scattering rate Γ= 1/τ= v F /l, where t is the (momentum) relaxation) time, l is the mean-free path, and v F is the Fermi velocity, can be expressed ("Matthiessen's Rule") as a sum of contributions from various scattering channels: Γ = Γ 0 + Γ el-el + Γ el-ph + …, where Γ 0 describes scattering on lattice imperfections, Γ el-el the electron-electron scattering, Γ el-ph the electron-phonon scattering, etc. Of these, Γ 0 = v F /l 0 , where l 0 is the average distance between the defects, is temperature-independent. Γ el-el should scale as T 2 because of Fermi statistics; for two electrons to scatter on one another, both must come from the "Debye shell" of the width k B T/E F , where k B is the Boltzmann constant and E F is the Fermi energy. Γ el-ph typically grows as T 5 , so we expect this to overwhelm the other channels at a high enough T. However, since l cannot be shorter than the distance between the atoms, the total Γ saturates at Mott-Ioffe-Regel (MIR) limit, roughly v F /a 0 , where v F is the Fermi velocity and a 0 is the lattice constant. The resistivity should also saturate at low T, at ρ 0 = m*v F /ne 2 l 0 T. FL theory also describes other electronic properties; e.g., it predicts that in the magnetic field B, the resistivity of the metal should increase with B 2 , because σ(B) = σ(B=0)/(1 + (ω c /Γ) 2 ), where ω c = eB/m*.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Measuring Galactic dark matter through unsupervised machine learning

ABSTRACT Measuring the density profile of dark matter in the Solar neighbourhood has important implications for both dark matter theory and experiment. In this work, we apply autoregressive flows to stars from a realistic simulation of a Milky Way-type galaxy to learn – in an unsupervised way – the stellar phase space density and its derivatives. With these as inputs, and under the assumption of dynamic equilibrium, the gravitational acceleration field and mass density can be calculated directly from the Boltzmann equation without the need to assume either cylindrical symmetry or specific functional forms for the galaxy’s mass density. We demonstrate our approach can accurately reconstruct the mass density and acceleration profiles of the simulated galaxy, even in the presence of Gaia-like errors in the kinematic measurements.

Astronomy & Astrophysics↗
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 OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗