Comparison of Shape Derivatives using CutFEM for Ill-posed Bernoulli Free Boundary Problem

by   Erik Burman, et al.

In this paper we discuss a level set approach for the identification of an unknown boundary in a computational domain. The problem takes the form of a Bernoulli problem where only the Dirichlet datum is known on the boundary that is to be identified, but additional information on the Neumann condition is available on the known part of the boundary. The approach uses a classical constrained optimization problem, where a cost functional is minimized with respect to the unknown boundary, the position of which is defined implicitly by a level set function. To solve the optimization problem a steepest descent algorithm using shape derivatives is applied. In each iteration the cut finite element method is used to obtain high accuracy approximations of the pde-model constraint for a given level set configuration without re-meshing. We consider three different shape derivatives. First the classical one, derived using the continuous optimization problem (optimize then discretize). Then the functional is first discretized using the CutFEM method and the shape derivative is evaluated on the finite element functional (discretize then optimize). Finally we consider a third approach, also using a discretized functional. In this case we do not perturb the domain, but consider a so-called boundary value correction method, where a small correction to the boundary position may be included in the weak boundary condition. Using this correction the shape derivative may be obtained by perturbing a distance parameter in the discrete variational formulation. The theoretical discussion is illustrated with a series of numerical examples showing that all three approaches produce similar result on the proposed Bernoulli problem.



There are no comments yet.


page 1

page 2

page 3

page 4


Maximizing Vortex for the Navier–Stokes Flow with a Convective Boundary Condition: A Shape Design Problem

In this study, a shape optimization problem for the two-dimensional stat...

φ-FEM, a finite element method on domains defined by level-sets: the Neumann boundary case

We extend a fictitious domain-type finite element method, called ϕ-FEM a...

Shape optimization for superconductors governed by H(curl)-elliptic variational inequalities

This paper is devoted to the theoretical and numerical study of an optim...

Model Hierarchy for the Shape Optimization of a Microchannel Cooling System

We model a microchannel cooling system and consider the optimization of ...

Vorticity Maximization of a Linear Fluid Flow via Volume Constrained and Perimeter Regularized Shape Optimization

We study an optimization problem that aims to determine the shape of an ...

Defect reconstruction in a 2D semi-analytical waveguide model via derivative-based optimization

This paper considers the reconstruction of a defect in a two-dimensional...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.