Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “utility rate analysis”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Genetic and Structural Analysis of SARS-CoV-2 Spike Protein for Universal Epitope Selection

Evaluation of immunogenic epitopes for universal vaccine development in the face of ongoing SARS-CoV-2 evolution remains a challenge. Herein, we investigate the genetic and structural conservation of an immunogenically relevant epitope (C662–C671) of spike (S) protein across SARS-CoV-2 variants to determine its potential utility as a broad-spectrum vaccine candidate against coronavirus diseases. Comparative sequence analysis, structural assessment, and molecular dynamics simulations of C662–C671 epitope were performed. Mathematical tools were employed to determine its mutational cost. We found that the amino acid sequence of C662–C671 epitope is entirely conserved across the observed major variants of SARS-CoV-2 in addition to SARS-CoV. Its conformation and accessibility are predicted to be conserved, even in the highly mutated Omicron variant. Costly mutational rate in the context of energy expenditure in genome replication and translation can explain this strict conservation. These observations may herald an approach to developing vaccine candidates for universal protection against emergent variants of coronavirus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The Distribution of U.S. Electric Utility Revenue Decoupling Rate Impacts from 2005 to 2017

Investments in energy efficiency and distributed generation reduce electric utility retail sales. Since electric utilities have historically collected a large portion of revenues from volumetric energy rates (¢/kWh), such reductions in sales can impact the utility’s ability to sufficiently recover non-production costs. Regulatory mechanisms that “decouple” utility revenues from sales were implemented to help make the utility indifferent to energy efficiency and distributed generation by ensuring the utility is able to collect an allowed level of revenue each year regardless of its sales. However, noticeable and consistent surcharges over time may create the perception of an incorrectly designed or implemented decoupling mechanism. To date, there are limited quantitative analyses of the rate adjustments due to decoupling mechanisms implemented among U.S. electric utilities (see Morgan, 2013). Given the recent rapid increase in distributed energy resource (DER) adoption in some states and utility service territories, consumers may be facing more prevalent and ongoing decoupling surcharges. This, in turn, may undermine stakeholder support for implementation of decoupling mechanisms and the associated utility support for energy efficiency and DERs. To determine the size of retail rate adjustments associated with decoupling mechanisms and whether they have a tendency towards bill surcharges or credits, we analyzed a large dataset of historical annual decoupling rate adjustments for 21 electric utilities in 11 states between 2005 and 2017. We found that decoupling mechanisms adjusted rates, both up and down, between rate cases, and the majority of those adjustments (54 percent) are small (within a range of -1 to 1 percent). However, we also found that 64 percent of the rate adjustment observations in our sample showed a positive rate adjustment. Importantly, once a surcharge is applied there is an 86 percent chance that there will be a surcharge in the next year as well. While our analysis did not seek to understand the root causes for such results, some possible factors include the accuracy of revenue requirement forecasts, emerging structural changes in customer use and production of energy, misaligned financial motivation, and other factors that influence sales (e.g., economic recession).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Real-time data reduction at 100 Tbps: Challenge and opportunity for AI-based data reduction for next-generation large-scale nuclear physics collider experiment

The modern large-scale nuclear physics (NP) experiments in high-energy particle colliders utilize streaming-readout electronics to digitize detector response at O(100) Tbps bandwidth. Prominent examples at Brookhaven National Lab (BNL) include the sPHENIX experiment at Relativistic Heavy Ion Collider (RHIC), which is under construction, and the experiments proposed for the Electron-Ion Collider (EIC), planned for the 2030s . One of the main challenges for these streaming readout systems is to manage the data rate with sufficient data reduction in real time so the end-data fit persistent storage for offline analysis, which is typically at O(1000) times smaller and O(100) Gbps. Such data reduction traditionally is achieved via real-time high level triggers, which select and save a small subset of collisions of interest. Although triggering is applicable to high energy collider experiments such as those at the Large Hardron Collider at CERN, it is insufficient for these nuclear physics experiments which study diverse collision topologies. And traditional triggering approach is inefficient to preserve the max information harvested from the operation of colliders that costs O(100)M per year to DOE. Meanwhile, in recent years, ML-based high-throughput data reduction has emerged as a promising approach to efficiently preserve max information for a given space of persistent storage, e.g. via AI data compression, feature extraction, and noise filtering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Unveiling Spatial and Temporal Dynamics of Plasmon-Enhanced Localized Fields in Metallic Nanoframes through Ultrafast Electron Microscopy

Plasmonic nanomaterials, particularly noble metal nanoframes (NFs), are important for applications such as catalysis, biosensing, and energy harvesting due to their ability to enhance localized electric fields and atomic efficiency via localized surface plasmon resonance (LSPR). Yet the fundamental structure-function relationships and plasmonic dynamics of the NFS are difficult to study experimentally and thus far rely predominately on computational methodologies, limiting their utilization. This study leverages the capabilities of ultrafast electron microscopy (UEM), specifically photon-induced near-field electron microscopy (PINEM), to probe the light-matter interactions within plasmonic NF structures. Here, the effects of shape, size, and plasmonic coupling of Pt@Au core-shell NFs on spatial and temporal characteristics of plasmon-enhanced localized electric fields are explored. Importantly, time-resolved PINEM analysis reveals that the plasmonic fields around hexagonal NF prisms exhibit a spatially dependent excitation and decay rate, indicating a nuanced interplay between the spatial geometry of the NF and the temporal evolution of the localized electric field. These results and observations uncover nanophotonic energy transfer dynamics in NFs and highlight their potential for applications in biosensing and photocatalysis.

Plasmonics↗
DOE OSTI2022

Additional Analysis of Carbon Capture at Industrial Facilities

This poster provides summary cost details for a technoeconomic analysis of CO2 capture from flue gas generated by cement plants, conducted by U.S. Department of Energy National Energy Technology Laboratory and presented at the 2022 Carbon Management Project Review Meeting held in Pittsburgh, PA (August 15-19, 2022). A range of modern cement plant configurations was considered, along with different fuel types (coal, coke, and natural gas) and two levels of CO2 capture (95% and 99%). The sensitivity of capture cost to a range of modeling assumptions (such as fuel price, financing assumptions, utilization rate, and plant size) was also conducted, to show possible cost ranges that could be expected for decarbonization of typical cement plants in the U.S.

Hughes, Sydney↗
DOE OSTI2023

Techno-economic life cycle assessment of CO 2 -EOR operations towards net negative emissions at farnsworth field unit

Optimizations of CO 2 Water Alternating Gas(WAG)- systems with multi-objectives of incremental recovery and maximization of CO 2 storage are challenging. Here, the incorporation of a total Greenhouse gas (GHG) life cycle assessment is mostly ignored leading to inaccurate estimation of overall net carbon emissions of their operations. In this study, the effect of a total GHG life cycle assessment on a multi-objective CO 2 -WAG optimization with integrated techno-economic assessment (TEA) which factors carbon tax credit is conducted. A life cycle assessment (LCA) was conducted utilizing a 20 -year optimized post history matched data from a high fidelity reservoir simulation model. Using data generated from the optimum result, a techno-economic life cycle analysis was further conducted. The first scenario classified as the base model had an estimated 81% of purchased CO 2 sequestered. The results through a comprehensive techno-economic LCA model yielded a net estimate of 73% of purchased CO 2 . The optimized forecasted model which considered key operational and reservoir factors such as WAG ratio, injection rates and periods, and well specification resulted in an improved sequestration of 92% of purchased CO 2 . However, this also dropped to 84% after taking it through LCA. These results clearly indicate a significant amount of net CO 2 is not accounted for when operations are not analyzed through LCA. From the LCA, direct flaring volumes of CO 2 , energy consumption and efficiency of unit equipment were noticed to be the major causes of these reductions. Considering ten main sources of energy as source of energy generation, a comparative techno-eco LCA was conducted. The results confirmed a lower net volume and NPV for energy sources with higher carbon footprints and vice versa. Thus, a total LCA of CO 2 -WAG greatly influences net storage factor of purchased CO 2 and hence project NPV where tax credit/incentives per ton of CO 2 sequestered is considered. Although operational conditions are optimized for best results, there are significant factors that leads to minimization of net storage factor. This study therefore provides an insightful information for optimizing CO 2 -WAG multi-objectives to achieve minimum GHG emission.

02 PETROLEUM↗
DOE OSTI2024

Analysis of neoclassical tearing mode stabilization experiment by electron cyclotron injection in KSTAR

We report the neoclassical tearing mode (NTM) stabilization experimental results by injection of the electron cyclotron (EC) beam in KSTAR, and its analysis with integrated modelling. In the KSTAR experiment, NTM was intentionally induced and the EC beam was injected by the plasma control system (PCS) to stabilize it. Here, the EC angle was controlled based on the minimum seeking algorithm for the island growth rate. Although the EC angle was not changed significantly, the island width significantly decreased towards the end of the experiment followed by the recovery of plasma performance. To assess the NTM stabilization experiment from a stability perspective, we performed integrated modelling with TRIASSIC utilizing the reconstructed equilibrium with magnetic diagnostics and motional stark effect diagnostics. As a result, we observed that the NTM is stabilized by EC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Boride-Carbon Hybrid Technology to Produce Ultra-Wear and Corrosion Resistant Surfaces for Applications in Harsh Conditions (Final Technical Report)

Engineered functional surfaces play an important role to enable new products and manufacturing processes that can endure harsh service conditions such as high impact and contact loads, highly abrasive wear, extreme temperatures, and corrosive environments. Engineered surfaces can also be instrumental to improving the efficient use of energy by reducing frictional losses and extending service life. The main objective of this project was to develop a hybrid surface engineering technology that combines the advantages of a novel ultrafast boriding process with the next generation of superhard carbon coatings. The hypothesis was that this hybrid process will offer an unprecedented combination of wear and corrosion resistance, low frictional losses and affordability for treated parts so that it can be utilized in many applications. During this project, a duplex process was developed that combines the advantages of ultra-fast electrochemical boriding with those of hard tetrahedral amorphous carbon coatings. Both technologies can be combined to form a hybrid technology that is characterized by low friction and wear properties combined with corrosion and fatigue resistance. Good adhesion of both layers to each other was one main goal of this project, that has been achieved with HF1 adhesion through the Rockwell-C adhesion test. In this project, the mechanical properties of the hybrid coating were modeled through a finite-element analysis approach. We can conclude that the FEA model resembles the actual samples and be utilized to predict mechanical behavior under impacts. Based on this model, application-oriented load conditions can be simulated for optimal layer design regarding thickness and mechanical properties. To exemplify, one conclusion that can be drawn from the nanoindentation model is a boride layer thickness of 50 µm is sufficient to effectively support the carbon coating on the identified AISI 1045 low carbon steel substrate material. The duplex treatment yields wear rates as low as 6 x 10 -8 mm 3 N -1 m -1 and a coefficient of friction of 0.14 when tested against a steel counter face in a ball-on-disk test setup. On the other hand, the wear rate of the only-borided AISI 1045 steel was 5 x 10 -5 mm 3 N -1 m -1 , about three orders of magnitude higher than the duplex coating. At the same time, duplex treated samples experience corrosion resistance, which could not be achieved with single-layer carbon coatings. The developed surface treatment withstands a 3-hour exposure to 15% HCl, while the only carbon coated counter sample shows severe delamination of the coating due to pin hole corrosion. The boride layer is chemically stable and pin hole free because it is formed through an electrochemical process under high current densities (700 mA/cm 2 ) and high temperature. Additionally, the hybrid coating led to at least 3x increase in fatigue strength of the steel substrate, which exceeds the target performance of 30% improvement. There are numerous potential applications for the duplex coatings. A representative application is bearing ball coatings for off-shore windmills. Compared to currently employed surface technologies in this field, the initial costs of applying our technology might be higher due to more process steps but the performance benefits lead to an increased life time of treated parts, which will lower the maintenance and replacement costs in the long-term. To validate the technology for this specific application, the team is currently investigating the process of white etching crack initiation of the duplex coating in collaboration with ANL. Overall, this project successfully validated that the boride-carbon hybrid technology can withstand harsh conditions. One possible approach to commercialization under consideration is to transfer the technology to a startup or an existing coatings company.

36 MATERIALS SCIENCE↗
DOE OSTI2024

2024 Geothermal Rising Conference Poster

As part of a U.S. DOE Geothermal Technologies Office funding opportunity, Geysers Power Company, LLC, an indirect subsidiary of Calpine Corporation, partnered with Sandia National Labs, EGI at the University of Utah, and Texas A&M University to demonstrate increased drilling performance at The Geysers Geothermal Field. The performance target in the drilling demonstrations is at least a 25% improvement in rates of penetration with increased time on bottom for each bit. The planned drilling demonstrations are being conducted as part of an existing drilling campaign intended to enhance reservoir utilization. A major goal of the project is to assess the effectiveness of implementing mechanical specific energy (MSE) and drilling dysfunction diagnosis and remediation in these challenging environments, as well as alternate bit technologies. The first demonstration well has been completed, with a total of 15 PDC bit runs in the 17.5”, 12.25” and 8.5” sections. Initial analysis shows ROP gains in all three hole sections, especially in the 17.5” and 12.25” sections, compared with conventional roller cone bit runs in the demonstration well and offset wells. However in the 8.5” hole, wear and damage to the PDC bits resulted in relatively short bit runs. Analysis is underway to take advantage of the positive results and remediate the challenges.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

OH Roaming and Beyond in the Unimolecular Decay of the Methyl-Ethyl-Substituted Criegee Intermediate: Observations and Predictions

Alkene ozonolysis generates short-lived Criegee intermediates that are a significant source of hydroxyl (OH) radicals. Here this study demonstrates that roaming of the separating OH radicals can yield alternate hydroxycarbonyl products, thereby reducing the OH yield. Specifically, hydroxybutanone has been detected as a stable product arising from roaming in the unimolecular decay of the methyl-ethyl-substituted Criegee intermediate (MECI) under thermal flow cell conditions. The dynamical features of this novel multistage dissociation plus a roaming unimolecular decay process have also been examined with ab initio kinetics calculations. Experimentally, hydroxybutanone isomers are distinguished from the isomeric MECI by their higher ionization threshold and distinctive photoionization spectra. Moreover, the exponential rise of the hydroxybutanone kinetic time profile matches that for the unimolecular decay of MECI. A weaker methyl vinyl ketone (MVK) photoionization signal is also attributed to OH roaming. Complementary multireference electronic structure calculations have been utilized to map the unimolecular decay pathways for MECI, starting with 1,4 H atom transfer from a methyl or methylene group to the terminal oxygen, followed by roaming of the separating OH and butanonyl radicals in the long-range region of the potential. Roaming via reorientation and the addition of OH to the vinyl group of butanonyl is shown to yield hydroxybutanone, and subsequent C–O elongation and H-transfer can lead to MVK. A comprehensive theoretical kinetic analysis has been conducted to evaluate rate constants and branching yields (ca. 10–11%) for thermal unimolecular decay of MECI to conventional and roaming products under laboratory and atmospheric conditions, consistent with the estimated experimental yield (ca. 7%).

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

Reliability-Informed Economic and Energy Evaluation for Bi-Level Design for Remanufacturing: A Case Study of Transmission and Hydraulic Manifold

Abstract Design for remanufacturing (DfRem) is one attractive strategy that encourages the reuse of a product and extends the product's life cycle. Traditional design processes often only consider product reliability at an early design stage. However, from the perspective of environmental sustainability, it is becoming increasingly important to evaluate the long-term economic and environmental impacts of design decisions during early-stage design. We propose a bi-level DfRem framework consisting of system-level reusability allocation and component-level design tradeoff analysis, considering reliability and product warranty policy. First, a system-level reusability allocation problem aims at a theoretical exploration of the design space where all the components comprising the system are allocated certain reuse rates to achieve target energy savings with minimum cost. Following the theoretical exploration at the system level, a component-level analysis looks at practical design options for each component and trades-off between the overall cost and energy consumption for multiple remanufacturing cycles. Both levels of the framework require modeling component reuse for multiple remanufacturing cycles, which we achieve by using a branched power-law model that provides probabilistic scenarios of reusing the component or replacing it with a new part. We demonstrate the utility of this framework with the case study of an infinitely variable transmission (IVT) used by some agricultural machines manufactured by John Deere and show snapshots of a prototype software tool that we developed for easy use by designers.

Engineering↗
DOE OSTI2021

Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning

Various modeling techniques are used to predict the capacity fade of Li-ion batteries. Algebraic reduced-order models, which are inherently interpretable and computationally fast, are ideal for use in battery controllers, technoeconomic models, and multi-objective optimizations. For Li-ion batteries with graphite anodes, solid-electrolyte-interphase (SEI) growth on the graphite surface dominates fade. This fade is often modeled using physically informed equations, such as square-root of time for predicting solvent-diffusion limited SEI growth, and Arrhenius and Tafel-like equations predicting the temperature and state-of-charge rate dependencies. In some cases, completely empirical relationships are proposed. However, statistical validation is rarely conducted to evaluate model optimality, and only a handful of possible models are usually investigated. This article demonstrates a novel procedure for automatically identifying reduced-order degradation models from millions of algorithmically generated equations via bi-level optimization and symbolic regression. Identified models are statistically validated using cross-validation, sensitivity analysis, and uncertainty quantification via bootstrapping. On a LiFePO 4 /Graphite cell calendar aging data set, automatically identified models utilizing square-root, power law, stretched exponential, and sigmoidal functions result in greater accuracy and lower uncertainty than models identified by human experts, and demonstrate that previously known physical relationships can be empirically "rediscovered" using machine learning.

25 ENERGY STORAGE↗
DOE OSTI2021

Optimal Sizing of an Electric Vehicle Charging Station with Integration of PV and Energy Storage

This paper proposes an optimization model for the optimal configuration of an grid-connected electric vehicle (EV) extreme fast charging station considering integration of photovoltaic (PV) and energy storage. The proposed model minimizes the annualized net cost (i.e., maximizes the annualized net profit) of the extreme fast charging station, including investment and maintenance cost of charging ports, PV and energy storage, net cost of purchasing energy from utility and selling energy to EV customers, degradation cost of energy storage and demand charge. The decision variables are number of charging ports, capacity of invested PV and the power and energy ratings of invested energy storage. The Erlang-loss system is adopted to model the EV mobility. Results of numerical simulations indicate that investment of PV and energy storage could increase the annualized profit of the extreme fast charging station. In addition, the impacts of various parameters on the optimal solution are investigated by sensitivity analysis.

Liu, Guodong↗
DOE OSTI2025

SAM Two-Phase Flow Model Development and Applications for Operational Transients in Advanced Reactors

As advanced nuclear technologies continue to develop, the need for the flexible operation and generation of these advanced reactors becomes necessary to maximize economic potential. As large-scale experiments are not always feasible, modeling and simulations of advanced reactors play a crucial role in design optimization and analysis. The SAM (System Analysis Module) code developed at Argonne National Laboratory is a state-of-the-art system-level thermal-hydraulic code aimed at simulating advanced reactor systems. Recent code developments have implemented two-phase flow modeling using the homogeneous equilibrium model, and a new steam generator component has been developed to utilize the two-phase flow implementation. In addition to verification tests, a load-following simulation was performed to model a realistic load-following transient in a proposed integrated system consisting of a conceptual advanced reactor known as the Advanced Burner Test Reactor (ABTR) and thermal energy storage (TES) tanks. The integrated system model uses two large TES tanks designed for sodium and a model helical coil steam generator to simulate the operational load-following transient. The flow rates of the feedwater and secondary loops are regulated to meet a prescribed steam generator load consistent with the electricity demand over a 24-h period. In conclusion, the results found the ABTR system was able to maintain stable reactor conditions and primary- and secondary-side characteristics over the course of the load-following transient.

Advanced Burner Test Reactor (ABTR)↗
DOE OSTI2021

Two-dimensional retrieval methods for ultrafast imaging of molecular structure using laser-induced electron diffraction

Molecular structural retrieval based on electron diffraction has been proposed to determine the atomic positions of molecules with sub-angstrom spatial and femtosecond temporal resolutions. Given its success on small molecular systems, in this work, we point out that the accuracy of structure retrieval is constrained by the availability of a wide range of experimental data in the momentum space in all molecular systems. To mitigate the limitations, for laser-induced electron diffraction, here we retrieve molecular structures using two-dimensional (energy and angle) electron momentum spectra in the laboratory frame for a number of small molecular systems, which have previously been studied with 1D methods. Compared to the conventional single-energy or single-angle analysis, our 2D methods effectively expand the momentum range of the measured data. Besides utilization of the 2D data, two complementary methods are developed for consistency check on the retrieved results. Further, the 2D nature of our methods also offers a way of estimating the error from retrieval, which has never been explored before. Comparing with results from prior experiments, our findings show evidence that our 2D methods outperform the conventional 1D methods. Paving the way to the retrieval of large molecular systems, in which their tunneling ionization rates are challenging to obtain, we estimate the error of using the isotropic model in place of including the orientation-dependent ionization rate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Research Needs and New Capabilities for Retail Electricity Rate Analysis

Retail electricity rates are at the center of supply- and demand-side changes in future power systems and influence system costs and DER adoption decisions, all with important implications for system reliability, resiliency, and energy affordability. Building on decades of research experience and analytical insights at Lawrence Berkeley National Laboratory (LBNL) and the National Renewable Energy Laboratory (NREL), the report identifies key near-term analysis questions, enhancements to existing capabilities, and longer-term new capabilities to provide actionable insights for electricity system decision-makers, including state utility regulators, electric and gas utilities, ratepayer advocates, and DER solution providers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Assessment of frequency and amplitude dependence on the cyclic degradation of polyurethane foams

Many energy absorption applications utilize flexible polymeric foams for their viscoelastic properties. It is desired that the material will perform consistently across repeated compression cycles. This study examines the effect of fatigue at low strain rates on the viscoelasticity of open-cell polyurethane foam. Six polyurethanes of the same base composition with two porosities (70% and 80%) and three chemical indexes (79i, 100i, and 121i) are tested. Large deformation cyclic compression of the foams is conducted on a universal testing system (UTS). These data are then post-processed leveraging dynamic mechanical analysis Fourier transform rheology to characterize changes in the viscoelasticity of the materials over fatigue cycles. Results show that foams can increase or decrease in stiffness up to 10% over 10 4 cycles. Specifically, higher chemical index, higher excitation frequency, and larger excitation amplitude correlate with a more pronounced decrease in stiffness. Damping can also change by 15% and correlates with chemical index and excitation frequency. Consequently, the findings suggest that internal foam structure and bulk material properties as well as applied loading parameters affect the viscoelastic fatigue response of flexible polymeric foams.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Industrial battery operation and utilization in the presence of electrical load uncertainty using Bayesian decision theory

Behind the meter battery storage is becoming increasing popular in all sectors, though enthusiasm has recently lagged in the industrial sector. Even though there may be many factors contributing to this including lack of innovation, prohibitive costs, and undesirable rate structures, a difficulty arises in accounting for uncertainty of electrical load in industrial facilities while still attempting to utilize battery storage as much as possible all while trying to achieve fiscal profitability. Here this study utilizes Gaussian process regression and Bayesian decision theory to organize load data and quantify electrical load uncertainty to properly and effectively discharge industrial battery storage. The study employs a simulation model to set battery load setpoints for the span of the utility billing period according to the degree of risk aversion. This combination of economic analysis according to utility billing period and utilization of degree of risk aversion to make decisions on the uncertainty of the data has not before been applied to battery storage. The method resulted in an annual average reduction of peak demand by 3.8 % at the lowest amount of savings and lowest risk aversion. The highest risk aversion resulted in an annual average reduction of peak demand of 7.5 %. The maximum reduction of peak load in any month was 13.8 % in the month of December with a relatively high risk aversion. With a the highest amount risk aversion tested, the model reduced demand ten of the twelve months of the year.

25 ENERGY STORAGE↗