The Fukushima Daiichi Nuclear Power Plant accident showed a significant effect on the environment due to the release of a large amount of radioactive material. To prevent the damage of the reactor pressure vessel (RPV), an in-vessel melt retention (IVR) by external reactor vessel cooling is considered to be a key severe accident management strategy. For its success, investigation of the thermal-mechanical behavior of the RPV lower head is of importance. The main objective of this study is the development of a thermal structural calculation tool for simulating the failure process and the visco-plastic behavior of the RPV lower head wall during the late phase of a core-melt severe accident. OpenFOAM, an open-source toolbox for developing numerical solvers, was used for the calculation platform. Thermal behavior of the molten pool in the RPV lower head is simulated by the phasechange effective convectivity model (PECM). One of the LIVE (Late In-Vessel Phase Experiments) experiments was analyzed for an evaluation purpose and reasonable results were obtained. The solver was then extended to couple PECM with a structural analysis model that considers thermal expansion, plasticity, creep and material damage. Two FOREVER (Failure Of Reactor Vessel Retention) tests, each of which uses different type of steel, were analyzed: EC-FOREVER-2 with French RPV steel 16MND5 and EC-FOREVER-4 with American RPV steel SA533B1. The deformation and failure time agreed reasonably well with the measured data. In addition, the failure mode was also well predicted qualitatively. The numerical analysis showed that the developed tool has a capability of simulating the lower head thermal structural behavior.