The periodically fully developed laminar heat transfer and fluid flow characteristics of an array of uniform plate length, aligned at angles to the air flow direction have been investigated numerically in the range Re=50-2300 using body-fitted coordinate system. Multisurface transformation is used in the grid generation in order to retain the periodic corresponding relations of the surfaces or lines in both physical and computational spaces. The computed results exhibit good agreement with the experimental data available in the literature. The results also show that both the intensity of heat transfer and the pressure drop increase with the increases in oblique angle and plate length.

Three practices of treating outflow boundary condition were adopted in computations for convective heat transfer of a two-dimensional jet impinging in a rectangular cavity. The three practices were local mass conservation method, local one-way method and fully developed assumption. The numerical solutions of the three methods were compared with test data obtained via. naphthalene sublimation technique. It was found that the fully developed assumption was inappropriate, and the local one-way method could provide reasonably good results for the cavity bottom, while for the lateral wall the results with this method qualitatively differed from the test data. The solution with the local mass conservation method was the best. It thus suggested that for a problem expected with a strong recirculating flow at the exit of the computation domain, the local mass conservation method be adopted to treat the outflow boundary condition.

Fluid flow and heat transfer in two-dimensional wavy channels of uniform cross-sectional area were analyzed numerically for constant-property Laminar flow. The channel walls were maintained at constant temperature. Streamwise pertodic geometry caused the flow and heat transfer to repeat periodically far downstream of the channel entrance. Computations were performed for different Reynolds numbers and geometric parameters. Prandtl number was kept constant (0.707). The whole solution domain was divided into five regions, and different coordinates were used for each region. The periodic boundary conditions were implemememted by mutual repiacements of field values at the two end regions. It was found that both Nussel: number and friction factor were greatly affected by Reynolds number and two geometric parameters.