“…Due to the complex THM coupled responses considering two-phase flow in porous media, intensive numerical coupled models have been performed with different codes, including OpenGeoSys (Bilke et al, 2019;Kolditz et al, 2012), Tough-Flac (Rutqvist et al, 2001(Rutqvist et al, , 2014, Code-bright (Li et al, 2013;Olivella et al, 1996), Comes-Geo (Cao et al, 2016;Sanavia et al, 2006), Compass (Cleall et al, 2006;Thomas et al, 1998) and others (Bond and Watson, 2008;Chen et al, 2012;Garitte et al, 2017;Rutqvist, 2015;Wu et al, 2004). Through combing advantages of Tough and Flac, respectively doing well with problems of multiphase-multicomponent flow / heat transfer and processes of soil deformation, it will be more efficient to establish complicated coupled models in Tough-Flac (Rutqvist et al, 2014).…”