Abstract: ased on the universal k-ε turbulent model,we trace free surface with the volume of fluid(VOF) method for taking the character of distinct surface elevation,and establish a two-dimensional model of vertical plane turbulent jet discharged from bottom of shallow water.Many numerical simulation experiments are implemented for differential groups of water depth,velocity and width of jets and the whole flow field is analyzed.We find that the distribution characters of velocity vector field and that of surface elevation are corresponding.A self-preservation law of surface elevation distribution around jet axis is found and an empirical formula is established.Due to the surface retardation,the attenuation rate of velocity becomes quicker,which depends on water depth and jet width.