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

A Two-Level Finite Element Discretization of the Streamfunction Formulation of the Stationary Quasi-Geostrophic Equations of the Ocean

In this paper we proposed a two-level finite element discretization of the nonlinear stationary quasi-geostrophic equations, which model the wind driven large scale ocean circulation. Optimal error estimates for the two-level finite element discretization were derived. Numerical experiments for the two-level algorithm with the Argyris finite element were also carried out. The numerical results verified the theoretical error estimates and showed that, for the appropriate scaling between the coarse and fine mesh sizes, the two-level algorithm significantly decreases the computational time of the standard one-level algorithm.Comment: Computers and Mathematics with Applications 66 201

Refined a posteriori error estimation for classical and pressure-robust Stokes finite element methods

Recent works showed that pressure-robust modifications of mixed finite element methods for the Stokes equations outperform their standard versions in many cases. This is achieved by divergence-free reconstruction operators and results in pressure independent velocity error estimates which are robust with respect to small viscosities. In this paper we develop a posteriori error control which reflects this robustness. The main difficulty lies in the volume contribution of the standard residual-based approach that includes the $L^2$-norm of the right-hand side. However, the velocity is only steered by the divergence-free part of this source term. An efficient error estimator must approximate this divergence-free part in a proper manner, otherwise it can be dominated by the pressure error. To overcome this difficulty a novel approach is suggested that uses arguments from the stream function and vorticity formulation of the Navier--Stokes equations. The novel error estimators only take the $\mathrm{curl}$ of the right-hand side into account and so lead to provably reliable, efficient and pressure-independent upper bounds in case of a pressure-robust method in particular in pressure-dominant situations. This is also confirmed by some numerical examples with the novel pressure-robust modifications of the Taylor--Hood and mini finite element methods

Applying a phase field approach for shape optimization of a stationary Navier-Stokes flow

We apply a phase field approach for a general shape optimization problem of a stationary Navier-Stokes flow. To be precise we add a multiple of the Ginzburg--Landau energy as a regularization to the objective functional and relax the non-permeability of the medium outside the fluid region. The resulting diffuse interface problem can be shown to be well-posed and optimality conditions are derived. We state suitable assumptions on the problem in order to derive a sharp interface limit for the minimizers and the optimality conditions. Additionally, we can derive a necessary optimality system for the sharp interface problem by geometric variations without stating additional regularity assumptions on the minimizing set