The computation of higher order QCD corrections to hard scattering processes is made extremely complicated by the presence of infrared (soft and collinear) singularities that appear at intermediate stages of the calculation and prevent a straightforward implementation of numerical techniques. We describe an extension of the subtraction formalism that can be used to perform fully exclusive NNLO QCD computations for a wide variety of important processes at hadron colliders. We present selected numerical results for vector and Higgs boson production, associated production of the Higgs boson with a W boson, and prompt photon pair production at the LHC.