SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Temporal dynamics of large-scale structures for turbulent Rayleigh–Bénard convection in a moderate aspect-ratio cylinder

We investigate the spatial organization and temporal dynamics of large-scale, coherent structures in turbulent Rayleigh–Bénard convection via direct numerical simulation of a 6.3 aspect-ratio cylinder with Rayleigh and Prandtl numbers of 9.6×10 7 and 6.7, respectively. Fourier modal decomposition is performed to investigate the structural organization of the coherent turbulent motions by analysing the length scales, time scales and the underlying dynamical processes that are ultimately responsible for the large-scale structure formation and evolution. We observe a high level of rotational symmetry in the large-scale structure in this study and that the structure is well described by the first four azimuthal Fourier modes. Two different large-scale organizations are observed during the duration of the simulation and these patterns are dominated spatially and energetically by azimuthal Fourier modes with frequencies of 2 and 3. Studies of the transition between these two large-scale patterns, radial and vertical variations in the azimuthal energy spectra, as well as the spatial and modal variations in the system's correlation time are conducted. Rotational dynamics are observed for individual Fourier modes and the global structure with strong similarities to the dynamics that have been reported for unit aspect-ratio domains in prior works. It is shown that the large-scale structures have very long correlation time scales, on the order of hundreds to thousands of free-fall time units, and that they are the primary source for a horizontal inhomogeneity within the system that can be observed during a finite, but a very long-time simulation or experiment.

42 ENGINEERING↗

DOE OSTI2026

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗

DOE OSTI2024

Urbanization exacerbates continental- to regional-scale warming

Urbanization is usually ignored when estimating past changes in large-scale climate and for future climate projections since cities historically covered a small fraction of the Earth’s surface. Here, by combining global land surface temperature observations with historical estimates of urban area, we demonstrate that the urban contribution to continental- to regional-scale warming has become non-negligible, especially for rapidly urbanizing regions and countries in Asia. Consequently, expected urban expansion over the next century suggests further increased urban influence on large-scale surface climate in the future (approximately 0.16 K for North America and Europe for high-emission scenario in 2100). Based on these results, also seen for air temperature, we argue that, in line with other forms of land use/land cover change, urbanization should be explicitly included in climate change assessments. This requires incorporation of dynamic urban extent and biophysics in current-generation Earth system models to quantify potential urban feedback on the climate system across scales.

Chakraborty, TC↗

DOE OSTI2024

Exploring the Frontiers of Energy Efficiency using Power Management at System Scale

In the face of surging power demands for exascale HPC systems, this work tackles the critical challenge of understanding the impact of software-driven power management techniques like Dynamic Voltage and Frequency Scaling (DVFS) and Power Capping. These techniques have been actively developed over the past few decades. By combining insights from GPU benchmarking to understand application power profiles, we present a telemetry data-driven approach for deriving energy savings projections. This approach has been demonstrably applied to the Frontier supercomputer at scale. Our findings based on three months of telemetry data indicate that, for certain resource-constrained jobs, significant energy savings (up to 8.5%) can be achieved without compromising performance. This translates to a substantial cost reduction, equivalent to 1438 MWh of energy saved. The key contribution of this work lies in the methodology for establishing an upper limit for these best-case scenarios and its successful application. This work enables HPC professionals to optimize the power-performance trade-off within constrained power budgets, not only for the exascale era but also beyond.

Karimi, Ahmad Maroof↗

DOE OSTI2023

Field-scale dynamics of planting dates in the US Corn Belt from 2000 to 2020

Crop planting dates are a dynamic feature of agricultural systems that respond to short- and long-term climate signals, crop and cultivar selection, and technology changes. Planting date records are essential for yield gap analyses, accurate crop modeling, and tracking farmer adaptations to weather and climate change. Although planting dates have high variation at local scales due to heterogeneity in farm resources and decision-making, available long-term data on planting dates is largely restricted to aggregated regional statistics or, at best, satellite-derived datasets with limited spatiotemporal extent and at resolutions unable to distinguish individual fields (> 250 m). Here, we generated retrospective annual field-scale (30 m) planting date maps for both maize and soybeans spanning 2000-2020 across a 12 state region in the United States Corn Belt based on Landsat satellite data and a large ground sample of over 28,000 maize and soybean fields. Using training data from 2015-2020 for model selection, we found that planting date predictions improved with harmonic regression of Landsat data and additional annual weather covariates. The preferred random forests model approximately doubled performance compared to a null model based on state median planting dates, capturing 47% of field-level variation for maize (mean absolute error, MAE = 7.4 days) and 44% for soybeans (MAE = 7.5 days) against held-out ground truth test data for 2008-2014. We also evaluated the full 2000-2020 dataset with state agricultural statistics, finding strong agreement with median planting dates for maize (R 2 = 0.76, MAE = 4.4 days) and slightly lower agreement for soybeans (R 2 = 0.65, MAE = 5.4 days) when aggregated to the state level. We then used this new dataset to analyze environmental determinants of planting dates at a finer-scale than previously possible, controlling for unobserved variation at the sub-state district level. We found that during 2000-2020, each standard deviation increase in rainfall delayed planting by ~ 2.5 days, and fields with higher soil productivity ratings tended to be planted earlier. We did not find meaningful trends over the last two decades in planting dates for maize or soybeans, in contrast to trends towards earlier planting dates late last century and predicted for this period in climate adaptation studies. We hypothesize increases in early season rainfall may have inhibited these shifts towards earlier planting. Remotely sensed planting dates will be a useful tool for yield gap analyses, crop simulation modeling, and ongoing assessment of climate adaptation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Leveraging National Laboratory Assets to Address Stability Challenges due to Declining Grid Inertia Using Geographically Distributed Electrical–Thermal Co-Emulation

Due to increased penetration of low-inertia resources into the electric grid, challenges are increasing for maintaining wide-area system stability. Grid stability assessment requires a faithful representation of the multiple-physics interaction at the system level, and timescales of interaction varying in orders of magnitude, from microseconds to seconds to several minutes. Along with the simulation-based techniques, hardware-in-the-loop (HIL), controller HIL, and power HIL techniques have been developed to better understand the emergent behavior of the system with emerging technologies. US National Laboratories have played a vital role in research and development to understand the behavior of individual technologies and devices integrated to the electric grid. Each national laboratory forwards a technological and strategic initiative tied core and enabling capabilities. Due to strategic, efficiency, and economic reasons, not all the labs have assets to conduct research on all technologies concomitantly, so it becomes crucial to integrate the labs across geographies to understand the interplay of different technologies together at the system level. This approach avoids duplication of the assets at different lab facilities and helps understand the integrated system behavior of various technologies representative of actual grid conditions by connecting multiple national labs. This article talks about techniques of connecting three national laboratories to enable co-emulation of electrical–mechanical–thermal characteristics of devices and systems. Such an approach can be used to understand the dynamic and transient interaction of multi-physics in a system level, at-scale emulation using real-time simulation tools and techniques.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

High-Temperature Latent Heat Storage System using Transportable Heat Pipes for Versatile Integration with Emerging Microreactors

With the growing interest in the development of micro nuclear reactors and their application to microgrid, high-temperature thermal energy storage (TES) has also been receiving great attention as a technology to reinforce the energy utilization and profitability of this class of reactors. Given the importance and research needs of the TES technology, a research project is underway at Idaho National Laboratory (INL) to develop high-temperature TES that can be easily interfaced with emerging micro nuclear reactors. Currently under development is the Heat pipe-Integrated Thermal Battery (HITB), a novel latent heat storage design that enables flexible and robust integration with any design of currently considered micro nuclear reactors. The development of detailed design based on the proposed TES concept, HITB, requires not only development and validation of modeling and simulation capabilities, but also research that spans thermal performance analysis, component- and system-scale design optimization, and experimental demonstrations including high-temperature feasibility testing. The current research efforts (and plans) are divided into the following four areas: (i) Thermal and structural modeling and analysis at both system- and component-scales for evaluating performance and conducting design optimization, (ii) Small-scale experimental demonstrations for the key components of the proposed HITB design, (iii) Modelica-based dynamic model development for developing a control strategy under dynamic operation within integrated energy systems configuration, (iv) Integrated experimental testing and demonstration. This paper introduces the overall research efforts that have been and are being conducted to develop and demonstrate the novel design of the high-temperature latent heat storage device, HITB. The key design features, ongoing research efforts including the development of modeling and simulation capabilities, small-scale demonstrations of the HITB’s key functionalities prior to scaled-up experimental demonstration, are discussed. Then, this paper concludes with a remark on the future research plans for this project.

25 ENERGY STORAGE↗

DOE OSTI2022

Statics and dynamics of skyrmions interacting with disorder and nanostructures

Magnetic skyrmions are topologically stable nanoscale particlelike objects that were discovered in 2009. Since that time, intense research interest in the field has led to the identification of numerous compounds that support skyrmions over a range of conditions spanning from cryogenic to room temperatures. Skyrmions can be set into motion under various types of driving, and the combination of their size, stability, and dynamics makes them ideal candidates for numerous applications. At the same time, skyrmions represent a new class of system in which the energy scales of the skyrmion-skyrmion interactions, sample disorder, temperature, and drive can compete. A growing body of work indicates that the static and dynamic states of skyrmions can be influenced strongly by pinning or disorder in the sample; thus, an understanding of such effects is essential for the eventual use of skyrmions in applications. The current state of knowledge regarding individual skyrmions and skyrmion assemblies interacting with quenched disorder or pinning is reviewed. Additionally, the microscopic mechanisms for skyrmion pinning, including the repulsive and attractive interactions that can arise from impurities, grain boundaries, or nanostructures, are outlined. This is followed by descriptions of depinning phenomena, sliding states over disorder, the effect of pinning on the skyrmion Hall angle, the competition between thermal and pinning effects, the control of skyrmion motion using ordered potential landscapes such as one- or two-dimensional periodic asymmetric substrates, the creation of skyrmion diodes, and skyrmion ratchet effects. Highlighted are the distinctions arising from internal modes and the strong gyrotropic or Magnus forces that cause the dynamical states of skyrmions to differ from those of other systems with pinning, such as vortices in type-II superconductors, charge density waves, or colloidal particles. Throughout this review future directions and open questions related to the pinning and dynamics in skyrmion systems are also discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2024

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗

DOE OSTI2021

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗

DOE OSTI2021

Improved Time-Stepping Methods in Global to Regional Ocean Modeling (Annual Status Report 2020)

Time stepping algorithms are an important part of ocean models, and strongly influence both the accuracy of solution and performance. There have been a number of projects investigating various improvements for ocean time-stepping schemes in the Model for Prediction Across Scales-Ocean (MPAS-Ocean), a component of the DOE Energy Exascale Earth System Model. Ocean dynamics include fast surface gravity waves, which are two-dimensional, and slower internal waves, which are three-dimensional, so ocean models use a split time-stepping scheme that separates these barotropic and baroclinic modes for efficiency. MPAS-Ocean runs on variable-resolution horizontal meshes, and must scale to tens of thousands of cores and millions of horizontal gridcells. Ocean models require time stepping algorithms that are customized to these needs, and which are tuned for performance on various resolutions and architectures.

58 GEOSCIENCES↗

DOE OSTI2022

Improved Time-Stepping Methods in Global to Regional Ocean Modeling (Annual Status Report)

58 GEOSCIENCES↗

DOE OSTI2021

A probabilistic graphical model foundation for enabling predictive digital twins at scale

A unifying mathematical formulation is needed to move from one-off digital twins built through custom implementations to robust digital twin implementations at scale. This work proposes a probabilistic graphical model as a formal mathematical representation of a digital twin and its associated physical asset. We create an abstraction of the asset–twin system as a set of coupled dynamical systems, evolving over time through their respective state spaces and interacting via observed data and control inputs. The formal definition of this coupled system as a probabilistic graphical model enables us to draw upon well-established theory and methods from Bayesian statistics, dynamical systems and control theory. The declarative and general nature of the proposed digital twin model make it rigorous yet flexible, enabling its application at scale in a diverse range of application areas. Here, we demonstrate how the model is instantiated to enable a structural digital twin of an unmanned aerial vehicle (UAV). The digital twin is calibrated using experimental data from a physical UAV asset. Its use in dynamic decision-making is then illustrated in a synthetic example where the UAV undergoes an in-flight damage event and the digital twin is dynamically updated using sensor data. The graphical model foundation ensures that the digital twin calibration and updating process is principled, unified and able to scale to an entire fleet of digital twins.

42 ENGINEERING↗

DOE OSTI2026

Hierarchical Testing of a Hybrid Machine Learning‐Physics Global Atmosphere Model

Machine learning (ML)-based models have demonstrated high skill and computational efficiency, often outperforming conventional physics-based models in weather and subseasonal predictions. While prior studies have assessed their fidelity in capturing synoptic-scale atmospheric dynamics, their performance across timescales and under out-of-distribution forcing, such as +3K or +4K uniform-warming forcings, and the sources of biases remain elusive, to establish the model's reliability for Earth science. Here, we design three sets of experiments targeting synoptic-scale phenomena, interannual variability, and out-of-distribution uniform-warming forcings. We evaluate the Neural General Circulation Model (NeuralGCM), a hybrid model integrating a dynamical core with ML-based component, against observations and physics-based Earth system models (ESMs). At the synoptic scale, NeuralGCM captures the evolution and propagation of extratropical cyclones with performance comparable to ESMs. At the interannual scale, when forced by El Niño-Southern Oscillation sea surface temperature (SST) anomalies, NeuralGCM successfully reproduces associated teleconnection patterns but exhibits deficiencies in capturing nonlinear response. Under out-of-distribution uniform-warming forcings, NeuralGCM simulates similar responses in global-average temperature and precipitation and reproduces large-scale tropospheric circulation features similar to those in ESMs. Notable weaknesses include overestimating the tracks and spatial extent of extratropical cyclones, biases in the teleconnected wave train triggered by tropical SST anomalies, and differences in upper-level warming and stratospheric circulation responses to SST warming compared to physics-based ESMs. The causes of these weaknesses were explored. Despite the noted weaknesses, NeuralGCM reproduces responses across experiments reasonably and performs comparably to ESMs. By integrating a dynamical core with ML, NeuralGCM shows potential for developing ML-based ESMs.

global warming↗

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↗