Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic State Estimation”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stress field measurements using quantitative schlieren

Quantitative schlieren analysis is extended here to optically transparent solids in quasi-static and dynamic experiments to measure stress distributions. The quasi-static experiments in polymethyl methacrylate (PMMA) compared refraction angles and stress gradients calculated from schlieren images to the analytical Flamant solution of a line load on a half-space. The quantitative schlieren measurements of the stress field in the thin sample with a load compared well to the analytical solution. The analysis method was then extended to explosive induced shock waves in PMMA. The explosive induced response of PMMA was experimentally studied using high-speed schlieren to visualize the shock propagation in conjunction with Photon Doppler Velocimetry (PDV) to record surface velocity histories. The stress state estimated from the schlieren images was compared to the stress calculated from the PDV measurements. High-speed imaging limitations caused the shock wave to not be fully resolved in the images, but was resolved in the PDV measurement. The stress state behind the shock calculated from the high-speed images followed a similar trend to the stress calculated from the PDV measurements.

Physics↗
DOE OSTI2025

Assessing High Burnup U-19Pu-10Zr Fuel Performance against Historical and Modeled Behavior

Advancing the deployment of sodium-cooled fast reactors (SFRs) requires thorough testing of metallic fuel pins under accident conditions to establish safe operational limits of high burnup fuel. To conduct transient testing, a comprehensive understanding of steady-state fuel behavior obtained through both experimental characterization and accurate predictive capabilities is needed. This study comparatively assesses the steady-state irradiation performance of two high burnup U-19Pu-10Zr fuel pins, DP-36 and DP-40, irradiated under prototypic fast reactor conditions in preparation for planned safety testing at the Transient Reactor Test Facility. Since DP-40 was designated for use in the test and DP-36 serves as its sibling pin, non-destructive, engineering-scale post-irradiation examinations (PIE) were conducted on both pins while destructive examinations were performed exclusively on DP-36. The results were then assessed against historical performance data from similar fuel pins irradiated in the Experimental Breeder Reactor-II. Additionally, the steady-state irradiation of each pin was modeled using the BISON fuel performance code to assess the accuracy of current modeling capabilities in predicting the baseline irradiation behavior. Non-destructive examinations included neutron radiography to measure fuel column elongation, gamma scanning to verify pin integrity and fission product migration, and profilometry to assess dimensional changes. Benchmarking against existing PIE data revealed consistent patterns in axial fuel column growth and cladding diametral strain, though both pins exhibited longer low-density “fluff” structures, which can have implications for core reactivity and source term calculations. Destructive examinations on DP-36 included fission gas release analysis and sectioning for optical microscopy, which showed more complex constituent redistribution patterns than the traditionally accepted 3-ring model. The axial evolution of fractional areas and porosities of each of the redistributed zones were quantified and presented. Modeling comparisons showed agreement in fractional fission gas release but consistently overestimated axial and radial swelling and disagreed with measured axial porosity patterns. These conservative overpredictions suggested that the pins would appear closer to failure or operational limits at the start of transient tests, potentially leading to higher strain accumulation during the transient. While conservative estimates provide safety margins, they can negatively impact fuel economics. A review of the swelling models identified areas for improvement in the gaseous swelling, solid swelling, and fuel hot-pressing models when applied to ternary fuel. The results of this study highlight the critical importance of conducting pre-test characterization on both test and sibling pins to accurately capture steady-state fuel behavior, providing a precise baseline for post-test evaluations and essential inputs for transient modeling of the planned experiments. The analysis also revealed significant data gaps that require further investigation to enhance the understanding and prediction of fuel swelling and pore dynamics. Collecting comprehensive data across different irradiation conditions, burnup levels, and fuel compositions are essential for refining existing models and developing mechanistic models for both binary and ternary metallic fuels, ultimately improving the integration of modeling and experimental approaches in accident testing.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Dynamic Properties of FeCrMnNi, a High Entropy Alloy

High entropy alloys (HEAs) represent a new class of materials that potentially possess an ideal combination of high strength and ductility and may be candidate materials that can withstand high amounts of damage before failure, especially under extreme conditions. To this end, a series of shock/spall loading experiments were conducted on an Fe/Cr/Mn/Ni HEA to assess not only its equation of state (EOS) but also its resistance to damage and failure. A single-stage light gas gun was used to conduct a series of flyer plate symmetric impact experiments coupled with recovery of the samples for post-mortem microstructural analysis to obtain fundamental understanding regarding the damage morphology. Our results show that the measured EOS for this alloy was in reasonable agreement with an estimated Hugoniot value. Additionally, the spall strength of these alloys was found to vary between 1.4 and 1.9 GPa within the same plate but was comparable to other iron-based conventional alloys. DOE/NV/03624--1059

36 MATERIALS SCIENCE↗
DOE OSTI2025

Predictions of Cell-to-Cell Propagation and Vent Gas Production in the Thermal Runaway of Lithium-Ion Battery Stacks

This work presents the thermal runaway propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway) as an efficient tool to predict different cell-to-cell thermal runaway propagation scenarios. Here, we explored the vent gas volume production and reaction duration highlighting the relationship between these parameters and thermal runaway propagation due to convection by the vented gases. Two metrics based on gas production rate and heating rate are utilized as good indicators of the start and end of thermal runaway. LIM1TR results are compared with and validated by experiments from the literature for single-cell and multicell array experiments of 5 Ah and 10 Ah cells. By accounting for intraparticle diffusion of reacting species in the electrodes, we were able to capture the general dynamics of thermal runaway propagation and estimate acceptable reaction durations compared with the experimental values. Simulation results further demonstrated that varying heating modes lead to distinct reaction durations, consistent with experimental observations. Vent gas volume predictions indicate the need to consider both full and partial oxidation of the electrolyte. The outcomes of this work are building blocks for further investigations of module-to-module propagation by vented gases through convective heat transfer.

25 ENERGY STORAGE↗
DOE OSTI2025

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗
DOE OSTI2022

Deblurring for nuclei: 3D characteristics of heavy-ion collisions

Observables from nuclear and high-energy experiments can be degraded by detector performance and/or methodology in extracting the observables, such as of the final-state characteristics of heavy-ion collisions in relation to a coarsely estimated reaction-plane direction. We propose the use of deblurring methods, such as in optics, to correct for observable degradation. Our main focus is the restoration of triple-differential particle distributions in heavy-ion collisions. We demonstrate that these could be extracted from collision measurements following the Richardson-Lucy deblurring method from optics. We illustrate basic features of the restoration methodology in a schematic model assuming either ideal or more realistic particle detection. Here, the inferred three-dimensional (3D) distributions for collisions may be easier to interpret in terms of collision dynamics and sought properties of bulk matter than the currently employed Fourier coefficients, that combine information from different azimuthal angles relative to the reaction plane.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Better practices for inferring ecosystem water use strategy from eddy covariance data

Eddy covariance data are critical for inferring ecosystem water use strategies. Yet, such inferences are sensitive to a range of assumptions applied across studies, hindering our understanding of water use strategies within and across eddy covariance sites. A recent analysis across 151 FLUXNET2015 and AmeriFlux-FLUXNET datasets found that poor model performance was the key driver of non-robust inferences of ecosystem water use strategies. Here, we leverage this previous analysis to (i) identify the specific assumptions that improve inference model performance across most sites, (ii) explain the mechanisms behind the performance improvements, and (iii) check whether better performance improves water use inference. We find that the common practice of fitting a model to canopy conductance (G c ) derived from the evapotranspiration (ET) observations, rather than to observed ET itself, artificially amplifies data errors and degrades the model performance. Next, accounting for vegetation dynamics by applying a growing season filter or incorporating satellite LAI data improves performance, but the former practice may remove soil water stress periods. Lastly, using the leaf-to-air vapor pressure deficit (VPD l ) derived from ET observations as a model input may artificially inflate performance. Based on these results, we recommend selecting observed ET (rather than derived G c ) as the response variable, carefully accounting for vegetation dynamics, and avoiding derived VPD l as a model input; these best practices improve model performance by c. 20% and robustness by c. 80% across all eddy covariance sites. Nevertheless, the performance improvements do not always correspond to more robust inference of water use strategies, as model parameter selection and surface energy budget closure corrections still strongly influence the ecosystem water use parameter estimation in a site-specific manner.

AmeriFlux↗
DOE OSTI2023

Efficient Super-Resolution of Near-Surface Climate Modeling Using the Fourier Neural Operator

Downscaling methods are critical in efficiently generating high-resolution atmospheric data. However, state-of-the-art statistical or dynamical downscaling techniques either suffer from the high computational cost of running a physical model or require high-resolution data to develop a downscaling tool. Here, we demonstrate a recently proposed zero-shot super-resolution method, the Fourier neural operator (FNO), to efficiently perform downscaling without the need for high-resolution data. Because the FNO learns dynamics in Fourier space, FNO is a resolution-invariant emulator; it can be trained at a coarse resolution and produces emulation at any high resolution. We applied FNO to downscale a 4-km resolution Weather Research and Forecasting (WRF) Model simulation of near-surface heat-related variables over the Great Lakes region. The FNO is driven by the atmospheric forcings, and topographic features used in the WRF model at the same resolution. We incorporated a physics-constrained loss in FNO by using the Clausius-Clapeyron relation to better constrain the relations among the emulated states. Trained on merely 600 WRF snapshots at 4-km resolution, the FNO shows comparable performance with a widely used convolutional network, U-Net, achieving averaged modified Kling-Gupta Efficiency of 0.88 and 0.94 on the test dataset for temperature and pressure, respectively. We then employed the FNO to produce 1-km emulations to reproduce the fine climate features. Further, by taking the WRF simulation as ground truth, we show consistent performances at the two resolutions, suggesting the reliability of FNO in producing high-resolution dynamics. Our study demonstrates the potential of using FNO for zero-shot super-resolution in generating first-order estimation on atmospheric modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Challenges of COVID-19 Case Forecasting in the US, 2020–2021

During the COVID-19 pandemic, forecasting COVID-19 trends to support planning and response was a priority for scientists and decision makers alike. In the United States, COVID-19 forecasting was coordinated by a large group of universities, companies, and government entities led by the Centers for Disease Control and Prevention and the US COVID-19 Forecast Hub ( https://covid19forecasthub.org ). We evaluated approximately 9.7 million forecasts of weekly state-level COVID-19 cases for predictions 1–4 weeks into the future submitted by 24 teams from August 2020 to December 2021. We assessed coverage of central prediction intervals and weighted interval scores (WIS), adjusting for missing forecasts relative to a baseline forecast, and used a Gaussian generalized estimating equation (GEE) model to evaluate differences in skill across epidemic phases that were defined by the effective reproduction number. Overall, we found high variation in skill across individual models, with ensemble-based forecasts outperforming other approaches. Forecast skill relative to the baseline was generally higher for larger jurisdictions (e.g., states compared to counties). Over time, forecasts generally performed worst in periods of rapid changes in reported cases (either in increasing or decreasing epidemic phases) with 95% prediction interval coverage dropping below 50% during the growth phases of the winter 2020, Delta, and Omicron waves. Ideally, case forecasts could serve as a leading indicator of changes in transmission dynamics. However, while most COVID-19 case forecasts outperformed a naïve baseline model, even the most accurate case forecasts were unreliable in key phases. Further research could improve forecasts of leading indicators, like COVID-19 cases, by leveraging additional real-time data, addressing performance across phases, improving the characterization of forecast confidence, and ensuring that forecasts were coherent across spatial scales. In the meantime, it is critical for forecast users to appreciate current limitations and use a broad set of indicators to inform pandemic-related decision making.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

A scalable variational method for estimating the latent infection-rate field of an outbreak

In this paper, we explore whether the infection-rate of a disease can serve as a robust monitoring variable in epidemiological surveillance algorithms. The infection-rate is dependent on population mixing patterns that do not vary erratically day-to-day; in contrast, daily case-counts used in contemporary surveillance algorithms are corrupted by reporting errors. The technical challenge lies in estimating the latent infection-rate from case-counts. Here we devise a Bayesian method to estimate the infection-rate across multiple adjoining areal units, and then use it, via an anomaly detector, to discern a change in epidemiological dynamics. We extend an existing model for estimating the infection-rate in an areal unit by incorporating a Markov random field model, so that we may estimate infection-rates across multiple areal units, while preserving spatial correlations observed in the epidemiological dynamics. To carry out the high-dimensional Bayesian inverse problem, we develop an implementation of mean-field variational inference specific to the infection model and integrate it with the random field model to incorporate correlations across counties. The method is tested on estimating the COVID-19 infection-rates across all 33 counties in New Mexico using data from the summer of 2020, and then employing them to detect the arrival of the Fall 2020 COVID-19 wave. We perform the detection using a temporal algorithm that is applied county-by-county. We also show how the infection-rate field can be used to cluster counties with similar epidemiological dynamics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Reducing heat load density with asymmetric and inclined double-crystal monochromators: principles and requirements revisited

Asymmetric double-crystal monochromators (aDCMs) and inclined DCMs (iDCMs) can significantly expand the X-ray beam footprint and consequently reduce the heat load density and gradient. Based on rigorous dynamical theory calculations, the major principles and properties of aDCMs and iDCMs are presented to guide their design and development, particularly for fourth-generation synchrotrons. In addition to the large beam footprint, aDCMs have very large bandwidths (up to ∼10 eV) and angular acceptance, but the narrow angular acceptance of the second crystal requires precise control of the relative orientations and strains. Based on Fourier coupled-wave diffraction theory calculations, it is rigorously proved that the iDCM has almost the same properties as the conventional symmetric DCM, including the efficiency, angular acceptance, bandwidth, tuning energy range and sensitivity to misalignment. The exception is that, for the extremely inclined geometry that can achieve very large footprint expansion, the iDCM has (beneficially) a larger bandwidth and wider angular acceptance. Inclined diffraction has the `rho-kick effect' that can be cancelled by the second reflection of the iDCM (even with misalignment), except that inhomogeneous strains may cause non-uniform rho-kick angles. At present, fabrication/mounting-induced strains pose low risk since they can be controlled to <0.5 µrad over large areas. The only uncertain challenge is the thermally induced strains, yet it is estimated that these strains are naturally lowered by the large footprint and may be further mitigated by optimized cryogenic cooling to the 1–2 µrad level. Overall, aDCMs and iDCMs have more stringent requirements than normal DCMs, but they are feasible schemes in practice.

asymmetric monochromator↗
DOE OSTI2023

Structural characterization of tin in toothpaste by dynamic nuclear polarization enhanced 119 Sn solid-state NMR spectroscopy

Stannous fluoride (SnF 2 ) is an effective fluoride source and antimicrobial agent that is widely used in commercial toothpaste formulations. The antimicrobial activity of SnF 2 is partly attributed to the presence of Sn(II) ions. However, it is challenging to directly determine the Sn speciation and oxidation state within commercially available toothpaste products due to the low weight loading of SnF 2 (0.454 wt% SnF 2 , 0.34 wt% Sn) and the amorphous, semi-solid nature of the toothpaste. Here, we show that dynamic nuclear polarization (DNP) enables 119 Sn solid-state NMR experiments that can probe the Sn speciation within commercially available toothpaste. Solid-state NMR experiments on SnF 2 and SnF 4 show that 19 F isotropic chemical shift and 119 Sn chemical shift anisotropy (CSA) are highly sensitive to the Sn oxidation state. DNP-enhanced 119 Sn magic-angle turning (MAT) 2D NMR spectra of toothpastes resolve Sn(II) and Sn(IV) by their 119 Sn chemical shift tensor parameters. Fits of DNP-enhanced 1D 1 H → 119 Sn solid-state NMR spectra allow the populations of Sn(II) and Sn(IV) within the toothpastes to be estimated. This analysis reveals that three of the four commercially available toothpastes contained at least 80% Sn(II), whereas one of the toothpaste contained a significantly higher amount of Sn(IV).

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

A model-independent data assimilation (MIDA) module and its applications in ecology

Abstract. Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

58 GEOSCIENCES↗
DOE OSTI2026

Differences in cluster and internal wake effects from mesoscale and large-eddy simulations off the US East Coast

Mesoscale simulations are increasingly used to estimate wake effects within and between large wind farms, despite limited validation for large-scale wake effects. This study evaluates the capabilities and limitations of mesoscale simulations in capturing wake-induced impacts on wind turbine power production through a direct comparison with large-domain large-eddy simulations (LESs) for three planned offshore wind farms under realistic atmospheric conditions and a range of atmospheric stabilities. We assess mesoscale performance in replicating wake characteristics behind single and multiple turbine clusters and quantify the resulting variability in mean turbine power. Results show that mesoscale Weather Research and Forecasting simulations with the Fitch wind farm parameterization capture key features of the velocity deficit downstream of both single and multiple wind farms, with mean root-mean-square errors near 5 % and good agreement with stability-driven wake behavior. However, in these simulations, the mesoscale Fitch parameterization underestimates power losses from internal wake effects, particularly when turbines align with the prevailing wind direction or under stable stratification. In these conditions, individual wakes persist and dominate downstream power deficits. The coarse resolution of the mesoscale simulations limits their ability to resolve individual wind turbine wakes that drive power fluctuations within wind farms. Nonetheless, mesoscale simulations can yield accurate estimates of combined wake losses from internal and cluster effects across some wind direction sectors, where errors in wake representation may cancel each other out. These findings underscore the strengths of mesoscale simulations for capturing broader wake patterns while highlighting their limitations for modeling turbine-level power losses. Future work should explore hybrid modeling approaches to capture both long-range cluster wake propagation and localized internal wake dynamics.

17 WIND ENERGY↗
DOE OSTI2023

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗
DOE OSTI2021

Dynamic thermal behavior of polycrystalline LaB 6 hollow cathodes

Lanthanum hexaboride (LaB6) hollow cathodes have demonstrated a capability for long life operation, which is critical to many space exploration missions. Thermal characterization of LaB6 hollow cathodes has revealed lower than expected electron emitter temperatures when the cathode reaches a steady state. This phenomenon is observed at discharge currents ranging from 5 to 35 A and xenon mass flow rates of 5–25 SCCM in cathodes with three different orifice diameters. Thus, the currently accepted value of the work function for polycrystalline LaB 6 , 2.67 eV, does not describe well the emission characteristics of LaB6 hollow cathodes operating with internal gas discharges at a steady state. We use empirically measured temperatures combined with a model of the hollow cathode emitter and xenon discharge to estimate the value of the work function, yielding a value ranging from 2.1 to 2.44 eV. This lower work function value implies that LaB6 hollow cathodes are expected to have even longer lifetimes than previously anticipated, further establishing them as a more suited alternative to other conventional cathode technologies for the task of long duration travel. Direct measurements of the work function as a function of depth on a hollow cathode emitter using x-ray photoelectron spectroscopy and ion beam milling indicate that the work function decreases with depth. We postulate several mechanisms that could explain the observed work function enhancement. Altogether, our results have important implications to the design, study approach, and operation of LaB 6 cathodes and potentially other cathodes with hollow configuration. Finally, our work opens the question of why the work function is reduced upon interaction with Xe plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Impacts of High Variable Renewable Energy Futures on Electric-Sector Decision Making: Demand-Side Effects

Previous work by the Berkeley Lab describes how high shares of variable renewable energy (VRE) such as wind and solar power could change wholesale electricity price dynamics. These include the timing of when electricity is cheap or expensive, locational differences in the cost of electricity, and the degree of regularity or predictability in those costs. Many decentralized decision-makers on the demand-side may not yet have considered the implications of these possible future changes. In this report, we evaluate the sensitivity of a set of demand-side decisions to different levels of VRE penetration ranging from a low of 5-20% to a high of 40-50%. The analysis builds on hourly wholesale energy and capacity prices in different VRE scenarios for four wholesale markets in the United States for the year 2030 (CAISO, ERCOT, NYISO, and SPP). The principal question for this exploration is whether private and public electric-sector decisions that are made based on assumptions reflecting low VRE levels still achieve their intended objective in a high VRE scenario with 40-50% wind and solar? This scoping report evaluates the impacts of changing patterns of peak system needs on the benefits of demand reductions by examining the altered value of different energy efficiency (EE) measures. Similarly, we investigate new opportunities for large energy consumers that may arise from periods with very low wholesale electricity prices. We calculate the value of new process investments (e.g., hydrogen production and other generalized electro-commodities), showcase the varying value of new product storage investments (such as reservoir extensions at a desalination plant), and estimate the benefits of increased process flexibility that uses electricity as a process-input in addition to traditional fossil fuels (e.g., district energy systems). Finally, many decentralized decision-makers and end-use customers are not directly exposed to wholesale electricity prices but instead receive price signals from their retail electricity rates. As wind and solar shares increase, we compare the economic efficiency of flat retail rates relative to more dynamic time-of-use tariffs with and without critical peak-pricing events.

14 SOLAR ENERGY↗
DOE OSTI2020

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗