Numerical discretization of a Brinkman-Darcy-Forchheimer model under singular forcing