Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical thermodynamics”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Finite-temperature vibrational full configuration interaction

Thermodynamic functions of an ideal molecular gas due to its anharmonic vibrations are evaluated in a wide range of temperature (T) by the vibrational full-configuration-interaction (FCI) method using a quartic force field and a finite number (N) of harmonic-oscillator basis functions along each normal mode. The thermodynamic functions considered are the grand potential (Ω), internal energy (U), and entropy (S). They are compared with those obtained from the Bose–Einstein theory with or without truncation of the harmonic-oscillator basis functions after quantum number N–1. The comparison reveals that the finite-basis-set errors in Ω and U are, respectively, k B Tln(k B T/Nℏω) and k B T per mode in the high-T limit, obscuring anharmonic effects when k B T > ℏω, where ω is the lowest mode frequency. Here, the benchmark data for several low-order perturbation corrections to Ω, U, and S are also obtained as the numerical derivatives of their FCI values with respect to dimensionless perturbation strength, and the domain of T and N in which these data are reliable (for the N → ∞ limits) is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Bubble Point Measurements of Three Binary Mixtures of Refrigerants: R-32/1234yf, R-32/1234ze(E), and R-1132a/1234yf

In this article, bubble point pressures of three binary mixtures at two compositions each have been measured utilizing a static method. The mixtures studied were difluoromethane + 2,3,3,3-tetrafluoroprop-1-ene (R-32/1234yf), difluoromethane + (E)-1,3,3,3-tetrafluoroprop-1-ene (R-32/1234ze(E)), and 1,1-difluoroethene + 2,3,3,3-tetrafluoroprop-1-ene (R-1132a/1234yf). For all of the mixtures, the minimum temperature measured was 270 K. The maximum temperatures measured for the R-32/1234ze(E), R-32/1234yf, and R-1132a/1234yf mixtures were 360, 345, and 325 K, respectively. A total of 100 bubble point pressures are reported and compared to available literature data relative to the REFPROP database. Lemmon, E. W. [ NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, 10.0; National Institute of Standards and Technology, Standard Reference Data Program: Gaithersburg, 2018]. Additionally, the bubble point pressures obtained in this study are modeled using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state. The PC-SAFT model parameters, as well as deviations of this work and literature data from the model, are reported.

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

Ionization state and dielectric constant in cold rarefied hydrocarbon plasmas of inertial confinement fusion

A combined approach to study cold rarefied matter is introduced that includes a semianalytical method based on the free-energy minimization and ab initio calculations based on the finite-temperature density-functional theory. The approach is used to calculate the ionization state of hydrocarbon (CH) under the shock-release conditions in inertial confinement fusion. Furthermore, the dielectric constant of CH is calculated using the Kubo-Greenwood formulation and contribution from atomic polarizabilities is found to be as important as the free-electron contribution. Using the ionization state and dielectric constant, the electron density profile in the rarefaction wave of the shock-release plasma is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kohn-Sham density functional theory of Abelian anyons

In this work, we develop a density functional treatment of noninteracting Abelian anyons, which is capable, in principle, of dealing with a system of a large number of anyons in an external potential. Comparison with exact results for a few particles shows that the model captures the behavior qualitatively and semiquantitatively, especially in the vicinity of the fermionic statistics. We then study anyons with statistics parameter 1+1/n, which are thought to condense into a superconducting state. An indication of the superconducting behavior is the mean-field result that, for uniform density systems, the ground-state energy increases under the application of an external magnetic field independent of its direction. Our density functional theory based analysis does not find that to be the case for finite systems of anyons, which can accommodate a weak external magnetic field through density transfer between the bulk and the boundary rather than through transitions across effective Landau levels, but the “Meissner repulsion” of the external magnetic field is recovered in the thermodynamic limit as the effect of the boundary becomes negligible. We also consider the quantum Hall effect of anyons, and show that its topological properties, such as the charge and statistics of the excitations and the quantized Hall conductance, arise in a self-consistent fashion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Fluctuating hydrodynamics of chiral active fluids

Active materials are characterized by continuous injection of energy at the microscopic level and typically cannot be adequately described by equilibrium thermodynamics. Here we study a class of active fluids in which equilibrium-like properties emerge when fluctuating and activated degrees of freedom are statistically decoupled, such that their mutual information is negligible. We analyse three paradigmatic systems: chiral active fluids composed of spinning frictional particles that are free to translate, oscillating granular gases and active Brownian rollers. In all of these systems, a single effective temperature generated by activity parameterizes both the equation of state and the emergent Boltzmann statistics. The same effective temperature, renormalized by velocity correlations, relates viscosities to steady-state stress fluctuations via a Green–Kubo relation. To rationalize these observations, we develop a theory for the fluctuating hydrodynamics of these non-equilibrium fluids and validate it through large-scale molecular dynamics simulations. Our work sheds light on the microscopic origin of odd viscosities and stress fluctuations characteristic of parity-violating fluids, in which mirror symmetry and detailed balance are broken.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Understanding the Recent Increase in Landfalling Tropical Cyclones Over Florida's Gulf Coast

Unlike Florida's Atlantic Coast, the Gulf Coast of Florida has seen heightened tropical cyclone (TC) activity in recent decades with several destructive landfalls. Here, we attempt to understand this regional contrast using a suite of observations for the period 1979–2024. First, we demonstrate that while the El Niño Southern Oscillation (ENSO), the Atlantic Multidecadal Oscillation (AMO) and the North Atlantic Oscillation (NAO) can explain ~23% of the interannual variability in landfalls over the Gulf Coast, the variance explained by them for the Atlantic Coast is statistically insignificant. Next, we show that this striking difference may be attributed to the regional patterns of wind shear, steering flow and air-sea thermodynamic state excited by those modes of variability. The differential control exerted by ENSO, AMO, and NAO on landfalling Florida TCs, in combination with decadal trends in those modes, is likely responsible for the observed increases in landfalls over Florida's Gulf Coast.

Florida↗
DOE OSTI2020

Machine learning enabled discovery of application dependent design principles for two-dimensional materials

Abstract The unique electronic and mechanical properties of two-dimensional (2D) materials make them promising next-generation candidates for a variety of applications. Large-scale searches for high-performing 2D materials are limited to calculating descriptors with computationally demanding first-principles density functional theory. In this work, we alleviate this issue by extending and generalizing crystal graph convolutional neural networks to systems with planar periodicity and train an ensemble of models to predict thermodynamic, mechanical and electronic properties. We carry out a screening of nearly 45,000 structures for two separate applications: mechanical strength and photovoltaics. By collecting statistics of the screened candidates, we investigate structural and compositional design principles that impact the properties of the structures surveyed. Our approach recovers some well-accepted design rules: hybrid organic-inorganic perovskites with lead and tin tend to be good candidates for solar cell applications and titanium based MXenes usually have high stiffness coefficients. Interestingly, other members of the group 4 elements also contribute to increasing the mechanical strength of MXenes. For all-inorganic perovskites, we discover some compositions that have not been deeply studied in the field of photovoltaics and thus open up paths for further investigation. We open-source the code-base to spur further development in this space.

Venturi, Victor (ORCID:0000000156727549)↗
DOE OSTI2025

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗
DOE OSTI2021

Interfacial Effects on Transport Coefficient Measurements in Li-ion Battery Electrolytes

Development of Li + -containing electrolytes with improved transport properties requires reliable, reproducible, and ideally low volume techniques to rigorously understand ion-transport with varying composition. Precisely measuring the complete set of transport coefficients in liquid electrolytes under battery-relevant operating conditions is difficult and the reliability of these methods are sparsely described in electrolyte transport literature. In this work, we apply a potentiostatic polarization-based transport characterization approach typically used for polymer electrolytes to liquid electrolyte systems in an attempt to fully measure all transport coefficients (conductivity, total salt diffusion coefficient, thermodynamic factor and transference number) for the model system of LiPF 6 in an ethylene carbonate—ethyl methyl carbonate (EC:EMC) mixture. Using systematic timescale and statistical analyses, we find that transport coefficients measured using potentiostatic polarization of Li-Li symmetric cells exhibit strong correlation to Li electrode interfacial resistance, indicating that such methods are probing both bulk and interfacial phenomena. This reveals a major roadblock in characterizing electrolyte systems where the interfacial resistance is significantly larger than ohmic electrolyte resistance. As a result, we find that methods that rely on potentiostatic Li metal stripping/plating do not readily result in reliable liquid electrolyte transport coefficients, unlike similar methods for solid polymer electrolytes, where interfacial resistances are typically smaller than electrolyte resistances at the elevated temperatures typically of interest for such electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2022

Activated relaxation in supercooled monodisperse atomic and polymeric WCA fluids: Simulation and ECNLE theory

Here, we combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks–Chandler–Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, a c , that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, N c , that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of N c ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.

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

Turbulent burning velocity of lean premixed hydrogen/air flames at engine conditions: Effects of turbulence intensity and length scale

For turbulent lean premixed hydrogen flames with strong thermodiffusively instabilities, most previous studies have focused on the influence of turbulence intensity, whereas the role of turbulence length scale is less well understood. Here, this study addresses this gap by conducting direct numerical simulations (DNS) of statistically planar turbulent premixed flames for a lean (ϕ=0.35) hydrogen/air mixture under independently varied turbulence intensity (u') and length scale (l T ) at engine-relevant thermodynamics conditions. Results show that as u' increases, the flame front becomes increasingly wrinkled, forming smaller cellular structures. In contrast, l T variations do not significantly alter the size of these structures. For the turbulent burning velocity (S T ), the normalized S T (i.e., S T /S L , where S L is the laminar flame speed) increases linearly with u', driven by both enhanced flame surface wrinkling (i.e., increased A T /A L ) and enhanced local burning rate (i.e., increased I 0 ). However, increasing l T reduces I 0 , despite a continued increase in A T /A L , resulting in only a marginal increase in S T /S L . To reveal the underlying mechanisms, especially the decreasing trend of I 0 with l T , local flame dynamics analyses are performed. It is found that as l T increases, the interaction between thermodiffusive effects and turbulence weakens due to the reduced tangential strain rate, while the flame curvature remains largely unchanged. This suppresses local reactivity enhancement and thus decreases I 0 , In contrast, an increase in u' enhances the interaction by amplifying both curvature fluctuation and tangential strain rate, leading to increased local reactivity (increased I 0 ). Finally, based on the DNS data, several new scaling models are proposed for the three global properties, S T /S L , A T /A L , and I 0 , and showed improvements compared to existing models. These findings provide new insights into the flame-turbulence interactions in thermodiffusively unstable hydrogen flames. The DNS dataset is also useful for the development of turbulent combustion models applicable to practical engine simulations.

Engine-relevant condition↗
DOE OSTI2025

Temperature measurement of Quark-Gluon plasma at different stages

In a Quark-Gluon Plasma (QGP), the fundamental building blocks of matter, quarks and gluons, are under extreme conditions of temperature and density. A QGP could exist in the early stages of the Universe, and in various objects and events in the cosmos. The thermodynamic and hydrodynamic properties of the QGP are described by Quantum Chromodynamics (QCD) and can be studied in heavy-ion collisions. Despite being a key thermodynamic parameter, the QGP temperature is still poorly known. Thermal lepton pairs (e + e − and μ + μ − ) are ideal penetrating probes of the true temperature of the emitting source, since their invariant-mass spectra suffer neither from strong final-state interactions nor from blue-shift effects due to rapid expansion. Here we measure the QGP temperature using thermal e+e− production at the Relativistic Heavy Ion Collider (RHIC). The average temperature from the low-mass region (in-medium ρ 0 vector-meson dominant) is (2.01 ± 0.23) × 10 12 K, consistent with the chemical freeze-out temperature from statistical models and the phase transition temperature from Lattice QCD. The average temperature from the intermediate mass region (above the ρ0 mass, QGP dominant) is significantly higher at (3.25 ± 0.60) × 10 12 K. This work provides essential experimental thermodynamic measurements to map out the QCD phase diagram and understand the properties of matter under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

CASM — A software package for first-principles based study of multicomponent crystalline solids

CASM is a software package that enables first-principles based studies of crystalline materials. It has been designed to treat coupled chemical, mechanical, vibrational, and magnetic degrees of freedom to determine ground state and finite temperature properties of crystals. The symmetry of the underlying parent crystal structure is used to enumerate perturbations of the parent crystal structure and generate derivative structures which can be input to first-principles calculations in order to explore the ground state energy landscape. CASM algorithmically constructs cluster expansions that fully couple discrete and continuous degrees of freedom, and generates highly efficient code to evaluate the cluster expansion basis functions. Widely used machine learning methods are integrated for fitting expansion coefficients to first-principle calculations. The fully parameterized cluster expansions can be combined with (kinetic) Monte Carlo methods to calculate finite temperature thermodynamic and kinetic properties. CASM Alloy Manager identifies distinct parent crystal structures in an alloy system, creating individual projects for each, and integrating the results. The integrated infrastructure facilitates the linkage between first-principles statistical mechanics predictions with higher length scale computational methods, such as phase field simulations. In conclusion, CASM projects are designed to be easy to share in repositories, to re-use, and to extend to include additional chemical species or types of degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Role of Subcloud Mesoscale Convergence in Sculpting Convective Updraft Width and Depth

The initiation of deep moist convection is governed in part by the horizontal width of updrafts near cloud base, which limits the deleterious effects of entrainment-driven dilution on buoyant thermals ascending through the free troposphere. However, the factors controlling cloud-base updraft width, which in turn dictates cloud depth, are not well understood. We track the evolving three-dimensional structure of the mesoscale subcloud forcing for vertical motion and near-cloud thermodynamic ingredients within a high-resolution ensemble of simulations of seven realistic daytime orographic convection initiation events to determine their relative roles in controlling cloud width and depth. Statistical analysis of approximately 5000 cloudy updraft samples indicates that the most important contributors to the width of cloudy updrafts across the ensemble are the depth and magnitude of the subcloud mesoscale ascent. However, the depth achieved by clouds is more consistently predicted by the near-cloud ambient relative humidity within the lower to middle free troposphere and convective available potential energy. Therefore, although the width of cloudy updrafts may be partly set at low levels by the mesoscale vertical mass and moisture flux, the likelihood of deep moist convection is governed by the generation of positive buoyancy within cumulus thermals and entrainment-driven dilution that reduces it. The persistence of the low-level mesoscale vertical forcing locally consolidates and vertically transports boundary layer moisture, helping to reduce updraft dilution. However, these factors vary in relative impacts on cloudy updrafts across individual cases, indicating multiple pathways for deep convection initiation.

Convective storms↗
DOE OSTI2024

Statistical mechanics and pressure of composite multimoded weakly nonlinear optical systems

Statistical mechanics can provide a versatile theoretical framework for investigating the collective dynamics of weakly nonlinear-wave settings that can be utterly complex to describe otherwise. In optics, composite systems arise due to interactions between different frequencies and polarizations. The purpose of this work is to develop a thermodynamic theory that takes into account the synergistic action of multiple components. We find that the type of the nonlinearity involved can have important implications in the thermalization process and, hence, can lead to different thermal equilibrium conditions. Importantly, we derive closed-form expressions for the actual optomechanical pressure that is exerted on the system. In particular, the total optomechanical pressure is the sum of the partial pressures due to each component. Our results can be applied to a variety of weakly nonlinear optical settings such as multimode fibers, bulk waveguides, photonic lattices, and coupled microresonators. We present two specific examples, where two colors interact in a one-waveguide array with either a cubic or quadratic nonlinearity.

Efremidis, Nikolaos K. (ORCID:0000000298300268)↗
DOE OSTI2024

Air separation and N 2 purification with Ba 0.15 Sr 0.85 FeO 3-$δ$ via a two-step thermochemical process

Thermochemical air separation to produce high-purity N 2 was demonstrated in a vertical tube reactor via a two-step reduction–oxidation cycle with an A-site substituted perovskite Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585). BSF1585 particles were synthesized and characterized in terms of their chemical, morphological, and thermophysical properties. A thermodynamic cycle model and sensitivity analysis using computational heat and mass transfer models of the reactor were used to select the system operating parameters for a concentrating solar thermal-driven process. Thermal reduction up to 800 °C in air and temperature-swing air separation from 800 °C to minimum temperatures between 400 and 600 °C were performed in the reactor containing a 35 g packed bed of BSF1585. The reactor was characterized for dispersion, and air separation was characterized via mass spectrometry. Gas measurements indicated that the reactor produced N 2 with O 2 impurity concentrations as low as 0.02 % for > 30 min of operation. Additionally, a parametric study of air flow rates suggested that differences in observed and thermodynamically predicted O 2 impurities were due to imperfect gas transport in the bed. Temperature swing reduction/oxidation cycling experiments between 800 and 400 °C in air were conducted with no statistically significant degradation in N 2 purity over 50 cycles.

14 SOLAR ENERGY↗