2,111 research outputs found
Numerical approximation of phase field based shape and topology optimization for fluids
We consider the problem of finding optimal shapes of fluid domains. The fluid
obeys the Navier--Stokes equations. Inside a holdall container we use a phase
field approach using diffuse interfaces to describe the domain of free flow. We
formulate a corresponding optimization problem where flow outside the fluid
domain is penalized. The resulting formulation of the shape optimization
problem is shown to be well-posed, hence there exists a minimizer, and first
order optimality conditions are derived.
For the numerical realization we introduce a mass conserving gradient flow
and obtain a Cahn--Hilliard type system, which is integrated numerically using
the finite element method. An adaptive concept using reliable, residual based
error estimation is exploited for the resolution of the spatial mesh.
The overall concept is numerically investigated and comparison values are
provided
Optimal stability polynomials for numerical integration of initial value problems
We consider the problem of finding optimally stable polynomial approximations
to the exponential for application to one-step integration of initial value
ordinary and partial differential equations. The objective is to find the
largest stable step size and corresponding method for a given problem when the
spectrum of the initial value problem is known. The problem is expressed in
terms of a general least deviation feasibility problem. Its solution is
obtained by a new fast, accurate, and robust algorithm based on convex
optimization techniques. Global convergence of the algorithm is proven in the
case that the order of approximation is one and in the case that the spectrum
encloses a starlike region. Examples demonstrate the effectiveness of the
proposed algorithm even when these conditions are not satisfied
Optimal Control of the Thermistor Problem in Three Spatial Dimensions
This paper is concerned with the state-constrained optimal control of the
three-dimensional thermistor problem, a fully quasilinear coupled system of a
parabolic and elliptic PDE with mixed boundary conditions. This system models
the heating of a conducting material by means of direct current. Local
existence, uniqueness and continuity for the state system are derived by
employing maximal parabolic regularity in the fundamental theorem of Pr\"uss.
Global solutions are addressed, which includes analysis of the linearized state
system via maximal parabolic regularity, and existence of optimal controls is
shown if the temperature gradient is under control. The adjoint system
involving measures is investigated using a duality argument. These results
allow to derive first-order necessary conditions for the optimal control
problem in form of a qualified optimality system. The theoretical findings are
illustrated by numerical results
- …