Theoretical method for generating solitary waves using plunger-type wavemakers and its Smoothed Particle Hydrodynamics validation
Introduction
Tsunamis are huge ocean waves triggered by the displacement of a large volume of water, mainly caused by submarine earthquakes [1, 2], volcanic eruptions [3, 4], and landslides [5, 6]. Tsunamis, once reaching coastal areas, can destroy the vast majority of marine and coastal structures and more seriously bring about enormous loss of life. A tsunami is generally represented by a solitary wave, but it should be clear that it is not a thoroughly realistic representation. A tsunami is a wave train, of which the leading wave can be either positive (i.e. surface elevation above the still water level (SWL)) or negative (i.e. surface elevation below the SWL) [7], [8], [9], while a solitary wave is single and always positive. Moreover, there is usually insufficient ocean space for an initial displacement of sea surface evolving into a solitary wave [10, 11]. Notwithstanding the differences, the solitary wave approximation has been widely used for tsunami studies [12], [13], [14], [15], [16], [17], [18] mainly due to its simplicity and high repeatability [19, 20].
It is common practice to generate a solitary wave using the method of Goring [21]. That is to match the velocity of piston-type wavemaker with the depth-averaged horizontal velocity of water particles adjacent to the piston. However, the generated wave exhibits a rapid decrease in height during its propagation, exceeding the theoretical level attributed to the boundary layer effect of the water tank [22]. To enhance the wave invariance, Katell and Eric [23] replaced the Boussinesq assumption [24] adopted by Goring [21] with the Rayleigh one [25]. Malek-Mohammadi and Testik [26] considered the evolving nature of wave generation and derived a time-dependent wave celerity based on the conservation of momentum. Both two measures help to generate a higher-quality wave with a shorter transition distance and less attenuation of wave height. The redundant trailing waves following the main pulse are another focus of the solitary wave generation. Although their heights under Boussinesq assumption are lower than 10% of the height of the main pulse, the ratios can be further suppressed to below 3% if Rayleigh assumption is adopted [23]. Farhadi et al. [27] compared the solitary wave generation under six wave theories, namely Boussinesq assumption [24], 1st-order and 2nd-order shallow water equations [28, 29], Rayleigh assumption [25], 3rd-order Grimshaw solution [30], and 9th-order Fenton solution [31]. Their results supported Katell and Eric [23] that Rayleigh assumption yields the least attenuation of wave height but also pointed out an unfavorable lead of wave phase. By comparison, the generated waves based on 3rd- and 9th-order solutions have better overall performance.
Although substantial theories have been established for piston-type solitary wave generation, they cannot be applied in laboratory tanks equipped with plunger-type wavemakers. However, there are indeed a large number of plunger-type wavemaker being used around the world (see our survey results in Appendix A). A plunger only disturbs near-surface water, thereby performing higher efficiency in a deep-water circumstance than a piston. It also saves horizontal space in the wave direction. Therefore, one cannot rule out the possibility that in the future more plunger-type wavemakers would be installed. Unfortunately, owing to the complex geometry, plunger-type wavemaker is merely designed to generate regular waves and multi-frequency random waves so far [32], [33], [34], [35]. There is hardly any research on plunger-type solitary wave generation. Thus, functionalities of the existing and future plunger-type wavemakers are discounted.
In fact, in the earliest study of solitary wave, Russell [36] has put forward the idea of generating a solitary wave by dropping a heavy object into the water. However, a freely falling object generally results in an unsteady wave train [37, 38]. A desired waveform can only be obtained after a long-distance evolution or it does not emerge at all. Considering the precious laboratory space or computational domain, the generated wave is hoped to achieve its target form in the shortest distance. Moreover, the trailing waves should be as weak as possible to enhance the wave pureness. Therefore, manual control of the movement of the plunger-type wavemaker is necessary. The mechanism of plunger-type solitary wave generation is the real-time equivalence between the immersion volume of the plunger and the volume of the wave. Since a solitary wave lays entirely above the SWL, its volume increases monotonically during its generation and the movement of the plunger should always be downward. A faster or slower descent, however, both results in inappropriate water volume and downgrades the quality of the generated wave. Thus, the primary objective of this study is to propose a new theory that can be used to precisely control the plunger descent so as to produce precise target solitary waves.
Unlike the piston-type wavemaker which has a simple geometry and can displace an arbitrary quantity of water by adjusting the piston stroke, the complex plunger-type wavemaker may encounter an incompatibility problem at the interface between plunger and wave (see Figs. 3 and 6). It also displaces limited water quantity due to the finite water depth or the finite volume of the plunger (see Fig. 2). Thus, to obtain the constraints on the generated wave and understand the consequences of exceeding the constraints are the secondary objective of this work.
There are two steps in validating the proposed theory: first, to generate solitary waves based on the proposed theory; second, to examine the accuracy of the generated wave. As for the first step, conducting a physical experiment is a reliable approach, however, it depends on the financial resources, laboratory space, facility availability, etc. Then, numerical simulation can be used alternatively due to its increased accuracy and popularity. The Smoothed Particle Hydrodynamics (SPH) [39], [40], [41], [42] is a Lagrangian meshfree numerical method, which has been successfully applied in various fields of coastal and ocean engineering, including wave/current-structure interaction [43], [44], [45], [46], wave propagation over topography [47], [48], [49], [50], [51], fluid resonance [52, 53], renewable energy utilization [54], [55], [56], sediment dynamics [57, 58], and many others as mentioned by Gotoh and Khayyer [59, 60]. The SPH method is adopted here to simulate the plunger water entry and solitary wave generation because it is naturally good at simulating highly nonlinear wave motion, its Lagrangian characteristic facilitates the treatment of the water-plunger interface, and the latest SPH method is able to capture the vortex dynamics [61, 62, 49] in the water entry process of a plunger with sharp corners. Over the past two decades, there have already been many SPH works on the water entry and the consequent wave generation. For example, Monaghan and Kos [63] and Monaghan et al. [64] simulated the vertical sinking and inclined slide of a box respectively, and established relations between the water entry speed and the wave height. Lin et al. [65] modelled the entire process of the wedge-shaped landslide, surge wave generation, propagation, and overtopping of a dam. By analysing subaerial landslides into various water body geometries, Heller et al. [66] provided a novel insight into the slide and wave kinematics. Capone et al. [67] and Xenakis et al. [6] treated sliding masses as non-Newtonian fluid so as to investigate the deformable landslide process. These representative studies are all circumstantial evidence that the SPH method is applicable to the present study.
Based on the SPH method, several open-source packages have been launched such as the classic SPHysics [68] and the advanced DualSPHysics [69]. SPHysics is a FORTRAN code that includes both serial and parallel versions and both 2D and 3D versions. It is easy to read and modify but is generally applicable to small and low-resolution domains. The code we use here is developed from the 2D serial SPHysics of version 2.2 but incorporating several beneficial functionalities, for instance, the OpenMP parallelization, modified dynamic boundary condition [70], δ-SPH method [71, 72], and particle shifting technique [73]. DualSPHysics is a set of C++, CUDA, and Java codes that can be executed either on a CPU or a GPU. It is characterized by high computational efficiency, which makes large-domain and high-resolution simulations ordinary. Moreover, it constantly incorporates the latest SPH achievements. Owing to the above features, DualSPHysics has been attracting a rapidly growing number of users.
All the simulations in this work are performed in 2D space based on the following considerations. In the laboratory, the width of a plunger-type wavemaker is generally close to that of the water tank. Thus, fluid motion is dominated by the wave propagation parallel to the sidewalls of the tank. There can be some fluid separation and wave diffraction due to the narrow gaps between the plunger and side walls or among multi-segment plungers, however, they are relatively weaker than the dominated wave motion. In addition, although 3D vortex motion can be triggered during the fluid-plunger interaction, their effects on the wave motion should also be small.
As for the second step in validation of the proposed theory, namely, to examine the accuracy of the generated wave, the SPH-simulated wave profile, velocity field, and pressure field can be examined by either comparing with experimental data or analytical solutions. Herein, the latter is preferred mainly due to three reasons. First, existing experiments usually focus on partial wave information [74], [75], [76], while analytical solutions are more comprehensive. Second, existing experiments only involve a limited number of wave conditions and are unlikely to cover all the conditions of this study, while analytical solutions are unlimited, provided their prerequisites are satisfied. Third, solitary wave theory has been almost mature. The experimental data are not necessarily more reliable than the analytical solutions if considering the instrument and measurement errors.
This paper is organized in the following manner. After the introduction, Sec. 2 derives the generic theory of plunger-type solitary wave generation and the generic constraints on the generated wave. Sec. 3 introduces the SPH method. In Sec. 4 to Sec. 6, three representative plungers, namely wedge, box, and cylinder are tested, respectively. In each section, formulae for the plunger descent and constraints on the produced wave height are first specified. Then, the reliability of the SPH model is validated by reproducing an available experiment. Finally, the validated SPH model is used to simulate the solitary wave generation based on the derived formulae, and physical quantities of the generated wave are examined by comparing SPH results with the analytical solutions. In Sec. 7, plunger-type and piston-type solitary wave generation and Rayleigh-based and Boussinesq-based solitary wave generation are compared. The main conclusions are drawn in Sec. 8, together with our perspectives on future research.
Section snippets
Wave-making theory and constraints
This section innovatively derives the generic formulae for the plunger-type solitary-wave generation. In Subsec. 2.1, theory of solitary wave generation regarding an arbitrary-geometry plunger will be proposed. In Subsec. 2.2, constraints on the generated waves will be investigated.
The SPH model
SPH method can be generally categorized into weakly compressible and incompressible modes. The weakly compressible SPH (WCSPH) assumes a variable density field and solves the fluid motion explicitly by establishing a relation between pressure and density [79, 80]. Thus, it is easy to program and is especially suitable for parallel computation [69]. However, WCSPH is characterized by unphysical pressure fluctuations, so that techniques such as the artificial viscosity [81], density
Case I:Wedge-shaped plunger
Taking a wedge-shaped plunger as the first example, specific formulae for the displacement of the plunger and constraints on the produced wave height will be given in Subsec. 4.1. Then, the SPH model will be validated in Subsec. 4.2 by reproducing the wedge-induced regular waves reported in the literature. In Subsec. 4.3, the validated SPH model will be applied to simulate the wedge-induced solitary waves by using the new formulae derived in Subsec. 4.1 and various aspects of the generated
Case II:Box-shaped plunger
Taking a box-shaped plunger as the second example, specific formulae for the displacement of the plunger and constraints on the produced wave height will be given in Subsec. 5.1. Then, the SPH model will be validated in Subsec. 5.2 by reproducing the box-induced impulse waves reported in the literature. In Subsec. 5.3, the validated SPH model will be applied to simulate the box-induced solitary waves by using the new formulae derived in Subsec. 5.1 and various aspects of the generated waves
Case III:Cylinder-shaped plunger
Taking a cylinder-shaped plunger as the third example, specific formulae for the displacement of the plunger and constraints on the produced wave height will be given in Subsec. 6.1. Then, the SPH model will be validated in Subsec. 6.2 by reproducing the cylinder-induced regular waves reported in the literature. In Subsec. 6.3, the validated SPH model will be applied to simulate the cylinder-induced solitary waves by using the new formulae derived in Subsec. 6.1 and various aspects of the
Two comparisons
Through cases I, II, and III, it has been confirmed that high-quality solitary waves can be generated using plunger-type wavemakers based on the proposed wave-making theory, as long as the constrains on the produced wave height are satisfied. In Subsec. 7.1, a further comparison will be made between plunger-type and piston-type solitary wave generation. In addition, since Boussinesq assumption is adopted throughout Secs. 4, 5, and 6, plunger-induce solitary waves under Rayleigh assumption will
Conclusions
Solitary wave generation ought to be a fundamental functionality of a wavemaker. However, despite the worldwide existence of plunger-type wavemakers, they are incompetent in generating accurate solitary waves because of lack of theory. This work solved this problem by deriving new formulae for plunger descent and constraint conditions on produced wave height. A fully validated SPH model was then used to examine the quality of the generated wave by taking wedge-shaped, box-shaped, and
CRediT authorship contribution statement
Ming He: Conceptualization, Investigation, Validation, Visualization, Writing - original draft. Abbas Khayyer: Writing - review & editing, Methodology, Conceptualization. Xifeng Gao: Funding acquisition, Project administration. Wanhai Xu: Supervision. Bijin Liu: Resources.
Declaration of Competing Interest
The authors declare that this manuscript has not been previously published in whole or in part nor is it under consideration for publication elsewhere. The authors have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. The publication of this work is approved by all authors and explicitly by the responsible authorities where the work was carried out.
Acknowledgements
This work was supported by the National Natural Science Foundation of China [grant number 51709201] and the 2019 Subsidy Fund Project for Marine Economy Development in Fujian Province [grant number FJHJF-L-20198]. The authors are grateful to three anonymous reviewers for their careful reviews and insightful suggestions that have enhanced the quality and clarity of the presented work.
References (111)
- et al.
Numerical investigation of transient harbor oscillations induced by N-waves
Coast. Eng.
(2017) - et al.
Topographic influences on transient harbor oscillations excited by N-waves
Ocean Eng.
(2019) - et al.
Numerical study on transient harbor oscillations induced by solitary waves
Ocean Eng.
(2016) - et al.
Effects of offshore fringing reefs on the transient harbor resonance excited by solitary waves
Ocean Eng.
(2019) - et al.
Interaction of solitary waves with emergent, rigid vegetation
Ocean Eng.
(2011) - et al.
Boussinesq modelling of solitary wave and N-wave runup on coast
Appl. Ocean Res.
(2013) - et al.
Reflection and transmission of tsunami waves by coastal structures
Appl. Ocean Res.
(2000) - et al.
Comparative study on the accuracy of solitary wave generations in an ISPH-based numerical wave flume
Appl. Ocean Res.
(2016) - et al.
Interaction between a Stokes wave packet and a solitary wave
Eur. J. Mech. B-Fluids
(1999) - et al.
Numerical simulation of nonlinear wave and force generated by a wedge-shape wave maker
Ocean Eng.
(2006)
Smoothed particle hydrodynamics and its applications in fluid-structure interactions
J. Hydrodyn.
Numerical simulation of a damaged ship cabin flooding in transversal waves with Smoothed Particle Hydrodynamics method
Ocean Eng.
Improved SPH simulation of spilled oil contained by flexible floating boom under wave–current coupling condition
J. Fluids Struct.
An improved solid boundary treatment for wave–float interactions using ISPH method
Int. J. Nav. Archit. Ocean Eng.
Consistent Particle Method simulation of solitary wave impinging on and overtopping a seawall
Eng. Anal. Bound. Elem.
An ISPH with k–ε closure for simulating turbulence under solitary waves
Coast. Eng.
Numerical analysis of wave-induced current within the inhomogeneous coral reef using a refined SPH model
Coast. Eng.
A set of canonical problems in sloshing, Part I: Pressure field in forced roll—comparison between experimental results and SPH
Ocean Eng.
Towards simulating floating offshore oscillating water column converters with Smoothed Particle Hydrodynamics
Coast. Eng.
Potential application of submerged horizontal plate as a wave energy breakwater: A 2D study using the WCSPH method
Ocean Eng.
An improved SPH model for turbulent hydrodynamics of a 2D oscillating water chamber
Ocean Eng.
Development of a two-phase SPH model for sediment laden flows
Comput. Phys. Commun.
An SPH study of driven turbulence near a free surface in a tank under gravity
Eur. J. Mech. B-Fluids
The δplus-SPH model: Simple procedures for a further improvement of the SPH scheme
Comput. Meth. Appl. Mech. Eng.
Composite modelling of subaerial landslide–tsunamis in different water body geometries and novel insight into slide and wave kinematics
Coast. Eng.
SPHysics – development of a free-surface fluid solver – Part 1: Theory and formulations
Comput. Geosci.
DualSPHysics: Open-source parallel CFD solver based on Smoothed Particle Hydrodynamics (SPH)
Comput. Phys. Commun.
Nonlinear simulations of wave-induced motions of a freely floating body using WCSPH method
Appl. Ocean Res.
Numerical diffusive terms in weakly-compressible SPH schemes
Comput. Phys. Commun.
Free-surface flows solved by means of SPH schemes with numerical diffusive terms
Comput. Phys. Commun.
Comparative study on accuracy and conservation properties of two particle regularization schemes and proposal of an optimized particle shifting scheme in ISPH context
J. Comput. Phys.
A simple procedure to improve the pressure evaluation in hydrodynamic context using the SPH
Comput. Phys. Commun.
Simulating free surface flows with SPH
J. Comput. Phys.
Shock simulation by the particle method SPH
J. Comput. Phys.
Numerical simulation of interfacial flows by smoothed particle hydrodynamics
J. Comput. Phys.
Numerical modeling of water waves with the SPH method
Coast. Eng.
Modified Moving Particle Semi-implicit methods for the prediction of 2D wave impact pressure
Coast. Eng.
Corrected incompressible SPH method for accurate water-surface tracking in breaking waves
Coast. Eng.
A higher order Laplacian model for enhancement and stabilization of pressure calculation by the MPS method
Appl. Ocean Res.
Enhancement of stability and accuracy of the moving particle semi-implicit method
J. Comput. Phys.
Smoothed particle hydrodynamics: Some recent improvements and applications
Comput. Meth. Appl. Mech. Eng.
Optimal sponge layer for water waves numerical models
Ocean Eng.
A measure of spatial disorder in particle methods
Comput. Phys. Commun.
Multi-resolution Delta-plus-SPH with tensile instability control: Towards high Reynolds number flows
Comput. Phys. Commun.
Incompressible smoothed particle hydrodynamics for free-surface flows: A generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves
J. Comput. Phys.
Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach
J. Comput. Phys.
Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method
J. Comput. Phys.
Survey of 2011 Tohoku earthquake tsunami inundation and run-up
Geophys. Res. Lett.
The September 16, 2015 Mw8.3 Illapel, Chile Earthquake: characteristics of tsunami wave from near-field to far-field
Acta Oceanol. Sin.
Source mechanisms of volcanic tsunamis
Philos. Trans. R. Soc. A-Math. Phys. Eng. Sci.
Cited by (42)
A robust δ-SPH<sup>C</sup> model for nonlinear water wave interactions with structures under complex wave conditions
2024, Engineering Analysis with Boundary ElementsExperimental validation of the plunger-type flow model for irregular waves
2023, Ocean EngineeringA novel Eulerian SPH shallow water model for 2D overland flow simulations
2023, Journal of HydrologyEnhanced resolution of the continuity equation in explicit weakly compressible SPH simulations of incompressible free-surface fluid flows
2023, Applied Mathematical ModellingCitation Excerpt :To resolve these issues, the incorporation of the so-called Particle Shifting (PS) scheme [17] that regularizes local particle distributions, has been gaining continued popularity. The δ-SPH method with the PS scheme has been applied to a wide range engineering applications and has been shown to provide enhanced robustness and accuracy (e.g., [18–23]). Despite enhanced stability and accuracy of the δ-SPH method with the PS scheme, there still exists a crucial drawback in the simulation of incompressible free-surface fluid flows; namely, the imprecise numerical resolution of the continuity equation resulting from the departure from incompressibility conditions corresponding to divergence-free velocity and invariant density field, especially in long-term simulations.