Abstract
The penalty method used to solve Darcy equations permits the decoupling of the velocity and the pressure. In this work, we propose the spectral element discretization to implement an algorithm based on the optimization of the penalty parameter. We present some numerical results that confirm the adequation of the used penalty method.