2. Problem formulation

Under the assumption of linearized theory, we study a two-dimensional irrotational motion in water due to scattering of a normally incident wave train by a thin porous circular arc shaped barrier submerged in water. A rectangular Cartesian coordinate system is chosen where the y-axis is taken vertically downwards into the fluid and the x-axis is along the undisturbed free surface. A two-dimensional motion is justified due to the reason that motion in every cross section $z=c, c$ being an arbitrary constant, is similar [Reference Mandal and Chakrabarti6]. Here, we study the motion for two cases when the water region is infinitely deep and also when water is of finite depth H. The water occupies the region $-\infty <x < \infty; \ 0<y<\infty $ for deep water (DW) and the region $-\infty <x<\infty; \ 0<y<H$ when the depth of water is finite (FDW). The circular arc shaped barrier C of radius u is placed symmetrically about the y-axis, with centre at $(0,v+u)$ , and the radius through end points makes an angle $\theta $ with the y-axis. A schematic diagram of the problem is given by Figure 1.

For a train of time harmonic surface waves from negative infinity, with circular frequency $\sigma $ represented by $Re\{\phi ^{inc}(x,y) e^{-i\sigma \tau }\}$ , when incident on the porous thin circular arc shaped barrier C, a part of it is reflected, and a part is transmitted above and below the barrier, as shown in Figure 1.

The resulting motion is described by the velocity potential $Re\{\phi (x,y)e^{-i\sigma \tau }\}$ , where $\phi (x,y)$ satisfies the following boundary value problem:

(2.1) $$ \begin{align} &\nabla^2 \phi=0 \quad \text{in the fluid region,} \nonumber\\ &K \phi+ \phi_y=0 \quad \text{on } y=0, \\ &r^{1/2}\nabla \phi \quad \text{is bounded as}\ r \rightarrow 0 , \nonumber \end{align} $$

where r is the distance of any point in the fluid region from either of the sharp ends of the barrier and

$$ \begin{align*} &\nabla \phi \rightarrow \ 0\quad \text{as } y \rightarrow\infty \ \text{for DW,}\nonumber \\ &\frac{\partial \phi}{\partial y}=0 \quad \text{on } y=H \ \text{for FDW.} \end{align*} $$

Figure 1

Schematic diagram of the problem.

The boundary condition on the curved plate surface C is given by

(2.2) $$ \begin{align} \frac{\partial \phi(q+)}{\partial n}=\frac{\partial \phi(q-)}{\partial n}= - i \kappa \lambda(q)\ [\phi](q) \quad \text{on } C, \end{align} $$

where

$$ \begin{align*}\kappa=\begin{cases} K={\sigma^2}/{g} & \text{for DW,}\\ k_0 & \text{for FDW,} \end{cases}\end{align*} $$

and $k_0$ is the unique real positive root of the transcendental equation $k\tanh (kH)=K$ .

Here, $[\phi ](q)=\phi (q+)-\phi (q-) \ (q\in C)$ is the difference of the potential function across the curved barrier C, where

$$ \begin{align*} \begin{aligned} q+ = \{(x,y)\in C\ | \ x^2+ (y-v-u)^2>u^2\}, \\ q- = \{(x,y)\in C\ | \ x^2+ (y-v-u)^2<u^2\}, \end{aligned} \end{align*} $$

${\partial }/{\partial n}$ denotes the normal derivative at a point on C, and $\lambda (q)=\lambda _r(q)+i \lambda _i(q)$ represents the nonuniform porosity parameter [Reference Singh, Gayen and Kundu16] which varies along the arc length of the barrier C. Here, $\lambda _r(q)$ denotes the resistance force coefficient and $\lambda _i(q)$ denotes the inertial force coefficient of the porous barrier, and if $\lambda _r>> \lambda _i,$ then $\lambda $ is taken to be real. Note that the resistance force coefficient $\lambda _r(q)$ resists the passage of water through the pores, while the inertial force coefficient $\lambda _i(q)$ allows the flow of water through the pores [Reference Yu20].

The far-field condition is

(2.3) $$ \begin{align} \phi(x,y)\sim \begin{cases} \phi^{inc}(x,y)+R\phi^{inc}(-x,y)& \text{as}\ x \rightarrow -\infty,\\ T\phi^{inc}(x,y)& \text{as}\ x \rightarrow \infty, \end{cases} \end{align} $$

where R and T are the reflection and transmission coefficients, respectively, which are to be determined. Here, the barrier C is represented parametrically as

$$ \begin{align*} C: \ x=u \sin s\theta; \ y= v+u-u\cos s\theta; \ \ -1 \leq s \leq 1; ~ -\pi<\theta<\pi \end{align*} $$

and

$$ \begin{align*} \phi^{inc}(x,y)= g(y) e^{i\kappa x}, \end{align*} $$

where

$$ \begin{align*} g(y)=\begin{cases} e^{-Ky}& \text{for DW,}\\[4pt] \dfrac{\cosh k_0(H-y)}{\cosh k_0 H} & \text{for FDW.} \end{cases} \end{align*} $$

3. Method of solution

Let $G(x,y;\xi ,\eta )$ be the Green’s function which is the velocity potential due to the presence of a line source at $(\xi ,\eta )$ in the fluid region whose expressions are given below [Reference Mandal and Chakrabarti6].

For DW,

(3.1)

For FDW,

(3.2)

where $r_{\mp }=\sqrt {(x-\xi )^2+(y\mp \eta )^2}$ , and the infinite integrals in (3.1) and (3.2) can be evaluated by standard procedure [Reference Parsons and Martin12].

We apply the Green’s integral theorem to the scattered potential $(\phi -\phi ^{inc})(x,y)$ and the source potential $G(x,y;\xi ,\eta )$ on the region (as shown in Figure 2), bounded externally by the lines $y=0, -X_1\leq x \leq X_1; \ x=-X_1, 0\leq y \leq Y_1; \ y=Y_1, \ -X_1\leq x \leq X_1; \ x=X_1, 0\leq y \leq Y_1$ ; and internally by a very small circle of radius $\delta $ centring $(\xi ,\eta )$ and a contour enclosing curve barrier C. It may be noted that for FDW, $Y_1=H$ and for DW, $Y_1$ will tend to infinity. Making $X_1 \rightarrow \infty ,~ \delta \rightarrow 0$ and $Y_1 \rightarrow \infty \text {for DW and} \ Y_1 \rightarrow ~H \ \text {for FDW, respectively},$ a representation of the potential function in terms of the unknown function $[\phi (p)]$ is obtained as

(3.3) $$ \begin{align} \phi(q)=\phi^{inc}(q)-\frac{1}{2\pi}\int_{C}{[\phi](p)\frac{\partial}{\partial n_p}G(p,q)}\,ds_p. \end{align} $$

Figure 2

Contour for Green’s integral theorem.

In obtaining the above representation, we have used (2.1) and the fact that $\nabla ^2 G= 0$ except at $(\xi , \eta )$ . Also note that the line integrals over all outer boundaries are zero because of the fact that both $(\phi -\phi ^{inc})(x,y)$ and G satisfy similar boundary conditions on all the outer boundaries. The only contribution comes from the contour C enclosing the barrier and the small circle of radius $\delta $ enclosing the point $(\xi , \eta )$ (Figure 2 and [Reference Mandal and Chakrabarti6]). In (3.3), $q=(\xi ,\eta )$ is a field point, $p\equiv (x,y)$ is a point on barrier C, $[\phi ](p)$ is the discontinuity of potential function $\phi (x,y)$ across C and ${\partial }/{\partial n_p}$ is the normal derivative at the point p on C. It may be noted that the unknown function $[\phi ](p)$ vanishes at the end points of the arc-shaped barrier [Reference Mandal and Chakrabarti6].

Using the boundary condition (2.2) on the barrier in (3.3), we obtain the following second kind of hypersingular integral equation:

(3.4) $$ \begin{align} \frac{1}{2\pi}\int_C\hspace{-0.4250cm}\times~{[\phi](p)\frac{\partial^2}{\partial n_q \partial n_p}{ G}(p,q)}\,ds_p \ - \ i \kappa \ \lambda(q)\ [\phi](q)\ =\frac{\partial \phi^{inc}}{\partial n_q}(q),\quad q\in C , \end{align} $$

where ${\partial }/{\partial n_q}$ is the normal derivative at the point q on C. The parametric equation of C at the points $p=(x,y)$ and $q=(\xi , \eta )$ in (3.4) is given by

(3.5) $$ \begin{align}\begin{aligned} x(s)&=u\sin s\theta,\quad y(s)=v+u-u\cos s\theta ;\quad -1\leq s\leq 1, \\ \xi(t)&=u\sin t\theta,\quad \eta(t)=v+u-u\cos t\theta;\quad -1\leq t\leq 1. \end{aligned}\end{align} $$

Thus, with the help of the parametrization given by (3.5) along with the introduction of a new unknown function $f(s)=[\phi ](p(s)),$ the integral equation (3.4) can be rewritten as

(3.6) $$ \begin{align} \int_{-1}^1 f(s)\bigg[\frac{1}{(s-t)^2} + u^2\theta^2 \mathcal{L}(s,t)\bigg]\,ds +2\pi i u \theta \kappa \lambda(t) f(t) =2\pi u \theta g_1(t) , \end{align} $$

where

$$ \begin{align*} g_1(t)=\begin{cases} Ke^{-K\eta(t)+i(K\xi(t)+t\theta)} & \text{for DW},\\ \dfrac{k_0e^{i k_0\xi(t)}\sinh{(k_0(H-\eta(t))+\imath t\theta)}}{\cosh{k_0 H}} & \text{for FDW}. \end{cases} \end{align*} $$

The expression of $\mathcal {L}(s,t)$ is given as follows.

For DW,

$$ \begin{align*} \mathcal{L}(s,t)&=\bigg[\frac{1}{4u^2 \sin^2({(s-t) \theta}/2)}-\frac{1}{u^2 \theta^2 (s-t)^2}\bigg] \nonumber\\ &\quad +\cos (s-t)\theta \bigg[\frac{Y^2-X^2}{(X^2+Y^2)^2}+\frac{2KY}{X^2+Y^2}+2K^2 \Phi_0(X,Y) \bigg]\nonumber \\ &\quad +2\sin (s-t)\theta\bigg[K \frac{\partial \Phi_0}{\partial X}-\frac{XY}{(X^2+Y^2)^2}\bigg] \end{align*} $$

with

$$ \begin{align*}X=x-\xi=u(\sin{s \theta}-\sin{t \theta}); \quad Y=y+\eta=2v+2u-u(\cos{s \theta}+\cos{}t \theta), \end{align*} $$

$$ \begin{align*}\Phi_0(X,Y)= \imath \pi e^{-KY+\imath K|X|} +\int_0^\infty{\frac{e^{-k|X|}}{k^2+K^2}\{k\cos{kY}-K\sin{kY}\}\,dk}.\end{align*} $$

For FDW,

$$ \begin{align*} \mathcal{L}&(s,t)=\bigg[\frac{1}{4u^2 \sin^2({(s-t) \theta}/{2})}-\frac{1}{u^2 \theta^2 (s-t)^2}\bigg] \nonumber \\&\ \qquad +\frac{X^2-Y^2}{(X^2+Y^2)^2}\ \cos(s-t)\theta-\frac{2XY}{(X^2+Y^2)^2}\ \sin(t-s)\theta\nonumber \\&\ \qquad +2\sin{t\theta}\sin{s\theta}\int_0^{\infty}\mu(k)\sinh{k\eta}\sinh{ky}\cos{kX}\,dk \nonumber \\&\ \qquad -2\sin{t\theta}\cos{s\theta}\int_0^{\infty}\mu(k)\sinh{k\eta}\cosh{ky}\sin{kX}\,dk \nonumber \\&\ \qquad +2\cos{t\theta}\sin{s\theta}\int_0^{\infty}\mu(k)\cosh{k\eta}\sinh{ky}\sin{kX}\,dk \nonumber \\&\ \qquad +2\cos{t\theta}\cos{s\theta}\int_0^{\infty}\mu(k)\cosh{k\eta}\cosh{ky}\cos{kX}\,dk \nonumber\\&\ \qquad+\sin t\theta \sin s\theta \bigg[\zeta i \cosh k_0(H-y)\cosh k_0(H-\eta)\ e^{ik_0|X|}\nonumber\\&\ \qquad-\sum_{n=1}^{\infty}\omega_n\cos k_n(H-y)\cos k_n(H-\eta)\ e^{-k_n|X|}+\frac{2\pi}{H}\sum_{n=0}^{\infty}\alpha_n\sin \alpha_ny\sin \alpha_n\eta\ e^{-\alpha_n|X|}\bigg] \nonumber \\&\ \qquad+\cos t\theta \cos s\theta \bigg[\zeta i \sinh k_0(H-y)\sinh k_0(H-\eta)\ e^{ik_0|X|} \nonumber\\&\ \qquad+\sum_{n=1}^{\infty}\omega_n\sin k_n(H-y)\sin k_n(H-\eta)\ e^{-k_n|X|}-\frac{2\pi}{H}\sum_{n=0}^{\infty}\alpha_n\cos \alpha_ny\cos \alpha_n\eta\ e^{-\alpha_n|X|}\bigg] \nonumber \\&\ \qquad+\sin t\theta \cos s\theta\bigg[\zeta \sinh k_0(H-y)\cosh k_0(H-\eta)\ e^{ik_0|X|} \nonumber \\&\ \qquad-\sum_{n=1}^{\infty}\omega_n\sin k_n(H-y)\cos k_n(H-\eta)\ e^{-k_n|X|}+\frac{2\pi}{H}\sum_{n=0}^{\infty}\alpha_n\cos \alpha_ny\sin \alpha_n\eta\ e^{-\alpha_n|X|}\bigg]\nonumber \\&\ \qquad-\cos t\theta \sin s\theta \bigg[\zeta \cosh k_0(H-y)\sinh k_0(H-\eta)\ e^{ik_0|X|} \nonumber \\&\ \qquad-\sum_{n=1}^{\infty}\omega_n\cos k_n(H-y)\sin k_n(H-\eta)\ e^{-k_n|X|}+\frac{2\pi}{H}\sum_{n=0}^{\infty}\alpha_n\sin \alpha_ny\cos \alpha_n\eta\ e^{-\alpha_n|X|}\bigg], \end{align*} $$

where $k_0$ and $ik_n$ are roots of the transcendental equations $k\tanh kH-K=0$ ; $i\alpha _n$ terms are roots of the equation $\cosh kH=0$ and

$$ \begin{align*}\mu(k)=\frac{k\ e^{-kH}}{\cosh kH}, \quad \zeta=\frac{4\pi k_0^2}{2k_0 H+\sinh(2k_0H)}, \quad \omega_n=\frac{4\pi k_n^2}{2k_nH+\sin(2k_nH)}.\end{align*} $$

It may be noted here that if $\lambda (t)=0$ , then the hypersingular integral equation of second kind (3.6) reduces to a hypersingular integral equation of first kind. Also for $\lambda (t)=0$ , the condition (2.2) on the porous barrier reduces to the condition on a rigid barrier. This shows that wave interaction with porous barrier produces the second kind of hypersingular integral equation, while that with a rigid barrier produces the first kind of hypersingular integral equation. It is also interesting to observe that the kernel of the hypersingular integral equation for water of finite depth is quite complicated in comparison to the kernel in the case of deep water.

The unknown physical quantities R and T are evaluated by making $\xi \rightarrow \mp \infty $ in (3.3) and comparing with the far-field condition for $\phi (\xi ,\eta )$ as given by (2.3). This yields

(3.7) $$ \begin{align} R=\begin{cases}\displaystyle -iKu\theta\int_{-1}^{1} f(s) e^{-K y(s)+iKx(s)s\theta}\, ds & \text{for DW,}\\[.2cm] \displaystyle-\dfrac{2u\theta i k_0 \cosh k_0H}{2k_0H+\sinh 2k_0H}\int_{-1}^{1} f(s) \sinh (k_0(H-y(s))+is\theta)\ e^{ik_0x(s)}\, ds &\text{for FDW,} \end{cases} \end{align} $$

and

(3.8) $$ \begin{align} T=\begin{cases}\displaystyle 1-iKu\theta\int_{-1}^{1} f(s) e^{-K y(s)-iKx(s)s\theta}\, ds & \text{for DW,}\\[.2cm] \displaystyle 1-\dfrac{2u\theta i k_0 \cosh k_0H}{2k_0H+\sinh 2k_0H}\int_{-1}^{1} f(s) \sinh (k_0(H-y(s))-is\theta)\ e^{-ik_0x(s)}\, ds & \text{for FDW.} \end{cases} \end{align} $$

Thus, from (3.7) and (3.8), we observe that R and T are expressed in terms of the unknown function $f(s)$ which can be obtained by solving the integral equation (3.6).

3.1. Solution of the integral equation

This section illustrates the boundary element method and the collocation method of the solution of the second kind of hypersingular integral equation (3.6). We may mention here that the numerical solution of the hypersingular integral equation is limited in the literature due to the reason that hypersingular integrals are not amenable to the numerical methods. Application of the collocation method using Chebychev’s polynomial in solving the hypersingular integral equation was initiated by Parsons and Martin [Reference Parsons and Martin12]. Later, Samanta et al. [Reference Samanta, Chakraborty and Banerjea13] applied the boundary element method to solve the hypersingular integral equation. The boundary element method is a simple method computationally, but the convergence of the collocation method is faster than the boundary element method. The present analysis demonstrates the application of two methods.

Boundary element method. We will solve the integral equation (3.6) numerically using the boundary element method (BEM). The method of solution of the hypersingular integral equation using BEM and the convergence study is described in detail in [Reference Samanta, Chakraborty and Banerjea13].

Since $[\phi ]=0$ at the ends of the barrier, we have $f(s)=0$ at $s=\pm 1$ . This behaviour asserts that we may assume

(3.9) $$ \begin{align} f(s)=\sqrt{1-s^2}\psi(s) , \end{align} $$

where $\psi (s)$ is a regular function in $[-1,1]$ . Noting (3.9), the hypersingular equation (3.6) takes the form

(3.10) $$ \begin{align} &\int_{-1}^1\sqrt{1-s^2}\bigg[\frac{1}{(s-t)^2}+u^2\theta^2\mathcal{L}(s,t)\bigg]~\psi(s)\,ds +2\pi i u \theta \kappa \lambda(t) \sqrt{1-t^2}\psi(t) \nonumber \\ &\quad=2\pi u \theta g_1(t), \quad -1\leq t \leq 1. \end{align} $$

As per the BEM, the range of integration $[-1,1]$ is divided into m number of equal line elements, that is, $[-1,1]=\bigcup _{j=1}^{m}[a_{j-1},a_j]$ , where $a_{j-1}$ and $a_j$ are the end points of the jth line element, and $a_0=-1$ , $a_m=1$ and $a_j=a_0+jr'$ , $r'=({a_m-a_0})/{m}$ . Writing $s=s_j$ where $s_j\in [a_{j-1},a_j],$ and $t=t_i$ where $t_i\in [a_{i-1},a_i], i,j=1,2,\ldots ,m$ ,

$$ \begin{align*} s_j=(1-\tau)a_{j-1}+\tau a_j, \quad 0\leq \tau \leq 1; \quad t=t_i=(1-\gamma)a_{i-1}+\gamma a_i, \quad 0\leq \gamma\leq 1. \end{align*} $$

Thus, (3.10) takes the form

(3.11) $$ \begin{align} &\sum_{j=1}^{m}\int_{0}^1 \sqrt{1-s_j^2}\bigg{\{}\frac{1}{(s_j-t_i)^2}+u^2 \theta^2 \mathcal{L}(s_j,t_i)\bigg{\}}\ \psi(s_j)\ r' \ d \tau+2\pi i u \theta \kappa \lambda(t_i) \sqrt{1-t_i^2}\psi(t_i)\nonumber\\ &\quad=2\pi u \theta g_1(t_i)\ , \quad i=1,2,\ldots ,m.\\ \end{align} $$

Assuming that the unknown function $\psi (s_j)=\psi _j$ is a constant for the jth line element, $j=1,2,\ldots ,m$ , the integral equation (3.11) is reduced to the following system of linear algebraic equations [Reference Samanta, Mondal and Banerjea14]:

(3.12) $$ \begin{align} \sum_{j=1}^{m}\ a_{ij}\ \psi_j=2\pi u \theta g_{1i}, \quad i=1,2,\ldots ,m, \end{align} $$

where

$$ \begin{align*} a_{ij}&=\int_{0}^1 \sqrt{1-s_j^2}\bigg{\{}\frac{1}{(s_j-t_i)^2}+u^2 \theta^2 \mathcal{L}(s_j,t_i)\bigg{\}}\ r' \ d \tau \ +\delta_{ij}2\pi i u \theta \kappa \lambda(t_i) \sqrt{1-t_i^2} , \nonumber\\ &\qquad \qquad \quad i=1,2,\ldots,m,\quad j=1,2,\ldots,m, \\ g_{1i}&=g_1(t_i), \quad\ i=1,2,\ldots,m. \end{align*} $$

Now solving the system of linear equations (3.12), we obtain the unknown coefficients $\psi _j,\ j=1,2,\ldots ,m$ , and then we approximate the solution of $f(s)$ in each line interval to evaluate the value of R and T from (3.7) and (3.8).

Collocation method. In this method, the unknown function $f(s)$ is approximated as [Reference Parsons and Martin12]

(3.13) $$ \begin{align} f(s)\approx \sqrt{1-s^2}\sum_{n=0}^N{d_n U_n(s)} , \end{align} $$

where N is an integer, $U_n(s)$ is the nth-order Chebyshev polynomial of second kind, and $d_n$ are unknown complex constants for each $n=0,1,2,\ldots ,N$ .

Substituting $f(s)$ from (3.13) into hypersingular integral equation (3.6), we obtain the following system of algebraic equations [Reference Parsons and Martin12]:

(3.14) $$ \begin{align} \sum_{n=0}^{N}\ d_n\ B_n(t)\ = 2\pi u \theta g_1(t) ,\quad -1\leq t \leq 1 , \end{align} $$

where

$$ \begin{align*} B_n(t)&= -\ \pi(n+1)U_n(t) \ +\ u^2\theta^2\int_{-1}^{1}\sqrt{1-s^2}\ U_n(s)\ \mathcal{L}(s,t)\, ds \nonumber \\ &\quad + 2\pi i u \theta \kappa \lambda(t) \sqrt{1-t^2}\ U_n(t). \end{align*} $$

The collocation points $t=t_j$ can be taken as

$$ \begin{align*} t_j=\cos \bigg(\frac{2j+1}{2N+2}\bigg) \pi , \quad j=0,1,2,\ldots,N , \end{align*} $$

so that (3.14) yields the following system of linear equations:

(3.15) $$ \begin{align} \sum_{n=0}^{N}\ d_n\ B_n(t_j)\ = 2\pi u \theta \ g_1(t_j),\quad j=0,1,2,\ldots,N , \end{align} $$

where $d_n$ terms are unknown constants to be determined. Now we can approximate the solution $f(s)$ of the integral equation (3.6) by substituting $d_n$ terms in (3.13), and then using it in (3.7) and (3.8), we can evaluate the value of reflection coefficient R and transmission coefficient T, respectively.

3.2. Energy identity

The energy dissipation due to porosity yields $|R|^2+|T|^2<1$ , and this can be mathematically justified by applying Green’s integral theorem to the functions $\phi $ and $\bar {\phi }$ in the suitable chosen region which results in the energy identity as

$$ \begin{align*}|R|^2+|T|^2=1-J ,\end{align*} $$

where

(3.16) $$ \begin{align} J=\begin{cases} \displaystyle 2K u\theta \ \int_{-1}^1 |f(t)|^2 \ \lambda_r (t) \, dt & \text{for DW,}\\[.2cm] \displaystyle\dfrac{2K u\theta \cosh k_0H}{k_0H+(1/2)\sinh 2k_0H}\ \int_{-1}^1 |f(t)|^2 \ \lambda_r (t) \, dt & \text{for FDW} \end{cases} \end{align} $$

is the energy dissipation coefficient.

4. Numerical results

In this section, the numerical estimates for reflection coefficients $|R|$ and energy dissipation coefficients J are studied graphically for deep water as well as finite depth of water. For DW, $|R|$ and J are depicted graphically against the dimensionless wave number $Ku$ for different values of nondimensional parameters $v/u, \theta , \lambda (t)$ . When the water region is of finite depth H (FDW), $|R|$ and J are plotted graphically against the wave number $KH$ for different values of nondimensional parameters $v/H,~ u/H ,\theta , \lambda (t)$ .

The porosity parameter $\lambda (t)=\lambda _r(t)+i \lambda _i(t), -1 \leq t \leq 1$ , is chosen here as:

1. (i) $\lambda (t)= (1+t)/2;$

2. (ii) $(1+i/2)(1+t)/2;$

3. (iii) $\lambda (t)= t^2;$

4. (iv) $t^2 (1+i/2),\quad 1 \leq t \leq 1. $

We may mention here that when the inertial force coefficient is negligible as compared with the resistance force coefficient, then $\lambda (t)$ is taken to be real. The distribution of pores in each of the cases is described below.

(a) When $\lambda (t)= (1+t)/2,$ then $\lambda _i (t)=0$ and $\lambda _r(t)=(1+t)/2$ , so that $\lambda (t)=\lambda _r(t)$ is real. In this case, the distribution of pores in the barrier is such that $\lambda (t)=0$ for $t=-1$ , that is, at one end of the barrier, and increases as t increases till $\lambda _r(t)$ reaches maximum value $\lambda _r(t)=1$ at $t=1$ (Figure 3(a)).

Figure 3

Porosity distribution of curved barrier (a) $\lambda (t)= (1+t)/2$ , (b) $\lambda (t)= t^2$ .

Table 1

Reflection coefficient for $\theta = {\pi }/{10}$ , ${v}/{u}=0.1$ , $\lambda (t)=t^2$ .

(b) When $\lambda (t)=(1+i/2)(1+t)/2,$ then $\lambda _r(t)=(1+t)/2$ and $\lambda _i(t)=(1+t)/4$ . In this case, at $t=-1$ , $\lambda _r(t)=\lambda _i(t)=0$ , that is, at one end of the plate, porosity is zero. As t increases, $\lambda _r(t)$ and $\lambda _i(t)$ increase and reach maximum values $\lambda _r(t)=1$ and $\lambda _i(t)=1/2,$ respectively. So $\lambda (t)=(1+i/2)(1+t)/2$ suggests that the distribution of pores in the barrier is such that the barrier is rigid at one end, and the porosity increases along the barrier in the sense that the resistance force coefficient $\lambda _r(t)$ and inertial force coefficient $\lambda _i(t)$ increase linearly along the barrier, until at the other end, the resistance force coefficient and the inertial force coefficient reach a maximum value. Similarly, the distribution of pores in the barrier for $\lambda (t)= t^2, t^2 (1+i/2)$ is such that $\lambda (t)=0$ for $t=0$ , that is, the barrier is rigid at the centre. The resistance force coefficient $\lambda _r(t)$ and inertial force coefficient $\lambda _i(t)$ increase towards the two ends of the barrier until they reach maximum values at $t=\pm 1$ , that is, at the two ends of the barrier (Figure 3(b)).

The reflection coefficient $|R|$ and the amount of wave energy dissipated J can be computed numerically from (3.7) and (3.16), respectively, by solving the integral equation (3.6). For numerical computation, the value of m (number of elements) in (3.12) is chosen as $m=45$ and the value of N (number of collocation points) in (3.15) is taken as 15. We may say here that the convergence of the collocation method is faster than BEM, but computation execution of BEM is simpler than the collocation method.

Tables 1 and 2 show a comparison of the values of $|R|$ obtained from BEM and the collocation method for DW with

$$ \begin{align*}Ku = \{ 0.5, 1.0, 1.5, 2.0 \},\quad \theta = \pi /10, \quad v/u=0.1, \quad \text{for } \lambda(t)= \{ t^2, t^2(1+i/2) \}.\end{align*} $$

Similar comparisons of $|R|$ are presented in Tables 3 and 4 for FDW with

$$ \begin{align*}KH &= \{ 0.5, 1.0, 1.5, 2.0 \}, \theta =\pi /10, ~v/H=0.3, ~u/H=0.5,\\ &\qquad\text{for } \lambda(t) = \{ t^2, t^2(1+i/2) \}.\end{align*} $$

Table 2

Reflection coefficient for $\theta = {\pi }/{10}$ , ${v}/{u}=0.1$ , $\lambda (t)=t^2(1+{i}/{2})$ .

Table 3

Reflection coefficient for $\theta = {\pi }/{10}$ , ${v}/{H}=0.3$ , ${u}/{H}=0.5$ , $\lambda (t)=t^2$ .

Table 4

Reflection coefficient for $\theta = {\pi }/{10}$ , ${v}/{H}=0.3$ , ${u}/{H}=0.5$ , $\lambda (t)=t^2(1+{i}/{2})$ .

It is observed from Tables 1, 2, 3 and 4 that values of $|R|$ agree with each other up to 3 or 4 decimal places. The accuracy can be improved by increasing the number of line elements in the boundary element method and the number of collocation points N in the collocation method. Here, note that the energy identity, $|R|^2+|T|^2=1$ , has been verified for the rigid curve barrier, that is, $\lambda =0$ . Also, the energy identity $|R|^2+|T|^2+J=1$ for a porous curve barrier has been verified.

The reflection coefficient $|R|$ for a rigid barrier, that is, for $\lambda =0$ in deep water, is compared with the results of Parsons and Martin [Reference Parsons and Martin12] in Figure 4(a) and in Figure 4(b), $|R|$ for rigid barrier, that is, for $\lambda =0$ in finite depth water are compared with [Reference Liu and Li5, Table 1]. It is observed that both results are in good agreement with each other. In Figure 5, $|R|$ for deep water is compared with $|R|$ for water of finite depth, for large depth, that is, $H/u=3,~v/u=0.1, ~\theta =\pi /10$ . A good matching of $|R|$ is observed for both infinite and finite depth of water region.

Figure 4

(a) Comparison between our result and the results of Parsons and Martin [Reference Parsons and Martin12] for $v/u=0.1, \theta =\pi /2$ . (b) Comparison between our result and the results of Liu and Li [Reference Liu and Li5] for $v/H=0.5, u/H=0.5$ and $\theta =\pi /2$ .

Figure 5

$|R|$ against $Ku$ for deep water and for water of finite depth where $v/u=0.1$ , $\theta =\pi /10$ , $H/u=3$ .

Behaviour of  $|R|$  and J for deep water. In Figures 6(a)–6(c) and 7(a)–7(c), the reflection coefficient $|R|$ and the wave energy dissipation coefficient J are computed numerically for a barrier in deep water (DW) and depicted graphically against the wave number $Ku$ , for $v/u=0.1$ , the porosity parameter $\lambda (t) =\{1,1+i/2,(1+t)/2, (1+i/2)(1+t)/2 ,t^2, t^2(1+i/2) \}$ and $\theta = \{\pi /10, 3\pi /10, \pi /2\}.$ In Figures 8(a)–8(c) and 9(a)–9(c), the reflection coefficient $|R|$ and the wave energy dissipation coefficient J are computed numerically for a barrier in water of finite depth (FDW), and depicted graphically against the wave number $KH$ , for $v/H=0.3, u/H=0.5$ , the porosity parameter $\lambda (t) =\{1,1+i/2,(1+t)/2, (1+i/2)(1+t)/2 ,t^2, t^2(1+i/2) \}$ and $\theta = \{\pi /10, 3\pi /10, \pi /2\}.$

Figure 6

$|R|$ against $Ku$ for $v/u=0.1$ , (a) $\theta =\pi /10$ , (b) $\theta =3\pi /10$ , (c) $\theta =\pi /2$ .

Figure 7

J against $Ku$ for $v/u=0.1$ : (a) $\theta =\pi /10$ ; (b) $\theta =3\pi /10$ ; (c) $\theta =\pi /2$ .

Figure 8

$|R|$ against $KH$ for $v/H=0.3$ , $u/H=0.5$ , for (a) $\theta =\pi /10$ , (b) $\theta =3\pi /10$ , (c) $\theta =\pi /2$ .

Figure 9

J against $KH$ for $v/H=0.3$ , $u/H=0.5$ , for (a) $\theta =\pi /10$ , (b) $\theta =3\pi /10$ , (c) $\theta =\pi /2$ .

In Figures 6(a)–6(c), the reflection coefficient $|R|$ is plotted against wave numbers $Ku$ for DW, for $v/u=0.1$ , the porosity parameter $\lambda (t) =\{1,1+i/2,(1+t)/2, (1+i/2)(1+t)/2 ,t^2, t^2(1+i/2) \}$ and $\theta = \{\pi /10, 3\pi /10, \pi /2\}.$

It is observed from Figures 6(a)–6(c) that for each $\theta = \{\pi /10, 3\pi /10, \pi /2\}$ ,

$$ \begin{align*}|R|_{\lambda(t)=t^2}>|R|_{\lambda(t)=(1+t)/2}>|R|_{\lambda(t)=1}, \quad |R|_{\lambda(t)=t^2(1+i/2)}>|R|_{\lambda(t)=(1+t)(1+i/2)/2}>|R|_{\lambda(t)=1+i/2}\end{align*} $$

for $-1<t<1 $ in general. This shows that a barrier with these particular choices of nonuniform porosity induces more reflection than that of a barrier with constant porosity.

In Figures 7(a)–7(c), the amount of wave energy dissipated J is plotted against wave numbers $Ku$ for DW, for $v/u=0.1$ , the porosity parameter $\lambda (t) =\{1,1+i/2,(1+t)/2, (1+i/2)(1+t)/2 ,t^2, t^2(1+i/2) \}$ and $\theta = \{\pi /10, 3\pi /10, \pi /2\}.$ From Figures 7(a)–7(c), the following observations are made.

(i) For a short barrier, $\theta =\pi /10,~~ J_{\lambda (t) =(1+t)/2}> J_{\lambda (t) =1}$ for $Ku>1.4$ . This shows that a nonuniform porous barrier whose porosity increases circumferentially from one end to the other end dissipates energy of short crested surface waves more than for a barrier with a uniform porosity distribution. This confirms the efficacy of a nonuniform porous barrier with this particular porosity distribution over a uniform porous barrier as a model for breakwater. As the length of the barrier increases, the energy dissipation for short-crested waves is more for a barrier with uniform porosity than for a barrier with nonuniform porosity.

(ii) For any length of the barrier, the energy dissipation coefficient J for a barrier which is rigid at the centre and porous towards the end is least.

(iii) $J_{\lambda _{\text {real}}}>J_{\lambda _{\text {complex}}}$ showing that the inertial force coefficient of the porous barrier allows the passage of water through the pores of the barrier and thereby prevents dissipation of wave energy.

Behaviour of $|R|$ and J for water of finite depth. In Figures 8(a)–8(c), the reflection coefficient $|R|$ for FDW is plotted against the wave number $KH$ , for $v/H=0.3, u/H=0.5$ , the porosity parameter $\lambda (t) =\{1,1+i/2,(1+t)/2, (1+i/2)(1+t)/2 ,t^2, t^2(1+i/2) \}$ and $\theta = \{\pi /10, 3\pi /10, \pi /2\}.$

From these figures, we make the following observations.

(1) From Figures 8(a)–8(c), it is observed in general that $|R|_{\lambda _r}>|R|_{\lambda _r+i \lambda _i}$ for all values of $\lambda (t) =\{1,1+i/2,(1+t)/2, (1+i/2)(1+t)/2 ,t^2, t^2(1+i/2) \}$ , which shows that the presence of an inertial force coefficient in the porous barrier allows the passage of water through it and thereby reduces the reflection.

(2) Figure 8(a) shows that for $KH>1$ ,

$$ \begin{align*} |R|_{\lambda(t)=t^2}> |R|_{\lambda(t)=(1+t)/2}>|R|_{\lambda(t)=1}, \quad |R|_{\lambda(t)=t^2(1+i/2)}> |R|_{\lambda(t)=(1+i/2)(1+t)/2}>|R|_{\lambda(t)=(1+i/2)}, \end{align*} $$

which infers that for a small arc length of the barrier ( $\theta =\pi /10,~ {3\pi }/10$ ), the reflection of waves with moderate to short wavelength is more for the barrier which is rigid at the centre and porous towards the ends ( $\lambda (t)= t^2,~ t^2 (1+i/2)$ ) than for the barrier which is rigid at one end and becomes porous towards the other end ( $\lambda (t)= (1+t)/2, (1+i/2)(1+t)/2$ ). Moreover, the reflection coefficient for these waves for a barrier with constant porosity is less than the barrier with variable porosity. Also for $KH<1$ , that is, the long waves which are towards the bottom, they do not feel the presence of the barrier. Almost similar behaviour is observed from Figure 8(b) when $\theta = 3\pi /10.$ It may be mentioned that for $KH>2.2$ , $|R|_{\lambda (t)=1}>|R|_{\lambda (t)=(1+t)/2}.$

However, when the length of the barrier increases to $\pi /2$ , the change in the behaviour of $|R|$ can be observed from Figure 8(c). For $KH<0.3,~ |R|$ coincides for all values of $\lambda (t)$ showing that the long waves do not feel the presence of the barrier. For $0.3<KH<1.5,~ |R|_{\lambda (t)=(1+t)/2}>|R|_{\lambda (t)=1}>|R|_{\lambda (t)=t^2}$ . For $1.5<KH<2.3, |R|_{\lambda (t)=(1+t)/2}>|R|_{\lambda (t)=t^2}>|R|_{\lambda (t)=1}$ and for $KH>2.3, |R|_{\lambda (t)=t^2}>|R|_{\lambda (t)=(1+t)/2}>|R|_{\lambda (t)=1}.$ Almost similar behaviour in $|R|$ is observed when the inertial force coefficient is present in the porous material of the barrier.

Thus, for these choices of porosity distributions of the barrier, the behaviour of $|R|$ in the case of finite depth of water region (FDW) is different from the behaviour of $|R|$ for deep water (DW), which may be attributed to the interaction of waves with the porous barrier and the bottom of the water region.

Figures 9(a)–9(c) capture the behaviour of wave energy dissipation coefficient J for FDW for various values of the wave number $KH$ , depth of the barrier below the mean free surface, that is, $v/H=0.3,$ radius of the barrier, that is, $u/H=0.5$ , and the porosity parameter $\lambda (t) = \{1,1+i/2,(1+t)/2, (1+i/2)(1+t)/2 ,t^2, t^2(1+i/2) \}$ and $\theta = \{\pi /10, 3\pi /10, \pi /2\}.$

(1) It is observed from these figures that

$$ \begin{align*}&J_{\lambda(t)=t^2}<J_{\lambda(t)=(1+t)/2}<J_{\lambda(t)=1}, \\ &J_{\lambda(t)=t^2(1+i/2)}<J_{\lambda(t)=(1+t)(1+i/2)/2}\leq J_{\lambda(t)=1+i/2}, \, -1<t<1 \end{align*} $$

for $\theta = \{\pi /10, 3\pi /10, \pi /2\} $ . This shows that energy dissipation for a barrier with constant porosity is more than that of a barrier with variable porosity. Also the energy dissipation is least for $\lambda (t)= t^2,~ t^2(1+i/2)$ .

(2) From these figures, it is seen that for a particular value of $KH$ , J increases as $\theta $ increases. This shows that the dissipation of wave energy increases with length of the curved barrier.

(3) Also from these figures, $J_{\lambda _r}>J_{\lambda _r+i \lambda _i}$ for any length of the barrier. So the inertial force coefficient of the porous barrier increases transmission and prevents dissipation of the wave energy.

Figure 10

(a) $|R|$ and J against $Ku$ for $v/u=0.1$ , $\theta =\pi /5$ for deep water, (b) $|R|$ and J against $Ku$ for $v/u=0.1$ , $H/u=1.2$ $\theta =\pi /5$ for water of finite depth.

Comparison of uniform and nonuniform porosity distribution. Figure 10(a) depicts $|R|$ and J against $Ku$ for $v/u=0.1$ , $\theta =\pi /5$ for deep water while Figure 10(b) depicts $|R|$ and J against $Ku$ for $v/u=0.1$ , $H/u=1.2$ $\theta =\pi /5$ for water of finite depth for uniform porosity distribution given by $\lambda (t)=1$ and for nonuniform porosity distribution $\lambda (t)=(0.5+t^2)$ . This particular choice of nonuniform porosity distribution shows that the porosity is less near the middle of the barrier and increases towards the ends of the barrier. From both figures, it is seen that energy dissipation J for nonuniform porosity is more than that of uniform porosity. In particular, the J for a barrier with nonuniform porosity is significantly more than J for a barrier with uniform porosity when the water is of finite depth. This shows that a barrier with porosity distribution $\lambda (t)=(0.5+t^2)$ reduces the wave power by dissipating the energy, and thus serves as a better model for breakwater than a barrier with uniform porosity. This phenomenon is also observed in Figure 7(a).