A numerical model of seawater intrusion in aquifers