Shape optimization for Stokes flows using sensitivity analysis and finite element method

In the context of structural optimization in fluid mechanics we propose a numerical method based on a combination of the classical shape derivative and Hadamard's boundary variation method. Our approach regards the viscous flows governed by Stokes equations with the objective function of energy dissipation and a constrained volume. The shape derivative is computed by Lagrange's approach via the solutions of Stokes and adjoint systems. The programs are written in FreeFem++ using the Finite Element method.