Refine
Year of publication
Keywords
- Bilevel optimization (16)
- Gas networks (8)
- Optimal control (8)
- Robust optimization (8)
- Uniqueness (7)
- optimal control (7)
- Mixed-Integer Nonlinear Optimization (6)
- Mixed-integer optimization (6)
- Convergence (5)
- Gas Networks (5)
We consider uncertain robust electricity market equilibrium problems including transmission and generation investments. Electricity market equilibrium modeling has a long tradition but is, in most of the cases, applied in a deterministic setting in which all data of the model are known. Whereas there exist some literature on stochastic equilibrium problems, the field of robust equilibrium models is still in its infancy. We contribute to this new field of research by considering Γ-robust electricity market equilibrium models on lossless DC networks with transmission and generation investments. We state the nominal market equilibrium problem as a mixed complementarity problem as well as its variational inequality and welfare optimization counterparts. For the latter, we then derive a Γ-robust formulation and show that it is indeed the counterpart of a market equilibrium problem with robustified player problems. Finally, we present two case studies to gain insights into the general effects of robustification on electricity market models. In particular, our case studies reveal that the transmission system operator tends to act more risk-neutral in the robust setting, whereas generating firms clearly behave more risk-averse.
In these lecture notes, we address the problem of large-time asymptotic behaviour of the solutions to scalar convection-diffusion equations set in [katex]\mathbb{R}^N[/katex]. The large-time asymptotic behaviour of the solutions to many convection-diffusion equations is strongly linked with the behavior of the initial data at infinity. In fact, when the initial datum is integrable and of mass [katex]M[/katex], the solutions to the equations under consideration oftentimes behave like the associated self-similar profile of mass [katex]M[/katex], thus emphasising the role of scaling variables in these scenarios. However, these equations can also manifest other asymptotic behaviors, including weakly non-linear, linear or strongly non-linear behavior depending on the form of the convective term. We give an exhaustive presentation of several results and techniques, where we clearly distinguish the role of the spatial dimension and the form of the nonlinear convective term.
This paper is devoted to analysing the explicit slow decay rate and turnpike in the infinite-horizon linear quadratic optimal control problems for hyperbolic systems. Assume that some weak observability or controllability are satisfied, by which, the lower and upper bounds of the corresponding algebraic Riccati operator are estimated, respectively. Then based on these two bounds, the explicit slow decay rate of the closed-loop system with Riccati-based optimal feedback control is obtained. The averaged turnpike property for this problem is also further discussed.
We then apply these results to the LQ optimal control problems constraint to networks of onedimensional wave equations and also some multi-dimensional ones with local controls which lack of GCC (Geometric Control Condition).
In real life decision problems, one almost always is confronted with uncertainty and risk. For practical optimization problems this is manifested by unknown parameters within the input data, or, an inexact knowledge about the system description itself. In case the uncertain problem data is governed by a known probability distribution, stochastic programming offers a variety of models hedging against uncertainty and risk. Most widely employed are two-stage models, who admit a recourse structure: The first-stage decisions are taken before the random event occurs. After its outcome, a recourse (second-stage) action is made, often but not always understood as some "compensation".
In the present thesis, the optimization problems that involve parameters which are not known with certainty are semidefinite programming problems. The constraint sets of these optimization problems are given by intersections of the cone of symmetric, positive semidefinite matrices with either affine or more general equations. Objective functions, formally, may be fairly general, although they often are linear as in the present thesis.
We consider risk neutral and risk averse two-stage stochastic semidefinite programs with continuous and mixed-integer recourse, respectively. For these stochastic optimization problems we analyze their structure, derive solution methods relying on decomposition, and finally apply our results to unit commitment in alternating current (AC) power systems.
Furthermore, deterministic unit commitment in AC power transmission systems is addressed. Beside traditional unit commitment constraints, the physics of power flow are included. To gain globally optimal solutions a recent semidefinite programming (SDP) approach is used which leads to large-scale semidefinite programs with discrete variables on top. As even the SDP relaxation of these programs is too large for being handled in an all-at-once manner by general SDP solvers, it requires an efficient and reliable method to tackle them. To this end, an algorithm based on Benders decomposition is proposed.
With power demand (load) and in-feed from renewables serving as sources of uncertainty, two-stage stochastic programs are set up heading for unit commitment schedules which are both cost-effective and robust with respect to data perturbations. The impact of different, risk neutral and risk averse, stochastic criteria on the shapes of the optimal stochastic solutions will be examined. To tackle the resulting two-stage programs, we propose to approximate AC power flow by semidefinite relaxations. This leads to two-stage stochastic mixed-integer semidefinite programs having a special structure. To solve the latter, the L-shaped method and dual decomposition have been applied and compared.
Finite Time Blow-up of Traveling Wave Solutions for the Flow of Real Gas through Pipeline Networks
(2016)
In the context of gas transportation, analytical solutions are essential
for the understanding of the underlying dynamics described
by a system of partial differential equations. We derive traveling wave
solutions for the 1-d isothermal Euler equations. A non-constant compressibility
factor is used to describe the correlation between density
and pressure. The blow-up of the traveling wave solution in finite time
is proven. We then extend our analysis to networks under appropriate
coupling conditions and derive compatibility conditions to fulfill these
coupling conditions.
We consider optimization problems with a joint probabilistic constraint under
normally distributed uncertain parameters. The parametric constraints
are replaced by one constraint stating that the probability of being feasible
shall exceed or be equal to a prescribed threshold. In order to apply the
concept to gas network optimization under uncertain boundary flows, which
corresponds to the demand of customers, we derive an analytic gradient formula.
The integral corresponding to the probability can be parameterized by
spherical radial decomposition. For this parameterization gradient formulas
are known under convexity assumptions of the parametric constraints in the
parameter. For the application in gas networks that we have in mind, the
convexity assumption of the parametric constraints is not satisfied. Therefore,
we weaken it to convexity of the region of feasible parameters for a fixed
optimization variable. We proceed to show that the assumptions needed for
the gradient formula are met in the gas network optimization problem on
a tree. For the numerical implementation we propose a multilevel sampling
algorithm that uses a coarse approximation of the chance constraint to generate
a warm start for the expensive approximation with fine sampling. The
numerical results illustrate that this approach significantly reduces the computation
time.
High-dimensional interpolation problems appear in various applications of uncertainty quantification, stochastic optimization and machine learning. Such problems are computationally expensive and request the use of adaptive grid generation strategies like anisotropic sparse grids to mitigate the curse of dimensionality. However, it is well known that the standard dimension-adaptive sparse grid method converges very slowly or even fails in the case of non-smooth functions. For piecewise smooth functions with kinks, we construct two novel hp-adaptive sparse grid collocation algorithms that combine low-order basis functions with local support in parts of the domain with less regularity and variable-order basis functions elsewhere. Spatial refinement is realized by means of a hierarchical multivariate knot tree which allows the construction of localised hierarchical basis functions with varying order. Hierarchical surplus is used as an error indicator to automatically detect the non-smooth region and adaptively refine the collocation points there. The local polynomial degrees are optionally selected by a greedy approach or a kink detection procedure. Three numerical benchmark examples with different dimensions are discussed and comparison with locally linear and highest degree basis functions are given to show the efficiency and accuracy of the proposed methods.
We build up a decomposition for the flow generated by the heat equation with a real analytic memory kernel. It
consists of three components: The first one is of parabolic nature; the second one gathers the hyperbolic component
of the dynamics, with null velocity of propagation; the last one exhibits a finite smoothing effect. This decomposition reveals the hybrid parabolic-hyperbolic nature of the flow and clearly illustrates the significant impact of the memory term on the parabolic behavior of the system in the absence of memory terms.
Compressor machines are crucial elements in a gas transmission network, required to compensate for the pressure loss caused by friction in the pipes. Modelling all physical and technical details of a compressor machine involves a large amount of nonlinearity, which makes it hard to use such models in the optimization of large-scale gas networks. In this paper, we are going to describe a modelling approach for the operating range of a compressor machine, starting from a physical reference model and resulting in a polyhedral representation in the 3D space of mass flow throughput as well as in- and outlet pressure.
Gas networks are an important application area for optimization. When considering long-range transmission, compressor stations play a crucial role in these applications. The purpose of this report is to collect and systematize the models used for compressor stations in the literature. The emphasis is on recent work on simple yet accurate polyhedral models that may replace more simplified traditional models without increasing model complexity. The report also describes an extension of the compressor station data available in GasLib (http://gaslib.zib.de/) with the parameters of these models.