Abstract: Las ecuaciones de Navier-Stokes promediadas en la vertical se resuelven numéricamente mediante volúmenes finitos en un caso de flujo formado por dos fluidos de diferente densidad y viscosidad sobre un lecho irregular fijo. Obtenido el modelo analítico como un sistema de ecuaciones hiperbólico con términos fuente, se aplica un esquema upwind de primer orden de Roe y un esquema FORCE para su resolución numérica programando en lenguaje C. La simulación de casos test y experimentales presentes en la literatura se utiliza para la evaluación y validación del modelo. Este trabajo se ha desarrollado con una beca JAE del CSIC de iniciación a la investigación en el marco del proyecto: Desarrollo y validación de modelos de simulación avanzados con base física para procesos geofísicos e hidráulicos. SIMGEOAGUA. Ministerio de Ciencia e Innovación. CGL2015-66114-R.