Hydrodynamic design optimization using non stationary porous media model

Abstract In this paper, we focus on the penalty finite element method for the non stationary porous media model. We begin by showing the existence and uniqueness of the solution for the initial problem. Error estimates for the velocity and the pressure are obtained via the energy method. We introduce a time discretization by the use of a backward Euler scheme combined with fully discrete finite element method to approximate the penalized problem and establish an error estimate for the velocity and the pressure which will be used to show the convergence of the approximate solution to the solution of the initial problem. The shape optimization problem is to find the shape which is optimal in that it minimizes a cost functional related to a comfort fish populations. We derive the adjoint system associated to the penalized problem. We compute the gradient in terms of state and adjoint variables. The optimization procedure is implemented using the continuous adjoint method and the finite element method. Numerical simulations are presented to show the efficiency and the robustness of the proposed method.