A computational program ［TDET］ of the particle transport equation is developed on radiation shielding problem in two-dimensional cartesian geometry based on the discrete element method. Not like the ordinary discrete ordinates method, the quadrature set of angles is not fixed but steered by the spatially dependent angular fluxes. The angular dependence of the scattering source term in the particle transport equation is described by series expansion in spherical harmonics, and the energy dependence of the particles is considered as well. Three different benchmark tests are made for verification of TDET : For the ray effect analysis on a square absorber with a flat isotropic source, the results of TDET calculation are quite well conformed to those of MORSE-CG calculation while TDET ameliorates the ray effect more effectively than S$_{N}$ calculation. In the analysis of the streaming leakage through a narrow vacuum duct in a shield, TDET shows conspicuous and remarkable results of streaming leakage through the duct as well as MORSE-CG does, and quite better than S$_{N}$ calculation. In a realistic reactor shielding situation which treats in two cases of the isotropic scattering and of linearly anisotropic scattering with two groups of energy, TDET calculations show local ray effect between neighboring meshes compared with S$_{N}$ calculations in which the ray effect extends broadly over several meshes.eshes.