Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex loading”

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 757 records · Page 42

A Markov framework for generalized post-event systems recovery modeling: From single to multihazards

State-dependent models can be used to represent the system recovery process as a series of stochastic transitions from lower to higher functional states. However, the applications of these models have been limited in scope and there is a lack of a generalized recovery modeling framework. A generalized framework would permit a robust forecasting of systems and system-of-systems recovery under multiple hazards, and more broadly, would contribute to community disaster preparedness. This paper develops a generalized post hazard-event recovery modeling framework based on state-dependent Markov-type processes. We then apply the proposed framework to solve a spectrum of problems that range from hind-casting single-system recovery following a single hazard event to forecasting post-event trajectories under multiple hazards and modeling the recovery of a system-of-systems. First, Markov chains are used to hind-cast the observed recovery for a portfolio of buildings affected by the 2014 South Napa, California, earthquake. Next, Markov processes are used to formulate a parametric post hazard-event recovery model, which can be updated using Bayesian statistics when relevant datasets become available. Semi-Markov processes are then used to develop a more general model of single hazard recovery, which accounts for the intensity of the loading and level of damage caused by the event. Semi-Markov processes with non-renewal features are then used to account for multihazard interactions in a post-event recovery model, and applied to a case study that involves a community in Charleston, South Carolina. Lastly, Markov-type processes are combined with Bayesian networks to model the recovery of residential, commercial, educational, and industrial buildings (system-of-systems) following a hazard event. Overall, these applications demonstrate the versatility of the Markov framework towards handling recovery problems with varying levels of complexity.

42 ENGINEERING↗

Synergistic Solvent-Surface Interactions Enable Alkyne Semihydrogenation at Palladium

Enabling higher yield and better selectivity for fine-chemical synthesis through heterogeneous catalysis is intricately linked to the interplay of active sites, reaction conditions, and mass transfer influence provided by the catalyst. Alkyne semihydrogenation is ubiquitous in the production of bulk chemicals in the pharmaceutical, polymer, or fine-chemical industries, but product selectivity remains a major challenge. Here, in this study, we demonstrate that the design of catalysts encompassing nickel (Ni) foams as contiguous monolith supports, decorated with ultralow loading of Pd/PdO x nanoparticles on a carbonized polydopamine interface and tuned with a thin layer of Al 2 O 3 , in conjunction with an optimized reaction environment leads to highly selective alkyne semihydrogenation. The reactions demonstrate good functional group tolerance and applicability to flow reactor systems. Combined computational and experimental studies are presented to describe the synergistic effect between the solvent-surface interaction and the degree of Pd surface reduction that are necessary to promote this selectivity. The system highlights the opportunity for catalyst-solvent codesign as a benign alternative to more complex reactants featuring extrinsic poisons or less-favored dopants.

atomic layer deposition↗

Scalable Deep Learning-Based Microarchitecture Simulation on GPUs

Cycle-accurate microarchitecture simulators are essential tools for designers to architect, estimate, optimize, and manufacture new processors that meet specific design expectations. However, conventional simulators based on discrete-event methods often require an exceedingly long time-to-solution for the simulation of applications and architectures at full complexity and scale. Given the excitement around wielding the machine learning (ML) hammer to tackle various architecture problems, there have been attempts to employ ML to perform architecture simulations, such as Ithemal and SimNet. However, the direct application of existing ML approaches to architecture simulation may be even slower due to overwhelming memory traffic and stringent sequential computation logic. This work proposes the first graphics processing unit (GPU)-based microarchitecture simulator that fully unleashes the potential of GPUs to accelerate state-of-the-art ML-based simulators. First, considering the application traces are loaded from central processing unit (CPU) to GPU for simulation, we introduce various designs to reduce the data movement cost between CPUs and GPUs. Second, we propose a parallel simulation paradigm that partitions the application trace into sub-traces to simulate them in parallel with rigorous error analysis and effective error correction mechanisms. Combined, this scalable GPU-based simulator outperforms by orders of magnitude the traditional CPU-based simulators and the state-of-the-art ML-based simulators, i.e., SimNet and Ithemal.

97 MATHEMATICS AND COMPUTING↗

Evaluation of Structural Lattices for a Davis Gun Earth Penetrator Impact Experiment

The advanced materials team investigated the use of additively manufactured metallic lattice structures for mitigating impact response in a Davis gun earth penetrator impact experiment. High-fidelity finite element models were developed and validated with quasistatic experiments. These models were then used to simulate the response of such lattices when subjected to the acceleration loads expected in the Davis gun experiment. Results reveal how the impact mitigation performance of lattices can change drastically at a certain relative density. Based on these observations, an experiment deck was designed to probe the response of lattices with different relative densities during the Davis gun phase 2 shots. The expected performance of these lattices is predicted before testing based on simulation results. The results of the Davis gun phase 2 shots are expected to provide data which will be used to assess the predictive capability of the finite element simulations in such a complex impact environment.

36 MATERIALS SCIENCE↗

Consolidated Hydropower Data Repository: Value and Opportunities

Hydropower is one of several types of generating assets that provides energy, capacity, and services to electric power systems. It does so under rubrics and objectives—market driven and regulated, internal and external to asset and fleet owners—that address reliability, cost, price, and, increasingly, flexibility of output. The aggregation of data from multiple hydropower units can provide insights into asset operations and maintenance practices and needs and assist in meeting hydropower objectives. This paper examines the concept and potential benefits of aggregating hydropower asset data—primarily supervisory control and data acquisition (SCADA) information—with examples of insights developed from data aggregated by the Hydropower Research Institute (HRI). Data aggregation as discussed herein, and as implemented by the HRI, extends beyond multiple units in a powerhouse and beyond multiple hydropower facilities in an electric utility fleet or river system. Examples of research and analytics from such aggregated datasets range from unit load dependency analyses to modeling sensor measurements to detect and diagnose anomalies in assets. These examples showed the benefits of utilizing the entire dataset for insights into how the sensor layout of a single unit or set of units compares to the hydropower industry overall. Such insights include whether additional sensors are needed to complete analyses or to make decisions. In addition, utilizing multiple sensors of the same kind within a unit can provide an indication of possible current or upcoming problems with equipment. Although other analyses are possible, their use requires the development of complex models and, potentially, access to types of data that are currently not included with the example dataset used in this study. However, the examples studied herein confirmed the value of data aggregation in the fleet and unit contexts, and the value extends beyond multiple units in a powerhouse and beyond multiple hydropower facilities in an electric utility fleet or river system. The assessments also provided insights into potential extensions to the data aggregation concept that could further add to their value to the hydropower community, and these are included in this document as a set of recommendations.

13 HYDRO ENERGY↗

Large Area Detection of Microstructural Defects with Multi-Mode Ultrasonic Signals

Cyclic loading or other stresses can lead to development of cracks and crack growth in mechanical structures, leading to eventual failure. While ultrasound imaging can be used for non-destructive testing of such structures, conventional ultrasound techniques are often limited by crack size, density, and areal coverage. An effective characterization of real-world, large-area structures is required at an early damage stage to prevent catastrophic failure and predict remaining life. In this study, a new nonlinear ultrasonic testing (NUT) method is proposed for large-area monitoring of practical structures with arbitrary complexity by using multiple-mode guided-wave ultrasonic signals. The proposed guided-wave NUT technique requires single-element transducers, simple electronics, and a mixed time-frequency domain signal processing. As a proof-of-concept demonstration, numerical simulations and experiments are performed on an A36 carbon steel beam assembly with previously formed microstructural defects that cause nonlinearities in ultrasonic response. The quadratic dependence of the nonlinear wave excitation on the input ultrasonic signal amplitude is shown by numerical simulations, and such a nonlinear ultrasonic response is experimentally observed in the zone with a high density of microstructural defects.

36 MATERIALS SCIENCE↗

Polypropylene Composites Reinforced With Recycled Waste Cellulosic Fiber/Fine Mixture: The Impact of Cellulose Sieving on Performance

This study explores how a sieving step of waste cellulosic fiber and fine (WCFF) mixture affects the performance of WCFF‐loaded polypropylene (PP) composites and whether the separation of fines from fibers offers an added benefit. The WCFF samples were downsized, and four different filler size ranges were sieved using a series of mesh sizes from 4 to 0.85 mm. The WCFF/PP composites were then compounded at 20 wt.% loading of WCFF using a twin‐screw extruder. Incorporating WCFF increased the tensile strength to 41.28 MPa and the modulus to 3207 MPa, accounting for 28% and 38% enhancements, respectively. Interestingly, the greatest improvements were associated with the nonsieved WCFF case, and the sieved WCFF fibers provided only marginal enhancements over virgin PP. The outperformance of nonsieved WCFF was attributed to the synergistic reinforcement of hybrid fibers and fines as well as the maintenance of longer fibers in the system. However, the strain at break and impact strength of PP decreased after introducing WCFF. Moreover, the complex viscosity and storage modulus increased with an increase in the filler size, due to the formation of a more effective percolative network. The PP's crystallinity exhibited a relatively strong dependency on the sieving, where WCFF samples with short‐aspect‐ratio fillers promoted the crystallinity significantly. It was also found that the WCFF degradation onset temperature increased once it was incorporated into PP. This study suggests that waste cellulosic feedstocks can be utilized as a reinforcement without additional sieving to manufacture high‐performance and cost‐effective composites.

36 MATERIALS SCIENCE↗

A Green's Function Wind Turbine Induction Model That Incorporates Complex Inflow Conditions

ABSTRACT In this work, we develop a new analytical turbine induction model that can incorporate complex inflow conditions including cases where the wind velocity and temperature profiles can vary as functions of height. This induction model is derived from the linearized Navier–Stokes and leads to a second‐order ODE that can be solved using a Green's function formulation. The corresponding Green's function for several configurations are found including the infinite domain, semi‐infinite domain with ground plane, and a power law velocity inflow profile. The results of this approach are then compared with simulations of the turbine induction field using the AMR‐Wind CFD solver with a uniformly loaded actuator disk model. These comparisons show that the Green's function approach captures the centerline blockage, three‐dimensional blockage flow field, and streamwise velocity slow down, with very good agreement for lower thrust conditions and at larger distances away from rotor disk. The effects of shear on the turbine blockage were also compared using a power law inflow profile, and we show that this approach matches the CFD predictions for the cases considered.

17 WIND ENERGY↗

Co-optimization of nuclear reactor flexible power operation and maintenance scheduling

As flexible power operation of nuclear power plants becomes more attractive due to the reduction in fossil-fueled dispatchable generation on energy grids, finding optimal power production strategies that balance revenue generation with operational concerns becomes more complex. This article presents a general framework to aid operators in designing economically optimal long term dispatch strategies for nuclear power plants. The principal novelty is the linking of estimated system remaining useable life (RUL) to strategic operational decisions. It is shown that, depending on the relationship between the fixed costs from maintenance and the associated lost revenue from an outage, it can be economically optimal in the long term to delay a maintenance outage and not perform this alongside refueling. For a given relationship between power ramping and degradation, optimal strategies were found that discouraged load following in some situations while minimizing unnecessary maintenance. It is shown that heavy load following can cause maintenance and refueling outages to diverge due to their inverse relationships with respect to load following, potentially leading to a significant loss in capacity factor. As a result, this general framework can be applied to specific reactor dispatch allowing operators to adapt operational strategies as future grid conditions change.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Using FLAC{sup R} Modeling Software to Analyze the Stability of a LLRW Disposal Facility - 20031

Waste Control Specialists (WCS) wished to amend its existing LLRW disposal license such that it might use gravel aggregate instead of concrete as backfill around large component (LC) parts from decommissioned nuclear reactors. A classical mechanics evaluation of the strength of the LCs was presented to the Texas Commission on Environmental Quality (TCEQ), demonstrating that the LCs had adequate capacity to withstand the loads within the facilities, and that concrete was unnecessary for structural support. TCEQ requested that additional analyses of the LCs within the facilities be performed using FLAC{sup R} finite difference modeling software. Because FLAC{sup R} has the capability of modeling the support provided by the lateral earth pressures, the factors of safety against failure for the LCs increased. TCEQ found the FLAC{sup R} analyses presented in this report to be satisfactory, and granted WCS' license amendment request. This work demonstrates that FLAC{sup R} may be used advantageously to more accurately represent the complexities of the physical situation being modeled. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗

A time‐varying formulation of the curled wake model within the FAST.Farm framework

Abstract In this article, we present a time‐varying formulation of the curled wake model that we implemented in FAST.Farm. The curled wake model, originally developed for steady‐state conditions, is used to produce realistic wake profiles behind a wind turbine in yawed (or skewed) conditions. We begin by introducing the key elements of the FAST.Farm framework. Then, after briefly summarizing the original wake dynamics formulation of FAST.Farm based on a polar wake profile, we present the new time‐varying formulation of the curled wake model, compare the two, and highlight the differences with the original curled wake model. After discussing some implementation details, we present different applications with increasing levels of complexity: single turbine with uniform and turbulent inflow, fixed and transient yaw, and multiple turbines. We verify our results using the original FAST.Farm implementation and large‐eddy simulations. The results with the new curled wake model are improved compared to the original implementation, as they include cross‐flow velocities and wake asymmetry. Yet, large‐eddy simulation results show a more pronounced lateral convection of the wake and a stronger concentration of vorticity at the top vortex. The new curled wake implementation in FAST.Farm should enable the calculation of not only generator power but also wind turbine structural loads for applications involving intentional or unintentional skewed flow and wind‐farm control involving wake steering.

17 WIND ENERGY↗

Recent advances and challenges in optimization models for expansion planning of power systems and reliability optimization

Optimization models for expansion planning of power systems aim to determine capacities, investment timing, and location of power systems to satisfy the power demands while minimizing the total cost. The models have become complex in recent years to reflect both regulations on conventional energy sources and the increasing penetration of renewable energy sources (RES). This paper reviews the basic concepts and optimization models for expansion planning of power systems. We first explain the definition and features of generation expansion planning (GEP), transmission expansion planning (TEP), and generation and transmission expansion planning (GTEP). To address the computational challenges of large-scale expansion planning problems, we review several simplifications including temporal and spatial aggregation, and decomposition methods. This paper also addresses power system reliability defined as the probability of satisfying the load demand while withstanding failures of components. Finally, the goal of this paper is to provide a research overview, discuss trends in expansion planning of power systems, and suggest directions for future research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic compression effects of H 2 ⁡O in a dynamic diamond anvil cell: Origin of metastable ice VII and its crystal growth kinetics

We report on the structural verification of metastable ice VII solidifying in the phase space of ice VI at 1.80 GPa at room temperature. Using time-resolved (TR) x-ray diffraction and TR ruby luminescence paired with high-speed microphotography utilizing a dynamic diamond anvil cell, an initial compression rate range from 0.12 to 95.84 GPa/s was explored. The solidification pressure of metastable ice VII has a potential sigmoidal dependence upon compression rate with a turnover compression rate of ∼80 GPa/s. The preferred crystallization of ice VII in the stability field of ice VI is due to the increased nucleation rate of ice VII over ice VI at 1.77 GPa that is driven by the surface energy difference between the liquid and solid phases along with the change in Gibbs free energy of solidification. The dynamic pressure-volume–compression behaviors of ice phases (VI and VII) show a lattice stiffening in both phases, especially during the compression loading. It is also found that the compression rate greatly affects the solid-solid phase transition between ice VI and VII but does not affect the liquid-solid transition between water and ice VI as much. Lastly, a third phase transition was found to occur after metastable ice VII transforms into high-density amorphous (HDA) ice, which could be a disordered hydrogen-bonded network configuration of ice VII forming out of HDA ice facilitated by the decoupling of the oxygen movement and reorientation of the H 2⁡ O molecule. These results demonstrate the complexity of a seemingly simple molecule H 2⁡ O, how it can readily change its static properties with the modification of (de)compression rate, and highlight the need to use multiple TR structural and spectroscopic probes at higher time resolutions to realize the most comprehensive understanding.

Chemical bonding↗

Smart Charging for Electric Ride-Hailing Vehicles using Renewables: A San Francisco Case Study

Charging large fleets of electric ride-hailing vehicles (ERVs) is a complex matter that could serve different objectives: lower carbon dioxide emissions, lower monetary expenditures, or maximize solar photovoltaics (PV) energy consumption. Currently, it is unclear how each of those objectives could impact the business and performance of a ride-hailing fleet. In order to fill this gap, this article employs a dynamic transportation model: a smart charging simulation that combines agent-based, discrete-event, and system dynamic modelling by comparing the above-mentioned objectives in separate scenarios. The results show that each scenario successfully manages to shift between 34% and 87% of all load to hours of the day when the objectives of those scenarios are met. Therefore, in comparison to the baseline, smart charging can save between 5% and 26% of monthly emissions and between 4% and 57% of monthly expenditures. The solar PV scenario, however, results in the highest savings, while ensuring profitable economics via net metering in the short- as well as long term. Finally, the sensitivity analysis points to important trade-offs between several fleet performance metrics. The article concludes by giving business and policy recommendations for maximising the economic, energy and environmental efficiency of large ERV fleets.

ADVANCED PROPULSION SYSTEMS,SOLAR ENERGY↗

Investigation into Reductive Separation of Tc-99 from Low Activity Waste Off-Gas Condensate by Zero Valent Iron Materials - 20311

This study is investigating the potential for reductive separation of technetium-99 (Tc) from off-gas condensate generated from vitrification of legacy tank waste at the Hanford Tank Waste Treatment and Immobilization Plant. Management of Tc is a high-priority activity for the U.S. Department of Energy (DOE) Office of Environmental Management (EM) complex, which must be addressed based on knowledge gaps and technology development needs outlined in the DOE-EM Technetium Management Program Plan [1]. Technetium is one of the main risk drivers in the permanent disposal of low activity waste (LAW) by vitrification at the Hanford Site due to its high volatility and only fractional incorporation into borosilicate glass. Volatilized Tc will be captured by an off-gas treatment system, and the current plan is to continuously recycle generated condensate back to the LAW glass melter feed. Off-gas recycling is effective at increasing Tc loading in the LAW glass, but it also disproportionately increases concentrations of sulfate, halides, and other volatile compounds impeding overall LAW processing leading to significantly greater waste glass volume requirements. Removal of Tc from secondary liquid waste offers an alternative disposal route to continuous off-gas recycling which may result in significant reduction of the volumes of LAW waste glass and substantial cost savings. The objective of this research funded by the DOE-EM Minority Serving Institutions Partnership Program (MSIPP) is to determine the suitability of zero valent iron (ZVI) as a separation technology to remove Tc from the off-gas condensate stream to minimize its recycle. While ZVI has been proposed for various remediation applications with Tc, it has not been studied for its suitability for removing Tc from the LAW off-gas condensate. The goal is to identify the optimal conditions for which Tc{sup (VII)} will be reduced to Tc{sup (IV)} species and incorporated into an iron matrix. This paper presents a programmatic overview of our evaluation of commercial ZVI materials for their efficiency towards Tc{sup (VII)} reductive removal from simplified off-gas simulant solutions under aerobic conditions. Tc concentrations in solutions were monitored for up to 30 days to determine efficiency and kinetics of reductive removal of Tc. The best performing candidate materials were selected for further experiments evaluating the effects of (1) variable ionic strength, (2) the presence of chromium and corrosive anions, and (3) ZVI aging on Tc{sup (VII)} reductive removal in batch experiments. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Counter Data Paucity through Adversarial Invariance Encoding: A Case Study on Modeling Battery Thermal Runaway

Lithium-ion batteries, widely used for their durability and high energy storage, face the risk of internal short circuits leading to catastrophic thermal runaway events. These events, triggered by external stimuli like mechanical loads, pose safety concerns in applications such as electric vehicles. Detecting and understanding thermal runaway events is crucial, but physics-driven models struggle to explain the non-linear evolution of battery temperature during these events, considering factors like material composition and state-of-charge. Due to the rarity of these events and the cost of data collection, we propose a deep learning (DL) model to predict battery temperature responses during thermal runaway. The challenge lies in the scarcity of data, making traditional DL models prone to overfitting and learning low-quality representations of the complex process.Our approach introduces a novel few-shot architecture that incorporates an adversarially governed invariant encoding process. This architecture aims to distill "invariant" relationships by addressing distributional shifts in data across various battery properties, facilitating the detection of thermal runaway events. Specifically, our results demonstrate that deep learning models conditioned on these "invariant" representations outperform state-of-the-art baselines, achieving a remarkable 96.8% performance improvement in terms of the popular metric MAPE. This framework presents a promising direction for enhancing battery safety modeling, particularly in the context of rare and complex events like thermal runaway. Our code and code and dataset used for the paper are public1.

Tabassum, Anika [ORNL] (ORCID:0000000254600955)↗

Site demonstration and performance evaluation of MPC for a large chiller plant with TES for renewable energy integration and grid decarbonization

Thermal energy storage (TES) for a cooling plant is a crucial resource for load flexibility. Traditionally, simple, heuristic control approaches, such as the storage priority control which charges TES during the nighttime and discharges during the daytime, have been widely used in practice, and shown reasonable performance in the past benefiting both the grid and the end-users such as buildings and district energy systems. However, the increasing penetration of renewables changes the situation, exposing the grid to a growing duck curve, which encourages the consumption of more energy in the daytime, and volatile renewable generation which requires dynamic planning. The growing pressure of diminishing greenhouse gas emissions also increases the complexity of cooling TES plant operations as different control strategies may apply to optimize operations for energy cost or carbon emissions. This paper presents a model predictive control (MPC), site demonstration and evaluation results of optimal operation of a chiller plant, TES and behind-meter photovoltaics for a campus-level district cooling system. The MPC was formulated as a mixed-integer linear program for better numerical and control properties. Compared with baseline rule-based controls, the MPC results show reductions of the excess PV power by around 25%, of the greenhouse gas emission by 10%, and of peak electricity demand by 10%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗