Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power optimizer”

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 559 records · Page 31

Joint scheduling of energy, fast and primary frequency response reserves in integrated transmission–distribution networks

Inverter-based distributed energy resources (DERs) connected to distribution networks (DNs) can provide fast frequency support, but their reserve deliverability depends on feeder constraints and differs from synchronous primary frequency response (PFR). Existing transmission–distribution coordination studies usually treat reserve generically or neglect feeder-level feasibility, while frequency-security scheduling studies rarely represent distribution feeders explicitly. This paper develops a bi-level day-ahead scheduling framework for integrated transmission–distribution networks that jointly clears energy, transmission-side PFR, and distribution-side fast frequency response (FFR) under exogenous hourly inertia and largest-loss inputs from an external unit commitment (UC) schedule. The transmission problem is modeled with DC-optimal power flow (OPF) and closed-form second-order cone (SOC) frequency-security constraints, whereas each DN is represented by a reserve-aware branch-flow AC-OPF so that scheduled fast reserves remain deliverable during activation. The bi-level problem is reformulated through Karush–Kuhn–Tucker (KKT) conditions into a mixed-integer SOC program, and a penalty term is used to tighten the distribution-network relaxation. In the reduced test system, lower exogenous inertia increased the required primary response from 179.64 MW to 191.08 MW, distribution-side fast response reduced total frequency-response procurement by up to 4.9%, and neglecting distribution constraints overstated the combined distribution-side energy and reserve award by up to 18%. In the expanded study, the largest case was solved in 2.02 s with a 0.00% optimality gap. Time-domain simulations kept the frequency nadir above 59.0 Hz in all tested hours. These results demonstrate the value of fast-response modeling and distribution-feasible reserve delivery in coordinated market clearing.

Noh, Seung-Gil↗

Power Prediction of Airborne Wind Energy Systems Using Multivariate Machine Learning

Kites can be used to harvest wind energy at higher altitudes while using only a fraction of the material required for conventional wind turbines. In this work, we present the kite system of Kyushu University and demonstrate how experimental data can be used to train machine learning regression models. The system is designed for 7 kW traction power and comprises an inflatable wing with suspended kite control unit that is either tethered to a fixed ground anchor or to a towing vehicle to produce a controlled relative flow environment. A measurement unit was attached to the kite for data acquisition. To predict the generated tether force, we collected input–output samples from a set of well-designed experimental runs to act as our labeled training data in a supervised machine learning setting. We then identified a set of key input parameters which were found to be consistent with our sensitivity analysis using Pearson input–output correlation metrics. Finally, we designed and tested the accuracy of a neural network, among other multivariate regression models. The quality metrics of our models show great promise in accurately predicting the tether force for new input/feature combinations and potentially guide new designs for optimal power generation.

17 WIND ENERGY↗

Area-preserving scheme for efficiency enhancement in single-pass tapered free electron lasers

Efficiency enhancement of a single-pass short-wavelength high-gain tapered free electron laser (FEL) has recently been intensively studied. The goal is to sustain the growth of radiation power in the post-saturation regime. Among the various schemes, the undulator tapering is considered an effective route to achieve higher power extraction efficiency. The tapering strategy can be of constant or varying resonant phase along the undulator axis. In this paper we propose an efficiency-enhancement scheme based on preservation of the longitudinal phase space area which ensures trapping of resonant particles in the ponderomotive bucket as long as possible along the undulator axis before significant particle depletion occurs. In the meanwhile such a scheme takes advantage of the increase of the radiation field amplitude to precipitate the particle deceleration process at the middle stage of undulator tapering. We analyze such an area-preserving scheme of undulator tapering by formulating the post-saturation FEL interaction in a one-dimensional (1-D) model via introduction of the particle trapping fraction. The output performance is evaluated through numerical iteration and confirmed with 1-D particle tracking simulations. The results show that the optimal power extraction efficiency based on the proposed scheme, together with a prebunched beam, can be greatly improved within relatively short taper length compared with other schemes before radiation diffraction effect becomes significant. Besides, the undesired sideband effects are found to be effectively suppressed. For the proposed area-preserving taper scheme, we also derive an analytical approximate formula for the resonant phase as a function of undulator axis. In conclusion, we expect that the analysis can shed light on the aim to further enhance the power extraction efficiency in single-pass tapered FELs.

47 OTHER INSTRUMENTATION↗

Ensuring Transient Stability With Guaranteed Region of Attraction in DC Microgrids

DC microgrids have promising applications in renewable integration due to their better energy efficiency when connecting DC components. However, they might be unstable since many loads in a DC microgrid are regulated as constant power loads (CPLs) that have a destabilizing negative impedance effect. As a result, the state trajectory displacement caused by abrupt load changes or contingencies can easily lead to instability. Many existing works have been devoted to studying the region of attraction (ROA) of a DC microgrid, in which the system is guaranteed to be asymptotically stable. Nevertheless, existing work either focuses on using numerical methods for ROA approximations that generally have no performance guarantees or cannot ensure a desired ROA for a general DC microgrid. To close this gap, this paper develops an innovative control synthesis algorithm to make a general DC microgrid have a theoretically guaranteed ROA, for example, to cover the entirety of its operating range regarding state trajectories. Here, we first study the nonlinear dynamics of a DC microgrid to derive a novel transient stability condition to rigorously certify whether a given operating range is a subset of the ROA; then, we formulate a control synthesis optimization problem to guarantee the condition’s satisfaction. This condition is a linear constraint, and the optimization problem resembles an optimal power flow problem and has a good computational behavior. Simulation case studies verify the validity of the proposed work.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A nonsmooth nonconvex optimization algorithm for two-stage optimization problems

An optimization algorithm for a group of nonsmooth nonconvex problems inspired by two-stage stochastic programming problems is proposed. The main challenges for these problems include (1) the problems lack the popular lower-type properties such as prox-regularity assumed in many nonsmooth nonconvex optimization algorithms, (2) the objective can not be analytically expressed and (3) the evaluation of function values and subgradients are computationally expensive. To address these challenges, this report first examines the properties that exist in many two-stage problems, specifically upper-C 2 objectives. Then, we show that quadratic penalty method for securityconstrained alternating current optimal power flow (SCACOPF) contingency problems can make the contingency solution functions upper-C 2 . Based on these observations, a simplified bundle algorithm that bears similarity to sequential quadratic programming (SQP) method is proposed. It is more efficient in implementation and computation compared to conventional bundle methods. Global convergence analysis of the algorithm is presented under novel and reasonable assumptions. The proposed algorithm therefore fills the gap of theoretical convergence for smoothed SCACOPF problems. The inconsistency that might arise in our treatment of the constraints are addressed through a penalty algorithm whose convergence analysis is also provided. Finally, theoretical capabilities and numerical performance of the algorithm are demonstrated through numerical examples.

97 MATHEMATICS AND COMPUTING↗

Real-Time Alpine Measurement System Using Wireless Sensor Networks

Monitoring the snow pack is crucial for many stakeholders, whether for hydro-power optimization, water management or flood control. Traditional forecasting relies on regression methods, which often results in snow melt runoff predictions of low accuracy in non-average years. Existing ground-based real-time measurement systems do not cover enough physiographic variability and are mostly installed at low elevations. We present the hardware and software design of a state-of-the-art distributedWireless Sensor Network (WSN)-based autonomous measurement system with real-time remote data transmission that gathers data of snow depth, air temperature, air relative humidity, soil moisture, soil temperature, and solar radiation in physiographically representative locations. Elevation, aspect, slope and vegetation are used to select network locations, and distribute sensors throughout a given network location, since they govern snow pack variability at various scales. Three WSNs were installed in the Sierra Nevada of Northern California throughout the North Fork of the Feather River, upstream of the Oroville dam and multiple powerhouses along the river. The WSNs gathered hydrologic variables and network health statistics throughout the 2017 water year, one of northern Sierra’s wettest years on record. These networks leverage an ultra-low-power wireless technology to interconnect their components and offer recovery features, resilience to data loss due to weather and wildlife disturbances and real-time topological visualizations of the network health. Data show considerable spatial variability of snow depth, even within a 1 km 2 network location. In conclusion, combined with existing systems, these WSNs can better detect precipitation timing and phase in, monitor sub-daily dynamics of infiltration and surface runoff during precipitation or snow melt, and inform hydro power managers about actual ablation and end-of-season date across the landscape.

47 OTHER INSTRUMENTATION↗

Robe Development for Electrical Conductivity Analysis in an Electron Gun Produced Helium Plasma

The use of magnetohydrodynamic (MHD) power conversion systems, potentially coupled with a fission power source, is currently being investigated as a driver for an advanced propulsion system, such as a plasma thruster. The efficiency of a MHD generator is strongly dependent on the electrical conductivity of the fluid that passes through the generator; power density increases as fluid conductivity increases. Although traditional MHD flows depend on thermal ionization to enhance the electrical conductivity, ionization due to nuclear interactions may achieve a comparable or improved conductivity enhancement while avoiding many of the limitations inherent to thermal ionization. Calculations suggest that nuclear-enhanced electrical conductivity increases as the neutron flux increases; conductivity of pure He-3 greater than 10 mho/m may be achievable if exposed to a flux greater than 10(exp 12) neutrons/cm2/s.) However, this remains to be demonstrated experimentally. An experimental facility has been constructed at the Propulsion Research Center at the NASA Marshall Space Flight Center, using helium as the test fluid. High energy electrons will be used to simulate the effects of neutron-induced ionization of helium gas to produce a plasma. These experiments will be focused on diagnosis of the plasma in a virtually static system; results will be applied to future tests with a MHD system. Initial experiments will utilize a 50 keV electron gun that can operate at up to a current of 200 micro A. Spreading the electron beam over a four inch diameter window results in an electron flux of 1.5x 10(exp 13) e/sq cm/s. The equivalent neutron flux that would produce the same ionization fraction in helium is 1x10(exp 12) n/sq cm/s. Experiments will simulate the neutron generated plasma modeled by Bitteker, which takes into account the products of thermal neutron absorption in He-3, and includes various ion species in estimating the conductivity of the resulting plasma. Several different probes will be designed and implemented to verify the plasma kinetics model. System parameters and estimated operating ranges are summarized. The predicted ionization fraction, electron density, and conductivity levels are provided in for an equivalent neutron flux of 1x10(exp 12) n/cm2/s. Understanding the complex plasma kinetics throughout a MHD channel is necessary to design an optimal power conversion system for space propulsion applications. The proposed experiments seek to fully characterize the helium plasma and to determine the reliability of each measurement technique, such that they may be applied to more advanced MHD studies. The expected value of each plasma parameter determined from theoretical models will be verified experimentally by several independent techniques to determine the most reliable method of obtaining each parameter. The results of these experiments will be presented in the final paper.

Bragg-Sitton, Shannon M.↗

Decision Analysis Methods For the Analysis of Nuclear Terrorism Threats with Imperfect Information

Protection of commercial nuclear fuel cycles against adversary groups who may attempt to steal weapons-usable material is an issue of interest to the U.S. Department of Energy, Office of Nuclear Energy. The wide array of potential nuclear fuel cycle configurations and the security measures that could be deployed to protect them pose a daunting system design problem. In general, there are billions of potential design configurations. Moreover, the design choices are not independent – the effectiveness of one design feature or security measure can affect the desirability of other design choices. To explore this issue, we develop a mixed-integer optimization problem that solves a complex decision tree to allocate security resources across a commercial nuclear fuel cycle configuration. A nominal UREX+ reprocessing facility is used to demonstrate how the model allocates resources across a particular fuel cycle facility and shows how the problem can be scaled up to include other facilities and design features in the fuel cycle. In our model, the defender first chooses a set of safeguards to implement in order to mitigate the risk of a terrorist theft from various unit operations in the reprocessing facility. The attacker observes the defender’s strategy choice and executes an attack to acquire material, including the possibility of declining to attack at all. The model considers the terrorist as an intelligent adversary with a set of attack strategies known to the defender. The strategic choices made by both the defender and attacker affect the probability of a successful attack by the terrorist and in turn affect the distribution of consequences to the defender. We use a decision tree to illustrate the interactions among costs, decisions, correlated uncertainties, and outcomes, and then show how to formulate the problem as a mixed-integer program with an equivalent solution. The novel features of our formulation include the following. First, as discussed above, we use a sequential game-theoretic formulation (a Stackelberg game). This formulation assumes the attacker can observe the defender’s actions before the attack is executed. Many previous game-theoretic studies have assumed attacker and defender act simultaneously without knowledge of the other’s strategy choice. We believe the sequential game formulation is more applicable for knowledgeable insider adversaries, the focus of this study. Second, we build on our previous work that has developed a methodology to reflect the correlation among probability distributions represented as discrete chance nodes in a decision tree. We illustrate how this methodology can be applied to two chance nodes in a decision tree that relate the quality of the materials that might be stolen or diverted by an adversary with the yield of a weapon that might be constructed from that material. In a decision tree, the dependence of the yield on the quality of the material could be modeled with conditional probabilities of the weapon yield given the quality of the materials obtained. For three outcome discretizations, this would require the assessment of nine conditional probabilities from domain experts as part of the model building effort. In contrast, our methodology will allow this same representation of the tree based on the assessment of the correlation between these two chance nodes by the experts, which we believe is a much easier and more reliable assessment task. The advantage of this approach increases as the number of dependent chance nodes grows in a decision tree representation of this problem. Third, we extend this approach to the approximation of the dependence between discrete binary events such as success or failure represented by binary chance nodes in a decision tree or in the special case of a probabilistic risk analysis (PRA) based on an event tree. This extension is a novel application in this domain, and uses Monte Carlo simulation to determine the appropriate conditional probabilities on the binary chance nodes when the number of correlated uncertainties is more than two. Fourth, we show how these decision trees or probability trees can be formulated and solved as equivalent integer programming (IP) problems. The advantages of the IP representation include the opportunity to use the computing power associated with powerful optimization solver algorithms available on high speed computers to solve very large versions of these trees, and the ability to include several such models of different facilities into a much larger IP with connecting constraints that represent budget restrictions, number of available inspection opportunities in a given time window, or other dependencies among these facilities. Policy makers can use the resulting model to inform their decisions about how to assess terrorism risk and to safeguard the commercial nuclear fuel cycle.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thiol‐Functionalized Conjugated Metal–Organic Frameworks for Stable and Efficient Perovskite Photovoltaics

Abstract Metal–organic frameworks (MOFs) have been investigated recently in perovskite photovoltaics owing to their potential to boost optoelectronic performance and device stability. However, the impact of variations in the MOF side chain on perovskite characteristics and the mechanism of MOF/perovskite film formation remains unclear. In this study, three nanoscale thiol‐functionalized UiO‐66‐type Zr‐based MOFs (UiO‐66‐(SH) 2 , UiO‐66‐MSA, and UiO‐66‐DMSA) are systematically employed and examined in perovskite solar cells (PSCs). Among these MOFs, UiO‐66‐(SH) 2 , with its rigid organic ligands, exhibited a strong interaction with perovskite materials with more efficient suppression of perovskite vacancy defects. More importantly, A detailed and in‐depth discussion is provided on the formation mechanism of UiO‐66‐(SH) 2 ‐assisted perovskite film upon in situ GIWAXS performed during the annealing process. The incorporation of UiO‐66‐(SH) 2 additives substantially facilitates the conversion of PbI 2 into the perovskite phase, prolongs the duration of stage I, and induces a delayed phase transformation pathway. Consequently, the UiO‐66‐(SH) 2 ‐assisted device demonstrates reduced defect density and superior optoelectronic properties with optimized power conversion efficiency of 24.09% and enhanced long‐term stability under ambient environment and continuous light illumination conditions. This study acts as a helpful design guide for desired MOF/perovskite structures, enabling further advancements in MOF/perovskite optoelectronic devices.

14 SOLAR ENERGY↗

Conjugation‐Curtailing of Benzodithionopyran‐Cored Molecular Acceptor Enables Efficient Air‐Processed Small Molecule Solar Cells

Abstract Small molecule solar cells (SMSCs) lag a long way behind polymer solar cells. A key limit is the less controllable morphology of small molecule materials, which can be aggravated when incorporating anisotropic nonfullerene acceptors. To fine‐tune the blending morphology within SMSCs, a π‐conjunction curtailing design is applied, which produces a efficient benzodithionopyran‐cored molecular acceptor for nonfullerene SMSCs (NF‐SMSCs). When blended with a molecular donor BDT3TR‐SF to fabricate NF‐SMSCs, the π‐conjunction curtailed molecular acceptor NBDTP‐M obtains an optimal power conversion efficiency (PCE) of up to 10.23%, which is much higher than that of NBDTTP‐M of longer π‐conjunction. It retains 93% of the PCE of devices fabricated in a glove box when all spin‐coating and post‐treating procedures are conducted in ambient air with relative humidity of 25%, which suggests the good air‐processing capability of π‐conjunction curtailed molecules. Detailed X‐ray scattering investigations indicate that the BDT3TR‐SF:NBDTP‐M blend exhibits a blend morphology featuring fine interpenetrating networks with smaller domains and higher phase purity, which results in more efficient charge generation, more balanced charge transport, and less recombination compared to the low‐performance BDT3TR‐SF:NBDTTP‐M blend. This work provides a guideline for molecular acceptors' design toward efficient, low‐cost, air‐processed NF‐SMSCs.

Wu, Hao↗

Perfecting one-loop BCJ numerators in SYM and supergravity

We take a major step towards computing D-dimensional one-loop amplitudes in general gauge theories, compatible with the principles of unitarity and the color-kinematics duality. For n-point amplitudes with either supersymmetry multiplets or generic non-supersymmetric matter in the loop, simple all-multiplicity expressions are obtained for the maximal cuts of kinematic numerators of n-gon diagrams. At n = 6, 7 points with maximal supersymmetry, we extend the cubic-diagram numerators to encode all contact terms, and thus solve the long-standing problem of simultaneously realizing the following properties: color-kinematics duality, manifest locality, optimal power counting of loop momenta, quadratic rather than linearized Feynman propagators, compatibility with double copy as well as all graph symmetries. Color-kinematics dual representations with similar properties are presented in the half-maximally supersymmetric case at n = 4, 5 points. The resulting gauge-theory integrands and their supergravity counterparts obtained from the double copy are checked to reproduce the expected ultraviolet divergences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A simple way to integrate distributed storage into a wholesale electricity market

Abstract Current plans to decarbonize the electric supply system imply that the generation from wind and solar sources will grow substantially. This growth will increase the uncertainty of system operations due to the inherent variability of these renewable sources, and as a result, more reserve capacity will be required to provide the ramping (flexibility) needed for reliable operations. This paper assumes that all of the increased uncertainty comes from wind farms on the grid, and it shows how distributed storage managed locally by aggregators can provide the ramping needed without introducing a separate market for flexibility. This can be accomplished when the aggregators minimize the expected daily cost of the energy purchased from the grid for their customers by submitting optimal bids into the wholesale market with high and low price thresholds for discharging and charging the storage. This model is illustrated using a stochastic multi-period security constrained optimal power flow together with realistic data for a reduction of the network in the Northeast Power Coordinating Council region of the United States. The results show that the bidding strategy for distributed storage provides ramping to the grid just as effectively as storage managed by a system operator.

Lamadrid, Alberto J.↗

Mechanochemical method of producing triethoxysilane

A mechanochemical method for synthesis of triethoxysilane from silicon-copper contact mass and ethyl alcohol in the developed vibration reactor is presented. It is shown that the process of a direct alkoxysilane synthesis in the vibro-boiling layer is affected by a series of control parameters such as the ratio between the contact mass and the mass of grinding bodies, the grinding body sizes and their ratios in a polydisperse mixture, power density. Optimization of these parameters allowed us to obtain HSi(OEt){sub 3} with a selectivity of 50% at a silicon conversion of 90% without the use of promoters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of substitutional doping and disorder on the phase stability, magnetism, and half-metallicity of Heusler alloys

Spintronics is the fast growing field that will play a key role in optimizing power consumption, memory, and processing capabilities of nanoelectronic devices. Heusler alloys are potential candidates for application in spintronics due to their room temperature (RT) half-metallicity, high Curie temperature, low lattice mismatch with most substrates, and strong control on electronic density of states at Fermi level. In this work, we investigate the effect of substitutional doping and disorder on the half-metallicity, phase stability, and magnetism of Heusler alloys using density functional theory methods. Our study shows that electronic and magnetic properties of half/full-Heusler alloys can be tuned by changing electron-count through controlled variation of chemical compositions of alloying elements. We provide a detailed discussion on the effect of substitutional doping and disorder on the tunability of half-metallic nature of Co 2 MnX and NiMnX based Heusler alloys, where X represents group 13–16 and period 3–6 elements of the periodic table. Based on the idea of electron count and disorder, we predicted a possible existence of thermodynamically stable half-metallic multicomponent bismuthides, for example, (CuNi 3 )Mn 4 Bi 4 and (ZnNi 7 )Mn 8 Bi 8 , through substitution doping at Ni site by specific Cu and Zn composition in half-Heusler NiMnBi. We believe that the design guide based on electron-counts presented for half-metals will play a key role in electronic-structure engineering of novel Heusler alloys for spintronic application, which will accelerate the development and synthesis of novel materials.

36 MATERIALS SCIENCE↗

Combining model-based and model-free methods for stochastic control of distributed energy resources

Modern distribution systems are experiencing a fast transformation with the growing penetration of distributed energy resources (DERs). Along with the economic and environmental benefits of DERs, challenges arise to address the uncertainties caused by their inherent volatility. If properly coordinated, however, DERs have the potential to provide the controllability that grid operators need. Here, we propose a hierarchical control framework that combines the model-based and model-free methods for stochastic DER control in distribution systems. The upper-level scheduler considers a chance-constrained optimal power flow problem (model-based) that schedules DER setpoints to minimize the operational cost and maintain the operating reserve. The lower-level distributed DER controllers absorb real-time disturbances and uncertainties using the extremum seeking control (model-free) to achieve grid objectives. The combination of model-based and model-free methods allows us to take the advantages of both methods to effectively manage the uncertainty in distribution systems. The proposed work is demonstrated on the IEEE 13-node feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Self-aligned concentrating immersion-lens arrays for patterning and efficiency recovery in scaffold-reinforced perovskite solar cells

In contrast to existing silicon, CIGS, and multi-junction cells that exhibit remarkable mechanical durability, perovskite solar cells have been shown to delaminate when subjected to the mechanical loads that occur during processing and field exposures, limiting their potential as a reliable solar technology. Mechanical reinforcement is thus essential for durable and reliable perovskite technologies. One current strategy to overcome the thermomechanical fragility of perovskite solar cells is to extrinsically shield them by introducing reinforcing scaffolds which partition the cell into many distinct microcells, but improvements in mechanical stability coincide with a reduced device efficiency due to parasitic absorption by the scaffold. We address this reduced efficiency by integrating concentrating immersion-lens arrays (CILAs) into scaffold-reinforced solar cells. These scaffolds are lithographically patterned by a maskless process, in which ultraviolet light is used to pattern the scaffolds through the CILAs, ensuring self-alignment and optical contact with the microcells. Perovskite devices are deposited into the microcells, and during operation, light is concentrated into the microcells and away from the insulating scaffolds, resulting in mechanically resilient solar cells with efficiencies comparable to planar devices. The CILAs also exhibit passive tracking of incident light—verified experimentally and with ray-tracing simulations—which is critical for optimal power output as the sun moves across the sky during the day. The simplicity of the fabrication process and the efficiency of the resulting scaffolded devices show that a lens-integrated scaffold-reinforced structure is a potential pathway to high-performance, robust, commercially viable perovskite solar cells.

14 SOLAR ENERGY↗

Scalable Preconditioning of Block-Structured Linear Algebra Systems using ADMM

In this paper, we study the solution of block-structured linear algebra systems arising in optimization by using iterative solution techniques. These systems are the core computational bottleneck of many problems of interest such as parameter estimation, optimal control, network optimization, and stochastic programming. Our approach uses a Krylov solver (GMRES) that is preconditioned with an alternating method of multipliers (ADMM). We show that this ADMM-GMRES approach overcomes well-known scalability issues of Schur complement decomposition in problems that exhibit a high degree of coupling. The effectiveness of the approach is demonstrated using linear systems that arise in stochastic optimal power flow problems and that contain up to 2 million total variables and 4,000 coupling variables. We find that ADMM-GMRES is nearly an order of magnitude faster than Schur complement decomposition. Moreover, we demonstrate that the approach is robust to the selection of the augmented Lagrangian penalty parameter, which is a key advantage over the direct use of ADMM.

97 MATHEMATICS AND COMPUTING↗