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.