Next Article in Journal
A Simplified Mathematical Formulation for Water Quality Index (WQI): A Case Study in the Kelani River Basin, Sri Lanka
Next Article in Special Issue
A Numerical Study on Distant Tsunami Propagation Considering the Strong Nonlinearity and Strong Dispersion of Waves, or the Plate Elasticity and Mantle Fluidity of Earth
Previous Article in Journal
Modeling and Parameterization of the Evaporation and Thermal Decomposition of an Iron(III) Nitrate Nonahydrate/Ethanol Droplet for Flame Spray Pyrolysis
Previous Article in Special Issue
Impact of the Dissipation on the Nonlinear Interactions and Turbulence of Gravity-Capillary Waves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Waves Passing over Rectangular Obstacles: Multimodal Method and Experimental Validation

1
Laboratoire de Physique et Mécanique des Milieux Hétérogènes, UMR CNRS 7636, ESPCI-Paris, PSL Research University, Sorbonne Université, Paris Diderot University, CEDEX 5, 75231 Paris, France
2
Institut Langevin, UMR CNRS 7587, ESPCI-Paris, 75005 Paris, France
3
Laboratoire d’Acoustique de l’Université du Maine, UMR CNRS 6613, CEDEX 9, 72085 Le Mans, France
*
Author to whom correspondence should be addressed.
Fluids 2022, 7(5), 145; https://doi.org/10.3390/fluids7050145
Submission received: 14 March 2022 / Revised: 19 April 2022 / Accepted: 21 April 2022 / Published: 23 April 2022
(This article belongs to the Special Issue Nonlinear Wave Hydrodynamics, Volume II)

Abstract

:
We report a theoretical and experimental investigation of the propagation of nonlinear waves passing over a submerged rectangular step. A multimodal method allows calculating the first- and second-order reflected and transmitted waves. In particular, at the second order, the propagation of free and bound waves is theoretically presented. A detailed analysis of the convergence of the second-order problem shows that a finite truncation of the series of evanescent bound waves is necessary to obtain a smooth and convergent solution. The computed coefficients of the first and second harmonics are experimentally validated via a complete space-time-resolved measurements of the wave propagation, which permits us to verify the relative amplitude, phase and spatial interference (beating) of the free and bound waves at the second order. This result can be useful in future multimodal models since it not only keeps the accuracy of the model with the inclusion of the first part of the evanescent bound terms (being also the dominants) but also ensures the convergence of the multimodal computation with an error that decreases as a function of the number of modes.

1. Introduction

Depth transitions are ubiquitous in coastal environments, being a crucial zone for understanding the hydrodynamics and morphology. For instance, recent investigations showed that natural hazards such as extreme waves have a place preferentially in large-depth discontinuities (Li et al. [1]). Due to widespread interest, a large variety of methods can be found in the mathematical modeling of surface water waves over variable bathymetry. In particular for the step problem, which tries to solve the wave propagation over an abrupt depth discontinuity, the scientific community’s interest has remained high due to its canonical character. Among the earliest models that solve the linear problem, we find Mei and Black [2], who developed a variational formulation, as well as Miles [3], who calculated the scattering matrix with the transmission and reflection coefficients in both directions. Moreover, the model proposed by Newman [4] considered an infinite depth before the step, which can be a reasonable simplification under certain conditions. In the same way, Boussinesq-type equations have been used in the step problem, as developed by Grue [5]. In the case of multimodal methods, the model developed by Massel [6] solves the problem of propagation of nonlinear waves over a submerged vertical step. The appeal of this model is the possibility of obtaining the reflection and transmission coefficients at the first and second order. Similarly, in the weakly nonlinear regime, more general bathymetries have been studied in Belibassakis and Athanassoulis [7] and Belibassakis and Athanassoulis [8]. Furthermore, in the linear regime, supplementary modes were proposed by Athanassoulis and Belibassakis [9], obtaining faster convergence. Moreover, oblique incident waves were considered by Rhee [10],who obtained interesting results for the phase shift of the transmitted and reflected waves. More recently, in the linear regime, Porter and Porter [11] used a conformal mapping of the fluid domain to transfer the steep deformations of the bottom topography into smooth functions applied to the modified free surface boundary conditions. More recently, some models have studied the propagation of surface wave packets (Li et al. [12]) and solitons (Ducrozet et al. [13]) over submerged steps.
In the propagation of water waves, depth transitions usually generate nonlinearities at different levels, going from weakly nonlinear waves up to breaking waves. The behavior of these nonlinear waves over submerged obstacles can be found in the theory of propagation of waves in shallow water by Massel [6] or Bryant [14]. These theories predict the evolution of waves in constant depth for different harmonics, where the second harmonic is of particular interest because of the presence of waves at different velocities. These waves are separated into free waves, which follow the dispersion relation in all the harmonics, and bound waves which are directly created by nonlinearities and are not dispersive. Thus, for a given forcing frequency ω , at the n-th harmonic, the free wave has a wavenumber given by k n = D ( n ω ) , where D represents the linear dispersion relation of water waves. On the other side, the bound waves at the n-th harmonic has a wavenumber n k 1 = n D ( ω ) , being merely an integer multiple of the fundamental mode wavenumber. The mechanism of decomposition of monochromatic incident waves into free and bound waves can be found in the characterization by Huang and Dong [15]. They described the mechanism involved when the primary wave passes over the dike and generates a secondary small wave that propagates with a lower velocity than the main crest. This difference in celerity corresponds to the deviation of the dispersion relation curve with respect to the non-dispersive harmonics (bound waves). Although the presence of free and bound waves is directly observable in the ( ω , k ) plane, which can be obtained from a two-dimensional Fourier transform of a space-time plane, most of the studies in nonlinear waves so far have considered punctual measurements, as carried out by Chapalain et al. [16], Beji and Battjes [17], Li and Ting [18], Benoit et al. [19] and Ohyama and Nadaoka [20]. These experimental works were mostly designed to validate numerical models by comparing temporal signals. Space-time measurements have been carried out by Brossard and Chagdali [21], Brossard et al. [22] with moving probes, where the Doppler effect permits the identification of free and bound waves in the k-spectrum. More recently, experiments performed by Li and Ting [18] and Ting et al. [23] with a vertical laser sheet measured space-time data along a channel, achieving an interesting separation of free and bound waves. However, this measurement is still in 2D, and it does not take into account the transverse effects that are always present in wave flume experiments.
In this study, we revisit a multimodal model developed by Massel [6] that gives accurate results for weakly nonlinear conditions, as verified experimentally by Ohyama et al. [24]. The numerical calculations are analyzed in terms of the convergence. We give special attention to the second order, whose source term has difficult convergence due to the abrupt depth discontinuity. Regarding the experimental validation, we report on a complete space-time-resolved measurement of nonlinear waves using Fourier Transform Profilometry (FTP) (Takeda et al. [25], Takeda and Mutoh [26]), which was adapted to surface water waves by Cobelli et al. [27]. A 3D field permits the verification of important experimental aspects such as wavefront flatness and transverse modes, making the measurements more reliable.

2. Multimodal Model

2.1. Statement of the Problem

Let us consider an infinite domain in the plane ( x , z ) where we place a depth discontinuity at x = 0 . Thus, as is shown in Figure 1, the depth h ( x ) is the piecewise function
h ( x ) = h , x < 0 h + < h , x > 0 ,
which describes a bathymetry composed of a deep water region in x < 0 and a shallow water region in x > 0 .
For simplicity, all the development is presented in the configuration called deep-to-shallow water, where an incident right-going wave comes from x = . Nevertheless, the model is perfectly applicable to the shallow-to-deep water case.
In both sides of the step, the velocity potential function ϕ satisfies the Laplace equation
2 x 2 + 2 z 2 ϕ = 0
and the potential ϕ and surface displacement η satisfy the free surface kinematic boundary conditions
ϕ z = η t + ϕ x η x , z = η
and the free surface dynamic boundary condition
ϕ t + 1 2 ϕ x 2 + ϕ z 2 + g η = 0 , z = η .
Eventually, on the bottom and on the step wall, we impose impervious boundary conditions
ϕ z = 0 , z = h ( x )
ϕ x = 0 , h < z < h + , x = 0
We start the solution of the problem via the perturbation method (Hsu et al. [28], Ohyama et al. [24] and Mei et al. [29]). To do so, let us consider a small parameter ϵ , which, in this dimensional form, corresponds to the wave amplitude a. This parameter is used in the formal expansion of the wave potential ϕ and the surface displacement η
ϕ ( x , z , t ) = ϵ ϕ 1 ( x , z , t ) + ϵ 2 ϕ 2 ( x , z , t ) +
η ( x , t ) = ϵ η 1 ( x , t ) + ϵ 2 η 2 ( x , t ) +
Next, a Taylor series approximated the velocity potential around z = 0 in the form
ϕ ( x , η , t ) = ϕ z = 0 + η ϕ z z = 0 + 1 2 η 2 2 ϕ z 2 z = 0 +
The power expansion of Equations (7) and (8) and the Taylor series of Equation (9) are replaced in the Laplace Equation (2) and in the free surface and bottom boundary conditions (Equations (3) to (6)). Then, we truncate the expressions at the order ϵ to obtain the first-order system with a Laplace Equation that reads
2 x 2 + 2 z 2 ϕ 1 = 0
with the boundary conditions at the free surface
2 ϕ 1 t 2 + g ϕ 1 z = 0 , z = 0 ,
and
g η 1 + ϕ 1 t = 0 , z = 0 ,
and the impervious bottom boundary conditions
ϕ 1 z = 0 , z = h ,
and
ϕ 1 x = 0 , h < z < h + , x = 0 .
We state the diffraction problem of water waves by setting at infinity, in x < 0 , that there are one incident wave and one reflected wave from the discontinuity. On the other hand, at the positive infinity in x > 0 , there is only one transmitted wave (outgoing wave). Eventually, at the depth transition, the potential ϕ and its spatial derivative are continuous in the whole domain, that is, [ ϕ ] = 0 and [ x ϕ ] = 0 , at x = 0 and at h + < z < 0 .

Solution of First Order

We define the space-time-dependent functions corresponding to the incident and reflected waves as
I 0 ( 1 ) ( x , t ) = g a ω e i ( k 0 x ω t ) , x < 0 ,
A n ( 1 ) ( x , t ) = g a ω e i ( k n x + ω t ) , x < 0 ,
B m ( 1 ) ( x , t ) = g a ω e i ( k m + x ω t ) , x > 0 ,
where k 0 ± corresponds to the unique real propagating solution and k n ± = i k 1 ± , i k 2 ± , i k 3 ± , to the infinite imaginary evanescent solutions of the dispersion relation of surface gravity waves
ω 2 = g k n ± tanh ( k n ± h ± ) .
Here, in Equations (16) and (17) and in what follows, for the sake of clarity, we shall use the variable n = 0 , 1 , 2 , 3 , as the index of the functions in the deep water region ( x < 0 , h ) and the variable m = 0 , 1 , 2 , 3 , as the index of the functions in the shallow water region ( x > 0 , h + ). In addition, the orthogonal basis function that satisfies the boundary conditions (11) to (13) is a z-dependent function in the form
F n ( 1 ) ( z ) = cosh k n ( z + h ) cosh k n h , x < 0 ,
G m ( 1 ) ( z ) = cosh k m + ( z + h + ) cosh k m + h + , x > 0 .
Hence, the linear potential for both deep and shallow water regions can be written as
ϕ 1 = I 0 ( 1 ) F 0 ( 1 ) + n R n ( 1 ) A n ( 1 ) F n ( 1 ) , x < 0 ,
ϕ 1 + = m T m ( 1 ) B m ( 1 ) G m ( 1 ) , x > 0 ,
where R n ( 1 ) and T m ( 1 ) are the unknown reflection and transmission coefficients of the first-order problem. From the continuity conditions at x = 0 for the potential ϕ and its spatial derivative, we obtain the following equations
ϕ 1 x = 0 = ϕ 1 + x = 0 +
ϕ 1 x x = 0 = ϕ 1 + x x = 0 + .
For faster convergence, we consider a number of modes N for the deep water region and N + h + / h N for the shallow water region. By replacing Equations (21) and (22) in the conditions (23) and (24), we obtain
I 0 ( 1 ) F 0 ( 1 ) + n R n ( 1 ) A n ( 1 ) F n ( 1 ) = m T m ( 1 ) B m ( 1 ) G m ( 1 ) ,
x I 0 ( 1 ) F 0 ( 1 ) + n R n ( 1 ) x A n ( 1 ) F n ( 1 ) = m T m ( 1 ) x B m ( 1 ) G m ( 1 ) .
Next, we project Equation (25) onto the orthogonal basis G m ( 1 )
I 0 ( 1 ) h + 0 F 0 ( 1 ) G m ( 1 ) + n R n ( 1 ) A n ( 1 ) h + 0 F n ( 1 ) G m ( 1 ) = T m ( 1 ) B m ( 1 ) h + 0 ( G m ( 1 ) ) 2 ,
where we have used the orthogonality of the basis G m ( 1 ) to eliminate the sum in m on the right-hand side, i.e., Equation (27) is written for each m with 0 m ( N + 1 ) . Before projecting Equation (26), we use the boundary condition from Equation (14) to impose the following equivalence of the integrals of the derivative
h 0 ϕ 1 x F n ( 1 ) d z = h + 0 ϕ 1 x F n ( 1 ) d z , x = 0 ,
which we use to project Equation (26) onto the orthogonal basis F n ( 1 ) , and obtain
x I 0 ( 1 ) h 0 F 0 ( 1 ) F n ( 1 ) + R n ( 1 ) x A n ( 1 ) h 0 ( F n ( 1 ) ) 2 = m T m ( 1 ) x B m ( 1 ) h + 0 F n ( 1 ) G m ( 1 ) ,
where we have used the orthogonality of the basis F n ( 1 ) to eliminate the sum in n on the left-hand side, i.e., Equation (29) is written for each n with 0 n ( N 1 ) . Therefore, we obtain N + N + projected equations.
Eventually, the system of Equations (27) and (29) can be solved for R n ( 1 ) , as is detailed in Appendix A.1, and then the transmission coefficients T m ( 1 ) can be easily obtained from Equation (27).

2.2. Second Order Problem

2.2.1. Statement of the Problem

Now after replacing Equations (7) to (9) in the system of Equation (5) to (2) to (4), we keep the terms at the order ϵ 2 . This leads to the nonlinear wave problem at second order where the solution ϕ 2 should satisfy the Laplace equation
2 x 2 + 2 z 2 ϕ 2 = 0
with the following boundary conditions at the free surface
2 ϕ 2 t 2 + g ϕ 2 z = t ϕ 1 x 2 + ϕ 1 z 2 η 1 z 2 ϕ 1 t 2 + g ϕ 1 z , z = 0 ,
g η 2 + ϕ 2 t = 1 2 ϕ 1 x 2 + ϕ 1 z 2 η 1 2 ϕ 1 z t , z = 0 ,
and impervious boundary conditions at the bottom
ϕ 2 z = 0 , z = h ( x ) ,
ϕ 2 x = 0 , h < z < h + , x = 0 .
This represents a non-homogeneous surface wave equation. Considering that the non-homogeneous wave problem has one particular and one homogeneous solution ( ϕ 2 = ϕ 2 b + ϕ 2 f ), we can state the boundary conditions separately at x = ± . First, the right side of Equation (31) depends on ϕ 1 , which implies that the particular solution ( ϕ 2 b ) has the same boundary conditions as for ϕ 1 . Thus, we have at the negative infinity, in x < 0 , there are one incident wave and one reflected wave, and at the positive infinity, in x > 0 , there is only one transmitted outgoing wave. Moreover, because the homogeneous second-order problem does not consider a source term with energy coming from the infinity, we have the Sommerfeld radiation condition (Sommerfeld [30]), which only considers outgoing waves in both sides of the domain; it is x = ± . Similarly to the first-order problem, the second-order problem should satisfy the continuity of the potential and its spatial derivative. Thus, we have jump conditions as [ ϕ 2 ] = 0 and [ x ϕ 2 ] = 0 at x = 0 and at h + < z < 0 .

2.2.2. Second-Order Particular Solution

The kinematic and dynamic boundary conditions expressed in Equations (31) and (32) have a source term depending on the solution of the linear problem. The solution of this equation is composed by one particular solution ϕ 2 b (bound wave) and one homogeneous solution ϕ 2 f (free wave). Thus, the complete solution is
ϕ 2 = ϕ 2 b + ϕ 2 f .
We start by finding the bound wave terms, which depend on ϕ 1 . We can find and write the incident second-order potential in terms of the following ( x , t ) and z-dependent functions
I 0 , 0 ( 1 , 1 ) ( x , t ) = C 0 , 0 ( 1 , 1 ) e 2 i k 0 x 2 i ω t + c . c .
F 0 , 0 ( 1 , 1 ) ( z ) = cosh 2 k 0 ( z + h ) cosh 2 k 0 h
where C 0 , 0 ( 1 , 1 ) is the constant of the second-order propagating bound wave (as in Equation (36), we shall use in what follows the super-index ( ) ( 1 , 1 ) to indicate the constants and functions resulting from replacing the solution of the first-order ϕ 1 and η 1 in the right side of Equations (31) and (32)). This second-order wave, usually called Stokes harmonic, is the solution of Equation (31) when there is only one incident wave (Equation (15)). Thus, it always exists in the propagation of a single wave in constant depth. The simplified expression of the Stokes harmonic is
I 0 , 0 ( 1 , 1 ) ( x , t ) F 0 , 0 ( 1 , 1 ) ( z ) = 3 8 ω a 2 cosh 2 k 0 ( z + h ) sinh 4 ( k 0 h ) e 2 i k 0 x 2 i ω t + c . c .
Likewise, the bound reflected waves can be written in terms of the ( x , t ) -dependent function
A n , n ( 1 , 1 ) ( x , t ) = C n , n ( 1 , 1 ) e i ( k n + k n ) x 2 i ω t + c . c .
where the constant C n , n ( 1 , 1 ) is the bound waves’ constant depending on the linear reflection coefficients and is detailed in the Appendix B.2.
For the shallow water region, the bound waves depend on the transmitted linear wave and are obtained by replacing Equation (22) in the right-hand side of Equation (31). Thus, the ( x , t ) -dependent function, in x > 0 , can be expressed as
B m , m ( 1 , 1 ) ( x , t ) = C m , m ( 1 , 1 ) e i ( k m + + k m + ) x 2 i ω t + c . c .
where the coefficient C m , m ( 1 , 1 ) has the same form as Equation (A5) but uses wavenumbers and depth from the shallow water region ( k m + , h + ). In order to satisfy the source term of Equations (31) and (32), the z dependency in the bound waves is expressed in the following form
F n , n ( 1 , 1 ) ( z ) = cosh ( k n + k n ) ( z + h ) cosh ( k n + k n ) h , x < 0 ,
G m , m ( 1 , 1 ) ( z ) = cosh ( k m + + k m + ) ( z + h + ) cosh ( k m + + k m + ) h + , x > 0 .
In addition, the second-order evanescent terms resulting from the multiplication of the Stokes incident harmonic with the evanescent reflected bound waves can be expressed as
D 0 , n ( 1 , 1 ) ( x ) = C 0 , n ( 1 , 1 ) e ( i k 0 k n ) x 2 i ω t + c . c .
where the constant coefficient C 0 , n ( 1 , 1 ) is detailed in Appendix B.2. The z-dependent function corresponds to
H 0 , n ( 1 , 1 ) ( z ) = cosh ( k 0 k n ) ( z + h ) cosh ( k 0 k n ) h , x < 0 .
Eventually, the second-order bound wave potential is
ϕ 2 b = I 0 , 0 ( 1 , 1 ) F 0 , 0 ( 1 , 1 ) + n n R n ( 1 ) R n ( 1 ) A n , n ( 1 , 1 ) F n , n ( 1 , 1 ) + n R n ( 1 ) D 0 , n ( 1 , 1 ) H 0 , n ( 1 , 1 ) + c . c . , x < 0 ,
ϕ 2 b + = m m T m ( 1 ) T m ( 1 ) B m , m ( 1 , 1 ) G m , m ( 1 , 1 ) + c . c . , x > 0 .

2.2.3. Second-Order Homogeneous Solution

In order to satisfy the homogeneous solution (without source terms from ϕ 1 ), the free waves should follow the dispersion relation at the frequency 2 ω . Thus, second-order dispersive terms are generated only at depth inhomogeneity due to the absence of dispersive terms in the incident wave. As a consequence, the Sommerfeld radiation condition is imposed to the free waves, which implies that the ( x , t ) -dependent functions have the following form
A p ( 2 ) ( x , t ) = g a 2 ω e i ( κ p x + 2 ω t ) + c . c . , x < 0
B q ( 2 ) ( x , t ) = g a 2 ω e i ( κ q + x 2 ω t ) + c . c . , x > 0
where κ p ± = κ 0 ± , i κ 1 ± , i κ 2 ± , i κ 3 ± , correspond to the unique real solution ( p = 0 ) and to the infinite imaginary solutions ( p = 1 , 2 , , N ) of the dispersion relation
( 2 ω ) 2 = g κ p ± tanh ( κ p ± h ± )
For the sake of clarity, we shall use, as in Equations (47) and (48), the sub-index p to identify the constants and functions of the second-order homogeneous solution in the deep water region ( x < 0 , h ) and the sub-index q to identify the constants and functions of the second-order homogeneous solution in the shallow water region ( x > 0 , h + ).
The orthogonal basis of z-dependent functions that satisfy the boundary conditions (31) to (34) is expressed in the same way as for the linear problem
F p ( 2 ) ( z ) = cosh κ p ( z + h ) cosh κ p h , x < 0 ,
G q ( 2 ) ( z ) = cosh κ q + ( z + h + ) cosh κ q + h + , x > 0 .
Next, the free wave potential at the second order corresponds to the following sum
ϕ 2 f = p R p ( 2 ) A p ( 2 ) F p ( 2 ) + c . c . , x < 0
ϕ 2 f + = q T q ( 2 ) B q ( 2 ) G q ( 2 ) + c . c . , x > 0
where the unknowns R p ( 2 ) and T q ( 2 ) are the transmission and reflection coefficients that shall be obtained in the next calculation.

2.2.4. Construction of the Complete Second-Order Solution

Similarly to the linear problem, we impose the continuity of the second-order potential and its derivative at the step position. Thus, we have
ϕ 2 = 0 , x = 0 ,
x ϕ 2 = 0 , x = 0 ,
which means that the free and bound wave potential should satisfy
ϕ 2 b + ϕ 2 f = ϕ 2 b + + ϕ 2 f + ,
x ϕ 2 b + x ϕ 2 f = x ϕ 2 b + + x ϕ 2 f + ,
which is the set of equations of the unknowns R p ( 2 ) and T q ( 2 ) .
Next, we project Equation (56) onto the orthogonal basis G q ( 2 ) , which yields
I 0 , 0 ( 1 , 1 ) h + 0 F 0 , 0 ( 1 , 1 ) G q ( 2 ) + n n R n ( 1 ) R n ( 1 ) A n , n ( 1 , 1 ) h + 0 F n , n ( 1 , 1 ) G q ( 2 ) + n R n ( 1 ) D 0 , n ( 1 , 1 ) h + 0 H 0 , n ( 1 , 1 ) G q ( 2 ) + p R p ( 2 ) A p ( 2 ) h + 0 F p ( 2 ) G q ( 2 ) + c . c . = m m T m ( 1 ) T m ( 1 ) B m , m ( 1 , 1 ) h + 0 G m , m ( 1 , 1 ) G q ( 2 ) + T q ( 2 ) B q ( 2 ) h + 0 ( G q ( 2 ) ) 2 + c . c .
where we have used the orthogonality of the basis G q ( 2 ) to eliminate the sum in q in the right-hand side of the equation, i.e., the projected Equation (58) is written for each q-th term.
Before projecting Equation (57), we use the boundary condition from Equation (34) to impose the following equivalence of the integrals of the derivative
h 0 ϕ 2 x d z = h + 0 ϕ 2 x d z , x = 0 ,
which is used to project Equation (57) onto the orthogonal basis F p ( 2 ) and to obtain
x I 0 , 0 ( 1 , 1 ) h 0 F 0 , 0 ( 1 , 1 ) F p ( 2 ) + n n R n ( 1 ) R n ( 1 ) x A n , n ( 1 , 1 ) h 0 F n , n ( 1 , 1 ) F p ( 2 ) + n R n ( 1 ) x D 0 , n ( 1 , 1 ) h 0 H 0 , n ( 1 , 1 ) F p ( 2 ) + R p ( 2 ) x A p ( 2 ) h 0 ( F p ( 2 ) ) 2 + c . c . = m m T m ( 1 ) T m ( 1 ) x B m , m ( 1 , 1 ) h + 0 G m , m ( 1 , 1 ) F p ( 2 ) + q T q ( 2 ) x B q ( 2 ) h + 0 G q ( 2 ) F p ( 2 ) + c . c . ,
Here, we have used the orthogonality of the basis F p ( 2 ) to eliminate the sum in p in the left-hand side of the equation, i.e., the projected Equation (60) is written for each p-th term.
Finally, the system of Equations (58) and (60) can be solved for R p ( 2 ) , as is detailed in Appendix B.2, and then the transmission coefficients T q ( 2 ) can be easily obtained from Equation (58).

3. Numerical Calculation

3.1. Solution of the First-Order Problem

We show in this section a numerical implementation of the first-order system of Equation (A1). This numerical example has been built with the parameters a = 0.023 m, h = 0.3 m, h + = 0.15 m and ω = 2 π s 1 . In Figure 2, we present the near field of the velocity potential ϕ 1 as well as the horizontal velocity x ϕ 1 calculated for N = 64 modes. As we observe, the first-order problem, considering it is homogeneous and linear, has a robust convergence with a smooth solution that is independent of the wave amplitude and frequency.
In Figure 2b, in the neighborhood of the corner, at ( x = 0 , z = 0.15 ) m, we observe a local divergence of the horizontal velocity, which is the classical behavior of a potential flow around a tip. In more detail, the profile of the matching of the potential and horizontal velocity is presented in Figure 3. The real part of the functions ϕ 1 ( x = 0 , z ) and ϕ 1 + ( x = 0 , z ) are plotted in Figure 3a, where a very good agreement is observed in the vertical range 0.15 < z < 0 m. On the other hand, in Figure 3b, we present the vertical profile of the match of the horizontal velocity, that is, x ϕ 1 ( x = 0 , z ) and x ϕ 1 + ( x = 0 , z ) . As expected, the divergence of the velocity around the corner is difficult to match in this part, showing nevertheless a good match in the rest of the domain. The boundary condition of the impervious wall at the step in Equation (14) is also well satisfied by the function x ϕ 1 ( x = 0 , z ) in the range 0.3 < z < 0.15 m, especially far from the rectangular corner. This match, as shall be shown further on, is improved as a function of the number of modes ( N , N + ) .

3.2. Solution of the Second-Order Problem

In the second-order problem, the presence of a source term in the right-hand side of Equations (31) and (32) can lead to convergence problems. The first natural trial is the solution of the second order considering all the terms obtained from the first-order solution, i.e., making the limits of the double and single sums in Equations (45) and (46) up to n = N and m = N + . The first term of these sums correspond to a propagating wave, whereas the following terms correspond to evanescent bound waves. Therefore, by injecting the first-order solution shown in Figure 2 in the source term of Equations (31) and (32) and solving the system of Equation (A7), we obtain a solution ϕ 2 and x ϕ 2 , as shown in Figure 4a,b, that is not smooth and has a clear divergence in the matching zone. This divergence, observed when the evanescent bound waves are included in the source term, has been recently reported by Li et al. [12], where, due to convergence problems, the contribution of bound evanescent waves was not considered.
Although the first-order problem has a good convergence, the second order shows problems when we increase the number of evanescent modes in the total bound wave, which is the result of the sum of all the evanescent modes obtained from Equations (A5) and (A6). The second order solves the problem of finding the proper free wave that compensates for the bound wave field in order to obtain a smooth solution (satisfying the matching conditions). To do so, the bound wave field should have a limited amplitude, i.e., it should not grow as a function of the number of modes.
In this analysis, we focus on n = n , which is the fastest growth in the double sum of Equation (45). Thus, we consider the term ( R n ( 1 ) ) 2 C n , n ( 1 , 1 ) . In Figure 5a, we show the decay of the first-order reflection coefficient and the slope n 2 for comparison. We observe that for the case of N = 100 , the decay is slower than n 2 for n > 50 approximately. In contrast, as we observe in Figure 5b, the bound wave coefficient | C n , n ( 1 , 1 ) | has growth at the power n 4 . Therefore, the decay of the first-order coefficients ( R n ( 1 ) , T m ( 1 ) ) does not compensate for the growth of the bound evanescent wave constants, and some modification should be carried out to ensure the convergence of the second-order problem.
From the observation of the decay of Figure 5a, we have that the exponent of the decay decreases close to the end of the series. This suggests that a compensation of the growth of the coefficients C n , n ( 1 , 1 ) can be possible if we consider only part of the series R n ( 1 ) , especially the first part that has faster decay. Therefore, we propose a truncation of the evanescent waves terms into the bound wave solution. To do so, we rewrite Equations (45) and (46) as
ϕ 2 b = I 0 , 0 ( 1 , 1 ) F 0 , 0 ( 1 , 1 ) + n = 0 α N n = 0 α N R n ( 1 ) R n ( 1 ) A n , n ( 1 , 1 ) F n , n ( 1 , 1 ) + n = 0 α N R n ( 1 ) D 0 , n ( 1 , 1 ) H 0 , n ( 1 , 1 ) + c . c . , x < 0
ϕ 2 b + = m = 0 α N + m = 0 α N + T m ( 1 ) T m ( 1 ) B m , m ( 1 , 1 ) G m , m ( 1 , 1 ) + c . c . , x > 0
where α is an arbitrary truncation coefficient in the range 0 < α 1 .
A systematic exploration of the convergence of the second-order solution as a function of α is shown in Figure 6a,b. The numerical error of the multimodal matching is calculated for the potential ϕ 2 as
E ( ϕ 2 ) = h + 0 | ϕ 2 + ( x = 0 , z ) ϕ 2 ( x = 0 , z ) | d z ,
and for the horizontal velocity x ϕ 2 as
E ( x ϕ 2 ) = h + 0 | x ϕ 2 + ( x = 0 , z ) x ϕ 2 ( x = 0 , z ) | d z ,
and is plotted in Figure 6a as a function of α , where we observe an abrupt increase in E ( ϕ 2 ) and E ( x ϕ 2 ) at α > 0.5 . For this calculation, we have considered N = 100 and N + = 50 modes and increased the truncation coefficient α progressively. Likewise, we have verified the evolution of the propagating coefficients at the second order. As we observe in Figure 6b, the coefficients R 0 ( 2 ) and T 0 ( 2 ) show a stable behavior up to the limit α = 0.5 , beyond which the coefficients diverge and change rapidly as a function of α . This threshold, after several numerical calculations, seems to be independent of input parameters, such as geometry, wave amplitude and frequency. Therefore, we consider in what follows that the truncation coefficient α = 0.5 ensures the convergence of the second-order solution.
A numerical example calculated with N = 64 , N + = 32 modes and α = 0.5 is presented in Figure 7a,b, for the near field of ϕ 2 and x ϕ 2 , respectively. Both fields, at the depth discontinuity, show a smooth transition between the deep and shallow water regions, confirming the convergence of the solution. The agreement of the multimodal matching can be verified more in detail in the vertical profile shown in Figure 8a,b, where an excellent match between the potential ϕ 2 and the velocity x ϕ 2 at both sides of the step confirms the accuracy of the multimodal matching after applying the truncation method.
Eventually, we show in Figure 9 the numerical error calculated in Equations (63) and (64) for the case α = 0.5 . In Figure 9a, we show the numerical error of the velocity potential for the first- and second-order problems, which are in both cases decaying. Similarly, in Figure 9b, the numerical error of the horizontal velocities at the first and second order also decay monotonically as a function of the number of modes.

4. Experimental Measurements

4.1. Experimental Set-Up

In this section, we present experimental space-time-resolved measurements of the propagation of nonlinear waves passing over a submerged step. The main objective of this section is provide quantitative measurements in such a way that the complex coefficients obtained in the multimodal model, especially for the second order, could be compared and validated with experiments. To do so, we take advantage of the Fourier transfer profilometry technique (FTP) to measure in space and time the surface deformation. In Figure 10a, we see a top view and a side view of the experimental setup. Experiments were carried out in a rectangular tank of 180 × 60 × 15 cm. A flap-type wavemaker is located 40 cm before the step from which waves propagate through a deep water region of 65 mm deep before passing over a rectangular step that separate the deep region from the shallow water region of 20 mm deep. The wave propagates 80 cm in the shallow region before attaining a 8% slope that avoid spurious reflections. As we see in Figure 10b, a sinosuidal pattern is projected onto a diffusive surface of water colored with titanuim dioxide (TiO2), which does not change the water’s physical properties (Przadka et al. [31]). The post-treatment of the patter deformed by the wave motion allows measuring the surface height (see Figure 10c) in any point of the image. The measured field has a spatial resolution of 1700 × 400 pixels covering a surface of 1 m long and 25 cm large. The amplitude of the wavemaker was calibrated in such a way that the wave amplitude measured at the step was constant at a 3 mm. The wave frequency with a monochromatic forcing varied in the range ω [ 10 , 20 ] s 1 . Additionally, spurious capillary waves coming from the lateral boundaries and affecting the wave front were canceled by the use of a wire-mesh in the lateral walls, as conducted by Monsalve et al. [32,33]. As we shall observe in the following experimental results and considering the small scale of the experiments, the influence of the surface tension is non-negligible at the second order, where the wave number k > 100 m 1 . Despite this, the main results show that the multimodal model of gravity waves, presented in the theoretical section, fitted the experimental results well in terms of transmitted and reflected waves at the first and second order. Afterwards, along the spatial propagation, the second-order free and bound waves enhance the difference due to the surface tension as an extra restoring force in the dispersion relation.

4.2. Space-Time Spectra

The measured field η ˜ ( x , y , t ) , considering that is invariant in the transverse direction y, is averaged in order to obtain a one-dimension measurement η ( x , t ) , which is presented as a space-time graph in Figure 11. Here, the origin of the abscissa is located at the step position. In x < 0 , the deep water region shows high modulation of the amplitude, indicating a high reflection from the step. On the other hand, in x > 0 , a narrower wave crest and a lager wave trough indicate nonlinear wave propagation. More interestingly, two different wave celerities, indicated as V 1 and V 2 in Figure 11, are clearly observed in the shallow water region. As we shall see further on, V 2 corresponds to the propagation of free waves at the second order.
From the space-time matrix η ( x , t ) , considering it in a harmonic regime, we can compute different temporal modes as
η n ( x ) = 1 T 0 T η ( x , t ) e i n ω t d t ,
where T is the duration of the signal and is equal to an integer multiple of the wave period. The amplitude of the harmonics η 1 ( x ) , η 2 ( x ) and η 3 ( x ) , for the experiment carried out at ω = 4 π s 1 , is presented in Figure 12. As we observe, | η 1 | has a large modulation due to the step reflection in the region x < 0 . In the region x > 0 , we observe large attenuation due to experimental conditions, such as bottom friction or surface impurities. The second harmonic | η 2 | has a clear beating in the shallow region, suggesting the spatial interference of waves at different velocities. Eventually, the third-order harmonic | η 3 | has a non-negligible amplitude that can be interesting to take into account. However, considering that the multimodal model only solves the first and second order, we focus on these two problems, leaving the third order to a future study.
Once the harmonics are calculated, we compute a spatial Fourier transform as
η ^ n ( k ) = 1 L 0 L η n ( x ) e i k x d x ,
from which we can analyze the spatial spectral information at a given harmonic frequency n ω . As an example, regardless of the spatial spectrum of the first harmonic that does not give any additional information, we show in Figure 13 the spatial spectrum of the second harmonic η 2 ( k ) measured at ω = 4 π s 1 . Interestingly, we clearly distinguish two maxima corresponding to the free and the bound wave that propagates at this frequency. As a reference, two vertical lines have been added. First, on the left side, the wavenumber 2 k 0 + corresponds to the bound wave, which is twice the wavenumber satisfying the dispersion relation of Equation (18), and second, on the right side, the wavenumber κ 0 + corresponds to the free wave, which satisfies the dispersion relation of the second order of Equation (49). The experimental separation of free and bound waves is limited by the spectral resolution (total length of the spatial domain). Regarding the amplitude of both maxima, the dispersion of the waves represented by the relative depth of the experiment k 0 + h + determines the contribution of the free wave with respect to the bound wave. In the multimodal model, free waves are designed to ensure the matching of ϕ and x ϕ once the bound wave field is calculated. We observe that the higher the relative depth k 0 + h + is, the smoother the transition between both sides of the bound wave field ϕ 2 b is (particular solutions of the second order) and the smaller the contribution of the free waves field ϕ 2 f required to satisfy the matching conditions is. This trend shall be confirmed experimentally further on.
In order to give a general overview of the wave propagation, we show in Figure 14a,b the experimental dispersion relation at the first and second order, respectively. Here, we present the power spectrum P S ( ω , k ) (in colored contours) normalized by its maximum value corresponding to the forcing wave— P S ( ω , k ) / m a x [ P S ( ω , k ) ] —which means that all maxima are equal to 1 in Figure 14a. In this graph, we have gathered an experimental series keeping only one experiment per frequency (when the number of experiments was more than one). Every spectrum in the plane ( ω , k ) was normalized by its maximum value, which corresponds to the fundamental mode. More interestingly, at the second order, the dispersion relation of the bound and free waves follows different curves. Whereas bound waves have the same phase velocity as the fundamental mode v ϕ b = 2 ω / ( 2 k 0 ) , free waves propagate at a slower phase velocity v ϕ f = 2 ω / κ 0 . Additionally, from Figure 14b, we can observe the decay of the amplitude of the free wave as a function of the wave frequency. This decay shall be analyzed in comparison with the theoretical multimodal model.

4.3. Amplitude and Phase of Free and Bound Waves

We start by reconstructing the surface displacement of free and bound waves from the solution of the second-order potential ϕ 2 f and ϕ 2 b . These complex valued coefficients shall be compared further on with experiments. Therefore, from the dynamic boundary condition of Equation (32), we can compute the free wave as
η 2 f = 1 g ϕ 2 f t
and the bound wave as
η 2 b = 1 g ϕ 2 b t + 1 2 ϕ 1 x 2 + ϕ 1 z 2 1 g ϕ 1 t 2 ϕ 1 t z
where ϕ 2 b and ϕ 2 f are the propagating terms of the total second-order solution (see Equation (35)). The free surface wave displacement is easily obtained from Equation (67); thus, we only give an explicit solution of η 2 b . Considering an incident wave I 0 , 0 ( 1 , 1 ) ( x , t ) (as in Equation (36)), the velocity potential of the incident bound wave is
ϕ 2 b = 3 8 ω a 2 cosh 2 k 0 ( z + h ) sinh 4 ( k 0 h ) sin 2 ( k 0 x ω t ) ,
which can be replaced, together with the incident term of ϕ 1 = I 0 ( 1 ) F 0 ( 1 ) , in Equation (68) to obtain the surface displacement of the incident bound wave
η 2 b = 1 4 a 2 k 0 2 + cosh 2 k 0 h cosh k 0 h sinh 3 k 0 h cos 2 ( k 0 x ω t ) 1 2 a 2 k 0 sinh 2 k 0 h
where the last term, which is constant, corresponds to the shift in the mean water level due to the wave non-linearity (Stokes [34]).
Next, we pass to the experimental validation of the theoretical waves η 2 f and η 2 b . From the space-time Fourier transform, we can extract at the frequency of the second order ( 2 ω ), the complex values of the free and bound waves corresponding to the local maxima shown in Figure 13. The ratio between the amplitudes of free and bound waves ( | η 2 f | / | η 2 b | ) is plotted in Figure 15a as a function of the non-dimensional wavenumber k 0 + h + . The theoretical curve was calculated with the same experimental parameters: h = 0.065 m, h + = 0.02 m, a = 0.003 m and ω [ 10 , 20 ] s 1 . We observe the decay of the theoretical curve that is well verified by the experimental points that are, however, mostly below the theoretical curve. This overestimation of the theoretical ratio with respect to the experiments may be explained by the actual attenuation due to the experimental conditions. The fact that the theory only solves the second order may also overestimate the contribution of this order to the total solution. To complement, the phase difference ( θ ( η 2 f / η 2 b ) ) measured at the step position between the free and bound waves is plotted in Figure 15b. This phase difference, according to the multimodal model, should be close to π , i.e., waves have opposite phases evolving slightly with the wavenumber. However, experimentally, the fact that the relative amplitude of the free wave decays as a function of the wavenumber means the signal-to-noise ratio is not sufficiently high enough to detect the free wave phase; thus, the experimental points are only accurately measured in the range k 0 + h + < 0.8 .

4.4. Beating Length and the Influence of the Surface Tension

The scale of the experimental set-up, in particular the shallow water depth h + = 2 cm, makes the second-order wavenumber be in the order of k 100 m 1 . At this scale, the influence of the surface tension in the dynamics of wave propagation starts to be non-negligible, as recently reported by Alippi et al. [35]. This is especially the case at the second order, where the high wavenumbers enhance the effect of the surface tension on the wave propagation. Actually, the wavenumber at the second-order free waves follows the dispersion relation of gravity-capillary waves
ω 2 = g k 0 ± + σ ρ ( k 0 ± ) 3 tanh ( k 0 ± h ± ) ,
where σ is the surface tension and ρ the density. As we observe in Figure 16, one of the clearest dimensions to be measured is the beating length of the second harmonic. This distance produced by the spatial interference between the free and bound waves is equal to the inverse of the wavenumber difference and reads
L 2 = 2 π κ 0 + 2 k 0 + .
In Figure 16, the experimental measured η 2 has a beating longer than the theoretical one. This mismatch is then corrected when the wavenumber of the propagating free and bound waves takes into account the surface tension, as in Equation (71). This addition of the surface tension has been included arbitrarily in the propagating wave after the multimodal calculation, which has been computed considering only a pure gravity regime.
Eventually, in Figure 17, we gather the measurement of the beating length of the second-order L 2 for all the experiments. We plot the beating length from Equation (72) as a reference when k 0 + and κ 0 + are calculated from the dispersion relation of gravity waves in Equation (49) and when they are calculated from the dispersion relation of gravity capillary waves in Equation (71). As we see, the beating length follows different curves depending on the wavenumber. Whereas for small waveumbers the beating length seems to be dominated by the gravity regime, for large wavenumbers, the influence of the surface tension in the beating length is clear.

5. Conclusions

In this article, we have revisited the multimodal model of the propagation of nonlinear waves over a submerged rectangular step proposed by Massel [6]. In this work, we focused on the problem of convergence that the model shows at the second order, where the bound evanescent waves, when calculated in the source term (right hand) of the kinematic and dynamic surface boundary conditions, produce a growth of the source term at a power n 4 that is higher than the decay exponent of the reflection and transmission coefficients of the first order, which are supposed to compensate for the divergence in the source term. This problem of convergence, although directly detectable when the model is correctly implemented, was not reported by Massel [6], apparently because the calculation was performed with only 20 modes. Nowadays, the computational power means the numerical computation can use a large number of modes at a low cost, which needs a robust implementation to ensures a convergent model independently of the number of modes. Therefore, by analyzing in detail the terms with the fastest growth in the source term that dominate the divergence of the second-order solution, we found that a truncation of the source term at a number of modes smaller than the first order produces convergence of the second-order problem, resulting in a smooth field in both potential ϕ 2 and horizontal velocity x ϕ 2 . We found that the threshold of the number of modes in the source term α N ± , with α 0.5 , ensures the convergence of the second-order solution. We have verified that this truncation works independently of the input parameters of the problem (geometry and wave characteristics) and has an error in the multimodal matching that decreases as a function of N ± . This convergence analysis can be useful to researchers interested in solving rectangular submerged obstacles in nonlinear wave problems.
The second part of the article shows the space-time-resolved measurement of the nonlinear wave propagation. This kind of measurement allows computing the complex coefficient (amplitude and phase) of the propagating waves at the first and second order. In particular, we were interested in the propagation of free and bound waves at the second order, where we have experimentally verified that the relative contribution of free waves decreases as a function of the non-dimensional wavenumber k 0 + h + (the system becomes more linear), and the phase between free and bound waves is the opposite at the step. Additionally, due to the scale of the experimental set-up, surface tension effects are non-negligible at the second order, where the wavenumber is in the order of k 100 m 1 . A systematic measurement of the beating length of the second harmonic reveals that at high wavenumbers ( k 0 + h + > 0.7 ), the surface tension should be considered. An interesting continuation of this work may consider the development of the multimodal model with surface tension, which would make the solution more suitable to small-scale phenomena.

Author Contributions

Conceptualization, E.M., A.M., V.P. and P.P.; Methodology, E.M., A.M., V.P. and P.P.; Numerical calculation, E.M.; Investigation, E.M.; Data curation, E.M.; Writing, E.M., A.M., V.P. and P.P.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

E.M. acknowledges the support of Agencia Nacional de Investigacion y Desarrollo (ANID) Chile Becas Chile Doctorado.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. First Order

Appendix A.1. Solution of the System

The system of Equations (27) and (29) can be solved for R n ( 1 )
n R n ( 1 ) x A n ( 1 ) N n ( 1 ) δ n , γ m A n ( 1 ) x B m ( 1 ) L n , m ( 1 ) L γ , m ( 1 ) M m ( 1 ) B m ( 1 ) = x I 0 ( 1 ) N 0 ( 1 ) δ 0 , γ + m I 0 ( 1 ) x B m ( 1 ) L 0 , m ( 1 ) L γ , m ( 1 ) M m ( 1 ) B m ( 1 )
where δ k , γ is the Kronecker delta, and the following notations were used to express the integrals:
L n , m ( 1 ) = h + 0 F n ( 1 ) G m ( 1 )
M m ( 1 ) = h + 0 ( G m ( 1 ) ) 2
N n ( 1 ) = h 0 ( F n ( 1 ) ) 2
The transmission coefficients T m ( 1 ) can be easily obtained from Equation (27).

Appendix B. Second Order

Appendix B.1. Constants

The bound wave constant C n , n ( 1 , 1 ) obtained from the source term by replacing Equation (21) in the right-hand side of Equation (31) and performing lengthy algebraic manipulations
C n , n ( 1 , 1 ) = i g 2 a 2 2 ω ( k n ) 2 + 4 k n k n + ( k n ) 2 ( k n ) 2 tanh 2 k n h ( k n ) 2 tanh 2 k n h 4 k n k n tanh k n h tanh k n h 4 ω 2 g ( k n + k n ) tanh ( k n + k n ) h .
The constant C 0 , n ( 1 , 1 ) of the second-order evanescent terms resulting from the multiplication of the Stokes incident harmonic with the evanescent reflected bound waves takes the following form
C 0 , n ( 1 , 1 ) = i g 2 a 2 2 ω ( k n ) 2 4 k 0 k n + ( k n ) 2 ( k 0 ) 2 tanh 2 k 0 h ( k n ) 2 tanh 2 k n h 4 k 0 k n tanh k 0 h tanh k n h 4 ω 2 g ( k 0 k n ) tanh ( k 0 k n ) h .

Appendix B.2. Solution of the System

The system of Equations (58) and (60) can be solved for R p ( 2 ) :
p R p ( 2 ) x A p ( 2 ) N p ( 2 ) δ p , γ q A p ( 2 ) x B q ( 2 ) L p , q ( 2 ) L γ , q ( 2 ) M q ( 2 ) B q ( 2 ) + c . c . = m m T m ( 1 ) T m ( 1 ) x B m , m ( 1 , 1 ) h + 0 G m , m ( 1 , 1 ) F γ ( 2 ) x I 0 , 0 ( 1 , 1 ) h 0 F 0 , 0 ( 1 , 1 ) F γ ( 2 ) n n R n ( 1 ) R n ( 1 ) x A n , n ( 1 , 1 ) h 0 F n , n ( 1 , 1 ) F γ ( 2 ) n R n ( 1 ) x D 0 , n ( 1 , 1 ) h 0 H 0 , n ( 1 , 1 ) F γ ( 2 ) + q x B q ( 2 ) L γ , q ( 2 ) M q ( 2 ) B q ( 2 ) I 0 , 0 ( 1 , 1 ) h + 0 F 0 , 0 ( 1 , 1 ) G q ( 2 ) + n n R n ( 1 ) R n ( 1 ) A n , n ( 1 , 1 ) h + 0 F n , n ( 1 , 1 ) G q ( 2 ) + n R n ( 1 ) D 0 , n ( 1 , 1 ) h + 0 H 0 , n ( 1 , 1 ) G q ( 2 ) m m T m ( 1 ) T m ( 1 ) B m , m ( 1 , 1 ) h + 0 G m , m ( 1 , 1 ) G q ( 2 ) + c . c .
Here, we have expressed some known integrals with the following notation
L p , q ( 2 ) = h + 0 F p ( 2 ) G q ( 2 ) ,
M q ( 2 ) = h + 0 ( G q ( 2 ) ) 2 ,
N p ( 2 ) = h 0 ( F p ( 2 ) ) 2 .
The transmission coefficients T q ( 2 ) can be easily obtained from Equation (58).

References

  1. Li, Y.; Draycott, S.; Zheng, Y.; Lin, Z.; Adcock, T.A.; Van Den Bremer, T.S. Why rogue waves occur atop abrupt depth transitions. J. Fluid Mech. 2021, 919, R5. [Google Scholar] [CrossRef]
  2. Mei, C.C.; Black, J.L. Scattering of surface waves by rectagular obstacles in waters of finite depth. J. Fluid Mech. 1969, 38, 499–511. [Google Scholar] [CrossRef]
  3. Miles, J.W. Surface-wave scattering matrix for a shelf. J. Fluid Mech. 1967, 28, 755–767. [Google Scholar] [CrossRef]
  4. Newman, J. Propagation of water waves over an infinite step. J. Fluid Mech. 1965, 23, 399–415. [Google Scholar] [CrossRef]
  5. Grue, J. Nonlinear water waves at a submerged obstacle or bottom topography. J. Fluid Mech. 1992, 244, 455–476. [Google Scholar] [CrossRef] [Green Version]
  6. Massel, S.R. Harmonic generation by waves propagating over submerged step. Coast. Eng. 1983, 7, 357–380. [Google Scholar] [CrossRef]
  7. Belibassakis, K.A.; Athanassoulis, G.A. Extension of second-order Stokes theory to variable bathymetry. J. Fluids Mech. 2002, 464, 35–80. [Google Scholar] [CrossRef] [Green Version]
  8. Belibassakis, K.A.; Athanassoulis, G.A. A coupled-mode system with application to nonlinear water waves propagating in finite water depth and in variable bathymetry regions. Coast. Eng. 2011, 58, 337–350. [Google Scholar] [CrossRef]
  9. Athanassoulis, G.A.; Belibassakis, K.A. A consistent coupled-mode theory for the propagation of small amplitude water waves over variable bathymetry regions. J. Fluids Mech. 1999, 389, 275–301. [Google Scholar] [CrossRef] [Green Version]
  10. Rhee, J.P. On the transmission of water waves over a shelf. Appl. Ocean Res. 1997, 19, 161–169. [Google Scholar] [CrossRef]
  11. Porter, R.; Porter, D. Approximations to the scattering of water waves by steep topography. J. Fluids Mech. 2006, 562, 279–302. [Google Scholar] [CrossRef]
  12. Li, Y.; Zheng, Y.; Lin, Z.; Adcock, T.A.; van den Bremer, T.S. Surface wavepackets subject to an abrupt depth change. Part 1. Second-order theory. J. Fluid Mech. 2021, 915, A71. [Google Scholar] [CrossRef]
  13. Ducrozet, G.; Slunyaev, A.; Stepanyants, Y. Transformation of envelope solitons on a bottom step. Phys. Fluids 2021, 33, 066606. [Google Scholar] [CrossRef]
  14. Bryant, P.J. Periodic waves in shallow water. J. Fluid Mech. 1972, 59, 625–644. [Google Scholar] [CrossRef]
  15. Huang, C.J.; Dong, C.M. Wave deformation and vortex generation in water waves propagating over a submerged dike. Coast. Eng. 1999, 37, 123–148. [Google Scholar] [CrossRef]
  16. Chapalain, G.; Cointe, R.; Temperville, A. Observed and modeled resonantly interacting progressive water-waves. Coast. Eng. 1992, 16, 267–300. [Google Scholar] [CrossRef]
  17. Beji, S.; Battjes, J.A. Experimental investigation of wave propagation over a bar. Coast. Eng. 1992, 19, 151–162. [Google Scholar] [CrossRef]
  18. Li, F.; Ting, C. Separation of free and bound harmonics in waves. Coast. Eng. 2012, 67, 29–40. [Google Scholar] [CrossRef]
  19. Benoit, M.; Raoult, C.; Yates, M. Fully nonlinear and dispersive modeling of surf zone waves: Non-breaking tests. Coast. Eng. Proc. 2014, 1, 15. [Google Scholar] [CrossRef]
  20. Ohyama, T.; Nadaoka, K. Transformation of a nonlinear wave train passing over a submerged shelf without breaking. Coast. Eng. 1993, 24, 1–22. [Google Scholar] [CrossRef]
  21. Brossard, J.; Chagdali, M. Experimental investigation of the harmonic generation by waves over a submerged plate. Coast. Eng. 2001, 42, 277–290. [Google Scholar] [CrossRef]
  22. Brossard, J.; Perret, G.; Blonce, L.; Diedhiou, A. Higher harmonics induced induced by a submerged horizontal plate and a submerged rectangular step in a wave flume. Coast. Eng. 2009, 56, 11–22. [Google Scholar] [CrossRef]
  23. Ting, C.L.; Chao, W.T.; Young, C.C. Experimental investigation of nonliear regular waves transformation over submerged step: Harmonic generation and wave height modulation. Coast. Eng. 2016, 117, 19–31. [Google Scholar] [CrossRef]
  24. Ohyama, T.; Kioka, W.; Tada, A. Applicability of numerical models to nonlinear dispersive waves. Coast. Eng. 1995, 24, 297–313. [Google Scholar] [CrossRef]
  25. Takeda, M.; Ina, H.; Kobayashi, S. Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry. JosA 1982, 72, 156–160. [Google Scholar] [CrossRef]
  26. Takeda, M.; Mutoh, K. Fourier transform profilometry for the automatic measurement of 3-D object shapes. Appl. Opt. 1983, 22, 3977–3982. [Google Scholar] [CrossRef] [PubMed]
  27. Cobelli, P.; Maurel, A.; Pagneux, V.; Petitjeans, P. Global measurement on water waves by Fourier transform profilometry. Exp. Fluids 2009, 46, 1037–1047. [Google Scholar] [CrossRef]
  28. Hsu, J.R.C.; Tsuchiya, Y.; Silvester, R. Third order approximation to short crested waves. J. Fluids Mech. 1979, 90, 179–196. [Google Scholar] [CrossRef]
  29. Mei, C.C.; Stiassnie, M.; Yue, D.K.P. Theory and Applications of Ocean Surface Waves. Part 1: Linear Aspects; World Scientific: Singapore, 2005; Volume 23. [Google Scholar]
  30. Sommerfeld, A. Partial Differential Equations in Physics; Academic Press: Cambridge, MA, USA, 1949; Volume 1. [Google Scholar]
  31. Przadka, A.; Cabane, B.; Pagneux, V.; Maurel, A.; Petitjeans, P. Fourier transform profilometry for water waves: How to achieve clean water attenuation with diffusive reflection at the water surface? Exp. Fluids 2012, 52, 519–527. [Google Scholar] [CrossRef]
  32. Monsalve, E.; Maurel, A.; Petitjeans, P.; Pagneux, V. Perfect absorption of water waves by linear or nonlinear critical coupling. Appl. Phys. Lett. 2019, 114, 013901. [Google Scholar] [CrossRef] [Green Version]
  33. Monsalve, E.; Maurel, A.; Pagneux, V.; Petitjeans, P. Space-time-resolved measurements of the effect of pinned contact line on the dispersion relation of water waves. Phys. Rev. Fluids 2022, 7, 014802. [Google Scholar] [CrossRef]
  34. Stokes, G.G. On the theory of the oscillatory waves. Trans. Camb. Phil. Soc. 1847, 8, 441–455. [Google Scholar]
  35. Alippi, A.; Bettucci, A.; Germano, M. Model and experiments for resonant generation of second harmonic capillary–gravity waves. Phys. D Nonlinear Phenom. 2019, 396, 12–17. [Google Scholar] [CrossRef]
Figure 1. Scheme of the propagation problem of a nonlinear wave over a submerged rectangular obstacle. The incident wave comes from the left side.
Figure 1. Scheme of the propagation problem of a nonlinear wave over a submerged rectangular obstacle. The incident wave comes from the left side.
Fluids 07 00145 g001
Figure 2. (a) Velocity potential at the first-order ϕ 1 near the step. (b) Horizontal velocity at the first-order x ϕ 1 near the step.
Figure 2. (a) Velocity potential at the first-order ϕ 1 near the step. (b) Horizontal velocity at the first-order x ϕ 1 near the step.
Fluids 07 00145 g002
Figure 3. (a) Vertical profile at x = 0 of the velocity potential at the first-order ϕ 1 . (b) Vertical profile at x = 0 of the horizontal velocity at the first-order x ϕ 1 .
Figure 3. (a) Vertical profile at x = 0 of the velocity potential at the first-order ϕ 1 . (b) Vertical profile at x = 0 of the horizontal velocity at the first-order x ϕ 1 .
Fluids 07 00145 g003
Figure 4. Non-convergent solution of the second-order problem, α = 1 (without truncation of forcing evanescent waves). (a) Velocity potential at the second-order ϕ 2 near the step. (b) Horizontal velocity at the second-order x ϕ 2 near the step.
Figure 4. Non-convergent solution of the second-order problem, α = 1 (without truncation of forcing evanescent waves). (a) Velocity potential at the second-order ϕ 2 near the step. (b) Horizontal velocity at the second-order x ϕ 2 near the step.
Fluids 07 00145 g004
Figure 5. (a) Reflection coefficient R n ( 1 ) of the first order for N = 100 modes. (b) Evanescent forcing coefficients C n , n ( 1 , 1 ) at the second order for N = 100 modes.
Figure 5. (a) Reflection coefficient R n ( 1 ) of the first order for N = 100 modes. (b) Evanescent forcing coefficients C n , n ( 1 , 1 ) at the second order for N = 100 modes.
Fluids 07 00145 g005
Figure 6. (a) Numerical error of the multimodal matching for the second-order functions ϕ 2 and x ϕ 2 as a function of the truncation parameter α . (b) Reflection and transmission coefficients of the free waves at the second order as a function of the truncation parameter α .
Figure 6. (a) Numerical error of the multimodal matching for the second-order functions ϕ 2 and x ϕ 2 as a function of the truncation parameter α . (b) Reflection and transmission coefficients of the free waves at the second order as a function of the truncation parameter α .
Fluids 07 00145 g006
Figure 7. Convergent solution of the second-order problem for α = 0.5 (with truncation of forcing evanescent waves). (a) Velocity potential ϕ 2 at the second order near the step. (b) Horizontal velocity x ϕ 2 at the second order near the step.
Figure 7. Convergent solution of the second-order problem for α = 0.5 (with truncation of forcing evanescent waves). (a) Velocity potential ϕ 2 at the second order near the step. (b) Horizontal velocity x ϕ 2 at the second order near the step.
Fluids 07 00145 g007
Figure 8. (a) Vertical profile at x = 0 of the velocity potential at the second order ϕ 2 ( 0 , z ) . (b) Vertical profile at x = 0 of the horizontal velocity at the first-order x ϕ 2 ( 0 , z ) .
Figure 8. (a) Vertical profile at x = 0 of the velocity potential at the second order ϕ 2 ( 0 , z ) . (b) Vertical profile at x = 0 of the horizontal velocity at the first-order x ϕ 2 ( 0 , z ) .
Fluids 07 00145 g008
Figure 9. Numerical error of the multimodal matching as a function of the number of modes for the case α = 0.5 . (a) Numerical error of the velocity potentials ϕ 1 and ϕ 2 . (b) Numerical error of the horizontal velocities x ϕ 1 and x ϕ 2 .
Figure 9. Numerical error of the multimodal matching as a function of the number of modes for the case α = 0.5 . (a) Numerical error of the velocity potentials ϕ 1 and ϕ 2 . (b) Numerical error of the horizontal velocities x ϕ 1 and x ϕ 2 .
Fluids 07 00145 g009
Figure 10. Experimental set-up. (a) Sketch of the top and side views of the experimental tank. The measurement area is highlighted in green. (b) Experimental tank with the sinusoidal pattern of the FTP technique projected onto the diffusive surface. (c) Typical experimentally measured field of a nonlinear-wave.
Figure 10. Experimental set-up. (a) Sketch of the top and side views of the experimental tank. The measurement area is highlighted in green. (b) Experimental tank with the sinusoidal pattern of the FTP technique projected onto the diffusive surface. (c) Typical experimentally measured field of a nonlinear-wave.
Fluids 07 00145 g010
Figure 11. Space-time representation of the wave propagation at ω = 4 π s 1 after transverse average in the y direction. Dashed and dashed-dotted lines shows two different celerities in the shallow water region.
Figure 11. Space-time representation of the wave propagation at ω = 4 π s 1 after transverse average in the y direction. Dashed and dashed-dotted lines shows two different celerities in the shallow water region.
Fluids 07 00145 g011
Figure 12. Amplitude of the temporal harmonics | η n ( x ) with n = 1 , 2 , 3 measured with Fourier transform profilometry at a frequency ω = 4 π s 1 .
Figure 12. Amplitude of the temporal harmonics | η n ( x ) with n = 1 , 2 , 3 measured with Fourier transform profilometry at a frequency ω = 4 π s 1 .
Fluids 07 00145 g012
Figure 13. Amplitude of the second harmonic η 2 ( k ) as a function of wave number. Non-dimensional wavenumber of the linear mode is k 0 + h + = 0.73. Vertical lines indicate the wavenumber expected from the dispersion relation of water waves in Equation (49).
Figure 13. Amplitude of the second harmonic η 2 ( k ) as a function of wave number. Non-dimensional wavenumber of the linear mode is k 0 + h + = 0.73. Vertical lines indicate the wavenumber expected from the dispersion relation of water waves in Equation (49).
Fluids 07 00145 g013
Figure 14. Experimental dispersion relation in the plane ( ω , k ) . (a) Dispersion relation of the linear mode η ( ω ) . (b) Dispersion relation of the second-order η ( 2 ω ) .
Figure 14. Experimental dispersion relation in the plane ( ω , k ) . (a) Dispersion relation of the linear mode η ( ω ) . (b) Dispersion relation of the second-order η ( 2 ω ) .
Fluids 07 00145 g014
Figure 15. Experimental measurement and theoretical model of the second-order free and bound waves. (a) Amplitude ratio between free and bound waves | η 2 f | / | η 2 b | . (b) Phase different between free and bound waves measured at the step position θ ( η 2 f | / | η 2 b ) .
Figure 15. Experimental measurement and theoretical model of the second-order free and bound waves. (a) Amplitude ratio between free and bound waves | η 2 f | / | η 2 b | . (b) Phase different between free and bound waves measured at the step position θ ( η 2 f | / | η 2 b ) .
Fluids 07 00145 g015
Figure 16. Comparison of the experimentally measured amplitude of the first (blue solid line) and second (red solid line) harmonics ( | η 1 | , | η 2 | ) with the propagating wave amplitudes computed in the multimodal model. At the second order, the theoretical beating is computed with (black dashed line) and without (magenta dash-dotted line) surface tension. Wave frequency ω = 4 π s 1 .
Figure 16. Comparison of the experimentally measured amplitude of the first (blue solid line) and second (red solid line) harmonics ( | η 1 | , | η 2 | ) with the propagating wave amplitudes computed in the multimodal model. At the second order, the theoretical beating is computed with (black dashed line) and without (magenta dash-dotted line) surface tension. Wave frequency ω = 4 π s 1 .
Fluids 07 00145 g016
Figure 17. Beating length of the second harmonic η 2 as a function of the non-dimensional wavenumber k 0 + h + . Two theoretical curves of L 2 are calculated from Equation (72): First, the blue dash-dotted line shows L 2 calculated with wavenumbers from Equation (71), and second, the red dashed line shows L 2 with wavenumbers from Equation (18).
Figure 17. Beating length of the second harmonic η 2 as a function of the non-dimensional wavenumber k 0 + h + . Two theoretical curves of L 2 are calculated from Equation (72): First, the blue dash-dotted line shows L 2 calculated with wavenumbers from Equation (71), and second, the red dashed line shows L 2 with wavenumbers from Equation (18).
Fluids 07 00145 g017
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Monsalve, E.; Maurel, A.; Pagneux, V.; Petitjeans, P. Nonlinear Waves Passing over Rectangular Obstacles: Multimodal Method and Experimental Validation. Fluids 2022, 7, 145. https://doi.org/10.3390/fluids7050145

AMA Style

Monsalve E, Maurel A, Pagneux V, Petitjeans P. Nonlinear Waves Passing over Rectangular Obstacles: Multimodal Method and Experimental Validation. Fluids. 2022; 7(5):145. https://doi.org/10.3390/fluids7050145

Chicago/Turabian Style

Monsalve, Eduardo, Agnès Maurel, Vincent Pagneux, and Philippe Petitjeans. 2022. "Nonlinear Waves Passing over Rectangular Obstacles: Multimodal Method and Experimental Validation" Fluids 7, no. 5: 145. https://doi.org/10.3390/fluids7050145

APA Style

Monsalve, E., Maurel, A., Pagneux, V., & Petitjeans, P. (2022). Nonlinear Waves Passing over Rectangular Obstacles: Multimodal Method and Experimental Validation. Fluids, 7(5), 145. https://doi.org/10.3390/fluids7050145

Article Metrics

Back to TopTop