Friedrich-Alexander-Universität Erlangen-Nürnberg
Refine
Year of publication
Keywords
- Bilevel optimization (8)
- Optimal control (7)
- Gas networks (6)
- Uniqueness (6)
- Convergence (5)
- Networks (5)
- Robust optimization (5)
- robust optimization (5)
- Branch-and-cut (4)
- Mixed-Integer Nonlinear Optimization (4)
Joint model of probabilistic/robust (probust) constraints applied to gas network optimization
(2017)
Optimization tasks under uncertain conditions abound in many
real-life applications. Whereas solution approaches for probabilistic constraints
are often developed in case the uncertainties can be assumed to follow a
certain probability distribution, robust approaches are usually used in case
solutions are sought that are feasible for all realizations of uncertainties within
some pre-defined uncertainty set. As many applications contain different types
of uncertainties that require robust as well as probabilistic treatments, we deal with a class of joint probabilistic/robust constraints as its appears in
optimization problems under uncertainty. Focusing on complex uncertain gas
network optimization problems, we show the relevance of this class of problems
for the task of maximizing free booked capacities in an algebraic model for a
stationary gas network. We furthermore present approaches for their solution.
Finally, we study the problem of controlling a transient system that is governed
by the wave equation. The task consists in determining controls such that a
certain robustness measure remains below some given upper bound, with high
probability.
Pareto efficiency for robust linear programs was introduced by
Iancu and Trichakis in [9]. We generalize their approach and theoretical results
to robust optimization problems in Euclidean spaces with linear uncertainty.
Additionally, we demonstrate the value of this approach in an exemplary
manner in the area of robust semidefinite programming (SDP). In particular,
we prove that computing a Pareto robustly optimal solution for a robust SDP
is tractable and illustrate the benefit of such solutions at the example of the
maximal eigenvalue problem. Furthermore, we modify the famous algorithm of
Goemans and Williamson [8] in order to compute cuts for the robust max cut
problem that yield an improved approximation guarantee in non-worst-case
scenarios.
We present a solution framework for general alternating current optimal power flow (AC OPF) problems that include discrete decisions.
The latter occur, for instance, in the context of the curtailment of renewables or the
switching of power generation units and transmission lines.
Our approach delivers globally optimal solutions and is provably convergent.
We model AC OPF problems with discrete decisions as mixed-integer nonlinear programs.
The solution method starts from a known framework that uses piecewise linear relaxations.
These relaxations are modeled as as mixed-integer linear programs and adaptively refined until some termination criterion is fulfilled.
In this work, we extend and complement this approach by problem-specific as well as very general algorithmic enhancements.
In particular, these are mixed-integer second-order cone programs as well as primal and dual cutting planes.
For example objective cuts and no-good-cuts help to compute good feasible solutions as where outer approximation constraints tighten the relaxations.
We present extensive numerical results for various AC OPF problems where discrete decisions play a major role.
Even for hard instances with a large proportion of discrete decisions, the method is able
to generate high quality solutions efficiently.
Furthermore, we compare our approach with state-of-the-art MINLP.
Our method outperforms all other algorithms.
Stochastic Optimization (SO) is a classical approach for optimization under uncertainty that typically requires knowledge about the probability distribution of uncertain parameters. As the latter is often unknown, Distributionally Robust Optimization (DRO) provides a strong alternative that determines the best guaranteed solution over a set of distributions (ambiguity set). In this work, we present an approach for DRO over time that uses online learning and scenario observations arriving as a data stream to learn more about the uncertainty. Our robust solutions adapt over time and reduce the cost of protection with shrinking ambiguity. For various kinds of ambiguity sets, the robust solutions converge to the SO solution. Our algorithm achieves the optimization and learning goals without solving the DRO problem exactly at any step. We also provide a regret bound for the quality of the online strategy which converges at a rate of $ O(\log T / \sqrt{T})$, where $T$ is the number of iterations. Furthermore, we illustrate the effectiveness of our procedure by numerical experiments on mixed-integer optimization instances from popular benchmark libraries and give practical examples stemming from telecommunications and routing. Our algorithm is able to solve the DRO over time problem significantly faster than standard reformulations.
We propose a mathematical optimization model and its solution for joint chance constrained DC Optimal Power Flow. In this application, it is particularly important that there is a high probability of transmission limits being satisfied, even in the case of uncertain or fluctuating feed-in from renewable energy sources. In critical network situations where the network risks overload, renewable energy feed-in has to be curtailed by the transmission system operator (TSO). The TSO can reduce the feed-in in discrete steps at each network node. The proposed optimization model minimizes curtailment while ensuring that there is a high probability of transmission limits being maintained. The latter is modeled via (joint) chance constraints that are computationally challenging. Thus, we propose a solution approach based on the robust safe approximation of these constraints. Hereby, probabilistic constraints are replaced by robust constraints with suitably defined uncertainty sets constructed from historical data. The uncertainty sets are calculated by encompassing randomly drawn scenarios using the scenario approach proposed by Margellos et al. (IEEE Transactions on Automatic Control, 59 (2014)). The ability to discretely control the power feed-in then leads to a robust optimization problem with decision-dependent uncertainties, i.e. the uncertainty sets depend on decision variables. We propose an equivalent mixed-integer linear reformulation for box uncertainties with the exact linearization of bilinear terms. Finally, we present numerical results for different test cases from the Nesta archive, as well as for a real network. We consider the discrete curtailment of solar feed-in, for which we use real-world weather and network data. The experimental tests demonstrate the effectiveness of this method and run times are very fast. Moreover, on average the calculated robust solutions lead only to a small increase in curtailment, when compared to nominal solutions.
Robust DC Optimal Power Flow with Modeling of Solar Power Supply Uncertainty via R-Vine Copulas
(2021)
We present a robust approximation of joint chance constrained DC Optimal Power Flow in combination with a model-based prediction of uncertain power supply via R-vine copulas.
It is applied to optimize the discrete curtailment of solar feed-in in an electrical distribution network and guarantees network stability under fluctuating feed-in.
This is modeled by a two-stage mixed-integer stochastic optimization problem proposed by Aigner et al. (European Journal of Operational Research, (2021)).
The solution approach is based on the approximation of chance constraints via robust constraints using suitable uncertainty sets.
The resulting robust optimization problem has a known equivalent tractable reformulation.
To compute uncertainty sets that lead to an inner approximation of the stochastic problem, an R-vine copula model is fitted to the distribution of the multi-dimensional power forecast error, i.e., the difference between the forecasted solar power and the measured feed-in at several network nodes.
The uncertainty sets are determined by encompassing a sufficient number of samples drawn from the R-vine copula model.
Furthermore, an enhanced algorithm is proposed to fit R-vine copulas which can be used to draw conditional samples for given solar radiation forecasts.
The experimental results obtained for real-world weather and network data demonstrate the effectiveness of the combination of stochastic programming and model-based prediction of uncertainty via copulas.
We improve the outcomes of previous work by showing that the resulting uncertainty sets are much smaller and lead to less conservative solutions while maintaining the same probabilistic guarantees.
Ongoing policy discussions on the reconfiguration of bidding zones in European electricity markets induce uncertainty about the future market design. This paper deals with the question of how this uncertainty affects market participants and their long-run investment decisions in generation and transmission capacity. Generalizing the literature on pro-active network expansion planning, we propose a stochastic multilevel model which incorporates generation capacity investment, network expansion, and market operation, taking into account uncertainty about the future bidding zone configuration. Using a stylized two-node network, we disentangle different effects that uncertainty has on market outcomes. If there is a possibility that future bidding zone configurations provide improved regional price signals, welfare gains materialize even if the change does not actually take place. As a consequence, welfare gains of an actual change of the bidding zone configuration are substantially lower due to those anticipatory effects. Additionally, we show substantial distributional effects in terms of both expected gains and risks, between producers and consumers and between different generation technologies.
In the course of the energy transition, load and supply centers are growing apart in electricity markets worldwide, rendering regional price signals even more important to provide adequate locational investment incentives. This paper focuses on electricity markets that operate under a zonal pricing market design. For a fixed number of zones, we endogenously derive the optimal configuration of price zones and available transfer capacities on a network in order to optimally govern investment and production decisions in the long run. In a multilevel mixed-integer nonlinear model that contains a graph partitioning problem on the first level, we determine welfare-maximizing price zones and available transfer capacities for a given electricity market and analyze their impact on market outcomes. Using a generalized Benders decomposition approach developed in Grimm et al. (2019) and a problem-tailored scenario clustering for reducing the input data size, we are able to solve the model to global optimality even for large instances. We apply the approach to the German electricity market as an example to examine the impact of optimal zoning on key performance indicators such as welfare, generation mix and locations, or electricity prices. It turns out that even for a small number of price zones, an optimal configuration of zones induces a welfare level that almost approaches the first best.
We study gas network problems with compressors and control valves under uncertainty that can be formulated as two-stage robust optimization problems. Uncertain data are present in the physical parameters of the pipes as well as in the overall demand. We show how to exploit the special decomposable structure of the problem in order to reformulate the two-stage robust problem as a standard single-stage optimization problem. Since this structure is present in similar problems on e.g., water or direct current electricity networks, we investigate the consequences of the decomposable structure in an abstract setting: The right-hand side of the single-stage problem can be precomputed by solving a series of optimization problems and multiple elements of the right-hand side can be combined into one optimization task. In order to apply our results to gas network problems, we extend piecewise relaxations and preprocessing techniques to incorporate uncertain input data. The practical feasibility and effectiveness of our approach is demonstrated with benchmarks on realistic gas network instances. We observe large speedups due to the described aggregation method together with the developed preprocessing strategies. Furthermore, we are able to solve even comparably large gas network instances quickly for the price of slightly more conservative solutions.
In this paper we study feasibility and infeasibility of nonlinear two-stage
fully adjustable robust feasibility problems with an empty first stage. This
is equivalent to deciding set containment of a projection of the feasible
region and the uncertainty set. For answering this question, two very general
approaches using methods from polynomial optimization are presented --- one
for showing feasibility and one for showing infeasibility. The developed
methods are approximated through sum of squares polynomials and solved using
semidefinite programs.
Deciding robust feasibility and infeasibility is important for gas network
operations, which is a \nonconvex quadratic problem with absolute values
functions. Concerning the gas network problem, different topologies are
considered. It is shown that a tree structured network can be decided exactly
using linear programming. Furthermore, a method is presented to reduce a tree
network with one additional arc to a single cycle network. In this case,
removing the absolute values and solving the problem can be decided with
linearly many polynomial optimization problems.
Lastly, the effectivity of the methods is tested on a variety of small cyclic
networks. For instances where robust feasibility or infeasibility can be
decided, level~2 or level~3 of the Lasserre relaxation hierarchy is typically
sufficient.
This paper deals with the averaged dynamics for heat equations in the degenerate case where the diffusivity coefficient, assumed to be constant, is allowed to take the null value. First we prove that the averaged dynamics is analytic. This allows to show that, most often, the averaged dynamics enjoys the property of unique continuation and is approximately controllable. We then determine if the averaged dynamics is actually null controllable or not depending on how the density of averaging behaves when the diffusivity vanishes. In the critical density threshold the dynamics of the average is similar to the \frac{1}{2}-fractional Laplacian, which is wellknown to be critical in the context of the controllability of fractional diffusion processes. Null controllability then fails (resp. holds) when the density weights more (resp. less) in the null diffusivity regime than in this critical regime.
The aim of this work is to give a broad panorama of the control properties of fractional diffusive models from a numerical analysis and simulation perspective. We do this by surveying several research results we obtained in the last years, focusing in particular on the numerical computation of controls, though not forgetting to recall other relevant contributions which can be currently found in the literature of this prolific field. Our reference model will be a non-local diffusive dynamics driven by the fractional Laplacian on a bounded domain ΩΩΩ. The starting point of our analysis will be a Finite Element approximation for the associated elliptic model in one and two space-dimensions, for which we also present error estimates and convergence rates in the L2L^2L2 and energy norm. Secondly, we will address two specific control scenarios: firstly, we consider the standard interior control problem, in which the control is acting from a small subset ω⊂Ωω ⊂ Ωω⊂Ω. Secondly, we move our attention to the exterior control problem, in which the control region O⊂ΩcO ⊂ Ω cO⊂Ωc is located outside ΩΩΩ. This exterior control notion extends boundary control to the fractional framework, in which the non-local nature of the models does not allow for controls supported on ∂Ω∂Ω∂Ω. We will conclude by discussing the interesting problem of simultaneous control, in which we consider families of parameter-dependent fractional heat equations and we aim at designing a unique control function capable of steering all the different realizations of the model to the same target configuration. In this framework, we will see how the employment of stochastic optimization techniques may help in alleviating the computational burden for the approximation of simultaneous controls. Our discussion is complemented by several open problems related with fractional models which are currently unsolved and may be of interest for future investigation.
We discuss the multilevel control problem for linear dynamical systems, consisting in designing a piece-wise constant control function taking values in a finite-dimensional set. In particular, we provide a complete characterization of multilevel controls through a duality approach, based on the minimization of a suitable cost functional. In this manner we build optimal multi-level controls and characterize the time needed for a given ensemble of levels to assure the controllability of the system. Moreover, this method leads to efficient numerical algorithms for computing multilevel controls.
Auctions are widely used to determine the remuneration for renewable energies. They typically induce a high concentration of renewable energy plants at very productive sites far-off the main load centres, leading to an inefficient allocation as transmission line capacities are restricted but not considered in the allocation, resulting in an inefficient system configuration in the long run. To counteract these tendencies effectively, we propose a combinatorial auction design that allows to implement regional target capacities, provides a simple pricing rule and maintains a high level of competition between bidders by permitting package bids. By means of extensive numerical experiments we evaluate the combinatorial auction as compared to three further RES auction designs, the current German nationwide auction design, a simple nationwide auction, and regional auctions. We find that if bidders benefit from high enough economies of scale, the combinatorial auction design implements system-optimal target capacities without increasing the average remuneration per kWh as compared to the current German auction design. The prices resulting from the combinatorial auction are linear and anonymous for each region whenever possible, while minimal personalised markups on the linear prices are applied only when necessary. We show that realistic problem sizes can be solved in seconds, even though the problem is computationally hard.
We consider equilibrium problems under uncertainty where firms
maximize their profits in a robust way when selling their output. Robust
optimization plays an increasingly important role when best guaranteed objective
values are to be determined, independently of the specific distributional
assumptions regarding uncertainty. In particular, solutions are to be determined
that are feasible regardless of how the uncertainty manifests itself within
some predefined uncertainty set. Our analysis adopts the robust optimization
perspective in the context of equilibrium problems. First, we consider a singlestage,
nonadjustable robust setting. We then go one step further and study the
more complex two-stage or adjustable case where a part of the variables can
adjust to the realization of the uncertainty. We compare equilibrium outcomes
with the corresponding centralized robust optimization problem where the
sum of all profits are maximized. As we find, the market equilibrium for
the perfectly competitive firms differs from the solution of the robust central
planner, which is in stark contrast to classical results regarding the efficiency of
market equilibria with perfectly competitive firms. For the different scenarios
considered, we furthermore are able to determine the resulting price of anarchy.
In the case of non-adjustable robustness, for fixed demand in every time step
the price of anarchy is bounded whereas it is unbounded if the buyers are
modeled by elastic demand functions. For the two-stage adjustable setting,
we show how to compute subsidies for the firms that lead to robust welfare
optimal equilibria.
Linear complementarity problems are a powerful tool for modeling many practically relevant situations such as market equilibria. They also connect many sub-areas of mathematics like game theory, optimization, and matrix theory. Despite their close relation to optimization, the protection of LCPs against uncertainties - especially in the sense of robust optimization - is still in its infancy. During the last years, robust LCPs have only been studied using the notions of strict and Γ-robustness. Unfortunately, both concepts lead to the problem that the existence of robust solutions cannot be guaranteed. In this paper, we consider affinely adjustable robust LCPs. In the latter, a part of the LCP solution is allowed to adjust via a function that is affine in the uncertainty. We show that this notion of robustness allows to establish strong characterizations of solutions for the cases of uncertain matrix and vector, separately, from which existence results can be derived. Our main results are valid for the case of an uncertain LCP vector. Here, we additionally provide sufficient conditions on the LCP matrix for the uniqueness of a solution. Moreover, based on characterizations of the affinely adjustable robust solutions, we derive a mixed-integer programming formulation that allows to solve the corresponding robust counterpart. If, in addition, the certain LCP matrix is positive semidefinite, we prove polynomial-time solvability and uniqueness of robust solutions. If the LCP matrix is uncertain, characterizations of solutions are developed for every nominal matrix, i.e., these characterizations are, in particular, independent of the definiteness of the nominal matrix. Robust solutions are also shown to be unique for positive definite LCP matrix but both uniqueness and mixed-integer programming formulations still remain open problems if the nominal LCP matrix is not positive definite.
Multi-modal distributed energy system planning is applied in the context of smart grids, industrial energy supply,and in the building energy sector. In real-world applications, these systems are commonly characterized by existing system structures of different age where monitoring and investment are conducted in a closed-loop, with the iterative possibility to invest. The literature contains two main approaches to approximate this computationally intensive multiperiod
investment problem. The first approach simplifies the temporal decision-making process collapsing the multistage decision to a two-stage decision, considering uncertainty in the second stage decision variables. The second approach considers multi-period investments under the assumption of perfect foresight. In this work, we propose a
multi-stage stochastic optimization problem that captures multi-period investment decisions under uncertainty and solves the problem to global optimality, serving as a first-best benchmark to the problem. To evaluate the performance of conventional approaches applied in a multi-year setup and to solve the multi-period problem at lower computational effort, we propose a rolling horizon heuristic that on the one hand reveals the performance of conventional approaches
applied in a multi-period set-up and on the other hand enables planners to identify approximate solutions to the original
multi-stage stochastic problem. Additionally, we consider an open-loop version of the rolling horizon algorithm to evaluate how single-period investments perform with respect to the entire scenario tree and compared to multi-period investments.
We conduct a real-world case study and investigate solution quality as well as the computational performance of the proposed approaches. Our findings indicate that the approximation of multi-period investments by two-stage stochastic approaches yield the best results regarding constraint satisfaction, while deterministic multi-period approximations yield
better economic and computational performance.
To support the uprise of demand response, especially in the context of industrial processes, we propose a new approach to integrally determine the production-inventory plan and the cost-minimizing bids to participate in sequential reserve and energy-only markets. In particular, our approach considers time-coupling constraints which occur in the context of a production-inventory planning problem. We extend this problem with a comprehensive bidding formulation, which allows evaluating revenues and potential cost from the market participation, considering price uncertainties and uncertain activations of committed reserve capacity. This results in a multistage stochastic mixed-integer linear program, which explicitly considers the stage-wise revelation of information in our setup. To illustrate the capabilities of our approach, we apply our model to a real-world case study in which we investigate the participation of a cement plant in the German energy-only and reserve markets. The results of our case study indicate significant revenues for flexible industrial processes when participating in German spot and reserve markets.
We extend the framework of sequential action control to systems of partial differential equations which can be posed as abstract linear control problems in a Hilbert space. We follow a late-lumping approach and show that the control action can be explicitly obtained from variational principles using adjoint information. Moreover, we analyze the closed-loop system obtained from the SAC feedback for quadratic stage costs. We apply this theory prototypically to an unstable heat equation and verify the results numerically.
The aim of this paper is to study the derivation of appropriate meso- and macroscopic
models for interactions as appearing in social processes. There are two main characteristics
the models take into account, namely a network structure of interactions, which we treat
by an appropriate mesoscopic description, and a different role of interacting agents. The
latter differs from interactions treated in classical statistical mechanics in the sense that
the agents do not have symmetric roles, but there is rather an active and a passive agent.
We will demonstrate how a certain form of kinetic equations can be obtained to describe
such interactions at a mesoscopic level and moreover obtain macroscopic models from
monokinetics solutions of those.
The derivation naturally leads to systems of nonlocal reaction-diffusion equations (or
in a suitable limit local versions thereof), which can explain spatial phase separation
phenomena found to emerge from the microscopic interactions. We will highlight the
approach in three examples, namely the evolution and coarsening of dialects in human
language, the construction of social norms, and the spread of an epidemic.