SUMMARYLow Reynolds number turbulence stress and heat flux equation models (LRSFM) have been developed to enhance predictive capabilities. A new method is proposed for providing the wall boundary condition for dissipation rate of turbulent kinetic energy, , to improve the model capability upon application of coarse meshes for practical use. The proposed method shows good agreement with accepted correlations and experimental data for ows with various Reynolds and Prandtl numbers including transitional regimes. Also, a mesh width about 5 times or larger than that used in existing models is applicable by using the present boundary condition. The present method thus enhanced computational e ciency in applying the complex turbulence model, LRSFM, to predictions of complicated ows. Unsteady Reynolds averaged Navier-Stokes (URANS) computations are conducted for an oscillatory non-isothermal quasiplanar triple-jet. Comparisons are made between an experiment and predictions with the LRSFM and the standard k-model. A water test facility with three vertical jets, the cold in between two hot jets, simulates temperature uctuations anticipated at the outlet of a liquid metal fast reactor core. The LRSFM shows good agreement with the experiment, with respect to mean proÿles and the oscillatory motion of the ow, while the k-model under-predicts the mixing due to the oscillation, such that a transverse mean temperature di erence remains far downstream.