For the simulation of the nonlinear wave propagation in coastal areas with complex boundaries,a numerical model is developed in curvilinear coordinates. In the model,the Boussinesq-type equations including the dissipation terms are employed as the governing equations. In the present model,the dependent variables of the transformed equations are the free surface elevation and the utility velocity variables,instead of the usual primitive velocity variables. The introduction of utility velocity variables which are the products of the contravariant components of the velocity vector and the Jacobi matrix can make the transformed equations relatively concise,the treatment of lateral boundary conditions easier and the development of the program simpler. The predictor-corrector method and five-point finite-difference scheme are employed to discretize the time derivatives and the spatial ones,respectively. The numerical model is tested for three cases. It is found that the numerical results are in good agreement with the analytical results and experimental data.
展开▼