We need to calculate multiphase stream of the one substance in the air atmosphere (this is why we using interFoam), but also we need to get heat transfer data between the substance, solid obstacle (moving reference frame) and atmosphere.
You might have to deal with the source code. My suggestion is to add 'alpha' to chtMultiRegionFoam. You may also try to modify 'multiphaseInterFoam' to solve 'T' and set one of the phase as a solid, which seems easier. Good luck in your new chtMultiRegionInterFoam!
I've done some modification on OpenFOAMv1912 solvers including a modification of chtMutliRegionFoam with interFoam. You may have a look here: https://github.com/c-schubert/OpenFOAM--MHD--VoF-Solver-Collection (no guarantee given that everything is 100% correct of course ...)
I guess the solver speeds by now could further be improved by doing the same for the new OpenFOAM8 release, which comes with some new convergence checking ability for the multiRegion solver and improved PLIC methods.