相变蓄热的数值模拟与热力学优化