With Managed Pressure Drilling (MPD) becoming increasingly important to oil and gas well construction, there has been a growing demand for comprehensive models that can accurately describe the well hydraulic dynamics, especially for complex well geometries, High-Pressure High-Temperature (HPHT) conditions, during well control events, etc. To meet this demand, new multi-phase hydraulic models have been developed. However, despite their improvement over the traditional lumped parameter models, these models typically do not consider detailed temperature effects, which markedly influence well hydraulics and need to be explicitly considered for accurate MPD control.The development of a new integrated multi-phase thermal and hydraulic modeling framework is presented in this paper. This integrated model can estimate the mud temperature in the drillstring and the annulus, as well as the temperature variations of the formation during complex well control situations. The model uses a semi-implicit finite volume approach and solves the mixture energy equation for the wellbore fluids, assuming that all the phases are in thermal equilibrium. Heat transfer between the drillstring and wellbore fluids, and between the wellbore fluids and the formation is calculated using thermal resistance networks.The steady-state and transient temperature behaviors during mud circulation are compared to the steady-state analytical model by Hasan and Kabir and commercial software. These comparisons show good agreement for both steady-state and transient cases. To demonstrate the importance of accurate temperature estimation of the drillstring and annulus fluids in HPHT conditions, offshore kick scenarios for non-aqueous drilling fluids and dynamic kick control methods are simulated. Simulation results indicate that the model enables real-time estimation of crucial parameters during well control, such as the wellbore pressure and temperature profiles, increased outflow and pit gain during kicks, gas thermodynamic behavior including solubility and unloading behavior at low pressure conditions, gas rising velocity, and temperature-dependent formation strength.