# The Computation of the Magnitude of the Far Field for an Eccentric Circular Cylinder

Authors: Bülent Yılmaz
Journal: Advances in Mathematical Physics
 (2015)

---

## Abstract

The specific case of scattering of a plane wave by a two-layered penetrable eccentric circular cylinder has been considered and it is about the validity of the on surface radiation condition method and its applications to the scattering of a plane wave by a two-layered penetrable eccentric circular cylinder. The transformation of the problem of scattering by the eccentric circular cylinder to the problem of scattering by the concentric circular cylinder by using higher order radiation conditions, is observed. Numerical results presented the magnitude of the far field.


---

## Body

## 1. Introduction

Approximate techniques have been introduced to study the scattering of waves by obstacles. The OSRC method has been devised by Kriegsmann et al. to investigate electromagnetic scattering problems involving cylindrical convex objects [1]. The main concept of this method is the application of a radiation boundary condition (RBC), connecting the field and its normal derivative, directly onto the surface of the scatterer, to determine approximately the surface field or its derivative in terms of the given field.Two main approaches have been employed to derive RBCs [2, 3]. One of the methods is based on the idea of killing the terms of the expansion of the scattering field satisfying the Helmholtz equation and Sommerfeld radiation condition. An nth-order RBC operator which annihilates the first n terms in the expansion is obtained either on a large circular cylinder enclosing a cylindrical convex object or on a large sphere enclosing a finite convex object, depending on the geometrical dimensions of the problem. These RBCs can be generalized so that they can be used in the OSRC method for constructing the approximate solution of a scattering problem involving an arbitrary convex object.The second approach which is capable of supplying SRCs of any order on a smooth scattering surface of arbitrary form, thus avoiding the difficulty mentioned above, was introduced by Kriegsmann and Moore [3]. This method is based on an asymptotic expansion, similar to the Luneburg-Kline expansion, made in the neighbourhood of a phase front by assuming that the field has a well defined phase front in the regions of interest. Then the assumption that the surface of the scattering is a phase front yields the SRC of the OSRC method. Although the method of Jones may, in principle, furnish SRCs to any desired order, only the second-order condition has been produced and applied to various problems to examine the predictions of the method. From these investigations it has been observed that, for hard objects, as frequency increases, second-order SRCs begin to fail to account for creeping wave physics adequately.Considering these results, it has been conjectured that, by introducing higher-order SRCs in the method, creeping wave physics may be modelled more accurately; that is, the approximation of the OSRC method may be improved. In [1, 3–7] only the first- and second-order RBCs have been produced and used in conjunction with the OSRC method. Later [2, 8, 9] third- and fourth-order RBCs have been used to examine whether the use of higher-order SRCs in the OSRC method models creeping wave physics more accurately than a second-order SRC.In this work we applied the higher-order SRCs to scattering of plane waves by an eccentric penetrable circular cylinder. The results are compared with those of second- and fourth-order SRCs and with those of the scattering of plane waves by a concentric penetrable cylinder [9].This paper is organized as follows. The formulation of the problem together with the exact and the approximate solutions of these equations with OSRC method is presented and calculated in Sections2 and 3, respectively. In Section 4 the comparisons and some concluding remarks take place.


## 2. Formulation and Exact Solution

We denote the outer and the inner layers byB
2, with sound speed c
2 (m/s) and constant density ρ
2 (kg/m3), and by B
3, with sound speed c
3 (m/s) and constant density ρ
3 (kg/m3), respectively. The interface between them is the circular cylinder Σ
2 of radius r
2. The radius of the circular cylinder between the exterior region B
1 which is a homogeneous isotropic medium with sound speed c
1 (m/s) and constant density ρ
1 (kg/m3) and the region B
2 is r
1,  B
1
=
r
>
r
1,  B
2
=
r
2
<
r
<
r
1, and  B
3
=
r
<
r
2.Let us denote byO
1 the center of the outer cylinder Σ
1 and by O
2 the center of the inner cylinder Σ
2 and take O
1 as the origin of the rectangular Cartesian coordinate system (
x
,
y
). We suppose that the x-axis passes through the points O
1 and O
2. In addition to this, we consider two polar coordinate systems centered at O
1 and O
2. Radial and angular variables are denoted by (
r
,
ϕ
) and (
ρ
,
θ
), respectively. ϕ and θ are measured from the positive x-axis. Under these assumptions the geometry of the problem is defined in Figure 1, where r
1, r
2 are the outer and inner cylinder radii, while the distance between the centers of the cylinders is represented by e.Figure 1The scattering problem is described by the following equations and boundary conditions:(1)
∇
2
u
l
+
k
l
2
u
l
=
0
x
,
y
∈
B
l
for
k
l
=
w
c
l
,
l
=
1,2
,
3
,
(2)
u
1
+
u
i
=
u
2
,
∂
∂
r
u
1
+
u
i
=
ξ
1
∂
u
2
∂
r
on
⁡
r
=
r
1
,
(3)
u
2
=
u
3
,
∂
u
2
∂
r
=
ξ
2
∂
u
3
∂
r
on
⁡
ρ
=
r
2
,
where ζ
1
=
ρ
1
/
ρ
2,  ζ
2
=
ρ
2
/
ρ
3,  k
l
=
w
/
c
l is a wave number, and w is the angular frequency. In addition, at infinity the scattered field u
1 must have the form of a radiating wave; that is, the following Sommerfeld radiation condition must be satisfied:
(4)
lim
r
→
∞
r
1
/
2
∂
u
1
∂
r
+
i
k
1
u
1
=
0
.Sinceϕ
0 is the direction of propagation of the incident wave, the propagation vector of the incident wave is assumed to be p
=
(
-
cos
⁡
ϕ
0
,
-
sin
ϕ
0
). ϕ
,
r are the polar coordinates and the assumed time dependence is e
-
i
w
t. All the field quantities are then independent of z. If x
=
r
cos
⁡
ϕ,  y
=
r
sin
ϕ then the incident wave can be expressed as
(5)
u
i
=
e
-
i
k
1
p
·
x
=
e
i
k
1
r
cos
⁡
(
ϕ
-
ϕ
0
)
=
∑
n
=
-
∞
∞
‍
J
n
(
k
1
r
)
e
i
n
(
ϕ
-
ϕ
0
+
π
/
2
)
[4]. Solutions for u
1, u
2, and u
3 can be expressed in polar coordinates as follows:
(6)
u
1
=
∑
n
=
-
∞
∞
‍
a
n
H
n
(
k
1
r
)
e
i
n
(
ϕ
-
ϕ
0
+
π
/
2
)
,
(7)
u
2
=
∑
n
=
-
∞
∞
‍
b
n
J
n
(
k
2
ρ
)
+
c
n
H
n
(
k
2
ρ
)
e
i
n
(
θ
-
π
/
2
)
,
(8)
u
3
=
∑
n
=
-
∞
∞
‍
d
n
J
n
(
k
3
ρ
)
e
i
n
(
θ
-
π
/
2
)
.
u
2 solution is expressed in polar coordinates centered at O
2. In order to be able to use this solution for r
=
r
1 it has to be written first in polar coordinates r
,
ϕ. We know that J
n
(
k
2
ρ
) and H
n
(
k
2
ρ
) can be expressed as follows:
(9)
J
n
k
2
ρ
e
i
n
θ
=
∑
q
=
-
∞
∞
‍
J
q
k
2
r
J
q
-
n
k
2
e
e
i
q
ϕ
,
H
n
k
2
ρ
e
i
n
θ
=
∑
q
=
-
∞
∞
‍
H
q
k
2
r
J
q
-
n
k
2
e
e
i
q
ϕ
.
If (9) is inserted in (7) the solution for u
2 in polar coordinates (
r
,
ϕ
) takes the following mentioned form:
(10)
u
2
=
∑
n
=
-
∞
∞
‍
b
n
∑
q
=
-
∞
∞
‍
J
q
-
n
(
k
2
e
)
J
q
(
k
2
r
)
e
i
q
ϕ
e
-
i
n
π
/
2
+
∑
n
=
-
∞
∞
‍
c
n
∑
q
=
-
∞
∞
‍
J
q
-
n
(
k
2
e
)
H
q
(
k
2
r
)
e
i
q
ϕ
e
-
i
n
π
/
2
.
Now, if (5), (6), and (10) are used under the boundary conditions (2), in order for these conditions to be satisfied
(11)
J
q
(
k
1
r
1
)
e
-
i
q
(
ϕ
0
-
π
/
2
)
+
a
q
H
q
(
k
1
r
1
)
e
-
i
q
(
ϕ
0
-
π
/
2
)
=
J
q
(
k
2
r
1
)
∑
n
=
-
∞
∞
‍
b
n
J
q
-
n
(
k
2
e
)
e
-
i
n
π
/
2
+
H
q
(
k
2
r
1
)
∑
n
=
-
∞
∞
‍
c
n
J
q
-
n
(
k
2
e
)
e
-
i
n
π
/
2
(12)
k
1
J
q
′
(
k
1
r
1
)
e
-
i
q
(
ϕ
0
-
π
/
2
)
+
k
1
a
q
H
q
′
(
k
1
r
1
)
e
-
i
q
(
ϕ
0
-
π
/
2
)
=
k
2
ξ
1
J
q
′
(
k
2
r
1
)
∑
n
=
-
∞
∞
‍
b
n
J
q
-
n
(
k
2
e
)
e
-
i
n
π
/
2
+
k
2
ξ
1
H
q
′
(
k
2
r
1
)
∑
n
=
-
∞
∞
‍
c
n
J
q
-
n
(
k
2
e
)
e
-
i
n
π
/
2
(13)
e
=
0
J
q
-
n
(
k
2
e
)
=
1
;
q
=
n
0
;
q
≠
n
.
Thus, if ϕ
0
=
π is taken it is seen that in case of e
=
0 or k
2
e
→
0  (11) and (12) are, respectively, transformed to (2.8) and (2.9) in [9]. Then if (7) and (8) are used with the boundary conditions (3), the requirements for these conditions to be satisfied are
(14)
b
n
J
n
k
2
r
2
+
c
n
H
n
k
2
r
2
=
d
n
J
n
k
3
r
2
,
(15)
k
2
b
n
J
n
′
k
2
r
2
+
k
2
c
n
H
n
′
k
2
r
2
=
ξ
2
k
3
d
n
J
n
′
k
3
r
2
.
Equations  (11), (12), (14), and (15) constitute a system of algebraic equations with an infinite dimension for a
n, b
n, c
n, and d
n. First the infinite sums in (11) and (12) are converted to finite sums by choosing a naturel number N
1. The summation is implemented up to N
1. Therefore a system of infinite dimension is approximated by this finite system and is numerically solved. Then the same process is repeated by choosing a greater naturel number N
2. This is repeated until a solution of required precision is obtained. In this process, in order to diminish the load of the numerical calculations, instead of dealing with all (11), (12), (14), and (15) an algebraic system of equations with a finite dimension which includes only b
n’s can be derived as follows. By eliminating a
n from (11) and (12) the below mentioned system of equations which includes only b
n and c
n is obtained:
(16)
∑
l
=
-
∞
∞
‍
Λ
l
(
n
)
b
l
+
Φ
l
(
n
)
c
l
=
Ψ
n
.
Here the elements of matrices Λ
l
(
n
), Φ
l
(
n
) and vector Ψ
n are defined as follows:
(17)
Λ
l
n
=
J
n
-
l
(
k
2
e
)
e
-
i
l
π
/
2
+
i
n
(
ϕ
0
-
π
/
2
)
·
J
n
k
2
r
1
H
n
k
1
r
1
-
k
2
ξ
1
J
n
′
k
2
r
1
k
1
H
n
′
k
1
r
1
,
Φ
l
n
=
J
n
-
l
(
k
2
e
)
e
-
i
l
π
/
2
+
i
n
(
ϕ
0
-
π
/
2
)
·
H
n
k
2
r
1
H
n
k
1
r
1
-
k
2
ξ
1
H
n
′
k
2
r
1
k
1
H
n
′
k
1
r
1
,
Ψ
n
=
J
n
(
k
1
r
1
)
H
n
(
k
1
r
1
)
-
J
n
′
(
k
1
r
1
)
H
n
′
(
k
1
r
1
)
.
On the other side if d
n is eliminated from (14) and (15) the relation
(18)
c
n
=
ζ
n
η
n
b
n
is obtained. ζ
n and η
n are defined as follows:
(19)
ζ
n
=
k
3
ξ
2
J
n
′
(
k
3
r
2
)
J
n
(
k
2
r
2
)
k
2
H
n
(
k
2
r
2
)
H
n
′
(
k
2
r
2
)
-
J
n
k
3
r
2
J
n
′
k
2
r
2
H
n
k
2
r
2
H
n
′
k
2
r
2
,
(20)
η
n
=
-
k
3
ξ
2
J
n
′
(
k
3
r
2
)
k
2
H
n
′
(
k
2
r
2
)
+
J
n
(
k
3
r
2
)
H
n
(
k
2
r
2
)
.
If c
n in (16) is replaced by the right hand side of (18), the below system of equations is obtained for b
n:
(21)
∑
l
=
-
∞
∞
‍
A
l
n
b
l
=
Ψ
n
.
Here the elements of the A
l
(
n
) matrix are defined as follows:
(22)
A
l
n
=
Λ
l
n
+
Φ
l
n
ζ
l
η
l
.
By replacing l and n by the elements of the increasing whole number sequence N
1
,
N
2
,
N
3
,
… and by making (21) finite, finite number of b
n is found with the same process as described above. By using b
n’s in (18), (14) (or (15)), and (11) (or (12)) the same numbers of c
n, d
n, and a
n’s are calculated. Thus, the scattered field is approximated by a finite sum.


## 3. Approximate Solution

By using on surface radiation condition method, the problem is converted to an internal one together with the boundary conditions of the impedance [9] type on Σ
1, boundary conditions (3) on Σ
2, and (1) (for l
=
2,3). For this reason in the boundary condition of the impedance type [9] θ is replaced by ϕ.When the incident wave is a plane wave,u
i
(
a
,
ϕ
)
=
v
i
(
ϕ
) and 1
/
k
1
∂
u
i
/
∂
r
(
a
,
ϕ
)
=
w
i
(
ϕ
)
(23)
v
i
ϕ
=
∑
n
=
-
∞
∞
‍
J
n
(
ε
)
e
i
n
ϕ
-
ϕ
0
+
π
/
2
,
ε
=
k
1
r
1
,
w
i
(
ϕ
)
=
∑
n
=
-
∞
∞
‍
J
n
′
(
ε
)
e
i
n
ϕ
-
ϕ
0
+
π
/
2
are written. Solutions of (1) (for l
=
2,3) can be expressed in the following forms:
(24)
u
2
=
∑
n
=
-
∞
∞
‍
b
~
n
J
n
(
k
2
ρ
)
+
c
~
n
H
n
(
k
2
ρ
)
e
i
n
θ
-
π
/
2
,
(25)
u
3
=
∑
n
=
-
∞
∞
‍
d
~
n
J
n
(
k
3
ρ
)
e
i
n
θ
-
π
/
2
.
The expression of u
2 in (24) is written in (
ρ
,
θ
) coordinates. From (10) this solution can be expressed in (
ρ
,
ϕ
) coordinates as follows:
(26)
u
2
=
∑
n
=
-
∞
∞
‍
∑
q
=
-
∞
∞
‍
b
~
q
J
n
-
q
(
k
2
e
)
e
-
i
q
π
/
2
J
n
(
k
2
r
)
e
i
n
ϕ
+
∑
n
=
-
∞
∞
‍
∑
q
=
-
∞
∞
‍
c
~
q
J
n
-
q
(
k
2
e
)
e
-
i
q
π
/
2
H
n
(
k
2
r
)
e
i
n
ϕ
.
When u
2 is taken as given in (26) and (23) are chosen for the incident wave, in order for (2.3) condition and (2.21) of [9] to be satisfied
(27)
Θ
(
l
)
(
ε
)
J
n
(
k
2
r
1
)
-
k
2
ξ
1
k
1
J
n
′
(
k
2
r
1
)
·
∑
q
=
-
∞
∞
‍
b
~
q
J
n
-
q
(
k
2
e
)
e
-
i
q
π
/
2
+
Θ
(
l
)
(
ε
)
H
n
(
k
2
r
1
)
-
k
2
ξ
1
k
1
H
n
′
(
k
2
r
1
)
·
∑
q
=
-
∞
∞
‍
c
~
q
J
n
-
q
(
k
2
e
)
e
-
i
q
π
/
2
=
Θ
(
l
)
(
ε
)
J
n
(
ε
)
-
J
n
′
(
ε
)
e
-
i
n
(
ϕ
0
-
π
/
2
)
is required. Here the function Θ
(
l
)
(
ε
) is defined in (2.32) of [9].Fore
=
0 two cylinders have the same axis and relation (13) holds. Thus if we take ϕ
0
=
π at the k
2
e
→
0 limit we can easily observe that (27) is converted to that of the problem which is the same as in [9]. If solutions (24) and (25) are used under conditions (3) in order for those conditions to be satisfied it is required that
(28)
b
~
n
J
n
k
2
r
2
+
c
~
n
H
n
k
2
r
2
=
d
~
n
J
n
k
3
r
2
,
(29)
k
2
b
~
n
J
n
′
(
k
2
r
2
)
+
k
2
c
~
n
H
n
′
(
k
2
r
2
)
=
ξ
2
k
3
d
~
n
J
n
′
(
k
3
r
2
)
hold. As a result (27), (28), and (29) constitute an algebraic system of equations of infinite dimension, for b
~
n, c
~
n, and d
~
n.We can write (27) as follows:
(30)
∑
q
=
-
∞
∞
‍
Λ
~
q
(
n
)
b
~
l
+
Φ
~
q
(
n
)
c
~
l
=
Ψ
~
n
.
Here Λ
~
q
(
n
), Φ
~
q
(
n
), and Ψ
~
n are defined as follows:
(31)
Λ
~
q
n
=
J
n
-
q
(
k
2
e
)
e
-
i
q
π
/
2
+
i
n
ϕ
0
-
π
/
2
·
Θ
l
,
m
ε
J
n
k
2
r
1
-
k
2
ξ
1
k
1
J
n
′
k
2
r
1
,
Φ
~
q
n
=
J
n
-
q
(
k
2
e
)
e
-
i
q
π
/
2
+
i
n
ϕ
0
-
π
/
2
·
Θ
l
,
m
ε
H
n
k
2
r
1
-
k
2
ξ
1
k
1
H
n
′
k
2
r
1
Ψ
~
n
=
Θ
l
,
m
ε
J
n
ε
-
J
n
′
ε
.
System (30) is the approximation prescribed by the method to system (16) that is derived for the exact solution. It is observed that the replacement of Θ
(
m
)
(
ε
) by H
n
′
(
ε
)
/
H
n
(
ε
) in (31) yields the exact ones. It is easily seen that the expressions Λ
~
q
(
n
), Φ
~
q
(
n
), and Ψ
~
n are transformed to H
n
(
2
)
′
(
ε
)
Λ
q
(
n
), H
n
(
2
)
′
(
ε
)
Φ
q
(
n
), and H
n
(
2
)
′
(
ε
)
Ψ
n. Therefore (30) is transformed to (16). Thus, for the problem under consideration the surface radiation conditions method is equivalent to introducing the approximation Θ
(
l
)
(
ε
)
=
H
n
(
2
)
′
(
ε
)
/
H
n
(
2
)
(
ε
) and, therefore, this result is independent of the boundary conditions prescribed on the surface of the circular cylinder Σ
1. Hence, the accuracy of the method for the cylinder problems will depend on the accuracy of the approximation in H
n
(
2
)
′
(
ε
)
/
H
n
(
2
)
(
ε
).Ifd
~
n is eliminated from (28) and (29) as in the exact solution, the following relation is obtained:
(32)
c
~
n
=
ζ
n
η
n
b
~
n
,
where ζ
n and η
n are given, respectively, by (19) and (20). The system of equations for b
~
n
(33)
∑
l
=
-
∞
∞
‍
A
~
l
(
n
)
b
~
l
=
Ψ
~
n
is obtained by (32) and (30). Here A
~
l
(
n
) is
(34)
A
~
l
n
=
Λ
~
l
n
+
Φ
~
l
n
ζ
l
η
l
.
First b
~
n is obtained as a result of this system and then c
~
n and d
~
n are obtained, respectively, from (32) and (28) or (29). After having calculated b
~
n and c
~
n from the above relations if at r
=
r
1, namely, over Σ
1, relations u
t
(
r
1
,
θ
)
=
v
2
(
θ
), ∂
u
t
/
∂
r
(
a
,
θ
)
=
k
2
ξ
1
w
2
(
θ
), and X
(
x
,
y
)
=
-
1
/
4
i
H
0
(
2
)
(
k
1
|
x
-
y
|
) are used in the integral representation below:
(35)
u
i
x
+
∫
Σ
1
∂
u
t
y
∂
n
y
X
x
,
y
p
p
p
p
p
p
p
p
p
-
u
t
y
∂
∂
n
y
X
x
,
y
d
s
y
=
u
t
x
,
x
∈
B
1
.
The scattered field in any point in the region B
1 is obtained by calculating the following integral [4]:
(36)
u
1
r
,
θ
=
i
4
∫
0
2
π
v
2
(
θ
′
)
∂
∂
r
1
H
0
p
p
p
p
p
p
·
k
1
r
2
+
r
1
2
-
2
r
r
1
cos
⁡
θ
-
θ
′
1
/
2
p
p
p
p
p
p
-
k
2
ξ
1
w
2
(
θ
′
)
H
0
p
p
p
p
p
p
∂
∂
r
1
·
k
1
r
2
+
r
1
2
-
2
r
r
1
cos
⁡
θ
-
θ
′
1
/
2
r
1
d
θ
′
.
Thus the amplitude of the far field is obtained as follows:
(37)
P
θ
=
i
ε
4
∫
0
2
π
v
2
(
θ
′
)
i
cos
⁡
(
θ
-
θ
′
)
k
2
k
1
p
p
p
p
p
p
p
-
k
2
k
1
ξ
1
w
2
θ
′
e
i
ε
cos
⁡
⁡
θ
-
θ
′
d
θ
′
.Here if (24) is used and the expressions ∫
0
2
π
e
i
ε
cos
⁡
(
θ
-
θ
′
)
+
i
n
(
θ
′
-
π
/
2
)
d
θ
′
=
2
π
J
n
(
ε
)
e
i
n
θ and ∫
0
2
π
cos
⁡
(
θ
-
θ
′
)
e
i
ε
cos
⁡
(
θ
-
θ
′
)
+
i
n
(
θ
′
-
π
/
2
)
d
θ
′
=
-
i
2
π
J
n
′
(
ε
)
e
i
n
θ are considered, after the calculations, the amplitude of the far field is obtained as follows:
(38)
P
θ
=
i
ε
π
2
∑
n
=
-
∞
∞
‍
f
~
n
e
i
n
θ
.
Here again f
~
n is expressed as follows: f
~
n
=
J
n
′
(
k
1
r
1
)
J
n
(
k
2
r
1
)
-
k
2
ξ
1
/
k
1
J
n
(
k
1
r
1
)
J
n
′
(
k
2
r
1
)
b
~
n
+
J
n
′
(
k
1
r
1
)
H
n
(
k
2
r
1
)
-
k
2
ξ
1
/
k
1
J
n
(
k
1
r
1
)
H
n
′
(
k
2
r
1
)
c
~
n. For the exact solution, the amplitude of the far field is calculated from P
e
(
θ
)
=
∑
n
=
-
∞
∞
‍
a
n
e
i
n
θ. If b
~
n and c
~
n are replaced by b
n and c
n, f
~
n
=
J
n
′
(
k
1
r
1
)
J
n
(
k
2
r
1
)
-
k
2
ξ
1
/
k
1
J
n
(
k
1
r
1
)
J
n
′
(
k
2
r
1
)
b
n
+
J
n
′
(
k
1
r
1
)
H
n
(
k
2
r
1
)
-
k
2
ξ
1
/
k
1
J
n
(
k
1
r
1
)
H
n
′
(
k
2
r
1
)
c
n is obtained. Therefore i
ε
π
/
2
f
~
n
=
a
n is obtained.


## 4. Comparison

Comparisons are made between the exact answer of the problem and the surface radiation condition solutions. It is observed that the introduction of higher-order radiation conditions improves the approximation considerably in comparison with the results obtained by the use of a second-order radiation condition, especially in cases where creeping waves are less pervasive.Various graphics are generated from analytic solutions obtained from the on surface radiation condition method and from the exact solution of the problem in question in order to examine the consequences of using radiation boundary conditions having an order higher than two in the on surface radiation condition method. For simplicity and readability of graphics only the curves corresponding to solutions of the second and fourth order are shown. Curves corresponding to the third-order condition behave in between those two. In the figures the variations of the modulus of far field, namely, of scattering functionP with respect to θ, are given. These studies are done depending on the density of the medium and the velocity of the wave.From the comparisons made between the exact solution of the problem and the approximate analytic solution obtained by on surface radiation condition method, it can be observed that using the fourth order instead of the second order highly ameliorates the approximation of the method.(A) Whenever the regions are, respectively, air, water, and a region denser than water from outside to inside the results of the second- and fourth-order surface radiation conditions for the following parameters are depicted together with the exact curve forr
1
=
10,  r
2
=
5, ε
=
10,  k
3
=
c
1
/
c
3
k
1,  k
2
=
c
1
/
c
2
k
1,  ξ
1
=
ρ
1
/
ρ
2,  ξ
2
=
ρ
2
/
ρ
3,  ρ
1
=
1.2,  c
1
=
340,  ρ
2
=
1000,  c
2
=
1480,  ρ
3
=
1200, and c
3
=
1600. As for the graphics of the solutions concerning the far fields and under the parameter values, the variation of the scattered field with e
=
0 and ϕ
0
=
π is given in Figure 2(a)(a)
(b)
(c)Here the attitude of the curve is like in the case of a hard cylinder. Withe
=
0.5 and ϕ
0
=
π
/
3 the scattered field gives acceptable results especially for the fourth order as it is seen in Figure 2(b). As it is seen in Figure 2(c) the scattered field gives good results in shadowed and luminous regions especially for the fourth order when the eccentricity e
=
2.0 and ϕ
0
=
π
/
3 are taken.(B) Whenever the regions are, respectively, water, air, and a region denser than air from outside to inside the results of the second- and fourth-order surface radiation conditions for the parameters arer
1
=
10,  r
2
=
5,  ε
=
10,  k
3
=
c
1
/
c
3
k
1,  k
2
=
c
1
/
c
2
k
1,  ξ
1
=
ρ
1
/
ρ
2,  ξ
2
=
ρ
2
/
ρ
3,  ρ
1
=
1000,  c
1
=
1480,  ρ
2
=
1.2,  c
2
=
340,  ρ
3
=
1.0,  c
3
=
280. As for the graphics of the solutions concerning the far fields and under the parameter values, the variation of the scattered field with e
=
0 and ϕ
0
=
π is given in Figure 3(a).(a)
(b)
(c)Here the attitude of the curve is like in the case of a soft cylinder. In Figures3(b) and 3(c) the results of the scattered field in shadowed and luminous regions are compared as the eccentricity increases. With e
=
0.5 and ϕ
0
=
π
/
3 the scattered field gives good results especially for the fourth order as it is seen in Figure 3(b).As it is seen in Figure3(c) the scattered field gives good results in shadowed and luminous regions especially for the fourth order when the eccentricity e
=
2.0 and ϕ
0
=
π
/
3 are taken.(C) Whenever the regions are, respectively, a region denser than water, water, and air from outside to inside the results of the second- and fourth-order surface radiation conditions for the parameters are,r
1
=
10,  r
2
=
5,  ε
=
10,  k
3
=
c
1
/
c
3
k
1,  k
2
=
c
1
/
c
2
k
1,  ξ
1
=
ρ
1
/
ρ
2,  ξ
2
=
ρ
2
/
ρ
3,  ρ
1
=
1200,  c
1
=
1600,  ρ
2
=
1000,  c
2
=
1480,  ρ
3
=
1.2,  c
3
=
340.As for the graphics of the solutions concerning the far fields and under the parameter values, the variation of the scattered field withe
=
0 and ϕ
0
=
π is given in Figure 4(a). Here also the attitude of the curve is almost like in the case of a soft cylinder.(a)
(b)
(c)Withe
=
0.5 and ϕ
0
=
π
/
3 the scattered field gives good results for all orders in both shadowed and luminous regions as it is seen in Figure 4(b).As it is seen in Figure4(c) the scattered field gives good results everywhere for the fourth order, when the eccentricity e
=
2.0 and ϕ
0
=
π
/
3 are taken.



---

*Source: 282601-2015-02-18.xml*