Engineering topics

Maldonado, Daniel Adrian

Publications and source records attributed to Maldonado, Daniel Adrian.

8 records

DOE OSTI2024

Parallel interior-point solver for block-structured nonlinear programs on SIMD/GPU architectures

Here, we investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses SIMD/GPU accelerators locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations. We demonstrate our method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the reduction method is 50x faster than the sparse linear solver HSL MA57 running in serial on the CPU, and 6x faster than Pardiso running in parallel on CPU on the same number of processes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Distributed Frequency Divider for Power System Bus Frequency Online Estimation Considering Virtual Inertia From DFIGs

In this work, a distributed frequency divider is proposed to estimate power system bus frequency with a limited number of PMUs while considering the inertial contributions from double-fed induction generators (DFIGs). The key idea is to reformulate the original frequency divider by modeling the contributions of DFIGs with inertia emulation and external regional power system through Thevenin equivalents. The distributed frequency divider is general and able to estimate local bus frequencies in a distributed manner. Besides, only the knowledge of the rotor speeds of synchronous generator, boundary bus frequencies and terminal bus frequencies of DFIGs as well as the admittance matrix are required to estimate the frequencies at all buses. This is drastically different from existing approaches that require system observability by PMUs to monitor all bus frequencies. Numerical results carried out on the IEEE 39-bus and modified 118-bus power systems with DFIGs demonstrate and effectiveness and robustness of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Trust-Region Approximation of Extreme Trajectories in Power System Dynamics

In this work we present a novel technique, based on a trust-region optimization algorithm and second-order trajectory sensitivities, to compute the extreme trajectories of power system dynamic simulations given a bounded set that represents parametric uncertainty. Furthermore, we show how this method, while remaining computationally efficient compared with sampling-based techniques, overcomes the limitations of previous sensitivity-based techniques to approximate the bounds of the trajectories when the local approximation loses validity because of the nonlinearity. We present several numerical experiments that showcase the accuracy and scalability of the technique, including a demonstration on the IEEE New England test system.

42 ENGINEERING↗

DOE OSTI2021

Extended Frequency Divider for Bus Frequency Estimation Considering Virtual Inertia from DFIGs

Accurate estimation of local bus frequency is important for effectively controlling both synchronous and nonsynchronous generators. As the power grid evolves to accommodate essential reliability services such as virtual inertia from nonsynchronous generators, conventional techniques to estimate frequency face challenges. This paper proposes a new frequency estimation method that can effectively include the inertia contributions from double-fed wind generators (DFIGs). This is achieved through the proposed extended frequency divider formula (FDF) to include the contributions of DFIGs via Thevenin equivalents. The proposed extended FDF does not suffer from numerical issues as compared to existing phasor angle derivative-based approaches. Moreover, the knowledge of the rotor speeds of the synchronous machine and DFIGs as well as of the network admittance matrix allows estimating the frequencies of all buses in the grid, thereby significantly improving system situational awareness with a limited number of measurements. Numerical results on the IEEE 39-bus power system with DFIGs show that the proposed method achieves more accurate bus frequency estimations than the original FDF formula and other approaches based on the numerical derivation of the bus voltage phase angles.

double-fed induction generator↗

DOE OSTI2020

Sequential Bayesian Parameter Estimation of Stochastic Dynamic Load Models

In this paper we focus on the parameter estimation of dynamic load models with stochastic terms-in particular, load models where protection settings are uncertain, such as in aggregated air conditioning units. We show how the uncertainty in the aggregated protection characteristics can be formulated as a stochastic differential equation with process noise. We cast the parameter inversion within a Bayesian parameter estimation framework, and we present methods to include process noise. We demonstrate the benefits of considering stochasticity in the parameter estimation and the risks of ignoring it.

Bayesian Statistics↗

DOE OSTI2020

Probabilistic analysis of masked loads with aggregated photovoltaic production

In this paper we present a probabilistic analysis framework to estimate behind-the-meter photovoltaic generation in real time. We develop a forward model consisting of a spatiotemporal stochastic process that represents the photovoltaic generation and a stochastic differential equation with jumps that represents the demand. Here, we employ this model to disaggregate the behind-the-meter photovoltaic generation using net load and irradiance measurements.

14 SOLAR ENERGY↗

DOE OSTI2020

PETSc DMNetwork: A Library for Scalable Network PDE-Based Multiphysics Simulations

We present DMNetwork, a high-level package included in the PETSc library for the simulation of multiphysics phenomena over large-scale networked systems. The library aims at applications that have networked structures such as those in electrical, gas, and water distribution systems. DMNetwork provides data and topology management, parallelization for multiphysics systems over a network, and hierarchical and composable solvers to exploit the problem structure. DMNetwork eases the simulation development cycle by providing the necessary infrastructure through simple abstractions to define and query the network components. This article presents the design of DMNetwork, illustrates its user interface, and demonstrates its ability to solve multiphysics systems, such as an electric circuit, a network of power grid and water subnetworks, and transient hydraulic systems over large networks with more than 2 billion variables on extreme-scale computers using up to 30,000 processors.

extreme-scale computing↗

DOE OSTI2019

A Framework for Distributed Approximation of Moments with Higher-Order Derivatives through Automatic Differentiation

We present a framework for the distributed approximation of moments, enabling the evaluation of the uncertainty in a dynamical system. The first and second moment, mean, and variance are computed with up to third-order Taylor series expansion. The required derivatives for the expansion are generated automatically by automatic differentiation and propagated through an implicit time stepper. The computational kernels are the accumulation of the derivatives (Jacobian, Hessian, tensor) and the covariance matrix. We apply distributed parallelism to the Hessian or third-order tensor, and the user merely has to provide a function for the differential equation, thus achieving similar ease of use as Monte Carlo-based methods. We demonstrate our approach using with benchmarks on Theta, a KNL-based system at the Argonne Leadership Computing Facility.

Automatic Differentiation↗