In this paper, a coupled thermo-hydro-mechanical (THM) simulation in a faulted deformable porous medium is presented. This model involves solving the mass conservation, linear momentum balance, and energy balance equations which are derived from the Biot's consolidation theory. Fluid pore pressure, solid displacement, and temperature are chosen as initial variables in these equations, and the finite element method in combination with the interface element is used for spatial discretization of continuous and discontinuities (fault) parts of the medium to solve the equations. The main purpose of this study is providing precise formulations, applicability, and ability of the triple-node zero-thickness interface element in THM modeling of faults. It should be noted that the system of equations is solved using a computer code written in Matlab program. In order to verify the developed method, simulations of index problems such as Mandel's problem, and coupled modeling of a faulted porous medium and a faulted aquifer are presented. The modeling results obtained from the developed method show a very good agreement with those by other modeling methods, which indicates its accuracy.