If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.
With the method of separation of variables and the eigenfunction expansion employed, an analytical solution is presented for the radiation and diffraction of a rectangular structure with an opening near a vertical wall in oblique seas, in which the unknown coefficients are determined by the boundary conditions and matching requirement on the interface. The effects of the width of the opening and the angle of incidence on the hydrodynamic characteristics of a rectangular structure with an opening near a vertical wall are mainly studied. The comparisons of the calculation results with wall-present and with wall-absent are also made. The results indicate that the variation trends of the heave added mass and excitation force with wall-present are almost the same as those with wall-absent, and that the peak values in the former case are obviously larger than those in the latter due to the reflection of the vertical wall.
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 em- ployed 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 ma- trix can make the transformed equations relatively concise, the treatment of lateral boundary conditions easier and the de- velopment 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.
