Hlavním cílem této diplomové práce je formulace, aproximace a numerické řešení Darcyovy-Forchheimerovy rovnice. Pro numerické řešení byly vybrány tři linearizační metody; Picardova iterace, Newtonova metoda a metoda L-schéma. Všechny tři linearizační metody byly naimplementovány do softwaru Flow123d. Metody byly nejdříve otestovány na jednoduchém analytickém řešení a na problému dvou kolmých puklin. Nakonec byla vybrána jedna linearizační metoda a otestována na benchmarkovém problému.