Exploring the cost implications of increased renewable energy for the U.S. power system
Not Available
SEARCH · Engineering Papers
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.
Not Available
Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005–1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2–3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.
For this work, we calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r 2 ) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee–Yang–Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee–Yang–Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.
The dipole moment is the molecular property that most directly indicates molecular polarity. The accuracy of computed dipole moments depends strongly on the quality of the calculated electron density, and the breakdown of single-reference methods for strongly correlated systems can lead to poor predictions of the dipole moments in those cases. Here, we derive the analytical expression for obtaining the electric dipole moment by multiconfiguration pair density functional theory (MCPDFT), and we assess the accuracy of MC-PDFT for predicting dipole moments at equilibrium and nonequilibrium geometries. We show that MC-PDFT dipole moment curves have reasonable behavior even for stretched geometries, and they significantly improve upon the CASSCF results by capturing more electron correlation. Herein, the analysis of a dataset consisting of 18 first-row transition metal diatomics and 6 main-group polyatomic molecules with multireference character suggests that MC-PDFT and its hybrid extension (HMC-PDFT) perform comparably to CASPT2 and MRCISD+Q methods and have a mean unsigned deviation of 0.2–0.3 D with respect to the best available dipole moment reference values. We explored the dependence of the predicted dipole moments upon the choice of the on-top density functional and active space, and we recommend the tPBE and hybrid tPBE0 on-top choices for the functionals combined with the moderate correlated participating orbital scheme for selecting the active space. With these choices, the mean unsigned deviations (in debyes) of the calculated equilibrium dipole moments from the best estimates are 0.77 for CASSCF, 0.29 for MC-PDFT, 0.24 for HMCPDFT, 0.28 for CASPT2, and 0.25 for MRCISD+Q. These results are encouraging because the computational cost of MC-PDFT or HMC-PDFT is largely reduced compared to the CASPT2 and MRCISD+Q methods.
Deep eutectic solvents (DESs) are emerging as environmentally friendly designer solvents for mass transport and heat transfer processes in industrial applications; however, the lack of accurate tools to predict and thus control their viscosities under both a range of environmental factors and formulations hinders their general application. While DESs may serve as designer solvents, with nearly unlimited combinations, this unfortunately makes it experimentally infeasible to comprehensively measure the viscosities of all DESs of potential industrial interest. To assist in the design of DESs, we have developed several new machine learning (ML) models that accurately and rapidly predict the viscosities of a diverse group of DESs at different temperatures and molar ratios using, to date, one of the most comprehensive data sets containing the properties of over 670 DESs over a wide range of temperatures (278.15–385.25 K). Three ML models, including support vector regression (SVR), feed forward neural networks (FFNNs), and categorical boosting (CatBoost), were developed to predict DES viscosity as a function of temperature and molar ratio and contrasted with multilinear and two-factor polynomial regression baselines. Further, quantum chemistry-based, COSMO-RS-derived sigma profile (σ-profile) features were used as inputs for the ML models. The CatBoost model is excellent at externally predicting DES viscosity, as indicated by high R 2 (0.99) and low root-mean-square-error (RMSE) and average absolute relative deviations (AARD) (5.22%) values for the testing data sets, and 98% of the data points lie within the 15% of AARD deviations. Furthermore, SHapley additive explanation (SHAP) analysis was employed to interpret the ML results and rationalize the viscosity predictions. The result is an ML approach that accurately predicts viscosity and will aid in accelerating the design of appropriate DESs for industrial applications.
Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.
In this work, we developed a chemical synthesis protocol to prepare both spherical and rod shaped YbH x nanoparticles and infiltrate a mesoporous CMK-3 structure with YbH x nanoparticles. The phase evolution of the YbH x nanoparticles was monitored as a function of the pressure in a diamond anvil cell. The cell volume of the Yb atom showed an unexpected positive deviation from the calculated equation of state (EoS) for the α-phase (space group Pnma) of YbH x . Additionally, the β-phase (space group Fm3̅m) YbH x nanoparticles appeared to become amorphous at approximately 27 GPa, although the same behavior was not observed in the Pnma (α-phase). The Pmna phase remained stable (measured up to 59.3 GPa) and showed a significantly higher than expected volume per Yb atom relative to the calculated EoS. This deviation could be attributed to the nanoscale effect on the thermodynamics of ytterbium hydride formation, which could provide a path forward to accessing higher order ytterbium hydrides (YbH x , x > 3) at pressures much lower than in the bulk metal.
This study systematically evaluates a vanadium–niobium (V 94.1 Nb 5.9 ) null-matrix alloy as a reference material for neutron spectrum normalization and compares its performance with that of pure vanadium under identical experimental conditions. Neutron diffraction experiments are conducted on the VULCAN Engineering Materials Diffractometer at the Spallation Neutron Source, Oak Ridge National Laboratory, over a temperature range from room temperature to 1200 °C under vacuum. Pure vanadium exhibited distinct Bragg peaks across all temperatures, with its diffraction behavior influenced by both sample orientation and temperature. As the temperature increased, the diffraction peaks shifted to larger d-spacings and decreased in intensity, while spectral deviation near d ≈ 2.8 Å exceeded 10% at 1200 °C. In contrast, the V–Nb alloy produced a nearly featureless spectrum over the full d-spacing range, confirming near-complete cancellation of coherent scattering over the wide temperature range. Its spectra were insensitive to sample orientation, temperature, and microstructural evolution, with spectral deviation around d ≈ 2.8 Å exceeded 5% at 1200 °C. In conclusion, these results demonstrate that the V–Nb null-matrix alloy provides a thermally stable, efficient, and reliable normalization standard for time-of-flight diffractometers or other instrument where it is needed, enabling reduced data acquisition time and improved data quality in high-temperature neutron diffraction experiments.
Volatile organic compounds (VOCs) play important roles throughout the atmosphere, many of which are altitude dependent. This highlights the need for easily deployable devices to sample VOCs across different atmospheric layers. To address this, we present the design and initial application of a Time Resolved Automated Volatile organIc compounds Sampling system (TRAVIS). VOCs are collected on sorbent tubes, which are subsequently analyzed by a thermal desorption gas chromatography mass spectrometry pipeline. TRAVIS leverages a piezoelectric pump with an integrated pressure sensor for precise (0.1% flow rate relative standard deviation) and accurate (−3 ± 2% error in VOC quantitation) measurements. Via deployment on a tethered balloon system over an agricultural area, TRAVIS is used to show consistent vertically resolved VOC concentrations in a well-mixed (i.e., turbulent) atmosphere (e.g., 5% relative standard deviation for isoprene) and vertically dependent concentrations for a stratified atmosphere (e.g., prior to boundary layer development). Furthermore, we also show VOC information from an intermittent plume via both targeted and untargeted analysis, highlighting future applications for spurious events in agriculture, air quality monitoring, and environmental impact. Overall, the development of TRAVIS represents a lightweight, accurate, sensitive, and precise VOC sampling module for the scientific community.
Magnetite nanocrystals show promise for electrically small gigahertz frequency applications, which could lead to miniaturizing transformer cores and new sensing technologies. This work presents a rigorous radiofrequency characterization of these nanocrystals using vector network analyzer (VNA) ferromagnetic resonance (FMR) measurements. For the first time, two different average diameters of Fe 3 O 4 nanocrystals are investigated (7.3 and 20.2 nm). When the VNA–FMR results were compared to micromagnetic simulations, the magnetic anisotropy (K 1 ) deviated from the ideal mean orientation value of K 1 /2 to K 1 /80 for the 7.3 nm nanocrystal. In contrast, the obtained magnetic anisotropy in 20.2 nm nanocrystals slightly deviated to K 1 /11 due to less structural deformations. These findings resulted in a newly proposed methodology for an approximate simulation based on VNA–FMR measurements. In addition, this work estimates the approximate amount of nanocrystals needed to measure a useful VNA–FMR spectrum.
Topography exerts major control on land surface processes. To improve representation of topographic impacts on land surface processes, a new topography-based subgrid structure has been introduced to the Energy Exascale Earth System Model to represent the subgrid heterogeneity of surface elevation. To take advantage of the new subgrid structure for improving land surface modeling, this study explores four variations of topography-based methods for downscaling grid precipitation to the corresponding subgrids. In the first three methods, the deviation of the subgrid precipitation from the grid’s is equal to the grid precipitation multiplied by the ratio of the elevation difference between the subgrid and grid mean to a specified elevation equals to the grid elevation, the difference between the maximum and minimum subgrid elevation, and the maximum subgrid elevation, respectively. The second method limits the ratio to 0.5 to avoid extreme values on mountains and the third method accounts for the slope effect. The fourth method is similar to the third method except that the Froude Number is used to limit the numerator in a blocking regime of the ambient flow. The downscaled precipitation is evaluated using the PRISM precipitation data over the U.S. using statistical metrics. Results show that by accounting for topographic slope besides elevation, the third and fourth methods show clear advantages over the first and second methods. Furthermore, introducing the Froude Number in the fourth method improves downscaling skill and shows consistent advantages over the third method in areas with larger subgrid heterogeneity across different grids sizes.
Hydrological parameters are scale dependent. Efficient monitoring techniques capable of measuring hydrological parameters, such as soil moisture content (θ), over a wide range of spatial scales are essential for understanding the complexity of water and energy movement across the landscape. Techniques to measure θ over spatial scales in the range from centimeters to thousands of meters, however, are sorely lacking. Recent improvements in the distributed temperature sensing (DTS) technology supported the development of novel techniques to fill that gap. However, improvements in the accuracy and applicability of DTS techniques are still needed. This study investigates the possibility of improving the accuracy of the fiber optics dual-probe heat-pulse (FO-DPHP) DTS technique by using a new design to maintain the spacing between the FO-DPHP probes and by introducing a novel data interpretation approach. The accuracy of the novel FO-DPHP design was tested at different θ in a sand column experiment. The FO-DPHP measurements obtained using traditional and novel data interpretation approaches were compared against independent measurements from several calibrated soil water content (EC5) sensors. Monte-Carlo analyses were also performed to assess the impact of DTS measurement errors on the accuracy achieved using the data interpretation approaches. The novel design and data interpretation approach allowed for accurate measurements of soil thermal properties and θ without the need to perform a hard-to-achieve soil-specific calibration. Measured θ had mean errors and standard deviations <0.03 and <0.01 m 3 m -3 , respectively, for moisture conditions ranging from dry to near saturation. The standard deviation in the measured heat capacity was <0.01 MJ m -3 K -1 .
Abstract Spectral width of the cloud droplet spectrum is important for radiative properties and drizzle/rain development in warm ice‐free clouds. We use an adiabatic rising parcel model to study activation and diffusional growth of cloud droplets, focusing on the spectral width evolution, and contrasting clean and polluted environments. A comprehensive droplet growth equation is used that includes kinetic, solute, and surface tension effects. We show that those effects have an appreciable impact on the spectral width evolution above the cloud base. Without those effects, the droplet area standard deviation should not change once activation is completed. In contrast, simulation results show that the area standard deviation does increase with height, especially for weak updrafts and polluted environments. Implications of those results for cloud modeling, especially applying conventional bin microphysics, are discussed.
Abstract Relative permeability is a key parameter for characterizing the multiphase flow dynamics in porous media at macroscopic scale while it can be significantly impacted by wettability. Recently, it has been reported in microfluidic experiments that wettability is dependent on the pore size (Van Rooijen et al., 2022). To investigate the effect of pore‐size‐dependent wettability on relative permeability, we propose a theoretical framework informed by digital core samples to quantify the deviation of relative permeability curves due to wettability change. We find that the significance of impact is highly dependent on two factors: (i) the function between contact angle and pore size (ii) overall pore size distribution. Under linear function, this impact can be significant for tight porous media with a maximum deviation of 1,000%.
The measured variations in the turbulent static pressure structure function D pp (r) with scale $r$ in the roughness sublayer above a subarctic forest are empirically shown to exhibit exponents that are smaller than $r$ 4/3 predicted for the inertial subrange (ISR). Three hypotheses are offered to explain these deviations. The first is based on conventional intermittency correction to the averaged turbulent kinetic energy dissipation rate, the second is based on shearing introducing deviations from locally isotropic state that must be sensed by both velocity and pressure structure functions, and the third is based on large and inertial scale pressure interactions that persist at values of $r$ within the resolvable ISR. The third hypothesis is shown to yield superior results, which allows a new formulation for D pp (r) to be derived that accommodates such finite interactions.
Relative permeability is a crucial two-phase property in porous media that can be significantly impacted by wettability conditions. While traditional research has predominantly examined homogeneous wettability, this work explores the less studied pore-size dependent (PSD) wettability, featured by a pore-size dependent wettability distribution. Leveraging high-fidelity Lattice Boltzmann simulations on CT-scanned porous samples, we demonstrate how PSD wettability would impact relative permeability at the pore scale. Our findings reveal that the deviation of relative permeability from the homogeneous wettabilitity induced by PSD wettability can be 5%–20%. The deviation of relative permeability curves increases as the spanning range of the contact angle increases. We also find that this impact is less pronounced as the capillary number increases. By adopting a pore-size-dependent contact angle relationship, our approach provides a more accurate and nuanced understanding of how PSD wettability would impact two-phase flow. These findings discovered at the pore scale may also provide valuable insights on relative permeability at the core to reservoir scales.
Design and processing of advanced lightweight structural alloys based on magnesium and titanium rely critically on a control over twinning which, however, remains elusive to date and is dependent on an explicit understanding on the twinning nucleation mechanism in hexagonal close-packed (HCP) crystals. Here, by using in-situ high resolution transmission electron microscopy, we directly reveal a “dual-step” twinning nucleation mechanism in HCP rhenium nanocrystals. It is found that nucleation of the predominant {10 -12} twinning is initiated by disconnections on the Prismatic¦Basal interfaces which establish the lattice correspondence of the twin with a minor deviation from the ideal orientation. Subsequently, the minor deviation is corrected by the formation of coherent twin boundaries through rearrangement of the disconnections on the Prismatic¦Basal interface; thereafter, the coherent twin boundaries propagate by twinning disconnections. The findings provide much-needed high-resolution evidences to the twinning nucleation mechanism in HCP crystals, with crucial implications for engineering advanced structural alloys.
There are now two single measurements of precision observables that have major anomalies in the Standard Model: the recent CDF measurement of the W mass shows a 7σ deviation and the Muon g – 2 experiment at FNAL confirmed a long-standing anomaly, implying a 4.2σ deviation. Doubts regarding new physics interpretations of these anomalies could stem from uncertainties in the common hadronic contributions. We demonstrate that these two anomalies pull the hadronic contributions in opposite directions by performing electroweak fits in which the hadronic contribution was allowed to float. The fits show that including the g – 2 measurement worsens the tension with the CDF measurement and conversely that adjustments that alleviate the CDF tension worsen the g – 2 tension beyond 5σ. This means that if we adopt the CDF W mass measurement, the case for new physics in either the W mass or muon g – 2 is inescapable regardless of the size of the SM hadronic contributions. Lastly, we demonstrate that a mixed scalar leptoquark extension of the Standard Model could explain both anomalies simultaneously.