Fully developed incompressible turbulent flow in a two-dimensional channel has been simulated numerically. Turbulence promoters were arranged intermittently on one side of the channel walls. The two walls and promoters were heated individually. A large-scale flow field and a temperature field were calculated by integrating the filtered, three-dimensional time-dependent Navier-Stokes equations and time-dependent energy equation. Motions of turbulence field with small scales were simulated using eddy-viscosity models. The Reynolds and Peclet numbers were set at 1.1×l0⁴ and 7.7×l0^, respectively. Only 9, 000 grid points were distributed and three kinds of averaging were used in order to complement the grid points. Mean-velocity profiles, turbulence statistics, time-dependent features of flow and mean-temperature profiles are in good agreement with the experimental data. Calculations were carried out using a HITAC M 200/280 H system computer of the University of Tokyo.