In this paper we are concerned with numerical methods for nonlinear time-dependent problem coupled by electron, ion and photon temperatures in two dimensions, which is called the 2D-3T heat conduction equations. We propose discontinuous Galerkin (DG) methods for the discretization of the equations. For solving the resulting discrete system, we employ two domain decomposition (DD) preconditioners, one of which is associated with the non-overlapping DDM and the other is based on DDM with small overlap. The preconditioners are constructed by dropping the couplings between particles and each preconditioner consists of three preconditioners with smaller matrix size. To gauge the efficiency of the preconditioners, we test two examples and make different settings of parameters. Numerical results show that the proposed preconditioners are very effective to the 2D-3T problem.