  We consider four turbulent models to simulate the boundary mixing layer of
the ocean. We show the existence of solutions to these models in the
steady-state case then we study the mathematical stability of these solutions.
