Source author record

Karl Kunisch

Karl Kunisch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Shape Derivative for Penalty-Constrained Nonsmooth-Nonconvex Optimization: Cohesive Crack Problem

A class of non-smooth and non-convex optimization problems with penalty constraints linked to variational inequalities (VI) is studied with respect to its shape differentiability. The specific problem stemming from quasi-brittle fracture describes an elastic body with a Barenblatt cohesive crack under the inequality condition of non-penetration at the crack faces. Based on the Lagrange approach and using smooth penalization with the Lavrentiev regularization, a formula for the shape derivative is derived. The explicit formula contains both primal and adjoint states and is useful for finding descent directions for a gradient algorithm to identify an optimal crack shape from a boundary measurement. Numerical examples of destructive testing are presented in 2D.

preprint2020arXiv

An inverse problem involving a viscous Eikonal equation with applications in electrophysiology

In this work we discuss the reconstruction of cardiac activation instants based on a viscous Eikonal equation from boundary observations. The problem is formulated as an least squares problem and solved by a projected version of the Levenberg Marquardt method. Moreover, we analyze the wellposeness of the state equation and derive the gradient of the least squares functional with respect to the activation instants. In the numerical examples we also conduct an experiment in which the location of the activation sites and the activation instants are reconstructed jointly based on an adapted version of the shape gradient method from https://link.springer.com/article/10.1007/s00285-019-01419-3. We are able to reconstruct the activation instants as well as the locations of the activations with high accuracy relative to the noise level.

preprint2020arXiv

Neural Network Based Nonlinear Observers

Nonlinear observers based on the well-known concept of minimum energy estimation are discussed. The approach relies on an output injection operator determined by a Hamilton-Jacobi-Bellman equation and is subsequently approximated by a neural network. A suitable optimization problem allowing to learn the network parameters is proposed and numerically investigated for linear and nonlinear oscillators.

preprint2020arXiv

Optimal Control of the 2D Evolutionary Navier-Stokes Equations with Measure Valued Controls

In this paper, we consider an optimal control problem for the two-dimensional evolutionary Navier-Stokes system. Looking for sparsity, we take controls as functions of time taking values in a space of Borel measures. The cost functional does not involve directly the control but we assume some constraints on them. We prove the well-posedness of the control problem and derive necessary and sufficient conditions for local optimality of the controls.

preprint2020arXiv

Policy iteration for Hamilton-Jacobi-Bellman equations with control constraints

Policy iteration is a widely used technique to solve the Hamilton Jacobi Bellman (HJB) equation, which arises from nonlinear optimal feedback control theory. Its convergence analysis has attracted much attention in the unconstrained case. Here we analyze the case with control constraints both for the HJB equations which arise in deterministic and in stochastic control cases. The linear equations in each iteration step are solved by an implicit upwind scheme. Numerical examples are conducted to solve the HJB equation with control constraints and comparisons are shown with the unconstrained cases.

preprint2020arXiv

Semiglobal optimal feedback stabilization of autonomous systems via deep neural network approximation

A learning approach for optimal feedback gains for nonlinear continuous time control systems is proposed and analysed. The goal is to establish a rigorous framework for computing approximating optimal feedback gains using neural networks. The approach rests on two main ingredients. First, an optimal control formulation involving an ensemble of trajectories with 'control' variables given by the feedback gain functions. Second, an approximation to the feedback functions via realizations of neural networks. Based on universal approximation properties we prove the existence and convergence of optimal stabilizing neural network feedback controllers.

preprint2020arXiv

Total Deep Variation for Linear Inverse Problems

Diverse inverse problems in imaging can be cast as variational problems composed of a task-specific data fidelity term and a regularization term. In this paper, we propose a novel learnable general-purpose regularizer exploiting recent architectural design patterns from deep learning. We cast the learning problem as a discrete sampled optimal control problem, for which we derive the adjoint state equations and an optimality condition. By exploiting the variational structure of our approach, we perform a sensitivity analysis with respect to the learned parameters obtained from different training datasets. Moreover, we carry out a nonlinear eigenfunction analysis, which reveals interesting properties of the learned regularizer. We show state-of-the-art performance for classical image restoration and medical image reconstruction problems.

preprint2020arXiv

Total Deep Variation: A Stable Regularizer for Inverse Problems

Various problems in computer vision and medical imaging can be cast as inverse problems. A frequent method for solving inverse problems is the variational approach, which amounts to minimizing an energy composed of a data fidelity term and a regularizer. Classically, handcrafted regularizers are used, which are commonly outperformed by state-of-the-art deep learning approaches. In this work, we combine the variational formulation of inverse problems with deep learning by introducing the data-driven general-purpose total deep variation regularizer. In its core, a convolutional neural network extracts local features on multiple scales and in successive blocks. This combination allows for a rigorous mathematical analysis including an optimal control formulation of the training problem in a mean-field setting and a stability analysis with respect to the initial values and the parameters of the regularizer. In addition, we experimentally verify the robustness against adversarial attacks and numerically derive upper bounds for the generalization error. Finally, we achieve state-of-the-art results for numerous imaging tasks.

preprint2016arXiv

Hybrid Optimal Control Problems for a Class of Semilinear Parabolic Equations

A class of optimal control problems of hybrid nature governed by semilinear parabolic equations is considered. These problems involve the optimization of switching times at which the dynamics, the integral cost, and the bounds on the control may change. First- and second-order optimality conditions are derived. The analysis is based on a reformulation involving a judiciously chosen transformation of the time domains. For autonomous systems and time-independent integral cost, we prove that the Hamiltonian is constant in time when evaluated along the optimal controls and trajectories. A numerical example is provided.

preprint2016arXiv

Sparse Control of Quantum Systems

A new class of cost functionals for optimal control of quantum systems which produces controls which are sparse in frequency and smooth in time is proposed. This is achieved by penalizing a suitable time-frequency representation of the control field, rather than the control field itself, and by employing norms which are of $L^1$ or measure form with respect to frequency but smooth with respect to time. We prove existence of optimal controls for the resulting nonsmooth optimization problem, derive necessary optimality conditions, and rigorously establish the frequency-sparsity of the optimizers. More precisely, we show that the time-frequency representation of the control field, which a priori admits a continuum of frequencies, is supported on only \textit{ finitely many} frequencies. These results cover important systems of physical interest, including (infinite-dimensional) Schrödinger dynamics on multiple potential energy surfaces as arising in laser control of chemical reactions. Numerical simulations confirm that the optimal controls, unlike those obtained with the usual $L^2$ costs, concentrate on just a few frequencies, even in the infinite-dimensional case of laser-controlled chemical reactions.

preprint2015arXiv

Local minimization algorithms for dynamic programming equations

The numerical realization of the dynamic programming principle for continuous-time optimal control leads to nonlinear Hamilton-Jacobi-Bellman equations which require the minimization of a nonlinear mapping over the set of admissible controls. This minimization is often performed by comparison over a finite number of elements of the control set. In this paper we demonstrate the importance of an accurate realization of these minimization problems and propose algorithms by which this can be achieved effectively. The considered class of equations includes nonsmooth control problems with $\ell_1$-penalization which lead to sparse controls.

preprint2015arXiv

Optimal control of the bidomain system (IV): Corrected proofs of the stability and regularity theorems

In a series of papers on optimal control problems for the monodomain as well as for the bidomain equations of cardiac electrophysiology, the authors studied existence of minimizers and derived first-order necessary optimality conditions. The analysis of these control problems is based on a regularity discussion for weak solutions, resulting in a stability estimate and a uniqueness theorem for the monodomain and bidomain system, respectively. Unfortunately, the authors recognized a serious error within the proof of these theorems. However, the present investigation shows that the assertions from [Kunisch/Wagner 12] and [Kunisch/Wagner 13a] can be maintained (with minor changes only) while the proofs must be subjected to considerable alterations. As a consequence, the optimization theorems from [Kunisch/Wagner 13b] allow for substantial improvements. Therefore, in the present paper we provide a refined regularity discussion of the bidomain system together with corrected proofs.