Higher order energy decay for damped wave equations with variable coefficients