Oblique Long Wave Scattering by an Array of Bottom-Standing Non-Smooth Breakwaters

: Bragg scattering of surface gravity waves by an array of submerged bottom-standing non-smooth breakwaters is studied under the assumption of linearized long wave theory. The closed-form long-wave analytical solutions are derived and validated by comparing them with the results available in the literature. The role of various physical parameters such as breakwaters friction coefﬁcient, depth, width and gap between the adjacent breakwaters are investigated by analyzing the reﬂection and transmission coefﬁcients. Further, the time-domain simulation for the scattering of long gravity waves over multiple breakwaters is analysed for different values of parameters of breakwaters. The results reveal that the rough surface of the breakwater plays a vital role in reducing wave reﬂection and transmission. Moreover, it is observed that the transmitted wave dissipates completely for larger values of friction parameters. For certain critical angles, change in wave dissipation becomes maximum due to the variation of phase of the incident wave. Various ﬁndings can be considered as benchmark results for the design of the non-smooth structures to attenuate the waves based on the Bragg reﬂection.


Introduction
Shoaling, refraction, diffraction, and reflection have an impact on the wave patterns that are created on the water's surface, causing them to undergo numerous changes. Bragg resonance is one such mechanism that is essential for spotting natural occurrences like earthquakes or free surface oscillations. This phenomenon occurs as a rise and fall of the water surface and often leads to an increase in wave amplitude over time. The mechanism of Bragg resonance was originally investigated by Bragg and Bragg [1] in crystallography. A few decades later, multiple groundbreaking research experimentally, theoretically, and computationally validated the Bragg resonance between surface water waves and undulating bottom topography (e.g., Heathershaw [2], Mei et al. [3], Cho et al. [4]). Bragg resonance could cause destruction when its period matches with the period of another external source, such as the wind, an earthquake, or rapid changes in air pressure. By adding friction to the top surface of the breakwater, one might lessen the resonance amplitude during Bragg resonance. Thus, there is a need to examine the effects of surface friction of a series of breakwaters on Bragg reflection.
There has been significant advancement in the understanding of Bragg scattering of long waves over submerged structures by a number of studies. Additionally, several analytical solutions have been developed based on linear long-wave equation (LWE) or modified mild-slope equation (MMSE). Liu et al. [5] analysed the reflection of waves by four different kinds of Bragg breakwaters such as triangular, rectified sinusoidal, trapezoidal, and parabolic and produced the associated optimal collocation curves for maximum wave reflection. Liu et al. [6] obtained the analytical MMSE solutions for Bragg reflection by a series of sinusoidal ripples and trapezoidal bars. Kar et al. [7] studied the Bragg scattering of long waves by a series of trenches using the solutions of Bessel equations and the findings were corroborated using numerical data produced based on the boundary element approach (BEM). Kar et al. [8] examined the Bragg scattering of long waves by an array of flexible floating plates in the presence of several submerged trenches. Gao et al. [9] analysed the Bragg resonance produced by deploying a sequence of sinusoidal bars with varying amplitudes and numbers outside a harbour. The interplay of waves, currents, and surges in a region of the shallow continental shelf was explored by Sahoo et al. [10]. Barman et al. [11] investigated the Bragg scattering of long flexural gravity waves by an array of submerged trenches in the perspective of blocking dynamics. Kar et al. [12] discovered the Bragg scattering of surface gravity waves due to surface wave interaction with multiple bottom undulations and a semi-infinite floating flexible structure. Boral et al. [13] investigated the role of viscous damping in the flexural gravity wave interaction in shallow water waves. Guo et al. [14] explored the phase downshift or upshift of Bragg resonance for water wave reflection by an array of cycloidal bars or trenches using the mild slope equation approach. Liu et al. [15] developed MMSE-based analytical modelling to study wave reflection by a single cycloidal geo tube or trench. Xie et al. [16] studied the long wave reflection by an array of submerged trapezoidal breakwaters on a sloping seabed.
The Bragg resonant reflections over several arrays of smooth, impermeable breakwaters were studied by a number of researchers (see Kar et al. [17] and Vijay et al. [18]). However, in actual use, breakwaters are often porous and not smooth. These structures withstand wave resonance and absorb and dissipate wave energy due to its non-smooth/friction nature. Losada et al. [19] explored monochromatic oblique wave propagation across and through porous beds or on a submerged rectangular structure. Few researchers have investigated water flows and resonant waves across the basin that entail the presence of wind and friction (see Visser et al. [20]). Additionally, Mullarney et al. [21] employed analytical technique to analyze wave propagation in a basin in the presence of bottom friction, while Reef et al. [22] investigated using numerical methodology. Ni et al. [23,24] examined the Bragg resonant reflection of water waves by Bragg breakwaters of two forms of porous bars, namely rectangular and trapezoidal, on a sloping permeable seabed. Magdalena et al. [25] analysed the attenuation of waves by mangroves which are modelled as porous structures. Magdalena et al. [26] performed both analytical and numerical studies for seiches in a closed basin with bottom friction. Magdalena et al. [27] investigated the resonance phenomena in lakes using a mathematical model. Subsequently, Magdalena et al. [28] examined the resonant periods of earthquakes in semi-closed basins with complicated bottom topography.
Furthermore, the effect of the rough surface of a rectangular submerged breakwater on resonance phenomena was studied by Magdalena et al. [29] using staggered finite volume method. They determined the optimal friction coefficient and shows that friction can prevent the resonance, even reduce the wave amplitude. Using a mathematical model based on the shallow water equations model, Magdalena et al. [30] simulated the resonance phenomena in basins of various shapes. They demonstrated that if either one of the mode parameters of rectangular basin is zero, a lower value of friction will be needed to prevent the wave resonance otherwise the value of the rougher bottom must be increased to prevent the incoming wave from resonating. Magdalena et al. [31] reported the effect of bottom friction on harbour oscillation for the three different harbour geometry (rectangular, triangular, and semi-parabolic). The results show that a relatively small value of friction is sufficient to prevent resonance in rectangular harbours whereas for the triangular harbour, a larger value of friction requires to prevent resonance. As rough surface plays a vital role on reduction of incidence wave energy, Bragg resonance occurring due to the multiple structures can be reduced by considering rough surface of the breakwaters. This has motivated the authors to study the role of rough surface of multiple bottom-standing structures on Bragg scattering using transform matrix method (TMM).
The rest of the paper is organised as follows: The detailed formulation of the physical model is given in Section 2. The details of the used Transfer Matrix Method and solution of the mathematical model are described in Section 3. The impact of roughness of the breakwaters on Bragg resonance and the role of various physical parameters such as breakwater width and depth, and the gap between the adjacent sides of breakwaters on Bragg scattering are analysed in Section 4. The shape of wave packets in time-domain interacting with rectangular smooth/non-smooth breakwaters is demonstrated in Subsection 4.1 for various values of the breakwater width, depth, breakwater friction, and initial width of the wave packet. Finally, the important observations are highlighted in Section 5.

Mathematical Formulation
The scattering of surface gravity waves by an array of non-smooth breakwaters is investigated analytically using the technique of Transfer Matrix Method (TMM) under the assumption of long wave theory. The problem is studied in the two-dimensional Cartesian coordinate system with x − y plane being assumed as the horizontal plane and z-axis being directed vertically upward. The position of the series of breakwaters associated with water depths is shown in Figure 1. It is assumed that the motion is simple-harmonic in time with angular frequency ω (see Dingemans [32], Behera and Sahoo [33]). Thus, the free surface elevation ζ(x, y, t) is assumed to be of the form ζ(x, y, t) =Re{η(x)e −ik y y−iωt } where, k y = k 1 sin θ, η(x) is the spatial component of the free surface elevation with k 1 being the incident wave number. Under the assumption that u, and v are the velocity component of the fluid in two dimensions along x, and y− directions respectively, the equation of motion can be expressed (see as Sahoo [34] for details) as with g being the acceleration due to gravity. Moreover, the continuity equation associated with the linearized long wave equation is given by Thus, from Equations (1) and (2), the two-dimensional long wave equation over the flat bottom is expressed as To account for the roughness of the breakwater, the two-dimensional linearized long waves equations are modified by adding friction component C f ωu and C f ωv that represents bottom friction where C f being the friction coefficient. This correction is included into the equation of motion Equation (1). As a result, the enhanced model may be expressed (see [26,29,30]) as ∂v ∂t + g ∂ζ ∂y After eliminating u, v from Equations (4)-(6), the two-dimensional long wave equation in terms of ζ can be expressed as Thus, the two-dimensional long wave equation over the breakwater region can be expressed in terms of η as where k is the solution of the dispersion equation ω 2 + iω 2 C f = gk 2 h 2 .

Closed-Form Analytical Solution
The small amplitude long gravity waves propagating over the series of submerged rectangular breakwaters as shown in Figure 1 are studied analytically using the technique of matrix multiplication. The form of solutions of Equations (3) and (8) in smooth and non-smooth regions having constant depths h 1 and h 2 respectively are expressed as where and A R and A T the unknown complex amplitudes of the reflected and the transmitted waves respectively. The unknowns constants A j 1 , A j 2 and B j 1 , B j 2 which are related to the amplitudes of the waves propagating along rough and smooth sea bed region respectively. Moreover, the reflection K r and transmission coefficients K t are defined as

Transfer Matrix Method (TMM)
Various unknown constants in Equation (9) are obtained using the boundary conditions on the breakwater boundaries along with the interface boundaries following the transfer matrix method (see Zeng et al. [35]) which is briefly discussed in the subsequent discussion.
Assuming x =x as an interface boundary, continuity of pressure and flow fluxes across the interface boundaries of the breakwater and open water regions yield Using the matching conditions at x = x 1,1 as in Equations (9) and (11) yields where Similarly, using the matching conditions at x = x j,2 for j = 1, 2, 3, ... . . , N − 1, where Finally, using the matching conditions as in Equation (11) Thus, the matrix multiplication to Equations (12)-(15) gives where whilst N represents the number of breakwaters. Therefore, Equation (16) gives rise to the reflection and transmission coefficients as

Results and Discussions
To analyze the effects of different waves and structural parameters on wave scattering over a non-smooth array of breakwaters for numerical computation, numerical codes are written using MATLAB software. It is pertinent to mention that certain physically realistic values of the wave and structural parameters, which are listed in Table 1, are used in the computation for an efficient breakwater system. Furthermore, in the subsequent numerical results and discussion, the legend and caption of the figures highlights the values of different physical parameters which are different from the parametric values mentioned in Table 1. The variation of the reflection coefficient K r against non-dimensional wave number k 1 h 1 for double smooth surface bottom-standing rigid breakwaters and trenches Figure 2a,b, respectively. The primary objective of these two figures is to validate the present result with published results in two different special cases. Thus, in Figure 2a,b, two different sets of parameters are taken from Zeng et al. [35] and Kar et al. [7], respectively for the purpose of validation. Figure 2a reveals that the results obtained by present theory is well matched with the results of Zeng et al. [35] for long wave scattering by multiple breakwaters with N = 2, C f = 0, h 1 = 4 m, h 2 = 3 m, w = 90 m, d = 180 m. In addition, in Figure 2b, the variation of reflection coefficient K r against wave number k 1 h 1 is plotted in the case of normal incident waves with double smooth trenches, which agree well with Figure 3a in Kar et al. [7] for In both the figures, Bragg resonance occurs due to the presence of a pair of smooth structures where C f = 0. In the subsequent figures, the effect of friction parameter C f on Bragg resonance through the reflection, transmission, and dissipation coefficients is analysed for various wave and structural parameters. In Figure 3, the variation of (a) reflection coefficient Figure 3a shows that Bragg reflection occurs for wave scattering by an array of multiple smooth breakwaters (C f = 0), which is similar to that of Kar et al. [7,12]. In Figure 3a,b, the occurrence of zero reflection and full transmission in the absence of friction happens due to constructive and destructive interference of the trapped waves. The peaks in each bandwidth of Bragg reflection coefficient K r remain periodic with an increase in wave number k 1 h 1 . In the presence of friction C f , the periodic pattern in the reflection coefficient K r is not observed which happens due to the dissipation of the incident wave energy by the non-smooth breakwater. These similar phenomena were observed by Gayathri et al. [36]. Further, it is observed that for certain range of wavenumber i.e., k 1 h 1 < 2.5, wave reflection K r decreases as friction parameter C f increases. Figure 3b reveals that no wave is transmitted in the presence of friction C f for k 1 h 1 > 2.5. Moreover, for wave number k 1 h 1 < 2.5, transmission coefficient K t decreases as the value of friction parameter C f increases which is similar as observed in Figure 3a. Figure 3c reveals that for wave number k 1 h 1 > 2.5, there is nearly 90% wave energy dissipation happens for irrespective of the value of friction parameter C f . Figure 3a shows that in the presence of friction, the oscillatory pattern diminishes, which may be due to the dissipation of wave energy. In Figure 3c, variation of energy dissipation in the absence of friction (with C f = 0) is zero, which is be due to the distribution of wave energy into the reflected and transmitted waves over the smooth surface of the breakwater.  In Figure 4, the variation of (a) reflection coefficient K r , (b) transmission coefficient K t , and (c) dissipation coefficient K d versus wavenumber k 1 h 1 are plotted for different values of N with C f = 1.5, w/h 1 = 1, d/h 1 = 2, θ = 20 • , h 2 /h 1 = 0.25, h 1 = 5 m. In Figure 4a, the reflection K r decreases as the number of breakwaters N increases, which is similar as observed in Behera and Sahoo [33]. Furthermore, the reflection coefficient K r remains constant beyond k 1 h 1 > 2.5 irrespective of number of breakwaters N which is not the case for smaller wave number. This is due to the fact that the effect of bottom undulation in shallow water region is more as compared to that of the deep water region. Figure 4b depicts that transmission coefficient K t decays faster as the number of breakwaters N increases. Moreover, in the presence of friction, no Bragg resonance occurs for wave transmission whereas Bragg resonance in wave reflection is observed for wave number k 1 h 1 < 2. A comparison of Figure 4c with Figure 3c reveals that in both the figures, nearly 90% incident wave energy dissipation occurs for wavenumber k 1 h 1 > 2.5 irrespective of values of C f or number of breakwaters N. In Figure 4c, the energy dissipation is due to the increase in the number of non-smooth breakwaters.
In Figure 5, the variation of (a) reflection coefficient K r and (b) transmission coefficient Figure 5a depicts that reflection coefficient K r decreases with an increase in breakwater depth h 2 /h 1 . This is due to the decrease in the height of the non-smooth breakwater which will transmit more wave energy.  Moreover, in the presence of friction C f , the reflection coefficient K r follows an oscillatory pattern with an increase of wave number k 1 h 1 . The peaks in the wave reflection K r shift forward with an increase in wavenumber k 1 h 1 . The shifting of peaks is due to the phase change in the reflected and transmitted waves. In Figure 5a, zero reflection does not occur as the height of breakwaters increases, which is similar to the observation as studied in [33]. Furthermore, it is observed that the number of peaks remains the same irrespective of the values of breakwaters depth. Figure 5b depicts that there is negligible variation of transmission coefficient K t observe for smaller values of wavenumber i.e., k 1 h 1 < 1.5. Moreover, the transmission coefficient K t increases with an increase in breakwater depth h 2 /h 1 . The zero transmission occurs for higher values of wavenumber k 1 h 1 . Both figures reveal that breakwater height plays a major role to reflect the wave energy toward the seaside. Figure 6 describes the variation of (a) reflection coefficient K r , (b) transmission coefficient K t , and (c) dissipation coefficient K d against wavenumber k 1 h 1 for different breakwater width w/h 1 with C f = 0.5, h 2 /h 1 = 0.25, d/h 1 = 2, θ = 20 • , N = 4, h 1 = 5. It is found that reflection coefficient K r follows an oscillatory pattern with an increase in wavenumber k 1 h 1 . Moreover, the oscillatory pattern of reflection coefficient K r demises with an increase in width of the breakwater w/h 1 . Further, it is observed that after a certain width, wave reflection remains constant. The number of resonating peaks for smaller breakwater width w/h 1 is more compared to the higher values of breakwater width w/h 1 . Moreover, the oscillatory pattern in wave transmission K t is not observed. This is due to the dissipation of incident wave energy by the rough surface of the multiple breakwaters. A comparison between the Figure 6a,b reveal that an oscillatory pattern in wave reflection K r occurs whereas no oscillatory peaks are observed in wave transmission. In general, dissipation of incident wave energy increases with an increase in width of the breakwaters as shown in Figure 6c. Additionally, dissipation coefficient K d varies significantly for smaller wave numbers with k 1 h 1 < 2.5 which ensures that the present physical model is more suitable for the long waves as compared to that of short waves. A comparison of Figures 5 and 6 reveals that in the presence of friction C f , the changes in transmission coefficient K t are observed more as compared to that of the reflected coefficient K r . In Figure 7, the variation of (a) reflection coefficient K r , (b) transmission coefficient K t , and (c) dissipation coefficient K d against wavenumber k 1 h 1 is plotted for different gap between the adjacent breakwaters d/h 1 with C f = 0.5, h 2 /h 1 = 0.25, d/h 1 = 2, θ = 20 • , N = 4, h 1 = 5 m. Figure 7a reveals that the peaks in the reflection coefficient K r decreases as gap between the adjacent sides of breakwaters d/h 1 increases which is same as observed in Figure 6a. Moreover, the mutual interaction between incident, reflected and transmitted waves decrease for higher gap gap between the structures d/h 1 thus, the oscillation in wave reflection K r demises as gap d/h 1 increases. Furthermore, due to this reason, the number of resonance peaks in the reflection coefficient K r decreases as the distance d/h 1 widens. Figure 7b shows the changes in transmission coefficient K t are very little as the gap between the breakwaters d/h 1 increases. Moreover, a comparison between the Figure 7a,b reveals that oscillatory trend is observed in wave reflection which is not the case with wave transmission. Figure 7c shows that the peaks in the dissipation coefficient K d decrease as the gap between the adjacent side of breakwaters d/h 1 increases. A comparison of Figure 7c with Figures 4c and 6c and reveals that 90% of the incident wave energy dissipates irrespective of gap between the adjacent breakwaters d/h 1 or number of breakwaters N.
In Figure 8, the variation of (a) reflection coefficient K r , (b) transmission coefficient K t , and (c) dissipation coefficient K d versus angle θ is plotted for different values of C f with d/h 1 = 2, width w/h 1 = 1.0, h 2 /h 1 = 0.25, and N = 4. Figure 8a shows that the reflection coefficient K r increases rapidly after a certain angle i.e., nearly θ = 65 • as the value of the friction parameter increases. Moreover, in the absence of friction (C f = 0), zero reflection occurs for certain countably many angles θ which is similar as observed in [33]. The occurrence of zero reflection is due to the destructive interference of reflected and incident waves. However, zero reflection does not happen in the case of breakwaters with rough surfaces (C f = 0) because of the dissipation caused by the structures. Further, a comparison with Figure 4a exhibits that the oscillatory pattern diminishes as the friction coefficient C f increases. Figure 8b evinces that the transmission coefficient K t decreases with an increase in angle θ which is obvious as reflection coefficient K r increases with an increase in angle θ. Moreover, in the absence of friction, full transmission occurs for certain values of angle θ • , which is due to the constructive interference of the resonating waves. Figure 8c shows that for higher values of friction parameter C f , more energy dissipation occurs and in the absence of friction, no dissipation occurs. Therefore friction parameter C f plays a major role in the dissipation of wave energy.   Figure 9a shows that reflection coefficient K r increases with an increase in angle θ • . Moreover, Figure 9a shows that for greater values of a certain angle, the reflection coefficient increases with an increase in breakwater depth h 2 /h 1 whereas for smaller values of a certain angle, reflection follows the reverse pattern. Figure 9b reveals the transmission coefficient K t decreases with an increases in angle θ. Moreover, it is observed that transmission coefficient K t decreases as breakwater depth h 2 /h 1 increases. Figure 9c shows that dissipation coefficient K d increases as breakwater depth h 2 /h 1 decreases. Further, Figure 9c reveals that in the presence of friction C f , for a certain angle i.e., nearly at 65 • , the maximum wave dissipation occurs. In Figure 10, the variation of (a) reflection coefficient K r and (b) transmission coefficient K t versus angle θ is plotted for different values of N with breakwater width w/h 1 = 1.0, h 2 /h 1 = 0.25, d/h 1 = 2,C f = 0.5, and N = 4. Figure 10a reveals that the reflection coefficient K r decreases with an increase in the number of breakwaters N. Moreover, the reflection coefficient K r increases with an increase in the number of breakwaters N as observed in [26]. Further, it is observed that common minima occur in the case of an even number of breakwaters, whereas common maxima occur in the case of an odd number of breakwaters as the angle of propagation θ increases. Moreover, the bandwidth of each harmonic peak increases with an increase in the number of breakwaters N. Further, Figure 10b reveals that the position of the band of each cycle remains fixed with an increase of oblique angle θ irrespective of the number of breakwaters N. Moreover, for a certain angle, θ, the wave reflection K r becomes completely zero which may be due to the phase mismatch of the resonating waves. Figure 11 describes the changes of (a) reflection coefficient K r , and (b) transmission coefficient K t against gap between the breakwater d/h 1 for different values of C f with breakwater width w/h 1 = 1.0, h 2 /h 1 = 0.5, N = 3. Figure 11a reveals that reflection coefficient K r follows a periodically oscillatory pattern with an increase in the gap between the adjacent breakwaters d/h 1 which is similar as observed in [17]. Further, the amplitude of the peaks in each band decreases with an increase in friction coefficient C f . Further, the sub-harmonic peaks demise with an increase in friction coefficient C f . Further, Figure 11b reveals that in the absence of friction, full transmission occurs for finitely many d/h 1 . The position of the band of each cycle remains fixed with an increase of oblique angle θ irrespective of the number of breakwaters N. Moreover, for certain angle θ, the wave reflection K r becomes completely zero which may be due to the phase mismatch of the resonating waves. In Figure 11a, wave reflection decreases due to the damping of wave energy.   Figure 11. Variation of (a) reflection coefficient K r , and (b) transmission coefficient K t versus gap d/h 1 for different values of C f with breakwater width w/h 1 = 1.0, h 2 /h 1 = 0.5 and N = 3. Figure 12 describes the variation of (a) reflection coefficient K r , and (b) transmission coefficient K t versus breakwater width w/h 1 for different values of C f with breakwater width w/h 1 = 1.0, h 2 /h 1 = 0.5 and N = 3. Figure 12a reveals that in the absence of friction, Bragg reflection occurs and follows periodically oscillatory pattern with an increase in breakwater width w/h 1 which is observed in our previous study (see Kar et al. [7,12]). Moreover, the reflection coefficient K r decreases with an increase in the friction coefficient C f . Further, the reflection coefficient K r remains constant after certain values of breakwater width (w/h 1 > 1) irrespective of friction coefficient C f . Further, it is observed that in the absence of breakwater friction, zero reflection occurs for the countably many values of breakwater width w/h 1 . Moreover, the bandwidth of each harmonic peak remains fixed with an increase in the breakwater width w/h 1 . Further, Figure 12b reveals that in the absence of friction, the amplitude of peaks in each band of the cycle remains the same with an increase in width of the breakwater w/h 1 . Moreover, in the absence of friction, the transmission coefficient K t follows the periodically oscillatory. Further, full transmission occurs for the countably many values of breakwater width w/h 1 . Moreover, in the presence of friction, zero transmission occurs for greater values of certain breakwater width w/h 1 > 1. It is observed that with an increase in breakwater friction C f , wave transmission coefficient K t decreases gradually. The periodic pattern in wave transmission vanishes in the presence of friction C f which is the same as observed in Figure 6b which may be due to phase mismatch of the resonating waves. The results in Figure 12 imply that in the presence of friction, the reduction of wave amplitude happens which may be due to the occurrence of wave energy dissipation.

Time-Dependent Simulations and Numerical Results
The time-dependent results are present for wave packets, which are incident from infinity. The free surface elevation, defined in Section 2 in the time domain, is computed using the formula withf (ω) being the Fourier transform of the incident wave pulse. In Equation (18), η(x, ω) is the frequency-dependent surface displacement of the open water region. The incident wave pulse is assumed to be a Gaussian wave packet of the formf (ω) = b π e −b(ω−ω 0 ) 2 where b is the spreading function and ω 0 is the central frequency of the incident wave pulse.
In Figures 13-15, the free surface elevation η(x, t) is plotted against spatial variable x for different values of time t, breakwater depth h 2 /h 1 and incident angle θ in the (a) absence and (b) presence of the breakwater friction C f with number of breakwaters N = 3. Figures 13a-15a depict that the amplitudes of wave pulses remain the same along the breakwater as observed in [12], whereas Figures 13b-15b show that the amplitude of wave pulse reduces in the presence of friction and demises with an increase in time. The amplitude of wave pulse reduces with an increase in time which is due to the dissipation of energy over the rough surface of breakwaters. The occurrence of this phenomena by time domain simulation is well verified with the results observed in Figure 3b. The reduction of the wave pulse happens due to the presence of bottom friction which plays an important role in restricting the unlimited growth of the wave height.

Conclusions
Under the assumption of linearized long wave theory, Bragg reflection of surface gravity waves by an array of submerged bottom-standing non-smooth breakwaters is investigated in this study. The linearized shallow water equations are solved analytically using the transfer matrix method to calculate the reflection, transmission and dissipation coefficients. The study reveals that wave transmission decays linearly in the presence of friction, whereas wave reflection follows an oscillating pattern similar to that observed for a pair of breakwaters in the literature as in Magdalena et al. [25]. In the presence of friction, the common minima or maxima of the reflection curve shift upward as the wave number increases. In addition, zero reflection or full transmission does not occur in the presence of friction parameters. Moreover, the transmission coefficient decays linearly in the presence of friction in the case of multiple breakwaters, similar to what has been observed in the literature. In contrast, the amplitude of surface elevation decreases in the presence of friction, which is due to wave dissipation. The present study will provide a helpful basis in the design of an array of submerged breakwaters for coastal protection and can be used as bench mark solution in case of long wave propagation. However, these type of submerged structures would not be very effective in the case of short-period small amplitude waves where wave energy concentration is higher near the free surface than that of long waves as discussed in the present study.