Search Results
Guaranteed error control for the pseudostress approximation of the Stokes equations
2015, Bringmann, Philipp, Carstensen, Carsten, Merdon, Christian
The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H (div) and the velocity in L2. Any standard mixed finite element function space can be utilized for this mixed formulation, e.g. the Raviart-Thomas discretization which is related to the Crouzeix-Raviart nonconforming finite element scheme in the lowest-order case. The effective and guaranteed a posteriori error control for this nonconforming velocity-oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf-sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy.
A goal-oriented dual-weighted adaptive finite element approach for the optimal control of a nonsmooth Cahn-Hilliard-Navier-Stokes system
2016, Hintermüller, Michael, Hinze, Michael, Kahle, Christian, Tobias, Keil
This paper is concerned with the development and implementation of an adaptive solution algorithm for the optimal control of a time-discrete Cahn-Hilliard-Navier-Stokes system with variable densities. The free energy density associated to the Cahn-Hilliard system incorporates the double-obstacle potential which yields an optimal control problem for a family of coupled systems in each time instant of a variational inequality of fourth order and the Navier-Stokes equation. A dual-weighted residual approach for goal-oriented adaptive finite elements is presented which is based on the concept of C-stationarity. The overall error representation depends on primal residuals weighted by approximate dual quantities and vice versa as well as various complementarity mismatch errors. Details on the numerical realization of the adaptive concept and a report on numerical tests are given.