Heat liberation calculation