Five versions of low Reynolds number k-epsilon models including the origina
l one have been tested against the DNS data for 1D oscillatory boundary lay
ers; one with sinusoidal and another with flat-crested free-stream velocity
variation. The comparison has been made for the cross-stream velocity, tur
bulent kinetic energy, Reynolds stress, and wall shear stress. It is found
that the original model is superior to the more recent models by virtue of
its numerically desirable boundary conditions and accuracy of prediction fo
r mean flow properties. A brief description about the wall-limiting behavio
r of important turbulent quantities is presented to recognize physically co
rrect k-epsilon models.