Search Results

Now showing 1 - 10 of 22
  • Item
    A boundary control problem for the pure Cahn–Hilliard equation with dynamic boundary conditions
    (Berlin ; Boston, Mass. : de Gruyter, 2015) Colli, Pierluigi; Gilardi, Gianni; Sprekels, Jürgen
    A boundary control problem for the pure Cahn–Hilliard equations with possibly singular potentialsand dynamic boundary conditions is studied and rst-order necessary conditions for optimality are proved.
  • Item
    Sparse optimal control of a phase field system with singular potentials arising in the modeling of tumor growth
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Sprekels, Jürgen; Tröltzsch, Fredi
    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 [13]. 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 L1-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.
  • Item
    Optimal control problems with sparsity for phase field tumor growth models involving variational inequalities
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2021) Colli, Pierluigi; Signori, Andrea; Sprekels, Jürgen
    This paper treats a distributed optimal control problem for a tumor growth model of Cahn--Hilliard type including chemotaxis. The evolution of the tumor fraction is governed by a variational inequality corresponding to a double obstacle nonlinearity occurring in the associated potential. In addition, the control and state variables are nonlinearly coupled and, furthermore, the cost functional contains a nondifferentiable term like the $L^1$--norm in order to include sparsity effects which is of utmost relevance, especially time sparsity, in the context of cancer therapies as applying a control to the system reflects in exposing the patient to an intensive medical treatment. To cope with the difficulties originating from the variational inequality in the state system, we employ the so-called ``deep quench approximation'' in which the convex part of the double obstacle potential is approximated by logarithmic functions. For such functions, first-order necessary conditions of optimality can be established by invoking recent results. We use these results to derive corresponding optimality conditions also for the double obstacle case, by deducing a variational inequality in terms of the associated adjoint state variables. The resulting variational inequality can be exploited to also obtain sparsity results for the optimal controls.
  • Item
    Solving joint chance constrained problems using regularization and Benders decomposition
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Adam, Lukás; Branda, Martin; Heitsch, Holger; Henrion, René
    In this paper we investigate stochastic programms with joint chance constraints. We consider discrete scenario set and reformulate the problem by adding auxiliary variables. Since the resulting problem has a difficult feasible set, we regularize it. To decrease the dependence on the scenario number, we propose a numerical method by iteratively solving a master problem while adding Benders cuts. We find the solution of the slave problem (generating the Benders cuts) in a closed form and propose a heuristic method to decrease the number of cuts. We perform a numerical study by increasing the number of scenarios and compare our solution with a solution obtained by solving the same problem with continuous distribution.
  • Item
    On M-stationarity conditions in MPECs and the associated qualification conditions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Adam, Lukáš; Henrion, René; Outrata, Jir̆í
    Depending on whether a mathematical program with equilibrium constraints (MPEC) is considered in its original or its enhanced (via KKT conditions) form, the assumed constraint qualifications (CQs) as well as the derived necessary optimality conditions may differ significantly. In this paper, we study this issue when imposing one of the weakest possible CQs, namely the calmness of the perturbation mapping associated with the respective generalized equations in both forms of the MPEC. It is well known that the calmness property allows one to derive socalled M-stationarity conditions. The strength of assumptions and conclusions in the two forms of the MPEC is strongly related with the CQs on the lower level imposed on the set whose normal cone appears in the generalized equation. For instance, under just the Mangasarian-Fromovitz CQ (a minimum assumption required for this set), the calmness properties of the original and the enhanced perturbation mapping are drastically different. They become identical in the case of a polyhedral set or when adding the Full Rank CQ. On the other hand, the resulting optimality conditions are affected too. If the considered set even satisfies the Linear Independence CQ, both the calmness assumption and the derived optimality conditions are fully equivalent for the original and the enhanced form of the MPEC. A compilation of practically relevant consequences of our analysis in the derivation of necessary optimality conditions is provided in the main Theorem 4.3. The obtained results are finally applied to MPECs with structured equilibria.
  • Item
    Optimal boundary control of a viscous Cahn-Hilliard system with dynamic boundary condition and double obstacle potentials
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Colli, Pierluigi; Farshbaf Shaker, Mohammad Hassan; Gilardi, Gianni; Sprekels, Jürgen
    In this paper, we investigate optimal boundary control problems for Cahn--Hilliard variational inequalities with a dynamic boundary condition involving double obstacle potentials and the Laplace--Beltrami operator. The cost functional is of standard tracking type, and box constraints for the controls are prescribed. We prove existence of optimal controls and derive first-order necessary conditions of optimality. The general strategy, which follows the lines of the recent approach by Colli, Farshbaf-Shaker, Sprekels (see Appl. Math. Optim., 2014) to the (simpler) Allen--Cahn case, is the following: we use the results that were recently established by Colli, Gilardi, Sprekels in the preprint arXiv:1407.3916 [math.AP] for the case of (differentiable) logarithmic potentials and perform a so-called ``deep quench limit''. Using compactness and monotonicity arguments, it is shown that this strategy leads to the desired first-order necessary optimality conditions for the case of (non-differentiable) double obstacle potentials.
  • Item
    A boundary control problem for the pure Cahn-Hilliard equation with dynamic boundary conditions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Colli, Pierluigi; Gilardi, Gianni; Sprekels, Jürgen
    A boundary control problem for the pure Cahn-Hilliard equations with possibly singular potentials and dynamic boundary conditions is studied and first-order necessary conditions for optimality are proved.
  • Item
    A boundary control problem for the viscous Cahn-Hilliard equation with dynamic boundary conditions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Colli, Pierluigi; Gilardi, Gianni; Sprekels, Jürgen
    A boundary control problem for the viscous Cahn-Hilliard equations with possibly singular potentials and dynamic boundary conditions is studied and first order necessary conditions for optimality are proved.
  • Item
    An enumerative formula for the spherical cap discrepancy
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Heitsch, Holger; Henrion, René
    The spherical cap discrepancy is a widely used measure for how uniformly a sample of points on the sphere is distributed. Being hard to compute, this discrepancy measure is typically replaced by some lower or upper estimates when designing optimal sampling schemes for the uniform distribution on the sphere. In this paper, we provide a fully explicit, easy to implement enumerative formula for the spherical cap discrepancy. Not surprisingly, this formula is of combinatorial nature and, thus, its application is limited to spheres of small dimension and moderate sample sizes. Nonetheless, it may serve as a useful calibrating tool for testing the efficiency of sampling schemes and its explicit character might be useful also to establish necessary optimality conditions when minimizing the discrepancy with respect to a sample of given size.
  • Item
    Topology optimization subject to additive manufacturing constraints
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Ebeling-Rump, Moritz; Hömberg, Dietmar; Lasarzik, Robert; Petzold, Thomas
    In Topology Optimization the goal is to find the ideal material distribution in a domain subject to external forces. The structure is optimal if it has the highest possible stiffness. A volume constraint ensures filigree structures, which are regulated via a Ginzburg-Landau term. During 3D Printing overhangs lead to instabilities, which have only been tackled unsatisfactorily. The novel idea is to incorporate an Additive Manufacturing Constraint into the phase field method. A rigorous analysis proves the existence of a solution and leads to first order necessary optimality conditions. With an Allen-Cahn interface propagation the optimization problem is solved iteratively. At a low computational cost the Additive Manufacturing Constraint brings about support structures, which can be fine tuned according to engineering demands. Stability during 3D Printing is assured, which solves a common Additive Manufacturing problem.