The open-source code DualSPHysics is used in this work to compute the wave run-up in an existing dike in the Chinese coast using realistic dimensions, bathymetry and wave conditions. The GPU computing power of the DualSPHysics allows simulating real-engineering problems that involve complex geometries with a high resolution in a reasonable computational time. The code is first validated by comparing the numerical free-surface elevation, the wave orbital velocities and the time series of the run-up with physical data in a wave flume. Those experiments include a smooth dike and an armored dike with two layers of cubic blocks. After validation,the code is applied to a real case to obtain the wave run-up under different incident wave conditions. In order to simulate the real open sea, the spurious reflections from the wavemaker are removed by using an active wave absorption technique. 相似文献
In this paper, some recent developments and new results concerning the trapping of waves by arrays of vertical circular cylinders is presented. In particular, the cases are examined when there is a circular arrangement of cylinders and both finite and infinite periodic linear arrays of identical cylinders. Only for the infinite array is there pure trapping of waves – known as Rayleigh–Bloch or edge waves – which, for particular dominant wavenumbers, reduce to the well-known trapped-mode solutions for a cylinder between two parallel walls having either Neumann or Dirichlet conditions upon them. This latter case is considered separately and some new results are presented. In the circular array and finite linear array the concept of near-trapping is introduced where large resonant motions are found to occur at certain frequencies of the incident wave field. In the case of the finite linear array, these near-trapping frequencies are related to the Rayleigh–Bloch trapped-wave frequencies for the infinite array. Finally, the case when there are two or more lines of cylinders in the linear array is examined. 相似文献
The numerical model for nonlinear wave propagation in the physical space, developed by Grilli, et al.12,13, uses a higher-order BEM for solving Laplace's equation, and a higher-order Taylor expansion for integrating in time the two nonlinear free surface boundary conditions. The corners of the fluid domain were modelled by double-nodes with imposition of potential continuity. Nonlinear wave generation, propagation and runup on slopes were successfully studied with this model. In some applications, however, the solution was found to be somewhat inaccurate in the corners and this sometimes led to wave instability after propagation in time.
In this paper, global and local accuracy of the model are improved by using a more stable free surface representation based on quasi-spline elements and an improved corner solution combining the enforcement of compatibility relationships in the double-nodes with an adaptive integration which provides almost arbitrary accuracy in the BEM numerical integrations. These improvements of the model are systematically checked on simple examples with analytical solutions. Effects of accuracy of the numerical integrations, convergence with refined discretization, domain aspect ratio in relation with horizontal and vertical grid steps, are separately assessed. Global accuracy of the computations with the new corner solution is studied by solving nonlinear water wave flows in a two-dimensional numerical wavetank. The optimum relationship between space and time discretization in the model is derived from these computations and expressed as an optimum Courant number of 0.5. Applications with both exact constant shape waves (solitary waves) and overturning waves generated by a piston wavemaker are presented in detail. 相似文献