A large eddy simulation (LES) of a plane turbulent channel flow is performed at a Reynolds number Re? = 590 based on the channel half width, ? and wall shear velocity, u? by approximating the near wall region using differential equation wall model (DEWM). The simulation is performed in a computational domain of 2?? x 2? x ??. The computational domain is discretized by staggered grid system with 32 x 30 x 32 grid points. In this domain the governing equations of LES are discretized spatially by second order finite difference formulation, and for temporal discretization the third order low-storage Runge-Kutta method is used. Essential turbulence statistics of the computed flow field based on this LES approach are calculated and compared with the available Direct Numerical Simulation (DNS) and LES data where no wall model was used. Comparing the results throughout the calculation domain we have found that the LES results based on DEWM show closer agreement with the DNS data, especially at the near wall region. That is, the LES approach based on DEWM can capture the effects of near wall structures more accurately. Flow structures in the computed flow field in the 3D turbulent channel have also been discussed and compared with LES data using no wall model.