.
Our site needs your help!
Site categories

Understanding Parabolic Approximation in Wave Propagation: Analytical Methods and Applications

Join Our Telegram (Seaman Community)

Discover the intricacies of wave propagation through stratified and random media with our in-depth analysis of parabolic approximation methods. This article covers the application of the wave equation in turbulent tropospheres, highlighting the impact of refractive index fluctuations.

Learn about the Green function for parabolic equations in stratified environments and enhance your knowledge of wave behavior in complex settings.

Analytical Methods in the Problems of Wave Propagation in a Stratified and Random Medium

The analogy between the non-stationary problem of quantum mechanics and the parabolic approximation in a problem wave propagation was likely first explored by Fock. It is well known that the methods of classic wave theory were utilised in quantum mechanics in the early stages of its development. At present, the situation is rather the reverse, and the mathematical methods of quantum mechanics have become widely adopted in the latest developments of wave theory.

The asymptotic solution to the problem of the diffraction of a plane monochromatic wave over a sphere of large radius, compared with the wavelength of radiation, was obtained by Fock in 1945. In obtaining this solution, he developed an asymptotic theory of diffraction on the basis of Leontovich’s boundary conditions and the concept of the “local” field. The approach utilised several large parameters involved in problem formulation, such as |η| ≫ 1, m = (kα/2)1/3 ≫ 1, η = εg+j4πσω, εg, σ, the dielectric permittivity and conductivity of the earth, respectively, ω is the cyclic frequency of the radiated wave, k = 2p/λ is a wavenumber, λ is a wavelength, α is the radius of the earth’s sphere. The same results were obtained by Fock later by means of the parabolic approximation to the wave equation which is based on significantly different parameters, with m ≫ 1, and the characteristic scales of the wave oscillations in two different directions: along the earth’s surface and normal to it. With regards to radio wave propagation at a frequency above 1 GHz the important development was a solution of the parabolic equation for the stratified troposphere, where the refractive index depends solely on the height above the ground surface.

The analysis of the wave field at low altitudes z, which are small compared with the radius of the earth’s curvature α, is convenient to perform by means of the introduction of the modified refractivity:

nm = n(z)+z/α.

The solution obtained is suitable for the very general case of the dependence n(z) and impedance boundary conditions at the earth’s surface. The solution is represented by a contour integral which in the shadow region is calculated by a sum over residues in the poles of the integrand, i. e. a series of normal waves. In the line-of-sight (LOS) region the contour integral is commonly calculated using stationary phase methods similar to a geometric optic approximation. Further study of the wave propagation in a stratified tropo- sphere was dedicated to the development of the effective methods of the field representation for various height profiles of the refractivity and combinations of the relative location of the transmitting and receiving antennae. An alternative analytical method to the solution of wave propagation in stratified medium, the method of “invariant submergence”, was developed by Klatskin. Using that method, a boundary problem, which is, in fact, a problem of diffraction, can be transformed into a problem of evolution. A significant advantage of this method is its inherent capability to study the wave propagation through random media, and the method is especially effective when applied to a randomly stratified medium.

Together with a ducting mechanism, the wave scattering on random inhomogeneities of the refractive index plays a significant role in the study of radio wave propa- gation in the atmospheric boundary layer at frequencies above 1 GHz. The straight-forward application of either the contour integral or normal mode series to a problem of wave propagation in a random troposphere will face significant difficulties. The major problem will be related to the divergence of the matrix elements of the normal wave transformation due to a scattering on random fluctuations in the refractive index. The deployment of the method of “invariant submergence” also does not lead to an analytical solution to the problem though numerical calculation is possible.

In this category we use the representation of the Green function in the form of an expansion over the eigen functions of a continuous spectrum which allows one to avoid the above mentioned difficulties of divergence of matrix elements and then to separate the matrix coefficients with a “nearly discrete” and continuum spectrum with normalised height gain functions. This approach led to some advances in solving several problems, in particular, the problem of wave propagation in a tropospheric duct filled with random irregularities of the refractive index in addition to a regular supper-refractive gradient of the refractivity. The disadvantage of the method is that the separation of the discrete spectrum of eigen functions is performed in the “unperturbed”, i. e. regular, part of the Schrödinger operator and, therefore, the spectrum of eigen functions is left unchanged by the random component of the potential term. The study of the spectrum of Schrödinger’s equation with random potential is a complex problem by itself and most advances so far have been achieved in the one-dimensional problem.

The theory of wave scattering in a random medium has been developed during the last decades. The most significant progress has been achieved in the solution to the problems of wave propagation in either a random but unbound medium or a uniform medium with a random boundary interface. In the solution of these problems a statistical approach is commonly used. This approach is targeted at obtaining the statistical characteristics of the scattered field:

  • probability function or,
  • statistical moments of the field.

Among them the first two moments: average field (coherent component) and either coherence function or field intensity, are of significant interest in many practical applications. There are three known methods of obtaining the closed equations for the field moments: Feynman diagrams, local perturbation and the Markov process approximation.

The study of a coherent component can be done using the modified perturbation theory, a similar result is obtained from the diagram methods. The advances in the solution of the higher order moments of the field have been obtained after application of the parabolic approximation to the wave equation. Neglecting back scattering in the parabolic approximation allows one to employ the concept of the Markov process and all the following analogies from the statistical theory of Markov processes.

The advantage of the functional approach, is that the formal solution to the problem is written via quadrature. The moment of the field of any order can, in principle, be obtained by averaging the respective dynamic solution which is described by a multiple functional integral in the space of the virtual trajectories. The solutions to the first two moments of the field obtained from both the method of closed equations for the moments and the method of the path integral provide exactly the same result, under conditions of the Markov process approximation. The study of the field moments of order higher than two is easier with the use of the path integral, since this approach does not formally require a solution of the equations for the field moments. In fact, this advantage became apparent in the calculation of the moments of the intensity (signal strength), which are not described by a closed equation.

The current limitations of the path integral method could be rather associated with the still limited mathematical methods of analytical calculation of the trajectory integrals, despite the serious studies in that area.

Finally, we comment on studies concerning scattering on a rough sea surface. The basic conclusion is that, until now, there is no adequate theory that takes into account the combined effect of refraction, diffraction and scattering on random inhomegeneities of the refractive index in the volume of the tropospheric layer as well as scattering on a rough sea surface. The impact from a rough surface cannot be regarded as negligible in the general case, at least from a theoretical point of view. Known attempts are limited to the deterministic model of the sea surface or a semi-empirical approach in treating the scattering on the sea surface in the Kirchhoff approximation that is applicable to the coherent component of the scattered field. While the theory of wave scattering on a random surface is an established science by itself, the studies known to the author were concerned with scattering theory on a random surface while treating the propagation media in a very simplistic way, basically, as deterministic and uniform.

In applying the scattering theory to a ducting mechanism, the approach should eventually be modified in order to take into account the multiple effects of refraction, diffraction and scattering. As discussed in Atmospheric Boundary Layer and Basics of the Propagation Mechanisms“Atmospheric Boundary Layer and Key Propagation Mechanisms”, the sea waves could be responsible for the modulation of the turbulence spectrum in the near-surface layer, where the major action of radio propagation takes place. On the other hand, all previous advances in theory have been driven by some unresolved problems in experiments on the radio wave propagation phenomenon. So far, the currently available radio coverage prediction system is reported to be in reasonably good agreement with the results of existing models, at least for measured signal strength.

Parabolic Approximation to a Wave Equation in a Stratified Troposphere Filled with Turbulent Fluctuations of the Refractive Index

Let us consider the vertical electric dipole in a spherical coordinate system (r, ϑ, ϕ) with the spherical axis coming through the dipole. Assume that the dielectric permeability ε is a function of the radius only, e. g. ε = ε(r). The classical solution to the electromagnetic field can be derived by introduction of Debye’s potentials U and V.

The field components are given by the following equations:

Er=1r* U,
Eϑ=1εr2rϑεrU+jksin ϑVϕ,              Eq. 1
Eϕ=1εr sin ϑ2rϕεrUjkVϑ,
Hr=1r * V,
Hϑ=jkεsin ϑUϕ+1r2(rV)rϑ,                Eq. 2
Hϕ=jkεUϑ+1r sin ϑ2(rV)rϕ.

We assume the field to be a harmonic function of time – exp (-jωt), ω is a cyclic frequency of the electromagnetic radiation, k = ω/c is a wave number, and c is the speed of light. The notation Δ* represents a Laplace operator on a sphere:

*=1sin ϑϑsin ϑϑ+1sin2ϑ2ϕ2.             Eq. 3

The Maxwell equations:

rotH=jkεE,rotE=jkH                 Eq. 4

will be satisfied if potentials U and V obey the following equations:

1rr1εεrUr+1r2*U+k2εU=0,                 Eq. 5
1r2(rV)r2+1r2*V+k2εV=0.                  Eq. 6

The approximate Leontovitch’s boundary conditions for U and V have the form:

εrUr=jkεηrU,              Eq. 7
rVr=jkηrV              Eq. 8

with r = α, where a is the earth’s radius.

Let us assume that dielectric permittivity
εr
is stratified on average, i. e.:
εr=ε0r.
This means that
εr
is a random function of the space coordinates:
r=r, ϑ, ϕ: εr=ε0r+εr
where εr
is a random component present in each realisation of
εr
, and
εr=0.
Assume also that |δε/ε| ≪ 1. The angular brackets here and below mean averaging over the statistical ensemble of the fluctuations δε.
Let us assume that, the relationships between field components
E, H
and potentials U, V are still given by Eqs. (1) and (2). Strictly speaking, the Maxwell equations (4) cannot be satisfied for an arbitrary function
δεr
, and one can seek approximation to Eq. (4) where the divergent terms are small.
Introduce vector
x=αϑ, ϕr sin ϑ.
The divergent terms in Eq. (4) will be small when:
εxik                  Eq. 9
where xi is any of component of the vector
x.ε=ε0r+δεr.

Let us examine the vertically polarised field for which V = 0, U ≠ 0. For high frequencies when >> 1, we can select a “dedicated” direction of wave propagation x = αϑ aligned with the direction along the arc of the earth’s radius. With >> 1, the radial component of the electric field Er is related to the potential U via: Er ≈ -kα2 U.

Introduce a new function:

U1=εrU                   Eq. 10

which is governed by the equation:

r1εU1r+1εr22U1ϑ2U1ϑ2εεϑcos ϑcos ϑ+U11ε2εω21ε2εϑ2cos ϑε sin ϑεϑ+1εr2 sin ϑ2U1φ22εU1φεϕU1ε22εφ22εεφ2+k2U1=0             Eq. 11

and obeys the boundary condition:

U1r=jkU1η+1,    with r=α.             Eq. 12

Let us isolate the slow varying complex amplitude W1 in the wave field by:

U1=W1ejkαϑ                 Eq. 13

and introduce coordinate γ = (r sin ϑ)φ. The amplitude W1 is then given by equation:

2W1r2+2jkrW1ϑ+2W1γ2+k2εk2α2r2W1=εrW1r1r22W1ϑ2+jkrW12εεϑcos ϑsin ϑ+W1r21ε2εϑ21ε2εϑ2cos ϑε sin ϑεϑ+1r22εW1γεγW1ε2εγ22εεγ2            Eq. 14

at the right-hand side of which we have a correction term. Using qualitative arguments, it can be shown that in the case of a stratified medium ε = ε(r) the terms on the left have an order of magnitude of k2W1/m2 while on the right the terms containing derivatives are of the order of k2W1/m4. The right-hand side terms containing sin ϑ in the denominator will also be small if the following inequality holds:

xαm.              Eq. 15
The problem now is to evaluate the terms in the right-hand side of Eq. (14) for the case of random
δεr
caused by turbulence. Let us introduce coordinate z = r-α. The characteristic scales L of variations in W1 over x, y, z can be estimated as follows:
Lx~αm,   Lγ~λx,   Lz~mk.           Eq. 16

respectively. Having compared the terms in the left- and right-hand sides of Eq. (14) we obtain the inequalities:

εx2m2α2,   εγ21λx,   εz2k2m2             Eq. 17

which, being satisfied, allow one to treat the terms in the right-hand side of Eq. (14) as negligible. Therefore, with fulfilment of inequalities (17) the right-hand side of Eq. (14) can be replaced with zero.

The inequalities (17) put limitations chiefly on the intensity of small-scale fluctuations
δεr
. We can assume that the fluctuations in dielectric permittivity
δεr
are locally uniform and isotropic. Let us introduce a spatial spectrum Φε of fluctuations
δεr
by the equation:
Φεκ=116π3d3ρ·grad<δεr+ρδεr2>κ sin kρκ            Eq. 18
where
κ=κx, κγ, κz
is the wave vector of the fluctuations in the dielectric permittivity,
κ=κ.
We can evaluate the gradients in
δεr
by:
εxn2~ d3κ·κn2Φεκ              Eq. 19
where xn is any of the coordinates x, y, z. For the spatial spectrum of fluctuations
δεr
we can use the formula related to the equilibrium interval of the locally-isotropic turbulence:
Φεκ=0,033 Cε2κ11/3expκ2κm2                   Eq. 20

where cz is a structure constant, κm = 5,92/l0, l0 is an internal scale of turbulence. Substituting Eq. (20) into Eq. (17) we obtain:

0,033 Cε2 κm4/3kx,k2m2,m2α2.                  Eq. 21

Assuming reasonably short distances and high frequencies of radiation, e. g. m ≫ 1 and x ≤ 103km, we can see that the major impact comes from variations of the z-coordinate, and the inequality (22) thus reduces to:

0,033 Cε2πα22/3λl04/3 1.                Eq. 22
Here we took into account that k/m = k2/3 (α/2)1/3. Estimating the parameters
Cε2=1014cm2/3, α=8 500 km
, we observe that the coefficient with (λ/l0)4/3 has an order of 10–9.
Now we can conclude that with fulfilment of inequalities (15), (22) and z/α ≪ 1, the complex amplitude
W1r
obeys the approximate equation:
j2kW1x+2W1γ22W1z2+k2εmz1+δεr W1=0             Eq. 23

and boundary condition:

W1z=jkW1η+1 with z=0             Eq. 24

The potential U is given by:

Ur=ejkaϑεrW1r.                 Eq. 25

The next step is to figure out the correct representation of the field amplitude W1 at small distances from the source in order to obtain an initial condition. In particular, at the distance where the curvature of the earth can be neglected as well as the ray’s refraction we have to obtain a reflection formula in the parabolic approximation. Introducing the coordinates of the source:

  • x = 0;
  • γ = γ0;
  • z = z0;
  • and assuming that η ≫ 1 the reflection formula will be satisfied when:
W1x, ρ=2εx, ρxα exp jkxρρ02      Eq. 26
where
ρ=γ, z, ρ0=γ0, z0.
Taking into account the relationship:
lim|xxk2πxxexpjkxxρρ2=ρρ
where
δρ
is a Dirac’s delta-function defined by:
δρ=14π2d2κ exp jκρ,             Eq. 27

we make a transition in Eq. (27) with x → 0 and obtain an expression for singularity at the source:

lim|x0W1x, ρ=4πjε0, ρ0αkδρρ0.          Eq. 28
Equation (27) shall take place in a line-of-sight region, i. e. with
x2αz0.
Taking into account that the inequality (16) is to hold as well, we obtain the limitations at the height of the source z0 under which the expression (29) is valid:
z0α/m2.                 Eq. 29

It should be noted that in the case of very small heights of the source, another approach can be used and the solution to Eq. (24) should be tailored with the solution of Veil-Van-der-Paul at small distances x.

In conclusion, one should note that there is some ambiguity with the definition of the complex amplitude via Debye’s potential U. Indeed, equations similar to Eq. (24) can be obtained for the other functions W2 and W3 defined by the following expressions:

W2=εrsin ϑ·ejkαϑU,W3=εr sin ϑjkαϑU.                Eq. 30

The area of applicability of the parabolic approximation for W3 is defined by the same set of inequalities (16), (22) while for W3 the inequality (16) can be replaced by a less stringent one:

xm/k.              Eq. 31

Let us summarize the major steps undertaken in this section:

  • The components of the electromagnetic field
    E
    and
    H
    can be represented via Debye’s potentials U and V. Even in the presence of a random turbulent component of the refractive index the de-polarization effects can be neglected and, therefore, the field components are governed by independent equations for U and V.
  • The wave equations for both potentials U and V can be approximated by parabolic equations for a slow varying complex amplitude. In this approximation all radiated waves are travelling only in one direction, away from the source. The Leontovich boundary conditions are used along with the parabolic equation. The conditions in the source are also modified in a parabolic approximation.
  • During the transformation from the Helmoltz equation to a parabolic equation we actually replaced the earth sphere with a cylinder of the same radius and then flattened the cylinder by means of compensating the earth’s curvature by “ray’s curvature” in the opposite direction. This is a well-known “flat earth” approximation with modified refractivity and locally Cartesian coordinates.;

Green Function for a Parabolic Equation in a Stratified Medium

Let
δεr
= 0 in Eq. (24) and take into account that for the frequencies above 1 GHz the boundary condition (25) is reduced to
Wr
at z = 0. By definition, a Green function
Gr, r
is governed by the equation:
2jkGx+G+k2εmz1G=4πδrr*.            Eq. 32

where:

=2γ2+2x2,

obeys the boundary condition:

Gr, r1=0 with z=0,
and is a continuous and limited function for all
rr.
. Introducing a Fourier transform over the transverse coordinates x, y:
Gr, r=14π2dp  dqGp, q, z, z exp jpxx+jqγγ

and taking into account the definition of the δ-function, we obtain:

d2Gdz2+k2εmz12kpq2G=4πδzz*.               Eq. 33

We will seek a solution to Eq. (34) as a composition over functions ΨE(z) which are governed by the equation:

d2Ψdz2+kεmz1Ψ=Eψ.                 Eq. 34

and satisfy the boundary conditions:

ΨEz=0,    z=0;   ΨEz=0,   z.            Eq. 35

Equation (35) is similar to the one-dimensional Schrödinger equation for a particle in a field with potential energy which is proportional to the term (εm(z)-1). As known, in the case of unlimited potential (in Eq. (35) m(z)-1)→2z/α with z → ∞)), the motion of the particle in a stationary state is infinite to z → +∞ and the spectrum of the eigen values E is purely continuous.

Therefore, the function
Gp, q, z, z
can be thought of as a composition of the eigen functions of the continuous spectrum which obey Eq. (35) and the boundary conditions (36):
Gp, q, z, z= dEBp, q, E, z ΨEz.               Eq. 36

The eigen functions of a continuous spectrum satisfy the orthogonality conditions:

0ΨEzΨE1*zdE=δEE1                    Eq. 37

and “completeness”:

ΨEzΨE*z=δzz.                    Eq. 38

The conditions (38) and (39) are known from quantum mechanics where they are proven for potentials limited at infinite z. Using methods similar to those, one can prove Eqs. (38) and (39) for a potential with a linear increment at infinity. The sign * means a complex conjugate, here and below.

Substituting Eq. (37) into Eq. (34) and taking into account Eqs. (38) and (39) we obtain:

Bp, q, E, z=4πψE*zE2kpq2,              Eq. 39

and then performing the integration in Eq. (37) we obtain:

Gr, r=2πkxx1/2expjπ4+jk2γγ2xxdEΨEzΨE*zejE2kxx.                 Eq. 40

In all practical problems related to radio wave propagation and scattering in the troposphere, the non-uniformities of the refractive index are localised in the boundary layer of the atmosphere and εm(z)-1 can be approximated by a linear function 2z/α with z → ∞.

Thus Eq. (34) has two solutions:
χE+z
and
χEz
behaving asymptotically as Airy functions
w2Eμ2μz
and
w2E/μ2μz
, where μ = k/m. We define the function w1(t) via a contour integral:
w1t=1πΓdξ exp ξ33+tξ                Eq. 41

where contour Γ comes from infinity along the ray arg ξ = -2π/3 to zero and then from zero to infinity ξ → +∞ along the real axis ξ. The function w2(t) is defined by a complex conjugate of Eq. (42).

The functions
χE+z
and
χEz
describe eigen waves coming to and from infinity. The regular (with z = 0) solution for ΨE(z) given by the composition.
ΨEz=AEχEzSEχE+z                 Eq. 42
where
AE1/2πμ
is a coefficient defined from normalisation to the delta-function Eqs. (38) and (39). The term S(E) (normally called the S-matrix in quantum mechanics) is determined by a kind of non-uniformity εm(z)-1 and the boundary condition at z = 0. In the case of an ideal boundary condition at z = 0, Eq. (36), the law of conservation for a number of particles holds true and |S(E)| = 1. Therefore, for real E, S(E) = exp (-2jδ(E)), where δ(E) is the real phase of scattering.

The S-matrix S(E) has poles En in the upper half-space of E and residues in En completely define the field in the shadow region. The value of ΨE(z) in the pole E = En provides a normalised height function of normal wave with number n:

ΨEnz=χEnzχnz.                Eq. 43

The normal wave χn(z) with complex value of the propagation constant En(Im En > 0) is unlimited at infinite z. This unlimited growth is a consequence of the exponential decay of the field over the x-coordinate, since the field at z → ∞ is created by radiation propagating from points distant to -∞ over the x-coordinate.

Let us calculate integral (41) for normal refraction when εm(z) = 1+2z/α at z > 0 and the eigenfunctions ΨE(z) the continuous spectrum (43) are defined by:

ΨEz=12πμw2Eμ2μzw2Eμ2w1Eμ2w1Eμ2μz.                Eq. 44

One can expect that, in a line-of-sight region, expression (41) shall represent a “reflection formula” in accordance with geometric optics, i. e. represent the field as a sum of the direct wave and the wave reflected from the earth’s surface.

Consider ray equations corresponding to the waves radiated either in an upward direction from the source (z > z0) or downward from the source (z > z0), and define the relation between the grazing angle of the ray and stationary value of E. We define the hit-angle ϑ to be a sliding angle of the ray relative to the surface z > z0 at the point of the source location
r0=0, γ0, z0.

The wave propagating in the direction up from the source can be expressed as follows:

V+~ejxE/2kw1Eμ2μzw2Eμ2μz0.               Eq. 45

Using an asymptotic expression for the Airy function we obtain the trajectory equation by defining the extremum of the phase in Eq. (46):

z+x=x2μ34k2+xkμ3z0E+z0.                Eq. 46

A similar equation can be obtained for the ray in the direction down from the source:

zx=x2μ34k2xkμ3z0E+z0.                Eq. 47

Differentiating Eqs. (47) and (48) we can find the equation for angle ϑ and its relation to the stationary value of E:

ϑ=±α tan 1kμ3z0E,                  Eq. 48
E=μ3z0k2 tan2 ϑ.                   Eq. 49

As observed from Eq. (49) the values of E satisfying the inequality E ≤ μ3 z0 determine the real hit-angles ϑ, those, in turn, represent uniform plane waves with angle ϑ/m between the normal to their wavefront and tangential to the radius r = α+z0 at the point of source location. Trajectory (48) has a minimum in:

xmin=2kμ3μ3z0E and zminxmin=Eμ3.
From that it follows that the values E < 0 determine the waves reflected from the earth’s surface. The sector
α tankz0/m3<ϑ<a tankz0/m3
determines the “space” rays (direct wave) related to a stationary value of E from the interval:
0<E<μ3z0.

The ray trajectories described by Eqs. (47) and (48) are shown in Figure.

Trajectory of the rays
Ray trajectories in the case of normal refraction
Therefore, in the line-of-sight region, i. e.
x<2αz+z0
we can select two areas of integration in Eq. (2.43):
I1=0dEΨEzΨE*z0ejE2kx,              Eq. 50
I2=0dEΨEzΨE*z0ejE2kx,              Eq. 51

Let us substitute Eqs. (45)–(47) into Eq. (51) and assume z > z0 and:

kz0/m1.                    Eq. 52

The integration in Eq. (51) can be performed using the method of stationary phase. The result is given by:

I1=θ2αzz0xI11+I12,                  Eq. 53
I11=2μA22kπxexpjπ4+jk2xzz02+jkx2αz+z0j112mxα3,
I12=kαξ4π1zz04·expjkm3/2z3/2+z03/2+jkx2αz+z0j112mξα3                 Eq. 54
where
ξ=x2αz+z0.
Consider integral (52) and make use of the asymptotic expression for Airy functions w1,2 (E/μ2) for large positive arguments:
w1Eμ2w2Eμ2Eμ21/2exp23μ3E3/2.             Eq. 55

Substituting Eq. (55) into Eq. (45) we obtain S(E) = 1. Let us introduce the Airy function v(t) by the relationship:

2jνt=w1tw2t                  Eq. 56

and the variable τ = E/μ2-μz0. The integral (52) can then be transformed into:

I214πejx2kμ3z0μz0ejxτ2kμ2ντντμzz0dτ.              Eq. 57

Without significant error the low limit in the integral (58) can be extended to -∞. Substituting into Eq. (58) the integral representation for ν(t):

νt=12πexpjξ33+jξtdξ              Eq. 58

and performing the integration, we obtain:

I2k2πxexpjπ4+jk2xzz02jkx2αz+z0jmxα3.                Eq. 59

Equation (60) describes the direct wave in the range of distances:

2αzz0<x<2αz+z0
and transforms into Eq. (53) with
x=2αzz0
. Therefore in a whole line-of-sight region the Green function (41) can be represented by superposition of the direct wave and the wave reflected from the earth’s surface.
Let us consider the Green function’s representation (41) in a shadow region, i. e.
x>2αz+z0.
In this case the contour of integration over dE can be closed in the upper half-plane and the integral is represented by a residue sum over the poles of S(E). In the particular case of normal refraction and ideal boundary conditions, the poles are given by:
w1Eμ2=0 and En=k2m2τnejπ/3, τ1=2,338,    τ2=4,02, ...                  Eq. 60

As a result, we obtain the Green function representation in terms of the normal wave series:

Gr, r=2πkxxexpjπ4+jk2γγ2xx·n=1expjEn2kxxw1Eμ2μzw1Eμ2w1Eμ2μzw1Eμ2.                 Eq. 61

In conclusion, the Green function has been built as an expansion over the set of complete and orthogonal eigenfunctions of the continuous spectrum. In the case of normal refraction, such a representation produces final formulas for the field of the point source that are similar to those introduced by Fock on the basis of the classical approach of the contour integral.

It is worthwhile to emphasise the purpose of the above presentation for a Green function. In a problem of multiple scattering the whole path is involved in the scattering process and one needs to have a representation for the field (or Green function) which is equally applicable in the line-of-sight and shadow regions. Also, in the process of obtaining the closed set of equations for statistical moments of the field, the Green function has to possess some properties of orthogonality and completeness which can be achieved by using the above expansion over the set of eigenfunctions of the continuous spectrum.

Author
Author photo - Olga Nesvetailova
Freelancer
Literature
  1. Fock, V.A. Electromagnetic Diffraction and Propagation Problems, Pergamon Press, Oxford, 1965.
  2. Fock, V.A. Distribution of currents excited by a plane wave on the surface of conductor, J. Exp. Theor. Phys., (Nauka), 1945, 15 (12), 693–702.
  3. Leontovitch, M.A. Method of solution to the problems of em wave propagation over the boundary surface, Izv. Acad. Nauk, Ser. Fiz., 1944, 8 (1), 16–22.
  4. Leontovitch, M.A., Fock, V.A. Solution to a problem of em wave propagation over the Earth’s surface using the method of parabolic equation, J. Exp.Theor. Phys., (Nauka), 1946, 16 (7), 557–573.
  5. Fock, V.A. Theory of wave propagation in a non-uniform atmosphere for elevated source, Izv. Acad. Nauk, Ser. Fiz., 1960, 14 (1), 70–94.
  6. Krasnushkin, A.E. The expansion over normal waves in a spherically-stratified medium, Dokl. Acad. Nauk SSSR, 1969, 185 (6), 1262–1265.
  7. Kravtsov, Yu.A., Orlov, Yu.I., Geometric Optics of Non-uniform Media, Nauka, Moscow, 1980, 340 pp.
  8. Rotheram, S. Radiowave propagation in the evaporation duct, The Marconi Rev., 1974, 67 (12), 18–40.
  9. Booker, H.G., Walkinshaw W. The mode theory of tropospheric refraction and its relation to waveguides and diffraction, in Meteorological Factors in Radiowave Propagation, The Royal Society, London, 1946, pp. 80–127.
  10. Andrianov, V.A. Diffraction of UHF in a bilinear model of the troposphere over the Earth’s surface, Radiophys. Quantum Electron., 1977, 22 (2), 212–222.
Footnotes
Sea-Man

Did you find mistake? Highlight and press CTRL+Enter

Февраль, 28, 2025 61 0
Add a comment


Notes
Text copied