Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic 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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater dependent ecosystems in areas with industrial and military land use are at risk of direct exposure to a wide range of contaminants, including PFAS chemicals. Glaciated terrain often has mixed high and low permeability sediments coupled with groundwater flow-through lake features. These hydrogeologic attributes create highly complex ‘source to seep’ dynamics that make spatiotemporal contaminant transport patterns difficult to predict. We investigated one such system in detail using a suite of heat-tracing and chemical methods. Numerous (n=57) preferential groundwater discharge zones (vertical flux rates ranging 0.2 to 3.2 m/d) were identified across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, adjacent to an Air Force Base with several known PFAS source areas. Surface-water and groundwater samples were collected and analyzed (n=145) for precursors and terminal PFAS compounds between March and September 2022. Samples were collected at identified seeps along the Quashnet River (n=59), from wells upgradient from the stream-wetland system (n= 44), from contributing flow-through kettle lakes (n=8), and at multiple locations along the Quashnet River (n=34). Samples from seeps and wells had measured PFAS concentrations ranging from non-detect to approximate 3,500 ng/L (mean= 1,650 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to emergence at the discharge zones. Groundwater-lake interaction along flowpaths that sourced the sampled seeps was farther supported by significant correlations (p < 0.01) between deuterium excess and %PFAS precursors, and between %PFAS precursors and multiple terminal PFAS compounds (e.g., PFPeS, PFBS, PFHxS). However, some sampled seeps contributing groundwater to the stream-wetland system had much higher total PFAS concentrations (>1000 ng/L) than the upgradient kettle lakes, despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. PFAS compounds and water isotopic composition at sampled multilevel groundwater wells, rivers, lakes, and seeps suggest that a complex mixture of source groundwater and flowpath characteristics are responsible for diverse observed PFAS mixtures at preferential discharge zones across the stream-wetland system. Further, total PFAS loading patterns to the Quashnet River via groundwater discharge remained remarkably similar from winter to summer to fall conditions, despite a regional dry period in late summer 2022 with the upper river channel completely drying. This work addresses gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations in controlling contaminant mass loading in preferential groundwater discharge zones and presents a transferrable field toolkit for efficient characterization of spatially preferential PFAS transport dynamics.

Contaminant transport↗
DOE OSTI2020

Umbilical Fatigue Analysis for a Wave Energy Converter

Unique umbilical designs for wave energy converters (WECs), including the ability to handle significantly larger motions and loads over long deployments, are often required when conventional marine umbilical designs for offshore oil and gas and offshore wind may not meet the design and cost needs of wave energy technologies. This study details a fatigue analysis of a dynamic power umbilical attached to a two-body floating point absorber WEC system, using the sea states provided for the PacWave testing facilities. The 6 degrees of freedom motion time history for the WEC was simulated, and the motions of the attachment point for the umbilical on the WEC and respective sea states were used to analyze the dynamic motions and fatigue of the connected power umbilical to predict the fatigue life. The results show that the fatigue damage observed is more significant in shallow water, and extensive fatigue damage may occur because of the larger curvature response of the umbilical. The umbilical configurations departing at 90 deg off incoming waves were found to have the highest fatigue life attributed to less extension or compression of the umbilical. However, additional bend stiffener/limiter features may need to be incorporated into the buoyancy section and touchdown regions to minimize curvature-induced fatigue.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2021

Pseudoelastic deformation in Mo-based refractory multi-principal element alloys

Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal-element alloys. We present phase and property analysis of quinary (MoW) x Zr y (TaTi) 1-x-y refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo-W-Ta-Ti-Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress-strain. Our analysis shows that only Mo-W-rich compositions of Mo-W-Ta-Ti-Zr, i.e., Mo+W, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW) 85 Zr 7.5 (TaTi) 7.5 was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning-assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. Furthermore, this alloy shows identical cyclic deformation characteristics during uniaxial < 100 > loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1,500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tailoring the Interfacial Interactions of Porous Polymer Membranes to Accelerate Atomic Layer Deposition: The Latent Path to Antifouling Membranes

Atomic layer deposition (ALD) is a powerful strategy to engineer hybrid organic-inorganic membranes with emergent functionalities. The combination of atomic-level thickness control, a wide materials palette, and unprecedented conformality allows the physiochemical properties (e.g., hydrophilicity) of mesoporous polymer membranes to be precisely tuned. The nucleation of ALD material growth on polymer surfaces relies on Lewis acid-base interactions and remains an overlooked motif with tremendous potential to accelerate ALD nucleation and growth. Strategies to enhance these interactions could enable desirable properties such as antifouling behavior to be imparted on inert polymer surfaces that lack the necessary functional groups for ALD nucleation. Here, in this study, we demonstrate that the reactivity of polyacrylonitrile (PAN) membranes toward ALD metal oxide (MO) precursors with Lewis acid characteristics is enhanced by introducing strong Lewis base functional groups (amidoxime: Am) on the PAN backbone (Am-PAN). The resulting Lewis acid-base interactions accelerate the MO nucleation in Am-PAN and reduce the number of deposition cycles required to achieve hydrophilicity compared with the untreated PAN membrane. Unveiling the reaction mechanism, in situ Fourier transform infrared (FTIR) spectroscopy measurements established enhanced interaction dynamics between the ALD MO precursors and the Am-PAN membrane, unlike the PAN membrane. Spectroscopic ellipsometry and thermogravimetric analysis measurements revealed higher MO loadings in Am-PAN membranes compared to PAN membranes for the same number of ALD cycles. We found that strong Lewis acid-base interactions accelerated the ALD for a range of materials including Al 2 O 3 , TiO 2 , SnO 2 , and ZnO. More broadly, our work demonstrates that tailoring metal-precursor-polymer interactions is a powerful strategy to accelerate and modulate the ALD. We used this design strategy to fabricate Al 2 O 3 -Am-PAN hybrid membranes that showed 2-fold higher antifouling capability compared to pristine PAN membranes prepared with an equivalent number of Al 2 O 3 ALD cycles. Our approach expands the scope of design options for fouling-resistant porous hybrid inorganic-organic membranes and may ultimately reduce the operational costs of water treatment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Servo-aero-gravo-elastic (SAGE) scaling and its application to a 13-MW downwind turbine

Reduced scale wind turbines can be extremely cost-effective to test new rotor concepts since prototype costs are heavily dependent on the rotor diameter. Ideally, the scaled model would have the same non-dimensional deflections, dynamics, and control behavior as the full-scale model. This would provide a high-fidelity demonstration of the full-scale performance, which is ideal if the full-scale turbine has significant aeroelastic interactions. To this end, servo-aero-gravo-elastic (SAGE) scaling is developed and applied to a 13-MW turbine that is scaled to a 20% scale model. The scaling preserves the tip-speed ratio, the rotor speed normalized by the flapping frequency, and the tip deflections normalized by the blade length. In addition, the controller employs the same control structure (gain-scheduled pitch control and variable speed torque control) and is scaled dynamically (e.g., matching non-dimensional time constant of the pitch angle, etc.). Furthermore, the thrust, gravity, and centrifugal moments are scaled such that the load angles are preserved as a function of a non-dimensional wind speed. However, the environmental scaling must consider differences in Reynolds number (since this parameter cannot be held constant) and subsequent changes in the axial induction factor. While the presented results showcase these differences during operational conditions, the non-dimensional tip deflections remain comparable through all wind speed ranges, indicating the viability of the SAGE scaling method in matching full-scale aeroelastic responses.

Kaminski, Meghan (ORCID:0000000322334784)↗
DOE OSTI2024

Final Report on Predictive Analyses of PRD as Function of Anomalies

This report summarizes the work completed in FY-2024 to analyze the power reactivity decrement (PRD) concepts of the ARC-100 core. The PRD has been traditionally defined by the reactivity change from a hot zero power (HZP) to a particular power. Consequently, the PRD accounts for the core reactivity changes due to increase coolant temperature gradient axially and radially across the core, and increased fuel temperature. The coolant temperature gradient leads to sodium and structure density changes, to radial core expansion from assembly flowering and bowing (due to axial and radial temperature gradients within the assemblies), and to control rod driveline thermal expansion. The fuel temperature increase associated with coolant temperature and power increases leads to Doppler effect and axial thermal expansion. In this work, the normal operating Hot Full Power (HFP) state is the only power level of interest, so analyses focus on the PRD calculated from HZP to HFP. The PRD has been used to assess the reactor safety features asymptotically in unprotected accident scenarios, including the loss of heat sink (LOHS), loss of flow (LOF), and transient overpower (TOP) without scram. The PRD concept relies on the “global” reactivity coefficients A, B, and C that are estimated based on “individual” reactivity effects (Doppler, sodium density, etc.). The objectives of this work are: 1) to improve and verify the methodology used to calculate the ABC coefficients used in the PRD, 2) to assess if the PRD can be used to reliably identify abnormal events. This report fulfills the FY-2024 scope of WBS#1.15.8.3 activity, “ANL0120 – Predictive Analyses of PRD as function of deformation”. The PRD concept is described in Sections 2. Additional effort in refining the methodology for core bowing modeling is performed in Sections 3. Then, two verification exercises are proposed in Section 4 to benchmark these coefficients based on direct neutronic-only calculations and on dynamic core transient simulations. Finally, the PRD approach is assessed for the detection of several unexpected events, such as primary flow perturbation or improper fuel loading, in Section 5.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Tidal Turbine Benchmarking Project: Stage I - Steady Flow Blind Predictions: Preprint

This paper presents the first blind prediction stage of the Tidal Turbine Benchmarking Project being conducted and funded by the UK's EPSRC and Supergen ORE Hub. In this first stage, only steady flow conditions, at low and elevated turbulence (3.1%) levels, were considered. Prior to the blind prediction stage, a large laboratory scale experiment was conducted in which a highly instrumented 1.6m diameter tidal rotor was towed through a large towing tank in well-defined flow conditions with and without an upstream turbulence grid. Details of the test campaign and rotor design were released as part of this community blind prediction exercise. Participants were invited to use a range of engineering modelling approaches to simulate the performance and loads of the turbine. 26 submissions were received from 12 groups from across academia and industry using solution techniques ranging from blade resolved computational fluid dynamics through actuator line, boundary integral element methods, vortex methods to engineering Blade Element Momentum methods. The comparisons between experiments and blind predictions were extremely positive helping to provide validation and uncertainty estimates for the models, but also validating the experimental tests themselves. The exercise demonstrated that the experimental turbine data provides a robust data set against which researchers and design engineers can test their models and implementations to ensure robustness in their processes, helping to reduce uncertainty and provide increased confidence in engineering processes. Furthermore, the data set provides the basis by which modellers can evaluate and refine approaches.

benchmarking↗
DOE OSTI2024

Retrofittable Thermal Switches for Dynamic Building Envelopes Integrated with Thermal Energy Storage: Preprint

Buildings in the United States consume about 40 quadrillion BTU of primary energy annually, which accounts for the nation's 40% of total energy use, 75% of all electricity use, and 35% of the net carbon emissions. Deploying thermal energy storage in the form of phase change material (PCM) in building envelopes is an effective method to reduce space heating/cooling loads, provide load shedding, and shift demand to periods of lower energy cost. However, the full potential of PCM-integrated envelopes can only be realized if the PCM undergoes complete phase change using free ambient heating/cooling, and the stored energy is effectively transferred between the exterior and the interior environments. Conventional thermal insulation (with a fixed thermal resistance) limits PCM utilization, particularly with the increasing emphasis on higher R-value in building envelopes, which negatively affects the energy-saving potential of a PCM-integrated envelope. In contrast, dynamic building envelopes integrated with PCMs provide the option of varying the thermal resistance based on the indoor and outdoor conditions, thereby enhancing utilization of free ambient cooling and heating to charge/discharge the PCM thermal storage, reducing the buildings' heating and cooling load, and shifting the peak energy demand. In this study, we demonstrate innovative retrofittable thermal switches in the form of the insertable plugs inside an insulation to provide variable thermal resistance depending on the operating temperature and direction of temperature gradient, thus allowing preferential directional heat flow. Notably, they are passive in nature, requiring no external power, and work solely based on the ambient temperature.

buildings↗
DOE OSTI2020

Crystallizing Atomic Xenon in a Flexible MOF to Probe and Understand Its Temperature-Dependent Breathing Behavior and Unusual Gas Adsorption Phenomenon

Flexible metal–organic frameworks (MOFs) hold great promise as smart materials for specific applications such as gas separation. These materials undergo interesting structural changes in response to guest molecules, which is often associated with unique adsorption behavior not possible for rigid MOFs. Understanding the dynamic behavior of flexible MOFs is crucial yet challenging as it involves weak host–guest interactions and subtle structural transformation not only at the atomic/molecular level but also in a nonsteady state. Here we report an in-depth study on the adsorbate- and temperature-dependent adsorption in a flexible MOF by crystallizing atomic gases into its pores. Mn(ina)2 shows an interesting temperature-dependent response toward noble gases. Its nonmonotonic, temperature-dependent adsorption profile results in an uptake maximum at a temperature threshold, a phenomenon that is unusual. Full characterization of Xe-loaded MOF structures is performed by in situ single-crystal and synchrotron X-ray diffraction, IR spectroscopy, and molecular modeling. The X-ray diffraction analysis offers a detailed explanation into the dynamic structural transformation and provides a convincing rationalization of the unique adsorption behavior at the molecular scale. The guest and temperature dependence of the structural breathing gives rise to an intriguing reverse of Xe/Kr adsorption selectivity as a function of temperature. The presented work may provide further understanding of the adsorption behavior of noble gases in flexible MOF structures.

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

Optimal Control Strategy With Efficiency and Reliability Improvement for Offshore DC Microgrids

Offshore microgrids, due to their remote location and lack of external energy support, face significant challenges in wide-range load operation and maintenance. Consequently, efficiency and reliability are critical concerns for converters in offshore dc microgrids. This article presents an optimal control strategy aimed at enhancing both efficiency and reliability. A normalized nonlinear relationship between power loss and thermal stress of a paralleled converter is first established. Based on this, a dual-objective optimization function with an active weight function as well as a system overall performance index is established. The active weight function dynamically adjusts the control priority based on converter efficiency and switching device thermal stress. Then, the optimal power-sharing strategy is derived by the Lagrange multiplier method with the proposed optimal function. Additionally, to accommodate a wide load range, an optimal selection strategy for operating converter combinations is proposed, requiring only low-bandwidth communication. Experiment verification is given to validate the effectiveness of the proposed control strategy. The experiment results demonstrate that the proposed control strategy can improve the overall performance of offshore microgrids by optimizing efficiency and reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

A Data-Driven Democratized Control Architecture for Regional Transmission Operators

As probably the most complicated and critical infrastructure system, U.S. power grids become increasingly vulnerable to extreme events such as cyber-attacks and severe weather, as well as higher DER penetrations and growing information mismatch among system operators, utilities (transmission or generation owners), and end-users. This paper proposes a data-driven democratized control architecture considering two democratization pathways to assist transmission system operators, with a targeted use case of developing online proactive islanding strategies. Detailed discussions on load capability profiling at transmission buses and disaggregation of DER generations are provided and illustrated with real-world utility data. By Combining network and operational constraints, transmission system operators can be equipped with new tools built on top of this architecture, to derive accurate, proactive, and strategic islanding decisions to incorporate the wide range of dynamic portfolios and needs when facing extreme events or unseen grid contingencies.

Power System Reliability, Islanding, Democratized ↗
DOE OSTI2024

Dynamic Behavior of Pt Multimetallic Alloys for Active and Stable Propane Dehydrogenation Catalysts

Improving the use of platinum in propane dehydrogenation catalysts is a crucial aspect to increasing the efficiency and sustainability of propylene production. A known and practiced strategy involves incorporating more abundant metals in supported platinum catalysts, increasing its activity and stability while decreasing the overall loading. Here, using colloidal techniques to control the size and composition of the active phase, we show that Pt/Cu alloy nanoparticles supported on alumina (Pt/Cu/Al 2 O 3 ) displayed elevated rates for propane dehydrogenation at low temperature compared to a monometallic Pt/Al 2 O 3 catalyst. We demonstrate that the enhanced catalytic activity is correlated with a higher surface Cu content and formation of a Pt-rich core and Cu-rich shell that isolates Pt sites and increases their intrinsic activity. However, rates declined on stream because of dynamic metal diffusion processes that led to a more uniform alloy structure. This transformation was only partially inhibited by adding excess hydrogen to the feed stream. Instead, cobalt was introduced to provide trimetallic Pt/Cu/Co catalysts with stabilized surface structure and stable activity and higher rates than the original Pt/Cu system. In conclusion, the structure–activity relationship insights in this work offer improved knowledge of propane dehydrogenation catalyst development featuring reduced Pt loadings and notable thermal stability for propylene production.

Alloys↗
DOE OSTI2023

Computational Fluid Dynamics Study of a Cross-Flow Marine Hydrokinetic Turbine and the Combined Influence of Struts and Helical Blades: Preprint

A computational fluid dynamics study was performed for a cross-flow marine hydro-kinetic turbine. The analysis was done in three dimensions and used the unsteady Reynolds averaged Navier-Stokes solver in the commercial code STAR-CCM+. The base turbine configuration is the RivGen® Turbine, designed by the Ocean Renewable Power Company (ORPC). A convergence and uncertainty analysis was performed for both the spatial and temporal discretization; this was done using the base configuration which features support struts and helical foils. The proposed study aims to compare the impact of the struts on both power performance and blade loading for helical and straight blades.

CFD↗
DOE OSTI2021

Informing the planning of rotating power outages in heat waves through data analytics of connected smart thermostats for residential buildings

Abstract With climate change, heat waves have become more frequent and intense. Rotating power outages happen when the power supply is unable to meet the cooling demand increase resulting from extreme high temperatures. Power outages during heat waves expose residents to high risks of overheating. In this study, we propose a novel data-driven inverse modelling approach to inform decision makers and grid operators on planning rotating power outages. We first infer the building thermal characteristics using the connected smart thermostat data, and used the estimated thermal dynamics to simulate the thermal resilience during a heat wave event. Our proposed method was tested for the California power outage in August 2020 by using the open source Ecobee Donate Your Data dataset. We found in California the power outage should not last more than two hours during heat waves to avoid overheating risks. Informing the residents in advance so they can prepare for it through pre-cooling is a simple but effective strategy to expand the acceptable power outage duration. In addition to assisting power outage planning, the proposed method can be used for other applications, such as to evaluate a building energy efficiency policy, to examine fuel poverty, and to estimate the load shifting potential of building stocks.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Dynamic Shear Deformation of a Precipitation Hardened Al 0.7 CoCrFeNi Eutectic High-Entropy Alloy Using Hat-Shaped Specimen Geometry

Lamellar eutectic structure in Al 0.7 CoCrFeNi high-entropy alloy (HEA) is emerging as a promising candidate for structural applications because of its high strength-ductility combination. The alloy consists of a fine-scale lamellar fcc + B2 microstructure with high flow stresses > 1300 MPa under quasi-static tensile deformation and >10% ductility. The response to shear loading was not investigated so far. This is the first report on the shear deformation of a eutectic structured HEA and effect of precipitation on shear deformation. A split-Hopkinson pressure bar (SHPB) was used to compress the hat-shaped specimens to study the local dynamic shear response of the alloy. The change in the width of shear bands with respect to precipitation and deformation rates was studied. The precipitation of L1 2 phase did not delay the formation of adiabatic shear bands (ASB) or affect the ASB width significantly, however, the deformed region around ASB, consisting of high density of twins in fcc phase, was reduced from 80 μm to 20 μm in the stronger precipitation strengthened condition. We observe dynamic recrystallization of grains within ASBs and local mechanical response of individual eutectic lamellae before and after shear deformation and within the shear bands was examined using nano-indentation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗