777 research outputs found
Total variation regularization of multi-material topology optimization
This work is concerned with the determination of the diffusion coefficient
from distributed data of the state. This problem is related to homogenization
theory on the one hand and to regularization theory on the other hand. An
approach is proposed which involves total variation regularization combined
with a suitably chosen cost functional that promotes the diffusion coefficient
assuming prespecified values at each point of the domain. The main difficulty
lies in the delicate functional-analytic structure of the resulting
nondifferentiable optimization problem with pointwise constraints for functions
of bounded variation, which makes the derivation of useful pointwise optimality
conditions challenging. To cope with this difficulty, a novel reparametrization
technique is introduced. Numerical examples using a regularized semismooth
Newton method illustrate the structure of the obtained diffusion coefficient.
Space-Time Mixed System Formulation of Phase-Field Fracture Optimal Control Problems
In this work, space-time formulations and Galerkin discretizations for phase-field fracture optimal control problems are considered. The fracture irreversibility constraint is formulated on the time-continuous level and is regularized by means of penalization. The optimization scheme is formulated in terms of the reduced approach and then solved with a Newton method. To this end, the state, adjoint, tangent, and adjoint Hessian equations are derived. The key focus is on the design of appropriate function spaces and the rigorous justification of all Fréchet derivatives that require fourth-order regularizations. Therein, a second-order time derivative on the phase-field variable appears, which is reformulated as a mixed first-order-in-time system. These derivations are carefully established for all four equations. Finally, the corresponding time-stepping schemes are derived by employing a dG(r) discretization in time
Simulation and control of a nonsmooth Cahn--Hilliard Navier--Stokes system with variable fluid densities
We are concerned with the simulation and control of a two phase flow model governed by a coupled Cahn--Hilliard Navier--Stokes system involving a nonsmooth energy potential.We establish the existence of optimal solutions and present two distinct approaches to derive suitable stationarity conditions for the bilevel problem, namely C- and strong stationarity. Moreover, we demonstrate the numerical realization of these concepts at the hands of two adaptive solution algorithms relying on a specifically developed goal-oriented error estimator.In addition, we present a model order reduction approach using proper orthogonal decomposition (POD-MOR) in order to replace high-fidelity models by low order surrogates. In particular, we combine POD with space-adapted snapshots and address the challenges which are the consideration of snapshots with different spatial resolutions and the conservation of a solenoidal property
International Conference on Continuous Optimization (ICCOPT) 2019 Conference Book
The Sixth International Conference on Continuous Optimization took place on the campus of the Technical University of Berlin, August 3-8, 2019. The ICCOPT is a flagship conference of the Mathematical Optimization Society (MOS), organized every three years. ICCOPT 2019 was hosted by the Weierstrass Institute for Applied Analysis and Stochastics (WIAS) Berlin. It included a Summer School and a Conference with a series of plenary and semi-plenary talks, organized and contributed sessions, and poster sessions.
This book comprises the full conference program. It contains, in particular, the scientific program in survey style as well as with all details, and information on the social program, the venue, special meetings, and more
First order algorithms in variational image processing
Variational methods in imaging are nowadays developing towards a quite
universal and flexible tool, allowing for highly successful approaches on tasks
like denoising, deblurring, inpainting, segmentation, super-resolution,
disparity, and optical flow estimation. The overall structure of such
approaches is of the form ; where the functional is a data fidelity term also
depending on some input data and measuring the deviation of from such
and is a regularization functional. Moreover is a (often linear)
forward operator modeling the dependence of data on an underlying image, and
is a positive regularization parameter. While is often
smooth and (strictly) convex, the current practice almost exclusively uses
nonsmooth regularization functionals. The majority of successful techniques is
using nonsmooth and convex functionals like the total variation and
generalizations thereof or -norms of coefficients arising from scalar
products with some frame system. The efficient solution of such variational
problems in imaging demands for appropriate algorithms. Taking into account the
specific structure as a sum of two very different terms to be minimized,
splitting algorithms are a quite canonical choice. Consequently this field has
revived the interest in techniques like operator splittings or augmented
Lagrangians. Here we shall provide an overview of methods currently developed
and recent results as well as some computational studies providing a comparison
of different methods and also illustrating their success in applications.Comment: 60 pages, 33 figure
A Posteriori Error Control for the Binary Mumford-Shah Model
The binary Mumford-Shah model is a widespread tool for image segmentation and
can be considered as a basic model in shape optimization with a broad range of
applications in computer vision, ranging from basic segmentation and labeling
to object reconstruction. This paper presents robust a posteriori error
estimates for a natural error quantity, namely the area of the non properly
segmented region. To this end, a suitable strictly convex and non-constrained
relaxation of the originally non-convex functional is investigated and Repin's
functional approach for a posteriori error estimation is used to control the
numerical error for the relaxed problem in the -norm. In combination with
a suitable cut out argument, a fully practical estimate for the area mismatch
is derived. This estimate is incorporated in an adaptive meshing strategy. Two
different adaptive primal-dual finite element schemes, and the most frequently
used finite difference discretization are investigated and compared. Numerical
experiments show qualitative and quantitative properties of the estimates and
demonstrate their usefulness in practical applications.Comment: 18 pages, 7 figures, 1 tabl
Topology optimization for incremental elastoplasticity: a phase-field approach
We discuss a topology optimization problem for an elastoplastic medium. The
distribution of material in a region is optimized with respect to a given
target functional taking into account compliance. The incremental elastoplastic
problem serves as state constraint. We prove that the topology optimization
problem admits a solution. First-order optimality conditions are obtained by
considering a regularized problem and passing to the limit
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
- …