Next Article in Journal
Two-Scale Investigation of the Retention Behavior of a Well-Graded Mixed Soil
Next Article in Special Issue
Numerical Investigation of the Effectiveness of Vegetation-Embankment Hybrid Structures for Tsunami Mitigation Introduced after the 2011 Tsunami
Previous Article in Journal
The Geo-Hydro-Mechanical Properties of a Turbiditic Formation as Internal Factors of Slope Failure Processes
Previous Article in Special Issue
The 2011 Tohoku Tsunami from the Sky: A Review on the Evolution of Artificial Intelligence Methods for Damage Assessment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical and Numerical Investigations Applied to Study the Reflections and Transmissions of a Rectangular Breakwater Placed at the Bottom of a Wave Tank

1
State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116023, China
2
Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, LAMA, 73000 Chambéry, France
3
Faculty of Civil Engineering, University of Bab-Ezzouar, Algiers 16000, Algeria
4
Faculty of Civil Engineering, Technical University of Kosice, Vysokoškolská 4, 042 00 Košice, Slovakia
*
Author to whom correspondence should be addressed.
Geosciences 2021, 11(10), 430; https://doi.org/10.3390/geosciences11100430
Submission received: 25 August 2021 / Revised: 12 October 2021 / Accepted: 14 October 2021 / Published: 18 October 2021

Abstract

:
The purpose of the work presented in this paper is to study the reflection and transmission coefficients resulting from the interactions of regular waves with a rectangular breakwater sited at the bottom of a tank. The present investigation is devoted to the analysis of the reflection and transmission coefficients within the framework of linearized potential flow theory using two methods, a numerical method based on the improved version of the meshless singular boundary method, and the analytical approach within the plane wave model. The numerical method is first validated by studying the accuracy of the numerical computations with respect to the number of boundary nodes and the location of the vertical boundaries of the computational domain, for different immersion ratios (h/d) and different relative lengths (w/d) of the obstacle. To assess the limitations of the analytical approach, a comparison analysis is carried out between the analytical and numerical results. To improve the calculations and the effectiveness of the analytical model, slight adjustments are made to the analytical procedure, which is termed here the corrected analytical plane wave model. Finally, the effects of the immersion ratio (h/d) and the relative length (w/d) of the obstacle on the reflection and transmission coefficients are computed using the three methods, and discussed for several wave and structural conditions.

1. Introduction

The coast is a site where several phenomena can appear and affect coastal structures, such as sea-level rise caused by climate change, erosion due to wave actions, and the decrease of fluvial sediment supply caused by the construction of dams, etc. Consequently, researchers, engineers, and scientists have shown that the reflection of swell is one of the convincing and relevant solutions to overcome problems of erosion, surges, marine submersions, and all phenomena to which these structures can be exposed. Therefore, the maximization of wave reflection has become extremely important in coastal engineering, and arouses the interest of many researchers in the literature to investigate the numerical, experimental, and theoretical methods that can study the reflection of wave–structure interactions.
Considerable research has been carried out in the past on different types of breakwaters. Nevertheless, the operating conditions of breakwaters are very challenging and, hence, there is still a need for additional research to fill in the gaps. A short account of what has been delivered up till now is given here. Dean [1] studied the effect of the wave amplitudes on the reflection of surface waves by a submerged plane barrier. Takano [2] evaluated the passage effect of waves propagating under a rectangular breakwater. Patarapanich [3] studied the wave reflection and transmission by a submerged thin horizontal plate using the finite element method (FEM). Using the matched asymptotic method, Liu and Jiankng [4] investigated the transmitted wave intensity through a submerged slit on a vertical barrier. On the other hand, using a parametric experimental design, Stamos et al. [5] compared the reflection and transmission coefficients resulting from the interaction of waves with a variety of submerged water-filled breakwater models of hemi-cylindrical and rectangular shapes. Molin et al. [6] carried out laboratory experiments to investigate the interaction of waves with a rigid vertical plate. Shortly afterwards, Lui et al. [7] examined the Bragg reflections of water waves by multiple submerged semi-circular breakwaters.
Further, the submerged rectangular step is a structure mostly used as a breakwater to protect shorelines by diminishing the destructive effects of the wave actions, reducing the erosion of coasts, and protecting coastal structures from damage [8]. Recently, the rectangular submerged breakwater has started to receive more attention compared to the traditional emerged structures, due to its attractive aesthetics and ability to allow water circulation as well as the passage of fish. Several experimental, analytical, and numerical studies have been devoted to studying the reflection and transmission of swells by this type of structure. Mei and Black [9] studied the problem of scattering properties for bottom and surface obstacles using the variational method. Massel [10] investigated the interaction of waves with an infinite- and finite-length rectangular submerged breakwater. Andrew et al. [11] presented an experimental and numerical study based on the boundary element method (BEM) to investigate the propagation of waves over a submerged impermeable obstacle of a rectangular cross section. Recently, Szmidt [12] studied numerically using the finite difference method (FDM), the interaction of waves with a rectangular breakwater fixed at the bottom of numerical wave tank (NWT), and estimated the efficiency of the breakwater in protecting sea shelf zones from open sea waves.
We propose to investigate, in this paper, the capabilities of the improved version of the meshless singular boundary method (ISBM) [13,14,15] to analyze the reflection and transmission coefficients resulting from the interactions of regular waves with a rectangular breakwater sited at the bottom of a tank. The method is validated by studying the accuracy of the numerical results with respect to the number of boundary nodes and the location of the vertical boundaries of the computational domain for different immersion ratios (h/d) and different relative lengths (w/d) of the obstacle. Further, the analytical reflection and transmission coefficients within the plane wave model (see Appendix A) are compared with the results of the numerical model and discussed for several wave and structural conditions. To improve the shortfalls of the analytical model, slight modifications are introduced to the analytical procedure, which is termed here the corrected analytical plane wave model. Finally, a general discussion is made to highlight the strengths and limitations of the corrected plane wave model.
This paper is divided into five main sections. After presenting the introduction of the work and stating the objective of this research in Section 1, the formulation of the problem and the numerical method used in this work are presented in Section 2 and Section 3. In Section 4, the error sensitivity indicator, and comparisons of the numerical and analytical results and analysis, are presented to study the interactions of regular waves with a rectangular breakwater sited at the bottom of a tank. Finally, some conclusions and perspectives are illustrated in Section 5.

2. Formulation of the Problem

We consider, in this study, a submerged single impermeable rectangular step (breakwater) which is placed at the bottom, as shown in Figure 1.
The idealized geometry of the two-dimensional (2D) problem in a Cartesian system (x-y) is shown in Figure 2. Regular waves of small amplitude a, period T, and wavelength L impinge from the left in water of depth d. Assuming an irrotational flow and incompressible fluid motion, the problem is formulated using a velocity potential ϕ ( x , y , t ) = R e [ ϕ ( x , y ) · exp ( i σ t ) ] where Re denotes the real part, ϕ ( x , y ) is the time independent spatial velocity potential, i = 1   ,   σ = 2 π / T is the wave angular frequency, and t is the time. The wave number k = 2 π / L is the solution of the dispersion relation g k · tanh ( k d ) = σ 2 , where g is the gravitational acceleration. The wave field is totally specified if the two-dimensional velocity potential ϕ is known.
The breakwater is described by the immersion ratio h/d and relative length (w/d), where d is the water depth in the absence of the obstacle, h is the water depth above the obstacle, and w is the length of the obstacle.
The total fluid domain is divided into three regions as shown in Figure 2. Region I at (−∞) is the region where the waves are incoming (inflow), and region III at (+∞) is where the waves are transmitted (outflow). Region II is between regions I and III, and is delimited by the rigid (impermeable) walls of the breakwater ( Γ b 1 , Γ b 2 , and Γ b 3 ), the free surface boundary Γ f , the seabed boundary Γ s , and the radiation boundaries Γ and Γ + , respectively, of the inflow and outflow regions. The spatial velocity potential ϕ satisfies the following conditions:
2 ϕ x 2 + 2 ϕ y 2 = 0   in   the   fluid   region   II  
ϕ n σ 2 g · ϕ = 0   at   y = d   ( free   surface   boundary   Γ f )  
ϕ n = 0   at   y = 0   ( seabed   boundary   Γ s )  
ϕ n = 0   ( x , y )     Γ b 1   Γ b 2   Γ b 3   ( breakwater   boundaries   Γ b )  
where n is the normal to the boundary pointing out of the flow region, and Γ b = Γ b 1   Γ b 2   Γ b 3 denotes the total rigid (impermeable) boundary of the breakwater.
The radiation conditions at the inflow and outflow regions are expressed as
( ϕ ϕ I ) n i · k · ( ϕ ϕ I ) = 0   radiation   condition   at   x   ( boundary   Γ )  
( ϕ ) n i · k · ( ϕ ) = 0   radiation   condition   at   x +   ( boundary   Γ + )  
where ϕ I   is the incident velocity potential.
The radiation conditions in the infinite strip problem are treated by transferring the far field potentials at two fictitious vertical boundaries at finite distances x = x r and x = + x r , representing, respectively, the left boundary Γ and the right boundary Γ + of the fluid domain. The analytical series at these boundaries are given by:
ϕ = ϕ I + A · cosh ( k · y ) sinh ( k · d ) e i k · ( x + x r )   and     ϕ   n =   ϕ   x   for   x = x r   ( boundary   Γ )  
ϕ + = A + · cosh ( k · y ) sinh ( k · d ) e i k · ( x x r )   and     ϕ +   n =   ϕ +   x   for   x = + x r   ( boundary   Γ + )  
where A and A + are unknown complex coefficients to be determined. The disturbances are guaranteed to be out-going waves only (see for example [16,17]). The incident velocity potential is defined as:
ϕ I = a · L T · cosh ( k · y ) sinh ( k · d ) e i k · ( x + x r )
The special matching conditions at the interfaces Γ and Γ + of the flow regions ensure a smooth transfer of the mass flow from one region to the next. Once the potentials ϕ and ϕ + are calculated by satisfying the radiation boundary conditions of Equations (5) and (6), they are matched to those of Equations (7) and (8), then the unknown coefficients A and A + are evaluated following the method of Yueh and Chuang [18]:
A = ( a · L T ) + k N 0 · cosh ( k · d ) · 0 d ϕ ( x r , y ) · cosh ( k · y ) · d y
A + = k N 0 · cosh ( k · d ) · 0 d ϕ + ( + x r , y ) · cosh ( k · y ) · d y
where N 0 = 1 2 ( 1 + 2 k · d sinh ( 2 k · d ) )   .
The reflection and transmission coefficients (R and T r ) are determined from the following expressions (see [16,17] and further details in Appendix B):
R = | A 0 | T a · L   and   T r = | A 0 + | T a · L

3. Numerical Solution by the ISBM

For the numerical solution the total boundary of the whole computational domain is discretized as shown in Figure 3 for the single breakwater.
In the ISBM, the nodal values of the potentials and their fluxes are expressed as linear combinations of fundamental solutions and their derivatives [13,14,15],
ϕ ( x i ) = j = 1 , i j N α j · Q ( x i , s j ) + α i · ϕ i i
q ( x i ) = ϕ ( x i ) n x i = j = 1 , i j N α j · Q ( x i , s j ) n x i + α i · q i i
where α j are unknown coefficients to be determined, x i and s j are, respectively, the collocation points ( x i , y i ) and the source points ( x j , y j ) , and N is the total number of points. ϕ ( x i ) and q ( x i ) are the Dirichlet and Neumann values, and n x i is the normal at the collocation point x i . The coefficients ϕ i i and q i i are source intensity factors corresponding, respectively, to the fundamental solution and its derivative. Q ( x i , s j ) is the fundamental solution of the 2D Laplace equation. It depends only on the Euclidean distance r i j = | x i s j | between the collocation points x i and the source points s j , i.e., r i j = ( x i x j ) 2 + ( y i y j ) 2 , and is given together with its normal derivative as:
Q ( x i , s j ) = 1 2 π ln   ( 1 r i j )
Q ( x i , s j ) n x i = [ ( n x x i ) · ( x i x j ) + ( n y x i ) · ( y i y j ) ] 2 π ( r i j ) 2
n x x i and n y x i are the component of the normal at the collocation point x i .
The coefficients ϕ i i and q i i are the diagonal elements of the ISBM interpolation matrices. They arise when the collocation points and the source points coincide ( x i = s j ). Direct evaluation of these coefficients is unfeasible because of the singularities inherent in the fundamental solution and its derivative. In this study, the coefficients ϕ i i are evaluated simply by the integration of the fundamental solution on line segments leading to a simple analytical expression as [19,20]:
ϕ i i = 1 i Γ s Q ( x i , s ) d Γ s = 1 2 π [ ln ( 2 i ) + 1 ]
For the coefficient q i i , a simple expression is derived by Gu [14], using a regularization process of subtracting and adding-back to remove singularities:
q i i = 1 i j = 1 , i j N j · Q ( x i , s j ) n s j
where i and j are the half distances, respectively, between the collocations points x i 1 and x i + 1 , and the source points s j 1 and s j + 1 . n s j is the normal at the source point s j .
The boundary conditions given by Equations (2)–(6) are satisfied by a linear combination of Equations (13) and (14). The discretization process leads to:
For nodes x i     Γ f (free surface boundary):
j = 1 , i j N α j · ( Q ( x i , s j ) n x i σ 2 g · Q ( x i , s j ) ) + α i · ( q i i σ 2 g · ϕ i i ) = 0
For nodes x i     Γ (the radiation boundary at x = x r ):
j = 1 , i j N α j · ( Q ( x i , s j ) n x i i · k · Q ( x i , s j ) ) + α i · ( q i i i · k · ϕ i i ) = ϕ I ( x i ) n x i i · k · ϕ I ( x i )
For nodes x i     Γ + (the radiation boundary at x = + x r ):
j = 1 , i j N α j · ( Q ( x i , s j ) n x i i · k · Q ( x i , s j ) ) + α i · ( q i i i · k · ϕ i i ) = 0
For nodes x i     Γ s and Γ b (seabed and breakwater boundaries):
j = 1 , i j N α j · ( Q ( x i , s j ) n x i ) + α i · q i i = 0
The resulting discretized Equations (19)–(22) are written in a more compact matrix form as:
[   H i j ] N   x   N   · {   α i } N = { B i } N   i , j = 1 , 2 , , N  
where N is the total number of nodes on the whole domain boundaries, e.g., N = N f + N + N + + N b + N s , where N f , N , N + , N b , and N s are the number of nodes respectively on the boundaries Γ f , Γ , Γ + , Γ b , and Γ s . The algebraic system of equations expressed by Equation (23) is solved numerically using a Gaussian elimination algorithm to yield the vector of unknowns {   α i } . The potential and its derivative at the nodes are then computed using Equations (13) and (14).

4. Results and Discussion

The results of the present investigation are divided into two separate sections. The first section is dedicated to studying the sensitivity of the numerical results of the ISBM to ensure the validity of the proposed numerical wave model. In the Section 2, comparison is made between the results of the ISBM, the analytical approach within the plane wave model (see Appendix A), and the corrected plane wave model.

4.1. Error Indicators Analysis

To ensure that the numerical results of the ISBM are acceptable, the sensitivity of the numerical results of the ISBM to the total number of boundary nodes N, and to the positions x = ± x r of the fictitious vertical boundaries, are tested for different immersion ratios (h/d) and for different relative lengths (w/d) of the obstacle. To guarantee tests for small and large values of relative water depth kd, the sensitivity tests are presented for two values of the relative water depths, namely kd = 0.2 and kd = 3.
Energy conservation is expressed as | 1 R 2 T r 2 | (see [16] for further details), and can be regarded as an error indicator of the numerical solutions. The results are assumed acceptable when errors do not exceed a value of 10−2.
In Figure 4, we examine effects of the total number of boundary nodes N for a fixed relative length w/d = 2 and different immersion ratios, h/d = 0.25, 0.50, 0.75. In Figure 5 the immersion ratio is fixed at h/d = 0.50 and the relative length is varied, w/d = 2, 4, 6. These results show that, for narrow breakwaters (Figure 4), when the immersion ratio h/d is varied, a minimum number of boundary nodes N = 200 is required to ensure that computational errors are small and do not exceed 10−2. However, for wider breakwaters (Figure 5), the total number of boundary nodes is recommended to be at least N = 600 to guarantee that the computational errors remain below 10−2.
To show the results of the numerical errors with respect to the positions x = ± x r of the fictitious vertical boundaries, we plot the errors against the quantity 2 x r / w (relative width of the computational domain). The total number of boundary nodes used in these tests is N = 600. Figure 6 shows the effects of the quantity 2 x r / w for a fixed relative length w/d = 2 and different immersion ratios, h/d = 0.25, 0.50, 0.75; whereas, in Figure 7, the immersion ratio is fixed at h/d = 0.50 and the relative length is varied, w/d = 2, 4, 6. The results of Figure 6 for narrow breakwaters show that errors reach minimal values when the quantity 2 x r / w is in the range 2.5 < 2 x r / w < 3.5 . For these computations, a representative value of 2 x r / w = 3 is adopted. Thus, the positions of the fictitious vertical boundaries are fixed, respectively, at x = ± 3 w / 2 . For wider breakwaters, as shown in Figure 7, the results indicate clearly that errors are lowest for 1.5 < 2 x r / w < 2.5 . A typical value of 2 x r / w = 2 is assumed for these calculations, and therefore the positions of the fictitious vertical boundaries are fixed, respectively, at x = ± w .

4.2. Reflection and Transmission Coefficients Analysis

This subsection stands deeply on discussing the reflection and transmission coefficients. Three approaches are used in this work to analyze the reflection and transmission coefficients during the wave interactions with submerged rectangular breakwater. The first approach used is the improved version of the ISBM approach that is already validated in reference [21].
The second approach used is the analytical method within the plane wave model (Appendix A), and the third is the analytical method within the corrected plane wave model. Further, the intention is to figure out the possibility to study the reflection and transmission coefficients using the corrected plane wave model without the introduction of evanescent modes. Firstly, the purpose is to study the validity of reflection and transmission coefficients within the corrected plane wave model, the limitations associated with relative length and the immersion ratio of the submerged rectangular breakwater are studied versus the relative water depth kd, which represents essentially the ratio of water depth to wavelength. Secondly, the goal is to discuss the effect of relative length on reflection and transmission coefficients. Therefore, in Figure 8, the effect of relative length is investigated by fixing the immersion ratio at h/d = 0.50, and varying the relative length into w/d = 2, w/d = 4, and w/d = 6. Thirdly, we aim to study the effect of immersion ratios on reflection and transmission coefficients. Hence, in Figure 9, the effect of immersion ratios is studied by fixing the relative length at w/d = 2, and varying the immersion ratio into h/d = 0.25, h/d = 0.50, and h/d = 0.75.
To analyze the reflection and transmission coefficients during the wave interactions with submerged rectangular breakwater, three methods are employed; the first approach is the ISBM, which is already validated in the reference [21]. The second approach used is the analytical method within the plane wave model (Appendix A), and the third is the analytical method within the corrected plane wave model. The plane wave model is corrected by adding 5% to the relative length and subtracting 4% from the immersion ratio. Further, the results of Figure 8 and Figure 9 show that the analytical method within the plane wave model appears slightly far from the ISBM approach, and, as the relative water depth kd increases, the deviation of profiles of the analytical method increases. After several tests, our computational experiments indicate that descriptions of analytical profiles of reflection and transmission coefficients have corrected, and have the same profiles of the ISBM method by adding 5% to the relative length and subtracting 4% from the immersion ratio. Furthermore, the results of the corrected plane wave model are much better than those of the uncorrected plane wave model for all immersion ratios and relative lengths of the obstacle. For these reasons, the corrected plane wave model is strongly recommended. Moreover, the analytical model is based on the use of simple assumptions and simple mathematical manipulations, and the use of the corrected plane wave model presents the advantage of reducing the dimension of algebraic systems treated. In addition, the ISBM method necessitates the mesh generation and the numerical integration. Thus, the analytical model is much easier to use than the ISBM method. Further, the advantages of the corrected plane wave model over the ISBM approach are that it is presented in no condition to the total nodes number, or fictitious vertical boundaries.
On the other hand, the effects of the relative length and immersion ratio on reflection and transmission coefficients, are deeply studied versus the relative water depth kd in Figure 8 and Figure 9. For the sake of clarity, Figure 8, describing the effect of relative length on the variation of reflection and transmission coefficients, shows that the reflection coefficients increase when the relative length increases, whereas the transmission coefficients decrease when the relative length increases. Furthermore, the results show that when the relative length increases, the oscillatory aspect of reflection and transmission coefficients increase, with the appearance of maximum of reflection and transmission coefficients at long wavelengths. Moreover, as the relative water depth increases, the maximum of reflection and transmission coefficients decrease. In addition, Figure 9, describing the effect of immersion ratio on the variation of reflection and transmission coefficients, shows that the reflection coefficients decrease when the immersion ratio increases, whereas the transmission coefficients increase when the immersion ratio decreases.

5. Conclusions and Perspectives

In this research paper, the reflection and transmission coefficients during the interactions of regular wave-rectangular breakwater cited at the bottom are studied using the ISBM approach, the analytical approach within the plane wave model (Appendix A), and the corrected plane wave model. Firstly, to verify the capability of the proposed NWT, the error sensitivity to the total number of boundary nodes, and to positions of fictitious vertical boundaries, for different immersion ratios (h/d) and for different relative length (w/d) of the obstacle, are deeply investigated. Further, the results of this work show that the minimum number of boundary nodes N = 200 is required to ensure an accurate computation solution for narrow breakwaters, and N = 600 for wider breakwaters. Furthermore, to ensure minimal values of errors, the quantity 2 x r / w is recommended to be in the range 2.5 2 x r / w 3.5 for narrow breakwaters, and 1.5 2 x r / w 2.5 for wider breakwaters.
Next, the improved version of the meshless singular boundary method (ISBM), the analytical approach within the plane wave model (Appendix A), and the corrected plane wave model are compared to investigate the capacity of the plane wave and corrected plane wave models to study the reflection and transmission coefficients during interactions of regular wave-rectangular breakwater cited at the bottom of tank. For the sake of details, the plane wave model is efficient for small relative water depths kd; whereas, by correcting the plane wave model by adding 5% to the relative length, and subtracting 4% from the immersion ratio, the corrected plane wave model appears to be in an acceptable agreement with the ISBM approach. Afterwards, the results show that the corrected plane wave model is successful in deeply analyzing the effects of relative length and immersion ratio on reflection and transmission coefficients. Then, it is recommended to use this approach compared to the ISBM method that is conditioned to the nodes number and the position of fictitious vertical boundaries. As perspective, we endeavor to study the wave-current–structure interactions using the Generating-Absorbing Boundary Conditions (GABCs) approach [22] to meticulously study the reflection and transmission coefficients for different aspects of currents.

Author Contributions

Conceptualization, M.L., D.D., C.N., D.N. and K.K.; methodology, M.L., D.D., C.N., D.N. and K.K.; software, M.L., D.D., C.N., D.N. and K.K.; validation M.L., D.D., C.N., D.N. and K.K.; formal analysis M.L., D.D., C.N., D.N. and K.K.; investigation, M.L., D.D., C.N., D.N. and K.K.; resources, M.L., D.D., C.N., D.N. and K.K.; writing—original draft preparation, M.L., D.D., C.N., D.N. and K.K.; writing—review and editing, M.L., D.D., C.N., D.N. and K.K.; visualization, M.L., D.D., C.N., D.N. and K.K.; supervision, M.L., D.D., C.N., D.N. and K.K.; project administration, M.L., D.D., C.N., D.N. and K.K.; funding acquisition, M.L., D.D., C.N., D.N. and K.K. All authors have read and agreed to the published version of the manuscript.

Funding

The work of DD has been supported by the French National Research Agency, through the Investments for Future Program (ref. ANR-18-EURE-0016—Solar Academy).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The Authors would like to express their gratitude to the Referees who helped us to improve our manuscript quality.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a The amplitude of the incident wave
TThe period
RThe reflection coefficient
T r The transmission coefficient
dThe water depth
hThe water depth above the obstacle
wThe length of the obstacle
h/dThe immersion ratio
w/dThe relative length
k The wave number that verifies the dispersion relation ω 2 g = k   th ( k d )
σ The wave number above of the obstacle that verifies the dispersion relation ω 2 g = σ   th ( σ h )

Appendix A

In this section, the objective is to present an analytical approach to calculate the reflection and transmission coefficients during the interaction of regular wave-rectangular obstacle cited at the bottom of tank. Within the plane wave approach, Dingemans [23] has already presented the reflection and transmission coefficients based on matching conditions approach where the continuity of the velocity potential as long with the horizontal velocity over the subdomains is neglected. In this paper, we subdivide the domain of study (cf. Figure A1) into three sub-domains ( D 1 , D 2 , D 3 ), where the continuity of velocity potential and horizontal velocity are expressed at edges of the obstacle to determine the reflection and transmission coefficients. By expressing the connection conditions at edges of the obstacle, our approach appears to be simpler and mre accurate than presented by Dingemans [23].
Figure A1. Domain of study.
Figure A1. Domain of study.
Geosciences 11 00430 g0a1
At each subdomain the velocity potential is expressed as:
-
at the subdomain D 1
ϕ 1 = a [ e x p ( j k ( x + w / 2 ) ) + R e x p ( j k ( x + w / 2 ) ) ] c h ( k ( y + d ) )  
-
at the subdomain D 2
ϕ 2 = a T r e x p ( j k ( x w / 2 ) ) c h ( k ( y + d ) )
-
at the subdomain D 3
ϕ 3 = a [ C e x p ( j σ x ) + D e x p ( j σ x ) ] c h ( σ ( y + h ) )
The continuity of the velocity potential and horizontal velocity are expressed at two positions x = w / 2 and position x =   w / 2 as:
  • At the position x = w / 2
    h 0 ϕ 1 ( w / 2 , y ) c h ( σ ( y + h ) ) d y = h 0 ϕ 3 ( w / 2 , y ) c h ( σ ( y + h ) ) d y
    h 0 ϕ 1 ( w / 2 , y ) x c h ( k ( y + d ) ) d y = h 0 ϕ 3 ( w / 2 , y ) x c h ( k ( y + d ) ) d y
  • At the position x=   w / 2
    h 0 ϕ 1 ( w / 2 , y ) c h ( σ ( y + h ) ) d y = h 0 ϕ 3 ( w / 2 , y ) c h ( σ ( y + h ) ) d y
    h 0 ϕ 1 ( w / 2 , y ) x c h ( k ( y + d ) ) d y = h 0 ϕ 3 ( w / 2 , y ) x c h ( k ( y + d ) ) d y
    after expressing the connection conditions at the positions x =   w / 2 and x = w / 2 , we obtain a linear algebraic system as
    { I 1 ( 1 + R ) = I 2 ( Cz + D z ¯ )                 I 1 T r = I 2 ( C z ¯ + Dz )   k I 3 ( 1 R ) = σ I 1 ( Cz D z ¯ )   k I 3 T r = σ I 1 ( C z ¯ Dz )              
    where:
    I 1 = h 0 c h ( k ( y + d ) ) c h ( σ ( y + d ) ) d y
      I 2 = h 0 c h 2 ( σ ( y + h ) ) d y
    I 3 = h 0 c h 2 ( k ( y + d ) ) d y
    by combining the system of Equation (8), we obtain
    Cz   = 1 2 I 1 I 2 ( 1 + R ) + 1 2 k σ I 3 I 1 ( 1 R )
    D z ¯ = 1 2 I 1 I 2 ( 1 + R ) 1 2 k σ I 3 I 1 ( 1 R )
Then, by injecting the expressions (12) and (13) in the second and fourth equation of the system (8), we obtain
I 1 I 2 T r = [ I 1 I 2 ( 1 + R ) + k I 3 σ I 1 ( 1 R ) ] z 2 ¯ 2 + [ I 1 I 2 ( 1 + R ) k I 3 σ I 1 ( 1 R ) ] z 2 2
k I 3     σ I 1 T r = [ I 1 I 2 ( 1 + R ) + k I 3 σ I 1 ( 1 R ) ] z 2 ¯ 2 [ I 1 I 2 ( 1 + R ) k I 3 σ I 1 ( 1 R ) ] z 2 2
with z = exp ( j σ w / 2 ) , more explicitly,
( I 1 I 2 + k I 3 σ I 1 ) T r = [ I 1 I 2 ( 1 + R ) + k I 3 σ I 1 ( 1 R ) ] z 2 ¯
( I 1 I 2 k I 3 σ I 1 ) T r = [ I 1 I 2 ( 1 + R ) k I 3 σ I 1 ( 1 R ) ] z 2  
Then, we obtain a matrix system as
[ A B B A ] [ T r 0 ] = [ z 2 ¯ 0 0 z 2 ] [ A B B A ] [ 1 R ]
The resolution of the matrix system makes it possible to obtain the expression of the reflection and transmission coefficients as
R = z 2 z 2 ¯ z 2 ¯ B / A z 2 A / B
T r = B / A A / B z 2 ¯ B / A z 2 A / B
where:
A = I 1 I 2 + k I 3 σ I 1   and   B = I 1 I 2 k I 3 σ I 1

Appendix B

The potentials at the fictitious boundaries Γ and Γ + , defined at x = x r and x = + x r respectively, are written as:
ϕ = a L T cosh ( k y ) sinh ( k d ) e i k ( x + x r ) a r L T cosh ( k y ) sinh ( k d ) e i k ( x + x r )
ϕ + = a t L T cosh ( k y ) sinh ( k d ) e i k ( x x r )
where a , a r and a t are the incident, reflected and transmitted wave amplitudes respectively.
The reflection and transmission coefficients R and T r are defined as:
R = a r a   and   T r = a t a
Substituting the relations given by Equation (A24) into Equations (A22) and (A23) we get:
ϕ = a L T cosh ( k y ) sinh ( k d ) e i k ( x + x r ) a R L T cosh ( k y ) sinh ( k d ) e i k ( x + x r )
ϕ + = a T L T cosh ( k y ) sinh ( k d ) e i k ( x x r )
Now defining the following expressions as:
| A 0 | = a R L T   and   | A 0 + | = a T L T
and substituting into Equations (A25) and (A26) , we finally achieve the final expressions of the inflow and outflow potentials:
ϕ = a L T cosh ( k y ) sinh ( k d ) e i k ( x + x r ) + A cosh ( k y ) sinh ( k d ) e i k ( x + x r )
ϕ + = A + cosh ( k y ) sinh ( k d ) e i k ( x x r )
After calculating the potentials φ and φ + from the ISBM, they are matched to those of Equations (A28) and (A29) at x = x r and x = + x r respectively, then the unknown coefficients A and A + are evaluated as follows. When Equations (A28) and (A29) are evaluated at x = x r and x = + x r they become:
ϕ ( x r , y ) = a L T cosh ( k y ) sinh ( k d ) + A cosh ( k y ) sinh ( k d )
ϕ + ( + x r , y ) = A + cosh ( k y ) sinh ( k d )
We can now multiply both sides of Equations (A30) and (A31) by cos h ( k y ) to give:
ϕ ( x r , y ) cos h ( k y ) = a L T [ cos h ( k y ) ] 2 sinh ( k d ) + A [ cos h ( k y ) ] 2 sinh ( k d )
ϕ + ( + x r , y ) cos h ( k y ) = A + [ cos h ( k y ) ] 2 sinh ( k d )
Integrating y in Equations (A32) and (A33) between 0 and d, and after manipulating the resulting algebraic expressions, we get:
A = ( a · L T ) + k N 0 · cosh ( k · d ) 0 d ϕ ( x r , y ) cos h ( k y ) d y
A + = k N 0 cos h ( k d ) 0 d ϕ + ( + x r , y ) cos h ( k y ) d y
where N 0 = 1 2 ( 1 + 2 k d sin h ( 2 k d ) )   .
Referring to Equation (A27) , the reflection and transmission coefficients are Finally calculated as:
R = | A 0 | T a · L   and   T r = | A 0 + | T a · L

References

  1. Dean, W.R. On the reflexion of surface waves by a submerged plane barrier. In Mathematical Proceedings of the Cambridge Philosophical Society; Cambridge University Press: Cambridge, UK, 1945; Volume 41, pp. 231–238. [Google Scholar]
  2. Takano, K. Effets d’un obstacle parallélépipédique sur la propagation de la houle. Houille Blanche 1960, 3, 247–267. [Google Scholar] [CrossRef] [Green Version]
  3. Patarapanich, M. Forces and moment on a horizontal plate due to wave scattering. Coast. Eng. 1984, 83, 279–301. [Google Scholar] [CrossRef]
  4. Liu, P.L.F.; Wu, J. Wave transmission through submerged apertures. J. Waterw. Port Coast. Ocean Eng. 1987, 113, 660–671. [Google Scholar] [CrossRef]
  5. Stamos, D.J.; Hajj, M.R.; Telionis, D.P. Performance of hemi-cylindrical and rectangular submerged breakwaters. Ocean Eng. 2003, 30, 813–828. [Google Scholar] [CrossRef]
  6. Molin, B.; Kimmoun., O.; Liu., Y.; Remy., F.; Bingham, H.B. Experimental and numerical study of the wave run-up along a vertical plate. J. Fluid Mech. 2010, 654, 363–386. [Google Scholar] [CrossRef]
  7. Liu, Y.; Li, H.J.; Zhu, L. Bragg reflection of water waves by multiple submerged semi-circular breakwaters. Appl. Ocean Res. 2016, 11, 67–78. [Google Scholar] [CrossRef]
  8. Scyphers, S.B.; Powers, S.P.; Heck, K.L. Ecological Value of Submerged Breakwaters for Habitat Enhancement on a Residential Scale. Environ. Manag. 2015, 55, 383–391. [Google Scholar] [CrossRef] [PubMed]
  9. Mei, C.; Black, J. Scattering of surface waves by rectangular obstacles in waters of finite depth. J. Fluid Mech. 1969, 38, 499–511. [Google Scholar] [CrossRef]
  10. Massel, S. Harmonic generation by waves propagating over a submerged step. Coast. Eng. 1983, 7, 357–380. [Google Scholar] [CrossRef]
  11. Driscoll, A.; Dalrymple, R.; Grilli, S. Harmonic generation and transmission past a submerged rectangular obstacle. Coast. Eng. 1992, 1142–1152. [Google Scholar] [CrossRef]
  12. Szmidt, K. Finite difference analysis of surface wave scattering by underwater rectangular obstacles. Arch. Hydro-Eng. Environ. Mech. 2010, 57, 179–198. [Google Scholar]
  13. Chen, W.; Gu, Y. An improved formulation of singular boundary method. Adv. Appl. Math. Mech. 2012, 4, 543–558. [Google Scholar] [CrossRef]
  14. Gu, Y.; Chen, W.; Zhang, J.Y. Investigation on near-boundary solutions by singular boundary method. Eng. Anal. Bound. Elem. 2012, 36, 1173–1182. [Google Scholar] [CrossRef]
  15. Gu, Y.; Chen, W. Infinite domain potential problems by a new formulation of singular boundary method. Appl. Math. Model. 2013, 37, 1638–1651. [Google Scholar] [CrossRef]
  16. Chioukh, N.; Çevik, E.; Yüksel, Y. Reflection and transmission of regular waves from/through single and double perforated thin walls. China Ocean Eng. 2017, 31, 466–475. [Google Scholar] [CrossRef]
  17. Bakhti, Y.; Chioukh, N.; Hamoudi, B.; Boukhari, M. A multi-domain boundary element method to analyze the reflection and transmission of oblique waves from double porous thin walls. J. Mar. Sci. Appl. 2017, 16, 276–285. [Google Scholar] [CrossRef]
  18. Yueh, C.Y.; Chuang, S.H. A boundary element model for a partially piston-type porous wave energy converter in gravity waves. Eng. Anal. Bound. Elem. 2012, 36, 658–664. [Google Scholar] [CrossRef]
  19. Brebbia, C.A.; Dominguez, J. Boundary Elements: An Introductory Course. In Computational Mechanics Publications; WIT Press: Southampton, UK; Boston, MA, USA, 1992. [Google Scholar]
  20. Gu, Y.; Chen, W. Recent advances in singular boundary method for ultra-thin structural problems. In Boundary Elements and Other Mesh Reduction Methods XXXVI; WIT Transactions on Modelling and Simulation: Southampton, UK; Boston, MA, USA, 2014; Volume 56, pp. 233–243. [Google Scholar]
  21. Senouci, F.; Chioukh, N.; Dris, M.E.A. Performance Evaluation of Bottom-Standing Submerged Breakwaters in Regular Waves Using the Meshless Singular Boundary Method. J. Ocean Univ. China 2019, 18, 823–833. [Google Scholar] [CrossRef]
  22. Loukili, M.; Dutykh, D.; Kotrasova, K.; Ning, D. Numerical Stability Investigations of the Method of Fundamental Solutions Applied to Wave-Current Interactions Using Generating-Absorbing Boundary Conditions. Symmetry 2021, 13, 1153. [Google Scholar] [CrossRef]
  23. Dingemans, M. Water Wave Propagation Over Uneven Bottoms. Adv. Ser. Ocean Eng. 1997, 13, 1–1016. [Google Scholar] [CrossRef]
Figure 1. Breakwater system of this study.
Figure 1. Breakwater system of this study.
Geosciences 11 00430 g001
Figure 2. Problem definitions for the breakwater system.
Figure 2. Problem definitions for the breakwater system.
Geosciences 11 00430 g002
Figure 3. Domain discretization.
Figure 3. Domain discretization.
Geosciences 11 00430 g003
Figure 4. Error indicator versus the number of boundary nodes for different immersion ratios.
Figure 4. Error indicator versus the number of boundary nodes for different immersion ratios.
Geosciences 11 00430 g004
Figure 5. Error indicator versus the number of boundary nodes for different relative lengths.
Figure 5. Error indicator versus the number of boundary nodes for different relative lengths.
Geosciences 11 00430 g005aGeosciences 11 00430 g005b
Figure 6. Error indicator versus the location of vertical fictitious boundaries for different immersion ratios.
Figure 6. Error indicator versus the location of vertical fictitious boundaries for different immersion ratios.
Geosciences 11 00430 g006
Figure 7. Error indicator versus the location of vertical fictitious boundaries for different relative lengths.
Figure 7. Error indicator versus the location of vertical fictitious boundaries for different relative lengths.
Geosciences 11 00430 g007
Figure 8. Variations of reflection and transmission coefficients by fixing the immersion ratio at h/d = 0.50 and varying the relative length: (a) w/d = 2, (b) w/d = 4, and (c) w/d = 6.
Figure 8. Variations of reflection and transmission coefficients by fixing the immersion ratio at h/d = 0.50 and varying the relative length: (a) w/d = 2, (b) w/d = 4, and (c) w/d = 6.
Geosciences 11 00430 g008
Figure 9. Variations of reflection and transmission coefficients by fixing the relative length at w/d = 2, and varying values of immersion ratio: (a) h/d = 0.25, (b) h/d = 0.50, (c) h/d = 0.75.
Figure 9. Variations of reflection and transmission coefficients by fixing the relative length at w/d = 2, and varying values of immersion ratio: (a) h/d = 0.25, (b) h/d = 0.50, (c) h/d = 0.75.
Geosciences 11 00430 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Loukili, M.; Dutykh, D.; Nadjib, C.; Ning, D.; Kotrasova, K. Analytical and Numerical Investigations Applied to Study the Reflections and Transmissions of a Rectangular Breakwater Placed at the Bottom of a Wave Tank. Geosciences 2021, 11, 430. https://doi.org/10.3390/geosciences11100430

AMA Style

Loukili M, Dutykh D, Nadjib C, Ning D, Kotrasova K. Analytical and Numerical Investigations Applied to Study the Reflections and Transmissions of a Rectangular Breakwater Placed at the Bottom of a Wave Tank. Geosciences. 2021; 11(10):430. https://doi.org/10.3390/geosciences11100430

Chicago/Turabian Style

Loukili, Mohammed, Denys Dutykh, Chioukh Nadjib, Dezhi Ning, and Kamila Kotrasova. 2021. "Analytical and Numerical Investigations Applied to Study the Reflections and Transmissions of a Rectangular Breakwater Placed at the Bottom of a Wave Tank" Geosciences 11, no. 10: 430. https://doi.org/10.3390/geosciences11100430

APA Style

Loukili, M., Dutykh, D., Nadjib, C., Ning, D., & Kotrasova, K. (2021). Analytical and Numerical Investigations Applied to Study the Reflections and Transmissions of a Rectangular Breakwater Placed at the Bottom of a Wave Tank. Geosciences, 11(10), 430. https://doi.org/10.3390/geosciences11100430

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