Accurate Solutions to Water Wave Scattering by Vertical Thin Porous Barriers

Abstract

The water wave scattering by vertical thin porous barriers is accurately solved in this study. Two typical structures of a surface-piercing barrier and a submerged bottom-standing barrier are considered. The solution procedure is based on the multi-term Galerkin method, in which the pressure jump across a porous barrier is expanded in a set of basis functions involving the Chebychev polynomials. Then, the square-root singularity of fluid velocity at the edge of the porous barrier is correctly modeled. The present solutions have the merits of very rapid convergence. Accurate results for both the reflection and the transmission coefficients and wave forces are presented. This study not only gives a promising procedure to tackle wave interaction with vertical thin porous barriers but also provides a reliable benchmark for complicated numerical solutions.

Academic Editor:Erik Cuevas

Shandong Provincial Key Laboratory of Ocean Engineering, Ocean University of China, Qingdao 266100, China

Received 1 April 2015; Revised 11 September 2015; Accepted 20 September 2015; 12 October 2015

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Vertical thin plates have been used as simple breakwaters in coastal engineering due to their merits of simple structure, constructing convenience, and low engineering cost. In some cases, a vertical barrier (plate) can play an important role in protecting coast from erosion as well as various activities along the coastline from wave attack. The practical vertical barriers are often surface-piercing and supported on piles or submerged below the ocean surface and bottom-standing [1-5].

An abundance of analytical studies for linear water wave scattering by vertical thin solid barriers have been conducted. Ursell [6] carried out an early study on normally incident wave transmission by a vertical thin barrier in deep water. Evans [7] developed expressions for the first- and second-order forces and moments on a submerged vertical thin wall based on linear potential theory. Morris [8] developed a variational approach to normally incident water interaction with two nonsymmetric vertical thin barriers. Losada et al. [9] and Abul-Azm [10] used matched eigenfunction expansion method to develop analytical solutions for obliquely and normally incident wave scattering by vertical thin barriers with different configurations, respectively. Kriebel and Bollmann [11] compared the transmission coefficients of surface-piercing solid barriers calculated by the wave power theory (an approximate method) and the matched eigenfunction expansion method. Porter and Evans [1] investigated oblique wave scattering by various partial solid barriers using a multi-term Galerkin method and obtained the upper and lower boundaries of reflection and transmission coefficients with extremely high accuracy. Banerjea et al. [12] also used the multi-term Galerkin method to obtain accurate solutions for oblique wave scattering by single and double submerged vertical solid barriers with gaps. Moreover, the multi-term Galerkin method has been used by Evans and Porter [13], Martins-Rivas and Mei [14], and Chang et al. [15] for three-dimensional vertical thin walls. Compared with matched eigenfunction expansion method, the beauty of multi-term Galerkin method is that it can correctly model the square-root singularity of fluid velocity near the edge of vertical plate. Then, extremely accurate results of hydrodynamic quantities can be obtained.

Besides good shelter function with required transmission coefficient, smaller wave force acting on a barrier is expected in engineering design for ensuring the safety of the structure. Thus, porous barriers have been proposed and used as breakwaters in coastal engineering. After generalizing the study of Tuck [16] for permeable thin barriers in deep water, Macaskill [17] examined normally incident wave interaction with a permeable thin barrier in finite water depth by solving integral equations. Sahoo [18] applied a perturbation analysis to study water wave scattering by porous barriers with small porosity effect in deep water. Isaacson et al. [2] used matched eigenfunction expansion method to develop an analytical solution for normally incident wave interaction with a partially immersed porous barrier and validated their solution using experimental data. Lee and Chwang [19] examined wave scattering and radiation by vertical thin porous barriers with four different configurations based on the eigenfunction expansion and least square methods. They found that due to the energy dissipation by holes, the reflection coefficient and wave force of a porous barrier are smaller than those of a solid barrier. Sahoo et al. [20] generalized the solution of Lee and Chwang [19] to obliquely incident waves. Recently, Karmakar and Guedes Soares [21] developed analytical solution for wave transformation due to multiple bottom-standing flexible porous barriers using eigenfunction expansion method. Gayen and Mondal [22] studied normally incident wave reflection and transmission by a submerged inclined porous barrier by solving a second kind hypersingular integral equation.

As mentioned above, a variety of solutions with respect to vertical thin porous barriers have been developed. But rapidly convergent solutions for thin porous barriers are still lacking. In particular, for multiple barriers, a rapidly convergent solution is of great significance for parameter studies. Also, an accurate analytical solution can serve as a reliable benchmark for complicated numerical potential flow solutions. Thus, we will develop extremely accurate solutions for vertical thin porous barriers in this study. As a first step, we will consider water wave scattering by two simple and typical thin porous barriers, that is, a surface-piercing barrier and a submerged bottom-standing barrier. We use the multi-term Galerkin method [1] to incorporate the square-root singularity of fluid velocity near the edge of porous barrier into the solution. The present solution procedure can be easily extended to more complicated structures, such as various multiple porous barriers.

The governing equation and boundary conditions of the present problem are given in the following section. In Section 3, the solving procedure using multi-term Galerkin method is detailed. In Section 4, the rapid convergence and high accuracy of the present solution are shown by comparing with matched eigenfunction expansion solution. Then, we present some typical curves of hydrodynamic quantities to show the effects of major barrier parameters. Finally, the main conclusions of this study are drawn.

2. The Boundary Value Problem

The idealized sketch of water wave scattering by vertical porous barriers is given in Figure 1. The vertical porous barrier is surface-piercing (Type 1) or totally submerged and bottom-standing (Type 2). The thickness of the porous barrier is assumed to be zero, as it is very small compared with the incident wavelength. The water depth is [figure omitted; refer to PDF] . We use [figure omitted; refer to PDF] to denote the area of barrier and [figure omitted; refer to PDF] to denote the fluid area of gap. A Cartesian coordinate system with the [figure omitted; refer to PDF] -axis taking vertically upwards along the barrier and the [figure omitted; refer to PDF] -axis along the still water level is used for mathematical descriptions. The incident waves propagate along the positive [figure omitted; refer to PDF] -direction.

Figure 1: Idealized sketch of water wave scattering by vertical porous barriers.

(a) A surface-piercing porous barrier (Type 1)

[figure omitted; refer to PDF]

(b) A submerged bottom-standing porous barrier (Type 2)

[figure omitted; refer to PDF]

It is assumed that the fluid is inviscid and incompressible, the fluid motion is irrotational, and the incident harmonic wave is of small amplitude with angular frequency [figure omitted; refer to PDF] . Then, a velocity potential [figure omitted; refer to PDF] can be used to describe the fluid motion. We separate out the time factor [figure omitted; refer to PDF] and then write the velocity potential and the dynamic pressure as [figure omitted; refer to PDF] where Re denote the real part of function; [figure omitted; refer to PDF] ; [figure omitted; refer to PDF] is the gravitational acceleration; [figure omitted; refer to PDF] is the incident wave amplitude; [figure omitted; refer to PDF] is the time; [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are, respectively, the spatial velocity potential and dynamic pressure.

For the convenience of study, the whole fluid domain is divided into two sub-regions: left sub-region [figure omitted; refer to PDF] and right sub-region [figure omitted; refer to PDF] . The spatial velocity potential in each sub-region satisfies the Laplace equation: [figure omitted; refer to PDF] where the subscript [figure omitted; refer to PDF] denotes variables in sub-region [figure omitted; refer to PDF] . The velocity potentials also satisfy the following boundary conditions on the free surface, the water bottom, and the far fields: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the incident wave number and [figure omitted; refer to PDF] is the velocity potential of incident waves.

At the edge of the porous barrier, the gradient of the velocity potential has a square-root singularity ([23, Eq. (2.85)] and [24, Eq. (1e)]): [figure omitted; refer to PDF] where [figure omitted; refer to PDF] for Type 1 barrier and [figure omitted; refer to PDF] for Type 2 barrier. In order to develop an accurate solution for the present problem, the square-root singularity of fluid velocity must be correctly modeled. This will be detailed in the following section.

Besides preceding boundary conditions, the velocity potentials satisfy the transmission boundary conditions at the interface of two sub-regions: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is a complex effect parameter of the porous barrier [25]. Equations (8) and (9) denote, respectively, that the horizontal fluid velocities and the dynamic pressures are continuous at the gap. Equation (10) is the boundary condition on the porous barrier, which was developed by Yu [25] based on the classical porous medium model of Sollitt and Cross [26]. The first equals sign in (10) denotes that the normal fluid velocity is continuous at the porous barrier. The second equals sign in (10) denotes that the normal fluid velocity through the porous barrier is proportional to the pressure jump across the barrier. Then, wave energy dissipation is introduced by (10).

According to Yu [25], the complex porous effect parameter [figure omitted; refer to PDF] can be estimated by [figure omitted; refer to PDF] , where [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are the porosity, the linearized resistance coefficient, the inertial coefficient, and the thickness of the porous plate (porous medium), respectively. The value of the inertial coefficient [figure omitted; refer to PDF] may be simply treated as unity [27, 28]. The value of the linearized resistance coefficient [figure omitted; refer to PDF] must be determined by experimental tests. Suh et al. [28] proposed an empirical formula for the resistance coefficient: [figure omitted; refer to PDF] . The real and imaginary parts of [figure omitted; refer to PDF] denote the linearized resistance effect and the inertial effect of a porous plate, respectively. The resistance effect brings wave energy dissipation, and the inertial effect produces the phase shift of wave motion. When [figure omitted; refer to PDF] equals zero, the barrier is impermeable. If [figure omitted; refer to PDF] tends to infinity, the barrier becomes entirely transparent (no porous barrier). This can be observed by rewritten (10) as: [figure omitted; refer to PDF] , where [figure omitted; refer to PDF] . In our computations, a very large value of [figure omitted; refer to PDF] is adopted for entirely transparent wall.

3. Methods of Solution

3.1. Expressions of Velocity Potentials

By the separation of variables, the velocity potentials satisfying the governing equation in (2) and the relevant boundary conditions in (3)-(6) can be written as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are unknown expansion complex coefficients. The wave numbers [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are the positive real roots of the following dispersion relation: [figure omitted; refer to PDF] In (11), the vertical eigenfunctions [figure omitted; refer to PDF] are given by [figure omitted; refer to PDF] It is noted that the eigenfunctions [figure omitted; refer to PDF] are orthogonal over [figure omitted; refer to PDF] : [figure omitted; refer to PDF] and integrals of their own square over [figure omitted; refer to PDF] are [figure omitted; refer to PDF] Then, applying the transmission boundary conditions in (8)-(10), the unknown expansion coefficients in (11) can be determined. For convenience, these transmission conditions are rewritten as [figure omitted; refer to PDF]

Substituting the expressions for the velocity potentials in (11) into (16), we have [figure omitted; refer to PDF] Multiplying both sides of (19) by [figure omitted; refer to PDF] (z ) and integrating with respect to [figure omitted; refer to PDF] over [-d , 0] and then using (14), we obtain [figure omitted; refer to PDF]

We note that the two conditions in (17) and (18) are known as dual series relations [29]. They can be combined into a mixed boundary condition along the whole water depth (see (A.1)), which can be solved by a least square method [30-32]. Here we use a similar procedure as that for (16) to transform (17) and (18) into a system of linear equations, and then determine all the unknown coefficients. This is just the traditional matched eigenfunction expansion method and is introduced in the Appendix. However, the square-root singularity of fluid velocity is not incorporated into such a solution procedure. As a result, the convergence of the solution is slow. In order to obtain extremely accurate analytical solution, we use the multi-term Galerkin method [1], instead of traditional matched eigenfunction expansion method, to solve the present problem.

3.2. Multi-Term Galerkin Methods

Following Porter and Evans [1, Eqs. (2.52) and (2.63)], we expand the jump of velocity potentials, that is, the pressure jump, at the porous barrier as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] are unknown expansion coefficients, and [figure omitted; refer to PDF] is given by [figure omitted; refer to PDF] with [figure omitted; refer to PDF] for Type 1 barrier, and [figure omitted; refer to PDF] for Type 2 barrier. In (23) and (24), [figure omitted; refer to PDF] is the [figure omitted; refer to PDF] th order Chebychev polynomial of the second kind. According to the expressions of [figure omitted; refer to PDF] , (22) satisfies the free surface condition in (3) and the square-root singularity at the lower edge of Type 1 barrier. Equation (24) satisfies the water bottom condition in (4) and the square-root singularity at the upper edge of Type 2 barrier.

Substituting the expressions for the velocity potentials in (11) into (21) and multiplying both sides of the obtained new equation by [figure omitted; refer to PDF] and then integrating with respect to [figure omitted; refer to PDF] over [figure omitted; refer to PDF] and using (14) and (15), we have [figure omitted; refer to PDF] where [figure omitted; refer to PDF] for Type 1 barrier, and [figure omitted; refer to PDF] for Type 2 barrier. In (26) and (27), [figure omitted; refer to PDF] and [figure omitted; refer to PDF] denote the Bessel function and the modified Bessel function of order [figure omitted; refer to PDF] , respectively. Substituting (11) with (20) and (25) into (17), multiplying both sides of the obtained new equation by [figure omitted; refer to PDF] , and then integrating with respect to [figure omitted; refer to PDF] over Γ, we obtain [figure omitted; refer to PDF]

We solve (28) by truncating [figure omitted; refer to PDF] and [figure omitted; refer to PDF] after [figure omitted; refer to PDF] terms and obtain the values of [figure omitted; refer to PDF] . Then, all the expansion coefficients [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] in velocity potentials are determined using (20) and (25).

3.3. Hydrodynamic Quantities

The reflection and transmission coefficients of the vertical porous barrier are calculated, respectively, by [figure omitted; refer to PDF] The energy loss coefficient is defined as [figure omitted; refer to PDF] The energy dissipation is introduced by the porous boundary condition in (10), which brings pressure loss on the porous barrier. If the barrier is impermeable [figure omitted; refer to PDF] , the energy loss coefficient is always zero.

The dynamic pressure jump on the porous barrier is calculated by the linear Bernoulli equation [figure omitted; refer to PDF] , where [figure omitted; refer to PDF] is the water density. Integrating the dynamic pressure jump along the porous barrier, we get the magnitude of horizontal wave force acting on the barrier: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] for Type 1 barrier, and [figure omitted; refer to PDF] for Type 2 barrier. The dimensionless wave force on the barrier is defined as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the wave force acting on per-unit width in an impermeable vertical wall extending from seabed to the still water level.

4. Results

The second part in the left hand side of (28) is still an infinite series after truncating [figure omitted; refer to PDF] and [figure omitted; refer to PDF] . Thus, the values of [figure omitted; refer to PDF] and [figure omitted; refer to PDF] need to be estimated by truncating [figure omitted; refer to PDF] . Porter and Evans [1] have detailed the computation method of this type of series for solid barriers. In this study, all the series are first computed taking 800 terms, and then they are improved by applying the asymptotic values of Bessel functions as those used in Porter and Evans [1]. In order to obtain the required accuracy, we use 10⁸ terms in computations of the improvements.

4.1. Special Cases of Solid Barriers

We first consider the special cases of solid barriers when the porous effect parameter [figure omitted; refer to PDF] is zero. The calculated results of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] for two types of barriers are given in Tables 1 and 2. Here, the truncated number of [figure omitted; refer to PDF] is 7. It can be seen from Tables 1 and 2 that the wave energy conservation ( [figure omitted; refer to PDF] ) is well satisfied for solid barriers. For comparisons, the results obtained using Tables 1 and 2 and (2.25) and (2.26) in Porter and Evans [1] are also added in these tables. For solid barriers, the fluid velocity at the gap Λ can also be expanded using a similar method to that in (21). The two approaches based on expanding the pressure jump at G and the fluid velocity at Λ for solid barrier can give the lower and upper bounds of hydrodynamic quantities, respectively. Thus, the results of Porter and Evans [1] for solid barriers can be treated as the exact solution. It can be seen from Tables 1 and 2 that the present results are the same as that of Porter and Evans [1] for five-figure accuracy.

Table 1: Results of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] for Type 1 barrier at [figure omitted; refer to PDF] and [figure omitted; refer to PDF] .

[figure omitted; refer to PDF]	The present results	Porter and Evans (1995) [1]
[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
0.1	0.439284	0.898348	0	0.43928	0.89835
0.4	0.349212	0.937044	0	-	-
0.5	0.353987	0.935250	0	0.35399	0.93525
0.8	0.493434	0.869783	0	-	-
0.9	0.602548	0.798083	0	0.60255	0.79808

Table 2: Results of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] for Type 2 barrier at [figure omitted; refer to PDF] .

[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	The present results	Porter and Evans (1995) [1]
[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
0.1	0.1	0.005598	0.999984	0	-	-
0.2	0.05	0.007817	0.999969	0	0.00782	0.99997
0.5	0.2	0.082591	0.996584	0	-	-
0.6	0.3	0.150706	0.988579	0	0.15071	0.98858
0.9	1.8	0.542662	0.839951	0	0.54266	0.83995

4.2. Convergence Examination for Porous Barriers

The main purpose of the results presented below is to demonstrate the high accuracy and rapid convergence of the present multi-term Galerkin solution for porous barriers. The computed results of the reflection coefficient, the transmission coefficient, and the wave force, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , for the surface-piercing porous barrier (Type 1) and the submerged porous barrier (Type 2) at different truncated number [figure omitted; refer to PDF] are listed in Tables 3 and 4, respectively. It can be seen from these tables that for both barriers with different sizes, the convergence of the solution is very rapid. Results with six-figure accuracy can be generally obtained using the truncating number of [figure omitted; refer to PDF] . Results with higher accuracy can also be obtained using larger truncated numbers. Besides results shown in Tables 3 and 4, our extensive numerical experiments indicate that other values of wave numbers and barrier sizes produce equally accurate results for [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] . The high accuracy of the present solution is due to the fact that the nature of the square-root singularity near the edge of the porous barrier is correctly modeled. We note that results with six-figure accuracy are presented to show the rapid convergence of the solution. One may adopt less Galerkin terms to obtain results with two- or three-figure accuracy.

Table 3: Results of the multi-term Galerkin solution for Type 1 porous barrier at [figure omitted; refer to PDF] and [figure omitted; refer to PDF] .

Truncated number [figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
0	0.222895	0.808722	0.139259	0.201558	0.835137	0.173011	0.267398	0.768831	0.261473
1	0.223158	0.809090	0.139902	0.202015	0.835583	0.174000	0.271319	0.766307	0.266388
2	0.223145	0.809113	0.139885	0.201999	0.835609	0.173977	0.271375	0.766303	0.266405
3	0.223144	0.809114	0.139884	0.201998	0.835611	0.173976	0.271376	0.766305	0.266409
4	0.223144	0.809114	0.139884	0.201998	0.835611	0.173976	0.271376	0.766305	0.266408
5	0.223144	0.809114	0.139884	0.201998	0.835611	0.173976	0.271376	0.766306	0.266408
6	0.223143	0.809114	0.139884	0.201998	0.835611	0.173975	0.271376	0.766306	0.266408
7	0.223143	0.809114	0.139884	0.201998	0.835611	0.173975	0.271376	0.766306	0.266408
9	0.223144	0.809114	0.139884	0.201998	0.835611	0.173975	0.271376	0.766306	0.266408
10	0.223144	0.809114	0.139884	0.201998	0.835611	0.173975	0.271376	0.766306	0.266408

Table 4: Results of the multi-term Galerkin solution for Type 2 porous barrier at [figure omitted; refer to PDF] and [figure omitted; refer to PDF] .

Truncated number [figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
0	0.003551	0.998938	0.007621	0.069920	0.975171	0.071480	0.175771	0.905392	0.176347
1	0.003548	0.998943	0.007614	0.069879	0.975289	0.071435	0.177309	0.905034	0.177850
2	0.003548	0.998943	0.007614	0.069878	0.975290	0.071434	0.177348	0.905026	0.177889
3	0.003548	0.998943	0.007614	0.069878	0.975290	0.071434	0.177349	0.905026	0.177891
4	0.003548	0.998943	0.007621	0.069878	0.975290	0.071434	0.177349	0.905026	0.177891
5	0.003548	0.998943	0.007614	0.069878	0.975290	0.071434	0.177349	0.905026	0.177891
6	0.003548	0.998943	0.007614	0.069878	0.975290	0.071434	0.177349	0.905026	0.177891
7	0.003548	0.998943	0.007614	0.069878	0.975290	0.071434	0.177349	0.905026	0.177891
9	0.003548	0.998943	0.007614	0.069878	0.975290	0.071434	0.177349	0.905026	0.177891

For comparisons, the corresponding results calculated by the matched eigenfunction expansion solution (see the Appendix) for porous barriers are listed in Tables 5 and 6. It can be seen from these tables that, after using a large truncated number of [figure omitted; refer to PDF] , the matched eigenfunction expansion solution can only give results with three-figure accuracy at best. If we adopt larger [figure omitted; refer to PDF] , the convergence is still slow. The reason is that, for the matched eigenfunction expansion solution, fluid velocity singularity near the edge of the porous barrier is not considered.

Table 5: Results of the matched eigenfunction expansion solution for Type 1 porous barrier at [figure omitted; refer to PDF] and [figure omitted; refer to PDF] .

Truncated number [figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
20	0.252275	0.771745	0.152222	0.209811	0.826342	0.180059	0.280401	0.754224	0.274120
50	0.238389	0.789733	0.147991	0.205739	0.831183	0.177159	0.275419	0.760897	0.270155
100	0.232329	0.797420	0.145332	0.204131	0.833086	0.175876	0.273561	0.763391	0.268507
200	0.228312	0.802531	0.143226	0.203180	0.834212	0.175066	0.272607	0.764665	0.267320
400	0.226035	0.805435	0.141870	0.202652	0.834837	0.174594	0.272048	0.765411	0.267082
600	0.225193	0.806507	0.141332	0.202457	0.835067	0.174415	0.271847	0.765679	0.266885

Table 6: Results of the matched eigenfunction expansion solution for Type 2 porous barrier at [figure omitted; refer to PDF] and [figure omitted; refer to PDF] .

Truncated number [figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
20	0.005142	0.998193	0.009922	0.076300	0.971694	0.077025	0.193357	0.889137	0.192532
50	0.004270	0.998621	0.008808	0.072877	0.973634	0.074184	0.184363	0.898210	0.1814523
100	0.003963	0.998760	0.008348	0.071573	0.974358	0.073031	0.181097	0.901417	0.181488
200	0.003774	0.998844	0.008036	0.070812	0.974778	0.072333	0.179449	0.903013	0.179928
400	0.003672	0.998889	0.007855	0.070393	0.975008	0.071937	0.178491	0.903934	0.179007
600	0.003636	0.998905	0.007815	0.070239	0.975092	0.071790	0.178149	0.904262	0.178675

4.3. Theoretical Curves of Hydrodynamic Quantities

Figure 2 gives the variations of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] for a surface-piercing porous barrier (Type 1) with dimensionless wave number [figure omitted; refer to PDF] at a fixed porous effect parameter [figure omitted; refer to PDF] . It is observed that for a fixed relative barrier height [figure omitted; refer to PDF] , [figure omitted; refer to PDF] and [figure omitted; refer to PDF] increase and [figure omitted; refer to PDF] decreases with the increasing value of [figure omitted; refer to PDF] . When the value of [figure omitted; refer to PDF] is large enough, [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] all tend constants. It can also be observed from Figure 2 that, for a fixed [figure omitted; refer to PDF] , the reflection coefficient and the wave force both decease with the decreasing relative barrier height [figure omitted; refer to PDF] , while the transmission coefficient increases with the deceasing [figure omitted; refer to PDF] .

Figure 2: Variations of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] versus [figure omitted; refer to PDF] for a surface-piercing porous barrier (Type 1) at [figure omitted; refer to PDF] and different a /d .

(a) The reflection and transmission coefficients, [figure omitted; refer to PDF] and [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(b) The dimensionless wave force [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

Figure 3 shows the variations of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] for a submerged porous barrier (Type 2) with dimensionless wave number [figure omitted; refer to PDF] at a fixed porous effect parameter [figure omitted; refer to PDF] . It can be seen that the variations of the dimensionless wave force [figure omitted; refer to PDF] of Type 2 barrier with the relative barrier height [figure omitted; refer to PDF] are similar to those of Type 1 barrier. However, the reflection and transmission coefficients attain, respectively, a maximum value and a minimum value with the increasing wave number [figure omitted; refer to PDF] . This is in fact a common feature for wave motion over submerged barriers.

Figure 3: Variations of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] versus [figure omitted; refer to PDF] for a submerged porous barrier (Type 2) at [figure omitted; refer to PDF] and different b /d .

(a) The reflection and transmission coefficients, [figure omitted; refer to PDF] and [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(b) The dimensionless wave force [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

Finally, we examine the effects of the porous effect parameter [figure omitted; refer to PDF] on [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] . For Type 1 porous barrier, we fix the relative barrier height at [figure omitted; refer to PDF] and illustrate the results of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] and [figure omitted; refer to PDF] at different values of [figure omitted; refer to PDF] in Figure 4. When the porous effect parameter [figure omitted; refer to PDF] increases at a fixed [figure omitted; refer to PDF] , the reflection coefficient and the wave force on the barrier both decrease significantly. When [figure omitted; refer to PDF] and [figure omitted; refer to PDF] , the transmission and reflection coefficients of porous barrier are both smaller than those of impermeable barrier. This is due to the additional wave energy dissipation caused by the holes in porous barrier. For Type 2 porous barrier, we fix the relative barrier height at [figure omitted; refer to PDF] and give the results of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] at different values of [figure omitted; refer to PDF] in Figure 5. It can be seen that the effects of [figure omitted; refer to PDF] on the [figure omitted; refer to PDF] and [figure omitted; refer to PDF] of Type 2 barrier are similar to those of Type 1 barrier.

Figure 4: Variations of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] versus [figure omitted; refer to PDF] for a surface-piercing porous barrier (Type 1) at [figure omitted; refer to PDF] and different G .

(a) The reflection and transmission coefficients, [figure omitted; refer to PDF] and [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(b) The dimensionless wave force [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

Figure 5: Variations of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] versus [figure omitted; refer to PDF] for a submerged porous barrier (Type 2) at [figure omitted; refer to PDF] and different G .

(a) The reflection coefficient [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(b) The transmission coefficient [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

(c) The dimensionless wave force [figure omitted; refer to PDF]

[figure omitted; refer to PDF]

5. Conclusions

This study has developed accurate solutions for water wave scattering by two types of thin vertical porous barriers based on the linear potential theory. The multi-term Galerkin method has been used to take into account the square-root singularity near the edge of the porous barrier. We have shown that the multi-term Galerkin method worked very well for porous barriers. The convergence of the present solution was very rapid in comparison with the traditional matched eigenfunction expansion solution. Accurate results of various hydrodynamic quantities have been presented. The reflection and transmission coefficients of porous barrier may be both smaller than those of impermeable barrier. The wave force acting on a porous barrier is significantly reduced compared to the corresponding solid barrier. This should be favorable in some coastal engineering applications. The present study gives a reliable benchmark for other numerical solutions with respect to porous thin walls. Most of all, this study provides a promising procedure for solving wave interaction with vertical thin porous barriers.

Acknowledgment

This work was supported by the National Natural Science Foundation of China (Grant nos. 51322903 and 51279224).

Conflict of Interests

The authors declare no conflict of interests.

References

[1] R. Porter, D. V. Evans, "Complementary approximations to wave scattering by vertical barriers," Journal of Fluid Mechanics , vol. 294, pp. 155-180, 1995.

[2] M. Isaacson, S. Premasiri, G. Yang, "Wave interactions with vertical slotted barrier," Journal of Waterway, Port, Coastal and Ocean Engineering , vol. 124, no. 3, pp. 118-126, 1998.

[3] K.-D. Suh, S. Shin, D. T. Cox, "Hydrodynamic characteristics of pile-supported vertical wall breakwaters," Journal of Waterway, Port, Coastal and Ocean Engineering , vol. 132, no. 2, pp. 83-96, 2006.

[4] H. Oumeraci, G. Koether, P. Lynett, "Hydraulic performance of a submerged wave absorber for coastal protection," Nonlinear Wave Dynamics , pp. 31-65, World Scientific Publishing, Hackensack, NJ, USA, 2009.

[5] Z. H. Huang, Y. C. Li, Y. Liu, "Hydraulic performance and wave loadings of perforated/slotted coastal structures: a review," Ocean Engineering , vol. 38, no. 10, pp. 1031-1053, 2011.

[6] F. Ursell, "The effect of a fixed vertical barrier on surface waves in deep water," Proceedings of the Cambridge Philosophical Society , vol. 43, no. 3, pp. 374-382, 1947.

[7] D. V. Evans, "Diffraction of water waves by submerged vertical plate," Journal of Fluid Mechanics , vol. 40, no. 3, pp. 433-451, 1970.

[8] C. A. N. Morris, "A variational approach to an unsymmetric water-wave scattering problem," Journal of Engineering Mathematics , vol. 9, no. 4, pp. 291-300, 1975.

[9] I. J. Losada, M. A. Losada, A. J. Roldan, "Propagation of oblique incident waves past rigid vertical thin barriers," Applied Ocean Research , vol. 14, no. 3, pp. 191-199, 1992.

[10] A. G. Abul-Azm, "Wave diffraction through submerged breakwaters," Journal of Waterway, Port, Coastal and Ocean Engineering , vol. 119, no. 6, pp. 587-605, 1993.

[11] D. L. Kriebel, C. A. Bollmann, "Wave transmission past vertical wave barriers," in Proceeding of the 25th Coastal Engineering Conference, pp. 2470-2480, ASCE, Orlando, Fla, USA, September 1996.

[12] S. Banerjea, M. Kanoria, D. P. Dolai, B. N. Mandal, "Oblique wave scattering by submerged thin wall with gap in finite-depth water," Applied Ocean Research , vol. 18, no. 6, pp. 319-327, 1996.

[13] D. V. Evans, R. Porter, "Efficient calculation of hydrodynamic properties of OWC-type devices," Journal of Offshore Mechanics and Arctic Engineering , vol. 119, no. 4, pp. 210-218, 1997.

[14] H. Martins-Rivas, C. C. Mei, "Wave power extraction from an oscillating water column at the tip of a breakwater," Journal of Fluid Mechanics , vol. 626, pp. 395-414, 2009.

[15] K.-H. Chang, D.-H. Tsaur, L.-H. Huang, "Accurate solution to diffraction around a modified V-shaped breakwater," Coastal Engineering , vol. 68, pp. 56-66, 2012.

[16] E. O. Tuck, "Matching problems involving flow through small holes," Advances in Applied Mechanics , vol. 15, pp. 89-158, 1975.

[17] C. Macaskill, "Reflexion of water waves by a permeable barrier," Journal of Fluid Mechanics , vol. 95, no. 1, pp. 141-157, 1979.

[18] T. Sahoo, "On the scattering of water waves by porous barriers," Zeitschrift für Angewandte Mathematik und Mechanik , vol. 78, no. 5, pp. 364-370, 1998.

[19] M. M. Lee, A. T. Chwang, "Scattering and radiation of water waves by permeable barriers," Physics of Fluids , vol. 12, no. 1, pp. 54-65, 2000.

[20] T. Sahoo, A. T. Chan, A. T. Chwang, "Scattering of oblique surface waves by permeable barriers," Journal of Waterway, Port, Coastal and Ocean Engineering , vol. 126, no. 4, pp. 196-205, 2000.

[21] D. Karmakar, C. Guedes Soares, "Wave transformation due to multiple bottom-standing porous barriers," Ocean Engineering , vol. 80, pp. 50-63, 2014.

[22] R. Gayen, A. Mondal, "A hypersingular integral equation approach to the porous plate problem," Applied Ocean Research , vol. 46, pp. 70-78, 2014.

[23] C. M. Linton, P. McIver Handbook of Mathematical Techniques for Wave/Structure Interactions , Chapman & Hall/CRC, Boca Raton, Fla, USA, 2001.

[24] D. V. Evans, M. A. Peter, "Asymptotic reflection of linear water waves by submerged horizontal porous plates," Journal of Engineering Mathematics , vol. 69, no. 2-3, pp. 135-154, 2011.

[25] X. P. Yu, "Diffraction of water waves by porous breakwaters," Journal of Waterway, Port, Coastal and Ocean Engineering , vol. 121, no. 6, pp. 275-282, 1995.

[26] C. K. Sollitt, R. H. Cross, "Wave transmission through permeable breakwaters," in Proceedings of the 13th Coastal Engineering Conference, pp. 1827-1846, Vancouver, BC, Canada, 1972.

[27] Y. C. Li, Y. Liu, B. Teng, "Porous effect parameter of thin permeable plates," Coastal Engineering Journal , vol. 48, no. 4, pp. 309-336, 2006.

[28] K.-D. Suh, Y. W. Kim, C.-H. Ji, "An empirical formula for friction coefficient of a perforated wall with vertical slits," Coastal Engineering , vol. 58, no. 1, pp. 85-93, 2011.

[29] I. N. Sneddon Mixed Boundary Value Problems in Potential Theory , North Holland, Amsterdam, The Netherlands, 1966.

[30] R. A. Dalrymple, P. A. Martin, "Wave diffraction through offshore breakwaters," Journal of Waterway, Port, Coastal and Ocean Engineering , vol. 116, no. 6, pp. 727-741, 1990.

[31] A. J. Roldan, M. A. Losada, I. J. Losada, "Theoretical study of the wave attenuation in a channel with roughened sides," in Proceedings of the 23rd Coastal Engineering Conference, pp. 2786-2798, Venice, Italy, 1992.

[32] Y. Liu, Y.-C. Li, "Wave interaction with a wave absorbing double curtain-wall breakwater," Ocean Engineering , vol. 38, no. 10, pp. 1237-1245, 2011.

Appendix

Appendix

Matched Eigenfunction Expansion Solution

According to (17) and (18), we define a new function: [figure omitted; refer to PDF] Substituting the expressions for the velocity potentials in (11) with (20) into (A.1), we have [figure omitted; refer to PDF] Multiplying both sides of (A.2) by [figure omitted; refer to PDF] , integrating with respect to [figure omitted; refer to PDF] over [figure omitted; refer to PDF] , and then truncating [figure omitted; refer to PDF] and [figure omitted; refer to PDF] after [figure omitted; refer to PDF] terms, we have [figure omitted; refer to PDF] where [figure omitted; refer to PDF] The unknown coefficients [figure omitted; refer to PDF] are obtained after solving (A.3). Then, the velocity potentials and various hydrodynamic quantities are determined.