Numerical discretization of a Darcy-Forchheimer problem coupled with a singular heat equation