Computational Electromagnetics
Chen Aixin School of Electronic and Information Engineering
Chapter 2 Analytical Methods
2.1 Introduction
z
The mos
t commonly used analytical methods in solving EM-related problems include:
– – – –
(1) separation of variables (2) series expansion (3) conformal mapping (4) integral methods
z
Perhaps the most powerful analytical method is the separation of variables.
3
2.2 Separation of Variables
z
The method of separation of variables (sometimes called the method of Fourier)
It entails seeking a solution which breaks up into a product of functions, each of which involves only one of the variables. For example,
A solution of the form in Eq. is said to be separable in x, y, z, and t .
4
z
To determine whether the method of independent separation of variables can be applied to a given physical problem, we must consider
– – –
the PDE describing the problem the shape of the solution region the boundary conditions
z
The three elements that uniquely define a problem.
5
z
We shall always take these three major steps:
– –
–
(1) separate the (independent) variables (2) find particular solutions of the separated equations, which satisfy some of the boundary conditions (3) combine these solutions to satisfy the remaining boundary conditions
6
z
The homogeneous scalar wave equation
7
If we let an arbitrary constant ?k2 be the common value of the two sides
solution or
8
z
Notice that if k = 0, the time dependence disappears and Eq. becomes Laplace’s equation.
9
2.3 Separation of Variables in Rectangular Coordinates
z
Helmholtz equation In rectangular coordinates, becomes
z
z
We let
10
z
We obtain Each term must be equal to a constant since each term depends only on the corresponding variable. We conclude that so that
11
z
z
z
Suppose we choose Then Or Introducing the time dependence gives
z
z z
12
z z
where ω = kc is the angular frequency of the wave. The solution in Eq. represents a plane wave of amplitude A propagating in the direction of the wave vector k = kxax + kyay + kzaz with velocity c.
13
2.4 Separation of Variables in Cylindrical Coordinates
z z z
Helmholtz’s equation Let Substituting into Eq. and dividing by RFZ/ρ2 yields where n = 0, 1, 2, . . . and n2 is the separation constant.
14
z
z z
Thus And Dividing both sides of Eq. by ρ2 leads to
z
z
where μ2 is another separation constant.
15
z
Hence And
z
z
If we let
16
z
The three separated equations become
z
The solutions or or
17
z
Let x = λρ and replace R by y; the equation becomes Bessel’s equation It has a general solution of the form
where Jn(x) and Yn(x) are Bessel functions of the first and second kinds of order n and real argument x. Yn is also called the Neumann function.
z
18
z
If x in Eq. is imaginary so that we may replace x by jx, the equation becomes modified Bessel’s equation This equation has a solution of the form
where In(x) and Kn(x) are modified Bessel functions of the first and second kind of order n.
z
19
z
For small values of x, some typical Bessel functions (or cylindrical functions)
20
21
z
For the Bessel function of the first kind
where Γ(k + 1) = k! is the Gamma function.
z
This is the most useful of all Bessel functions.
22
z
Some of its important properties and identities (also hold for Yn(x))
23
z
For the Neumann function
–
when n > 0
–
If n = 0
where γ = 1.781 is Euler’s constant and
24
z
For the modified Bessel function of the first kind
z
For the modified Bessel function of the second kind
25
–
If n > 0
–
if n =0
26
z
Other functions closely related to Bessel functions are Hankel functions of the first and second kinds, defined respectively by
z
Hankel functions are analogous to functions exp(±jx) just as Jn and Yn are analogous to cosine and sine functions.
27
z
This is evident from asymptotic expressions
28
z
z
With the time factor ejωt, and represent inward and outward traveling waves, respectively, while Jn(x) or Yn(x) represents a standing wave. With the time factor e-jωt , the roles of and are reversed.
29
z
Any of the Bessel functions or related functions can be a solution to Eq. depending on the problem. If we choose R(ρ) = Jn(x) = Jn(λρ) with the solutions of and and apply the superposition theorem
30
z
z
The solution is
z
Introducing the time dependence, we get
where ω = kc.
31
2.5 Separation of Variables in Spherical Coordinates
z
The wave equation becomes
z
z
We substitute into the equation. Multiplying the result by r2 sin2θ/RHF gives
32
z
Since the left-hand side of this equation is independent of φ, we let
where m, the first separation constant, is chosen to be nonnegative integer such that U is periodic in φ.
z
This requirement is necessary for physical reasons that will be evident later.
33
z
Thus Eq. reduces to
where λ is the second separation constant.
z
Let λ = n(n+1), the separated equations are
34
z
The solution to Eq. Is or
z
35
z
If we let
,Eq. Becomes
which has the solution
36
z
Functions zn(x) are spherical Bessel functions and are related to ordinary Bessel functions Zn+1/2 according to
Zn+1/2(x) may be any of the ordinary Bessel functions of half-integer order, Jn+1/2(x), Yn+1/2(x), In+1/2(x), Kn+1/2(x), ,and while zn(x) may be any of the corresponding spherical Bessel functions jn(x), yn(x),in(x),kn(x), ,and .
37
z
Bessel functions of fractional order are, in general, given by
where J-ν and I-ν are, respectively, obtained from Eqs. by replacing ν with-ν.
38
z
Since Gamma function of half-integer order is needed in Eq., it is necessary to add that
39
z
Thus the lower order spherical Bessel functions are as follows:
40
z
Other zn(x) can be obtained from the series expansion in Eqs. of Bessel functions of fractional order or the recurrence relations and properties of zn(x) presented as below
41
z
By replacing H in Eq. with y, cosθ by x, and making other substitutions, we obtain Legendre’s associated differential equation
42
z
Its general solution is of the form
where and are called associated Legendre functions of the first and second kind, respectively.
43
z
and can be obtained from ordinary Legendre functions Pn(x) and Qn(x) using
where -1 < x < 1. We note that
44
z
Pn(x) and Qn(x) are Legendre functions of the first and second kind, respectively.
45
z
Typical associated Legendre functions are:
z z
Higher-order associated Legendre functions can be obtained using Eqs. of and is unbounded at x = ±1, and hence it is only used when x = ±1 is excluded.
46
z
Substituting and applying superposition theorem,we obtain
z
Note that the products H(θ)F(φ) are known as spherical harmonics.
47
2.6 Concluding Remarks
z z z
Analytic solutions are of major interest as test models for comparison with numerical techniques. The emphasis has been on the method of separation of variables, the most powerful analytic method. In the course of applying the method of separation of variables, we have encountered some mathematical functions such as Bessel functions and Legendre polynomials.
48
z
We conclude this chapter by remarking that the most satisfactory solution of a field problem is an exact analytical one. In many practical situations, no solution can be obtained by the analytical methods now available, and one must therefore resort to numerical approximation, graphical or experimental solutions. (Experimental solutions are usually very expensive, while graphical solutions are not so accurate).
49