Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit correlation”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Multifaceted Water Dynamics in Spherical Nanocages

We present here a new method to study position-dependent, anisotropic diffusion tensors inside spherically confined systems—a geometry that is common to many chemical nanoreactors. We use this method to elucidate the surprisingly rich solvent dynamics of confined water. The spatial variation of the strongly anisotropic diffusion predicted by the model agrees with the results of explicit molecular dynamics simulations. The same approach can be directly transferred to the transport of solutes to and from reaction sites located at nanoreactor interfaces. We complement our study by a detailed analysis of water hydrogen bond kinetics, which is intimately coupled to diffusion. Despite the inhomogeneity in structure and translational dynamics inside our nanocages, a single set of well-defined rate constants is sufficient to accurately describe the kinetics of hydrogen bond breaking and formation. We find that once system size effects have been eliminated, the residence times of water molecules inside the coordination shell of a hydrogen bond partner are well correlated to average diffusion constants obtained from the procedure above.

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

Biases to primordial non-Gaussianity measurements from CMB secondary anisotropies

ABSTRACT Our view of the last-scattering surface in the cosmic microwave background (CMB) is obscured by secondary anisotropies, sourced by scattering, extragalactic emission, and gravitational processes between recombination and observation. Whilst it is established that non-Gaussianity from the correlation between the integrated-Sachs–Wolfe (ISW) effect and gravitational lensing can significantly bias primordial non-Gaussianity (PNG) searches, recent work by Hill suggests that other combinations of secondary anisotropies can also produce significant biases. Building on that work, we use the WebSky and Sehgal et al. simulations to perform an extensive examination of possible biases to PNG measurements for the local, equilateral and orthogonal shapes. For a Planck-like CMB experiment, without foreground cleaning, we find significant biases from cosmic infrared background (CIB)-lensing and thermal Sunyaev–Zel’dovich (tSZ)-lensing bispectra for the local and orthogonal templates, and from CIB and tSZ bispectra for the equilateral template. For future experiments, such as the Simons Observatory, biases from correlations between the ISW effect and the tSZ and CIB will also become important. Finally, we investigate the effectiveness of foreground-cleaning techniques to suppress these biases. We find that the majority of these biases are effectively suppressed by the internal-linear combination method with a total bias below the $1\, \sigma$ statistical error for both experiments. However, the small total bias arises from the cancellation of several $1\, \sigma$ biases for Planck-like experiments and $2\, \sigma$ biases for SO-like. As this cancellation is likely sensitive to the modelling, to ensure robustness against these biases, we recommend that explicit removal methods should be used.

Astronomy & Astrophysics↗
DOE OSTI2016

Estimation of heat transfer coefficients for biomass particles by direct numerical simulation using microstructured particle models in the Laminar regime

Here, direct numerical simulation of convective heat transfer from hot gas to isolated biomass particle models with realistic morphology and explicit microstructure was performed over a range of conditions with laminar flow of hot gas (500 degrees C). Steady-state results demonstrated that convective interfacial heat transfer is dependent on the wood species. The computed heat transfer coefficients were shown to vary between the pine and aspen models by nearly 20%. These differences are attributed to the species-specific variations in the exterior surface morphology of the biomass particles. We also quantify variations in heat transfer experienced by the particle when positioned in different orientations with respect to the direction of fluid flow. These results are compared to previously reported heat transfer coefficient correlations in the range of 0.1 < Pr < 1.5 and 10 < Re < 500. Comparison of these simulation results to correlations commonly used in the literature (Gunn, Ranz-Marshall, and Bird-Stewart-Lightfoot) shows that the Ranz-Marshall (sphere) correlation gave the closest h values to our steady-state simulations for both wood species, though no existing correlation was within 20% of both species at all conditions studied. In general, this work exemplifies the fact that all biomass feedstocks are not created equal, and that their species-specific characteristics must be appreciated in order to facilitate accurate simulations of conversion processes.

09 BIOMASS FUELS↗
DOE OSTI2014

A stochastic extension of the explicit algebraic subgrid-scale models

The explicit algebraic subgrid-scale (SGS) stress model (EASM) of Marstorp et al. [“Explicit algebraic subgrid stress models with application to rotating channel flow,” J. Fluid Mech. 639, 403–432 (2009)] and explicit algebraic SGS scalar flux model (EASFM) of Rasam et al. [“An explicit algebraic model for the subgrid-scale passive scalar flux,” J. Fluid Mech. 721, 541–577 (2013)] are extended with stochastic terms based on the Langevin equation formalism for the subgrid-scales by Marstorp et al. [“A stochastic subgrid model with application to turbulent flow and scalar mixing,” Phys. Fluids 19, 035107 (2007)]. The EASM and EASFM are nonlinear mixed and tensor eddy-diffusivity models, which improve large eddy simulation (LES) predictions of the mean flow, Reynolds stresses, and scalar fluxes of wall-bounded flows compared to isotropic eddy-viscosity and eddy-diffusivity SGS models, especially at coarse resolutions. The purpose of the stochastic extension of the explicit algebraic SGS models is to further improve the characteristics of the kinetic energy and scalar variance SGS dissipation, which are key quantities that govern the small-scale mixing and dispersion dynamics. LES of turbulent channel flow with passive scalar transport shows that the stochastic terms enhance SGS dissipation statistics such as length scale, variance, and probability density functions and introduce a significant amount of backscatter of energy from the subgrid to the resolved scales without causing numerical stability problems. The improvements in the SGS dissipation predictions in turn enhances the predicted resolved statistics such as the mean scalar, scalar fluxes, Reynolds stresses, and correlation lengths. Moreover, the nonalignment between the SGS stress and resolved strain-rate tensors predicted by the EASM with stochastic extension is in much closer agreement with direct numerical simulation data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Solvation of furfural at metal–water interfaces: Implications for aqueous phase hydrogenation reactions

Metal–water interfaces are central to understanding aqueous-phase heterogeneous catalytic processes. However, the explicit modeling of the interface is still challenging as it necessitates extensive sampling of the interfaces’ degrees of freedom. Herein, we use ab initio molecular dynamics (AIMD) simulations to study the adsorption of furfural, a platform biomass chemical on several catalytically relevant metal–water interfaces (Pt, Rh, Pd, Cu, and Au) at low coverages. We find that furfural adsorption is destabilized on all the metal–water interfaces compared to the metal–gas interfaces considered in this work. This destabilization is a result of the energetic penalty associated with the displacement of water molecules near the surface upon adsorption of furfural, further evidenced by a linear correlation between solvation energy and the change in surface water coverage. To predict solvation energies without the need for computationally expensive AIMD simulations, we demonstrate OH binding energy as a good descriptor to estimate the solvation energies of furfural. Using microkinetic modeling, we further explain the origin of the activity for furfural hydrogenation on intrinsically strong-binding metals under aqueous conditions, i.e., the endothermic solvation energies for furfural adsorption prevent surface poisoning. In conclusion, our work sheds light on the development of active aqueous-phase catalytic systems via rationally tuning the solvation energies of reaction intermediates.

aqueous-phase hydrogenation↗
DOE OSTI2024

Dark matter from anomaly cancellation at the LHC

We discuss a class of theories that predict a fermionic dark matter candidate from gauge anomaly cancellation. As an explicit example, we study the predictions in theories where the global symmetry associated with baryon number is promoted to a local gauge symmetry. In this context, the symmetry-breaking scale has to be below the multi-TeV scale in order to be in agreement with the cosmological constraints on the dark matter relic density. The new physical “Cucuyo” Higgs boson in the theory has very interesting properties, decaying mainly into two photons in the low mass region, and mainly into dark matter in the intermediate mass region. We study the most important signatures at the Large Hadron Collider, evaluating the experimental bounds. We discuss the correlation between the dark matter relic density, direct detection, and collider constraints. We find that these theories are still viable and are susceptible to being probed in current, and future high-luminosity, running. Published by the American Physical Society 2024

Butterworth, Jon (ORCID:0000000259055394)↗
DOE OSTI2020

Cumulants of multiple conserved charges and global conservation laws

We analyze the behavior of cumulants of conserved charges in a subvolume of a thermal system with exact global conservation laws by extending a recently developed subensemble acceptance method (SAM) to multiple conserved charges. Explicit expressions for all diagonal and off-diagonal cumulants up to sixth order that relate them to the grand canonical susceptibilities are obtained. The derivation is presented for an arbitrary equation of state with an arbitrary number of different conserved charges. The global conservation effects cancel out in any ratio of two second order cumulants, in any ratio of two third order cumulants, as well as in a ratio of strongly intensive measures Σ and Δ involving any two conserved charges, making all these quantities particularly suitable for theory-to-experiment comparisons in heavy-ion collisions. We also show that the same cancellation occurs in correlators of a conserved charge, like the electric charge, with any non-conserved quantity such as net proton or net kaon number. The main results of the SAM are illustrated in the framework of the hadron resonance gas model. We also elucidate how net-proton and net-Λ fluctuations are affected by conservation of electric charge and strangeness in addition to baryon number.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Avoiding ergodicity problems in lattice discretizations of the Hubbard model

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. Here, we show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

First-principles calculations↗
DOE OSTI2021

Revisiting Understanding of Electrochemical CO 2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory

Copper (Cu) electrodes, as the most efficacious of CO 2 reduction reaction (CO 2 RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO 2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO 2 RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO 2 RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO 2 RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at –0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Furthermore, our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

High strain-rate compression behavior of polymeric rod and plate Kelvin lattice structures

The compressive high strain-rate behavior of polymeric Kelvin lattice structures with rod-based or plate-based unit cells was investigated through experimental techniques and finite element simulations. Polymeric lattice structures with 5x5x5 unit cell geometries were manufactured on the millimeter scale using vat polymerization additive manufacturing and tested at low (0.001/s) and high (1000/s) strain-rates. High strain-rate experiments were performed and validated for a viscoelastic split-Hopkinson (Kolsky) pressure bar system (SHPB) coupled with high-speed imaging and digital image correlation (DIC). Experimental results at both low and high strain-rates show the formation of a localized deformation band which was more prevalent in low relative density specimens and low strain-rate experiments. Strain-rate effects of lattice specimens strongly correlate with effects of the base polymer material; both bulk polymer and lattice specimen demonstrated strain-rate hardening, strain-rate stiffening, and decreased fracture strain under dynamic loading. Results show mechanical failure properties and energy absorption depended strongly on the relative density of the lattice specimen and exhibited distinct scaling between relative density and geometry type (rod, plate) and loading rate. High relative density plate-lattices demonstrated inferior mechanical properties to rod-lattices; however, there exists a critical relative density for a given mechanical property (17%- 28%) below which plate-lattices outperform rod-lattices of similar mass. As a result, high strain-rate explicit finite element simulations were performed and showed good agreement with the mechanical failure trends and deformation modes observed in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photometric Analyses of Saturn’s Small Moons: Aegaeon, Methone, and Pallene Are Dark; Helene and Calypso Are Bright

We examine the surface brightnesses of Saturn’s smaller satellites using a photometric model that explicitly accounts for their elongated shapes and thus facilitates comparisons among different moons. Analyses of Cassini imaging data with this model reveal that the moons Aegaeon, Methone, and Pallene are darker than one would expect given trends previously observed among the nearby mid-sized satellites. On the other hand, the trojan moons Calypso and Helene have substantially brighter surfaces than their co-orbital companions Tethys and Dione. These observations are inconsistent with the moons’ surface brightnesses being entirely controlled by the local flux of E-ring particles, and therefore strongly imply that other phenomena are affecting their surface properties. The darkness of Aegaeon, Methone, and Pallene is correlated with the fluxes of high-energy protons, implying that high-energy radiation is responsible for darkening these small moons. Meanwhile, Prometheus and Pandora appear to be brightened by their interactions with the nearby dusty F ring, implying that enhanced dust fluxes are most likely responsible for Calypso’s and Helene’s excess brightness. However, there are no obvious structures in the E ring that would preferentially brighten these two moons, so there must either be something subtle in the E-ring particles’ orbital properties that leads to asymmetries in the relevant fluxes, or something happened recently to temporarily increase these moons’ brightnesses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Origin of parameter degeneracy and molecular shape relationships in geometric-flow calculations of solvation free energies

Implicit solvent models are important tools for calculating solvation free energies for chemical and biophysical studies since they require fewer computational resources but can achieve accuracy comparable to that of explicit-solvent models. In past papers, geometric flow-based solvation models have been established for solvation analysis of small and large compounds. In the present work, the use of realistic experiment-based parameter choices for the geometric flow models is studied. We find that the experimental parameters of solvent internal pressure p = 172 MPa and surface tension γ = 72 mN/m produce solvation free energies within 1 RT of the global minimum root-mean-squared deviation from experimental data over the expanded set. Our results demonstrate that experimental values can be used for geometric flow solvent model parameters, thus eliminating the need for additional parameterization. We also examine the correlations between optimal values of p and γ which are strongly anti-correlated. Geometric analysis of the small molecule test set shows that these results are inter-connected with an approximately linear relationship between area and volume in the range of molecular sizes spanned by the data set. In spite of this considerable degeneracy between the surface tension and pressure terms in the model, both terms are important for the broader applicability of the model.

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

Investigating the magnetovolume effect in isotropic body-centered-cubic iron using spin-lattice dynamics simulations

The understanding of the magnetovolume effect lacks explicit consideration of spin-lattice coupling at the atomic level, despite abundant theoretical and experimental studies throughout the years. This research gap is filled by the recently developed spin-lattice dynamics technique implemented in this study, which investigates the magnetovolume effect of isotropic body-centered-cubic (BCC) iron, a topic that has previously been subject to macroscopic analysis only. This approach demonstrates the magnetic anomaly followed by the volumetric changes associated with the effect, each characterized by the corresponding field-induced inflection temperature. The temperature of the heat capacity peaks is useful in determining the temperature for retarding the atomic volume increase. Moreover, this work shows the correlation between the effects of temperature and field strength in determining the equilibrium atomic volume of a ferromagnetic material under a magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Examining the Characteristics of the Cropland Data Layer in the Context of Estimating Land Cover Change

The United States Department of Agriculture (USDA) Cropland Data Layer (CDL) provides spatially explicit information about crop production area and has served as a prevalent data source for characterizing cropland change in the U.S. in the last decade. Understanding the accuracy of the CDL is paramount because of the reliance on it for management and policy making. This study examined the characteristics of the CDL from 2007 to 2017 using comparisons to other USDA datasets. The results showed when examining the cropland area for the same year, the CDL produced comparable trends with other datasets (R 2 > 0.95), but absolute area differed. The estimated area of cropland changes from 2007 to 2012, 2008 to 2012 and 2012 to 2017 varied from weak to moderate correlation between the CDL and the tabular data (R 2 = 0.005~0.63). Differences in area of cropland change varied widely between data sources with the CDL estimating much larger change area. A series of image processing techniques designed to improve the confidence in cropland change estimated using the CDL reduced the area of estimated cropland change. The techniques also, unexpectedly, lowered the correlation in change estimated between the CDL and the tabular datasets. Estimated land cover change area varied widely based on analyses applied and could reverse from increasing to declining area in cropland. Further analyses showed unlikely change scenarios when comparing different year combinations. The authors recommend the CDL only be used for land cover change analysis if the error can be estimated and is within change estimates.

58 GEOSCIENCES↗
DOE OSTI2016

Redshift-space distortions with wide angular separations

Redshift-space distortions are generally considered in the plane parallel limit, where the angular separation between the two sources can be neglected. Given that galaxy catalogues now cover large fractions of the sky, it becomes necessary to consider them in a formalism which takes into account the wide angle separations. In this article we derive an operational formula for the matter correlators in the Newtonian limit to be used in actual data sets. In order to describe the geometrical nature of the wide angle RSD effect on Fourier space, we extend the formalism developed in configuration space to Fourier space without relying on a plane-parallel approximation, but under the extra assumption of no bias evolution. We then recover the plane-parallel limit not only in configuration space where the geometry is simpler, but also in Fourier space, and we exhibit the first corrections that should be included in large surveys as a perturbative expansion over the plane-parallel results. We finally compare our results to existing literature, and show explicitly how they are related.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

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

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗