SEARCH · Engineering Papers
Results for “water model”
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.
DeePKS-kit: A package for developing machine learning-based chemically accurate energy and density functional models
We introduce DeePKS-kit, an open-source software package for developing machine learning based energy and density functional models. DeePKS-kit is interfaced with PyTorch, an open-source machine learning library, and PySCF, an ab initio computational chemistry program that provides simple and customized tools for developing quantum chemistry codes. It supports the DeePHF and DeePKS methods. In addition to explaining the details in the methodology and the software, we also provide an example of developing a chemically accurate model for water clusters.
Using SWMM for emergency response planning: A case study evaluating biological agent transport under various rainfall scenarios and urban surfaces
To assist in emergency preparedness for a biological agent terrorist attack or accidental pathogen release, potential contaminant levels and migration pathways of spores spread by urban stormwater were evaluated using a Storm Water Management Model (SWMM) of U.S. Coast Guard Base Elizabeth City, North Carolina. The high temporal-spatial resolution SWMM model was built using spore concentrations in stormwater runoff from asphalt, grass, and concrete collected from a point-scale field study. The subsequent modeled contamination scenarios included a notional plume release and point releases mimicking the field study under three rainfall conditions. The rainfall scenarios included a 6-hour natural rainfall event on Dec. 8, 2021 and two design storms (2-year and 100-year events). The observed spore concentrations from asphalt and concrete from the actual field experiment were applied to calibrate the washoff parameters in the SWMM model, using an exponential washoff function. The calibrated washoff coefficient (c 1 ) and exponent (c 2 ) were 0.01 and 1.00 for asphalt, 0.05 and 1.45 for grass, and 2.45 and 1.00 for concrete, respectively. The calibrated SWMM model simulated spore concentrations in runoff at times and magnitudes similar to the field study data. In the point release modeled scenario, the concrete surface generated 55.6% higher average spore concentrations than asphalt. Similarly, in the field experiment, a 175% (p < 0.05) higher average spore concentration in surface runoff was observed from concrete than from asphalt. Here, this study demonstrates how SWMM may be used to evaluate spore washoff from urban surfaces under different precipitation amounts, intensities, and durations, and how visualized spatial migration pathways in stormwater runoff may be used for emergency planning and remediation.
Can artificial intelligence and data-driven machine learning models match or even replace process-driven hydrologic models for streamflow simulation?: A case study of four watersheds with different hydro-climatic regions across the CONUS
With recent developments in computational techniques, Data-driven Machine Learning Models (DMLs) have shown great potential in simulating streamflow and capturing the rainfall-runoff relationship in given watersheds, which are traditionally fulfilled by Process-based Hydrologic Models (PHMs). There are debates on whether the DMLs can outperform and possibly replace the classical PHMs for streamflow simulation and river forecasting, but no clear conclusions have been made. This study aims to investigate whether the newer DMLs have any potential in further improving the simulation accuracy of classical PHMs, and vice versa. To do this, we compared a few popular PHMs and DMLs over four watersheds across the Continental US (CONUS) that are associated with different input, climate, and regional conditions. A total of five hydrologic models were chosen, including (1) two classical lumped models, i.e., the Sacramento Soil Moisture Accounting (SAC-SMA) and Xinanjiang (XAJ); (2) one modern distributed model, termed Coupled Routing and Excess Storage (CREST); (3) and two DMLs including an Artificial Neural Networks (ANN) and a deep learning model, termed Long Short Term Memory (LSTM). Our results demonstrated that the DMLs still significantly biased when using the baseline input scenario with the PHMs. However, the DMLs fed with delayed input scenarios had great potential and can reach high simulation accuracy. The DMLs, especially the ANN, outperformed other employed models under the rainfall-runoff relationship in which rainfall dominantly drives. Furthermore, the DMLs also showed better performance in the high-flow regime, while the PHMs had a better performance for the low-flow regime, implying both PHMs and DMLs have their own merits and are worthy of joint development. In general, our study indicated a great potential of using DMLs to simulate streamflow, but further studies are still needed to verify the transferability and scalability of DMLs in large-scale experiments, such as the Distributed Model Intercomparison Projects 1&2 conducted by National Weather Services but to compare modern DMLs and PHMs.
First-Principles Modeling of Sodium Ion and Water Intercalation into Titanium Disulfide Interlayers for Water Desalination
Recent experiments revealed the possibility of using titanium disulfide (TiS 2 ) as the cathode material in capacitive deionization (CDI) devices for water desalination. Although it performed stably up to 70 cycles with a salt removal capacity of 14 mg/g (corresponding to a sodium ion removal capacity of 35.8 mg/g) at high molar concentration (600 mM NaCl), the maximum capacity of TiS 2 as a CDI electrode was much lower (~70 mAh/g) than as a supercapacitor (239 mAh/g). Understanding why the ion capacities of these two configurations of the same material differ will entail elucidating detailed charge/discharge mechanisms at the atomic scale. Here, we present density functional theory simulations of sodium intercalation into TiS 2 interlayers in order to gain such understanding. We systematically investigated TiS 2 stacking patterns and ion intercalation sites, energetics of the intercalated compounds, and phase transformation during sodium intercalation. The calculated structural evolution and capacitance–voltage curve agree quite well with previous measurements. We conclude that the different maximum capacities of TiS 2 measured in aqueous and dry environments originate from weaker interlayer interactions with respect to shear strain after ~33% intercalation of Na + –H 2 O pairs, which is detrimental to the mechanical stability of TiS 2 . In conclusion, this study sheds light on the underlying mechanisms of ion intercalation into layered materials and contributes to understanding requirements for future design and optimization of CDI electrode materials for water desalination.
Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite
Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.
The Contribution of Subtropical Moisture Within an Atmospheric River on Moisture Flux, Cloud Structure, and Precipitation Over the Salmon River Mountains of Idaho Using Moisture Tracers
The impact of an atmospheric river (AR) on the flux of subtropical moisture across Idaho's Salmon River Mountains and precipitation over the mountains is evaluated using the Weather, Research, and Forecasting model with water vapor tracers (WRF-WVT). The AR impacted Idaho between 17 and 19 January 2017 during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE) campaign. WRF-WVT is configured to isolate the subtropical moisture contribution to the AR, the moisture flux, and precipitation. Subtropical water vapor advected by the AR into Idaho is tagged and tracked in three-dimensional space throughout the run. This allows the contribution of the subtropical moisture to the vertical distribution of water vapor and the precipitation to be directly calculated. The simulated cloud structure is compared with airborne radar data collected during two SNOWIE intensive operation periods. This study found that more than 70% of the moisture flux and more than 80% of the precipitation across the Idaho Mountains during SNOWIE IOP 4 could be attributed to subtropical moisture within the AR. Nearly all of the moisture flux in the upper cloud and 50% of the moisture in the lower cloud was attributable to the subtropical moisture. The subtropical moisture contribution within the AR to precipitation ranged from 35% in northern Idaho to more than 90% in southern Idaho. Across the entire period of impact of the AR, more than 60% of precipitation in Idaho was attributable to the subtropical moisture within the AR, with this percentage increasing toward the south across the state.
Hydropower expansion in eco-sensitive river basins under global energy-economic change
Vast hydropower resources remain untapped globally, the deployment of which could provide energy-economic benefits but impact riverine ecosystems. Across eco-sensitive river basins, it is unclear how drivers of hydropower expansion, such as rapid economic growth and a low-carbon energy transition, could interact with countervailing forces, such as increasingly cost-competitive variable renewable energy (VRE). Using an integrated energy-water-economy model, we explore the effects of these forces on long-term hydropower expansion in the world’s 20 most eco-sensitive basins, which have high ecological richness and untapped hydropower potential. Here, we find that a low-carbon transition exerts the strongest development pressure, causing deployment exceeding 80% of exploitable potential in more than 72% of eco-sensitive basins by 2050, most of which have limited deployment today. Rapid economic growth induces such extensive deployment in only 44% of eco-sensitive basins. Enhanced integration of VRE reduces deployment, alleviating the impacts of rapid economic growth but not the low-carbon transition.
Ballistic Δ S = 2 intersystem crossing in a cobalt cubane following ligand-field excitation probed by extreme ultraviolet spectroscopy
We report femtosecond M 2,3 -edge X-ray absorption near-edge structure (XANES) spectroscopy is used to probe the excited-state dynamics of the cobalt cubane [Co III 4 O 4 ](OAc) 4 (py) 4 (OAc = acetate, py = pyridine), a model for water oxidation catalysts. After ligand-field excitation, intersystem crossing (ISC) to a metal-centered quintet occurs in 38 fs. 30% of the hot quintet undergoes ballistic back-ISC directly to the singlet ground state, with the remainder relaxing to a long-lived triplet.
p K a of alcohols dictates their reactivity with reduced uranium-substituted thiomolybdate clusters
The uranium-substituted thiomolybdate cluster, (Cp* 3 Mo 3 S 4 )UCp*, has been demonstrated as a model for water reduction by single uranium atoms supported on a molybdenum sulfide surface (U@MoS 2 ). In this study, the scope of O–H bond activation is expanded through the investigation of the reactivity of various alcohols with differing pKa values for the –OH proton. The reaction of (Cp* 3 Mo 3 S 4 )UCp* with stoichiometric amounts of methanol, phenol, 2,6-dichlorophenol, and nonafluoro-tert-butyl alcohol affords the corresponding mono-alkoxide species, (Cp* 3 Mo 3 S 4 )Cp*U(OR), via a uranium-metalloligand cooperative activation of the O–H bond. This observed reactivity is analogous to the O–H bond activation reported by (Cp* 3 Mo 3 S 4 )UCp* in the presence of water. However, addition of tert-butanol induces protonolysis of the Cp* ligand on uranium, resulting in the formation of a uranium tris-tert-butoxide cluster, (Cp* 3 Mo 3 S 4 )U(O t Bu) 3 . Independent synthesis of (Cp* 3 Mo 3 S 4 )Cp*U(O t Bu) was possible via an alternative pathway, eliminating sterics as a justification for the observed discrepancy in reactivity. Furthermore, these results offer insight into the role the –OH proton pKa plays in dictating the mechanism of O–H bond activation of alcohols by the uranium-substituted thiomolybdate cluster.
Study of self-interaction errors in density functional predictions of dipole polarizabilities and ionization energies of water clusters using Perdew-Zunger and locally scaled self-interaction corrected methods
Here, we studied the effects of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modelled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz. the local spin density approximation (LDA), the Perdew-Burke-Ernzherof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-Lowdin orbital SIC (FLO-SIC) framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.58 Bohr 3 with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but over-corrects, resulting in underestimated polarizabilities in SIC-LDA, -PBE, and -SCAN. Finally, we applied a recently proposed local-scaling SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute error of 0.08 Bohr 3 for LSIC-LDA and 0.06 Bohr 3 for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as an absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAE of 0.4 eV for LSIC-LDA and 0.04 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies while the LSIC-PBE ionization energies are more accurate than reported GW results.
Comparison of water nanodroplet properties on different graphite-based substrates
The molecular structure and dynamics of water differ considerably at various interfaces. We compare the interfacial water structure–property relationship on three different carbon substrates, namely, amorphous carbon, compressed expanded natural graphite, and pure graphite by utilizing atomistic molecular dynamics simulations. The effect of different substrates on the structural and dynamical properties of water can readily be observed. Here, the density distributions parallel and normal to the substrates show oblate droplet structures. The normal to the substrate water distribution shows a strong hydration layer at the interface that does not vary with substrates. However, the disparity in the structure and dynamics on three different substrates shows that the surface morphologies of the substrates are critical for determining nanoscale water properties. Furthermore, it is observed that the formation of an interfacial water layer or the hydration layer is a direct consequence of both water “confinement” at the nanoscale and “attraction” between water molecules and the carbon substrates.
Anomalous temperature dependence of the experimental x-ray structure factor of supercooled water
The structural changes of water upon deep supercooling were studied through wide-angle x-ray scattering at SwissFEL. The experimental setup had a momentum transfer range of 4.5 Å-1, which covered the principal doublet of the x-ray structure factor of water. The oxygen–oxygen structure factor was obtained for temperatures down to 228.5 ± 0.6 K. Similar to previous studies, the second diffraction peak increased strongly in amplitude as the structural change accelerated toward a local tetrahedral structure upon deep supercooling. We also observed an anomalous trend for the second peak position of the oxygen–oxygen structure factor (q 2 ). We found that q 2 exhibits an unprecedented positive partial derivative with respect to temperature for temperatures below 236 K. Based on Fourier inversion of our experimental data combined with reference data, we propose that the anomalous q 2 shift originates from that a repeat spacing in the tetrahedral network, associated with all peaks in the oxygen–oxygen pair-correlation function, gives rise to a less dense local ordering that resembles that of low-density amorphous ice. The findings are consistent with that liquid water consists of a pentamer-based hydrogen-bonded network with low density upon deep supercooling.
Isotope effects on the structural transformation and relaxation of deeply supercooled water
Herein, we have examined the structure of supercooled liquid D 2 O as a function of temperature between 185 and 255 K using pulsed laser heating to rapidly heat and cool the sample on a nanosecond timescale. The liquid structure can be represented as a linear combination of two structural motifs, with a transition between them described by a logistic function centered at 218 K with a width of 10 K. The relaxation to a metastable state, which occurred prior to crystallization, exhibited nonexponential kinetics with a rate that was dependent on the initial structural configuration. When the temperature is scaled by the temperature of maximum density, which is an isostructural point of the isotopologues, the structural transition and the non-equilibrium relaxation kinetics of D 2 O agree remarkably well with those for H 2 O.
Density functional theory of water with the machine-learned DM21 functional
The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn–Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. Furthermore, we find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.
Anomaly detection for MPC forecast in Fleet of Water Heaters
Among residential devices, water heaters consume 20% of home energy use in the United States. Water heaters possess the capability to store energy within their reservoirs, enabling the ability to decouple energy use from hot water use. This capability can be used to reduce energy usage and costs while also supporting grid services. This requires accurate forecasting of the parameters of the water heater such as upper and lower temperatures. In this study, we analyzed the performance and behavior of a water heater model used in the real-world to predict a control mechanism that is implemented in a smart residential neighborhood. The model forecasts are accurate in most cases but not all. In such scenarios, error correction of the model is necessary to further improve model predictive control accuracy. Anomaly detection is the first step of error correction. This study complements existing research by grouping time series data into two clusters one with anomalies and another without anomalies. To achieve this task, we explored and compared multiple unsupervised machine learning algorithms to perform clustering. Among these algorithms, Ward clustering has the lowest running time and identified the highest number of anomalies for the upper temperature limit. The proposed approach is tested based on the data collected in a neighborhood with 46 townhomes located in Atlanta, GA.
An Automatic Processing Framework for In Situ Determination of Ecohydrological Root Water Content by Ground-Penetrating Radar
This report concerns Root water content (RWC) as a vital component in water flux in soil-plant-atmosphere continuum. Knowledge of RWC helps to better understand the root function and the soil-root interaction and improves water cycle modeling. However, due to the lack of appropriate methods, field monitoring of RWC is seriously constrained. In this study, we used ground-penetrating radar (GPR), a common geophysical technique, to characterize RWC of coarse roots noninvasively. An automatic GPR data processing framework was proposed to (1) identify hyperbolic root reflections and locate roots in GPR images and (2) extract waveform parameters from the reflected wave of identified roots. These waveform parameters were then used to establish an empirical model and a semiempirical model to determine RWC. We validated the developed models using GPR root data at three antenna center frequencies (500 MHz, 900 MHz, and 2 GHz) that were produced from simulation experiments (with RWC ranging from 70% to 150%) and field experiments in sandy soils (with RWC ranging from 66% to 144%). Our results show that both the empirical and the semiempirical models achieved a good performance in estimating RWC with similar accuracy, i.e., the prediction error [root-mean-square error (RMSE)] was less than 8% for the simulation data and 12% for the field data. For both models, the accuracy of RWC estimation was the highest when applied to 2-GHz data. This study renders a new opportunity to determine RWC under field conditions that enhances the application of GPR for root study and the understanding and modeling of ecohydrology in the rhizosphere.