| Based on the Naviertokes equation, a numerical model to predict vertical two-dimensional turbulent wave field has been developed in the present thesis. The key technology introduced and implemented is as following: 1. A modified k?t turbulence model was derived to deal with the closure problem by using the Renolds-Stress formula proposed by Zhang Changgao. The finitefference method is used to approximate the predicting equations, in which the convective term of inner points is discretized by using third order upwind difference scheme and a linear combination of one order upwind difference scheme and centered difference scheme is applied to discretize the convective term of boundary points. 2. The modified VOF method is adopted to deal with free surface. The governing equation for the Volume of Fluid is solved by Donorcceptor method. The calculation results show that it can predict the surface transformation of wave over the irregular landform very well. 3. A partial cell parameter is introduced to treat irregular boundary, which is based on the finite ifference method. For the grid nodes immediately adjacent to the wall boundary, the viscosity sublayer is treated by wall function approximation. The transmittal of linear waves over flat bottom, the transformation of Cnoidal waves over a submerged bar and the transformation of regular waves over a slope are simulated by the model. The results are compared with theoretical and experimental values respectively. It shows that the calculation results agree reasonably with the theoretical and experimental values. |