The characteristics of the turbulent MHD channel flows are analysed within the elliptic blending    model. The evaluation of additional terms representing the MDH interactions in the transport equation for the turbulent kinetic energy and dissipation rate is carried out firstly. And then, the improvement of model coefficients for the additional terms is achieved. Regardless of the magnetic field orientation, perpendicular and parrel to the main flow, the model coefficients for the additional terms are not changed to maintain the generality of the present model. The prediction results are directly compared to the DNS data to assess the performance of the present model.