Modeling of multiphase transport of multicomponent organic contaminants and heat in the subsurface: Numerical model formulation