An explicit Loop Tree Duality (LTD) formula for two-loop Feynman integrals with integer power of propagators is presented and used for a numerical UV divergence subtraction algorithm. This algorithm proceeds recursively and it is based on the $\mathcal{R}$ operator and the Hopf algebraic structure of UV divergences. After a short review of LTD and the numerical evaluation of multi-loop integrals, LTD is extended to two-loop integrals with generalized powers of propagators. The $\mathcal{R}$ operator and the tadpole UV subtraction are employed for the numerical calculation of two-loop UV divergent integrals, including the case of quadratic divergences.