Finite element solution of Navier-Stokes equations in shallow domains

A finite element method for the numerical solution of the anisotropic Navier-Stokes equations in shallow domain is presented. This method take into account the low regularity of the vertical component of the velocity in the hydrostatic approximation of the Navier-Stokes equations [2, 3, 5]. A projection method [8] is used for the time discretization. The linear systems are solved via a some preconditioned conjugate algorithm, well adapted to massively parallel computers [4]. Some results are presented for the wind driven water circulation in lakes Geneva and Neuchâtel.