Dissolved oxygen convection and diffusion numerical simulation of stokes wave

In order to study the reoxygenation of the coastal wave, the two-dimensional renormalized-group (2-D RNG) k-epsilon (γ-ε) mathematical model and Volume of Fluid (VOF) method are employed to compute the motion of a stokes wave. Dissolved oxygen convection and diffusion model is established to simulate the dissolved oxygen concentration with user defined scalar method. The computational results of above model are compared with analytical solutions, and they agree with each other approximately. Based on the analysis of computational results, we can found dissolved oxygen concentration decreases with the increase of water depth. The average dissolved oxygen concentration of vertical line reaches the minimum at wave crest position, and the average dissolved oxygen concentration of vertical line reaches the maximum at wave trough position. Dissolved oxygen concentration varies much more near the wave surface than other positions.