AbstractGrid adaptive methods combined with domain adaptation are discussed for two‐dimensional seepage flow problems with free boundaries through porous media. Examples of grid and domain adaptive methods are presented to demonstrate several ways to predict grids and shapes of free boundaries using an iterative scheme. Finally, the combined adaptive methods are applied to obtain smooth non‐oscillatory shape of a free boundary of seepage flow through non‐homogeneous porous