Search Results
Generalized Nash equilibrium problems with partial differential operators: Theory, algorithms, and risk aversion
PDE-constrained (generalized) Nash equilibrium problems (GNEPs) are considered in a deterministic setting as well as under uncertainty. This includes a study of deterministic GNEPs with nonlinear and/or multivalued operator equations as forward problems and PDE-constrained GNEPs with uncertain data. The deterministic nonlinear problems are analyzed using the theory of generalized convexity for set-valued operators, and a variational approximation approach is proposed. The stochastic setting includes a detailed overview of the recently developed theory and algorithms for risk-averse PDE-constrained optimization problems. These new results open the way to a rigorous study of stochastic PDE-constrained GNEPs.
Strong stationarity conditions for the optimal control of a Cahn--Hilliard--Navier--Stokes system
This paper is concerned with the distributed optimal control of a time-discrete Cahn-Hilliard-Navier-Stokes system with variable densities. It focuses on the double-obstacle potential which yields an optimal control problem for a variational inequality of fourth order and the Navier-Stokes equation. The existence of solutions to the primal system and of optimal controls is established. The Lipschitz continuity of the constraint mapping is derived and used to characterize the directional derivative of the constraint mapping via a system of variational inequalities and partial differential equations. Finally, strong stationarity conditions are presented following an approach from Mignot and Puel.
Optimal control of geometric partial differential equations
Optimal control problems for geometric (evolutionary) partial differential inclusions are considered. The focus is on problems which, in addition to the nonlinearity due to geometric evolution, contain optimization theoretic challenges because of non-smoothness. The latter might stem from energies containing non-smooth constituents such as obstacle-type potentials or terms modeling, e.g., pinning phenomena in microfluidics. Several techniques to remedy the resulting constraint degeneracy when deriving stationarity conditions are presented. A particular focus is on Yosida-type mollifications approximating the original degenerate problem by a sequence of nondegenerate nonconvex optimal control problems. This technique is also the starting point for the development of numerical solution schemes. In this context, also dual-weighted residual based error estimates are addressed to facilitate an adaptive mesh refinement. Concerning the underlying state model, sharp and diffuse interface formulations are discussed. While the former always allows for accurately tracing interfacial motion, the latter model may be dictated by the underlying physical phenomenon, where near the interface mixed phases may exist, but it may also be used as an approximate model for (sharp) interface motion. In view of the latter, (sharp interface) limits of diffuse interface models are addressed. For the sake of presentation, this exposition confines itself to phase field type diffuse interface models and, moreover, develops the optimal control of either of the two interface models along model applications. More precisely, electro-wetting on dielectric is used in the sharp interface context, and the control of multiphase fluids involving spinodal decomposition highlights the phase field technique. Mathematically, the former leads to a Hele-Shaw flow with geometric boundary conditions involving a complementarity system due to contact line pinning, and the latter gives rise to a Cahn-Hilliard Navier-Stokes model including a non-smooth obstacle type potential leading to a variational inequality constraint.
First-order conditions for the optimal control of learning-informed nonsmooth PDEs
In this paper we study the optimal control of a class of semilinear elliptic partial differential equations which have nonlinear constituents that are only accessible by data and are approximated by nonsmooth ReLU neural networks. The optimal control problem is studied in detail. In particular, the existence and uniqueness of the state equation are shown, and continuity as well as directional differentiability properties of the corresponding control-to-state map are established. Based on approximation capabilities of the pertinent networks, we address fundamental questions regarding approximating properties of the learning-informed control-to-state map and the solution of the corresponding optimal control problem. Finally, several stationarity conditions are derived based on different notions of generalized differentiability.
Dualization and automatic distributed parameter selection of total generalized variation via bilevel optimization
Total Generalized Variation (TGV) regularization in image reconstruction relies on an infimal convolution type combination of generalized first- and second-order derivatives. This helps to avoid the staircasing effect of Total Variation (TV) regularization, while still preserving sharp contrasts in images. The associated regularization effect crucially hinges on two parameters whose proper adjustment represents a challenging task. In this work, a bilevel optimization framework with a suitable statistics-based upper level objective is proposed in order to automatically select these parameters. The framework allows for spatially varying parameters, thus enabling better recovery in high-detail image areas. A rigorous dualization framework is established, and for the numerical solution, two Newton type methods for the solution of the lower level problem, i.e. the image reconstruction problem, and two bilevel TGV algorithms are introduced, respectively. Denoising tests confirm that automatically selected distributed regularization parameters lead in general to improved reconstructions when compared to results for scalar parameters.
Uncertainty quantification in image segmentation using the Ambrosio--Tortorelli approximation of the Mumford--Shah energy
The quantification of uncertainties in image segmentation based on the Mumford-Shah model is studied. The aim is to address the error propagation of noise and other error types in the original image to the restoration result and especially the reconstructed edges (sharp image contrasts). Analytically, we rely on the Ambrosio-Tortorelli approximation and discuss the existence of measurable selections of its solutions as well as sampling-based methods and the limitations of other popular methods. Numerical examples illustrate the theoretical findings.
A generalized $Gamma$-convergence concept for a type of equilibrium problems
A novel generalization of Γ-convergence applicable to a class of equilibrium problems is studied. After the introduction of the latter, a variety of its applications is discussed. The existence of equilibria with emphasis on Nash equilibrium problems is investigated. Subsequently, our Γ-convergence notion for equilibrium problems, generalizing the existing one from optimization, is introduced and discussed. The work ends with its application to a class of penalized generalized Nash equilibrium problems and quasi-variational inequalities.
Constrained exact boundary controllability of a semilinear model for pipeline gas flow
While the quasilinear isothermal Euler equations are an excellent model for gas pipeline flow, the operation of the pipeline flow with high pressure and small Mach numbers allows us to obtain approximate solutions by a simpler semilinear model. We provide a derivation of the semilinear model that shows that the semilinear model is valid for sufficiently low Mach numbers and sufficiently high pressures. We prove an existence result for continuous solutions of the semilinear model that takes into account lower and upper bounds for the pressure and an upper bound for the magnitude of the Mach number of the gas flow. These state constraints are important both in the operation of gas pipelines and to guarantee that the solution remains in the set where the model is physically valid. We show the constrained exact boundary controllability of the system with the same pressure and Mach number constraints.
Differentiability properties for boundary control of fluid-structure interactions of linear elasticity with Navier--Stokes equations with mixed-boundary conditions in a channel
In this paper we consider a fluid-structure interaction problem given by the steady Navier Stokes equations coupled with linear elasticity taken from [Lasiecka, Szulc, and Zochoswki, Nonl. Anal.: Real World Appl., 44, 2018]. An elastic body surrounded by a liquid in a rectangular domain is deformed by the flow which can be controlled by the Dirichlet boundary condition at the inlet. On the walls along the channel homogeneous Dirichlet boundary conditions and on the outflow boundary do-nothing conditions are prescribed. We recall existence results for the nonlinear system from that reference and analyze the control to state mapping generaziling the results of [Wollner and Wick, J. Math. Fluid Mech., 21, 2019] to the setting of the nonlinear Navier-Stokes equation for the fluid and the situation of mixed boundary conditions in a domain with corners.
A class of second-order geometric quasilinear hyperbolic PDEs and their application in imaging science
In this paper, we study damped second-order dynamics, which are quasilinear hyperbolic partial differential equations (PDEs). This is inspired by the recent development of second-order damping systems for accelerating energy decay of gradient flows. We concentrate on two equations: one is a damped second-order total variation flow, which is primarily motivated by the application of image denoising; the other is a damped second-order mean curvature flow for level sets of scalar functions, which is related to a non-convex variational model capable of correcting displacement errors in image data (e.g. dejittering). For the former equation, we prove the existence and uniqueness of the solution. For the latter, we draw a connection between the equation and some second-order geometric PDEs evolving the hypersurfaces which are described by level sets of scalar functions, and show the existence and uniqueness of the solution for a regularized version of the equation. The latter is used in our algorithmic development. A general algorithm for numerical discretization of the two nonlinear PDEs is proposed and analyzed. Its efficiency is demonstrated by various numerical examples, where simulations on the behavior of solutions of the new equations and comparisons with first-order flows are also documented.