Several domain decomposition methods for approximating solutions of parabolic problems are given. These methods rely on implicit Galerkin procedures in the subdomains and explicit flux calculation on the inter-domain boundaries. The procedures are conservative both in the subdomains and across inter-domain boundaries. A priori error bounds and experimental results are presented.