Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Predictive Control”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗
DOE OSTI2025

Temperature and Strain Sensing Characteristics of a 128° YX-Cut LiNbO 3 Rayleigh-Mode SAW Sensor From Room to Cryogenic Temperatures

Accurate, passive, and wireless monitoring of cryogenic hardware is essential for high-energy physics, space propulsion, and biomedical instrumentation. Here, this study quantifies the coupled temperature-strain behavior of Rayleigh-mode surface acoustic-wave (SAW) delay-line sensors fabricated on 128° YX-cut LiNbO 3 . A nonlinear finite element (FE) model incorporating Varshni-based elastic constants, higher-order thermal expansion, and temperature-dependent piezo- and dielectric coefficients was developed and validated experimentally between 280K and 80K. Free-standing (first test condition) and bonded/wired (second test condition) devices exhibited indistinguishable thermal responses; the average temperature coefficient of delay (TCD) in the critical cryogenic range from 130K down to 80K differed by only 0.15 ppm/K (0.32%), confirming that bonding-induced stress is negligible. Over 280-80K the measured TCD was 61.77 ppm/K, while the FE model predicted an equivalent temperature coefficient of frequency (TCF) of −62.74 ppm/K with an overall coefficient of determination R2 = 0.998. In the critical cryogenic interval 130-80K the TCD fell to 47.66 ppm/K, indicating improved thermal stability at low temperature. Controlled loading (0-300 με) revealed a strain coefficient of delay (SCD) that rises from 0.53 ± 0.02 ppm/με at 300K to 1.05 ± 0.02 ppm/με at 80K. This modest sensitivity confirms that, for temperature sensing, strain is a second-order perturbation above 135K but must be compensated at deeper cryogenic levels. Overall, this work establishes a predictive multiphysics model together with repeatable wired measurements that confirm the suitability of SAW sensors for temperature and strain monitoring in extreme cryogenic environments, while also providing a baseline for future wireless implementations.

25 ENERGY STORAGE↗
DOE OSTI2020

Managing time-substitutable electricity usage using dynamic controls

A predictive-control approach allows an electricity provider to monitor and proactively manage peak and off-peak residential intra-day electricity usage in an emerging smart energy grid using time-dependent dynamic pricing incentives. The daily load is modeled as time-shifted, but cost-differentiated and substitutable, copies of the continuously-consumed electricity resource, and a consumer-choice prediction model is constructed to forecast the corresponding intra-day shares of total daily load according to this model. This is embedded within an optimization framework for managing the daily electricity usage. A series of transformations are employed, including the reformulation-linearization technique (RLT) to obtain a Mixed-Integer Programming (MIP) model representation of the resulting nonlinear optimization problem. In addition, various regulatory and pricing constraints are incorporated in conjunction with the specified profit and capacity utilization objectives.

Ghosh, Soumyadip↗
DOE OSTI2026

Demonstration of gold nanorod systems for enhanced total efficiency: Experimental and numerical analysis

This study investigates the photothermal performance of gold nanorods engineered to exhibit longitudinal plasmon resonances at 695 nm, 780 nm, and 970 nm. The work combines synthesis, structural characterization, extinction measurements, numerical modeling, and controlled temperature experiments to quantify how nanorod geometry, resonance tuning, concentration, and chamber shape jointly influence heat generation. Transmission electron microscopy confirms that increasing nanorod aspect ratio systematically shifts the longitudinal plasmon peak toward the near-infrared region. Extinction measurements show strong agreement with theoretical predictions, validating the numerical model across two independent datasets. Three chamber geometries were tested under laser excitation at 640 nm, 808 nm, and 980 nm: an ascending stepped base, a flat base, and a descending stepped base. Without nanorods, the ascending geometry produced the highest efficiency due to enhanced natural convection. After introducing gold nanorods, all geometries exhibited substantial thermal enhancement, with total efficiencies exceeding 20%. The strongest improvement was obtained for nanorods resonant at 780 nm with a mass concentration of 4.6 mg/mL implemented on the descending stepped-base geometry. This performance resulted from the combined effect of spectral overlapping with the 808 nm laser, the highest nanorod concentration, and localized heat accumulation that intensified buoyancy-driven flow. The findings demonstrate that total efficiency is governed by a synergistic interplay between optical resonance, nanoparticle concentration, and macroscopic chamber design, revealing the system-level coupling between nanoscale plasmonic absorption and macroscale heat-transfer phenomena. The results provide a validated framework for tuning nanoscale plasmonic absorbers and optimizing thermal systems for applications requiring efficient light-to-heat conversion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Autoignition of cyclohexane at low-to-intermediate temperatures: Rapid compression machine experiments and improved comprehensive chemical kinetic model

Cycloalkanes are an important class of hydrocarbons found in conventional and alternative transportation fuels. Being one of the simplest practical cycloalkanes, cyclohexane is often used as a surrogate in studies to understand the oxidation, ignition, and sooting characteristics of cycloalkanes more broadly. Although cyclohexane has been extensively studied, very limited datasets for cyclohexane oxidation at low-to-intermediate temperatures are available in the literature. Moreover, the literature autoignition data from rapid compression machines (RCMs) did not provide information on volume histories or experimental heat transfer characteristics which can significantly impact model predictions, thereby hindering the detailed validation of chemical kinetic models. It has also been observed that the cyclohexane kinetic models developed so far are not able to predict the literature data across a wide range of operating and mixture conditions. Hence, further studies are recommended to understand and improve model predictions for cyclohexane oxidation. In this study, RCM experiments were conducted using cyclohexane/air mixtures at varying compressed pressures (P C =15, 20, and 30 bar), equivalence ratios (ϕ=0.5, 1.0, and 2.0), and compressed temperatures (T C =670–960 K). The current RCM results of cyclohexane were found to complement well with the literature data, with the first-stage IDTs at high pressures and volume histories at all test conditions being first reported for better development and validation of chemical kinetic models. Additionally, a detailed chemical kinetic model for the cyclohexane oxidation has been developed and validated using both newly acquired experimental IDTs in this study and literature experimental studies (IDTs in RCMs and STs, species profiles in jet-stirred reactors, and laminar burning velocities). Unlike literature models, the proposed model was observed to provide an overall better prediction for experimental results over a wide range of conditions. Chemical kinetic analysis of the developed model was further conducted to identify the reaction pathways controlling the autoignition of cyclohexane.

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

Relative Impact of Sea Ice and Temperature Changes on Arctic Marine Production

We use a modern Earth system model to approximate the relative importance of ice versus temperature on Arctic marine ecosystem dynamics. We show that while the model adequately simulates ice volume, water temperature, air-sea CO 2 flux, and annual primary production in the Arctic, itunderestimates upper water column nitrate across the region. This nitrate bias is likely responsible for the apparent underestimation of ice algae production. Despite this shortcoming, the model appears to be a useful tool for exploring the impacts of environmental change on phytoplankton production and carbon dynamics over the Arctic Ocean. Our experiments indicate that under a warmer climate scenario, the percentage of ocean warming that could be apportioned to a reduction in ice area ranged from 11% to 100%, while decreasing ice area could account for 22–100% of the increase in annual ocean primary production. The change to CO2 air-sea flux in response to ice and temperature changes averaged an Arctic-wide 5.5 Tg C yr –1 (3.5%) increase, into the ocean. This increased carbon sink may be short-lived, as ice cover continues to decrease and the ocean warms. The change in carbon fixation from phytoplankton in response to increased temperatures and reduced ice was generally more than a magnitude larger than the changes to CO 2 flux, highlighting the importance of fully considering changes to the marine ecosystem when assessing Arctic carbon cycle dynamics. Our work demonstrates the importance of ice dynamics in controlling ocean warming and production and thus the need for well-behaved ice and BGC models within Earth system models if we hope to accurately predict Arctic changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fatty acid foams for nonselective physical removal of microplastics from aqueous solutions

Microplastics (MPs) are pervasive environmental contaminants whose removal from water remains a major challenge due to their small size, chemical diversity, and dynamic surface properties arising from environmental aging/weathering. Here, we present a concept of foam-based separation method that physically traps MPs in the foam phase using microtubular assemblies of 12-hydroxystearic acid. These foams are stabilized by anisotropic fatty acid microtubules formed in the presence of ethanolamine, which jam within the foam channels and suppress fluid drainage thereby enhancing MP retention and foam stability. MPs of different sizes, polymer compositions (including polystyrene, polypropylene, polyethylene terephthalate, and polytetrafluoroethylene), and weathered states were retained in the foam phase without requiring chemical modification or relying on chemical interactions between the fatty acid and MPs. Thermally induced transition of the fatty acid microtubules into nanomicelles above the characteristic phase transition temperature (∼35 °C) enables controlled foam collapse and recovery of trapped MPs. The cumulative removal efficiency can exceed 85% through multiple foaming cycles, matching predictions from a probabilistic retention model. This work shows that foams can provide a simple platform to trap MPs, thus providing a new physical-removal strategy that does not rely on the particles’ chemistry.

Guillot, Kennedy A. [Louisiana State Univ., Baton ↗
DOE OSTI2024

A Deep Learning Approach for In-Network Synchrophasor Missing Data Recovery Using Programmable Network Switches

Phasor measurement unit (PMU) networks deliver accurate and timely measurements, which is essential for managing today’s electric power systems. To ensure data quality and enhance the cyber-resilience of PMU networks against malicious attacks and data errors, this study presents an online PMU missing data recovery scheme by leveraging P4 programmable switches. The data plane incorporates a customized PMU protocol parser that abstracts the necessary payload data for recovery. Recovery processes are executed in the control plane using a pre-trained machine learning model. Both traditional and advanced ML models, such as transformer and TimeGPT, are explicitly employed for data prediction. This approach ensures rapid and precise data recovery. Performance evaluations focus on recovery speed and accuracy, using a real dataset from a campus microgrid. With 20% missing PMU data, the mean absolute percentage error for voltage magnitude is 0.0384%, and the phase angle error discrepancy is approximately 0.4064%.

Phasor Measurement Unit, Machine Learning, Program↗
DOE OSTI2020

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗
DOE OSTI2026

The two radiative states of the Arctic atmosphere and their impacts on the surface energy budget of sea ice

The surface energy budget (SEB) is a central regulator of Arctic climate and sea ice evolution, yet its processes remain poorly constrained due to sparse observations and complex, coupled surface-atmosphere interactions. This study leverages year-long measurements from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) to provide the most comprehensive assessment to date of the central Arctic SEB and its modulation by atmospheric variability. Ship- and ice-based observations from October 2019 to September 2020 were used to directly measure or tightly constrain each term of the SEB, leading to exceptional energetic closure with the seasonal snow and ice mass balance. The analysis reveals strong seasonal transitions in atmosphere-surface energy transfer that are modulated by the atmospheric state and constrained by the ability of the surface temperature to respond. Classification of the atmosphere into its two dominant radiative states—the semi-transparent (ST) and opaque (OP)—highlights the central role of synoptic-scale variability in clouds. The ST atmospheric state dominated the long winter ice growth season, with limited cloudiness supporting persistent surface radiative cooling and ice growth. The OP state, associated with liquid-containing or thick ice clouds, became dominant in spring, with the combination of increased solar heating and cloud surface longwave warming driving ice and snow melt. Eddy covariance versus bulk approaches for deriving surface turbulent heat fluxes provide vastly different perspectives on the role of turbulence in modulating the SEB. These results establish a high-quality benchmark dataset for Arctic SEB studies and demonstrate how the balance of atmospheric radiative states exerts a first-order control on the annual evolution of the sea ice. The findings have broad implications for advancing observing technologies, understanding Arctic amplification, improving climate models, and predicting future sea ice change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Planar Systems for Quantum Information

This project aims to develop two‐dimensional (2D) moiré materials as a quantum simulator to implement model Hamiltonians and their phase diagrams. Progress in quantum information science (QIS) requires the development of advanced quantum materials systems. The rich family of layered van der Waals materials and their heterostructures present opportunities to create previously unrealized types of applications for QIS. Specifically, when two layers of van der Waals materials are overlaid with a small twist angle or/and lattice mismatch, a moiré superlattice with a period of about ten nanometers is formed. This provides a periodic trapping potential for electrons. Electrons can tunnel between the traps and repel each other by their mutual Coulomb interactions. The platform of 2D moiré materials provides many attractive features, including tunability of length and energy scales, charge density, and even lattice symmetry. It presents new possibilities for realizing quantum simulation of the many-body physics in a solid-state platform. This integrated team of six investigators seeks to develop relevant theoretical treatments to link ab-initio studies of 2D moiré materials to model Hamiltonians and to evaluate correlated phases predicted by these model Hamiltonians in the relevant regimes. On the experimental side, the team aims to develop methods to realize a homogeneous and highly controlled potential landscape for the electrons and to initiate, protect, and measure their quantum many-body states.

36 MATERIALS SCIENCE↗
DOE OSTI2022

MODELLING FIBER ORIENTATION DURING ADDITIVE MANUFACTURING COMPRESSION MOLDING PROCESSES

The production of high-performance thermoplastic composites reinforced with short carbon fibers can be achieved by a novel “additive manufacturing-compression molding” technique. An advantage of such combination is two-fold: controlled fiber orientation in additive manufacturing and less void content by compression molding. In this study, a computational fluid dynamics model has been developed to predict the behavior of printed layers during fiber-reinforced thermoplastic extrusion and subsequent compression molding. The fiber orientation was modelled with a simple quadratic closure model. The interaction between the fibers is included using a rotary diffusion coefficient which becomes significant in concentrated regimes. Finally, the second order orientation tensor is coupled with the momentum equation as an anisotropic part of the stress term. The effect of processing parameters on the behavior of printed layers was investigated to determine the favorable printing scenarios. The developed numerical model enables design of high-performance composites with tunable mechanical properties.

Seta, Berin↗
DOE OSTI2025

Unveiling the Hidden Evolution of Crystal Defects and Disorder in Energy Materials

Control of point defects and disorder in functional thin films and 2D materials is critical to realizing their full potential in applications ranging from energy storage to advanced electronics. However, these phenomena are often poorly understood, difficult to characterize, and challenging to direct with precision. This presentation explores emerging multi-modal computer vision to decipher and predict order in materials across multiple length scales in the electron microscope, from the atomic to the nanoscale. By fusing data from diverse sources, these powerful models provide unprecedented insights into materials' lifecycles, enabling the control of defects and their associated properties at a fundamental level. This capability promises to transform materials design and accelerate the development of next-generation technologies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Optimization of Energy Storage System Economics and Controls by Incorporating Battery Degradation Costs in REopt

The use of stationary electrochemical energy storage systems utilizing lithium-ion batteries has increased rapidly as the production scale and price for lithium-ion batteries has decreased. These energy storage systems are crucial for maintaining grid resiliency, especially for grids operating with high penetration of renewable energy generation assets or for with a variety of distributed energy generation and storage systems. One challenging factor for the development of battery energy storage systems is estimating the proper sizing, in terms of both power and energy, that minimizes total costs over the lifetime of the systems; this calculation is difficult in simple cases, where a battery is costed independently, but is extremely challenging when building loads and electrical generation by photovoltaic resources are also considered. REopt is a techoeconomic optimization tool developed by NREL to address these challenges. Previously, battery degradation has been priced by simply assuming a 10-year replacement schedule for battery systems. However, this does not account for varying degradation trends observed across real-world batteries, or allow for batteries to be operated in a degradation-aware manner that optimizes battery dispatch based on operating costs. This work incorporates a battery life model into REopt. This battery life model is simple, so that it may be solvable within the constrains of a mixed-integer linear optimization problem, but is fit to accelerated aging data recorded in the lab. To achieve the best possible accuracy for lifetime estimates given these constraints, parameters for the battery life model in REopt are estimated by fitting 20-year simulations of battery life after identifying state-space battery degradation model from accelerated aging data. Comparisons of battery life predicted in REopt and from the state-space battery degradation model to ensure validity of lifetime estimates made by REopt. Battery life and cost is optimized by controlling three decision to minimize system life cost: battery sizing, daily state-of-charge, and daily energy-throughput. The cost of battery degradation as a function of these control variables is then estimated assuming two possible maintenance strategies: replacement, where the entire battery system is replaced if cell reach an end-of-life capacity threshold; and augmentation, which establishes a fund to pay for continual purchase of new batteries to maintain the initial energy capacity of the system. These two strategies offer conservative (for replacement) and optimistic (for augmentation) bounds for total system cost. The degradation cost incurred by these strategies is then used to control battery dispatch decisions, operating the battery in a degradation-aware manner that maximizes battery lifetime while also providing energy when favorable. Because the mixed-integer linear program has perfect foresight of future energy needs, batteries with degradation costs are always operated using 'just-in-time' charging, which is unrealistic, as no energy is left in the storage system to perform other energy services or to serve as emergency back-up power. To combat this, an inequality constraint on the average annual state-of-charge is imposed, and the sensitivity of system cost to average stored energy, e.g., the cost of system resiliency, can be quantified. Analysis of results has several conclusions, for instance, oversizing of battery storage systems is not a cost burden when battery storage is an optimal solution, as any additional battery capacity can simply be utilized to avoid costs of purchasing energy from a utility.

battery↗
DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗