Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “residual network”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Structural basis of transcription: RNA polymerase II substrate binding and metal coordination using a free-electron laser

Catalysis and translocation of multisubunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near-atomic resolution and precise arrangement of key active site components have been elusive. Here, we present the free-electron laser (FEL) structures of a matched ATP-bound Pol II and the hyperactive Rpb1 T834P bridge helix (BH) mutant at the highest resolution to date. The radiation-damage-free FEL structures reveal the full active site interaction network, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg 2+ , and, more importantly, a putative third (site C) Mg 2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structures indicate that the third Mg 2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/BH interactions induce conformational changes that could allow translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the T834P mutant reveal rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Science & Technology - Other Topics↗
DOE OSTI2024

Independent Mutation of Two Bridging Carboxylate Ligands Stabilizes Alternate Conformers of the Photosynthetic O 2 -Evolving Mn 4 CaO 5 Cluster in Photosystem II

The O 2 -evolving Mn 4 CaO 5 cluster in Photosystem II is ligated by six carboxylate residues. One of these is D170 of the D1 subunit. This carboxylate bridges between one Mn ion (Mn4) and the Ca ion. A second carboxylate ligand is D342 of the D1 subunit. This carboxylate bridges between two Mn ions (Mn1 and Mn2). D170 and D342 are located on opposite sides of the Mn 4 CaO 5 cluster. Recently, it was shown that the D170E mutation perturbs both the intricate networks of H-bonds that surround the Mn 4 CaO 5 cluster and the equilibrium between different conformers of the cluster in two of its lower oxidation states, S 1 and S 2 , while still supporting O 2 evolution at approx. 50% the rate of wild-type. In this study, we show that the D342E mutation produces much the same alterations to the cluster’s FTIR and EPR spectra as D170E, while still supporting O 2 evolution at approx. 20% the rate of wild-type. Furthermore, the double mutation, D170E + D342E, behaves similarly to the two single mutations. We conclude that D342E alters the equilibrium between different conformers of the cluster in its S 1 and S 2 states in the same manner as D170E and perturbs the H-bond networks in a similar fashion. This is the second identification of a Mn 4 CaO 5 metal ligand whose mutation influences the equilibrium between the different conformers of the S 1 and S 2 states without eliminating O 2 evolution. Finally, this finding has implications for our understanding of the mechanism of O 2 formation in terms of catalytically active/inactive conformations of the Mn 4 CaO 5 cluster in its lower oxidation states.

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

Three-dimensional reconstruction of inertial confinement fusion hot-spot plasma from x-ray and nuclear diagnostics on OMEGA

Multidimensional effects degrade the neutron yield and the compressed areal density of laser-direct-drive inertial confinement fusion implosions of layered deuterium–tritium cryogenic targets on the OMEGA Laser System with respect to 1D radiation-hydrodynamic simulation predictions. A comprehensive physics-informed 3D reconstruction effort is under way to infer hot-spot and shell conditions at stagnation from four x-ray and seven neutron detectors distributed around the OMEGA target chamber. Neutron diagnostics, providing measurements of the neutron yield, hot-spot flow velocity, and apparent ion-temperature distribution, are used to infer the mode-1 perturbation at stagnation. The x-ray imagers record the shape of the hot-spot plasma to diagnose mode-1 and mode-2 perturbations. A deep-learning convolutional neural network trained on an extensive set of 3D radiation-hydrodynamic simulations is used to interpret the x-ray and nuclear measurements to infer the 3D profiles of the hot-spot plasma conditions and the amount of laser energy coupled to the hot-spot plasma. A 3D simulation database shows that larger mode-1 asymmetries are correlated with higher hot-spot flow velocities and reduced laser-energy coupling and neutron yield. Three-dimensional hot-spot reconstructions from x-ray measurements indicate that higher amounts of residual kinetic energy are correlated with higher measured hot-spot flow velocities, consistent with 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Amino acid interacting network in the receptor-binding domain of SARS-CoV-2 spike protein

The relation between amino acid (AA) sequence and biologically active conformation controls the process of polypeptide chains folding into three-dimensional (3d) protein structures. The recent achievements in the resolution achieved in cryo-electron microscopy coupled with improvements in computational methodologies have accelerated the analysis of structures and properties of proteins. However, the detailed interaction between AAs has not been fully elucidated. Herein, we present a de novo method to evaluate inter-amino acid interactions based on the concept of accurately evaluating the amino acid bond pairs (AABP). The results obtained enabled the identification of complex 3d long-range interconnected AA interacting network in proteins. The method is applied to the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. We show that although nearest-neighbor AAs in the primary sequence have large AABP, other nonlocal AAs make substantial contribution to AABP with significant participation of both covalent and hydrogen bonding. Detailed analysis of AABP in RBD reveals the pivotal role they play in sequence conservation with profound implications on residue mutations and for therapeutic drug design. This approach could be easily applied to many other proteins of biomedical interest in life sciences.

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

Reconstructing cosmological initial conditions from late-time structure with convolutional neural networks

ABSTRACT We present a method to reconstruct the initial linear-regime matter density field from the late-time non-linearly evolved density field in which we channel the output of standard first-order reconstruction to a convolutional neural network (CNN). Our method shows dramatic improvement over the reconstruction of either component alone. We show why CNNs are not well-suited for reconstructing the initial density directly from the late-time density: CNNs are local models, but the relationship between initial and late-time density is not local. Our method leverages standard reconstruction as a preprocessing step, which inverts bulk gravitational flows sourced over very large scales, transforming the residual reconstruction problem from long-range to local and making it ideally suited for a CNN. We develop additional techniques to account for redshift distortions, which warp the density fields measured by galaxy surveys. Our method improves the range of scales of high-fidelity reconstruction by a factor of 2 in wavenumber above standard reconstruction, corresponding to a factor of 8 increase in the number of well-reconstructed modes. In addition, our method almost completely eliminates the anisotropy caused by redshift distortions. As galaxy surveys continue to map the Universe in increasingly greater detail, our results demonstrate the opportunity offered by CNNs to untangle the non-linear clustering at intermediate scales more accurately than ever before.

Astronomy & Astrophysics↗
DOE OSTI2026

Soil physical and chemical measurements for topsoils collected during NEON campaign in East River, CO (06/14/2018-06/28/2018)

The package is part of the DOE Watershed Function Science Focus Area (SFA) project and includes soil physical and chemical measurements from topsoils collected at the East River, Colorado, in conjunction with the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) survey conducted in June 2018. The soil measurements include soil bulk density, soil volumetric water content, soil microbial biomass C (Carbon), N (Nitrogen) and C:N (C to N ratio), soil DNA yield, soil total extractable organic C, soil total extractable N, soil extractable nitrate, soil extractable ammonium, soil dissolved inorganic N, soil dissolved organic N, soil pH, soil TOC400 (total organic carbon at 400°C), soil ROC (residual oxidizable carbon), soil TIC (total inorganic carbon), soil TOC (total organic carbon), soil TC (total carbon), soil N, soil OM (organic matter) loss on ignition. Additional associated site metadata can be found in the ESS-DIVE package 10.15485/1618130. The dataset includes (1) 2018_NEON_soil_physical_chemical_measurements.csv: soil physical and chemical measurements indexed by soil sample IGSNs; (2) samples.csv: sample metadata file used to register International Generic Sample Numbers (IGSNs); (3) flmd.csv: file level metadata file; and (4) dd.csv: data dictionary file. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS (Catchment Hydrology and ↗
DOE OSTI2024

Study of electron-induced chemical transformations in polymers

In extreme ultraviolet (EUV) photoresist exposure, the primary and secondary electrons drive chemistry rather than the EUV photons themselves. These electrons have a wide range of energies below approximately 80 eV, which are capable of complex network of reactions during exposure. To better understand the ability of electrons of different energies within the EUV primary and secondary electron range, we want to characterize and compare the chemistry induced in pure polymer films by direct exposure to electrons. Thin films of poly(tert-butyl methacrylate), poly(methyl methacrylate), and poly(4-hydroxystyrene) were exposed to a 20 to 80 eV electron beam. Outgassing during exposure was characterized in-situ using a quadrupole residual gas analyzer. The thickness changes were measured using ellipsometry and chemical bond structure data were collected using Fourier-transform infrared spectroscopy (FTIR) after exposure to compare different exposure conditions. Poly(4-hydroxystyrene) demonstrated stability during electron exposures. Exposures of the other two materials led to outgassing of protecting groups, the intensity of which decayed in time. Outgassing, FTIR, and thickness loss data exhibited approximately linear relationships to each other.

Mueller, Maximillian↗
DOE OSTI2022

GUTs, hybrid topological defects, and gravitational waves

The symmetry breaking of grand unified gauge groups in the early Universe often leaves behind relic topological defects such as cosmic strings, domain walls, or monopoles. For some symmetry breaking chains, hybrid defects can form where cosmic strings attach to domain walls or monopoles attach to strings. In general, such hybrid defects are unstable, with one defect "eating" the other via the conversion of its rest mass into the other's kinetic energy and, subsequently, decaying via gravitational waves. In this work, we determine the gravitational wave spectrum from 1) the destruction of a cosmic string network by the nucleation of monopoles which cut up and "eat" the strings, 2) the collapse and decay of a monopole-string network by strings that eat the monopoles, 3) the destruction of a domain wall network by the nucleation of string-bounded holes on the wall that expand and eat the wall, and 4) the collapse and decay of a string-bounded wall network by walls that eat the strings. We call the gravitational wave signals produced from the eating of one topological defect by another "gravitational wave gastronomy." We find that the four gravitational wave gastronomy signals considered yield unique spectra that can be used to narrow down the SO (10) symmetry breaking chain to the Standard Model and the scales of symmetry breaking associated with the consumed topological defects. Moreover, the systems we consider are unlikely to have a residual monopole or domain wall problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Record-Low Thermal Boundary Resistance between Diamond and GaN-on-SiC for Enabling Radiofrequency Device Cooling

The implementation of 5G-and-beyond networks requires faster, high-performance, and power-efficient semiconductor devices, which are only possible with materials that can support higher frequencies. Gallium nitride (GaN) power amplifiers are essential for 5G-and-beyond technologies since they provide the desired combination of high frequency and high power. These applications along with terrestrial hub and backhaul communications at high power output can present severe heat removal challenges. The cooling of GaN devices with diamond as the heat spreader has gained significant momentum since device self-heating limits GaN’s performance. However, one of the significant challenges in integrating polycrystalline diamond on GaN devices is maintaining the device performance while achieving a low diamond/GaN channel thermal boundary resistance. In this study, we achieved a record-low thermal boundary resistance of around 3.1 ± 0.7 m 2 K/GW at the diamond/Si 3 N 4 /GaN interface, which is the closest to theoretical prediction to date. The diamond was integrated within ~1 nm of the GaN channel layer without degrading the channel’s electrical behavior. Furthermore, we successfully minimized the residual stress in the diamond layer, enabling more isotropic polycrystalline diamond growth on GaN with thicknesses >2 μm and a ~1.9 μm lateral grain size. More isotropic grains can spread the heat in both vertical and lateral directions efficiently. Using transient thermoreflectance, the thermal conductivity of the grains was measured to be 638 ± 48 W/m K, which when combined with the record-low thermal boundary resistance makes it a leading-edge achievement.

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

Deep inference of simulated strong lenses in ground-based surveys

The large number of strong lenses discoverable in future astronomical surveys will likely enhance the value of strong gravitational lensing as a cosmic probe of dark energy and dark matter. However, leveraging the increased statistical power of such large samples will require further development of automated lens modeling techniques. We show that deep learning and simulation-based inference (SBI) methods produce informative and reliable estimates of parameter posteriors for strong lensing systems in ground-based surveys. We present the examination and comparison of two approaches to lens parameter estimation for strong galaxy-galaxy lenses — Neural Posterior Estimation (NPE) and Bayesian Neural Networks (BNNs). We perform inference on 1-, 5-, and 12-parameter lens models for ground-based imaging data that mimics the Dark Energy Survey (DES). We find that NPE outperforms BNNs, producing posterior distributions that are more accurate, precise, and well-calibrated for most parameters. For the 12-parameter NPE model, the calibration is consistently within <10% of optimal calibration for all parameters, while the BNN is rarely within 20% of optimal calibration for any of the parameters. Similarly, residuals for most of the parameters are smaller (by up to an order of magnitude) with the NPE model than the BNN model. This work takes important steps in the systematic comparison of methods for different levels of model complexity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Phonon-informed Neural Thermal Scattering (NeTS) Optimization for Crystalline Graphite and Beryllium Metal

Fast neutrons born from fission lose energy through scattering interactions in the process of slowing-down. As neutrons thermalize to the order of $k$ $b$ $T$ (where $k$ $b$ is the Boltzmann constant, and $T$ is the temperature of the medium), their de Broglie wavelength and energy approaches the order of inter-atomic spacing and quantized lattice vibrations, i.e., phonons. At thermal energies, the thermal scattering law (TSL), i.e., $S$($α, β$), captures crystal binding contributions to the total reaction rate, or cross section. This dimensionless material property describes the energy ($β$) and momentum ($α$) exchanges available in a medium. Currently, $S$($α, β$) is evaluated in the Full Law Analysis Scattering System Hub (FLASSH) code for discrete inputs and stored as ENDF/B File 7 for 0-phonon elastic (MT 2) and n-phonon inelastic (MT 4) processes. Further processing recasts $S$($α, β$) into cumulative distribution functions for sampling post-collision scattering kinematics. In practice, interpolation schemes are employed to access data between tabulated values. An improvement to this juncture of the nuclear data pipeline is supplying cross sections on-the-fly (OTF), as has been developed for the un-resolved resonance region to minimize non-physical interpolation errors. This capability may improve simulation accuracy for accident and transient analyses, where rapidly varying changes in temperature and pressure are difficult to predict beforehand. To do so, deep artificial neural networks (ANNs) can be employed which collapse non-linear, complex data into a lightweight dictionary of neural weights and biases. This has been successfully demonstrated for the hydrogen in light water $S$($α, β$) dataset in the form of a Neural Thermal Scattering (NeTS) module. In this work, the NeTS framework is extended to consider the impact of material-dependent dynamical features on optimal neural pre-processing and architecture design decisions, such as number of neurons per hidden layer, residual skip connections and neural depth. New NeTS modules for crystalline graphite and beryllium metal illuminate a novel correlation between dynamical nonlinearity and optimal neural parametrization when deploying $S$($α, β$) on-the-fly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Data-driven recovery of hidden physics in reduced order modeling of fluid flows

In this article, we introduce a modular hybrid analysis and modeling (HAM) approach to account for hidden physics in reduced order modeling (ROM) of parameterized systems relevant to fluid dynamics. The hybrid ROM framework is based on using first principles to model the known physics in conjunction with utilizing the data-driven machine learning tools to model the remaining residual that is hidden in data. This framework employs proper orthogonal decomposition as a compression tool to construct orthonormal bases and Galerkin projection (GP) as a model to build the dynamical core of the system. Our proposed methodology hence compensates structural or epistemic uncertainties in models and utilizes the observed data snapshots to compute true modal coefficients spanned by these bases. The GP model is then corrected at every time step with a data-driven rectification using a long short-term memory (LSTM) neural network architecture to incorporate hidden physics. A Grassmannian manifold approach is also adopted for interpolating basis functions to unseen parametric conditions. The control parameter governing the system's behavior is thus implicitly considered through true modal coefficients as input features to the LSTM network. The effectiveness of the HAM approach is then discussed through illustrative examples that are generated synthetically to take hidden physics into account. Furthermore, our approach thus provides insights addressing a fundamental limitation of the physics-based models when the governing equations are incomplete to represent underlying physical processes.

42 ENGINEERING↗
DOE OSTI2025

Energy metric prediction for double insertion mutants via the RoseNet deep learning framework

Studying the structural and functional implications of protein mutations is an important task in computational biology and bioinformatics. We leverage our previously proposed RoseNet neural network architecture to predict energy metrics of proteins with double amino acid insertions or deletions (InDels). We train models on previously generated benchmark datasets containing the exhaustive double InDel mutations for three proteins, as well as an additional three proteins for which ∼145k random mutants, each with two InDels, have been generated. We expand on our previous work by evaluating three additional proteins and analyzing domain features that impact the prediction capabilities of RoseNet. These features include InDels into secondary structures and the solvent accessible surface area (SASA) scores of the residues. We uncover further evidence to support that RoseNet has a higher proficiency of generalizing to unseen residue combinations than unseen insertion positions. We also observe that RoseNet produces higher-quality predictions when inserting into a β-sheet over an α-helix. Additionally, when the insertions fall in an area of high SASA, RoseNet often displays better performance than inserting into areas of low SASA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗
DOE OSTI2024

De novo design of pH-responsive self-assembling helical protein filaments

Abstract Biological evolution has led to precise and dynamic nanostructures that reconfigure in response to pH and other environmental conditions. However, designing micrometre-scale protein nanostructures that are environmentally responsive remains a challenge. Here we describe the de novo design of pH-responsive protein filaments built from subunits containing six or nine buried histidine residues that assemble into micrometre-scale, well-ordered fibres at neutral pH. The cryogenic electron microscopy structure of an optimized design is nearly identical to the computational design model for both the subunit internal geometry and the subunit packing into the fibre. Electron, fluorescent and atomic force microscopy characterization reveal a sharp and reversible transition from assembled to disassembled fibres over 0.3 pH units, and rapid fibre disassembly in less than 1 s following a drop in pH. The midpoint of the transition can be tuned by modulating buried histidine-containing hydrogen bond networks. Computational protein design thus provides a route to creating unbound nanomaterials that rapidly respond to small pH changes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Phosphoserine Charge State Drives Ion Condensation and Spatial Polyamine Presentation in Multirepeat Silaffin

Diatom silaffins direct silica biomineralization through heavily post-translationally modified repeat domains, yet how these modifications reshape the multirepeat conformational ensemble remains unknown. We report all-atom MD simulations of a 195-residue construct spanning repeats R1–R7 of Sil1p from Cylindrotheca fusiformis, carrying the full complement of native PTMs: phosphoserine (pSer), long-chain polyamines (LCPA), dimethyllysine (MLY), and trimethylhydroxylysine phosphate (TPL). We simulated three variants (Native, P1/singly deprotonated phosphate, and P2/doubly deprotonated phosphate) at two NaCl concentrations in triplicate for 500 ns each. All systems disorder from the AlphaFold 3 starting structure. Phosphate charge state, not ionic strength, is the dominant control of ensemble compaction and ion organization. Doubly deprotonated phosphate organizes an extensive Na + condensation shell (∼100 ions, 20% of box Na + ) and a heterogeneous bridging network that integrates both pSer and TPL phosphate groups. The resulting ensemble is compact with LCPA side chains exhibiting above-median solvent accessibility in 84% of simulation frames in P2 at 300 mM. This is higher than any other condition we simulated and indicates that polyamine groups are preferentially surface-presented in the most compact, ion-organized state. A charge-neutralization control confirms that this compact state is a structured intermediate maintained by the bridging network, not a simple collapsed globule. Here, this repeat-scale spatial organization is not captured by single-repeat peptide studies. Understanding the dynamics, mechanism, and spatial organization of PTM-rich silaffin at the repeat scale is a step closer to hierarchical biomimetic materials beyond simple silica morphologies.

Amines↗
DOE OSTI2021

Structure of the polysaccharide sheath from the B race of the green microalga Botryococcus braunii

The green microalga Botryococcus braunii forms a colony of cells held together by an intricate extracellular matrix. This extracellular matrix also acts as a storage location for liquid hydrocarbons initially produced inside the cells. These hydrocarbons can be used as feedstocks to produce transportation fuels. In the B race of B. braunii, the extracellular matrix is composed of several elements: a cross-linked hydrocarbon network filled with the liquid hydrocarbons, a retaining wall to hold the hydrocarbons in the extracellular matrix, an arabinogalactan polysaccharide fibrillar sheath that extends outward from the retaining wall, and a single protein known as polysaccharide associated protein is found at the base of the polysaccharide fibers where they connect to the retaining wall. During growth, colonies of the B race of B. braunii shed fragments of the retaining wall with polysaccharide fibrils attached into the media. Here we describe the isolation and structural analysis of these polysaccharides using a series of chemical degradations in combination with glycosyl residue composition and linkage analysis, as well as 2-dimensional nuclear magnetic resonance spectroscopy. The results demonstrate that the polysaccharide is highly branched, consisting of a galactan backbone with short branches comprised of Dgalactopyranose, L-arabinofuranose, 4-O-methyl-D-glucopyranuronic acid, and 6-deoxy-D-altropyranose. To our knowledge, this is the first description of 6-deoxyaltrose in algae or plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Scalable and Distributed Algorithm for Managing Residential Demand Response Programs using Alternating Direction Method of Multipliers (ADMM)

For effective engagement of residential demand-side resources and to ensure efficient operation of distribution networks, we must overcome the challenges of controlling and coordinating residential components and devices at scale. To overcome this challenge, we present a distributed and scalable algorithm with a three-level hierarchical information exchange architecture for managing the residential demand response programs. First, a centralized optimization model is formulated to maximize community social welfare. Then, this centralized model is solved in a distributed manner with alternating direction method of multipliers (ADMM) by decomposing the original problem to utility-level and house-level problems. The information exchange between the different layers is limited to the primary residual (i.e., supply-demand mismatch), Lagrangian multipliers, and the total load of each house to protect each customer’s privacy. Simulation studies are performed on the IEEE 33 bus test system with 605 residential customers. The results demonstrate that the proposed approach can save customers’ electricity bills and reduce the peak load at the utility level without much affecting customers’ comfort and privacy. Finally, a quantitative comparison of the distributed and centralized algorithms shows the scalability advantage of the proposed ADMM-based approach, and it gives benchmarking results with achievable value for future research works.

24 POWER TRANSMISSION AND DISTRIBUTION↗