Simulation of one-dimensional evolution of wind waves in a deep water

A direct wave model based on the one-dimensional nonlinear equations for potential waves is used for simulation of wave field development under the action of energy input, dissipation, and nonlinear wave-wave interaction. The equations are written in conformal surface-fitted nonstationary coordinate system. New schemes for calculating the input and dissipation of wave energy are implemented. The wind input is calculated on the basis of the parameterization developed through the coupled modeling of waves and turbulent boundary layer. The wave dissipation algorithm, introduced to prevent wave breaking instability, is based on highly selective smoothing of the wave surface and surface potential. The integration is performed in Fourier domain with the number of modes M = 2048, broad enough to reproduce the energy downshifting. As the initial conditions, the wave field is assigned as train of Stokes waves with steepness ak = 0.15 at nondimensional wavenumber k = 512. Under the action of nonlinearity and energy input the spectrum starts to grow. This growth is followed by the downshifting. The total time of integration is equal to 7203 initial wave periods. During this time the energy increased by 1111 times. Peak of the spectrum gradually shifts from wavenumber nondimensional k = 512 down to k = 10. Significant wave height increases 33 times, while the peak period increases 51 times. Rates of the peak downshift and wave energy evolution are in good agreement with the JONSWAP formulation