The small open area available at the slots of underdrains in pressurized granular bed filters for drip irrigation implies: (1) the existence of a region with non-uniform flow, and (2) local values of modified particle Reynolds number >500. These flow conditions may disagree with those accepted as valid for common pressure drop-flow rate correlations proposed for packed beds. Here, we carried out detailed computational fluid dynamics (CFD) simulations of a laboratory filter to analyze the results obtained with five different equations of head losses in porous media: (1) Ergun, (2) Darcy-Forchheimer, (3) Darcy, (4) Kozeny-Carman and (5) power function. Simulations were compared with experimental data at different superficial velocities obtained from previous studies. Results for two silica sand media indicated that all equations predicted total filter pressure drop values within the experimental uncertainty range when superficial velocities <38.3 m h−1. At higher flow rates, Ergun equation approximated the best to the observed results for silica sand media, being the expression recommended. A simple analytical model of the pressure drop along flow streamlines that matched CFD simulation results was developed.