Modeling coupled heat transfer and multiphase flow during the expanding solvent steam-assisted gravity drainage (ES-SAGD) process