Se presenta la formulación de un modelo matemático, hidrodinámico y de transporte por difusión, y su solución numérica por el método de los elementos finitos. el modelo parte de las ecuaciones básicas bidimensionales de Navier-Stokes integradas a lo largo de la profundidad, válidas en un principio para aguas poco profundas, y de la ecuación de difusión Euleriana, a la que se incorpora una ecuación de estado que expresa la relación entre la densidad y la sustancia a difundir.