Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice Boltzmann”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Kohn anomalies and phonon anharmonicity in iridium

Elemental iridium presents surprising challenges for both inelastic neutron scattering (INS) and theoretical thermal transport calculations due to its high neutron absorption cross-section and strong electron-phonon interactions, respectively. Here, in this study, we overcome these challenges to measure temperature-dependent phonon dispersion curves, compare these with calculations based on density functional theory (DFT), and ultimately examine the electron-phonon limited transport behaviors of this material. Our DFT calculations demonstrate Kohn anomalies, near the 𝐾 point of the iridium Brillouin zone, indicating coupling between electrons and phonons. Strong electron-phonon coupling can compete with anharmonic effects to determine electrical and thermal transport behaviors and make the Kohn anomalies challenging to observe. Nonetheless, our INS measurements map these anomalies and other dispersion features over the Brillouin zone from 100 to 700 K. These measurements also uncover unexpectedly large mode specific Grüneisen parameters obtained from the temperature-dependent phonon energies, highlighting strong anharmonicity in iridium. DFT-based Boltzmann transport calculations demonstrate how anharmonicity and electron-phonon couplings determine electronic and lattice transport behaviors. Furthermore, we correlate the Kohn anomalies with calculated electron-phonon nesting functions, Fermi surfaces, and DFT-derived coupling strengths. This study provides detailed insights into the temperature-dependent mode-resolved lattice dynamics and anharmonicity, transport behavior, and electron-phonon interactions.

DFT↗
NASA NTRS1981

Langmuir probe theory and the problem of anisotropic collection

A more general understanding of Langmuir probe theory by going beyond particular coordinate systems as far as possible was studied. To delineate the role of initial velocities in the problem and to characterize the interaction between the velocity flow field and the electrostatic potential field both of which enter into the determination of particle trajectories are desired. The probe is supposed to be negatively charged throughout and the ensuing region of positive space charge about it differs fundamentally from the plasma region beyond in the sense that the former represents a region in which deterministic trajectories dominate while the lattice is a region dominated by statistical interactions largely controlled by the collisionless stationary Boltzmann equation. A general solution of this Boltzmann equation is given and its relationship to the problem of saturation current density in 3D is elucidated.

Kyrala, A.↗
DOE OSTI2021

Accelerated discovery of a large family of quaternary chalcogenides with very low lattice thermal conductivity

The development of efficient thermal energy management devices such as thermoelectrics and barrier coatings often relies on compounds having low lattice thermal conductivity (κ l ). Here, we present the computational discovery of a large family of 628 thermodynamically stable quaternary chalcogenides, AMM'Q 3 (A = alkali/alkaline earth/post-transition metals; M/M' = transition metals, lanthanides; Q = chalcogens) using high-throughput density functional theory (DFT) calculations. We validate the presence of low κ l in these materials by calculating κ l of several predicted stable compounds using the Peierls–Boltzmann transport equation. Our analysis reveals that the low κ l originates from the presence of either a strong lattice anharmonicity that enhances the phonon-scatterings or rattler cations that lead to multiple scattering channels in their crystal structures. Our thermoelectric calculations indicate that some of the predicted semiconductors may possess high energy conversion efficiency with their figure-of-merits exceeding 1 near 600 K. Our predictions suggest experimental research opportunities in the synthesis and characterization of these stable, low κ l compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗
DOE OSTI2022

Thermalization of Fluorescent Protein Exciton–Polaritons at Room Temperature

Fluorescent proteins (FPs) have recently emerged as a serious contender for realizing ultralow threshold room temperature exciton–polariton condensation and lasing. Here, this contribution investigates the thermalization of FP microcavity exciton–polaritons upon optical pumping under ambient conditions. Polariton cooling is realized using a new FP molecule, called mScarlet, coupled strongly to the optical modes in a Fabry–Pérot cavity. Interestingly, at the threshold excitation energy (fluence) of ≈9 nJ per pulse (15.6 mJ cm –2 ), an effective temperature is observed, T eff ≈ 350 ± 35 K close to the lattice temperature indicative of strongly thermalized exciton–polaritons at equilibrium. This efficient thermalization results from the interplay of radiative pumping facilitated by the energetics of the lower polariton branch and the cavity Q-factor. Direct evidence for dramatic switching from an equilibrium state into a metastable state is observed for the organic cavity polariton device at room temperature via deviation from the Maxwell–Boltzmann statistics at k || = 0 above the threshold. Thermalized polariton gases in organic systems at equilibrium hold substantial promise for designing room temperature polaritonic circuits, switches, and lattices for analog simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles determination of the phonon-point defect scattering and thermal transport due to fission products in ThO 2

This work presents the first-principles calculations of the lattice thermal conductivity degradation due to point defects in thorium dioxide using an iterative solution of the Pierels-Boltzmann transport equation. We have used the non perturbative Green’s function methodology to compute the phonon- point defect scattering rates that consider the local distortion around the point defect, including the mass difference changes, interatomic force constants, and structural relaxation near the point defects. The point defects considered in this work include the vacancy of thorium ($V$ $TH$ ) and oxygen ($V$ $O$ ), substitutions of helium ($He$ $Th$ ), krypton ($Kr$ $Th$ ), zirconium ($Zr$ $Th$ ), iodine ($I$ $Th$ ) and xenon ($Xe$ $Th$ ) in the thorium site, and the three different configurations of the Schottky defects. The results of the phonon-defect scattering rate reveal that, among all the considered intrinsic defects, the thorium vacancy and helium substitution in the thorium site scatter the phonon most due to the substantial changes in the force constant and structural distortions. The scattering of phonons due to the substitutional defects unveils that the zirconium atom scatters phonons the least, followed by xenon, iodine, krypton, and helium. This is contrary to the intuition that the scattering strength follows $He$ $Th$ > $Kr$ $Th$ > $Zr$ $Th$ > $I$ $Th$ > $Xe$ $Th$ based on the mass difference. This striking difference in the zirconium phonon scattering is due to the local chemical environment changes. Zirconium is an electropositive element with valency similar to thorium and, therefore, can bond with the oxygen atoms, thus creating less force constant variance compared to iodine, an electronegative element noble gas helium, xenon, and krypton. In conclusion, these results can serve as the benchmark for the analytical models and help the engineering-scale modeling effort for nuclear design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Energy amplification in plasma crystals due to multiple torsions

Interacting torsions are examined within a two-dimensional monolayer crystal suspended in an argon complex plasma for 1–10 W discharge powers and pressures of 135–155 mTorr. Two torsions embedded in a lattice are shown to amplify the kinetic energy and range of motion of particles located between the torsions to nearly three times that observed in single torsion systems. It is also shown that multiple torsions can interact via amplified particle energy when separated by up to 14 interparticle distances (Δ). The torsion separation distance also showed a positive linear trend with power and a slightly positive correlation with the pressure. This amplification of energy is possible due to the fact that multiple torsions in a lattice increase the interparticle distance of the lattice by 16% more than single torsion systems, leading to additional freedom of motion in the lattice plane. These combined findings show that multiple torsions heat the lattice differently depending on their separation from the other torsion. The midpoint particles between torsions absorb the majority of energy from the two torsions, and energy addition at the midpoint is nonlinear. The addition of more torsions to the lattice may lead to melting of the plasma crystal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗
NASA NTRS1982

Modeling of thin film GaAs growth

A potential scaling Monte Carlo model of crystal growth is developed. The model is a modification of the solid-on-solid method for studying crystal growth in that potentials at surface sites are continuously updated on a time scale reflecting the surface events of migration, incorporation and evaporation. The model allows for B on A type of crystal growth and lattice disregistry by the assignment of potential values at various surface sites. The surface adatoms are periodically assigned a random energy from a Boltzmann distribution and this energy determines whether the adatoms evaporate, migrate or remain stationary during the sampling interval. For each addition or migration of an adatom, the surface potentials are adjusted to reflect the adsorption, migration or desorption potential changes.

Heinbockel, J. H.↗
DOE OSTI2020

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

Intrinsic Thermal Hall Effect in Mott Insulators

Here, in light of recent experimental data indicating a substantial thermal Hall effect in square lattice antiferromagnetic Mott insulators, we investigate whether a simple Mott insulator can sustain a finite thermal Hall effect. We verify that the answer is “no” if one performs calculations within a spin-only low-energy effective spin model with noninteracting magnons. However, by performing determinant quantum Monte Carlo simulations, we show the single-band 𝑡−𝑡′−𝑈 Hubbard model coupled to an orbital magnetic field does support a finite thermal Hall effect when 𝑡′ ≠0 and 𝐵 ≠0 in the Mott insulating phase. We argue that the (carrier agnostic) necessary conditions for observing a finite thermal Hall effect are time-reversal and particle-hole symmetry breaking. By considering magnon-magnon scattering using a semiclassical Boltzmann analysis, we illustrate a physical mechanism by which finite transverse thermal conductivity may arise, consistent with our symmetry argument and numerical results. Our results indicate that square and triangular lattices with SU(2) symmetry can support a finite thermal Hall effect and call for a critical reexamination of thermal Hall effect data in insulating magnets, as the magnon contribution should not be excluded a priori .

Mott insulators↗
DOE OSTI2023

Breaking Rayleigh’s Law with Spatially Correlated Disorder to Control Phonon Transport

Controlling thermal transport in insulators and semiconductors is crucial for many technological fields such as thermoelectrics and thermal insulation, for which a low thermal conductivity (κ) is desirable. A major obstacle for realizing low κ materials is Rayleigh’s law, which implies that acoustic phonons, which carry most of the heat, are insensitive to scattering by point defects at low energy. We demonstrate, with large scale simulations on tens of millions of atoms, that isotropic long-range spatial correlations in the defect distribution can dramatically reduce phonon lifetimes of important low-frequency heat-carrying modes, leading to a large reduction of κ—potentially an order of magnitude at room temperature. We propose a general and quantitative framework for controlling thermal transport in complex functional materials through structural spatial correlations, and we establish the optimal functional form of spatial correlations that minimize κ. Finally, we end by briefly discussing experimental realizations of various correlated structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Thermoelectric Figure of Merit of a Superatomic Crystal Re 6 Se 8 I 2 Monolayer

Superatomic materials are newly emerging candidates for high-performance thermoelectric (TE) devices due to their intrinsic ultralow thermal conductivities. However, the low TE power factor becomes a huge obstacle to reaching the required the dimensionless figure of merit (ZT) values for practical applications. Here, motivated by the recently synthesized superatomic Re 6 Se 8 I 2 monolayer [J. Am. Chem. Soc. 144, 74 (2022)], we study its superior thermoelectric properties by using density functional theory combined with phonon Boltzmann transport theory and deep potential molecular dynamics. We show that the large mass and anharmonic Re-I bonds introduce strong phonon scattering and result in low lattice thermal conductivity of 1.20 W m -1 K -1 at 300 K, while the strong and harmonic Re-Se network ensures the high TE power factor of 4344 μW m -1 K -2 in b direction for n-type doping. This is an order of magnitude higher than that of other cluster-based materials. Due to the low thermal conductivity and high TE power factor, Re 6 Se 8 I 2 exhibits high ZT values of 1.20 (500 K) and 1.43 (900 K) with n-type doping along the b direction among all the cluster-based thermoelectric materials reported so far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

First principles thermal transport modeling in GaN and related materials

Gallium nitride is a wide bandgap material utilized in a variety of technologies, including high-power electronics and light-emitting diodes, partly due to its favorable thermal properties. This chapter describes modern first-principles-based modeling of phonons and lattice thermal conductivity (k) of GaN, III-nitrides and related materials. In particular, we describe the theoretical underpinnings of calculating phonon dispersions, intrinsic phonon interactions, and other lattice dynamical properties from quantum perturbation theory and density functional theory (DFT) methods. Description of how these methods are then coupled with the Peierls-Boltzmann transport (PBT) equation to determine phonon distributions and lifetimes relevant for thermal transport is given. These theoretical and numerical methods have demonstrated quantitative accuracy and predictive power for calculating pristine k and defect-limited k from first principles for a variety of materials. We present a review of the literature utilizing DFT-PBT methods to understand novel k behaviors in III-nitrides and related materials.

Lindsay, Lucas↗
DOE OSTI2023

Wide-range continuous tuning of the thermal conductivity of La0.5Sr0.5CoO3-δ films via room-temperature ion-gel gating

Abstract Solid-state control of the thermal conductivity of materials is of exceptional interest for novel devices such as thermal diodes and switches. Here, we demonstrate the ability to continuously tune the thermal conductivity of nanoscale films of La 0.5 Sr 0.5 CoO 3- δ (LSCO) by a factor of over 5, via a room-temperature electrolyte-gate-induced non-volatile topotactic phase transformation from perovskite (with δ ≈ 0.1) to an oxygen-vacancy-ordered brownmillerite phase (with δ = 0.5), accompanied by a metal-insulator transition. Combining time-domain thermoreflectance and electronic transport measurements, model analyses based on molecular dynamics and Boltzmann transport equation, and structural characterization by X-ray diffraction, we uncover and deconvolve the effects of these transitions on heat carriers, including electrons and lattice vibrations. The wide-range continuous tunability of LSCO thermal conductivity enabled by low-voltage (below 4 V) room-temperature electrolyte gating opens the door to non-volatile dynamic control of thermal transport in perovskite-based functional materials, for thermal regulation and management in device applications.

42 ENGINEERING↗
DOE OSTI2022

Perturbation theory and thermal transport in mass-disordered alloys: Insights from Green's function methods

Lowest-order quantum perturbation theory (Fermi's golden rule) for phonon-disorder scattering has been used to predict thermal conductivities in several semiconducting alloys with surprising success given its underlying hypothesis of weak and dilute disorder. In this paper, we explain how this is possible by focusing on the case of maximally mass-disordered Mg 2 Si 1-x Sn x . In this work, we use a Chebyshev polynomials Green's function method, which allows a full treatment of disorder on very large systems (tens of millions of atoms) to probe individual phonon linewidths and frequency-resolved thermal transport. We demonstrate that the success of perturbation theory originates from the specific form of mass disorder terms in the phonon Green's function and from the interplay between anharmonic and disorder scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phonon-informed Neural Thermal Scattering (NeTS) Optimization for Crystalline Graphite and Beryllium Metal

Fast neutrons born from fission lose energy through scattering interactions in the process of slowing-down. As neutrons thermalize to the order of $k$ $b$ $T$ (where $k$ $b$ is the Boltzmann constant, and $T$ is the temperature of the medium), their de Broglie wavelength and energy approaches the order of inter-atomic spacing and quantized lattice vibrations, i.e., phonons. At thermal energies, the thermal scattering law (TSL), i.e., $S$($α, β$), captures crystal binding contributions to the total reaction rate, or cross section. This dimensionless material property describes the energy ($β$) and momentum ($α$) exchanges available in a medium. Currently, $S$($α, β$) is evaluated in the Full Law Analysis Scattering System Hub (FLASSH) code for discrete inputs and stored as ENDF/B File 7 for 0-phonon elastic (MT 2) and n-phonon inelastic (MT 4) processes. Further processing recasts $S$($α, β$) into cumulative distribution functions for sampling post-collision scattering kinematics. In practice, interpolation schemes are employed to access data between tabulated values. An improvement to this juncture of the nuclear data pipeline is supplying cross sections on-the-fly (OTF), as has been developed for the un-resolved resonance region to minimize non-physical interpolation errors. This capability may improve simulation accuracy for accident and transient analyses, where rapidly varying changes in temperature and pressure are difficult to predict beforehand. To do so, deep artificial neural networks (ANNs) can be employed which collapse non-linear, complex data into a lightweight dictionary of neural weights and biases. This has been successfully demonstrated for the hydrogen in light water $S$($α, β$) dataset in the form of a Neural Thermal Scattering (NeTS) module. In this work, the NeTS framework is extended to consider the impact of material-dependent dynamical features on optimal neural pre-processing and architecture design decisions, such as number of neurons per hidden layer, residual skip connections and neural depth. New NeTS modules for crystalline graphite and beryllium metal illuminate a novel correlation between dynamical nonlinearity and optimal neural parametrization when deploying $S$($α, β$) on-the-fly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Microscopic mechanism of unusual lattice thermal transport in TlInTe 2

We investigate the microscopic mechanism of ultralow lattice thermal conductivity (κ l ) of TlInTe 2 and its weak temperature dependence using a unified theory of lattice heat transport, that considers contributions arising from the particle-like propagation as well as wave-like tunneling of phonons. While we use the Peierls–Boltzmann transport equation (PBTE) to calculate the particle-like contributions (κ l (PBTE)), we explicitly calculate the off-diagonal (OD) components of the heat-flux operator within a first-principles density functional theory framework to determine the contributions (κ l (OD)) arising from the wave-like tunneling of phonons. At each temperature, T, we anharmonically renormalize the phonon frequencies using the self-consistent phonon theory including quartic anharmonicity, and utilize them to calculate κ l (PBTE) and κ l (OD). With the combined inclusion of κ l (PBTE), κ l (OD), and additional grain-boundary scatterings, our calculations successfully reproduce the experimental results. Our analysis shows that large quartic anharmonicity of TlInTe 2 (a) strongly hardens the low-energy phonon branches, (b) diminishes the three-phonon scattering processes at finite T, and (c) recovers the weaker than T –1 decay of the measured κ l .

36 MATERIALS SCIENCE↗