El presente articulo presenta la modelacion numerica de problemas de flujo en presas de tierra, empleando dos tipos de estrategias. La primera estrategia consiste en modelar el dominio saturado de la presa, delimitando una superficie libre inicial arbitraria o aproximada, por medio de la parabola de Koseni, y realizando iteraciones de la superficie libre mediante la adaptacion de un esquema de optimizacion que minimiza el error en el calculo y re-enmalla cada iteracion hasta alcanzar resultados adecuados. La segunda estrategia es el modelado transitorio del problema de flujo hasta alcanzar el estado estacionario; momento en el cual queda determinada la superficie libre en la malla, que permanece fija durante el proceso de modelado. La primera estrategia se implementa en el programa multifisico denominado Ansys, y, para ello, se desarrolla un codigo en lenguaje APDL, con el fin de adaptar un esquema de optimizacion que itere la posicion de la superficie libre. La segunda estrategia se implementa en el programa geotecnico de elementos finitos, Slide. Las notables diferencias entre la superficie libre obtenidas por Slide y Ansys, cuantitativa y cualitativamente, llaman la atencion de los disenadores, ya que la posicion de la superficie libre afecta el calculo de factores de seguridad en posibles superficies de falla.