Table of Contents
ISRN Applied Mathematics
VolumeΒ 2011Β (2011), Article IDΒ 341564, 16 pages
Research Article

Solution of Singular Integral Equations Involving Logarithmically Singular Kernel with an Application in a Water Wave Problem

1Department of Mathematics, Jadavpur University, Kolkata 700032, India
2Department of Mathematics, Netajinagar Day College, Regent Estate, NSC Bose Road, Kolkata 700023, India
3Department of Mathematics, Indian Institute of Science, Bangalore 560012, India

Received 10 March 2011; Accepted 1 April 2011

Academic Editors: A.Β Bellouquid and T. Y.Β Kam

Copyright Β© 2011 Sudeshna Banerjea et al. 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.


A direct function theoretic method is employed to solve certain weakly singular integral equations arising in the study of scattering of surface water waves by vertical barriers with gaps. Such integral equations possess logarithmically singular kernel, and a direct function theoretic method is shown to produce their solutions involving singular integrals of similar types instead of the stronger Cauchy-type singular integrals used by previous workers. Two specific ranges of integration are examined in detail, which involve the following: Case(i) two disjoint finite intervals (0,π‘Ž)βˆͺ(𝑏,𝑐) and (π‘Ž,𝑏,𝑐beingfinite) and Case(ii) a finite union of 𝑛 disjoint intervals. The connection of such integral equations for Case(i), with a particular water wave scattering problem, is explained clearly, and the important quantities of practical interest (the reflection and transmission coefficients) are determined numerically by using the solution of the associated weakly singular integral equation.

1. Introduction

The singular integral equationξ€œπ΅|||πœ™(𝑑)lnπ‘‘βˆ’π‘₯|||𝑑+π‘₯𝑑𝑑=𝑓(π‘₯),π‘₯∈𝐡,(1.1) where 𝐡 is a single or a finite union of finite intervals is encountered in many branches of mathematical physics (cf. Lewin [1], Ursell [2], and Mandal and Kundu [3]).This type of integral equation is usually solved in the literature by converting the weak singularity (logarithmic type) in the kernel to a Cauchy-type singularity which is strong in nature. For an integral equation with weakly singular kernel, the associated integral is defined in the standard Riemann sense, while for an integral equation with strongly singular kernel, the singular integral has to be defined appropriately so as to give a consistent mathematical sense. Consequently, the integrals with weak singularities are more amenable to the numerical techniques as compared to the integrals with strong singularity which requires artificial methods for its numerical evaluation. Noting this difficulty with the integrals with strong singularities, Chakrabarti and Manam [4], Manam [5], Banerjea and Dutta [6], and Banerjea and Sarkar [7] used function theoretic method to solve (1.1) for various 𝐡 without converting it to Cauchy-type integral equation. However, the solution function in [4–7] involves Cauchy-type singularity.

Recently, Chakrabarti [8] has developed a direct function theoretic method for (1.1) with 𝐡=[π‘Ž,𝑏] to obtain its solution. The novelty of the method is that it avoids the occurrence of strong singularity such as Cauchy-type singularity in the solution function of the integral equation.

The present paper is a generalization of the method in [8] with a little modification. We have applied the method of [8] to the integral equation (1.1), where consists of (i) two disjoint finite intervals and (ii) a finite union of 𝑛 disjoint intervals βˆͺ𝑛𝑖=1(π‘Žπ‘–,𝑏𝑖). As an application, we have applied the solution of (1.1) where 𝐡 consists of two disjoint finite intervals (0,π‘Ž)βˆͺ(𝑏,𝑐) to the problem of water wave scattering by a partially immersed vertical barrier with a gap present in deep water. The important quantities of physical interest, that is, the reflection and transmission coefficients, are evaluated. The numerical values of the reflection coefficient are presented graphically for various values of the wave number and various lengths of the barrier. It is observed that for certain small range of small values of wave number, there occurs almost total reflection. This phenomenon was also observed by Banerjea [9]. It is also observed that this phenomenon of almost total reflection of waves decreases when the length of lower part of the barrier increases and also when the length of the gap increases. We observe that by letting π‘Žβ†’0, the known results available in the literature are recovered. It is emphasized here that for numerical evaluation of the reflection and transmission coefficients, computation of singular integrals with strong singularities is avoided, and this is the main advantage in the present method over other methods existing in the literature.

In Section 2, we have presented the detailed analysis to determine the solution of the integral equation (1.1) in cases (i) and (ii) mentioned in the abstract. Then, in Section 3, we have shown the connection of (1.1) with a water wave scattering problem.

2. Solution of the Integral Equation (1.1)

2.1. Case(i)

Consider that 𝐁=(0,π‘Ž)βˆͺ(𝑏,𝑐).(2.1)

In this section, we present a method of solution of ξ€œπ΅|||πœ™(𝑑)lnπ‘‘βˆ’π‘₯|||𝑑+π‘₯𝑑𝑑=𝑓(π‘₯),π‘₯∈𝐡,(2.2)

with 𝐡≑(0,π‘Ž)βˆͺ(𝑏,𝑐) and ξ€·||||ξ€Έπœ™(𝑑)βˆΌπ‘‚π‘‘βˆ’π‘βˆ’1/2aπ‘ π‘‘βŸΆπ‘,(2.3)𝑝 being π‘Ž, 𝑏, or 𝑐.

We introduce a sectionally analytic function 𝐹(𝑧) in the complex 𝑧(=π‘₯+𝑖𝑦) plane, cut along the segment (βˆ’π‘,βˆ’π‘)βˆͺ(βˆ’π‘Ž,0)βˆͺ(0,π‘Ž)βˆͺ(𝑏,𝑐) on the real axis as given by ξ€œπΉ(𝑧)=π΅ξ‚€πœ™(𝑑)lnπ‘‘βˆ’π‘§ξ‚π‘‘+𝑧𝑑𝑑,(𝑧=π‘₯+𝑖𝑦,π‘₯∈𝐡).(2.4) If we use the following results, valid for the limiting values of Ξ¦(𝑧)=log(𝑧)≑log(π‘₯+𝑖𝑦):Φ±(π‘₯)∼ln|π‘₯|forπ‘₯>0,ln|π‘₯|Β±πœ‹π‘–forπ‘₯<0,(2.5) we find the following limiting values of 𝐹(𝑧):𝐹±⎧βŽͺβŽͺβŽͺβŽͺβŽͺ⎨βŽͺβŽͺβŽͺβŽͺβŽͺβŽ©ξ€œ(π‘₯)=𝐡|||πœ™(𝑑)lnπ‘‘βˆ’π‘₯|||ξ€œπ‘‘+π‘₯π‘‘π‘‘βˆ“πœ‹π‘–π‘₯0ξ€œπœ™(𝑑)𝑑𝑑,π‘₯∈(0,π‘Ž),π΅πœ™|||(𝑑)lnπ‘‘βˆ’π‘₯|||ξ€œπ‘‘+π‘₯π‘‘π‘‘βˆ“πœ‹π‘–π‘Ž0πœ™ξ€œ(𝑑)π‘‘π‘‘βˆ“πœ‹π‘–π‘₯π‘πœ™ξ€œ(𝑑)𝑑𝑑,π‘₯∈(𝑏,𝑐),π΅πœ™|||(𝑑)lnπ‘‘βˆ’π‘₯|||ξ€œπ‘‘+π‘₯π‘‘π‘‘Β±πœ‹π‘–0βˆ’π‘₯πœ™ξ€œ(𝑑)𝑑𝑑,π‘₯∈(βˆ’π‘Ž,0),π΅πœ™|||(𝑑)lnπ‘‘βˆ’x|||ξ€œπ‘‘+π‘₯π‘‘π‘‘Β±πœ‹π‘–0βˆ’π‘₯πœ™(𝑑)𝑑𝑑,π‘₯∈(βˆ’π‘,βˆ’π‘).(2.6) Equations (2.6) give rise to the following Plemelj-type alternative formulae: 𝐹+(π‘₯)βˆ’πΉβˆ’(π‘₯)=π‘Ÿ(π‘₯),(2.7) where ⎧βŽͺβŽͺβŽͺβŽͺβŽͺ⎨βŽͺβŽͺβŽͺβŽͺβŽͺβŽ©ξ€œπ‘Ÿ(π‘₯)=βˆ’2πœ‹π‘–π‘₯0ξ€œπœ™(𝑑)𝑑𝑑,π‘₯∈(0,π‘Ž),βˆ’2πœ‹π‘–π‘Ž0πœ™ξ€œ(𝑑)π‘‘π‘‘βˆ’2πœ‹π‘–π‘₯π‘πœ™ξ€œ(𝑑)𝑑𝑑,π‘₯∈(𝑏,𝑐),2πœ‹π‘–0βˆ’π‘₯πœ™ξ€œ(𝑑)𝑑𝑑,π‘₯∈(βˆ’π‘Ž,0),2πœ‹π‘–π‘Ž0πœ™ξ€œ(𝑑)𝑑𝑑+2πœ‹π‘–π‘βˆ’π‘₯πœ™(𝑑)𝑑𝑑,π‘₯∈(βˆ’π‘,βˆ’π‘),(2.8)𝐹+(π‘₯)+πΉβˆ’(ξ€œπ‘₯)=2𝐡|||πœ™(𝑑)lnπ‘‘βˆ’π‘₯|||𝑑+π‘₯𝑑𝑑,π‘₯∈𝐡βˆͺπ΅ξ…žξ€Έ,(2.9) where π΅ξ…ž=(βˆ’π‘,βˆ’π‘)βˆͺ(βˆ’π‘Ž,0).

By using the relation (2.9), the singular integral equation (2.2) can be expressed in the form of a Riemann-Hilbert problem as given by 𝐹+(π‘₯)+πΉβˆ’ξ€·(π‘₯)=2β„Ž(π‘₯),π‘₯∈𝐡βˆͺπ΅ξ…žξ€Έ,(2.10) where ξƒ―β„Ž(π‘₯)=𝑓(π‘₯),π‘₯∈𝐡,βˆ’π‘“(βˆ’π‘₯),π‘₯βˆˆπ΅ξ…ž.(2.11) Again noting the relation (2.8), we can write π‘Ÿξ…žξƒ―(π‘₯)=βˆ’2πœ‹π‘–πœ™(π‘₯),π‘₯∈𝐡,βˆ’2πœ‹π‘–πœ™(βˆ’π‘₯),π‘₯βˆˆπ΅ξ…ž.(2.12) Hence, 1πœ™(π‘₯)=βˆ’π‘Ÿ2πœ‹π‘–ξ…ž(π‘₯),π‘₯∈𝐡.(2.13)

Using the relation (2.13) in (2.4), we get the solution of special Riemann-Hilbert problem (2.7) as 1𝐹(𝑧)=βˆ’ξ€œ2πœ‹π‘–π΅π‘Ÿξ…ž(𝑑)lnπ‘‘βˆ’π‘§ξ‚π‘‘+𝑧𝑑𝑑.(2.14)

Utilizing the above idea, we can solve the Riemann-Hilbert problem (2.10) in the following manner.

Let 𝐹(𝑧)=𝐹0(𝑧)𝐹2(𝑧),(2.15) where 𝐹0(𝑧) is the solution of the corresponding homogeneous problem (2.10), that is, 𝐹+0(π‘₯)+πΉβˆ’0(π‘₯)=0.(2.16) Thus, (2.10) becomes 𝐹+2(π‘₯)βˆ’πΉβˆ’2(π‘₯)=2β„Ž(π‘₯)𝐹+0(π‘₯),π‘₯∈𝐡βˆͺπ΅ξ…ž.(2.17) This is a Riemann-Hilbert problem which is similar to (2.7). Noting (2.7), (2.12), and (2.14), the solution of this Riemann-Hilbert problem is given by 𝐹2(𝑧)=𝐹(𝑧)𝐹01(𝑧)=βˆ’ξ€œ2πœ‹π‘–π΅π‘‘ξƒ©ξƒ©π‘‘π‘‘2𝑓(𝑑)𝐹+0ξƒͺξ‚€(𝑑)lnπ‘‘βˆ’π‘§ξ‚ξƒͺ𝑑+𝑧𝑑𝑑+𝐸(𝑧),(2.18) where 𝐸(𝑧) is an entire function.

We can choose 𝐹0(𝑧𝑧)=2𝑧2βˆ’π‘Ž2𝑧2βˆ’π‘2𝑧2βˆ’π‘2ξ€Έξ€»βˆ’1/2(2.19) with 𝐹±0(⎧βŽͺ⎨βŽͺβŽ©Β±π‘–π‘₯)=πœŒβˆ“π‘–(π‘₯),π‘₯∈(0,π‘Ž),𝜌(π‘₯),π‘₯∈(𝑏,𝑐),(2.20) where ||π‘₯𝜌(π‘₯)=2ξ€·π‘₯2βˆ’π‘Ž2π‘₯ξ€Έξ€·2βˆ’π‘2π‘₯ξ€Έξ€·2βˆ’π‘2ξ€Έ||1/2.(2.21)

Knowing 𝐹(𝑧) from (2.18), (2.19), and (2.20), we can find the limiting values 𝐹±(π‘₯) of 𝐹(𝑧) using the formula (2.5). Substituting 𝐹±(π‘₯) in (2.7) and noting the relation (2.8), the solution of the integral equation (2.2) can be obtained as⎧βŽͺβŽͺ⎨βŽͺβŽͺβŽ©π‘‘πœ™(π‘₯)=𝑑π‘₯πœ“(π‘₯)ξ‚Άβˆ’π‘‘πœŒ(π‘₯),π‘₯∈(0,π‘Ž),𝑑π‘₯πœ“(π‘₯)ξ‚ΆπœŒ(π‘₯),π‘₯∈(𝑏,𝑐),(2.22) where 1πœ“(π‘₯)=βˆ’πœ‹2ξ€œπ‘Ž0𝑑|||𝑑𝑑(𝑓(𝑑)𝜌(𝑑))lnπ‘‘βˆ’π‘₯|||1𝑑+π‘₯𝑑𝑑+πœ‹2ξ€œπ‘π‘π‘‘|||𝑑𝑑(𝑓(𝑑)𝜌(𝑑))lnπ‘‘βˆ’π‘₯|||𝑑+π‘₯π‘‘π‘‘βˆ’πœ‹π΄π‘₯4+𝐡π‘₯3+𝐢π‘₯2ξ€Έ,+𝐷π‘₯+𝐸(2.23) and 𝐴, 𝐡, 𝐢, 𝐷, and 𝐸 are arbitrary constants.

It is observed here that the original integral equation (2.2) can be solved by the aid of formula (2.22) if and only if the following consistency conditions are satisfied: (𝑖)πœ“(0)=0,(𝑖𝑖)πœ“ξ…ž(0)=0,(𝑖𝑖𝑖)πœ“(π‘Ž)=0,(𝑖𝑣)πœ“(𝑏)=0,(𝑣)πœ“(𝑐)=0.(2.24) We can determine the arbitrary constants 𝐴, 𝐡, 𝐢, 𝐷, and 𝐸 by using the relation (2.24). Thus, πœ™(π‘₯) can be determined completely by utilizing the relation (2.22). It is observed that the solution as obtained by the relations (2.22) and (2.23) of (2.2) neither involves any strong singularity nor any arbitrary constant.

2.2. Case(ii)

Consider that 𝐁=𝐡1=βˆͺ𝑛𝑖=1ξ€·π‘Žπ‘–,𝑏𝑖.(2.25) In this section, we proceed to solve the integral equation ξ€œπ΅1πœ™1(|||𝑑)lnπ‘‘βˆ’π‘₯|||𝑑+π‘₯𝑑𝑑=𝑓1(π‘₯),π‘₯∈𝐡1,(2.26) where 𝐡1 consists of 𝑛 finite disjoint intervals 𝐡1=βˆͺ𝑛𝑖=1(π‘Žπ‘–,𝑏𝑖) and πœ™1(𝑑)βˆΌπ‘‚(|π‘‘βˆ’π‘π‘–|)βˆ’1/2 as 𝑑→𝑝𝑖, 𝑝𝑖 being π‘Žπ‘–, 𝑏𝑖, 𝑖=1,…,𝑛.

In this case, we introduce the following sectionally analytic function 𝐹1(ξ€œπ‘§)=𝐡1πœ™1(𝑑)lnπ‘‘βˆ’π‘§ξ‚ξ€·π‘‘+𝑧𝑑𝑑,𝑧=π‘₯+𝑖𝑦,π‘₯∈𝐡1ξ€Έ,(2.27) cut along .

Using the formula (2.5), we get following limiting values of : This gives the following Plemelj-type alternative formulae for the integral equation (2.26): where Hence, from the relation (2.30), we get Thus, noting the relation (2.32), the equation (2.27) becomes Utilizing the above idea, we can solve the Riemann-Hilbert problem (2.31) in the following form: where is an entire function, and we can choose , the solution of the homogeneous equation of (2.31) as with where

Using (2.35) in (2.34), is known. Noting the relation (2.5), one can find the limiting values of . Using the limiting values , we obtain, from (2.29) and (2.30), the solution of the integral equation (2.26) as where We observe that the original integral equation (2.26) can be solved by using the above relations if and only if the following consistency conditions are satisfied:πœ“1ξ€·π‘Žπ‘–ξ€Έ=πœ“1𝑏𝑖=0,𝑖=1,…,𝑛.(2.40) From these 2𝑛 conditions, we find the arbitrary constant Ξ 2𝑛+1𝑖=2𝐴𝑖, and from the order condition, the arbitrary constant 𝐴1 can be determined, and finally the solution of the singular integral equation (2.26) can be obtained from (2.38).

In the next section, we give the genesis of the integral equation (2.2) with 𝐡=(0,π‘Ž)βˆͺ(𝑏,𝑐) and use its solution to the problem of scattering of water waves by a partially immersed vertical barrier with gap.

3. Genesis of the Integral Equation (2.2) Connection with a Water Wave Problem

3.1. The Detailed Analysis and the Analytical Solution

In this section, we show the connection of integral equation (2.2) with the problem of scattering of surface water waves by a thin partially immersed vertical barrier with a gap in deep water occupying the position π‘₯=0, 0<𝑦<π‘Ž, and 𝑏<𝑦<𝑐, 𝑦-axis being directed downwards and π‘₯-axis along the mean free surface. We assume the motion to be two dimensional, irrotational, and time harmonic. The velocity potential is represented by Re(πœ™(π‘₯,𝑦)exp(βˆ’π‘–πœŽπ‘‘)), 𝜎 being circular frequency. Here, πœ™(π‘₯,𝑦) satisfies the following boundary value problem:βˆ‡2πœ™=0in𝑦β‰₯0,𝐾Φ+πœ™π‘¦=0on𝑦=0,(3.1) where 𝐾=𝜎2/𝑔, 𝑔 being the acceleration of gravity, πœ™π‘₯π‘Ÿ=0onπ‘₯=0,π‘¦βˆˆπ΅β‰‘(0,π‘Ž)βˆͺ(𝑏,𝑐),(3.2)1/2ξ€½βˆ‡πœ™isboundedasπ‘ŸβŸΆ0,π‘Ÿ=(π‘₯)2+(π‘¦βˆ’π‘£)2ξ€Ύ1/2πœ™ξƒ―,𝑣=π‘Ž,𝑏or𝑐,(3.3)βˆ‡πœ™βŸΆ0asπ‘¦βŸΆβˆž,(3.4)(π‘₯,𝑦)=𝑅exp(βˆ’πΎπ‘¦βˆ’π‘–πΎπ‘₯)+exp(βˆ’πΎπ‘¦+𝑖𝐾π‘₯)asπ‘₯βŸΆβˆ’βˆž,𝑇exp(βˆ’πΎπ‘¦+𝑖𝐾π‘₯)asπ‘₯⟢∞,(3.5) where 𝑅 and 𝑇 are reflection and transmission coefficients. A suitable representation of πœ™ satisfying the equations (3.1), (3.4), and (3.5) is given by ⎧βŽͺ⎨βŽͺβŽ©ξ€œπœ™(π‘₯,𝑦)=𝑅exp(βˆ’πΎπ‘¦βˆ’π‘–πΎπ‘₯)+exp(βˆ’πΎπ‘¦+𝑖𝐾π‘₯)+∞0ξ€œπ΅(π‘˜)𝐿(π‘˜,𝑦)exp(π‘˜π‘₯)π‘‘π‘˜,π‘₯<0,𝑇exp(βˆ’πΎπ‘¦+𝑖𝐾π‘₯)+∞0𝐴(π‘˜)𝐿(π‘˜,𝑦)exp(βˆ’π‘˜π‘₯)π‘‘π‘˜,π‘₯>0,(3.6) where, 𝐿(π‘˜,𝑦)=π‘˜cosπ‘˜π‘¦βˆ’πΎsinπ‘˜π‘¦.(3.7) By Havelock's expansion theorem, 𝐴(π‘˜)=βˆ’π΅(π‘˜),𝑅=1βˆ’π‘‡.(3.8) Then, using the condition (3.2) along with the condition of continuity of πœ™ across the gap 𝐺∢(π‘Ž,𝑏)βˆͺ(𝑐,∞), along π‘₯=0, we arrive at the following multiple integral equations: ξ€œβˆž0π‘˜π΄(π‘˜)𝐿(π‘˜,𝑦)π‘‘π‘˜=𝑖𝐾(1βˆ’π‘…)exp(βˆ’πΎπ‘¦),π‘¦βˆˆπ΅,(3.9)ξ€œβˆž0𝐴(π‘˜)𝐿(π‘˜,𝑦)π‘‘π‘˜=𝑅exp(βˆ’πΎπ‘¦),π‘¦βˆˆπΊ.(3.10) The Relation (3.9) can be recast in the form π‘‘ξ€œπ‘‘π‘¦βˆž0𝐴(π‘˜)π‘˜ξƒ―πΏ(π‘˜,𝑦)π‘‘π‘˜=𝑖(1βˆ’π‘…)exp(βˆ’πΎπ‘¦)+𝐷1,0<𝑦<π‘Ž,𝑖(1βˆ’π‘…)exp(βˆ’πΎπ‘¦)+𝐷2,𝑏<𝑦<𝑐,(3.11) where 𝐷1 and 𝐷2 are two arbitrary constants. Then, applying 𝑑/𝑑𝑦+𝐾 to the equations (3.10) and (3.11), we obtain ξ€œβˆž0𝐹(π‘˜)sinπ‘˜π‘¦π‘‘π‘˜=0,π‘¦βˆˆπΊ,(3.12)ξ€œβˆž0𝐹(π‘˜)π‘˜ξƒ―sinπ‘˜π‘¦π‘‘π‘˜=𝐾𝐷1𝑦+𝐢1,0<𝑦<π‘Ž,𝐾𝐷2𝑦+𝐢2,𝑏<𝑦<𝑐,(3.13) where πΉξ€·π‘˜(π‘˜)=βˆ’2+𝐾2𝐴(π‘˜).(3.14) It is observed from the relation (3.13) that we must choose 𝐢1=0 for consistency. Setting ξ€œβˆž0𝐹(π‘˜)sinπ‘˜π‘¦π‘‘π‘˜=𝑔(𝑦),π‘¦βˆˆπ΅,(3.15) where 𝑔(𝑦) is an unknown function and using sine inversion formula, to the relations (3.12) and (3.15), we obtain from (3.13) ξ€œπ΅π‘”||||(𝑑)logπ‘¦βˆ’π‘‘||||𝑦+𝑑=𝑓(𝑦),π‘¦βˆˆπ΅,(3.16) with 𝑓(𝑦)=βˆ’πœ‹πΎπ·1𝑦,0<𝑦<π‘Ž,βˆ’πœ‹πΎπ·2π‘¦βˆ’πœ‹πΆ2,𝑏<𝑦<𝑐.(3.17) The solution of the integral equation (3.16) is given by the equations (2.22), (2.23), and (2.24) with 𝑓(𝑦) given by relation (3.17).

Thus, knowing 𝑔(𝑑), we can obtain 𝐹(π‘˜) by using Fourier Sine inversion formulae to the relations (3.12) and (3.15). Hence, 𝐴(π‘˜) can be obtained from equation (3.14). Now, 𝑓(𝑦) contains unknown constants 𝐷1, 𝐷2, and 𝐢2. To obtain the unknowns 𝐷1, 𝐷2, 𝐢2, and 𝑅, we substitute 𝐴(π‘˜) in (3.10) and (3.11) to get the following four equations: ξ€œπ‘π‘πœ“πœŒπ‘’πΎπ‘‘πΎξ‚Έβˆ’ξ€œπ‘‘π‘‘=0,π‘Ž0πœ“πœŒπ‘’βˆ’πΎπ‘‘ξ€œπ‘‘π‘‘+π‘π‘πœ“πœŒπ‘’βˆ’πΎπ‘‘ξ€œπ‘‘π‘‘βˆ’π‘Ž0πœ“πœŒπ‘’πΎπ‘‘ξ‚Ήξ€œπ‘‘π‘‘=𝑅,βˆ’πΎπ‘π‘Žπœ“πœŒπ‘’πΎπ‘‘π‘‘π‘‘βˆ’πΎ2𝐷1ξ€œπ‘Ž0π‘₯𝑒𝐾π‘₯+𝑒𝐾𝑏𝐾𝑏𝐷2βˆ’π·2+𝐢2ξ€Έ+𝐷1πΎξ€œ=0,π‘π‘Žπœ“πœŒπ‘’βˆ’πΎπ‘‘ξ€œπ‘‘π‘‘βˆ’πΎβˆžπ‘πœ“πœŒπ‘’βˆ’πΎπ‘‘πΎπ‘‘π‘‘βˆ’πœ‹ξ€œπ‘Ž0ξ€·βˆ’πœ‹πΎπ‘₯𝐷1ξ€Έπ‘’βˆ’πΎπ‘₯βˆ’πΎπœ‹ξ€œπ‘π‘ξ€·βˆ’πœ‹πΎπ‘₯𝐷2βˆ’πœ‹πΆ2ξ€Έπ‘’βˆ’πΎπ‘₯βˆ’π·1=𝑖(1βˆ’π‘…),(3.18) where 1πœ“(π‘₯)=βˆ’πœ‹2ξ€œπ‘Ž0𝑑|||𝑑𝑑(𝑓(𝑑)𝜌(𝑑))lnπ‘‘βˆ’π‘₯|||1𝑑+π‘₯𝑑𝑑+πœ‹2ξ€œπ‘π‘π‘‘|||𝑑𝑑(𝑓(𝑑)𝜌(𝑑))lnπ‘‘βˆ’π‘₯|||𝑑+π‘₯π‘‘π‘‘βˆ’πœ‹π΄π‘₯4+𝐡π‘₯3+𝐢π‘₯2ξ€Έ,+𝐷π‘₯+𝐸(3.19) and 𝐴, 𝐡, 𝐢, 𝐷, and 𝐸 are arbitrary constants and ||π‘₯𝜌(π‘₯)=2ξ€·π‘₯2βˆ’π‘Ž2π‘₯ξ€Έξ€·2βˆ’π‘2π‘₯ξ€Έξ€·2βˆ’π‘2ξ€Έ||1/2.(3.20) Thus, the relations (3.18) can be solved to determine four unknowns 𝐷1, 𝐷2, 𝐢2, and 𝑅. This completes the description of the analytical solution of the water wave problem.

3.2. Numerical Results

The reflection coefficient |𝑅| has been numerically evaluated for various values of wave number 𝐾𝑏 for four cases, namely, π‘Ž/𝑏→0, 𝑐/𝑏=15; π‘Ž/𝑏=0.5, 𝑐/𝑏=4; π‘Ž/𝑏=0.5, 𝑐/𝑏=6; π‘Ž/𝑏=0.25, 𝑐/𝑏=6.

Numerical computation of |𝑅| involves the solution of four equations (3.18). It may be noted from the relation (3.19) that for the solution of (3.18), one does not need to compute any integral with strong singularity. The integrals occurring in (3.18) are defined in Riemann sense. This is the main advantage in the present method as compared to the other methods based on strongly singular integral equation formulation available in the literature. Since there is no integral with strong singularity, so for numerical computation, we have used Mathematica 5.2 to evaluate various integrals in (3.18) and (3.19) and to solve the system of (3.18).

In Table 1, we have presented the values of reflection coefficient |𝑅𝐴| computed from our result when π‘Ž/𝑏→0 and 𝑐/𝑏=15 and |𝑅𝐸|, the reflection coefficients corresponding to the problem of scattering by a submerged vertical barrier π‘₯=0, 𝑏<𝑦<𝑐(𝑐/𝑏=15) as given in the literature (cf. Banerjea and Dutta [6]) for various values of wave number 𝐾𝑏. It is observed that the values reflection coefficient |𝑅𝐴| coincides with |𝑅𝐸| upto two or three places of decimal. This shows that in limiting case, the reflection coefficient |𝑅𝐴| computed from the present paper agrees with the exact value of |𝑅𝐸| as given in [6].

Table 1

In Figures 1, 2(a), 2(b), and 3 (f), |𝑅|  is depicted against 𝐾𝑏 for various values of π‘Ž/𝑏 and 𝑐/𝑏=6. In each case, we have verified that |𝑅|2+|𝑇|2=1 for all values of 𝐾𝑏.

Figure 1
Figure 2
Figure 3

In Figure 1, we have plotted |𝑅| against various values of 𝐾𝑏 when π‘Ž/𝑏=0.5 and 𝑐/𝑏=4. For 0.15<𝐾𝑏<0.25, it is observed that |𝑅| rises steeply, becomes almost a unity, and then falls rapidly forming a sharp peak. This indicates that for this range of values of 𝐾𝑏, there occurs almost total reflection of wave energy. Also for 0.35<𝐾𝑏<0.4, there occurs a small oscillation in |𝑅|. In Figures 2(a) and 2(b),  |𝑅| is depicted against 𝐾𝑏 for π‘Ž/𝑏=0.5 and 𝑐/𝑏=6. In Figure 2(a), it is observed that for very small values of 𝐾𝑏, that is, for 𝐾𝑏<0.1, |𝑅| increases at first and then decreases almost to zero indicating that most of the wave energy is transmitted. If the values of 𝐾𝑏 are further increased, that is, 𝐾𝑏>0.1, it is observed from Figure 2(b) that |𝑅| increases. For 0.13<𝐾𝑏<0.15, |𝑅| sharply increases almost to unity and then decreases sharply. The peak value in this case is lower than the peak value when 𝑐/𝑏=4 as observed in Figure 1. Thus, when the length of the lower part of the barrier increases, the peak value in |𝑅| decreases. In Figure 3, |𝑅| is plotted against 𝐾𝑏 for π‘Ž/𝑏=0.25 and 𝑐/𝑏=6. In this figure, it is observed that when the length of the gap is increased keeping the length of the lower part of the barrier the same, the sharp peak in the values of |𝑅| occurs for 0.111<𝐾𝑏<0.1111. In this case the peak value is much lower than the peak value observed in Figures 1 and 2(b). Thus, an increase in the length of the gap decreases the peak value in |𝑅|. Thus, it is observed from Figures 1, 2, and 3 that for certain range of values of 𝐾𝑏, there occurs a sharp peak in |𝑅| showing that |𝑅| reaches a maximum value. In these three cases, the maximum value |𝑅| is almost unity, that is, the incident wave experiences almost total reflection for certain range of small values of 𝐾𝑏. This type of behaviour of |𝑅| is also observed by Banerjea in [9].

4. Conclusion

A direct function theoretic method is applied here to solve the integral equations (2.2) and (2.26) without converting them to the integral equation with strong singularity. It may be noted here that the solution functions given in (2.22) and (2.38) do not involve any strong singularity. Also using the presently detailed solution of (2.2), the problem of scattering of water waves by a partially immersed barrier with a gap is examined fully.


This work is also partially supported by DST through PURSE scheme (no. SR/S9/Z-23/2008/5) at Jadavpur University. A. Chakrabarti is thankful to NASI, Allahabad for offering Senior Scientist Platinum Jubilee Fellowship which helped him in pursuing this research.


  1. L. Lewin, Theory of Wave Guides, Newness Butterworths, London, UK, 1975.
  2. F. Ursell, β€œThe effect of a fixed vertical barrier on surface waves in deep water,” Proceedings of the Cambridge Philosophical Society, vol. 43, pp. 374–382, 1947. View at Google Scholar Β· View at Zentralblatt MATH
  3. B. N. Mandal and P. K. Kundu, β€œScattering of waterwaves by vertical barriers and associated Mathematical methods,” Proceedings of the Indian National Science Academy, vol. 53, pp. 514–530, 1987. View at Google Scholar
  4. A. Chakrabarti and S. R. Manam, β€œSolution of a logarithmic singular integral equation,” Applied Mathematics Letters, vol. 16, no. 3, pp. 369–373, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  5. S. R. Manam, β€œA logarithmic singular integral equation over multiple intervals,” Applied Mathematics Letters, vol. 16, no. 7, pp. 1031–1037, 2003. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  6. S. Banerjea and B. Dutta, β€œSolution of a singular integral equation and its application to water wave problems,” Indian Institute of Science, vol. 86, no. 3, pp. 265–278, 2006. View at Google Scholar
  7. S. Banerjea and C. Sarkar, β€œOn a singular integral equation with log kernel and its application,” International Journal of Mathematics and Mathematical Sciences, vol. 2006, Article ID 35630, 12 pages, 2006. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  8. A. Chakrabarti, β€œSolution of certain weakly singular integral equations,” IMA Journal of Applied Mathematics, vol. 71, no. 4, pp. 534–543, 2006. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  9. S. Banerjea, β€œScattering of water waves by a vertical wall with gaps,” Australian Mathematical Society, vol. 37, no. 4, pp. 512–529, 1996. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH