Step Approximation for Water Wave Scattering by Multiple Thin Barriers over Undulated Bottoms

: This paper investigates the scattering of oblique water waves by multiple thin barriers over undulation bottoms using the eigenfunction matching method (EMM). In the solution procedures of the EMM, the bottom topographies are sliced into shelves separated by steps. On each step, surface-piercing or/and bottom-standing barriers can be presented or not. For each shelf, the solution is composed of eigenfunctions with unknown coefﬁcients representing the wave amplitudes. Then applying the conservations of mass and momentum, a system of linear equations is resulted and can be solved by a sparse-matrix solver. If no barriers are presented on the steps, the proposed EMM formulation degenerates to the water wave scattering over undulating bottoms. The effects on the barrier lengths, barrier positions and oblique wave incidences by different undulated bottoms are studied. In addition, the EMM is also applied to solve the Bragg reﬂections of normal and oblique water waves by periodic barrier sinusoidal bottoms. The accuracy of the solution is demonstrated by comparing it with the results in the literature.


Introduction
The interaction of surface water waves with an undulation bottom topography is a fundamental interest in understanding wave diffraction along the coastal region. When water waves are generated in deep to a shallow area in natural environments, they experience bathymetric variations, underwater obstacles, nonlinear wave interactions and others. These are some of the main natural causes that lead to coastal hazards, which seriously affect coastline and human activities. Several previous studies [1][2][3][4] have focused on the water wave scattering by bathymetric variations and underwater obstacles.
Breakwaters are conventional structures that often rest on the bottom foundation or pierced on the free surface. They are specially designed to reduce the directed affectation of water waves from the open sea approach to the nearshore area. In addition, they are usually built to satisfy the requirements of maritime economic activities. As traditional breakwaters may require higher construction costs due to larger amount of material and complex installation, thin barriers could be alternative solutions. Thin barriers may create natural habitats and/or fishing farms for ecology conservation, such as oil exploration stations, offshore accommodations and oceanographic research centers. Engineers also use assembled floating thin barriers to constructing artificial islands on the sea to provide hospitable environments on the surface of the water. Thin barriers are also intended to protect the shoreline from the direct impacts of coastal erosion and storm surge. In the literature, several authors have studied the problems of water wave interaction by thin barriers. The issue of water wave scattering by a single surface-piercing or bottomstanding barrier over uniform bottom topography has been examined by Losada et al. [5][6][7] and Abul-Azm [8]. Porter and Evans [9] used the Chebychev polynomials to formulate a complementary approximation and obtained highly accurate results to improve the convergence of the solutions. Additionally, Das et al. [10] and Roy et al. [11] studied the water wave scattering by multiple barriers. Furthermore, Wang et al. [12] considered the combined effects of multiple barriers and a single step bottom.
Numerical solutions are important for solving water wave problems as analytic solutions are rare. Therefore, Berkhoff [13] derived the mild-slope equation (MSE) to solve water wave scattering problems. The MSE has also been employed by Kirby [4], Belibassakis et al. [14][15][16] and Toledo et al. [17], to examine the problems of wave scattering by rippled bottom, nonlinear waves and wave-current interactions, respectively. Even though MSE has been modified and improved several times, it still has some special requirements on the spatial derivatives of the eigenfunctions and needs to be improved. In addition, the MSE has seldom been applied for problems of water wave scattering by thin barriers.
Alternately, Takano [18] developed the eigenfunction matching method (EMM) for solving water wave scattering by an elevated sill and a fixed surface obstacle. Mei and Black [19] applied the EMM for solving surface waves normally incident on a rectangular obstacle in a channel of finite depth. Sequentially, Kirby and his co-authors [20,21] used the EMM to study the problems of oblique wave scattering by trenches. Newman [22,23] employed the EMM to transform the integral equation into an infinite set of algebraic equations in terms of the eigenfunctions describing the evanescent modes. The EMM has also been extended to water wave scattering by rectangular submerged bar [24] and surfacepiercing breakwaters [25]. For waves propagating over an arbitrary bottom topography, Devillard et al. [26], O'Hare and Davies [27] and Tsai et al. [28][29][30] decomposed the bottom profiles into a sequence of flat shelves separated by multiples steps and used the EMM to solved the resultant problem. The accuracy of the EMM solutions was shown by Tsai et al. [31] to be comparable with that of the MSE solutions.
In the present study, we consider the applications of EMM for solving problems of water wave scattering by multiple barriers over undulation bottom profiles. These barriers can be either surface-piercing or bottom-standing. By applying the conservations of mass and momentum to the eigen solutions on shelves, the problems can be converted into a system of linear equations with unknown coefficients representing the wave amplitudes on the shelves. The sparse-matrix solver SuperLU [32] is adopted to solve the resulted system. The present EMM formulation is the extension of the existing EMM for solving problems of water wave scattering solely by undulation bottom profiles [20,27,31] or by barriers [33]. The proposed EMM model is validated by comparisons with analytical solutions in the literatures of Losada et al. [5], Porter and Evans [9], Wang et al. [12], Porter and Porter [34], and Davies and Heathershaw [1] for some degenerated cases. This paper is organized as follows: the wave problem is mathematically modeled, and the EMM solution is developed in Section 2, and the EMM model is validated in Section 3. Discussions on Bragg reflections by multiple periodic thin barriers over sinusoidal rippled bottom are provided in Section 4. Finally, Section 5 presented the conclusions of the whole study.

The Mathematical Model
In this section, surface water wave scattering problems by thin impermeable barriers of submerged and bottom-standing configurations are formulated for waves over undulation bottom topography. Considering a train of monochromatic surface waves with incidence angle γ, amplitude a, angular frequency σ and wavelength λ propagate towards thin impermeable barriers over an undulating bottom. Figure 1 illustrates a schematic representation of the wave scattering problem induced by multiple thin barriers over the undulation bottom. In this figure, a Cartesian coordinate system is taken in which z-axis is chosen vertically upwards, and the x-plane is taken as the rest position of the free surface. The wave amplitude is assumed to be small enough that the linear wave theory is applicable. The wave motion is supposed to be time-harmonic by e −iσt , where σ is defined by 2π/T with T being the wave period, t is the time and i is the unit of complex numbers. To facilitate the approximate bottom topography process, thin barriers and sea bottom could be discretized into a series of M shelves in the interval of x m−1 ≤ x ≤ x m for m = 1, 2, . . . , M with a water depth h m . Furthermore, x 0 = −∞ and x M = ∞ are assumed. The i-th surface-piercing barrier is placed at a location x = v i with the submerged length being equal to a i . In addition, the i-th bottom-standing barrier is placed at a location x = w i with length b i . In the following of this article, a = a i , v = v i , b = b i or w = w i if there is only one surface-piercing or bottom-standing barrier. lation bottom topography. Considering a train of monochromatic surface waves with incidence angle , amplitude a , angular frequency  and wavelength  propagate towards thin impermeable barriers over an undulating bottom. Figure 1 illustrates a schematic representation of the wave scattering problem induced by multiple thin barriers over the undulation bottom. In this figure, a Cartesian coordinate system is taken in which z -axis is chosen vertically upwards, and the x -plane is taken as the rest position of the free surface. The wave amplitude is assumed to be small enough that the linear wave theory is applicable. The wave motion is supposed to be time-harmonic by if there is only one surface-piercing or bottom-standing barrier.
where ( / , / , / ) x y z         denotes the three-dimensional del operator with respect to the three-dimensional Cartesian co-ordinate   , , x y z and m  is the velocity potential.
As maintained by the linear wave theory, the velocity potential is governed by the Laplace equation: Considering the solution on the m-th shelf in the interval x m−1 ≤ x ≤ x m for m = 1, 2, . . . , M, the velocity of the fluid is defined by: where ∇ = (∂/∂x, ∂/∂y, ∂/∂z) denotes the three-dimensional del operator with respect to the three-dimensional Cartesian co-ordinate (x, y, z) and φ m is the velocity potential. As maintained by the linear wave theory, the velocity potential is governed by the Laplace equation: which subjected to the kinematic and dynamic free-surface boundary conditions, respectively as: and −iσφ m + gη m = 0 on z = 0, where η m is surface elevation. Equations (3) and (4) can be combined to obtain: Furthermore, the condition for bottom boundary can be written as: The diffraction of water waves at a step-type bottom is considered by Miles [35]. As described in Figure 2, between two shelves, the velocity potentials φ m and φ m+1 require connection conditions: and ∂φ m ∂x where L g and L b are the vertical intervals of the gap and barrier, respectively. Additionally, the condition for the vertical wall is described by:    In Equation (10), φ stands for either φ m or φ m+1 depending on the waterside of the barrier and L w is the vertical intervals of the wall.
In order to make the solution unique, the following far-field conditions are required: η = a e ik 1.0 x + K R e iθ R e −ik 1.0 x e ik y y as x → −∞ (11) and η = aK T e iθ T e ik M.0 x e ik y y as x → ∞, where K R , K T , θ R and θ T are real numbers, the reflection and transmission coefficients are defined by K R e iθ R and K T e iθ T , respectively. In Equation (11) and Equation (12),k M,0 , k y andk 1.0 are positive real wavenumbers expressed by:k and where k 1,0 = 2π/λ > 0 and k M,0 > 0 are the progressive wavenumbers obtained from the dispersion relation: In Equation (13), the indices vary as m = 1, 2, . . . , M and n = 0, 1, . . .. Additionally, the evanescent wavenumbers k m,n in Equation (13) with n = 1, 2, 3, . . . are defined by: where κ m,n is the n-th smallest positive root of the dispersion relation:

Eigenfuction Matching Method (EMM)
Constructed on the linear theory, the complete solution of the velocity potential for the m-th shelf can be defined by: m,n (x) + B m,n ξ (2) m,n (x) ζ m,n (z)e ik y y (18) for m = 1, 2, 3, . . . , M, where A m,n and B m,n are unknown coefficients to be determined. In Equation (18), N is the number of evanescent modes. To build up a complete solution by the method of the separation of variables, the eigenfunctions from Equation (18) can be defined by: and ξ with It can be noticed that the solutions defined by Equations (13)-(22) satisfy Equations (2), (5) and (6) analytically and the unknown coefficients A m,n and B m,n can be found by using Equations (7)- (12).
The conservation of mass in Equations (8)- (10) can be formulated in the EMM as: where the inner product of two depth eigenfunctions is defined by: where G 1 and G 2 are the depth eigenfunction of ζ m,n with arbitrary m and n, λ denotes the water depth of the depth eigenfunction G 1 . In Equation (23), the depth eigenfunction The conservation of momentum in Equation (7) as well as the barrier condition (9) can be expressed as: where and Based on the far-field conditions (11), (12) and the dynamic free-surface boundary condition (4), the far-field solutions of the velocity potential can be expressed as: and Combine Equation (29) and Equation (30) into Equation (18), we can obtain the following equations: and B M,n = 0 for n = 0, 1, . . . , N.
By using Equation (18), Equations (23) and (26) can be rewritten as: and for l = 0, 1, . . . , N and m = 1, 2, . . . , M − 1. Afterward, the EMM solution procedure can begin with Equations (33)-(37), which compose 2M(N + 1) linear equations and can be used to solve the 2M(N + 1) unknown A m,n and B m,n . Then, Equations (31) and (32) can be used to solve the reflection and transmission coefficients K R e iθ R and K T e iθ T . These complete the EMM solution procedures. In our study, the SuperLU library is used to solve the sparse matrix of the resulting linear Equation (21). It is noticeable that if L b = 0 or equivalently F m = 0, the proposed formulation reduces to the EMM formulation for cases of undulation bottom without barriers [20,27,31] and barriers on flat bottom [33].

Results
In order to validate the accuracy of EMM model, comparisons between the results by the proposed EMM and from the literature are made for cases of water wave scattering by barriers or undulated bottoms. The problem of convergence has been carefully taken into consideration by examining certain test cases.

Water Wave Scattering by a Single Barrier over Flat-Bottom Topography
First of all, let us consider the water waves scattering by a single surface-piercing thin barrier over a uniform bottom where h 1 = h 2 = 1.0m with the submerged length of barrier from the free surface being equal to a/h 1 = 0.5. In order to study the convergences on the number of evanescent modes N, the reflection |R| and transmission |T| coefficients are plotted against the dimensionless wavenumber kh in Figure 3. In the figure, it can be noted that the convergence is achieved for evanescent modes increase up to N = 50. Furthermore, the results are in good agreement with those in Porter and Evans [9] for the whole frequency range.  (31) and (32) can be used to solve the reflection and transmission coefficients These complete the EMM solution procedures. In our study, the SuperLU library is used to solve the sparse matrix of the resulting linear Equation (21). It is noticeable that if 0 b L  or equivalently 0 m F  , the proposed formulation reduces to the EMM formulation for cases of undulation bottom without barriers [20,27,31] and barriers on flat bottom [33].

Results
In order to validate the accuracy of EMM model, comparisons between the results by the proposed EMM and from the literature are made for cases of water wave scattering by barriers or undulated bottoms. The problem of convergence has been carefully taken into consideration by examining certain test cases.

Water Wave Scattering by a Single Barrier over Flat-Bottom Topography
First of all, let us consider the water waves scattering by a single surface-piercing thin barrier over a uniform bottom where 1 with the submerged length of barrier from the free surface being equal to . In order to study the convergences on the number of evanescent modes N , the reflection | | R and transmission | | T coefficients are plotted against the dimensionless wavenumber kh in Figure 3. In the figure, it can be noted that the convergence is achieved for evanescent modes increase up to 50 N  . Furthermore, the results are in good agreement with those in Porter and Evans [9] for the whole frequency range.  After considering the water wave scattering over a single surface-piercing barrier, the EMM is also implemented to simulates the case of scattered waves over a single bottomstanding thin barrier by uniform bottom. We assume a single barrier standing from the bottom with length b/h 1 = 0.3, where the water depth is h 1 = h 2 = 1.0 m as exhibited in Figure 4, where curves of |R| and |T| are plotted against dimensionless wavenumber kh 1 . The figure shows that convergence can be found for N = 50 and good agreements with those in Losada et al. [5] are clear. After considering the water wave scattering over a single surface-piercing barrier, the EMM is also implemented to simulates the case of scattered waves over a single bottom-standing thin barrier by uniform bottom. We assume a single barrier standing from the bottom with length as exhibited in Figure 4, where curves of | | R and | | T are plotted against dimensionless wavenumber 1 kh . The figure shows that convergence can be found for 50 N and good agreements with those in Losada et al. [5] are clear.

Water Wave Scattering by Barriers over Step-Type Bottom Topography
Then, we validate the proposed EMM by considering the scattering of oblique water waves by parallel thin barriers and a step bottom. Following Wang et al. [12], the two parallel surface-piercing barriers are of an equal submerged length , respectively. Additionally, the oblique incidence angle  is set to 30 . Figure 5 illustrates the relationship between the reflection | | R and dimensionless wavenumber 1 kh . Notably, numerical results of the present model with a number of modes 50 N  are in good agreement with those from literature [12]. Based on the results of the above examples, we will adopt 50 N  in the following of this study. Overall, this example demonstrates the proposed EMM can be applied for solving water wave scattering by multiple barriers over a step bottom.

Water Wave Scattering by Barriers over Step-Type Bottom Topography
Then, we validate the proposed EMM by considering the scattering of oblique water waves by parallel thin barriers and a step bottom. Following Wang et al. [12], the two parallel surface-piercing barriers are of an equal submerged length a/h 1 = 0.25 and the step separates the deeper and shallower waters with depths being h 1 = 1.0 m and h 2 = 0.25 m, respectively. Additionally, the oblique incidence angle γ is set to 30 • . Figure 5 illustrates the relationship between the reflection |R| and dimensionless wavenumber kh 1 . Notably, numerical results of the present model with a number of modes N = 50 are in good agreement with those from literature [12]. Based on the results of the above examples, we will adopt N = 50 in the following of this study. Overall, this example demonstrates the proposed EMM can be applied for solving water wave scattering by multiple barriers over a step bottom.

Water Wave Scattering over Undulated Slope Bottom
Then we applied the proposed method for solving water wave scattering over the undulated bottom. Following Porter and Porter [34], we consider the bottom defined by where α is parameter which defines the shape of the bottom. Figure 6 illustrates the reflection coefficients against dimensionless wavenumber kh 1 for the undulation bottom topography where h 1 = 1m is the upstream water depth. In the figure, it can be observed that the curves of |R| converge to the results computed by Porter and Porter [34] when step numbers increase to M = 50.

Water Wave Scattering over Undulated Slope Bottom
Then we applied the proposed method for solving water wave scattering over the undulated bottom. Following Porter and Porter [34], we consider the bottom defined by where  is parameter which defines the shape of the bottom. Figure 6 illustrates the reflection coefficients against dimensionless wavenumber

Water Wave Scattering over the Sinusoidal Rippled Bottom
Davies et al. [1,2] studied the Bragg reflection by the sinusoidal rippled bottoms, which provides the groundwork for considering Bragg refletions as one of the main impact factors to water scattering. Combining the proposed step approximation with these studies, the EMM is applied to calculate the resonant reflection coefficients by sinusoidal rippled bottom, which is defined by:

Water Wave Scattering over the Sinusoidal Rippled Bottom
Davies et al. [1,2] studied the Bragg reflection by the sinusoidal rippled bottoms, which provides the groundwork for considering Bragg refletions as one of the main impact factors to water scattering. Combining the proposed step approximation with these studies, the EMM is applied to calculate the resonant reflection coefficients by sinusoidal rippled bottom, which is defined by: where h 1 = 0.156 m is the water depth away from the ripples, K = 2π is the wavenumber of ripples and β is the number of sinusoidal ripples. When the case of two ripples (β = 2) are considered, Figure 7 depicts the reflection coefficients against 2k/K. In the figure, it is significant that the Bragg resonance happens when 2k/K ∼ 1. In the figure, it can be observed that the curves of |R| converge well with Davies and Heathershaw [1] as the step numbers increase to M = 100. These two examples validate the application of the EMM for solving problems of water wave scattering over undulated bottoms. Therefore, in this study, we will adopt for the undulated slope (38) and sinusoidal rippled (39) bottoms, respectively.

Discussion
In this section, the proposed EMM is applied to study the combined effects of multiple barriers and undulation bottom profiles for water wave scattering. We focus our attention on the undulated slope and sinusoidal rippled bottoms in the following.

A Surface-Piercing Barrier
In this examination, the bottom profile is defined by Equation (38)   These two examples validate the application of the EMM for solving problems of water wave scattering over undulated bottoms. Therefore, in this study, we will adopt M = 50 and M = 50β for the undulated slope (38) and sinusoidal rippled (39) bottoms, respectively.

Discussion
In this section, the proposed EMM is applied to study the combined effects of multiple barriers and undulation bottom profiles for water wave scattering. We focus our attention on the undulated slope and sinusoidal rippled bottoms in the following.

A Surface-Piercing Barrier
In this examination, the bottom profile is defined by Equation (38) as depicts in Figure 1. Additionally, a surface-piercing barrier is located at x = v with v/h 1 = 0.3. Four values of barrier length a/h 1 equal to 0, 0.1, 0.2 and 0.3 are considered. As illustrates in Figure 8, with the first three values of a/h 1 , with long waves which wavenumber kh 1 increases from 0 to 2, the reflection coefficient |R| firstly drops down before grow up with the increasing of the dimensionless wavenumber kh 1 . Different from the above three values, with a/h 1 = 0.3, |R| stably increases when water waves become shorter. Basically, the results indicate that the reflection is stronger when the barrier length becomes larger, especially for short waves. can be observed that the reflection becomes stronger when the barrier moves to the right of the undulated slope bottom, especially for medium waves, i.e., 1~1 .5 kh .  Then, we consider a surface-piercing barrier with length defined by a/h 1 = 0.2. Four values of barrier location v/h 1 equal to 0, 0.2, 0.4 and 0.6 are considered. Figure 9 illustrates the variation of |R| against dimensionless wavenumber kh 1 . From the figure, it can be observed that the reflection becomes stronger when the barrier moves to the right of the undulated slope bottom, especially for medium waves, i.e., kh 1 ∼ 1.5.

A Bottom-Standing Barrier
Then, another case of a bottom-standing barrier is considered. The barrier length is located at x = w with w/h 1 = 0.1. Four values of barrier length b/h 1 equal to 0, 0.2, 0.3 and 0.4 are considered. As demonstrates in Figure 10, the values of |R| are plotted against dimensionless wavenumber kh 1 . From the figure, It can be seen that the effects of barrier length are not significant when the barrier is located on the very left of the undulated slope bottom. Additionally, for the short wave limit kh 1 → ∞ , total transmissions occur for all cases.
The problem of barrier location is also considered. We assume that the barrier length equal to b/h 1 = 0.3. The barrier moves to the right of the undulated slope bottom with four barrier locations w/h 1 being equal to 0, 0.3, 0.4 and 0.5, respectively. Figure 11 illustrates the variation of |R| against dimensionless wavenumber kh 1 . From the figure, it can be observed that the reflection becomes stronger when the barrier moves to the right of the undulated slope bottom for almost all wavelengths.

Combined Effects of Surfacing-Piecing and Bottom-Standing Barriers
Sequentially, the combined effects of surface-piecing and bottom-standing barriers over the undulated slope bottom, defined by Equation (38), are studied. The surfacepiercing and bottom-standing barrier are located at x = v and x = w, respectively. Typically, we adopt v/h 1 = 0.4 and w/h 1 = 0.1. Additionally, their barrier lengths are set to be equal, as depicted in Figure 12. Basically, the results indicate that the reflection is stronger when the barrier length becomes larger, especially for short waves. In the short wave limit kh 1 → ∞ , total reflections occur when there are barriers (a/h 1 = b/h 1 > 0). can be observed that the reflection becomes stronger when the barrier moves to the right of the undulated slope bottom, especially for medium waves, i.e., 1~1 .5 kh .   The problem of barrier location is also considered. We assume that the barrier length equal to    Figure 12. Basically, the results indicate that the reflection is stronger when the barrier length becomes larger, especially for short waves. In the short wave limit 1 kh , total reflections occur when there are barriers (    Figure 12. Basically, the results indicate that the reflection is stronger when the barrier length becomes larger, especially for short waves. In the short wave limit 1 kh , total reflections occur when there are barriers (

Effect of Oblique Water Waves
In general, water waves could approach to coastal line from either normal incident or oblique incidence. To enlarge the working efficiency of structures, the main target of this section is focusing on the mutation of approach oblique incidences. The parameters of multiple barriers are typically chosen as a/h 1 = 0.2, v/h 1 = 0.4, b/h 1 = 0.3 and w/h 1 = 0.1 with a variety of approaches oblique water wave incidence γ. As depicted in Figure 13, reflection coefficients are drawn against dimensionless wavenumber kh 1 under the effect of incidences γ = 0 • , 30 • and 60 • , respectively. The resulting reflection coefficients are basically not significant with respect to the changes of the incident angles in this specific case. Further studies are still required.
In general, water waves could approach to coastal line from either normal incident or oblique incidence. To enlarge the working efficiency of structures, the main target of this section is focusing on the mutation of approach oblique incidences. The parameters of multiple barriers are typically chosen as with a variety of approaches oblique water wave incidence . As depicted in

Combined Effects of Thin Barriers and Sinusoidal Rippled Bottom
Then, the EMM is applied to solve water wave scattering problems by multiple barriers over the sinusoidal rippled bottom to study the Bragg resonance. The sinusoidal rippled bottom is set according to Equation (39) with the parameter defined therein if not otherwise mentioned.

A Single Surface-Piecing or Bottom-Standing Barrier
Let's begin with the water wave scattering by a surface-piercing or bottom-standing barrier over the prescribed sinusoidal rippled bottom. The surface-piercing or bottomstanding barrier is located at x v  and x w  , respectively. Firstly, a surface-piercing barrier is set at location  Figure 14a illustrates the variation of | | R against 2 / k K , which indicates that the effect of the surface-piercing barrier is to strengthen the reflection, especially for shorter waves by larger barriers. Secondly, a bottom-standing barrier is set at location  Figure 14b shows the variation of | | R Figure 13. Reflection coefficient |R| varying against kh 1 for a set of one bottom-standing barrier and one surface-piercing barrier by waves with different oblique incidence angles γ.

Combined Effects of Thin Barriers and Sinusoidal Rippled Bottom
Then, the EMM is applied to solve water wave scattering problems by multiple barriers over the sinusoidal rippled bottom to study the Bragg resonance. The sinusoidal rippled bottom is set according to Equation (39) with the parameter defined therein if not otherwise mentioned.

A Single Surface-Piecing or Bottom-Standing Barrier
Let's begin with the water wave scattering by a surface-piercing or bottom-standing barrier over the prescribed sinusoidal rippled bottom. The surface-piercing or bottomstanding barrier is located at x = v and x = w, respectively. Firstly, a surface-piercing barrier is set at location v/h 1 = 0.4 with four values of barrier length a/h 1 equal to 0, 0.2, 0.3 and 0.4. Figure 14a illustrates the variation of |R| against 2k/K, which indicates that the effect of the surface-piercing barrier is to strengthen the reflection, especially for shorter waves by larger barriers.
Secondly, a bottom-standing barrier is set at location w/h 1 = 0.1 with four values of barrier length b/h 1 equal to 0, 0.4, 0.5 and 0.6. Figure 14b shows the variation of |R| against 2k/K. In the figure, it is clear that the larger bottom-standing barrier gives a stronger reflection for all wavelengths.

Enhanced Bragg Resonance by Dual Periodic Barriers and Sinusoidal Rippled Bottom
Sequentially, we consider the enhanced Bragg resonance by dual barriers and sinusoidal rippled bottom. At first, we consider the dual periodic surface-piercing barriers located at x = v 1 and x = v 2 with v 1 /h 1 = π/2K and v 2 /h 1 = 5π/2K, respectively. The lengths of the barriers are equal as a 1 /h 1 = a 2 /h 1 and set to four values of 0, 0.1, 0.2 and 0.3. Figure 15 depicts the variation of |R| against 2k/K. In the figure, it is significant that the main effect of dual periodic barriers is to enhance both the primary (2k/K ∼ 1) and secondary (2k/K ∼ 2) Bragg resonances, especially for larger barriers. Here, the enhancement of the secondary Bragg resonances is especially observable for the case of a 1 /h 1 = a 2 /h 1 = 0.3.  Secondly, we study the effect of dual periodic bottom-standing barriers located exactly the same as the previous case. The lengths of the barriers are also equal as b 1 /h 1 = b 2 /h 1 and set to four values of 0, 0.4, 0.5 and 0.6. Figure 16 shows the variation of |R| against 2k/K, in which the enhanced primary and secondary Bragg resonances are also significant, especially for larger barriers.
Sequentially, we consider the enhanced Bragg resonance by dual barriers and sinusoidal rippled bottom. At first, we consider the dual periodic surface-piercing barriers located at 1 x v  and 2 x v   Figure 15 depicts the variation of | | R against 2 / k K . In the figure, it is significant that the main effect of dual periodic barriers is to enhance both the primary ( 2 /~1 k K ) and secondary ( 2 / 2 k K ) Bragg resonances, especially for larger barriers. Here, the enhancement of the secondary Bragg resonances is especially observable for the case of   From these two cases, it is clear that the effects of both the dual periodic surfacepiercing and bottom-standing barriers are to enhance the Bragg resonances since they are located with the same frequency as the sinusoidal rippled bottom. The enhancement is significant, especially for the second Bragg resonance for cases of larger barrier lengths.
Additionally, the enhanced Bragg resonances are obvious compared to the cases of single barriers in the previous subsection.

Enhanced Bragg Resonance by Four Periodic Barriers and Sinusoidal Rippled Bottom
In order to demonstrate the enhancements of the Bragg resonance, we consider the problems of water wave scattering by four periodic barries over sinusoidal rippled bottom with the number of ripples β = 4. Both the four periodic barries and the sinusoidal rippled bottom are set to be of the same frequency K. In other words, the four periodic surfacepiercing barriers are located at x = v 1 , x = v 2 , x = v 3 and x = v 4 with v 1 /h 1 = π/2K, v 2 /h 1 = 5π/2K, v 3 /h 1 = 9π/2K and v 4 /h 1 = 13π/2K, respectively. In addition, the four bottom-standing barriers are located exactly the same. Figure 17 depicts the variation of |R| against 2k/K, where the enhanced primary and secondary Bragg resonances are still significant for both the cases of surfacing-piecing and bottom-standing barriers over the sinusoidal rippled bottom. The enhancements are more considerable, especially for the issues of larger barriers.

Enhanced Bragg Resonance by Multiple Barriers and Sinusoidal Rippled Bottom under Effects of Oblique Water Waves
In order to study the effects of oblique incidences, the dual thin barrier cases with a 1 /h 1 = a 2 /h 1 = 0.3 and b 1 /h 1 = b 2 /h 1 = 0.5 (in Figure 16) are extended by considering different incident angles. As demonstrated in Figure 18a, |R| are plotted against 2k(cos γ)/K for the surface-piercing case, where the enhanced primary (2k(cos γ)/K = 1) and secondary Bragg resonances (2k(cos γ)/K = 2) are still significant for the oblique water wave scattering by dual surfacing-piecing barriers over the sinusoidal rippled bottom. Additionally, the numerical results indicated that the cases with larger incident angles result in less intensive primary Bragg resonances. Then the effects of oblique water waves on the dual bottom-standing barrier over the sinusoidal rippled bottom are studied as depicted in in Figure 18b. In the figure, the Bragg's law can also be found to be significant and the cases with larger incident angles give less intensive primary Bragg resonances.

Enhanced Bragg Resonance by Multiple Barriers and Sinusoidal Rippled Bottom under Effects of Oblique Water Waves
In order to study the effects of oblique incidences, the dual thin barrier cases with

Conclusions
In this paper, we proposed the EMM formulation to solve the scattering of oblique water waves by multiple thin barriers over undulation bottoms. The bottom topographies were sliced into shelves separated by steps in the solution procedures. Additionally, the solutions were composed of eigenfunctions and solved by the conservations of mass and momentum. The sparse-matrix solver, SuperLU, solved the resulted system of linear equations. The proposed EMM formulation can be reduced to the traditional EMM formulation for cases of undulation bottom without barriers and barriers on a flat bottom. Sequentially, the EMM was validated by the problems of water wave scattering by a single surface-piecing or bottom-standing barrier, dual barriers over a step, undulated slope and sinusoidal rippled bottoms. Then, the effects of barriers on the undulated slope bottom were discussed. Finally, the EMM was applied for studying the impact of barriers on the Bragg resonance by the sinusoidal rippled bottoms. It was found that the effects of the periodic barriers are to enhance the Bragg resonances since they are located with the same frequency as the sinusoidal rippled bottom. The enhanced Bragg resonances are obvious as compared with the cases of single barries over the sinusoidal rippled bottom. Additionally, the Bragg's law of oblique waves is validated for the prescribed enhanced resonances of periodic barriers and sinusoidal rippled bottoms. Numerical results indicated that the cases with larger incident angles result in less intensive Bragg resonances.

Conclusions
In this paper, we proposed the EMM formulation to solve the scattering of oblique water waves by multiple thin barriers over undulation bottoms. The bottom topographies were sliced into shelves separated by steps in the solution procedures. Additionally, the solutions were composed of eigenfunctions and solved by the conservations of mass and momentum. The sparse-matrix solver, SuperLU, solved the resulted system of linear equations. The proposed EMM formulation can be reduced to the traditional EMM formulation for cases of undulation bottom without barriers and barriers on a flat bottom. Sequentially, the EMM was validated by the problems of water wave scattering by a single surface-piecing or bottom-standing barrier, dual barriers over a step, undulated slope and sinusoidal rippled bottoms. Then, the effects of barriers on the undulated slope bottom were discussed. Finally, the EMM was applied for studying the impact of barriers on the Bragg resonance by the sinusoidal rippled bottoms. It was found that the effects of the periodic barriers are to enhance the Bragg resonances since they are located with the same frequency as the sinusoidal rippled bottom. The enhanced Bragg resonances are obvious as compared with the cases of single barries over the sinusoidal rippled bottom. Additionally, the Bragg's law of oblique waves is validated for the prescribed enhanced resonances of periodic barriers and sinusoidal rippled bottoms. Numerical results indicated that the cases with larger incident angles result in less intensive Bragg resonances.