The flow characteristics and the heat transfer rate on a surface by the interaction of a pair of vortices are studied numerically. To analyze the common flow up produced by vortex generators in a rectangular channel flow, the pseudo-compressibility viscous method is introduced into the Reynolds-averaged Navier-Stokes equation for 3-dimensional unsteady, incompressible viscous flows. To predict turbulence characteristics, a two-layer $k-varepsilon$ turbulence model is used on the flat plate 3-dimensional turbulence boundary The computational results predict accurately Reynolds stress, turbulent kinetic energy and flow field generated by the vortex generators. The numerical results, such as thermal boundary layers, skin friction characteristics and heat transfers, are also reasonably close to the experimental data.