In this paper, Forchheimer equation is used as the constitutive equation for flow through rockfill, and the non-linear two-dimensional governing equation with free surface is solved by a new finite element method in a fixed grid. The model is verified by applying it to different flow conditions. The first scenario, which is assumed to be one-dimensional with analytical solution available for it...