SPH modeling of tidal bore scenarios

Abstract

The paper presented a Smoothed Particle Hydrodynamics (SPH) method to study the three-dimensional (3D) tidal bore scenarios. The SPH method is a mesh free particle modeling technique that can track the large deformation of free surfaces in a straightforward and accurate way. Two benchmark cases of the tidal bore propagation were computed and compared with the experimental results. The first one is related to the undular and breaking bores in a regular open channel and the second one considers the undular bore passing through the contraction of bridge piers. Physical laboratory experiments have also been carried out to validate the numerical investigations. The comparisons of both the free-surface profile and velocity field demonstrated that the SPH technique could provide a very promising tool to simulate tidal bore phenomena in engineering practice. The work is the first to systematically explore the potentials of mesh-free SPH modeling approach in predicting the tidal bore features under 3D flow conditions.

Cylindrical Smoothed Particle Hydrodynamics Simulations of Water Entry

“…The difference is that the symmetric boundary here ensures the fluid particles to move following the tangential direction of the symmetric plane and thus complement the motion of fluid particles on another side. Similar treatment has also been done by Liu et al [28] for modeling the 3D tidal bores. Fig.…”

“…1 (27) By referring to the principles of density re-initialization under a Cartesian coordinate, the following representation can be derived from the division of above two equations as (28) This should be applied every 20 time steps in the computation, as recommended by [24]. Based on the previous computational experiences, the choice of time step seems not to influence this number.…”

Consideration of the Mechanisms for Tidal Bore Formation in an Idealized Planform Geometry

“…In the case of a positive surge generated by the rapid closure of a downstream gate in a laboratory flume with uniform flow, the bulk flow velocity downstream the wave front does not reverse (Viero et al, ). Experimental and numerical analyses showed some transient flow reversal next to the bed and close to the bore passage (Docherty & Chanson, ; Khezri & Chanson, ; Koch & Chanson, ; Leng & Chanson, , ; Lubin et al, ); however, this occurrence is relatively weak and related to flow separation and recirculation due to the upward deviation of the main flow at the foot of the bore (Liu et al, ; Lubin et al, ). In the case of a tidal bore, in which also the bulk flow velocity at the head of the front actually reverses (Furgerot et al, ; Masoud et al, ; Simpson et al, ), the reverse flow is typically much stronger.…”

SPH-FE-Based Numerical Simulation on Dynamic Characteristics of Structure under Water Waves

“…1/2 (12) Properly reducing the initial particle spacing can improve the accuracy of calculation, but too small particle spacing will lead to the instability of computation [32]. In this study, the water fluid deformation is solved by the SPH method.…”

“…For the SPH method, the particles move in Lagrangian coordinates and the advection is directly calculated by particle motion without numerical diffusion [9]. As SPH is a type of mesh-free method, it is inherently suitable for solving the procedures of moving discontinuities with large deformations [3], and it has been utilized to resolve a wide range of hydrodynamics problems [10][11][12][13]. Rafiee et al [14] used the pressure Poisson equation in the SPH method to solve the incompressible problem of fluid, and verified that the method can well simulate the dynamic response between the free surface of fluid and the structure; Dao et al [15] established a numerical wave flume based on SPH method, which was successfully applied to simulate the cases of tsunami wave propagating on the slope and wave impacting on coastal structures; and St-Germain et al [16] used the weak compressibility SPH method to compute the wave force of the rapidly developing tsunami wave acting on the square cylinder, and the computed results fit well with the results of large-scale physical experiments.…”