-
Second order analysis for the optimal selection of time delays
Authors:
Karl Kunisch,
Fredi Troeltzsch
Abstract:
For a nonlinear ordinary differential equation with time delay, the differentiation of the solution with respect to the delay is investigated. Special emphasis is laid on the second-order derivative. The results are applied to an associated optimization problem for the time delay. A first- and second-order sensitivity analysis is performed including an adjoint calculus that avoids the second deriv…
▽ More
For a nonlinear ordinary differential equation with time delay, the differentiation of the solution with respect to the delay is investigated. Special emphasis is laid on the second-order derivative. The results are applied to an associated optimization problem for the time delay. A first- and second-order sensitivity analysis is performed including an adjoint calculus that avoids the second derivative of the state with respect to the delay.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Optimality conditions for sparse optimal control of viscous Cahn-Hilliard systems with logarithmic potential
Authors:
Pierluigi Colli,
Jürgen Sprekels,
Fredi Tröltzsch
Abstract:
In this paper we study the optimal control of a parabolic initial-boundary value problem of viscous Cahn-Hilliard type with zero Neumann boundary conditions. Phase field systems of this type govern the evolution of diffusive phase transition processes with conserved order parameter. It is assumed that the nonlinear functions driving the physical processes within the spatial domain are double-well…
▽ More
In this paper we study the optimal control of a parabolic initial-boundary value problem of viscous Cahn-Hilliard type with zero Neumann boundary conditions. Phase field systems of this type govern the evolution of diffusive phase transition processes with conserved order parameter. It is assumed that the nonlinear functions driving the physical processes within the spatial domain are double-well potentials of logarithmic type whose derivatives become singular at the boundary of their respective domains of definition. For such systems, optimal control problems have been studied in the past. We focus here on the situation when the cost functional of the optimal control problem contains a nondifferentiable term like the L1-norm, which leads to sparsity of optimal controls. For such cases, we establish first-order necessary and second-order sufficient optimality conditions for locally optimal controls. In the approach to second-order sufficient conditions, the main novelty of this paper, we adapt a technique introduced by E. Casas, C. Ryll and F. Tröltzsch in the paper [SIAM J. Control Optim. 53 (2015), 2168-2202]. In this paper, we show that this method can also be successfully applied to systems of viscous Cahn-Hilliard type with logarithmic nonlinearity. Since the Cahn-Hilliard system corresponds to a fourth-order partial differential equation in contrast to the second-order systems investigated before, additional technical difficulties have to be overcome.
△ Less
Submitted 19 September, 2024; v1 submitted 28 February, 2024;
originally announced February 2024.
-
Second-order sufficient conditions in the sparse optimal control of a phase field tumor growth model with logarithmic potential
Authors:
Jürgen Sprekels,
Fredi Tröltzsch
Abstract:
This paper treats a distributed optimal control problem for a tumor growth model of viscous Cahn--Hilliard type. The evolution of the tumor fraction is governed by a thermodynamic force induced by a double-well potential of logarithmic type. The cost functional contains a nondifferentiable term like the $L^1$-norm in order to enhance the occurrence of sparsity effects in the optimal control, i.e.,…
▽ More
This paper treats a distributed optimal control problem for a tumor growth model of viscous Cahn--Hilliard type. The evolution of the tumor fraction is governed by a thermodynamic force induced by a double-well potential of logarithmic type. The cost functional contains a nondifferentiable term like the $L^1$-norm in order to enhance the occurrence of sparsity effects in the optimal control, i.e., of subdomains of the space-time cylinder where the controls vanish. In the context of cancer therapies, sparsity is very important in order that the patient is not exposed to unnecessary intensive medical treatment. In this work, we focus on the derivation of second-order sufficient optimality conditions for the optimal control problem. While in previous works on the system under investigation such conditions have been established for the case without sparsity, the case with sparsity has not been treated before.
△ Less
Submitted 10 June, 2023;
originally announced June 2023.
-
Second-order sufficient conditions for sparse optimal control of singular Allen--Cahn systems with dynamic boundary conditions
Authors:
Jürgen Sprekels,
Fredi Tröltzsch
Abstract:
In this paper we study the optimal control of a parabolic initial-boundary value problem of Allen--Cahn type with dynamic boundary conditions. Phase field systems of this type govern the evolution of coupled diffuse phase transition processes with nonconserved order parameters that occur in a container and on its surface, respectively. It is assumed that the nonlinear function driving the physical…
▽ More
In this paper we study the optimal control of a parabolic initial-boundary value problem of Allen--Cahn type with dynamic boundary conditions. Phase field systems of this type govern the evolution of coupled diffuse phase transition processes with nonconserved order parameters that occur in a container and on its surface, respectively. It is assumed that the nonlinear function driving the physical processes within the bulk and on the surface are double well potentials of logarithmic type whose derivatives become singular at the boundary of their respective domains of definition. For such systems, optimal control problems have been studied in the past. We focus here on the situation when the cost functional of the optimal control problem contains a nondifferentiable term like the $L^1$-norm leading to sparsity of optimal controls. For such cases, we derive second-order sufficient conditions for locally optimal controls.
△ Less
Submitted 29 March, 2023;
originally announced March 2023.
-
Space-time finite element methods for the initial temperature reconstruction
Authors:
Ulrich Langer,
Olaf Steinbach,
Fredi Tröltzsch,
Huidong Yang
Abstract:
This work is devoted to the reconstruction of the initial temperature in the backward heat equation using the space-time finite element method on fully unstructured space-time simplicial meshes proposed by Steinbach (2015). Such a severely ill-posed problem is tackled by the standard Tikhonov regularization method. This leads to a related optimal control for an parabolic equation in the space-time…
▽ More
This work is devoted to the reconstruction of the initial temperature in the backward heat equation using the space-time finite element method on fully unstructured space-time simplicial meshes proposed by Steinbach (2015). Such a severely ill-posed problem is tackled by the standard Tikhonov regularization method. This leads to a related optimal control for an parabolic equation in the space-time domain. In this setting, the control is taken as initial condition, whereas the terminal observation data serve as target. The objective becomes a standard terminal observation functional combined with the Tikhonov regularization. The space-time finite element method is applied to the space-time optimality system that is well-posed for a fixed regularization parameter.
△ Less
Submitted 30 March, 2021;
originally announced March 2021.
-
Sparse optimal control of a phase field system with singular potentials arising in the modeling of tumor growth
Authors:
Jürgen Sprekels,
Fredi Tröltzsch
Abstract:
In this paper, we study an optimal control problem for a nonlinear system of reaction-diffusion equations that constitutes a simplified and relaxed version of a thermodynamically consistent phase field model for tumor growth originally introduced in [12]. The model takes the effect of chemotaxis into account but neglects velocity contributions. The unknown quantities of the governing state equatio…
▽ More
In this paper, we study an optimal control problem for a nonlinear system of reaction-diffusion equations that constitutes a simplified and relaxed version of a thermodynamically consistent phase field model for tumor growth originally introduced in [12]. The model takes the effect of chemotaxis into account but neglects velocity contributions. The unknown quantities of the governing state equations are the chemical potential, the (normalized) tumor fraction, and the nutrient extra-cellular water concentration. The equation governing the evolution of the tumor fraction is dominated by the variational derivative of a double-well potential which may be of singular (e.g., logarithmic) type. In contrast to the recent paper [10] on the same system, we consider in this paper sparsity effects, which means that the cost functional contains a nondifferentiable (but convex) contribution like the $L^1-$norm. For such problems, we derive first-order necessary optimality conditions and conditions for directional sparsity, both with respect to space and time, where the latter case is of particular interest for practical medical applications in which the control variables are given by the administration of cytotoxic drugs or by the supply of nutrients. In addition to these results, we prove that the corresponding control-to-state operator is twice continuously differentiable between suitable Banach spaces, using the implicit function theorem. This result, which complements and sharpens a differentiability result derived in [10], constitutes a prerequisite for a future derivation of second-order sufficient optimality conditions.
△ Less
Submitted 22 August, 2020; v1 submitted 6 May, 2020;
originally announced May 2020.
-
Space-time finite element discretization of parabolic optimal control problems with energy regularization
Authors:
Ulrich Langer,
Olaf Steinbach,
Fredi Tröltzsch,
Huidong Yang
Abstract:
We analyze space-time finite element methods for the numerical solution of distributed parabolic optimal control problems with energy regularization in the Bochner space $L^2(0,T;H^{-1}(Ω))$. By duality, the related norm can be evaluated by means of the solution of an elliptic quasi-stationary boundary value problem. When eliminating the control, we end up with the reduced optimality system that i…
▽ More
We analyze space-time finite element methods for the numerical solution of distributed parabolic optimal control problems with energy regularization in the Bochner space $L^2(0,T;H^{-1}(Ω))$. By duality, the related norm can be evaluated by means of the solution of an elliptic quasi-stationary boundary value problem. When eliminating the control, we end up with the reduced optimality system that is nothing but the variational formulation of the coupled forward-backward primal and adjoint equations. Using Babuška's theorem, we prove unique solvability in the continuous case. Furthermore, we establish the discrete inf-sup condition for any conforming space-time finite element discretization yielding quasi-optimal discretization error estimates. Various numerical examples confirm the theoretical findings. We emphasize that the energy regularization results in a more localized control with sharper contours for discontinuous target functions, which is demonstrated by a comparison with an $L^2$ regularization and with a sparse optimal control approach.
△ Less
Submitted 20 April, 2020;
originally announced April 2020.
-
Unstructured space-time finite element methods for optimal control of parabolic equations
Authors:
Ulrich Langer,
Olaf Steinbach,
Fredi Tröltzsch,
Huidong Yang
Abstract:
This work presents and analyzes space-time finite element methods on fully unstructured simplicial space-time meshes for the numerical solution of parabolic optimal control problems. Using Babuška's theorem, we show well-posedness of the first-order optimality systems for a typical model problem with linear state equations, but without control constraints. This is done for both continuous and disc…
▽ More
This work presents and analyzes space-time finite element methods on fully unstructured simplicial space-time meshes for the numerical solution of parabolic optimal control problems. Using Babuška's theorem, we show well-posedness of the first-order optimality systems for a typical model problem with linear state equations, but without control constraints. This is done for both continuous and discrete levels. Based on these results, we derive discretization error estimates. Then we consider a semilinear parabolic optimal control problem arising from the Schlögl model. The associated nonlinear optimality system is solved by Newton's method, where a linear system, that is similar to the first-order optimality systems considered for the linear model problems, has to be solved at each Newton step. We present various numerical experiments including results for adaptive space-time finite element discretizations based on residual-type error indicators. In the last two examples, we also consider semilinear parabolic optimal control problems with box constraints imposed on the control.
△ Less
Submitted 4 April, 2020;
originally announced April 2020.
-
Unstructured Space-Time Finite Element Methods for Optimal Sparse Control of Parabolic Equations
Authors:
Ulrich Langer,
Olaf Steinbach,
Fredi Tröltzsch,
Huidong Yang
Abstract:
We consider a space-time finite element method on fully unstructured simplicial meshes for optimal sparse control of semilinear parabolic equations. The objective is a combination of a standard quadratic tracking-type functional including a Tikhonov regularization term and of the $L^1$-norm of the control that accounts for its spatio-temporal sparsity. We use a space-time Petrov-Galerkin finite el…
▽ More
We consider a space-time finite element method on fully unstructured simplicial meshes for optimal sparse control of semilinear parabolic equations. The objective is a combination of a standard quadratic tracking-type functional including a Tikhonov regularization term and of the $L^1$-norm of the control that accounts for its spatio-temporal sparsity. We use a space-time Petrov-Galerkin finite element discretization for the first-order necessary optimality system of the associated discrete optimal sparse control problem. The discretization is based on a variational formulation that employs piecewise linear finite elements simultaneously in space and time. Finally, the discrete nonlinear optimality system that consists of coupled forward-backward state and adjoint state equations is solved by a semismooth Newton method.
△ Less
Submitted 31 March, 2020;
originally announced March 2020.
-
First and Second Order Optimality Conditions for the Control of Fokker-Planck Equations
Authors:
M. Soledad Aronna,
Fredi Tröltzsch
Abstract:
In this article we study an optimal control problem subject to the Fokker-Planck equation \[ \partial_t ρ- νΔρ- {\rm div } \big(ρB[u]\big) = 0. \] The control variable $u$ is time-dependent and possibly multidimensional, and the function $B$ depends on the space variable and the control. The cost functional is of tracking type and includes a quadratic regularization term on the control. For this p…
▽ More
In this article we study an optimal control problem subject to the Fokker-Planck equation \[ \partial_t ρ- νΔρ- {\rm div } \big(ρB[u]\big) = 0. \] The control variable $u$ is time-dependent and possibly multidimensional, and the function $B$ depends on the space variable and the control. The cost functional is of tracking type and includes a quadratic regularization term on the control. For this problem, we prove existence of optimal controls and first order necessary conditions. Main emphasis is placed on second order necessary and sufficient conditions.
△ Less
Submitted 16 January, 2021; v1 submitted 10 February, 2020;
originally announced February 2020.
-
Optimal time delays in a class of reaction-diffusion equations
Authors:
Eduardo Casas,
Mariano Mateos,
Fredi Tröltzsch
Abstract:
A class of semilinear parabolic reaction diffusion equations with multiple time delays is considered. These time delays and corresponding weights are to be optimized such that the associated solution of the delay equation is the best approximation of a desired state function. The differentiability of the mapping is proved that associates the solution of the delay equation to the vector of weights…
▽ More
A class of semilinear parabolic reaction diffusion equations with multiple time delays is considered. These time delays and corresponding weights are to be optimized such that the associated solution of the delay equation is the best approximation of a desired state function. The differentiability of the mapping is proved that associates the solution of the delay equation to the vector of weights and delays. Based on an adjoint calculus, first-order necessary optimality conditions are derived. Numerical test examples show the applicability of the concept of optimizing time delays.
△ Less
Submitted 4 May, 2018;
originally announced May 2018.
-
Measure Control of a Semilinear Parabolic Equation with a Nonlocal Time Delay
Authors:
Eduardo Casas,
Mariano Mateos,
Fredi Tröltzsch
Abstract:
We study a control problem governed by a semilinear parabolic equation. The control is a measure that acts as the kernel of a possibly nonlocal time delay term and the functional includes a non-differentiable term with the measure-norm of the control. Existence, uniqueness and regularity of the solution of the state equation, as well as differentiability properties of the control-to-state operator…
▽ More
We study a control problem governed by a semilinear parabolic equation. The control is a measure that acts as the kernel of a possibly nonlocal time delay term and the functional includes a non-differentiable term with the measure-norm of the control. Existence, uniqueness and regularity of the solution of the state equation, as well as differentiability properties of the control-to-state operator are obtained. Next, we provide first order optimality conditions for local solutions. Finally, the control space is suitably discretized and we prove convergence of the solutions of the discrete problems to the solutions of the original problem. Several numerical examples are included to illustrate the theoretical results.
△ Less
Submitted 2 May, 2018;
originally announced May 2018.
-
On the switching behavior of sparse optimal controls for the one-dimensional heat equation
Authors:
Fredi Troeltzsch,
Daniel Wachsmuth
Abstract:
An optimal boundary control problem for the one-dimensional heat equation is considered. The objective functional includes a standard quadratic terminal observation, a Tikhonov regularization term with regularization parameter $ν$, and the $L^1$-norm of the control that accounts for sparsity. The switching structure of the optimal control is discussed for $ν\ge 0$. Under natural assumptions, it is…
▽ More
An optimal boundary control problem for the one-dimensional heat equation is considered. The objective functional includes a standard quadratic terminal observation, a Tikhonov regularization term with regularization parameter $ν$, and the $L^1$-norm of the control that accounts for sparsity. The switching structure of the optimal control is discussed for $ν\ge 0$. Under natural assumptions, it is shown that the set of switching points of the optimal control is countable with the final time as only possible accumulation point. The convergence of switching points is investigated for $ν\searrow 0$.
△ Less
Submitted 10 May, 2017; v1 submitted 9 May, 2017;
originally announced May 2017.
-
Control of traveling localized spots
Authors:
Steffen Martens,
Christopher Ryll,
Jakob Löber,
Fredi Tröltzsch,
Harald Engel
Abstract:
Traveling localized spots represent an important class of self-organized two-dimensional patterns in reaction-diffusion systems. We study open-loop control intended to guide a stable spot along a desired trajectory with desired velocity. Simultaneously, the spot's concentration profile does not change under control. For a given protocol of motion, we first express the control signal analytically i…
▽ More
Traveling localized spots represent an important class of self-organized two-dimensional patterns in reaction-diffusion systems. We study open-loop control intended to guide a stable spot along a desired trajectory with desired velocity. Simultaneously, the spot's concentration profile does not change under control. For a given protocol of motion, we first express the control signal analytically in terms of the Goldstone modes and the propagation velocity of the uncontrolled spot. Thus, detailed information about the underlying nonlinear reaction kinetics is unnecessary. Then, we confirm the optimality of this solution by demonstrating numerically its equivalence to the solution of a regularized, optimal control problem. To solve the latter, the analytical expressions for the control are excellent initial guesses speeding-up substantially the otherwise time-consuming calculations.
△ Less
Submitted 20 September, 2018; v1 submitted 13 March, 2017;
originally announced March 2017.
-
Analytical, Optimal, and Sparse Optimal Control of Traveling Wave Solutions to Reaction-Diffusion Systems
Authors:
Christopher Ryll,
Jakob Löber,
Steffen Martens,
Harald Engel,
Fredi Tröltzsch
Abstract:
This work deals with the position control of selected patterns in reaction-diffusion systems. Exemplarily, the Schlögl and FitzHugh-Nagumo model are discussed using three different approaches. First, an analytical solution is proposed. Second, the standard optimal control procedure is applied. The third approach extends standard optimal control to so-called sparse optimal control that results in v…
▽ More
This work deals with the position control of selected patterns in reaction-diffusion systems. Exemplarily, the Schlögl and FitzHugh-Nagumo model are discussed using three different approaches. First, an analytical solution is proposed. Second, the standard optimal control procedure is applied. The third approach extends standard optimal control to so-called sparse optimal control that results in very localized control signals and allows the analysis of second order optimality conditions.
△ Less
Submitted 22 April, 2016; v1 submitted 29 July, 2015;
originally announced July 2015.
-
Optimization of nonlocal time-delayed feedback controllers
Authors:
Peter Nestler,
Eckehard Schöll,
Fredi Tröltzsch
Abstract:
A class of Pyragas type nonlocal feedback controllers with time-delay is investigated for the Schlögl model. The main goal is to find an optimal kernel in the controller such that the associated solution of the controlled equation is as close as possible to a desired spatio-temporal pattern. An optimal kernel is the solution to a nonlinear optimal control problem with the kernel taken as control f…
▽ More
A class of Pyragas type nonlocal feedback controllers with time-delay is investigated for the Schlögl model. The main goal is to find an optimal kernel in the controller such that the associated solution of the controlled equation is as close as possible to a desired spatio-temporal pattern. An optimal kernel is the solution to a nonlinear optimal control problem with the kernel taken as control function. The well-posedness of the optimal control problem and necessary optimality conditions are discussed for different types of kernels. Special emphasis is laid on time-periodic functions as desired patterns. Here, the cross correlation between the state and the desired pattern is invoked to set up an associated objective functional that is to be minimized. Numerical examples are presented for the 1D Schlögl model and a class of simple step functions for the kernel.
△ Less
Submitted 8 May, 2015;
originally announced May 2015.