An Efficient Two-Layer Non-Hydrostatic Model for Investigating Wave Run-Up Phenomena

In this paper, we study the maximum run-up of solitary waves on a sloping beach and over a reef through a non-hydrostatic model. We do a modification on the non-hydrostatic model derived by Stelling and Zijlema. The model is approximated by resolving the vertical fluid depth into two-layer system. In contrast to the two-layer model proposed by Stelling, here, we have a block of a tridiagonal matrix for the hydrodynamic pressure. The equations are then solved by applying a staggered finite volume method with predictor-corrector step. For validation, several test cases are presented. The first test is simulating the propagation of solitary waves over a flat bottom. Good results in amplitude and shape preservation are obtained. Furthermore, run-up simulations are conducted for solitary waves climbing up a sloping beach, following the experimental set-up by Synolakis. In this case, two simulations are performed with solitary waves of small and large amplitude. Again, good agreements are obtained, especially for the prediction of run-up height. Moreover, we validate our numerical scheme for wave run-up simulation over a reef, and the result confirms the experimental data.