Refine
Year of publication
Language
- English (1103) (remove)
Keywords
- optimal control (27)
- stability (14)
- integer programming (11)
- Stochastic programming (9)
- finite elements (9)
- mixed integer programming (9)
- Hamiltonian matrix (8)
- finite element method (8)
- model reduction (8)
- state constraints (8)
Improved optimization models for potential-driven network flow problems via ASTS orientations
(2019)
The class of potential-driven network flow problems provides important models for a range of infrastructure networks. For real-world applications, they need to be combined with integer models for switching certain network elements, giving rise to hard-to-solve MINLPs. We observe that on large-scale real-world meshed networks the relaxations usu-
ally employed are rather weak due to cycles in the network. To address this situation, we introduce the concept of ASTS orientations, a generalization of bipolar orientations, as a combinatorial relaxation of feasible solutions of potential-driven flow problems, study their structure, and
show how they can be used to strengthen existing relaxations and thus provide improved optimization models. Our computational results indicate that ASTS orientations can be used to derive much stronger bounds on the flow variables than existing bound tightening methods. We also
show that our proposed model extensions yield significant performance improvements for an existing state-of-the-art MILP model for large-scale gas networks.
In commodity transport networks such as natural gas, hydrogen and water networks, flows arise from nonlinear potential differences between the nodes, which can be represented by so-called “potential-driven” network models. When operators of these networks face increasing demand or the need to handle more diverse transport situations, they regularly seek to expand the capacity of their network by building new pipelines parallel to existing ones (“looping”). The paper introduces a new mixed-integer non-linear programming (MINLP) model and a new non-linear programming (NLP) model and compares these with existing models for the looping problem and related problems in the literature, both theoretically and experimentally. On this basis, we give recommendations about the circumstances under which a certain model should be used. In particular, it turns out that one of our novel models outperforms the existing models.
Moreover, the paper is the first to include the practically relevant option that a particular pipeline may be looped several times.
We discuss drift-diffusion models for charge-carrier transport in organic semiconductor de- vices. The crucial feature in organic materials is the energetic disorder due to random alignment of molecules and the hopping transport of carriers between adjacent energetic sites. The for- mer leads to so-called Gauss-Fermi statistics, which describe the occupation of energy levels by electrons and holes. The latter gives rise to complicated mobility models with a strongly nonlin- ear dependence on temperature, density of carriers, and electric field strength. We present the state-of-the-art modeling of the transport processes and provide a first existence result for the stationary drift-diffusion model taking all of the peculiarities of organic materials into account. The existence proof is based on Schauder’s fixed-point theorem. Finally, we discuss the numerical discretization of the model using finite-volume methods and a generalized Scharfetter-Gummel scheme for the Gauss-Fermi statistics.
This paper deals with the analysis of an instationary drift-diffusion model for organic semicon- ductor devices including Gauss–Fermi statistics and application-specific mobility functions. The charge transport in organic materials is realized by hopping of carriers between adjacent ener- getic sites and is described by complicated mobility laws with a strong nonlinear dependence on temperature, carrier densities and the electric field strength.
To prove the existence of global weak solutions, we consider a problem with (for small den- sities) regularized state equations on any arbitrarily chosen finite time interval. We ensure its solvability by time discretization and passage to the time-continuous limit. Positive lower a priori estimates for the densities of its solutions that are independent of the regularization level en- sure the existence of solutions to the original problem. Furthermore, we derive for these solutions global positive lower and upper bounds strictly below the density of transport states for the densi- ties. The estimates rely on Moser iteration techniques.
This article deals with the spectra of Laplacians of weighted graphs. In this context, two objects are of fundamental importance for the dynamics of complex networks: the second eigenvalue of such a spectrum (called algebraic connectivity) and its associated eigenvector, the so-called Fiedler vector. Here we prove that, given a Laplacian matrix, it is possible to perturb the weights of the existing edges in the underlying graph in order to obtain simple eigenvalues and a Fiedler vector composed of only non-zero entries. These structural genericity properties with the constraint of not adding edges in the underlying
graph are stronger than the classical ones, for which arbitrary structural perturbations are allowed. These results open the opportunity to understand the impact of structural changes on the dynamics of complex systems.
An optimal control problem is studied for a quasilinear Maxwell equation of nondegenerate parabolic
type. Well-posedness of the quasilinear state equation, existence of an optimal control, and weak G\^ateaux-differentiability of the control-to-state mapping are proved. Based on these results, first-order necessary optimality conditions and
an associated adjoint calculus are derived.
A mathematical model is set up that can be useful for controlled voltage excitation in time-dependent electromagnetism.
The well-posedness of the model is proved and an associated optimal control problem is investigated. Here, the control
function is a transient voltage and the aim of the control is the best approximation of desired electric and magnetic fields in
suitable $L^2$-norms.
Special emphasis is laid on an adjoint calculus for first-order necessary optimality conditions.
Moreover, a {peculiar attention is devoted to propose a formulation for which the computational complexity of the finite element solution method is substantially reduced}.
We report the cancellation of the soliton self-frequency shift in nonlinear optical fibers. A soliton which interacts with a group velocity matched low intensity dispersive pump pulse, experiences a continuous blue-shift in frequency, which counteracts the soliton self- frequency shift due to Raman scattering. The soliton self-frequency shift can be fully compensated by a suitably prepared dispersive wave. We quantify this kind of soliton-dispersive wave interaction by an adiabatic approach and demonstrate that the compensation is stable in agreement with numerical simulations.
We consider the phenomenon of an optical soliton controlled (e.g. amplified) by a much weaker second pulse which is efficiently scattered at the soliton. An important problem in this context is to quantify the small range of parameters at which the interaction takes place. This has been achieved by using adiabatic ODEs for the soliton characteristics, which is much faster than an empirical scan of the full propagation equations for all parameters in question.
We consider scattering of small-amplitude dispersive waves at an intense optical soliton which constitutes a nonlinear perturbation of the refractive index. Specifically, we consider a single-mode optical fiber and a group velocity matched pair: an optical soliton and a nearly perfectly reflected dispersive wave, a fiber-optical analogue of the event horizon. By combining (i) an adiabatic approach that is used in soliton perturbation theory and (ii) scattering theory from Quantum Mechanics, we give a quantitative account for the evolution of all soliton parameters. In particular, we quantify the increase in the soliton peak power that may result in spontaneous appearance of an extremely large, so-called champion soliton. The presented adiabatic theory agrees well with the numerical solutions of the pulse propagation equation. Moreover, for the first time we predict the full frequency band of the scattered dispersive waves and explain an emerging caustic structure in the space-time domain.
Physics of short optical pulses is an important and active research area in nonlinear optics. In what follows we theoretically consider the most extreme representatives of short pulses that contain only several oscillations of electromagnetic field. Description of such pulses is traditionally based on envelope equations and slowly varying envelope approximation, despite the fact that the envelope is not “slow” and, moreover, there is no clear definition of such a “fast” envelope. This happens due to another paradoxical feature: the standard (envelope) generalized nonlinear Schrödinger equation yields very good correspondence to numerical solutions of full Maxwell equations even for few-cycle pulses, a thing that should not be. In what follows we address ultrashort optical pulses using Hamiltonian framework for nonlinear waves. As it appears, the standard optical envelope equation is just a reformulation of general Hamiltonian equations. In a sense, no approximations are required, this is why the generalized nonlinear Schrödinger equation is so effective. Moreover, the Hamiltonian framework greatly contributes to our understanding of “fast” envelope, ultrashort solitons, stability and radiation of optical pulses. Even the inclusion of dissipative terms is possible making the Hamiltonian approach an universal theoretical tool also in extreme nonlinear optics.
Chance constraints represent a popular tool for finding decisions that enforce a
robust satisfaction of random inequality systems in terms of probability. They
are widely used in optimization problems subject to uncertain parameters as they
arise in many engineering applications. Most structural results of chance constraints (e.g., closedness, convexity, Lipschitz continuity, differentiability etc.) have been formulated in a finite-dimensional
setting. The aim of this paper is to generalize some of these well-known semi-continuity and convexity properties to a setting of control problems
subject to (uniform) state chance constraints.
We present a globally convergent method for the solution of frictionless large deformation contact problems involving hyperelastic materials.
For the discretisation we apply the dual mortar method which is known to be more stable than node-to-segment approaches. The resulting non-convex constrained minimisation problems are solved using a filter–trust-region scheme.
This method combines several techniques from non-linear optimisation to achieve global convergence towards first-order optimal points. To speed up the method inexact linearisations of the non-penetration constraint are used
whenever the current iterate is far away from a critical point. A monotone multigrid method is applied for the fast solution of the constrained Newton
problems.
We derive an optoelectronic model based on a gradient formulation for the relaxation of electron-, hole- and photon-densities to their equilibrium state. This leads to a coupled system of partial and ordinary differential equations, for which we discuss the isothermal and the non-isothermal scenario separately.
We show the existence of solutions to a system of elliptic PDEs, that was recently introduced to describe the electrothermal behavior of organic semiconductor devices. Here, two diffculties appear: (i) the elliptic term in the current-flow equation is of p(x)-Laplacian type with discontinuous exponent p, which limits the use of standard methods, and (ii) in the heat equation, we have to deal with an a priori L1 term on the right hand side describing the Joule heating in the device. We prove the existence of a weak solution under very weak assumptions on the data. Our existence proof is based on Schauder’s fixed point theorem and the concept of entropy solutions for the heat equation. Here, the crucial point is the continuous dependence of the entropy solutions on the data of the problem.
Classical gradient systems have a linear relation between rates and driving forces. In generalized gradient systems we allow for arbitrary relations derived from general non-quadratic dissipation potentials. This paper describes two natural origins for these structures.
A first microscopic origin of generalized gradient structures is given by the theory of large-deviation principles. While Markovian diffusion processes lead to classical gradient structures, Poissonian jump processes give rise to cosh-type dissipation potentials.
A second origin arises via a new form of convergence, that we call EDP-convergence. Even when starting with classical gradient systems, where the dissipation potential is a quadratic functional of the rate, we may obtain a generalized gradient system in the evolutionary -limit. As examples we treat (i) the limit of a diffusion equation having a thin layer of low diffusivity, which leads to a membrane model, and (ii) the limit of difusion over a high barrier, which gives a reaction-diffusion system.
Large-area organic light-emitting diodes are thin-film multilayer devices that show pronounced self-heating and brightness inhomogeneities at high currents. As these high currents are typical for lighting applications, a deeper understanding of the mechanisms causing these inhomogeneities is necessary. We discuss the modeling of the interplay between current flow, self-heating, and heat transfer in such devices using a system of partial differential equations of thermistor type, that is capable of explaining the development of luminance inhomogeneities. The system is based on the heat equation for the temperature coupled to a p(x)-Laplace-type equation for the electrostatic potential with mixed boundary conditions. The p(x)-Laplacian allows to take into account non-Ohmic electrical behavior of the different organic layers. Moreover, we present analytical results on the existence, boundedness, and regularity of solutions to the system. A numerical scheme based on the finite-volume method allows for efficient simulations of device structures.
We show the existence of solutions to a system of elliptic PDEs, that was recently introduced to describe the electrothermal behavior of organic semiconductor devices. Here, two difficulties appear: (i) the elliptic term in the current-flow equation is of $p(x)$-Laplacian-type with discontinuous exponent $p$, which limits the use of standard methods, and (ii) in the heat equation, we have to deal with an a priori $L^1$ term on the right hand side describing the Joule heating in the device. We prove the existence of a weak solution under very weak assumptions on the data. Our existence proof is based on Schauder’s fixed point theorem and the concept of entropy solutions for the heat equation. Here, the crucial point is the continuous dependence of the entropy solutions on the data of the problem.
We study linear dissipative Hamiltonian (DH) systems with real constant coefficients that arise in energy based modeling of dynamical
systems. In this paper we analyze when such a system is on the boundary
of the region of asymptotic stability, i.e., when it has purely imaginary eigenvalues,
or how much the dissipation term has to be perturbed to be on this boundary. For unstructured systems the explicit construction of the \emph{real distance to instability (real stability radius)} has been a challenging problem. In this paper, we analyze this real distance under different structured perturbations to the dissipation term that preserve the DH structure and we derive explicit formulas for this distance in terms of low rank perturbations. We also show (via numerical examples) that under real structured perturbations to the dissipation the asymptotical
stability of a DH system is much more robust than for unstructured perturbations.