Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “uncertain systems”

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 271 records · Page 15

Integrating Maximum Entropy Production Theory and Machine Learning to Improve Global Evapotranspiration Modeling

Accurate estimation of terrestrial evapotranspiration (ET) is vital for understanding global water and energy cycles. However, current global ET estimations are not well constrained. This study introduces an integrated framework combining the Maximum Entropy Production (MEP) theory with Random Forest (RF) model to improve global ET estimation. Specifically, in contrast to direct ET estimation by the RF model, the integrated framework (MEP‐RF) trains to predict error of MEP‐simulated ET. MEP‐RF outperforms RF in spatiotemporal extrapolation. Attribution analysis with in situ observations reveals that the inputs of MEP are the most critical variables for the ET process, including net radiation, vegetated area, soil moisture, and surface temperature. We further drive MEP‐RF with global reanalysis and satellite data sets of these four inputs, yielding a global mean terrestrial ET of 548 mm/year, with 77% attributed to transpiration. The global ET increased at a rate of 0.85 mm/year per year during 2003–2021, primarily due to vegetation greening rather than rising temperature, while decreasing soil moisture led to decreasing regional ET. The integrated framework provides a novel approach for the estimation of global ET without the need for hard‐to‐obtain and thus uncertain inputs, such as wind speed, surface roughness, aerodynamic and canopy stomatal resistance. Therefore, MEP‐RF offers an independent method on existing global ET products. It represents a promising physically based approach that can be incorporated into Earth System Models to enhance water and energy cycle simulations.

54 ENVIRONMENTAL SCIENCES↗

Methods for Assessing Opportunities for Ring Dam Pumped Storage Hydropower

There is growing interest in new pumped storage hydropower (PSH) deployment to provide a range of grid flexibility, reliability, and resiliency services under an evolving and uncertain future power sector. The National Laboratory of the Rockies develops open PSH resource assessment and cost modeling tools to help evaluate PSH deployment opportunities, and this report describes expansions to those tools to consider an additional PSH system configuration - ring-dam reservoirs built on flat topographical features that are constructed from roller-compacted concrete material. This reservoir type is common among current PSH proposals and requires new methods to identify sites with this reservoir geometry throughout the United States and characterize the associated dam cost. Cost characterization for ring dam reservoirs required collecting historical dam cost data for earthen, rockfill, and roller-compacted concrete dams and regressing equations that relate costs between alternative materials. The ring dam site identification algorithm follows a 5-step procedure to identify circular geometry reservoirs. Once ring dam reservoirs are identified, they are then paired with potential dry-gully reservoirs, and the full set of potential paired reservoirs is cost-optimized to produce a least-cost set of potential PSH sites with no overlapping reservoirs. The resulting analysis found 1,663 ring-dam to dry-gully systems in the contiguous United States that are lower cost than any overlapping dry-gully to dry-gully systems, 29 in Alaska, and none in Hawaii or Puerto Rico. These systems constitute 1.5 TW of capacity in the contiguous United States and nearly 29 GW in Alaska, demonstrating that under suitable topography and head, ring-dam systems can provide cost-effective PSH opportunities. The greatest density of these opportunities are found in the intermountain west where there are mesas and flat land at bases of mountain ranges, but continued work could incorporate additional site characteristics or consider more complex reservoir shapes to find additional PSH deployment opportunities.

13 HYDRO ENERGY↗

Marine energy supported multi-energy system planning and operation optimization for sustainable coastal community

The growing need for sustainable energy solutions in coastal areas necessitates the development of integrated systems that leverage abundant marine resources. In this study, a standalone Marine Energy Supported Multi-Energy System (MRE-MES) is designed for sustainable coastal community development, utilizing renewable marine resources, including offshore wind, wave, and solar energy, to address the energy needs of electricity, heat, freshwater, and hydrogen. The proposed MRE-MES incorporates a co-optimization model that simultaneously balances capacity planning and operational efficiency to minimize costs and environmental impacts. The system is tested under different renewable energy penetration levels and demand uncertainties, using a two-stage stochastic programming to account for variability in renewable resources and consumption needs. The experimental results indicate that in the optimal system capacity configuration, the percentage of total renewable energy generation is around 80 %, with or without capacity limitation constraints on PV, water tank, and hydrogen storage. Compared to the worst-case scenario in Monte Carlo experiments, two-stage stochastic optimization results in a more robust decision that effectively mitigates the risks posed by future uncertain demand conditions. In conclusion, the findings highlight the viability of marine energy for providing a resilient, comprehensive energy solution to coastal communities.

Capacity planning↗

In situ Study on the Evolution of Atomic and Electronic Structure of LaTiO 3 /SrTiO 3 System

As with other transition-metal oxide interfaces, LaTiO 3 /SrTiO 3 interfaces exhibit interesting electronic properties, such as a two-dimensional electron gas. While the available Mott states and multiple pathways to metallicity in LaTiO 3 can lead to a variety of applications in oxide electronics, the origin of the different properties remains uncertain. Herein, utilizing real-time, in situ x-ray diffraction at the synchrotron, robust but thin LaTiO 3 /SrTiO 3 heterostructures were fabricated using oxide molecular beam epitaxy. A thickness of six unit cells was determined to be critical, at which point the LaTiO 3 /SrTiO 3 heterostructures begin to sustain sharp, well-defined interfaces. The layers of LaTiO 3 and SrTiO 3 were sufficiently thin to benefit from the finite escape length of electrons at resonant soft x-ray photoemission, and the heterostructures were transferred in vacuo after growth to another chamber for synchrotron-based x-ray photoemission spectroscopy studies. With the x-ray energy tuned to the Ti-2p resonance, the LaTiO 3 /SrTiO 3 interfaces could be probed with enough penetration depth and selectivity. It is shown that all of the heterointerfaces exhibit the two-dimensional electron gas.

36 MATERIALS SCIENCE↗

Multistage economic MPC for systems with a cyclic steady state: A gas network case study

Multistage model predictive control (MPC) provides a robust control strategy for dynamic systems with uncertainties and a setpoint tracking objective. Moreover, extending MPC to minimize an economic cost instead of tracking a pre-calculated optimal setpoint improves controller performance. This paper presents a novel multistage economic nonlinear model predictive control (E-NMPC) framework for dynamic systems operating under uncertainty, with specific application to natural gas transmission networks. A key innovation lies in the integration of cyclic steady-state (CSS) constraints within the multistage MPC formulation, enabling the controller to manage periodic operating conditions commonly observed in energy systems. A Lyapunov-based descent condition is enforced to ensure robust stability of the controller. The multistage economic MPC framework is validated on two gas pipeline case studies, where it successfully minimizes net energy consumption, respects operational constraints under uncertain demand profiles, and guides the network to optimal cyclic operation. The Lyapunov function remains bounded in both case studies, validating the robust stability of multistage E-NMPC.

03 NATURAL GAS↗

EMIS-CompetitiveEquilibrium.jl [SWR-19-56]

The Electricity Markets Investment Suite - Competitive Equilibrium (EMIS-CE) model is an optimization-based power system investment model developed at NREL. It attempts to approximate a market-driven investment equilibrium outcome under the assumption of idealized or perfect competition (e.g. no participants are able to exert market power or leverage exclusive cost advantages) and uncertain future conditions.

Stephen, Gord [National Renewable Energy Lab. (NRE↗

Optimization under uncertainty of a hybrid waste tire and natural gas feedstock flexible polygeneration system using a decomposition algorithm

Market uncertainties motivate the development of flexible polygeneration systems that are able to adjust operating conditions to favor production of the most profitable product portfolio. However, this operational flexibility comes at the cost of higher capital expenditure. A scenario-based two-stage stochastic nonconvex Mixed-Integer Nonlinear Programming (MINLP) approach lends itself naturally to optimizing these trade-offs. This work studies the optimal design and operation under uncertainty of a hybrid feedstock flexible polygeneration system producing electricity, methanol, dimethyl ether, olefins or liquefied (synthetic) natural gas. A recently developed C++ based software framework (named GOSSIP) is used for modeling the optimization problem as well as its efficient solution using the Nonconvex Generalized Benders Decomposition (NGBD) algorithm. Two different cases are studied: The first uses estimates of the means and variances of the uncertain parameters from historical data, whereas the second assesses the impact of increased uncertain parameter volatility. The value of implementing flexible designs characterized by the value of the stochastic solution (VSS) is in the range of 260–405 M$ for a scale of approximately 893 MW of thermal input. Increased price volatility around the same mean results in higher expected net present value and VSS as operational flexibility allows for asymmetric exploitation of price peaks.

42 ENGINEERING↗

Optimal Eco-Driving Control of Connected and Autonomous Vehicles Through Signalized Intersections

This article focuses on the speed planning problem for connected and automated vehicles (CAVs) communicating to traffic lights. The uncertainty of traffic signal timing for signalized intersections on the road is considered. The eco-driving problem is formulated as a data-driven chance-constrained robust optimization problem. Effective red-light duration (ERD) is defined as a random variable, and describes the feasible passing time through the signalized intersections. Usually, the true probability distribution for ERD is unknown. Consequently, a data-driven approach is adopted to formulate chance constraints based on empirical sample data. This incorporates robustness into the eco-driving control problem with respect to uncertain signal timing. Dynamic programming (DP) is employed to solve the optimization problem. The simulation results demonstrate that the proposed method can generate optimal speed reference trajectories with 40% less vehicle fuel consumption, while maintaining the arrival time at a similar level compared to a modified intelligent driver model (IDM). The proposed control approach significantly improves the controller's robustness in the face of uncertain signal timing, without requiring to know the distribution of the random variable a priori.

33 ADVANCED PROPULSION SYSTEMS↗

Risk-Averse Bi-Level Stochastic Network Interdiction Model for Cyber-Security Risk Management

Security of cyber networks is crucial; recent severe cyber-attacks have had a devastating effect on many large organizations. The attack graph, which maps the potential attack paths of a cyber network, is a popular tool for analyzing cyber system vulnerability. In this study, we propose a bi-level stochastic network interdiction model on an attack graph to enable a risk-averse, resource constrained cyber network defender to optimally deploy security countermeasures that protect against attackers with an uncertain budget. This risk- averse conditional-value-at-risk (CVaR) model minimizes a weighted sum of the expected maximum loss over all scenarios and the expected maximum loss from the most damaging attack scenarios. We develop a customized constraint and column generation algorithm to solve our model as well as several acceleration techniques to improve the computational efficiency. Numerical experiments demonstrate that the acceleration techniques enable the solution of relatively large problems within a reasonable amount of time: applying all the acceleration techniques also reduces the average computation time of the basic algorithm by 71% for 100-node graphs. Using metrics called mean-risk value of stochastic solution and value of risk-aversion, computational results suggest that our stochastic risk-averse model significantly outperforms deterministic and risk-neutral models when 1) the distribution of attacker budget is heavy-right-tailed and 2) the defender is highly risk-averse.

Bhuiyan, Tanveer H.↗

Shallow-Water Habitat in the Lower Columbia River Estuary: A Highly Altered System

Abstract Decreases in shallow-water habitat area (SWHA) in the Lower Columbia River and Estuary (LCRE) have adversely affected salmonid populations. We investigate the causes by hindcasting SWHA from 1928 to 2004, system-wide, based on daily higher high water (HHW) and system hypsometry. Physics-based regression models are used to represent HHW along the system as a function of river inflow, tides, and coastal processes, and hypsometry is used to estimate the associated SWHA. Scenario modeling is employed to attribute SWHA losses to levees, flow regulation, diversion, navigational development, and climate-induced hydrologic change, for subsidence scenarios of up to 2 m, and for 0.5 m fill. For zero subsidence, the system-wide annual-average loss of SWHA is 55 ± 5%, or 51 × 10 5 ha/year; levees have caused the largest decrease ( $${54}_{-14}^{+5}$$ 54 - 14 + 5 %, or ~ 50 × 10 5 ha/year). The loss in SWHA due to operation of the hydropower system is small, but spatially and seasonally variable. During the spring freshet critical to juvenile salmonids, the total SWHA loss was $${63}_{-3}^{+2}$$ 63 - 3 + 2 %, with the hydropower system causing losses of 5–16% (depending on subsidence). Climate change and navigation have caused SWHA losses of $${5}_{-5}^{+16}$$ 5 - 5 + 16 % and $${4}_{-6}^{+14}$$ 4 - 6 + 14 %, respectively, but with high spatial variability; irrigation impacts have been small. Uncertain subsidence causes most of the uncertainty in estimates; the sum of the individual factors exceeds the total loss, because factors interact. Any factor that reduces mean or peak flows (reservoirs, diversion, and climate change) or alters tides and along-channel slope (navigation) becomes more impactful as assumed historical elevations are increased to account for subsidence, while levees matter less.

59 BASIC BIOLOGICAL SCIENCES↗

SIT 45: An interacting, compact, and star-forming isolated galaxy triplet

The underlying scenario of the formation and evolution of galaxy triplets is still uncertain. Mergers of galaxies in isolated triplets give us the opportunity to study the already complex merging process, with minimal contamination of other environmental effects that potentially allow and accelerate galaxy transitions from active star-forming to passive galaxies. The merging system SIT 45 (UGC 12589) is one of 315 systems in the SDSS-based catalogue of Isolated Triplets (SIT); it is an unusual isolated galaxy triplet, consisting of three merging late-type galaxies. The main aims of this work are to study the dynamical evolution and star formation history (SFH) of SIT 45, as well as its dependence on its local and large-scale environment. Methods. To study its dynamics, parameters such as the velocity dispersion (σ v ), the harmonic radius (R H ), the crossing time (H 0 t c ), and the virial mass (M vir ), along with the compactness of the triplet (S) were considered. To investigate the possible dependence of these dynamical parameters on the environment, the tidal force Q parameters (both local and large-scale) and the projected local density (η k ) were used. To constrain the SFH, we used CIGALE to fit its observed spectral energy distribution using multiwavelength data from the ultraviolet to the infrared. SIT 45 is one of the most compact triplets in the SIT, and it is also more compact than triplets in other samples. According to its SFH, SIT 45 presents star formation, where the galaxies also present recent (~200 Myr) star formation increase, indicating that this activity may have been triggered by the interaction. Its dynamical configuration suggests that the system is highly evolved in comparison to the SIT. However, this is not expected for systems composed of star-forming late-type galaxies, based on observations in compact groups. We conclude that SIT 45 is a system of three interacting galaxies that are evolving within the same dark matter halo; its compact configuration is a consequence of the ongoing interaction rather than being due to a long-term evolution (as suggested from its H 0 t c value). We consider two scenarios for the present configuration of the triplet, one where one of the members is a tidal galaxy, and another where this galaxy arrives to the system after the interaction. Both scenarios need further exploration. The isolated triplet SIT 45 is therefore an ideal system to study short timescale mechanisms (~10 8 yr), such as starbursts triggered by interactions which are more frequent at higher redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparative Performance of Gaussian Plume and Backward Lagrangian Stochastic Models for Near-Field Methane Emission Estimation Using a Single Controlled Release Experiment

Methane (CH 4 ) is a major component of natural gas and a potent greenhouse gas. Increasing atmospheric methane concentrations are attributed to emissive anthropogenic activities by an average of 13 ppb per yr since 2020 and are linked to a changing global climate. Mitigating CH 4 emissions from oil and gas production sites has recently become a target to reduce overall greenhouse gas emissions; however, monitoring the efficacy of mitigation strategies depends on accurate quantification of CH 4 emissions at the facility-level. Near-field quantification of methane (CH 4 ) emissions from oil and gas (O&G) facilities remains challenging due to the effects of atmospheric variability and sensor configuration on atmospheric dispersion models. This study evaluates the performance of two atmospheric dispersion models, the Gaussian plume (GP) and backward Lagrangian stochastic (bLS), by comparing calculated CH 4 emissions to controlled single-point emissions between 0.4 and 5.2 kg CH 4 h −1 . Emissions were calculated by both models using 121 individual sets of measurements comprising five-minute averaged downwind methane mixing ratios and matching meteorological data. The comparison shows that the bLS approach achieved a higher proportion of emission estimates within a factor of two (FAC2) of the known emission rates compared to the GP approach. The emissions calculated by the bLS model also had a lower multiplicative error and reduced bias relative to GP. Other error-based metrics further confirmed the bLS model performed better, as it yielded lower RMSE and MAE than GP. Statistical analysis of the emission data shows that the lateral and vertical alignment of the source and the sensor plays a critical role in emission estimations, as measurements made closer to the plume centerline and at a distance between 40 and 80 m downwind yielded the best FAC2 agreement. High wind meander degraded the ability of both approaches to generate representative emissions, particularly with the GP approach, as it violates the modeling approach’s assumption of steady-state emissions. Data suggest emissions calculated by the bLS model are comprehensively in better agreement, but the computational demands of the modeling approach and integration into fenceline systems limit real-time applicability. While these results provide insight into model performance under controlled near-field conditions, their applicability to more complex or heterogeneous oil and gas production environments (e.g., the regions Marcellus or Unita Basins) remains limited and uncertain.

gaussian plume↗

Comparison of aircraft measurements during GoAmazon2014/5 and ACRIDICON-CHUVA

The indirect effect of atmospheric aerosol particles on the Earth's radiation balance remains one of the most uncertain components affecting climate change throughout the industrial period. The large uncertainty is partly due to the incomplete understanding of aerosol–cloud interactions. One objective of the GoAmazon2014/5 and the ACRIDICON (Aerosol, Cloud, Precipitation, and Radiation Interactions and Dynamics of Convective Cloud Systems)-CHUVA (Cloud Processes of the Main Precipitation Systems in Brazil) projects was to understand the influence of emissions from the tropical megacity of Manaus (Brazil) on the surrounding atmospheric environment of the rainforest and to investigate its role in the life cycle of convective clouds. During one of the intensive observation periods (IOPs) in the dry season from 1 September to 10 October 2014, comprehensive measurements of trace gases and aerosol properties were carried out at several ground sites. In a coordinated way, the advanced suites of sophisticated in situ instruments were deployed aboard both the US Department of Energy Gulfstream-1 (G1) aircraft and the German High Altitude and Long-Range Research Aircraft (HALO) during three coordinated flights on 9 and 21 September and 1 October. Here, we report on the comparison of measurements collected by the two aircraft during these three flights. Such comparisons are challenging but essential for assessing the data quality from the individual platforms and quantifying their uncertainty sources. Similar instruments mounted on the G1 and HALO collected vertical profile measurements of aerosol particle number concentrations and size distribution, cloud condensation nuclei concentrations, ozone and carbon monoxide mixing ratios, cloud droplet size distributions, and downward solar irradiance. We find that the above measurements from the two aircraft agreed within the measurement uncertainties. The relative fraction of the aerosol chemical composition measured by instruments on HALO agreed with the corresponding G1 data, although the total mass loadings only have a good agreement at high altitudes. Furthermore, possible causes of the discrepancies between measurements on the G1 and HALO are examined in this paper. Based on these results, criteria for meaningful aircraft measurement comparisons are discussed.

47 OTHER INSTRUMENTATION↗

Learning with Adaptive Conservativeness for Distributionally Robust Optimization: Incentive Design for Voltage Regulation: Preprint

Information asymmetry between the Distribution System Operator (DSO) and Distributed Energy Resource Aggregators (DERAs) obstructs designing effective incentives for voltage regulation. To capture this effect, we employ a Stackelberg game-theoretic framework, where the DSO seeks to overcome the information asymmetry and refine its incentive strategies by learning from DERA behavior over multiple iterations. We introduce a model-based online learning algorithm for the DSO, aimed at inferring the relationship between incentives and DERA responses. Given the uncertain nature of these responses, we also propose a distributionally robust incentive design model to control the probability of voltage regulation failure and then reformulate it into a convex problem. This model allows the DSO to periodically revise distribution assumptions on uncertain parameters in the decision model of the DERA. Finally, we present a gradient-based method that permits the DSO to adaptively modify its conservativeness level, measured by the size of a Wasserstein metric-based ambiguity set, according to historical voltage regulation performance. The effectiveness of our proposed method is demonstrated through numerical experiments.

distribution system operator↗

Characterizing the Potential for Sustainable Azelaic Acid Production from High-Oleic Vegetable Oil Using Two-Step Oxidative Cleavage

Azelaic acid is a renewable monomer conventionally produced via the energy-intensive ozonolysis of oleic acid. Recent advancements have enabled the use of high-oleic vegetable oils (rather than tallow-derived oleic acid) and replaced ozonolysis with two-step oxidative cleavage using hydrogen and oxygen. Although this shift would improve process safety, the financial viability and environmental implications remain uncertain. In this study, we characterized the sustainability of azelaic acid production from high-oleic vegetable oil using two-step oxidative cleavage. Process design, simulation, technoeconomic analysis (TEA), and life cycle assessment (LCA) were executed under uncertainty using BioSTEAM. The modeled system produces azelaic acid at a market-competitive minimum selling price (MSP) of 8.32 [4.93−13.34] $\$kg$ −1 (median 5th−95th percentiles), below the minimum estimated market price of 9.93 $\$kg$ −1 . Further, it has the potential to approach carbon neutrality (0.0 [−5.5 to 5.6] kg of CO 2 -eq kg −1 ) under displacement allocation. Improvements to dihydroxylation (86 to 99%) and oxidative cleavage conversions (93 to 99%) would reduce MSP to $\$5.24$ kg −1 and carbon intensity to −1.90 kg of CO 2 -eq kg −1 (displacement). Additionally, increasing the feedstock triolein content (75 to 85%) lowers MSP by $\$0.82$ kg −1 . Overall, this research demonstrates the potential for financially viable production of azelaic acid from vegetable oils and the utility of agile TEA/LCA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Radial-velocity Search for Binary RR Lyrae Variables

We report 272 radial velocities for 19 RR Lyrae variables. For most of the stars we have radial velocities for the complete pulsation cycle. These data are used to determine robust center-of-mass radial velocities that have been compared to values from the literature in a search for evidence of binary systems. Center-of-mass velocities were determined for each star using Fourier Series and template fits to the radial velocities. Our center-of-mass velocities have uncertainties from ±0.16 km s{sup −1} to ±2.5 km s{sup −1}, with a mean uncertainty of ±0.92 km s{sup −1}. We combined our center-of-mass velocities with values from the literature to look for deviations from the mean center-of-mass velocity of each star. Fifteen RR Lyrae show no evidence of binary motion (BK And, CI And, Z CVn, DM Cyg, BK Dra, RR Gem, XX Hya, SZ Leo, BX Leo, TT Lyn, CN Lyr, TU Per, U Tri, RV UMa, and AV Vir). In most cases this conclusion is reached due to the sporadic sampling of the center-of-mass velocities over time. Three RR Lyrae show suspicious variation in the center-of-mass velocities that may indicate binary motion but do not prove it (SS Leo, ST Leo, and AO Peg). TU UMa was observed by us near a predicted periastron passage (at 0.14 in orbital phase) but the absence of additional center-of-mass velocities near periastron makes the binary detection, based on radial velocities alone, uncertain. Two stars in our sample show Hγ emission in phases 0.9–1.0: SS Leo and TU UMa.

47 OTHER INSTRUMENTATION↗

Understanding the plasticity contributions during laser-shock loading and spall failure of Cu microstructures at the atomic scales

A hybrid atomic-scale and continuum modeling framework is used to study the microstructural evolution during the laser-induced shock deformation and failure (spallation) of copper microstructures. A continuum two-temperature model (TTM) is used to account for the interaction of Cu atoms with a laser in molecular dynamics (MD) simulations. The MD-TTM simulations study the effect of laser loading conditions (laser fluence) on the microstructure (defects) evolution during various stages of shock wave propagation, reflection, and interaction in single-crystal (sc) Cu systems. In addition, the role of the microstructure is investigated by comparing the defect evolution and spall response of sc-Cu and nanocrystalline (nc) Cu systems. The defect (stacking faults and twin faults) evolution behavior in the metal at various times is further characterized using virtual in situ selected area electron diffraction and x-ray diffraction during various stages of evolution of microstructure. Here, the simulations elucidate the uncertain relation between spall strength and strain-rate and the much stronger relation between the spall strength and the temperatures generated due to laser shock loading for the small Cu sample dimensions considered here.

36 MATERIALS SCIENCE↗

Time-of-Use and Demand Charge Battery Controller Using Stochastic Model Predictive Control

Stationary batteries in residential and commercial buildings are often used to smooth customer load profiles and to lower customer electricity bills. Controllers for these battery systems should account for customer energy consumption, rate structures, and high internal battery temperatures, which can lead to reduced performance over the battery lifetime. It is important to consider the uncertainty in forecasting energy consumption and temperature, especially for customers with highly variable and uncertain loads. We propose a novel battery controller using stochastic model predictive control that accounts for these uncertainties and can handle complex rate structures, including demand charges. We show that the controller performs better than standard model predictive control when there is significant uncertainty in the forecast. We also show improvements in the performance with more accurate forecasts and with a more aggressive control strategy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗