Open access

Fractional Operators Approach and Fractional Boundary Conditions

Written By

Eldar Veliev, Turab Ahmedov, Maksym Ivakhnychenko

Submitted: October 8th, 2010 Published: June 21st, 2011

DOI: 10.5772/16300

Chapter metrics overview

2,850 Chapter Downloads

View Full Metrics

1. Introduction

Tools of fractional calculus including fractional operators and transforms have been utilized in physics by many authors (Hilfer, 2000). Fractional operators defined as fractionalizations of some commonly used operators allow describing of intermediate states. For example, fractional derivatives and integrals (Oldham & Spanier, 1974; Samko et al., 1993) are generalizations of derivative and integral. Fractional curl operator defined in (Engheta, 1998) is a fractionalized analogue of conventional curl operator used in many equations of mathematical physics. A fractionalized operator generalizes the original operator. The idea to use fractional operators in electromagnetic problems was formulated by N. Engheta (Engheta, 2000) and named “fractional paradigm in electromagnetic theory”.

Our purpose is to find possible applications of the use of fractional operators in the problems of electromagnetic wave diffraction. In this paper two-dimensional problems of diffraction by infinitely thin surfaces are considered: a strip, a half-plane and a strip resonator (Fig.1). Assume that an incident field is an E-polarized plane wave, described by the function


Figure 1.

Geometry of the diffraction problems: a) strip, b) half-plane, c) two parallel strips.

where θis the incidence angle, k=2πλis the wavenumber. Here, the time dependence is assumed to be eiωtand omitted throughout the paper. There are three structures considered in this paper:

  1. a strip located in the plane y=0(x[a,a]) infinite along the axis z (Fig. 1a);

  2. a half-plane (y=0,x0) (Fig. 1b);

  3. two parallel strips infinite along the axis z (a strip resonator). The first strip is located aty=l, x[a,a], and the second one is aty=l, x[a,a](Fig. 1c).

One may ask what new features are that the fractional operators can bring to the theory of diffraction. The concept of intermediate states, obtained with the aid of fractional

derivatives and integrals, yields to various generalizations of commonly used models in electrodynamics such as:

  1. Intermediate waves. For instance, intermediate waves between plane and cylindrical waves (Engheta, 1996, 1999) can be obtained using fractional integral of scalar Green’s function:


where G2is two-dimensional Green’s function of the free space. Gαdescribes an intermediate case between one- and two-dimensional Green’s functions and have the following behavior in the far-zone (Engheta, 1999):

Gα~i4πcos(πα2)(ksin|φ|)α2πkρeikρiπ/4+i2kαΓ(α)eik|x|k|y|1α, kρ=kx2+y2,

This function consists of two waves: a cylindrical wave and a non-uniform plane wave propagating in the xdirection and behaving with yas|y|α1.

  1. Fractional Green’s function Gαdefined as a fractional derivative (integral) of the ordinary Green’s function of the free space -GαDkyαG. αdenotes the fractional order and varies from 0 to 1 (0α1). In two-dimensional case Gαis expressed as

  1. Fractional Green’s theorem which involves fractional derivatives of ordinary Green’s function and fractional derivatives of the considered function on a boundary of a domain (Veliev & Engheta, 2003). The corresponding equations will be presented later in this paper.

  2. Fractional boundary conditions (FBC) defined via fractional derivatives of the tangential electric field componentsU(x,y). For an infinitely thin boundary Slocated in the planey=d, FBC are defined asDyαU(x,y)|yS=0,y±d.

The order of the fractional derivative αis assumed to be between 0 and 1. Fractional derivative Dαis applied along the direction normal to the surfaceS. Fractional boundary conditions describe an intermediate boundary between the perfect electric conductor (PEC) and the perfect magnetic conductor (PMC), obtained from FBC if the fractional order equals to 0 and 1, respectively.

We will use the symbol Dyαfto denote operator of fractional derivative or integralDyαf, which is defined by the integral of Riemann-Liouville on semi-infinite interval (Samko et al., 1993):(Dxαf)(x)=1Γ(1α)ddxxf(t)dt(xt)α,0<α<1.

where Γ(1α)is Gamma function.

This paper is devoted to the problems of diffraction by a strip, a strip resonator and a half-plane characterized with fractional boundary conditions with 0α1expressed asDkyαEz(x,y)=0, y±0,xL.

where L=(a,a)for a strip and L=(0,)for a half-plane. For convenience, fractional derivative is applied with respect to dimensionless variableky. The function Ez(x,y)denotes z-component of the total electric field, Ez(x,y)=Ezi+Ezs, that is the sum of the incident plane wave Ezi(x,y)and the scattered waveEzs(x,y).

In case of a strip resonator we have two equations to impose fractional boundary conditions:DkyαEz(x,y)=0, yl±0, x(a,a),DkyαEz(x,y)=0, yl±0,x(a,a).

From the one hand, introduction of new boundary conditions should describe a new physical boundary world, and from the other hand they must allow to build an effective computational algorithm to solve the stated problems with a desired accuracy. Simple mathematical description of the scattering properties of surfaces is a common problem in modeling in diffraction theory.

One of the well-studied boundaries, which can be treated as an intermediate state between PEC and PMC, is an impedance boundary defined by the equationn×E(r)=ηn×(n×H(r)),rS.

where nis the normal to the surfaceS. The value of the impedance ηvaries from 0 for PEC to ifor PMC.

There are many papers devoted to diffraction by impedance boundaries. Impedance boundary conditions (IBC) have been used for the modeling of the scattering properties of good conductors, gratings, etc. In each case there are formulas to define the value of the impedance as a function of material parameters. IBC are approximate BC and therefore they have limitations in usage and cannot describe all diversity of boundaries.

Further approximation of IBC can be made with the aid of derivatives of higher but integer orders or generalized boundary conditions (Hope & Rahmat-Samii, 1995; Senior & Volakis, 1995). A general methodology to obtain exact IBC of higher order in spectral domain is presented in (Hope & Rahmat-Samii, 1995), where flat covers (and also surfaces with curvature) consisting of homogeneous materials with an arbitrary (linear, bi-anisotropic) constitutive equations. It is possible to obtain exact IBC in the spectral domain that can be often done in an analytical form very often. However, it is not always possible to get IBC in the spatial domain in an exact form. That is why it is necessary to approximate IBC in the spectral domain in order to apply inverse Fourier transform.

Another boundary condition that generalizes the perfect boundaries like PEC and PMC was introduced in (Lindell & Sihvola, 2005a). The corresponding surface was named perfect electromagnetic conductor (PEMC) and the mentioned condition is defined asH+ME=0.

ForM=0, PEMC defines a PEC boundary and for M=we get a PMC. The physical model of PEMC boundary was proposed in (Lindell & Sihvola, 2005b) where it was shown that the PEMC condition can simulate reflection from an anisotropic layer for the normal incidence of the plane wave. Diffraction by a PEMC boundary has not been considered yet. Further generalization of PEMC can be made using concept of the generalized soft-and-hard surface (GSHS) (Haninnen et al., 2006):aE=0,bH=0.

wherea, bare complex vectors that satisfy equations na=nb=0andab=1. GSHS can transform an incident plane wave with any given polarization into any other polarization of the reflected plane wave if the vectorsa, bare chosen appropriately (Haninnen et al., 2006).

Fractional boundary conditions (FBC) can be compared with impedance boundary conditions (IBC). First of all FBC are intermediate between PEC and PMC as well as IBC. The value of fractional order α=0(α=1) corresponds to the value of impedance η=0(η=i), respectively. For other values of 0<α<1the deeper analysis is needed.

Physical analysis of the strip with FBC shows that the induced surface currents behave similarly to the currents on an impedance strip. Due to specific properties the strip with FBC is compared with the well-known impedance strip. It can be shown that for a wide range of input parameters the “fractional strip” behaves similarly to the impedance strip if the fractional order is chosen appropriately (Veliev et al., 2008b). The proposed method used for a “fractional strip” has some advantages over the known methods applied to the analysis of the wave scattering by an impedance strip.

The purpose of this work is to build an effective analytic-numerical method to solve two-dimensional diffraction problems for the boundaries described by fractional boundary conditions with α ∈ [0,1]. The method will be applied to two canonical scattering objects: a strip and a half plane. The method is based on presenting the scattered field via fractional Green’s function,Ezs(x,y)Lf1α(x')Gα(xx',y)dx',

where f1α(x)is the unknown function and Gα(xx',y)=i4DkyαH0(1)(k(xx')2+y2)is the fractional derivative of the Green’s function defined by equation (2). This presentation leads to the following dual integral equations (DIE) with respect to the Fourier transform F1α(q)=Lf1α(ξ)eikqξdξof the functionf1α(x)

{F1α(q)eikdLξq(1q2)α1/2dq=4πeiπ/2(1α)sinαθeikdLξcosθ, ξL,F1α(q)eikdLξqdq=0,             ξL,E5

where dL=aforL=(a,a), dL=1forL=(0,).

In the case of a strip resonator, we obtain more complicated set of integral equations which will be presented later in this paper.

The method generalizes the known method used for the PEC and PMC strip and half plane. As will be shown later, this method allows obtaining a solution for the value α=0.5in the explicit analytical form. For other values of α[0,1]the scattering problems are reduced to solving of the infinite systems of linear algebraic equations (SLAE). In order to discretize the DIE the function f1α(x)is represented as a series in terms of orthogonal polynomials: Gegenbauer polynomials for the strip and Laguerre polynomials for the half-plane. These representations result in a special kind of the edge conditions for the fractional current density functionf1α(x). The physical characteristics of the considered scattering objects can be found with any desired accuracy by solving SLAE.


2. Diffraction by a strip with fractional boundary conditions

Assume that an E-polarized plane wave is characterized with the functionEi=zEzi(x,y)=zeik(xcosθ+ysinθ). The total field E=zEz(x,y)must satisfy fractional boundary conditions


where L=(a,a)for a strip. For convenience, fractional derivative Dkyαis applied with respect to a dimensionless variableky. The function Ez(x,y)denotes the z-component of the total electric field Ez(x,y)=Ezi+Ezsthat is the sum of the incident plane wavе Ezi(x,y)and the scattered fieldEzs(x,y). Solution to the diffraction by the screen S={(x,y):y=0,a<x<a}is to be sought under the following conditions:

  1. The total field Emust satisfy the Helmholtz equation everywhere outside the screen

  1. The scattered field Ezs(x,y)must satisfy Sommerfeld radiation condition at the infinity

  1. The total field Emust satisfy the edge condition, i.e. the finiteness of energy in every local area near the edges of the screen (Honl et al., 1961).

  2. The total field Ez(x,y)must satisfy the boundary conditions (3).

The method is based on representation of the scattered field with the aid of the fractional derivative of the Green’s function:


In (6), the function f1α(x)is the unknown function called the density of the fractional potential, and Gαis the fractional derivative of two-dimensional the Green’s function of the free space defined by equation (2).

For the limit cases of the fractional order with α=0and α=1representation (6) corresponds to the single-layer and double-layer potentials commonly used to present the scattered fields in diffraction problems:

Ezs(x,y)={i4aaf'(x')H0(1)(k(xx')2+y2)dx',          α=0i4aaf(x')yH0(1)(k(xx')2+y2)dx',        α=1

More general representations (6) can be derived from the fractional Green’s theorem (Veliev & Engheta, 2003) which generalizes the ordinary Green’s theorem.

2.1. Fractional Green’s theorem

Consider a functionψ(r), which satisfies inhomogeneous scalar Helmholtz equation with the source density given by the functionρ(r):


Besides, defineG(r,r0)as the Green’s function of the Helmholtz equation:


Here, δ(rr0)is the three-dimensional Dirac delta function, rand r0are the position vectors for the observation and source points, respectively, Δ=2x2+2y2+2z2is the Laplacian, and kis a scalar constant. After applying fractional derivatives to equations (7) and (8) with respect to the xvariable, multiplying the first equation withDxνG(r,r0), and the second withDxμψ(r), subtracting one from another, integrating this over all source coordinates x0,y0,z0insideS, and finally using the Green’s theorem, we obtain the following representation:

Dxβψ(r)={VDx0βνρ(r0)·Dx0νG(r,r0)dv0++14πS[Dx0νG(r,r0)·0Dx0βνψ(r0)Dx0βνψ(r0)·0Dx0νG(r,r0)]·ds0,  rV 0,          rVE12

whereμ+ν=β. Operator 0denotes the operator of gradient in respect of variabler0(x0,y0,z0). Here it was used the property of the fractional derivative of the Dirac delta function:


We use the uniform symbol Dxα(orDxα) to denote both fractional derivatives and fractional integrals, and it defines a fractional derivative for 0<α<1and a fractional integral forα<0.

Equation (9) is a generalization of well-known Green’s theorem for the case of fractional derivatives.

Consider some important particular cases, which can be obtained from (9).

In the case of excitation in a free space so that the volume Vis the whole space, the surface integrals in (9) vanish, and we have:


Originally function ψ(r)characterizes the field excited by the source with the volume densityρ(r). From the other hand, for β=0representation (11) means that the field ψ(r)is expressed through the distribution of fractional sources with density Dνρ(r0)inside the volume Vand by using fractional integral of conventional Green’s functionDνG(r0,r).

Assumingρ(r)=0, we can obtain some other important representations:

Dxβψ(r)={14πS[Dx0βG(r,r0)·0ψ(r0)ψ(r0)·0Dx0βG(r,r0)]ds0,   if  ν=β,μ=014πS[G(r,r0)·0Dx0βψ(r0)Dx0βψ(r0)·0G(r,r0)]ds0,   if  ν=0 E15

From this representation we see that the fractional derivative of function ψ(r)is expressed either via the value of the function and its first derivative at the boundary and the fractional derivatives of Green’s function, or by the fractional derivatives of the function at the boundary and the usual Green’s function.

Ifν=μ, i.e.β=0, we obtain a representation for the function ψ(r)itself:


This expression means that the function ψ(r)is represented through its fractional derivatives at the boundary and the fractional derivatives of Green's function. The equation (13) can be useful in scattering problems. If we have boundary conditions for the function ψ(r)on the surface Sas 0Dx0μψ(r0)|r0S=0(orDx0μψ(r0)|r0S=0) then one of the surface integrals in (13) vanishes and we get a simple presentation forψ(r). This fact will be used to present the scattered field in all diffraction problems considered in this paper (6). Equations (12), (13) generalize the Huygens principle in such a sense that the fractional derivative of the functionψ(r), which characterizes a wave process, is presented as a superposition of waves radiated by elementary "fractional" sources distributed on the given surface. “Fractional” potentials, SDx0βνψ(r0)·0Dx0νG(r,r0)·ds0, SDx0νG(r,r0)·0Dx0βνψ(r0)·ds0, can be treated as a generalization of well-known single and double layer potentials.

2.2. Solution to integral equations

Substituting the expression (6) for Ez(x,y)into fractional boundary conditions (3) we get the equation

limy0 DkyαLf1α(x')Gα(xx',y)dx'=limy0DkyαEzi(x,y)E17

It is convenient to use the Fourier transform of the fractional potential density f1α(x)


where a new function f˜1α(ξ)is introduced:


Then the scattered field is expressed via the Fourier transform F1α(q)as


where the upper (lower) sign is chosen for y>0 (y<0). Here, in (15), the following representation for the fractional Green’s function was used:


It can be shown that the equation (14) can be reduced to dual integral equations (DIE)

{F1α(q)eikaξq(1q2)α1/2dq=4πeiπ/2(1α)sinαθeikdLξcosθ, |ξ|<1,F1α(q)eikaξqdq=0,             |ξ|>1,E23

For the limit cases of the fractional order α =0 and α =1 the equations (17) are reduced to the well known integral equations used for PEC and PMC strips (Honl et al., 1961; Veliev & Veremey, 1993; Veliev & Shestopalov, 1988; Uflyand, 1977), respectively. In this paper the method to solve DIE (17) is proposed for arbitrary value of α∈[0,1].

DIE (17) can be solved analytically for one special case of α= 0.5. In this case we get the solutions for any value of kas


In the case of arbitrary αthe solutions can be obtained numerically. First, we modify the equations (17). After multiplying by eikaτξand integrating in ξfrom -1 to 1, the first equation in (17) can be rewritten in the following form:


In order to discretize this equation, we present the unknown function f˜1α(ξ)as a uniformly convergent series in terms of the orthogonal polynomials with corresponding weight functions which allow satisfying the edge conditions:


where Cnα(x)are the Gegenbauer polynomials and fnαare the unknown coefficients. Gegenbauer polynomials can be treated as intermediate polynomials between Chebyshev polynomials of the first and second kind:


The Fourier transform F1α(q)is expressed as the series


where Jn+α(kaq)is the Bessel function.

It must be noted that the edge conditions are chosen in the following form


For special cases of α=0and α=1the edge conditions have the form as

f˜1α(ξ)={O((1ξ2)1/2),     α=0O((1ξ2)1/2),      α=1,ξ±1E31

These are well-known Meixner edge conditions in diffraction problems (Honl et al., 1961).

Substituting (22) into (17) and taking into account the properties of discontinuous integrals of Weber-Shafheitlin (Bateman & Erdelyi, 1953) and the following formula (Prudnikov et al., 1986)


one can show that the homogenous equation in the set (17) is satisfied identically.

The first equation of (17) written in the form (20) can be reduced to an infinite system of linear algebraic equations (SLAE) with respect to the unknown coefficientsfnα:


where the matrix coefficients are expressed as


It can be shown that the SLAE (26) can be reduced to SLAE of the Fredholm type of the second kind (Veliev et al., 2008a). Then the coefficients fnαcan be found with any desired accuracy (within the machine precision) using the truncation of SLAE. The fractional density f1α(x)is computed by using (21) and the scattered field (6) and other physical characteristics can be obtained as series in terms of the found coefficientsfnα.

In order to solve the diffraction problem on a plane screen with fractional boundary conditions and obtain a convenient SLAE we applied several techniques. First of all, the fractional Green’s theorem presented above allowed searching the unknown scattered field as a potential with the fractional Green’s function. The order of the fractional Green’s function is defined from the fractional order of the boundary conditions. In general, the fractional derivative of Green’s function may have a complicated form, but we used the Fourier transform where application of the fractional derivative maps to a simple multiplication by(iq)α. Finally, utilization of the orthogonal Gegenbauer polynomials along with the specific form of the edge conditions allowed to reduce integral equations to SLAE in a convenient form. One can compare the method presented for fractional boundary conditions with the known methods applied to solve diffraction by an impedance strip. The impedance strip requires to consider two unknown densities in presentation of the scattered field as a sum of single- and double-layer potentials. The usage of two unknown functions leads to more complicated SLAE in spite of the SLAE obtained for fractional boundary conditions.

2.3. Physical characteristics

We consider such electrodynamic characteristics of the scattered field as the radiation pattern (RP), monostatic radar cross-section (MRCS) and surface current densities depending on the coefficientsfnα. The scattered field Ezs(x,y)in the far-zone krin the cylindrical coordinate system(r,ϕ), x=rcosφ,y=rsinφ, is expressed as


where the upper sign is chosen forφ[0,π], and the lower one whenφ[π,2π]. Using the stationary phase method for krwe present Ezs(x,y)as




The function Φα(φ)describes RP and can be expressed via the coefficients fnαas


In physical optics (PO) approximation (ka1) Φα(φ)has a simpler form. Using the following formula


in IE (20) we get the following expressions for Fα(q)andΦα(φ):


In the special case of α=0.5and arbitrary value of kawe get an analytical expression for the RP


Bi-static radar cross section (BRCS) is expressed from RP Φ(φ)asσ2dλ(φ)=2π|Φ(φ)|2. MRCS σ2Dmonois defined asσ2Dmono=σ2dλ(θ)=2π|Φ(θ)|2.

We have the following representations in PO approximation


It must be noted that the density function f1α(x)in the integral (6) does not describe the density of physical surface currents on the strip for0<α<1. The function f1α(x)is defined as the discontinuity of fractional derivatives of E-field at the planey=0:

f1α(x)= Dky1αEz(x,y)|y=+0Dky1αEz(x,y)|y=0,x(a,a)E46

For the limit cases of α=0and α=1the equation (29) is reduced to well-known presentations for electric and magnetic surface currents, respectively, i.e.

f1α(x)={Ez(x,y)y|y=+0Ez(x,y)y|y=0=Hx(x,+0)Hx(x,0),    α=0Ez(x,+0)Ez(x,0),                    α=1

In order to obtain physical surface currents from f1α(x)we have to apply additional integration. In case of E-polarized incident plane wave we have the following induced currents on a strip: electric current jα(e)=zjzα(e)and magnetic current jα(m)=xjxα(m)expressed from f1α(x)as


The detailed analysis of the scattering properties of the strip with fractional boundary conditions one can find in papers (Veliev et al., 2008a; Veliev et al., 2008b).

2.4. H-polarization

In the case of the H-polarized incident plane waveHi(0,0,Hzi), whereHzi(x,y)=eik(xcosθ+ysinθ), the method proposed above can be applied as well. We define fractional boundary conditions as


The case of α=0corresponds to diffraction of the H-polarized plane wave on a PEC strip, while the case of α=1describes diffraction of the H-polarized plane wave on a PMC strip. As before, we represent the scattered field via the fractional Green’s function


After substituting (18) into fractional boundary conditions (19) we get the equation


This equation can be solved by repeating all steps of the E-polarization case after changing αto1α.


3. Diffraction by a half-plane with fractional boundary conditions

Another problem studied in this paper is the diffraction by a half-plane with fractional boundary conditions. The method introduced to solve the dual integral equation (DIE) for a finite object (a strip) will be modified to solve DIE for semi-infinite scatterers such as half-plane. There are many papers devoted to the classical problem of diffraction by a half-plane. The method to solve the scattering problem for a perfectly conducting half-plane is presented in (Honl et al., 1961). Usually, it is solved using Wiener-Hopf method. The first application of the method to a PEC half-plane can be referred to the papers of Copson (Copson, 1946) and independently to papers of Carlson and Heins (Carlson & Heins, 1947). In 1952 Senior first applied Wiener-Hopf method to the diffraction by an impedance half-plane (Senior, 1952) and later oblique incidence was considered (Senior, 1959). Diffraction by a resistive and conductive half-plane and also by various types of junctions is analyzed in details in (Senior & Volakis, 1995). We propose a new approach for the rigorous analysis of the considered problem which generalizes the results of (Veliev, 1999) obtained for the PEC boundaries and includes them as special cases.

Let an E-polarized plane wave Ezi(x,y)=eik(xcosθ+ysinθ)(1) be scattered by a half-plane (y=0,x>0). The total field Ez=Ezi+Ezsmust satisfy fractional boundary conditions


and Meixner’s edge conditions must be satisfied forx0.

Following the idea used for the analysis of diffraction by a strip we represent the scattered field using the fractional Green’s function


where f1α(x)is the unknown function, Gαis the fractional Green’s function (2).

After substituting the representation (31) into fractional boundary conditions (30) we get the equation


The Fourier transform of f1α(x)is defined as


where f˜1α(ξ)f1α(ξ)for ξ>0and f˜1α(ξ)0forξ<0.

Then the scattered field will be expressed via the Fourier transform F1α(q)as


Using the Fourier transform the equation (32) is reduced to the DIE with respect toF1α(q):


The kernels in integrals (34) are similar to the ones in DIE (17) obtained for a strip if the constant dLis equal to 1 (L=(0,)in the case of a half-plane).

For the limit cases of the fractional order α=0and α=1these equations are reduced to well known integral equations used for the PEC and PMC half-planes (Veliev, 1999), respectively. In this paper the method to solve DIE (5) is proposed for arbitrary values ofα[0,1].

DIE allows an analytical solution in the special case of α=0.5in the same manner as for a strip with fractional boundary conditions. Indeed, for α=0.5we obtain the solution for any value of kas


The scattered field can be found in the following form:

Ezs(x,y)=i2ke±iπα/2eiπ/4sinα1/2θeik(cosθx+|y|sinθ),α=0.5, for
y>0   (y<0)E60

In the general case of 0<α<1the equations (34) can be reduced to SLAE. To do this we represent the unknown function f˜1α(ξ)as a series in terms of the Laguerre polynomials with coefficientsfnα:


Laguerre polynomials are orthogonal polynomials on the interval L=(0,)with the appropriate weight functions used in (35). It can be shown from (35) that f˜1α(ξ)satisfies the following edge condition:


For the special cases of α =0 and α =1, the edge conditions are reduced to the well-known equations (Honl et al., 1961) used for a perfectly conducting half-plane.

After substituting (35) into the first equation of (34) we get an integral equation (IE)


where R(ξ)=4πeiπ/2(1α)sinαθeikξcosθis known.

Using the representation for Fourier transform of Laguerre polynomials (Prudnikov et al., 1986) we can evaluate the integral over dtas


After some transformations IE (37) is reduced to


Then we integrate both sides of equation (38) with appropriate weight functions, as0()eξξα1/2Lmα1/2(2ξ)dξ. Using orthogonality of Laguerre polynomials we get the following SLAE:


with matrix coefficients


It can be shown that the coefficients fnαcan be found with any desired accuracy by using the truncation of SLAE. Then the function f˜1α(x)is found from (35) that allows obtaining the scattered field (33).


4. Diffraction by two parallel strips with fractional boundary conditions

The proposed method to solve diffraction problems on surfaces described by fractional boundary conditions can be applied to more complicated structures. The interest to such structures is related to the resonance properties of scattering if the distance between the strips varies. Two strips of the width 2ainfinite along the axis z are located in the planes y=landy=l. Let the E-polarized plane wave Ezi(x,y)=eik(xcosθ+ysinθ)(1) be the incident field. The total field Ez=Ezi+Ezssatisfies fractional boundary conditions on each strip:


and Meixner’s edge conditions must be satisfied on the edges of both strips (y=±l,x±a).

The scattered field Ezs(x,y)consists of two parts




Here, Gαis the fractional Green’s function defined in (2). y1,2 are the coordinates in the corresponding coordinate systems related to each strip,


Using Fourier transforms, defined as


the scattered field is expressed as


Fractional boundary conditions (30) correspond to two equations


After substituting expressions (41) and (42) into the equations (43) and (44) we obtain


Multiplying both equations with e–ikxτand integrating them in ζon the interval [–a,a], the system (45), (46) leads to


Similarly to the method described for the diffraction by one strip, the set (47) can be reduced to a SLAE by presenting the unknown functions fj1α(x)as a series in terms of the orthogonal polynomials. We represent the unknown functions f˜j1α(ξ)as series in terms of the Gegenbauer polynomials:


For the Fourier transforms Fj1α(q)we have the representations (22). Substituting the representations for Fj1α(q)into the (47), using the formula (25), then integrating (.)Jm+α(kaτ)mαdτform=0,1,2,.., we obtain the following SLAE:


where the matrix coefficients are defined as


Consider the case of the physical optics approximation, whereka1. In this case we can obtain the solution of (47) in the explicit form. Indeed, using the formula (28) we get


Finally, we obtain the solution as


Having expressions for Fj1α(q)we can obtain the physical characteristics. The radiation pattern of the scattered field in the far zone (27) is expressed as




5. Conclusion

The problems of diffraction by flat screens characterized by the fractional boundary conditions have been considered. Fractional boundary conditions involve fractional derivative of tangential field components. The order of fractional derivative is chosen between 0 and 1. Fractional boundary conditions can be treated as intermediate case between well known boundary conditions for the perfect electric conductor (PEC) and perfect magnetic conductor (PMC). A method to solve two-dimensional problems of scattering of the E-polarized plane wave by a strip and a half-plane with fractional boundary conditions has been proposed. The considered problems have been reduced to dual integral equations discretized using orthogonal polynomials. The method allowed obtaining the physical characteristics with a desired accuracy. One important feature of the considered integral equations has been noted: these equations can be solved analytically for one special value of the fractional order equal to 0.5 for any value of frequency. In that case the solution to diffraction problem has an analytical form. The developed method has been also applied to the analysis of a more complicated structure: two parallel strips. Introducing of fractional derivative in boundary conditions and the developed method of solving such diffraction problems can be a promising technique in modeling of scattering properties of complicated surfaces when the order of fractional derivative is defined from physical parameters of a surface.


  1. 1. BatemanH.ErdelyiA.1953Higher Transcendental Functions,2McGraw-Hill, New York
  2. 2. Carlson J.F. & Heins A.E.1947The reflection of an electromagnetic plane wave by an infinite set of plates.Quart. Appl. Math.,4313329
  3. 3. Copson E.T.1946On an integral equation arising in the theory of diffraction,Quart. J. Math.,171934
  4. 4. EnghetaN.1996Use of Fractional Integration to Propose Some ‘Fractional’ Solutions for the Scalar Helmholtz Equation. A chapter in Progress in Electromagnetics Research (PIER), Monograph Series, Chapter 5,12Jin A. Kong, ed.EMW Pub., Cambridge, MA,107132
  5. 5. EnghetaN.1998Fractional curl operator in electromagnetic.Microwave and Optical Technology Letters,1728691
  6. 6. EnghetaN.1999Phase and amplitude of fractional-order intermediate wave,Microwave and optical technology letters,215
  7. 7. EnghetaN.2000Fractional Paradigm in Electromagnetic Theory,a chapter in IEEE Press,chapter 12,523553
  8. 8. HanninenI.LindellI. V.SihvolaA. H.2006Realization of Generalized Soft-and-Hard Boundary,Progress In Electromagnetics Research, PIER 64,317333
  9. 9. HilferR.1999Applications of Fractional Calculus in Physics, World Scientific Publishing,9-81023-457-0
  10. 10. HonlH.A.MaueW.WestpfahlK.1961Theorie der Beugung, Springer-Verlag, Berlin
  11. 11. HopeD. J.Rahmat-SamiiY.1995Impedance boundary conditions in electromagnetic, Taylor and Francis, Washington, USA
  12. 12. Lindell I.V. & Sihvola A.H.2005Transformation method for Problems Involving Perfect Electromagnetic Conductor (PEMC) Structures.IEEE Trans. Antennas Propag.,5330053011
  13. 13. Lindell I.V. & Sihvola A.H.2005Realization of the PEMC Boundary.IEEE Trans. Antennas Propag.,5330123018
  14. 14. OldhamK. B.SpanierJ.1974The Fractional Calculus: Integrations and Differentiations of Arbitrary Order, Academic Press, New York
  15. 15. PrudnikovH. P.BrychkovY. H.MarichevO. I.1986Special Functions,Integrals and Series,2Gordon and Breach Science Publishers
  16. 16. SamkoS. G.KilbasA. A.MarichevO. I.1993Fractional Integrals and Derivatives, Theory and Applications,Gordon and Breach Science Publ., Langhorne
  17. 17. SeniorT. B. A.1952Diffraction by a semi-infinite metallic sheet.Proc. Roy. Soc. London, Seria A, 213,436458.
  18. 18. SeniorT. B. A.1959Diffraction by an imperfectly conducting half plane at oblique incidence.Appl. Sci. Res., B8,3561
  19. 19. SeniorT. B.VolakisJ. L.1995Approximate Boundary Conditions in Electromagnetics, IEE, London
  20. 20. UflyandY. S.1977The method of dual equations in problems of mathematical physics [in russian].Nauka, Leningrad
  21. 21. VelievE. I.ShestopalovV. P.1988A general method of solving dual integral equations.Sov. Physics Dokl.,336411413
  22. 22. VelievE. I.VeremeyV. V.1993Numerical-analytical approach for the solution to the wave scattering by polygonal cylinders and flat strip structures.Analytical and Numerical Methods in Electromagnetic Wave Theory, M. Hashimoto, M. Idemen, and O. A. Tretyakov (eds.), Chap. 10, Science House, Tokyo
  23. 23. VelievE. I.1999Plane wave diffraction by a half-plane: a new analytical approach.Journal of electromagnetic waves and applications,131014391453
  24. 24. VelievE. I.EnghetaN.2003Generalization of Green’s Theorem with Fractional Differintegration, IEEE AP-S International Symposium & USNC/URSI National Radio Science Meeting
  25. 25. VelievE. I.IvakhnychenkoM. V.AhmedovT. M.2008Fractional boundary conditions in plane waves diffraction on a strip.Progress In Electromagnetics Research,79443462
  26. 26. VelievE. I.IvakhnychenkoM. V.AhmedovT. M.2008Scattering properties of the strip with fractional boundary conditions and comparison with the impedance strip.Progress In Electromagnetics Research C,2189205

Written By

Eldar Veliev, Turab Ahmedov, Maksym Ivakhnychenko

Submitted: October 8th, 2010 Published: June 21st, 2011