A Lyapunov Function Approach to Energy Based Model Reduction