Next Article in Journal
Phototactic Behavioral Responses of Mesozooplankton in the Barents Sea as an Indicator of Anthropogenic Impact
Next Article in Special Issue
Size Distribution and Variation in Surface-Suspended Sediments in the North Passage Estuarine Turbidity Maximum of the Yangtze Estuary
Previous Article in Journal
Hydrogeochemical Characteristics and Evolution Processes of Karst Groundwater Affected by Multiple Influencing Factors in a Karst Spring Basin, Eastern China
Previous Article in Special Issue
Water Renewal Simulation in Two Flow-Through Water Bodies in Western Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Order Boussinesq Equations for Water Wave Propagation in Porous Media

1
State Environmental Protection Key Laboratory of Marine Ecosystem Restoration, National Marine Environmental Monitoring Center, Dalian 116023, China
2
College of Transportation Engineering, Dalian Maritime University, Dalian 116026, China
3
State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China
4
Liaoning Transportation Affairs Service Center, The Transportation Department of Liaoning Province, Shenyang 110003, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(22), 3900; https://doi.org/10.3390/w15223900
Submission received: 18 September 2023 / Revised: 2 November 2023 / Accepted: 3 November 2023 / Published: 8 November 2023
(This article belongs to the Special Issue Hydrodynamics in Coastal Areas)

Abstract

:
To accurately capture wave dynamics in porous media, the higher-order Boussinesq-type equations for wave propagation in deep water are derived in this paper. Starting with the Laplace equations combined with the linear and nonlinear resistance force of the dynamic conditions on the free surface, the governing equations were formulated using various independent velocity variables, such as the depth-averaged velocity and the velocity at the still water level and at an arbitrary vertical position in the water column. The derived equations were then improved, and theoretical analyses were carried out to investigate the linear performances with respect to phase celerity and damping rate. It is shown that Boussinesq-type models with Padé [4, 4] dispersion can be applied in deep water. A numerical implementation for one-dimensional equations expressed with free surface elevation and depth-averaged velocity is presented. Solitary wave propagation in porous media was simulated, and the computed results were found to be generally in good agreement with the measurements.

1. Introduction

Human-made porous structures, such as rock breakwaters and revetments, are often built close to the shore to protect coastal buildings and prevent severe coastal erosion. These types of coastal structures partially reflect wave energy and partially dissipate wave energy due to the drag and inertial resistance of the structures, resulting in only partial wave transmission. Therefore, it is important to have a good understanding of how the wave motion interacts with porous structures. Much research has been carried out since 1990, and an excellent review was conducted by Losada et al. [1].
Research on the interaction between waves and porous structures has been successfully carried out using differential approaches such as analytical solutions [2], laboratory measurements [1,3], and numerical modeling [4,5,6,7]. Among these works, numerical models based on the Navier–Stokes (RANS) equations and nonlinear shallow water (NLSW) equations are commonly used to calculate wave interaction with porous structures and appear to be more advantageous. RANS models can represent the details of the flow fields of porous structures, as shown in research conducted by Sasikumar et al. [8], but require longer computation times, making the application of RANS models to the interaction between waves and porous structures very challenging at present. In contrast to RANS models, NLSW models are much more computationally efficient as only depth-averaged parameters are calculated, as they can efficiently simulate long waves in a shallow water region but are not suitable for accurately describing short, dispersive waves and are inadequate for the interactions between waves and porous structures. Mohamed et al. [9] used NLSW models to simulate shallow waves with porosity structures based on the finite volume method. In recent years, smoothed-particle hydrodynamic (SPH) models have also been used to simulate wave propagation; these are meshless and achieve similar accuracy to RANS-type models, and their validity has been demonstrated by Ren et al. [10], but these models still require a huge computational effort, which limits their application in modeling wave propagation in porous structures.
In recent decades, Boussinesq-type models have been significantly enriched and widely applied in modeling the coastal waves and currents in complex bathymetric configurations, such as the FUNWAVE model [11,12,13] and the MIKE21 BW model [14]. The Boussinesq-type models have better computational efficiency, with their main advantage being the simplification of the three-dimensional problem into the two-dimensional one, as outlined in excellent reviews by Kirby [15], Madsen and Fuhrman [16], Brocchini [17], Kirby [18], and Sun et al. [19]. As stated by Sun et al. [19], the Boussinesq-type models can be classified into two groups: the quasi-three-dimensional (3D) Boussinesq-type models and the two-horizontal-dimension (2HD) models. Sun et al. [19] also summarized the use of Boussinesq-type models for wave propagation over a porous seabed and in multi-layer water with different densities. The use of Boussinesq-type models for wave propagation and evolution on the impermeable seabed [20,21] has been fully developed, and some representative types of these equations include the extended Boussinesq-type models, with improved linear dispersion/shoaling characteristics [22,23,24,25]; the second-order nonlinear Boussinesq-type models [26,27]; the so-called fully nonlinear Boussinesq-type models [28,29,30,31,32]; and the extremely dispersive and nonlinear Boussinesq-type models [33,34,35,36]. More recently, using the methods proposed in the literature [30,33,35,37], a Boussinesq-type model for interfacial waves over a two-density system has been proposed.
The above Boussinesq-type models with exact linear and nonlinear properties cannot be used directly for wave propagation in porous media. As can be seen, the inclusion of porosity in Boussinesq-type models has rarely been reported. Cruz et al. [38] first proposed a set of weakly nonlinear Boussinesq-type models using depth-averaged velocities in both free and permeable water surfaces to model wave propagation over a permeable seabed. Starting with the Euler equations, Hsiao et al. [39] and Chen [40] independently proposed two sets of higher-order Boussinesq-type equations expressed in terms of velocities at two arbitrary vertical locations. Liu and Sun [41] extended the Boussinesq equations of Hsiao et al. [39] to make them applicable to wave propagation in deeper water. More recently, the composite Boussinesq-type models for use over porous beds were proposed by Klonaris and Memos [42]. However, when the porous structures were sub-aerial with surface-piercing boundaries, the aforementioned Boussinesq-type models could not be applied. For such cases, Hsiao et al. [43] were the first to develop a second-order Boussinesq model for starting to simulate similar conditions, and more recently, Fang et al. [44] developed a Boussinesq-type wave model to simulate the interaction of coastal waves with bottom-mounted porous structures. However, these models (e.g., [43,44]) may not be valid for short-wave propagation from deeper water into shallower water.
To overcome these shortcomings and to fully account for wave dispersion, nonlinearity, and modeling efficiency, it is essential to develop Boussinesq-type models suitable for cases of porous structures that span from deeper water to shallower water. The aim of this study was to develop a Boussinesq-type model that is both valid for wave propagation in deeper water and suitable for wave interactions with a surface-piercing porous structure. By analyzing the above-mentioned Boussinesq-type models and comparing them with the models for waves in impermeable seabed cases, the Boussinesq-type models valid for wave interaction with a surface-penetrating porous structure have the linear and nonlinear resistance forces resulting from the water passing through the porous media. Therefore, we will derive the Boussinesq-type models following the main procedures as developed by Madsen and Schäffer [45].
This paper is organized as follows. The detailed derivation of four sets of Boussinesq-type equations for wave propagation in a porous structure is given in Section 2. In Section 3, the phase celerity and damping rate of the Boussinesq models are theoretically analyzed and compared with the analytical solutions. Numerical implementation of 1D equations and simulation of solitary wave interaction with porous breakwater are carried out to validate the present model in Section 4. Finally, conclusions are drawn in Section 5.

2. Derivation of High Order of Boussinesq-Type Equations

2.1. Governing Equations and Boundary Conditions

The sketch of wave propagation in a porous media is shown in Figure 1, where the computational domain is divided into non-porous media and porous media. In Figure 1, z = 0 represents the still water level, η x , y , t is the free surface elevation (measured at still water level), h b ( x , y ) is the depth of the porous media, and ψ x , y , z , t denotes the velocity potential in non-porous and porous media.
The flow is assumed to be impressible and irrotational, and the equation of motion inside and outside the porous media is derived by Cruz et al. [38] and can be written as:
n d U s d t + 1 ρ 3 p s + ρ g z + F r + F i = 0
where n is the porosity, n = 1 represents wave propagation in non-porous media and n < 1 represents porous media. U s ( u s , v s , w s ) is the three-dimensional velocity vector, p s is the pore pressure, ρ is the fluid density, g is the gravity acceleration, 3 = ( / x ,   / y ,   / z ) is the gradient operator, and d / d t = / t + U s 3 denotes the total derivative. F r and F i denote drag resistance and inertial force, respectively, and they are defined as
F r = α 1 U s + α 2 U s U s , F i = ( 1 n ) ( 1 + c m ) d U s d t
where α 1 and α 2 denote coefficients for laminar fluid and turbulence force, respectively, and c m is the added mass coefficient. Inserting Equation (2) into Equation (1), we can obtain the following equation
c r d U s d t + 1 ρ 3 ( p s + ρ g z ) + α 1 U s + α 2 U s U s = 0
where c r = n + ( 1 n ) ( 1 + c m / n ) is the inertial force coefficient.
The mass conservation equation in porous media is
3 ( n U s ) = 0
If the porosity is assumed to be a constant value, Equation (4) can be written as
3 U s = 0
At the free surface, the dynamic and kinematic conditions are
p s = 0 z = η ( x , y , t )
d d t ( z η ) = w s η t u s η = 0 z = η ( x , y , t )
where = x , y is the horizontal gradient operator and u s = ( u s , v s ) is the corresponding velocity vector at the free surface. At the bottom boundary, we have
u s h b + w s = 0 z = h b
To reduce the number of unknown variables in the derivation, the three-dimensional velocity vector can be replaced by the velocity potential, as shown in the following formula
U s 3 ψ
The mass conservation equation is thus transformed into the Laplace equations
2 ψ + ψ z z = 0 , h b < z < η
and Equation (3) can be written as
c r ( ψ t + 1 2 3 ψ 2 ) + p s + g z + α ψ = 0
According to Equation (11), the dynamic and kinematic conditions at the free surface can be written as follows
c r ( ψ t + 1 2 3 ψ 2 ) + g η + α ψ = 0 z = η ( x , y , t )
ψ z = η t ψ η = 0 z = η ( x , y , t )
And at the bottom, the boundary condition is
ψ h b + ψ z = 0 z = h b
Thus, Equations (10) and (12)–(14) are the motion equations and boundary conditions inside the porous media expressed by velocity potential.

2.2. Dimensionless Equations

To obtain the dimensionless equations, we normalize the variables using characteristic wave amplitude a, characteristic water depth h0, and characteristic wavelength l0 as follows
x = x l 0 , y = y l 0 , z = z h 0 , h b = h b h 0 , η = η a , t = t g h 0 l 0 , ψ = ψ a g h 0 l 0 h 0
Two parameters representing nonlinearity and dispersion are introduced in the derivation as
ε = a h 0 , μ = h 0 l 0
Omitting primes for simplicity, the Equations (10) and (12)–(14) can be written as
μ 2 2 ψ + ψ z z = 0 , h b < z < ε η
μ 2 ( c r ψ t + α ψ + η ) + ε c r 1 2 [ μ 2 ( ψ ) 2 + ψ z 2 ] = 0 z = ε η
μ 2 ( η t + ε ψ η ) = ψ z z = ε η
μ 2 ψ h b + ψ z = 0 , z = h b
Notice that by setting cr = 1 and α = 0, Equation (18) recovers to the expressions presented in the original paper of Madsen and Schäffer [45].

2.3. Power Series Solution to the Laplace Equations

The main advantage of the Boussinesq equations is to simplify the three-dimensional problem into a two-dimensional one. We first expand the velocity potential as
ψ ( x , y , z , t ) = n = 0 z n ψ ( n ) ( x , y , t )
The first and second order of the potential derivatives can be expressed as
ψ ( x , y , z , t ) = n = 0 z n ψ ( n ) ( x , y , t )
2 ψ ( x , y , z , t ) = n = 0 z n 2 ψ ( n ) ( x , y , t )
ψ z ( x , y , z , t ) = n = 0 ( n + 1 ) z n ψ ( n + 1 ) ( x , y , t )
ψ z z ( x , y , z , t ) = n = 0 ( n + 2 ) ( n + 1 ) z n ψ ( n + 2 ) ( x , y , t )
Inserting Equations (23) and (24) into Equation (17), we obtain a recurrence relation as
ψ ( n + 2 ) = 1 ( n + 1 ) ( n + 2 ) μ 2 2 ψ ( n ) n = 0 , 1 , 2 ,
The velocity potential is thus expressed as follows
ψ ( x , y , z , t ) = n = 0 ( 1 ) n μ 2 n z 2 n ( 2 n ) ! 2 n ψ ( 0 ) + z 2 n + 1 ( 2 n + 1 ) ! 2 n ψ ( 1 )
Equation (27) indicates that the velocity potential can be expressed by the expansion ψ ( 0 ) and ψ ( 1 ) . The derivatives of velocity potential over horizontal plane (x, y) and vertical plane z direction can be written as
u s ( x , y , z , t ) = n = 0 ( 1 ) n μ 2 n z 2 n ( 2 n ) ! ( 2 n 2 ( u s ) ) + z 2 n + 1 ( 2 n + 1 ) ! μ 2 ( 2 n w s )
w s ( x , y , z , t ) = n = 0 ( 1 ) n μ 2 n + 2 z 2 n + 1 ( 2 n + 1 ) ! 2 n ( u s ) + z 2 n ( 2 n ) ! 2 n w s
where u s ψ ( 0 ) , w s ( 1 / μ 2 ) ψ ( 1 ) , i.e.,
u s ( x , y , 0 , t ) = u s , w s ( x , y , 0 , t ) = μ 2 w s
Inserting Expression (29) into Equation (20), we obtain
w s + n = 0 ( 1 ) n μ 2 n h b 2 n + 1 ( 2 n + 1 ) ! ( 2 n 2 ( u s ) ) ) μ 2 h b 2 n + 2 ( 2 n + 2 ) ! ( 2 n w s ) = 0
Similarly, w s x , y , t = m = 0 μ 2 n w s n x , y , t is introduced and truncated to the order of O ( μ 6 ) , and we have
w s ( x , y , t ) = ( h b u s ) + μ 2 { 1 6 h b 3 ( u s ) 1 2 h b 2 [ ( h b u s ) ] } + μ 4 { 1 24 h b 4 [ 2 ( ( h b u s ) ) ] 1 120 h b 5 [ 2 ( u s ) ] + 1 2 h b 2 [ ( 1 6 h b 3 ( u s ) 1 2 h b 2 ( ( h b u s ) ) ) ] } + O ( μ 6 )

2.4. Boussinesq Equations in Terms of Velocity at the Free Surface u ^ s

Inserting Equations (27) and (28) into Equation (18), applying the horizontal gradient operator, and eliminating w s ( x , y , t ) using Expression (32) under the assumption that ε = O μ , we have
u ^ s t + α / c r u ^ s + η / c r + ε 1 2 ( u ^ s ) 2 + μ 2 { ε T 21 + ε 2 T 22 + ε 3 T 23 } + μ 4 ( ε T 41 ) = O ( μ 6 , ε 2 μ 4 )
T 21 = [ η ( h b u ^ s t ) α / c r η ( h b u ^ s ) + 1 2 ( ( h b u ^ s ) ) 2 ]
T 22 = [ 1 2 η 2 u ^ s t 1 2 α / c r η 2 u ^ s η u ^ s ( h b u ^ s ) + η ( h b u ^ s ) u ^ s ]
T 23 = [ 1 2 η 2 u ^ s u ^ s + 1 2 ( η u ^ s ) 2 ]
T 41 = [ η ( h b 2 Γ t ) + α / c r η ( h b 2 Γ ) ( h b 2 Γ ) ( h b u s ) ]
Γ = 1 6 h b ( u s ) 1 2 ( ( h b u s ) )
Similarly, using the same simplification method, we can obtain the depth-integrated continuity equations
η t + Q = 0
where
Q h b ε η ψ d z = u ^ s h b + ε η μ 2 { 1 6 ( ε 3 η 3 + h b 3 ) ( u ^ s ) + 1 2 ( ε 2 η 2 h b 2 ) [ ( h b u ^ s ) ] } + μ 4 { 1 2 ( ε 2 η 2 h b 2 ) [ ( h b 2 Γ ) ] + 1 24 ( ε 4 η 4 h b 4 ) [ 2 ( ( h b u ^ s ) ) ] + 1 120 ( ε 5 η 5 + h b 5 ) [ 2 ( u ^ s ) ] } + O ( μ 6 )

2.5. Boussinesq Equations in Terms of Depth-Averaged Velocity u ¯ s

The depth-averaged velocity u ¯ s is defined as
u ¯ s = 1 h b + ε η h b ε η ψ d z
and it can be expressed in terms of u ^ s as
u ¯ s = u ^ s + μ 2 { 1 2 h b ε η [ h b u ^ s ] 1 6 h b 2 ε h b η + ε 2 η 2 u ^ s } + μ 4 { 1 2 ( h b ε η ) [ ( h b 2 Γ ) ] 1 24 ( h b 3 ε h b 2 η + ε 2 η 2 h b ε 3 η 3 ) [ 2 ( ( h b u ^ s ) ) ] + 1 120 ( h b 4 ε h b 3 η + ε 2 η 2 h b 2 ε 3 η 3 h b + ε 4 η 4 ) [ 2 ( u ^ s ) ] } + O ( μ 6 )
From the expression, we can obtain
u ^ s = u ¯ s μ 2 { 1 2 ( h b ε η ) [ ( h b u ¯ s ) ] 1 6 ( h b 2 ε h b η + ε 2 η 2 ) ( u ¯ s ) } + μ 4 { 1 24 ( h b 3 ε h b 2 η + ε 2 η 2 h b ε 3 η 3 ) [ 2 ( ( h b u ¯ s ) ) ] 1 120 ( h b 4 ε h b 3 η + ε 2 η 2 h b 2 ε 3 η 3 h b + ε 4 η 4 ) [ 2 ( u ¯ s ) ] 1 2 ( h b ε η ) [ ( 1 2 ε η h b ( ( h b u ¯ s ) ) 1 6 ( ε η h b 2 ε 2 η 2 h b ) ( u ¯ s ) ) ] 1 6 ( h b 2 ε h b η + ε 2 η 2 ) [ ( 1 2 ( h b ε η ) ( ( h b u ¯ s ) ) 1 6 ( h b 2 ε h b η + ε 2 η 2 ) ( u ¯ s ) ) ] } + O ( μ 6 )
The continuity equation in terms of u ¯ s is
η t + [ ( h b + ε η ) u ¯ s ] = 0
Inserting Expression (39) into Equations (33) and (34), and introducing the mild-slope assumption n h b = O μ , n = 1 , 2 , 3 , we obtain
u ¯ s t + α / c r u ¯ s + η / c r + ε 1 2 ( u ¯ s ) 2 + μ 2 { T M 20 + ε T M 21 + ε 2 T M 22 + ε 3 T M 23 } + μ 4 ( T M 40 + ε T M 41 ) = O ( μ 6 , ε 2 μ 4 )
T M 20 = h b Γ t + α / c r h b Γ
T M 21 = η Γ t α / c r η Γ + [ u ¯ s ( h b Γ ) η ( h b u ¯ s t ) α / c r η ( h b u ¯ s ) + 1 2 ( ( h b u ¯ s ) ) 2 ]
T M 22 = 1 6 η 2 ( u ¯ s t ) + α / c r η 2 ( u ¯ s ) 1 3 η ( h b u ¯ s ) ( u ¯ s ) + ( η u ¯ s ) Γ + [ η u ¯ s Γ 1 2 η 2 u ¯ s t 1 2 α / c r η 2 u ¯ s η u ¯ s ( h b u ¯ s ) + η ( h b u ¯ s ) u ¯ s ]
T M 23 = 1 3 η ( η u ¯ s ) ( u ¯ s ) + [ 1 3 η 2 u ¯ s u ¯ s + 1 2 ( η u ¯ s ) 2 ]
T M 40 = 1 24 h b 3 ( 2 ( ( h b u ¯ s t ) ) ) 1 120 h b 4 ( 2 ( u ¯ s t ) ) + 1 6 h b 2 ( ( h b Γ t ) ) + α / c r [ 1 24 h b 3 ( 2 ( ( h b u ¯ s ) ) ) 1 120 h b 4 ( 2 ( u ¯ s ) ) + 1 6 h b 2 ( ( h b Γ ) ) ]
T M 41 = 1 45 h b 3 η ( 2 ( u ¯ s t ) ) 1 9 h b 3 ( ( η ( u ¯ s t ) ) ) + α / c r [ 1 45 h b 3 η ( 2 ( u ¯ s ) ) 1 9 h b 3 ( ( η ( u ¯ s ) ) ) ] 1 45 h b 4 u ¯ s ( 2 ( u ¯ s ) ) + 1 9 h b 4 [ ( u ¯ s ( ( u ¯ s ) ) ) ] + 1 18 h b 4 ( ( u ¯ s ) ) 2
Γ = 1 6 h b ( u ¯ s ) 1 2 ( ( h b u ¯ s ) )

2.6. Boussinesq Equations in Terms of Velocity at an Arbitrary Water Column u s a

Keeping the order accurate up to O μ 4 , the velocity at an arbitrary water depth can be written as
ψ ( x , y , z α h b , t ) = ψ ( 0 ) μ 2 z ( h b ψ ( 0 ) ) 1 2 μ 2 z 2 h b 2 2 ψ ( 0 ) + O ( μ 4 )
The expression of u s α in terms of u ^ s is given
u s α = u ^ s μ 2 { z α h b ( h b u ^ s ) + 1 2 z α 2 h b 2 u ^ s }
So, we can obtain the following expression
u ^ s = u s α + μ 2 { z α h b ( h b u s α ) + 1 2 z α 2 h b 2 u s α }
Inserting (45) into momentum Equation (33), we have
u s a t + α / c r u s a + η / c r + ε 1 2 ( u s a ) 2 + μ 2 T s 20 + ε T s 21 + ε 2 T s 22 + ε 3 T s 23 = O ( μ 4 )
T s 20 = Γ s a t + α / c r Γ s a
T s 21 = [ η ( h b u s α t ) α / c r η ( h b u s α ) + 1 2 ( ( h b u s α ) ) 2 + u s α Γ s α ]
T s 22 = [ 1 2 η 2 u s α t 1 2 α / c r η 2 u s α η u s α ( h b u s α ) + η ( h b u s α ) u s α ]
T s 23 = [ 1 2 η 2 u s α u s α + 1 2 ( η u s α ) 2 ]
Γ s α = z α h b ( h b u s α ) + 1 2 z α 2 h b 2 u s α
Inserting (45) into continuity Equation (35), we obtain
η t + [ ( h b + ε η ) u s a + μ 2 [ ( 1 2 h b 2 + z α h b 2 + ε z α η h b 1 2 ε 2 η 2 ) ( ( h b u s a ) ) ( 1 2 z α 2 h b 3 1 6 h b 3 + ε z α 2 η h b 2 1 6 ε 3 η 3 ) ( u s a ) ] = 0

2.7. Boussinesq Models Extended to Deeper Water Depth

In order to extend the applicability of the models to deeper water, the technique presented by Madsen and Schäffer [45] is used to improve the derived model. This method has been widely used in the literature to obtain the improved Boussinesq-type equations for open water problems [45,46].

2.7.1. The Improvement of Equations in Terms of Depth-Averaged Velocity

(1)
Keeping at order O ( μ 2 )
Keeping at order O ( μ 2 ) , the depth-averaged velocity Equations (40)–(42) have identical dispersion as the classical Boussinesq equations [47] and are only valid in the shallow water region. To extend its application range, we introduce the following Expression [48].
T e x t a 1 = β 2 ( h b + ε η ) μ 2 ( ( h b + ε η ) ( u ¯ s t + η / c r + α / c r u ¯ s + ε ( u ¯ s 2 ) / 2 ) ) + β 1 ( h b + ε η ) 2 μ 2 ( u ¯ s t + η / c r + α / c r u ¯ s + ε ( u ¯ s 2 ) / 2 ) = O ( μ 4 )
where β 1 and β 2 are dispersive parameters.
By adding Equation (49) to Equation (41), the resulting equation is written as
u ¯ s t + α / c r u ¯ s + η / c r + ε 1 2 ( u ¯ s ) 2 + μ 2 { T M 20 + ε T M 21 + ε 2 T M 22 + ε 3 T M 23 } + T e x t r a 1 = O ( μ 4 )
(2)
Keeping at order O ( μ 4 )
Following Madsen and Schäffer [45], we add the following terms into Equation (41)
T e x t a 2 = ( γ 2 γ 1 ) h b 2 μ 2 [ ( u ¯ s t + η / c r + α / c r u ¯ s + ε ( u ¯ s 2 ) / 2 ) + μ 2 ( T M 20 ) + ε μ 2 ( T M 21 ) ] γ 2 h b μ 2 { [ h b ( u ¯ s t + η / c r + α / c r u ¯ s + ε ( u ¯ s 2 ) / 2 ) ] + μ 2 ( h b T M 20 ) + ε μ 2 ( h b T M 21 ) } + γ 3 h b 4 μ 4 { [ 2 ( ( u ¯ s t + η / c r + α / c r u ¯ s + ε ( u ¯ s 2 ) / 2 ) ) ] } + γ 4 h b 3 h b μ 4 [ 2 ( ( u ¯ s t + η / c r + α / c r u ¯ s ) ] = O ( μ 6 , ε 2 μ 4 )
where γ 1 and γ 3 are dispersive parameters, γ 2 and γ 4 are shoaling parameters. The enhanced version of Equation (41) can be written as
u ¯ s t + α / c r u ¯ s + η / c r + ε 1 2 ( u ¯ s ) 2 + μ 2 { T M 20 + ε T M 21 + ε 2 T M 22 + ε 3 T M 23 } + μ 4 ( T M 40 + ε T M 41 ) + T e x t r a 2 = O ( μ 6 , ε 2 μ 4 )

2.7.2. The Improvement of Equations in Terms of u s α

Neglecting high-order terms, Equations (46)–(48) are written as follows
η t + [ ( h + ε η ) u s α ] = O ( μ 2 )
u s α t + ε u s α u s α + η + α / c r u s α = O ( μ 2 )
Also following Madsen and Schäffer [45], we have
T e x t r a 3 = ( δ 2 - δ 1 ) μ 2 { h b 2 [ η t + ( ( h b + ε η ) u s a ) ] } δ 2 μ 2 { h b 2 [ η t + ( ( h b + ε η ) u s a ) } = O ( μ 4 )
T e x t r a 4 = ( δ 4 - δ 3 ) μ 2 ( h b + ε η ) 2 2 [ u s α t + ε u s α u s α + η / c r + α / c r u s α ] δ 4 μ 2 ( h b + ε η ) 2 { ( h b + ε η ) [ u s α t + ε u s α u s α + η / c r + α / c r u s α } = O ( μ 4 )
Adding Expression (55) to Equation (46) and Expression (56) to Equation (48), we obtain
u s α t + α / c r u s α + η / c r + ε 1 2 ( u s α ) 2 + μ 2 { T s 0 + ε T s 1 + ε 2 T s 2 + ε 3 T s 3 } + T e x t r a 4 = O ( μ 4 )
η t + [ ( h b + ε η ) u s a + μ 2 [ ( 1 2 h b 2 + z α h b 2 + ε z α η h b 1 2 ε 2 η 2 ) ( ( h b u s a ) ) ( 1 2 z α 2 h b 3 1 6 h b 3 + ε z α 2 η h b 2 1 6 ε 3 η 3 ) ( u s a ) ] + T e x t r a 3 = O ( μ 4 )
In summary, four sets of higher-order Boussinesq models for wave propagation in a porous structure have been derived, i.e., Model 1, Equations (40) and (50); Model 2, Equations (40) and (52); Model 3, Equations (46) and (48); and Model 4, Equations (57) and (58). By neglecting the drag and inertial force induced by the porosity, the four models presented here recover those of Madsen and Schäffer [45]. Compared with the Boussinesq-type models in the literature, Models 2 and 4 are the new Boussinesq-type models that are not presented in any literature. Model 3 can be seen as a slightly different version of Hsiao et al. [43], while Model 1 can be seen as an alternative choice to the Boussiensq-type model of Hsiao et al. [43]. Meanwhile, Models 1 and 3 have more accuracy in regards to second-order nonlinear terms compared with Fang et al. [44].

3. Dispersive Analysis on a Horizontal Bottom

The dispersive analysis of the four sets of Boussinesq models is performed on a horizontal bed in one dimension, neglecting the nonlinear terms and focusing on the phase velocity and damping rate. Following a similar procedure in Hsiao et al. [43], we can obtain the following dispersion expressions for the different models.
For Model 1, the dispersive expression is
ω 2 ( c r + α 1 i / ω ) = g K h b 1 + ( β 1 + β 2 ) K 2 h b 2 1 + ( 1 / 3 + β 1 + β 2 ) K 2 h b 2
where ω denotes wave frequency and K = k r + i k i is the complex wave number ( i 2 = 1 ).
For Model 2, the dispersive expression is
ω 2 ( c r + α 1 i / ω ) = g K h b 1 + γ 1 K 2 h b 2 + γ 3 K 4 h b 4 1 + ( 1 / 3 + γ 1 ) K 2 h b 2 + ( γ 3 + γ 1 / 3 1 / 45 ) K 4 h b 4
For Model 3, we have the following expression
ω 2 ( c r + α 1 i / ω ) = g K h b 1 ( 1 / 3 + B ) K 2 h b 2 1 B K 2 h b 2
where B = z α + z α 2 / 2 .
For Model 4, the dispersive expression is
ω 2 ( c r + α 1 i / ω ) = g K h b 1 + ( δ 1 + δ 3 B 1 / 3 ) K 2 h b 2 + δ 1 ( δ 3 B 1 / 3 ) K 4 h b 4 1 + ( δ 1 + δ 3 B ) K 2 h b 2 + δ 3 ( δ 1 B ) K 4 h b 4
If we choose B = ( 1 / 3 + β 1 + β 2 ) , the dispersion expressions of Model 1 and Model 3 are equivalent. Similarly, the expression of Model 2 is equivalent to that of Model 4 if the parameters are chosen appropriately. Therefore, in the following section, we will focus on Model 1 and Model 2. As suggested by Madsen and Schäffer [45], these dispersion expressions can be closer to the Padé [2, 2] or Padé [4, 4] expansion of the Stokes analytic solutions without considering the porosity effect, if appropriate parameter values are chosen. The coefficients are summarized in Table 1, where two sets of parameter values are given for Model 1.
The following analytic dispersive Expression [43] is adopted to evaluate the accuracy of the models
ω 2 ( c r + α 1 i / ω ) = g K tanh ( K h b )
Phase celerity (c = ω/kr) of the different models is compared to the analytic solution for α1 = 0.2, 2, and 4 s−1 in Figure 2, Figure 3 and Figure 4 (where c is normalized by the analytic value), where Model 1 uses the two set group parameters in Table 1. It can be seen from the figures that Model 2 presents much more accurate results than Model 1, and the maximum error is only 2.5% when the applicable range of water depth is k r h = 6.0 The applicable range of water depth for Model 1 using the second set of parameter value is k r h 4.05 ,   4.57 ,   1.58 with a 5% tolerance error for three considered cases, which is higher than that using the first set of parameter value k r h 3.2 ,   3.75 ,   1.67 .
Defining c ( k i ) as the ratio of imaginary wave number between the model equations and analytic solutions, we could obtain the variation trend of c ( k i ) with water depth for Models 1 and 2, as shown in Figure 5, Figure 6 and Figure 7 for α 1 = 0.2, 2, and 4 s−1. Model 2 again presents better results than Model 1. The figures show that Model 2 is applicable for water depth k r h 4.19 , 4.41 , 4.98 for three cases within 2% error, whilst Model 1 with the first set of parameter values is only applicable for k r h 2.0 , 2.16 , 1.55 with 5% error, and with the second set of values only for k r h 2.62 , 2.68 , 1.52 . According to the above analysis, Model 2 is better adapted to deep water conditions, and its accuracy in deep water varies for different linear force coefficients.
The above comparisons are carried out without considering the effect of the inertial force coefficient c r , now let us discuss its effect. From the expression c r = 1 + ( 1 n ) c m , we can know that c r = 1 is the minimum value. For α 1 = 1, the phase celerity of Model 2 is compared against the analytical solutions for c r = 1, 1.1, 1.2, and 1.3, which is plotted in Figure 8. It shows that the applicable water depth decreases with the increase in c r , and the value of the inertial force coefficient c r greatly affects the results with the increase in the water depth.

4. Numerical Implementation and Validation

4.1. Numerical Implementation

Considering the overall accuracy of the derived model equations, the one-dimensional version of Model 2 is chosen for numerical implementation. The details of the numerical implementation are given in Kirby et al. [49], with a summary of the core procedures provided for clarity.
The one-dimensional equations of Model 2 are discretized on uniform grids. Higher-order finite difference formulations are used to approximate both the temporal and spatial derivatives. A composite fourth-order Adams–Bashforth–Moulton integration scheme is used for time marching. In the calculation of η and u values at time step n + 1, the third-order Adams–Bashforth time scheme is used for the prediction step, and the fourth-order Adams–Moulton scheme is used for the correction step, which is:
η n + 1 η n = Δ t 12 ( 23 f n 16 f n 1 + 5 f n 2 )
η n + 1 η n = Δ t 24 ( 9 f n + 1 + 19 f n 5 f n 1 + f n 2 )
where the superscript n denotes the time step corresponding to the parameter value at time n × Δ t and n + 1 corresponds to the parameter value at time n + 1 × Δ t , f represents the remaining terms in Equations (40) and (52) except the time derivative term.
The derivatives of the space in Equations (40) and (52) are discretized by the 4th order Taylor formats as Equations (66) and (67), which uses 5-point data on the spatial grid, so the velocity components are solved by a bandwidth of 5 broadband solution.
f x = f i + 2 + 8 f i + 1 8 f i 1 + f i + 2 / 12 Δ x + O ( Δ x ) 4
f x x = f i 2 + 16 f i + 1 30 f i + 16 f i 1 f i + 2 / 12 Δ x 2 + O ( Δ x ) 4
To enhance numerical stability, the model is solved using a predictor–corrector iterative approach. During the correction stage, the process is repeated until the difference in calculated variables between the two iterations reaches a pre-determined value, typically 0.0001. And the convergence of the iterations is accelerated by using an over-relaxation technique proposed by Kirby et al. [49].
Sponger layers are placed on either side of the numerical flume to absorb wave energy and wave generation occurs internally in the computational domain through the addition of a source function to the mass continuity equation. In this paper, we numerically simulate the propagation of a solitary wave in a porous media, due to the transient characteristics of the solitary wave, the surface elevation and velocities of the solitary wave are directly specified in the computational domain as the initial conditions to generate solitary waves, which is:
η t = H sec h 2 3 H 4 h 3 x c t
c = g H 1 + 1 2 H h g H + h
u = g h η h
where η is the water surface change; H is the solitary wave height; h is the water depth; t is the time; c is the wave speed; g is the acceleration of gravity. u is the solitary wave of any mass in the water body of the horizontal velocity.

4.2. Model Validation

The experimental data of the solitary wave interaction with the porous breakwater were obtained by Vidal et al. [50], which has been used to verify the numerical models of Lynett et al. [4] and Lin and Karunarathna [51]. The sketch of the solitary wave interaction with the porous breakwater is shown in Figure 9, where H and h are the wave height and still water depth respectively, and a, d50, and n are the porous structure width, mean diameter, and material porosity. Due to the transient nature of solitary waves, incident and reflected waves do not superimpose in the flume. Wave heights of reflected and transmitted waves on the left and right sides of the porous structure can be directly obtained. Therefore, the calculated reflection coefficient is defined as Kr = HL/H0, and the transmission coefficient is defined as Kt= HR/H0, where H0 is the input wave height of the solitary wave in the wave generation region, and HL is the wave height on the left side of the solitary wave after passing through the porous structure, HR is the wave height on the right side of the solitary wave. The numerical results were also compared with the experimental data.
In the experiment of Vidal et al. [50], the porous breakwater width is 20 cm or 40 cm, and the water depth ranges from 25 cm to 32 cm. The material porosity is 0.44 with two particle diameters d50 = 1.43 cm and 2.43 cm. In the numerical simulation, the space interval is 0.01 m, and the time interval is 0.005 s. By setting α 1 = 18 s−1 and α 2 = 234.45 m for the d50 = 1.43 cm case, and α 1 = 6.23 s−1 and α 2 = 137.97 m for the d50 = 2.43 cm case, the reflection coefficient and transmission coefficient are calculated numerically using Model 2, as shown in Figure 10. For the breakwater with d50 = 2.43 cm, both the calculated reflection and transmission coefficients are in good agreement with the experimental results. However, for the breakwater with d50 = 1.43 cm, the calculated reflection coefficients are slightly larger than the experimental data, and the calculated transmission coefficients are slightly smaller than the experimental data, but the trend with wave height (H/h) is consistent with the measured values. Considering the errors introduced by the modeling assumptions and parameters, the model developed in this paper can better simulate the wave interaction with porous structures which has surface-piercing boundaries.
The experiments of Lynett et al. [4] are nearly identical to those of Vidal et al. [50], where four types of porous breakwaters were tested. The gravels in two of the tests had mean diameters of 1.6 and 2.0 cm, and the breakwater length was 15 cm. The mean diameters of the others were 1.6 and 2.0 cm, with a breakwater length of 30 cm. The porosity of the porous media was n = 0.5. The water depth was 10 cm. The solitary waves had amplitudes between 1–4 cm. The reflection and transmission coefficient are computed using the free surface time series recorded at 1 m in front of and 1 m behind the breakwater. In the numerical simulation, the spatial interval is 0.01 m, and the time interval is 0.005 s. For d50 = 1.6 cm, α 1 = 9.5 s−1 and α 2 = 200 m are used, and for d50 = 2.0 cm, α 1 = 8 s−1 and α 2 = 160 m are used. Figure 11 shows the numerical reflection and transmission coefficients and the measurements from Lynett et al. [4]. For the 15 cm breakwater, both the calculated reflection and transmission coefficients are in good agreement with the experimental results. However, for the 30 cm breakwater, only the reflection coefficients are in good agreement with the experimental data. The numerical results for the thin breakwater are slightly better than those for the wide breakwater. In general, the agreements between the model results and the experiment data are acceptable after choosing the appropriate coefficients.

5. Conclusions

To accurately capture the wave motion in porous structures, four sets of Boussinesq-type equations with different linear and nonlinear performances are derived by including the drag resistance and inertial force induced by the media porosity. Two of them can be considered slightly different versions of Hsiao et al. [43], and the others are the new Boussinesq-type models, which have more accuracy of the second-order nonlinear terms compared to Fang et al. [44]. Two Boussinesq-type models (Mode 2 and Model 4) have Padé [4, 4] dispersion with accurate second or fourth-order nonlinearity, while the other Boussinesq-type models (Mode 1 and Model 3) have Padé [2, 2] dispersion with accurate second-order nonlinearity.
The phase celerity and damping rate of the four models are analyzed and compared with the analytical solutions. The results show that the fourth-order dispersive Boussinesq-type equations have accurate linear phase celerity and damping rate. In particular, the model works in the range of 0 < kh < 6.19 for damping rate and 0 < kh < 6.06 for phase celerity.
Some preliminary numerical results are presented for the solitary wave propagating through a vertical porous breakwater. The calculated reflection and transmission coefficients are compared with published experimental data, the variation trend with wave height (H/h) is consistent with the measured values, and the agreements between the model result and experiment data are acceptable after selecting the appropriate coefficients.
Porous structures have an effect on the nonlinearity of waves propagation. The analysis of the nonlinear performance of the models and the implementation of the two-dimensional numerical model are left for further study.

Author Contributions

Methodology, P.W., Z.L. and K.F.; Software, P.W.; Validation, P.W. and Z.L.; Formal Analysis, Z.L. and K.F.; Data Curation, P.W.; Writing—Original Draft, P.W., Z.L., K.F. and J.S.; Writing—Review and Editing, P.W., Z.L. and D.G.; Visualization, K.F. and J.S. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank the National Key Research and Development Program of China (2021YFB2600201 and 2022YFC3106101) and the National Natural Science Foundation of China (52171247 and 52071057) for their financial support.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Losada, I.J.; Patterson, M.D.; Losada, M.A. Harmonic generation past a submerged porous step. Coast. Eng. 1997, 31, 281–304. [Google Scholar] [CrossRef]
  2. Gu, Z.; Wang, H. Gravity waves over porous bottoms. Coast. Eng. 1991, 15, 497–524. [Google Scholar] [CrossRef]
  3. Lara, J.L.; Losada, I.J.; Liu, P.L.F. Breaking waves over a mild gravel slope: Experimental and numerical analysis. J. Geophys. Res. 2006, 111, C11019. [Google Scholar] [CrossRef]
  4. Lynett, P.J.; Liu, P.L.F.; Losada, I.J. Solitary wave interaction with porous breakwaters. J. Waterw. Port Coast. Ocean Eng. 2000, 126, 314–322. [Google Scholar] [CrossRef]
  5. Losada, I.J.; Lara, J.L.; Guanche, R.; Gonzalez-Ondina, J.M. Numerical analysis of wave overtopping of rubble mound breakwaters. Coast. Eng. 2008, 55, 47–62. [Google Scholar] [CrossRef]
  6. Cheng, Y.Z.; Jiang, C.B.; Wang, Y.Y. A coupled numerical model of wave interaction with porous medium. Ocean Eng. 2009, 36, 952–959. [Google Scholar] [CrossRef]
  7. Higuera, P.; Lara, J.L.; Losada, I.J. Three-dimensional interaction of waves and porous coastal structures using OpenFOAM®. Part I: Formulation and validation. Coast. Eng. 2014, 83, 243–258. [Google Scholar] [CrossRef]
  8. Sasikumar, A.; Kamath, A.; Bihs, H. Modeling porous coastal structures using a level set method based VRANS-solver on staggered grids. Coast. Eng. J. 2020, 62, 198–216. [Google Scholar] [CrossRef]
  9. Mohamed, K. A finite volume method for numerical simulation of shallow water models with porosity. Comput. Fluids 2014, 104, 9–19. [Google Scholar] [CrossRef]
  10. Ren, B.; Wen, H.; Dong, P.; Wang, Y. Improved SPH simulation of wave motions and turbulent flows through porous media. Coast. Eng. 2016, 107, 14–27. [Google Scholar] [CrossRef]
  11. Gao, J.L.; Zhou, X.J.; Zhou, L.; Zang, J.; Chen, H.Z. Numerical investigation on effects of fringing reefs on low-frequency oscillations within a harbor. Ocean Eng. 2019, 172, 86–95. [Google Scholar] [CrossRef]
  12. Gao, J.L.; Ma, X.Z.; Zang, J.; Dong, G.H.; Ma, X.J.; Zhu, Y.Z.; Zhou, L. Numerical investigation of harbor oscillations induced by focused transient wave groups. Coast. Eng. 2020, 158, 103670. [Google Scholar] [CrossRef]
  13. Gao, J.L.; Ma, X.Z.; Chen, H.Z.; Zang, J.; Dong, G.H. On hydrodynamic characteristics of transient harbor resonance excited by double solitary waves. Ocean Eng. 2021, 219, 108345. [Google Scholar] [CrossRef]
  14. Buccino, M.; Tuozzo, S.; Ciccaglione, M.C.; Calabrese, M. Predicting Crenulate Bay Profiles from Wave Fronts: Numerical Experiments and Empirical Formulae. Geosciences 2021, 11, 208. [Google Scholar] [CrossRef]
  15. Kirby, J.T. Boussinesq models and applications to nearshore wave propagation, surf zone processes and wave induced-current. In Advances in Coastal Modelling; Lakhan, V.C., Ed.; Elsevier Science: Amsterdam, The Netherlands, 2003; pp. 1–41. [Google Scholar]
  16. Madsen, P.A.; Fuhrman, D.R. High-order Boussinesq-type modeling of nonlinear wave phenomena in deep and shallow water. In Advances in Numerical Simulation of Nonlinear Water Waves; Ma, Q.W., Ed.; World Scientific Publishing Co. Pte. Ltd.: Singapore, 2010; pp. 245–285. [Google Scholar]
  17. Brocchini, M.A. reasoned overview on Boussinesq-type models: The interplay between physics, mathematics and numerics. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2013, 469, 20130496. [Google Scholar] [CrossRef]
  18. Kirby, J.T. Boussinesq models and their application to coastal processes across a wide range of scales. J. Waterw. Port Coast. Ocean Eng. 2016, 142, 03116005. [Google Scholar] [CrossRef]
  19. Sun, J.W.; Fang, K.Z.; Liu, Z.B.; Fan, H.X.; Sun, Z.C.; Wang, P. A review on the theory and application of Boussinesq-type equations for water waves. Haiyang Xuebao 2020, 42, 1–11. (In Chinese) [Google Scholar]
  20. Gao, J.L.; Shi, H.B.; Zang, J.; Liu, Y.Y. Mechanism analysis on the mitigation of harbor resonance by periodic undulating topography. Ocean Eng. 2023, 281, 114923. [Google Scholar] [CrossRef]
  21. Gao, J.L.; Ma, X.Z.; Dong, G.H.; Chen, H.Z.; Liu, Q.; Zang, J. Investigation on the effects of Bragg reflection on harbor oscillations. Coast. Eng. 2021, 170, 103977. [Google Scholar] [CrossRef]
  22. Madsen, P.A.; Murray, R.; Sørensen, O.R. A new form of the Boussinesq equations with improved linear dispersion characteristics. Coast. Eng. 1991, 15, 371–388. [Google Scholar] [CrossRef]
  23. Madsen, P.A.; Sørensen, O.R. A new form of the Boussinesq equations with improved linear dispersion characteristics. Part 2. A slowly-varying bathymetry. Coast. Eng. 1992, 18, 183–204. [Google Scholar] [CrossRef]
  24. Nwogu, O. Alternative form of Boussinesq equations for nearshore wave propagation. J. Waterw. Port Coast. Ocean Eng. 1993, 119, 618–638. [Google Scholar] [CrossRef]
  25. Schäffer, H.A.; Madsen, P.A. Further enhancements of Boussinesq-type equations. Coast. Eng. 1995, 26, 1–14. [Google Scholar] [CrossRef]
  26. Zou, Z.L. Higher-order Boussinesq equations for rapidly varying topography. Haiyang Xuebao 2001, 23, 109–119. (In Chinese) [Google Scholar]
  27. Liu, Z.B.; Sun, Z.C. Two sets of higher-order Boussinesq-type equations for water waves. Ocean Eng. 2005, 32, 1296–1310. [Google Scholar] [CrossRef]
  28. Wei, G.; Kirby, J.T.; Grilli, S.T.; Subramanya, R. A fully nonlinear Boussinesq model for surface waves. Part I. Highly nonlinear unsteady waves. J. Fluid Mech. 1995, 294, 71–92. [Google Scholar] [CrossRef]
  29. Gobbi, M.F.; Kirby, J.T.; Wei, G. A fully nonlinear Boussinesq model for surface waves. Part 2. Extension to O(kh)4. J. Fluid Mech. 2000, 405, 181–210. [Google Scholar] [CrossRef]
  30. Lynett, P.J.; Liu, P.L.F. Linear analysis of the multi-layer model. Coast. Eng. 2004, 51, 439–454. [Google Scholar] [CrossRef]
  31. Zou, Z.L.; Fang, K.Z. Alternative forms of the higher-order Boussinesq equations: Derivations and validations. Coast. Eng. 2008, 55, 506–521. [Google Scholar] [CrossRef]
  32. Liu, Z.B.; Fang, K.Z. Two-layer Boussinesq models for coastal water waves. Wave Motion 2015, 57, 88–111. [Google Scholar] [CrossRef]
  33. Madsen, P.A.; Bingham, H.B.; Liu, H. A new method for fully nonlinear waves from shallow water to deep water. J. Fluid Mech. 2002, 462, 1–30. [Google Scholar] [CrossRef]
  34. Chazel, F.; Benoit, M.; Ern, A.; Piperno, S. A double-layer Boussinesq-type model for highly nonlinear and dispersive waves. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2009, 465, 2319–2346. [Google Scholar] [CrossRef]
  35. Liu, Z.B.; Fang, K.Z. A new two-layer Boussinesq model for coastal waves from deep to shallow water: Derivation and analysis. Wave Motion 2016, 67, 1–14. [Google Scholar] [CrossRef]
  36. Liu, Z.B.; Fang, K.Z.; Cheng, Y.Z. A new multi-layer irrotational Boussinesq-type model for highly nonlinear and dispersive surface waves over a mildly sloping seabed. J. Fluid Mech. 2018, 842, 323–853. [Google Scholar] [CrossRef]
  37. Liu, Z.B.; Han, P.X.; Fang, K.Z.; Liu, Y. A high-order nonlinear Boussinesq-type model for internal waves over a mildly-sloping topography in a two-fluid system. Ocean Eng. 2023, 285, 115283. [Google Scholar] [CrossRef]
  38. Cruz, E.C.; Isobe, M.; Watanabe, A. Boussinesq equations for wave transformation on porous beds. Coast. Eng. 1997, 30, 125–156. [Google Scholar] [CrossRef]
  39. Hsiao, S.; Liu, P.L.F.; Chen, Y. Nonlinear water waves propagating over a permeable bed. Philos. Trans. R. Soc. Lond. Ser. A 2002, 458, 1291–1322. [Google Scholar] [CrossRef]
  40. Chen, Q. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. J. Eng. Mech. 2006, 132, 220–230. [Google Scholar] [CrossRef]
  41. Liu, Z.B.; Sun, Z.C. Wave propagating model over a porous seabed. China Sci. Pap. 2011, 6, 374–379. (In Chinese) [Google Scholar]
  42. Klonaris, G.T.; Memos, C.D. Compound Boussinesq-type modelling over porous beds. Appl. Ocean Res. 2020, 105, 102422. [Google Scholar] [CrossRef]
  43. Hsiao, S.; Hu, K.; Hwung, H. Extended Boussinesq Equations for Water-Wave Propagation in Porous Media. J. Eng. Mech. 2010, 136, 625–640. [Google Scholar] [CrossRef]
  44. Fang, K.Z.; Huang, M.H.; Chen, G.L.; Wu, J.; Wu, H.; Jiang, T. Boussinesq Simulation of Coastal Wave Interaction with Bottom-Mounted Porous Structures. J. Mar. Sci. Eng. 2022, 10, 1367. [Google Scholar] [CrossRef]
  45. Madsen, P.A.; Schäffer, H.A. Higher-order Boussinesq-type equations for surface gravity waves: Derivation and analysis. Philos. Trans. R. Soc. Lond. Ser. A 1998, 356, 3123–3184. [Google Scholar] [CrossRef]
  46. Zou, Z.L.; Liu, Z.B.; Fang, K.Z. Further improvements to the higher-order Boussinesq equations: Bragg reflection. Coast. Eng. 2009, 56, 672–687. [Google Scholar] [CrossRef]
  47. Peregrine, D.H. Long waves on a beach. J. Fluid Mech. 1967, 27, 815–827. [Google Scholar] [CrossRef]
  48. Liu, Z.B.; Fang, K.Z.; Zou, Z.L. Boussinesq wave equations with full nonlinear characteristics at order O(µ2). J. Harbin Eng. Univ. 2012, 33, 556–561. (In Chinese) [Google Scholar]
  49. Kirby, J.T.; Wei, G.; Chen, Q.; Kennedy, A.B.; Dalrymple, R.A. FUNWAVE 1.0 Fully Nonlinear BOUSSINESQ Wave Model Documentation and User’s Manual; Report, Center for Applied Coastal Research; University of Delawar: Newwark, NJ, USA, 1998. [Google Scholar]
  50. Vidal, C.; Losada, M.A.; Medina, R.; Rubio, J. Solitary wave transmission through porous breakwaters. In Proceedings of the 24th International Conference on Coastal Engineering, Costa del Sol-Malaga, Spain, 20–25 June 1988; pp. 1073–1083. [Google Scholar]
  51. Lin, P.; Karunarathna, S.A. Numerical study of Solitary wave interaction with porous breakwater. J. Waterw. Port Coast. Ocean Eng. 2007, 133, 352–363. [Google Scholar] [CrossRef]
Figure 1. The sketch of wave propagation in porous media.
Figure 1. The sketch of wave propagation in porous media.
Water 15 03900 g001
Figure 2. Non-dimensional wave celerity versus water depth ( α 1 = 0.2 , c r = 1 ).
Figure 2. Non-dimensional wave celerity versus water depth ( α 1 = 0.2 , c r = 1 ).
Water 15 03900 g002
Figure 3. Non-dimensional wave celerity versus water depth ( α 1 = 2 , c r = 1 ).
Figure 3. Non-dimensional wave celerity versus water depth ( α 1 = 2 , c r = 1 ).
Water 15 03900 g003
Figure 4. Non-dimensional wave celerity versus water depth ( α 1 = 4 , c r = 1 ).
Figure 4. Non-dimensional wave celerity versus water depth ( α 1 = 4 , c r = 1 ).
Water 15 03900 g004
Figure 5. Non-dimensional imaginary wave number versus water depth ( α 1 = 0.2 , c r = 1 ).
Figure 5. Non-dimensional imaginary wave number versus water depth ( α 1 = 0.2 , c r = 1 ).
Water 15 03900 g005
Figure 6. Non-dimensional imaginary wave number versus water depth ( α 1 = 2 , c r = 1 ).
Figure 6. Non-dimensional imaginary wave number versus water depth ( α 1 = 2 , c r = 1 ).
Water 15 03900 g006
Figure 7. Non-dimensional imaginary wave number versus water depth ( α 1 = 4 , c r = 1 ).
Figure 7. Non-dimensional imaginary wave number versus water depth ( α 1 = 4 , c r = 1 ).
Water 15 03900 g007
Figure 8. The effect of c r on wave phase celerity and imaginary wave number.
Figure 8. The effect of c r on wave phase celerity and imaginary wave number.
Water 15 03900 g008aWater 15 03900 g008b
Figure 9. The sketch of solitary wave interaction with a porous breakwater in the experiments by Vidal et al. [50] and Lynett et al. [4].
Figure 9. The sketch of solitary wave interaction with a porous breakwater in the experiments by Vidal et al. [50] and Lynett et al. [4].
Water 15 03900 g009
Figure 10. The comparison between computed (line) and measured (points) reflection coefficient (solid line and rhombus) and transmission coefficient (dashed line and circle) in Vidal et al.’s experiment [50].
Figure 10. The comparison between computed (line) and measured (points) reflection coefficient (solid line and rhombus) and transmission coefficient (dashed line and circle) in Vidal et al.’s experiment [50].
Water 15 03900 g010
Figure 11. The comparison between computed (line) and measured (points) reflection coefficient (solid line and rhombus) and transmission coefficient (dashed line and circle) in Lynett et al.’s experiment [4].
Figure 11. The comparison between computed (line) and measured (points) reflection coefficient (solid line and rhombus) and transmission coefficient (dashed line and circle) in Lynett et al.’s experiment [4].
Water 15 03900 g011
Table 1. The determined values for coefficient in different models.
Table 1. The determined values for coefficient in different models.
Model SetsParameter Values
1β1 = −0.0013, β2 = −0.0654 (I) or β1 = 0.0073, β2 = −0.064 (II)
2γ1 = 1/9, γ2 = 0.146488, γ3 = 1/945, γ4 = 0.00798359
3B = −0.4 or B = −0.395
4δ1= 0.101, δ3 = 0.039, B = −0.305, δ3 = 0.082 *, δ4 = 0.162 *
Note: * The two parameters are determined by reanalysis of the shoaling property, which differs from Madsen and Schäffer [45].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, P.; Liu, Z.; Fang, K.; Sun, J.; Gou, D. High-Order Boussinesq Equations for Water Wave Propagation in Porous Media. Water 2023, 15, 3900. https://doi.org/10.3390/w15223900

AMA Style

Wang P, Liu Z, Fang K, Sun J, Gou D. High-Order Boussinesq Equations for Water Wave Propagation in Porous Media. Water. 2023; 15(22):3900. https://doi.org/10.3390/w15223900

Chicago/Turabian Style

Wang, Ping, Zhongbo Liu, Kezhao Fang, Jiawen Sun, and Daxun Gou. 2023. "High-Order Boussinesq Equations for Water Wave Propagation in Porous Media" Water 15, no. 22: 3900. https://doi.org/10.3390/w15223900

APA Style

Wang, P., Liu, Z., Fang, K., Sun, J., & Gou, D. (2023). High-Order Boussinesq Equations for Water Wave Propagation in Porous Media. Water, 15(22), 3900. https://doi.org/10.3390/w15223900

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop