Electromagnetic Field Theory A Collection of Problems
Electromagnetic Field Theory A Collection of Problems
Springer Vieweg
© Vieweg+Teubner Verlag | Springer Fachmedien Wiesbaden 2013
This work is subject to copyright. All rights are reserved by the Publisher, whether the whole or part of
the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation,
broadcasting, reproduction on microfilms or in any other physical way, and transmission or information
storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology
now known or hereafter developed. Exempted from this legal reservation are brief excerpts in connection
with reviews or scholarly analysis or material supplied specifically for the purpose of being entered and
executed on a computer system, for exclusive use by the purchaser of the work. Duplication of this pub-
lication or parts thereof is permitted only under the provisions of the Copyright Law of the Publisher’s
location, in its current version, and permission for use must always be obtained from Springer. Permissions
for use may be obtained through RightsLink at the Copyright Clearance Center. Violations are liable to
prosecution under the respective Copyright Law.
The use of general descriptive names, registered names, trademarks, service marks, etc. in this publication
does not imply , even in the absence of a specific statement, that such names are exempt from the relevant
protective laws and regulations and therefore free for general use.
While the advice and information in this book are believed to be true and accurate at the date of publica-
tion, neither the authors nor the editors nor the publisher can accept any legal responsibility for any errors
or omissions that may be made. The publisher makes no warranty, express or implied, with respect to the
material contained herein.
1 Fundamental Equations 1
2 Electrostatic Fields 12
Appendix 260
References 265
Index 266
Symbols
A Vector potential
a Surface, Distance
B Magnetic flux density
C Capacitance, Integration path
cij Capacitance coefficients
D Electric flux density
E Electric field
F Vector potential, Force
H Magnetic field
I, i Electric current
J Current density
K Current sheet
Lik Inductance
M Dipole moment
m
Dipole moment density
P Power
pv Power loss density
pij Potential coefficients
Q Charge
r
R, Position vectors
S Poynting vector
t, T Time
u Voltage
U, V Potential
v Volume, Velocity
W Energy
w Energy density, Complex variable
Z Wave impedance
μ Permeability
κ Conductivity
Volume charge density, Coordinate
σ Surface charge density, Step function
λq Line charge density
λ Wavelength
Ψe,m Electric/Magnetic flux
ϕ, ψ Potential function
Ω Solid angle
ω Angular frequency
Φ Radiation pattern, Potential
A∗ conjugate-complex quantity.
B = μH
; = εE
D
a = 0;
Bd
Dda = dv ; a+ ∂
Jd dv = 0 ; J = κE.
∂t
a a v a v
All fields depend on the space variable r and the time t and have to be integrated over
surfaces a with a contour C or over surfaces a enclosing the volume v. The current
r , t) and the charge density (r, t) satisfy the continuity equation, which
density J(
implies the conservation of charge.
= J + ∂ D ;
rot H = − ∂B ;
rot E = μH
B ; = εE
D
∂t ∂t
= 0;
div B = ;
div D div J + ∂ /∂t = 0 ; J = κE
.
The force F of the electromagnetic field on a charge with velocity v , namely the Lorentz
force, is given by
dF = (E + v × B)
dv.
For a surface a with normal n directed from domain (1) with material properties
ε1 , μ1 , κ1 into domain (2) with material properties ε2 , μ2 , κ2 the following boundary
conditions can be deduced
n × E 1 = 0 ; n × H
2 − E 2 − H 1 = 0, K
a a
n B 1 = 0 ;
2 − B n D 1 = σ ;
2 −D
n J2 − J1 = −
∂σ
.
a a a ∂t
Here σ is a space- and time-dependent surface charge density within the boundary
layer a.
The following equations hold in regions with sectional homogeneous material properties,
that are characterized by scalar parameters ε, μ, and κ.
Electrostatic Fields
The equations for the electric field of static charges are
= 0;
rot E = ;
div D = εE
D
= −grad V ;
E ΔV = − /ε .
The electric field at coordinate rp is given by the gradient of the scalar electrostatic
potential V (rp ). Charge distributions are classified as space charge densities (rq )
in a volume v, surface charge densities σ(rq ) on a surface a, or line charge densities
λq (rq ) on a line with contour C. In a homogeneous space of permittivity ε the charge
distributions possess the following potentials and fields at point rp .
1 (rq ) 1 r
V (rp ) = dv ; E(rp ) = −grad V (rp ) = (rq ) 3 dv
4πε r 4πε r
v v
1 σ(rq ) rp ) = −grad V (rp ) = 1 r
V (rp ) = da ; E( σ(rq ) da
4πε r 4πε r3
a v
1 λq (rq ) rp ) = −grad V (rp ) = 1 r
V (rp ) = ds ; E( λq (rq ) ds
4πε r 4πε r3
C C
with r = rp − rq , |r| = r. Similarly the field of n point charges Qi at positions rqi ;
i = 1, 2, 3, . . . , n is given by
1
n n
Qi rp ) = −grad V (rp ) = 1 rp − rqi
V (rp ) = ; E( Qi .
4πε i=1 |rp − rqi | 4πε i=1 |rp − rqi |3
1 Fundamental Equations 3
Line- and point charges are quantities that were introduced to simplify the field compu-
tation, but due to the unbounded energy of the fields they are just theoretical models.
over a surface a leads to the electric Flux
The integral of the electric flux density D
Ψe = da .
D
a
If the electric flux that runs through a tube is constant Ψe = const , then the electric
field is parallel to the side surface of the tube. Thus the determination of these flux
tubes also defines electric lines of force. In this way often lines of force are obtained
easier than by solving the corresponding differential equation
= 0.
ds × E
The transition trough a dipole layer in the direction of the surface normal n(rq ) from
region (1) to region (2) induces a change of m/ε in the potential and for the electric
field applies
If the dipoles have a homogeneous density m 0 = n(rq ) m0 , then the potential is de-
termined by means of the solid angle Ω, which is given by the border contour of the
surface.
m0 3 m0
V (rp ) = n(rq )r/r da = − Ω ; Ω = − n(rq ) r/r3 da
4πε 4πε
a a
4 1 Fundamental Equations
In a system of n conducting and charged bodies a simple linear relation between the
potentials Vi and the charges Qi , i = 1, 2, . . . , n exists.
n
n
Vk = pki Qi ; Qk = cki Vi
i=1 i=1
n
Qk = Ck∞ Vk + Cki (Vk − Vi )
i=1
Here pki and cki are the potential and capacitance coefficients and Cki are mutual
capacitances.
1 Fundamental Equations 5
= 0;
rot E div J = 0 ; J = κE
; = −grad V ;
E ΔV = 0 .
As there is an analogy to the electrostatic field in source free regions, most relations of
the previous section remain valid.
Magnetic Field of Stationary Currents
The field equations are
= J ;
rot H = 0;
div B = μH
B
= rot A
B ; = 0;
div A = −μJ .
ΔA
is
Hence the solution for the magnetic field H
rp ) = B = 1
H( rq ) × r dv ;
J( r = rp − rq .
μ 4π r3
v
The magnetic Flux Ψm through a surface a with the contour C is defined by the
integrals
Ψm =
Bda =
rot A da = ds .
A
a a C
As the curl of the magnetic field is zero in space with no currents one can also use a
scalar potential Vm for the description of the magnetic field. For a line current the
solution is
1 r
H = −grad Vm ; ΔVm = 0 ; Vm (rp ) = m
m 3 da ; m m = μnI.
4πμ r
a
Plane magnetic fields can be calculated by a complex magnetic potential P m (z). With
ΔPmr = 0 and ΔPmi = 0 it follows
The second integral over the volume of the conductor vL assumes a regular behavior
of the magnetic field and the vector potential at large distances.
In a system of n current-carrying conductors with current densities Ji the energy of
the field is given by
n n
1
n n n n
Wm = H k dv = 1
iB ik dv = 1
Ji A Lik Ii Ik
2 i=1 2 i=1 2 i=1
k=1 v k=1vLi k=1
For two line currents Ii and Ik with the contours Ci and Ck it follows
1 ik dsi = Ψmik ;
Lik = A i = k .
Ik Ik
Ci
= JE + J + ∂ D ;
rot H = − ∂B ;
rot E = 0;
div B = E
div D
∂t ∂t
∂ E
div JE + = 0; B = μH
; D = εE ; J = κE.
∂t
r , t) and ϕ(r, t) satisfy the differential equations
The electrodynamic potentials A(
− κμ
∂A ∂ 2A
∂ϕ ∂2ϕ E
ΔA − με 2 = −μJE ; Δϕ − κμ − με 2 = −
∂t ∂t ∂t ∂t ε
+ κμϕ + με ∂ϕ = 0 ; B
div A = rot A; E = − ∂ A − grad ϕ .
∂t ∂t
In case of a sinusoidal time dependency of the field the corresponding equations for the
complex amplitudes are
E
+ ω 2 μεA
− jωκμA
ΔA = −μJ ; Δϕ − jωκμϕ + ω 2 μεϕ = −
E
ε
+ κμϕ + jωμεϕ = 0 ;
div A = rot A
B ; = −jω A
E − grad ϕ .
The Poynting theorem is the law of energy conservation for a volume v surrounded by
the surface a
Here P s is the complex power of all sources in v, W m and W e are the time-averaged
energies stored in the electric and magnetic field
⎧ ⎫
1 ⎨ ∗ ⎬ μ 2 dv
W m = Re H B dv = |H|
4 ⎩ ⎭ 4
v v
8 1 Fundamental Equations
⎧ ⎫
1 ⎨ ⎬ ε
We = Re E ∗ dv
D = 2 dv
|E|
4 ⎩ ⎭ 4
v v
holds.
Quasi Stationary Fields
The contribution of the displacement current ∂ D/∂t to the magnetic field is negligible in
case of slowly time-varying fields compared to that of the conduction current. Therefore
the field equations are
rot H = − ∂B
= J ; rot E
∂t
div B = 0 ; div J = 0 ; B = μH ; J = κE
∂A ∂ϕ
= rot A
B ; = κμ
ΔA ; Δϕ = κμ ; + κμϕ = 0
div A
∂t ∂t
= − ∂ A − grad ϕ .
E
∂t
Electromagnetic Waves
The fields of plane waves in a homogeneous space of permittivity ε and permeability μ
depend only on one linear coordinate, here z, and on the time t (wave impedance Z,
phase velocity v, wavelength λ).
1,2 (z, t) = Re H
H exp j(ωt ∓ βz)
1,2
1,2 and H
where E 1,2 are complex amplitudes with
1
H 1,2 = ±
;
ez × E 1,2
E
1,2 = ±Z H 1,2 ×
ez ; β 2 = ω 2 με = (2π/λ)2 .
Z
A F T E = eF T E
TM =
eAT M ;
ΔF T E + β 2 F T E = 0 ; ΔAT M + β 2 AT M = 0
=H
1 1
H TE + HTM = rot A TM + rot rot F T E
μ jωμε
1 1
E
=E
T E + E T M = − rot F T E +
rot rot A TM ; β 2 = ω 2 με ; κ = 0.
ε jωμε
A Hertzian dipole located at the point of origin with the moment M = ez M = ez Is
and s λ possesses in the far-field region r λ the field
r ) = jβ M × r exp(−jβr) ; √ 2π
H( β = ω με =
4π r r λ
2
r) = β
E(
r
× M × r exp(−jβr) = Z H × r .
jω4πε r r r r
The radiation pattern Φ is given by
2
Z |
|M r
Φ= sin2 ϑ ; cos ϑ = ez .
8 λ r
10 1 Fundamental Equations
A linear antenna with the given current distribution I(rq ) on a line with contour C
leads to the following field in the far-field region:
rp ) = jβ exp(−jβrp ) I(rq ) exp jβrq rp
H( dsq ×
rp
4π rp rp rp
C
×
rp ) = Z H rp μ
E( ; Z = ; r = rp − rq ; rp = |rp | .
rp ε
= 1 rot A
H = 1 grad Φ(x, y) × ez exp(∓jβz) = ± 1 ez × E
μ μ Z
√
ΔΦ(x, y) = 0 ; β = ω με .
The integrals of the fields bring out the voltage u(z, t) and the current i(z, t) of the
transmission line
∂ 2u ∂2u ∂2i ∂2i ∂u ∂i ∂i ∂u
2
= L C 2 ; 2
= L C 2 ; − = L ; − = C
∂z ∂t ∂z ∂t ∂z ∂t ∂z ∂t
where L is the inductance per unit length and C is the capacitance per unit length
in z-direction. The solutions are
1
u(z, t) = u1,2 (z ∓ vt) ; i(z, t) = i1,2 (z ∓ vt) ; v=√
L C
L
Z = ; u1 = Z i1 ; u2 = −Z i2 .
C
TE- and TM-waves within ideal cylindrical waveguides in z-direction are described by
the potentials A
T M and AT E .
A ez AT M = ez U T M (u1 , u2 ) exp ∓jβzT M z
TM =
Here (u1 , u2 ) are the orthogonal curvilinear coordinates in the plane z = const .
1 Fundamental Equations 11
At the perfect conducting boundaries of the waveguide with the contour C the boundary
conditions
∂U T E
U T M |C = 0 ; =0
∂n C
have to be satisfied. Thereby the phase constants βzT E and βzT M can be identified.
If the source distributions JE (rq , t) and E (rq , t) are given in a homogeneous space,
then the fields can be deduced from retarded potentials.
rp , t) = μ JE (rq , t∗ )
A( dv
4π |rp − rq |
v
1 E (rq , t∗ )
ϕ(rp , t) = dv
4πε |rp − rq |
v
∗ √
Therein t = t − r/vph is the retarded time with the phase velocity vph = 1/ με.
In case of a sinusoidal time dependency the complex amplitudes are given by
rp ) = μ exp(−jβ|rp − rq |)
A( JE (rq ) dv
4π |rp − rq |
v
1 exp(−jβ|rp − rq |)
ϕ(rp ) = E (rq ) dv
4πε |rp − rq |
v
√
with β = ω με.
2. Electrostatic Fields
Consider two conducting hollow spheres of radius r = b and r = c > b. The inner
sphere contains a spherical space charge with density [As/m3 ] and radius a < b − e,
which is eccentrically positioned at a distance e from the center. The homogeneous
permittivity in b < r < c is ε, while the permittivity in the remaining space is ε0 . The
potential on the outer sphere at r = c is V0 .
Find the potential on the sphere r = b, when it is uncharged.
Now the charge on the sphere r = b is set to Q1 . What is the Potential on the surface
and what is the charge Q0 of the sphere at r = c?
? >
A
=
In b < r < c the electric field is of spherical symmetry and is described by the potential
V (r).
V (r) = C1 + C2 /r ; E = −grad V = −er ∂V = er E(r)
∂r
E(r) = C2 /r2 ; D(r) = εE(r) = ε C2 /r2
If the sphere r = b is uncharged, the electrical flux through a surface b < r = const < c
is equal to the charge quantity in the region r < b.
Q
D(r) 4πr2 = 4/3 πa3 = Q = 4πε C2 ⇒ C2 =
4πε
Hence the potential is
Q 1 c
V (r) = C1 + ; b≤r≤c; V (r) = V0 ; r ≥ c.
4πε r r
Q0 = −(Q + Q1 ) + 4π ε0 cV0 .
Finally, it should be noted that the position of the charge density has no influence
on the potential in > b.
j y
a r
a
cd j dc
N
2 2 bb C 0 C C 0
e
where the real part vanishes on the cylinder = a. Rearranging the arguments of the
logarithm functions results in
λq1 z − a2 /b b a−b λq2 z + a2 /b
P e (z) = ln + ln + ln + ln(b/a) .
2πε z−b a b−a 2πε z+b
C 0
C n
0 .3 d = a = { 1 0 -3
; 2 ·1 0 -3
; 5 ·1 0 -3
; 1 0 -2
}
0 .2
d = a
0 .1
0 0 .2 0 .4 0 .6 0 .8 b = a
C
C n
-3 -3 -3 -2
d = a = { 1 0 ; 2 ·1 0 ; 5 ·1 0 ; 1 0 }
0 .0 8
0 .0 4
d = a
0
0 0 .2 0 .4 0 .6 0 .8 b = a
Fig. 2.2–2: Equipotential lines and electric lines of force for λq1 /λq2 = −2 (top) and
λq1 /λq2 = 2 (bottom)
2.3 Singular Points and Lines in the Field of Point Charges 17
With the capacitance coefficients cik the equations for the charges read as
V1 p0 − V2 p
λq1 = c11 V1 + c12 V2 =
p20 − p2
V2 p0 − V1 p
λq2 = c21 V1 + c22 V2 = ,
p20 − p2
Finally rewriting the equations of the the charges in terms of potential differences
On a straight line in the space of permittivity ε the point charges Q1 and Q2 are placed
in front of the point charge Q with the distances r1 and r2 .
What are the values of the charges Q1 and Q2 , in order that a singular point (a singular
line) occurs at arbitrary distance rs from Q ?
The condition for the electric field is:
Q rs Q1 rs − r1 Q2 rs − r2
+ + = 0; |rs | = rs
4πε rs3 4πε |rs − r1 |3 4πε |rs − r2 |3
18 2 Electrostatic Fields
Fig. 2.3–1: Equipotential lines and electric lines of force of the three point charges,
|r2 | = 2|r1 |
top: Q1 /Q = −1/4; Q2 /Q = 1/2
bottom: Q1 /Q = 2; Q2 /Q = −2/3
2.4 Force on a Point Charge by the Field of a Space Charge 19
r s
r 2
r
Q 1
Q 1 Q 2
Building the vector and scalar product with rs brings out the equations
and
Q r2 − r1 rs r2 − r2rs
+ Q1 s 3
+ Q2 s = 0.
rs |rs − r1 | |rs − r2 |3
Hence we get the solution:
rs − r1 3 3
|rs − r2 |3 r1
Q1 = −Q2
r2 = Q |rs − r1 | r2
; Q2 = Q .
r1 rs − r2 3
rs r1 − r2 rs3 r2 − r1
A hollow spherical charge with an inner and outer radius of a and b is placed in front
of the conducting plane z = 0 with its center at z = c > b. The charge density is
[As/m3 ]. On the axis of rotation, normal to the conducting plane and through the
center of the sphere, a point charge Q is positioned at z = d with c − b < d < c − a.
The permittivity ε is constant.
If there is no conducting plane, then the force F1 = ez F1 on the charge Q is given by
4 3
1 1
F1 = Q E1 = −Q π (c − d) − a3
3 4πε (c − d)2
3
Q a
= − (c − d) 1 − .
3ε c−d
Thus the force on the charge Q caused by the hollow spherical charge density is equiv-
alent to that of the total charge within the radii c − d and a located in the center z = c.
The remaining charge between the radii b and c − d does not exert a force on the point
charge.
20 2 Electrostatic Fields
"
d c z
a
F = ez (F1 + F2 ) .
Calculate the maximum of the charge density on the cylinder, if it is on the potential
V0 . The permittivity is ε.
For x > 0 outside the cylinder the field is equivalently described by two line charges
±λq of infinite extension at positions (x = ±c; y = 0).
2.5 Charge Density on a Conducting Cylinder in Front of a Conducting Plane 21
y "
¾ m a x
V 0
x = -c x = c x = b
-¸ q ¸ q
x
a
V = 0
⇒ b2 − (a − c)2 = (a + c)2 − b2 ; c2 = b 2 − a2
¾
" V 0
[ m1 ]
1 2 0 0
b { a = { 0 .3 3 3 ; 0 .5 ; 1 ; 2 } [c m ]
1 0 0 0
8 0 0
b { a
6 0 0
4 0 0
2 0 0
0
0 0 .2 0 .4 0 .6 0 .8 a = b 1
Fig. 2.5–1: Maximum value of σ in dependence on the cylinder radius a (the distance
b − a is constant)
∂V λq 1 1 λq 2c
= + =
∂x x=b−a;y=0 2πε c − x c + x x=b−a 2πε c − (b − a)2
2
√ √
λq b 2 − a2 λq b2 − a2 λq 1 + a/b
σ = = = .
π (b2 − a2 ) − (b − a)2 π 2a(b − a) 2πa 1 − a/b
What is the total charge of the conducting surfaces? On what terms does the total
charge of the surface r = b vanish?
In r ≤ a the potential is constant and the electric field vanishes. Outside of the surface
r = b the potential is given in terms of the reciprocal distance 1/r (compare with the
potential of a point charge or an homogeneous spherical charge). For this reasons the
2.6 Potential of Concentric Spheres 23
V 2
V 1
" 1
b
" 2
V (r = a) = V1 = C1 + C2 ; V (r = b) = V2 = C1 a/b + C2
V1 − V2 V1 − V2
C1 = ; C2 = V1 −
1 − a/b 1 − a/b
Thus the potential in a ≤ r ≤ b is given by
1 − a/r
V (r) = V1 − (V1 − V2 ) .
1 − a/b
In the following σ1 is the charge density on the sphere r = a. σ2i and σ2a are the inner
(r < b) and outer (r > b) charge densities on the conducting sphere r = b.
∂V
σ1 = −ε1
∂r a≤r≤b,r→a
∂V ∂V
σ2i = ε1 ; σ2a = −ε2
∂r a≤r≤b,r→b ∂r r>b,r→b
C1 a V2
σ1 = ε1 ; σ2i = −ε1 C1 ; σ2a = ε2
a b2 b
24 2 Electrostatic Fields
In this case the potential V2 relates to the value at points r = b, when only the sphere
r = a with potential V1 is specified (V (r) = V1 a/r for r ≥ a).
z
#
" a
" 0
Solutions are the functions V (r, u) = R(r)Θ(u) with Rn (r) = {rn ; r−(n+1) } and
Θn (u) = {Pn (u) ; Qn (u)}; n = 0, 1, 2, 3, . . . . The functions Pn (u) and Qn (u) are the
Legendre polynomials.
As the exciting potential
M Pn (u) M cos ϑ M a 2
Ve = = = cos ϑ
4πε r(n+1) n=1 4πε r2 4πεa2 r
is a solution for n = 1, the overall solution is given by
⎧ ⎫
M ⎨ (a/r)2 + C1 r/a ⎬ r≤a
V (r, ϑ) = cos ϑ ; .
4πεa2 ⎩ C2 (a/r)2 ⎭ r≥a
The constants C1 and C2 result from the boundary conditions for the potential and
the normal component of the electric flux density.
∂V ∂V
V (r, ϑ)| r>a = V (r, ϑ)| r<a ; ε0 = ε
r→a r→a ∂r r>a ∂r r<a
r→a r→a
1 + C1 = C2 ; C1 − 2 = −2 ε0 /ε C2
3 1 − ε0 /ε
C2 = ; C1 = 2
1 + 2 ε0 /ε 1 + 2ε0 /ε
= −grad V .
For the electric field it follows with E
z = a
0 .5
0
0 0 .5 1 x = a
Fig. 2.7–1: Equipotential lines and electric lines of force for ε/ε0 = 10
M " #
V (r, ϑ) = (a/r)2 − (r/a) cos ϑ ; r ≤ a.
4πεa2
Finally the surface charge σ is
∂Vc ∂Vc
σ = n Dc r>a − Dc r<a = −ε0 +ε =
r→a r→a
∂r r>a ∂r r<a
r→a r→a
ε + 2ε0 M 3M
= − 3
C2 cos ϑ = − cos ϑ .
4πε a 4πa3
exp(−αr)
V (r) = C .
r
α, C are positive constants.
What is the charge density according to the potential in homogeneous space of permit-
tivity ε ?
As the spheres r = const are equipotential surfaces, the direction of the electric field
is radial with
dV exp(−αr)
E(r) = er E(r) = −er = er C(1 + αr) .
dr r2
The space charge density satisfies the Laplace–Poisson equation:
d2 V 2 dV
ΔV = + = − ; r>0
dr2 r dr ε
dV exp(−αr)
= −C (1 + αr)
dr r2
d2 V exp(−αr)
= −C [αr + (1 + αr) (−αr − 2)]
dr2 r3
exp(−αr)
− = C (αr)2
ε r3
exp(−αr)
(r) = −ε C (αr)2 ; r > 0.
r3
Obviously this solution is not complete, because the space charge is negative in r > 0
although the whole potential is positive.
28 2 Electrostatic Fields
The calculation of the charge Qk inside a sphere of radius rk and surface ak = 4πrk2
results in
Qk (rk ) = Ψe (rk ) = da = 4πε r2 C (1 + αrk ) exp(−αrk )
D k
rk2
ak
Apparently a point charge Q0 is located in the center of the coordinate system, which
is not included in the calculation of the Laplace operator defined for r > 0.
If the charge Qk is derived by means of the integral of the incomplete solution over
a sphere with radius r = rk and volume vk we get
rk rk
2 2
Qk (rk ) = d v = 4π (r) r dr = −4πε C α r exp(−αr) dr
vk 0 0
rk
r exp(−αr) 1 exp(−αr)
= −Q0 α2 +
−α α −α 0
As the limit rk → 0 is zero in this case, the complete solution for the charge density
reads
exp(−αr)
(r) = Q0 δ(r) − ε C (αr)2 ,
r3
where the Dirac delta function is defined by
δ(r) = 0 with r > 0 and δ(r) dv = 1 .
v
A sphere with radius r = a has the permittivity ε1 and is surrounded by material with
a permittivity ε2 . A homogeneous line charge of density λq0 is positioned on the z-axis
at points |zq | < c < a inside the sphere or bounded by the interval a < b < zq < d
outside the sphere.
Calculate the potential and the field for all points r ≥ 0 respectively.
¸ q
P
b
# r p
¸ q
" 2
a " 1
-c
For 2
rp
+ 2 rp cos ϑ < 1
zq zq
the following expansion holds.
$
2
1 1 1 rp rp
= 1− − 2 cos ϑ
R zq 2 zq zq
2
2
1 3 rp rp
+ − 2 cos ϑ
2 4 zq zq
3 ⎫
1 3 5 rp
2
rp ⎬
− − 2 cos ϑ + − . . .
2 4 6 zq zq ⎭
P0 = 1,
P1 = cos ϑ ,
1
P2 = (3 cos2 ϑ − 1) ,
2
cos ϑ
P3 = (5 cos2 ϑ − 3) . . . .
2
They satisfy the orthogonality relation
⎧
+1 ⎨ 0 n = k
Pn (u)Pk (u)du = 2 .
⎩ n=k
−1 2n + 1
With the Legendre polynomials the reciprocal distance 1/R takes the form
⎧ n ⎫
⎪
⎪ 1 rp ⎪
⎪ rp ≤ zq
∞ ⎨⎪ ⎪
⎬
1 zq zq
= n Pn (cos ϑ) ;
R ⎪ 1 ⎪
n=0 ⎪ ⎪
⎪ zq ⎪
⎩ ⎭ zq ≤ rp ,
rp rp
where an analog expansion for zq ≤ rp has been used.
2.9 Dielectric Sphere Exposed to the Field of an Axial Line Source 31
For the interpretation of the reciprocal distance we can employ a Taylor expansion
around the point zq = 0:
1 " 2 #−1/2
= rp + zq2 − 2rp zq cos ϑ = f (zq )
R
f (zq )
= f (0) + f (zq ) zq + zq2 + . . . .
zq =0 2!
zq =0
Thus:
1 1 ∂ 1 zq2 ∂ 2 1 (−1) n ∂ n 1
= − zq + − ...+ zqn + ... .
R rp ∂zp rp 2! ∂zp2 rp n! ∂zpn rp
This result can be used to express the potential of the line charge λq (zq ) with |zq | < c
in the following way:
c
1 λq (zq )
V (rp , ϑ) = dzq
4πε0 R(zq , rp , ϑ)
−c
c ∞ n
1 1 zq
rp > c : V (rp , ϑ) = Pn (cos ϑ) λq (zq )dzq .
4πε0 rp n=0 rp
−c & '( )
1
R(zq , rp , ϑ)
Now if we evaluate the integral for the first term of the series we get
c
1 1 Q 1
V0 (rp ) = λq (zq )dzq =
4πε0 rp 4πε rp
−c
c
with Q = λq (zq )dzq .
−c
32 2 Electrostatic Fields
This is the result for the potential V0 (rp ) of the total charge Q located at the center
of the coordinate system. The expression for n = 1 yields
c
1 zq m1 cos ϑ
V1 (rp , ϑ) = λq (zq )P1 (cos ϑ)dzq =
4πε0 rp2 4πε0 rp2
−c
c
with m1 = zq λ(zq ) dzq
−c
and describes the potential of a dipole moment m1 . The series terms with n ≥ 2 define
axial multipoles of higher order
c
mn = zqn λq (zq )dzq .
−c
Their potentials
mn Pn (cos ϑ)
Vn =
4πε0 rpn+1
are in total ∞
V = Vn
n=0
the potential of the line charge.
Example: Homogeneous line charge:
λq0 = Q/2c
c
Q Q cn+1 − (−c)n+1
mn = zqn dzq =
2c 2c n+1
−c
Qc2k
m2k+1 = 0; m2k = ; k = 0, 1, 2, . . .
⎧ 2k + 1
⎨ 0 k>0
lim m2k =
c→0 ⎩ Q k=0
The initial problems are now solved by means of the expansion of the reciprocal dis-
tance.
A homogeneous line charge λq0 of range |zq | < c < a is given in the domain r < a
inside the sphere.
2.9 Dielectric Sphere Exposed to the Field of an Axial Line Source 33
The exciting potential Ve (r, ϑ) of the line charge in homogeneous space of permittivity
ε1 is
∞
∞
mn Pn (cos ϑ)
Ve (r, ϑ) = Vn = ; r ≥ c.
n=0 n=0
4πε 1 rn+1
An analog expression holds for r ≤ c, but is not directly needed for the following
calculation.
c2k+1
m2k = 2λq0 ; m2k+1 = 0
2k + 1
∞
λq0 1 c 2k+1
Ve (r, ϑ) = P2k (cos ϑ) ; r≥c
2πε1 2k + 1 r
k=0
V = Ve (r, ϑ) + Vs (r, ϑ) ; r ≥ c.
This expression already satisfies the continuity of the potential in r = a. The constants
Ak are determined by evaluating the boundary condition for the normal component of
the electric flux density.
∂V ∂V
ε1 − ε = 0
∂r r>a r→a
2
∂r r<a
2k + 1 c 2k+2 1 ε2 2k + 1 c 2k+2 2k + 1
− + Ak 2k = − − Ak
c a a ε1 c a a
In r > a the potential is independent of a. The limit case ε1 = ε2 provides directly the
sole exciting potential. For ε2 → ∞ it follows
c 2k+1
Ak |ε2 →∞ = − ; c≤r≤a
a
34 2 Electrostatic Fields
∞
λq0 1 c 2k+1 c 2k+1 r 2k
V |ε2 →∞ = −
P2k (cos ϑ)
2πε1 2k + 1 r '( a ) &a
k=0
2k+1
a/r cr/a2
a
a
= Ve (r, ϑ) − Ve a2 /r, ϑ = Ve (r) − Ve (r(a/r)2 ) .
r r
The resulting potential is derived directly from the exciting potential (Law of potential
mirroring).
∞ d
λq0 −(n+1)
r≥d: Ve = r zqn dzq Pn (cos ϑ)
4πε2 n=0
b
∞
n+1 n+1
λq0 1 b d
= 1− Pn (cos ϑ)
4πε2 n=0 n + 1 d r
2.9 Dielectric Sphere Exposed to the Field of an Axial Line Source 35
z = a
-1
-1 0 1 x = a
Fig. 2.9–1: Electric lines of force for b/a = 1.2, d/a = 2.2, and ε1 /ε2 = 10
⎧ ⎫
λq0
∞
⎨ (r/a)n ⎬ r≤a
V (r, ϑ) = Ve + Vs = Ve + Bn Pn (cos ϑ); .
4πε2 ⎩ (a/r)n+1 ⎭ r≥a
n=0
z = a
-1
-1 0 1 x = a
Fig. 2.9–2: Electric lines of force for b/a = 1.2, d/a = 2.2, and ε1 /ε2 = 0.5
∂V ∂V
ε1 − ε2 = 0
∂r r<a
∂r r>a r→a
1 − ε1 /ε2 a n n
b
Bn = 1− .
1 + n(1 + ε1 /ε2 ) b d
2.10 Concentric Cylinders With Given Potential 37
Once again the limit case ε1 = ε2 brings out the exciting potential. For large ε1 we get
n
1 a n b
Bn |ε1 →∞ = − 1−
n b d
and
n
λq0 1 a n
∞
b a n+1
V |ε1 →∞ ; r≥a = Ve − 1− Pn (cos ϑ)
4πε2 n=0 n b d r
& '( )
n n
1 a a2 /r b
1−
n r b d
a a
= Ve (r, ϑ) − Ve (a2 /r, ϑ) = Ve (r) − Ve (r(a/r)2 ) .
r r
Thus again the resulting potential is derived directly from the exciting potential.
Two concentric cylinders with radii = a and = b > a are on the potentials
V ( = a, z) = V0 sin(πz/h) and V ( = b, z) = V0 sin(2πz/h). The planes z = 0 and
z = h are equipotential surfaces V ( , z = 0) = 0 and V ( , z = h) = 0. The permittivity
ε is constant.
Calculate the potential V ( , z) in the domain (a ≤ ≤ b; 0 ≤ z ≤ h) and determine
the charge density σ( ) in the (conducting) plane z = 0 within a < < b.
z
V = 0 V 0 s in ( ¼ z = h )
h
V 0 s in ( 2 ¼ z = h )
"
a b %
V = 0
Because the specified potentials on the cylinders = const depend on the z-coordinate,
we have to choose orthogonal functions in this coordinate, e.g. trigonometric functions
38 2 Electrostatic Fields
z = h
0 .8
0 .6
0 .4
0 .2
0
1 1 .2 1 .4 1 .6 1 .8 % = h 2
The given potential on the surfaces = a and = b requires the parameters m = π/h
and m = 2π/h.
that obeys the property Rn (η, η) = 0, the solution for the potential reads
R1 ( , b) R2 ( , a) a≤ ≤b
V ( , z) = V0 sin(πz/h) + sin(2πz/h) ;
R1 (a, b) R2 (b, a) 0 ≤ z ≤ h.
The center of a spherical conducting surface of radius a and potential V0 lies in the
center of the coordinate system, where in addition a point charge Q is positioned.
Within the plane z = c > a a circular surface charge σ of radius b is centered around
the z-axis. The permittivity ε is constant.
Calculate the potential on the z-axis. What is the total charge of the sphere, when the
potential is V0 = 0?
z
2 b
c ¾
q
"
r q
"
r q( a = r q) 2
V 0
The fields inside and outside the conducting sphere are determined independently.
40 2 Electrostatic Fields
Inside the sphere |r| ≤ a the potential is the sum of a constant potential and the
potential of a point charge. As the potential on the surface |r| = a is V0 , the solution
is
Q 1 1
V (r) = V0 + − ; r≤a
4πε r a
and on the z-axis applies
Q 1 1
V (z) = V0 + − ; |z| ≤ a .
4πε |z| a
For the determination of the potential outside the sphere the method of images for
conducting spheres is used. However, the mirrored surface charge is inhomogeneously
distributed on a curvilinear surface. Its contribution to the potential leads to a com-
plicated integration. A considerably easier approach is to derive this contribution by
evaluating the expression for the potential of the initial charge at mirrored observation
points. Thus no additional calculation is needed (cp. problem 2.9).
For example, if a point charge Q is located in front of a grounded sphere of radius a at
position rq with |rq | > a, then the resulting potential V (rp ) at points rp with |rp | ≥ a
is the sum of the potential Ve (rp ) of the charge Q in homogeneous space and of the
potential of the mirrored charge −Qa/|rq | at position rq (a/rq )2 .
We get the same result if the total potential V (rp ) is derived only by means of the
exciting potential.
a
V (rp ) = Ve (rp ) − Ve (rp (a/rp )2 ) ; |rp | > a
rp
This approach is even independent of the initial charge distribution, and hence it is
preferable, if the charge distribution does not comprise simple point charges. This is
actually the case.
In order to apply this approach for the determination of the potential on the rotational
z-axis of the surface charge in front of the conducting sphere, we first have to calculate
the exciting potential Ve (z) of the circular surface charge in the homogeneous space.
1 σ(rq )
Ve (z) = da =
4πε |rp − rq |
b 2π b
σ q dϕq d q σ q d q
= 0 = 0
4πε 2q + (z − c)2 2ε 2q + (z − c)2
q =0 ϕq =0 q =0
σ 0 2 b σ 2
= q + (z − c)2 = b + (z − c)2 − |z − c|
2ε 0 2ε
2.12 Rectangular Cylinder with Given Potential 41
V
V 0
1 .7 5
1 .5
b = a
1 .2 5
0 .7 5
0 .5
b = a = { 0 .5 ; 1 ; 2 ; 3 }
0 .2 5
0
0 1 2 3 4 z = a 5
Fig. 2.11–1: Potential on the z-axis for Q = 10−12 [As], ε = ε0 , σ = 10−11 [As/m2 ],
and c/a = 3
Now with the presented method of images (for potentials) the resulting potential of
our problem is
a a
V (z) = Ve (z) − Ve (a2 /z) + V0 ; |z| ≥ a .
|z| |z|
The last term accounts for the given potential V0 on the surface of the sphere. This
term vanishes, i.e. V (a) = V0 = 0, if the additional total charge of Qz = −4πεaV0 is
brought onto the spherical surface.
The potentials on the boundary of a cylinder with rectangular cross section are
V (x = 0, y) = V (x = a, y) = 0 ; V (y = 0, x) = V2 cos(πx/a)
V (y = b, x) = V1 sin(πx/a) .
42 2 Electrostatic Fields
y
V 1 s in ( ¼ x = a )
V = 0
V = 0 "
a x
V 2 c o s(¼ x = a )
As the given potentials are functions of the x-coordinate, the solutions of Laplace’s
equation ΔV (x, y) = 0 require orthogonal functions for the description of the x-
dependency, accordingly trigonometric functions. Due to the homogeneous boundary
conditions in x = 0 and x = a the functions sin(mπx/a) have to be considered. In
consequence the dependency of the y-coordinate is given by hyperbolic or exponential
functions. With these functions a possible ansatz for the potential inside the rectangle
(0 ≤ x ≤ a; 0 ≤ y ≤ b) is
∞
V (x, y) = [Am sinh(mπy/a) + Bm cosh(mπy/a)] sin(mπx/a) .
m=1
It is useful to split the problem, because of the simple description of the boundary
values. For the first subproblem with solution VI only the boundary values at y = b
exist, whereas in the second subproblem with solution VII only the boundary values
at y = 0 are considered. All other boundaries have the potential V = 0, respectively.
The solution of the original problem is V = VI + VII .
The potential on the surface y = b suggests, that only the term with m = 1 is needed
in the summation above for the potential VI . With the requirement VI (x, y = 0) = 0
we get the solution
sinh(πy/a)
VI (x, y) = V1 sin(πx/a) .
sinh(πb/a)
Things are different for the second subproblem with solution VII . With the given
potential on the surface y = 0 the potential is discontinuous at x = 0 and x = a. As a
2.12 Rectangular Cylinder with Given Potential 43
y = a
1 .2
0 .8
0 .6
0 .4
0 .2
0
0 0 .2 0 .4 0 .6 0 .8 x = a 1
Here the arguments of the hyperbolic functions have been adapted to satisfy the bound-
ary condition VII (x, y = b) = 0. The division by the hyperbolic function with constant
argument simplifies the evaluation of the boundary condition and thus leads to
∞
VII (x, y = 0) = Cm sin(mπx/a) = V2 cos(πx/a) .
m=1
44 2 Electrostatic Fields
Now, making use of the orthogonality relations for trigonometric functions by multi-
plication with sin(nπx/a) and integration from x = 0 to x = a results in
a a
Cm sin(mπx/a) sin(nπx/a) dx = V2 cos(πx/a) sin(nπx/a) dx
m 0 0
a
cos((n + 1)πx/a) cos((n − 1)πx/a)
Cn a/2 = −V2 +
2(n + 1)π/a 2(n − 1)π/a 0
a 4k
= V2 2
; n = 2k ; k = 1, 2, . . .
π 4k − 1
8 k
⇒ C2k = V2 .
π 4k 2 − 1
Finally the solution of the subproblem is
∞
8 k sinh(2kπ(b − y)/a)
VII (x, y) = V2 sin(2kπx/a)
π 4k 2 − 1 sinh(2kπb/a)
k=1
A sphere of radius a is intersected by the plane z = 0 running through its center. The
two hemispheres z > 0 and z < 0 carry the space charge densities ± . The permittivity
ε is constant.
Calculate the potential and the field on the z–axis. Simplify the expression for the field
for the case z a and for z = 0.
The determination of the potential on the z-axis requires the solution of
1 (rq )
V (z) = dv ; V (−z) = −V (z) ; |ez z − rq | = r.
4πε |ez z − rq |
"
+ %
- %
a
a π/2 2π
rq n
∞
1
V (z) = [Pn (cos ϑq ) +
4π ε z n=0 z
rq =0 ϑq =0 ϕq =0
∞ 1
1 1 an+3
= [Pn (u) − Pn (−u)] du ; u = cos ϑq .
2ε z n=0 z n n + 3
u=0
n Pn+1 − Pn−1
With Pn (−u) = (−1) Pn (u) and Pn (u)du = it follows
2n + 1
1
a2 1 a 2(k+1)
∞
P2k+2 − P2k
V (z) = ; n = 2k + 1 .
2ε k+2 z 4k + 3 0
k=0
For points 0 < z < a the integration has to be split up because of the expansion of the
reciprocal distance 1/r.
⎡ ⎤
∞
z a n
⎢ 1 rq n 1 z ⎥
V (z) = ⎣ rq2 drq + rq2 drq ⎦ ·
2ε n=0 z z rq rq
rq =0 z
1
· [Pn (u) − Pn (−u)] du
0
∞
z 2k+1 P (0)
z2 a2 z 2 2k
= + −
2ε 2(k + 2) 2k − 1 a a (k + 1)
k=0
∞
a 2 4k + 3 z 2 z 2k+1 P2k (0)
V (z) = −
2ε 2(k + 2) a a (k + 1)(2k − 1)
k=0
holds, which is a continuous function for z → a, thus takes the value of Ez , that has
been derived from the potential in z > a.
The expression for the field simplifies for z → 0 to
a
Ez (z = 0) = − ,
2ε
because only the series term k = 0 contributes to the field. This result can also be
obtained by a direct integration for the point z = 0.
E = 1 r
3 dv = ez Ez ; r = −rq
4πε r
1 ez (−rq )
Ez = dv ; dv = rq2 sin ϑq drq dϑq dϕq
4πε rq3
⎡ 0 ⎤
a 1
⎣− u du + u du⎦ drq ;
= − u = cos ϑq
2ε
rq =0 −1 0
a
Ez = −
2ε
2.13 Potential of Hemispherical Charge Distributions 47
z = a
-1
-1 0 1 x = a
Fig. 2.13–1: Equipotential lines and electric lines of force inside and outside of the
charge distribution (In the charged area the flux inside a “flux tube” is
no longer constant like in the charge-free space, but decreases until the
electric field is zero.)
For great distances z |a| it is possible to calculate the potential by means of the
positive and negative charge center. The (positive) charge center of the hemispherical
charge in z > 0 is (Q is the total charge):
π/2 a
1 3
rqs = rq dv = ez rqs ; rqs = 3 ez rq rq2 drq sin ϑq dϑq
Q a
ϑq =0 rq =0
48 2 Electrostatic Fields
1
3 3
u = cos ϑq ; rqs = a udu = a.
4 8
0
6
Hence the distance of the the positive and negative charge center is s = 8a and for
points z |a| the potential Vs is defined by the dipole moment
4 3 6 π
m
= ez m ; m = Q·s = πa · a = a4 .
6 8 2
The potential on the z-axis is
a4
Vs = sign(z) ; |z| a.
8ε z 2
This expression agrees with the series term k = 0 in the series expansion of V (z a),
which is the dominant term for great distances.
Calculate the charge distributions on the plates. What is the force on the dielectric?
¾ 2 ¾ 1
d
F x
x p b x
2.14 Energy and Force inside a Partially Filled Parallel-Plate Capacitor 49
Before the dielectric has been inserted the charge density on the plates is ±σ0 = ±D0 =
= ±ε0 E0 and the total charge is ±Q0 = ±σ0 b2 . The energy stored in the electric field
W0 is given by the integral of the energy density.
1 1 U2 1 b2
W0 = E0 D0 · b2 d = ε0 20 b2 d = U02 C0 ; C0 = ε0
2 2 d 2 d
After inserting the dielectric the electric flux density D is different in rooms of different
materials and therefore the associated charge density σ is also different. This is because
the electric field is continuous at the boundary x = xp , thus the electric field is identical
everywhere between the plates. Its value depends on the position of the dielectric and
therewith on the value xp . As the total charge on the plates remains the same, we get
the following result with two different charge densities σ1 and σ2 for the domains with
permittivities ε0 and ε.
σ1 (b − xp ) + σ2 xp = σ0 b ; E1 = E2 = E
ε0 ε E1 = ε0 ε E2 = ε D1 = ε σ1 = ε0 D2 = ε0 σ2
ε σ ε0
σ2 = σ1 ; σ1 = 0 = σ2 = D1
ε0 xp ε ε
1+ −1
b ε0
ε0 E0 ε
D1 = ε0 E = = 0 D2
xp ε ε
1+ −1
b ε0
E
E = 0
xp ε
1+ −1
b ε0
The limits xp = 0, b and ε = ε0 provide the correct results.
1 1
W (xp ) = ε0 E 2 db(b − xp ) + ε E 2 dbxp
2 2
1 2 1 xp ε
= E db [ε0 (b − xp ) + εxp ] = ε0 E 2 db2 1 + −1
2 2 b ε0
W0
= .
xp ε
1+ −1
b ε0
50 2 Electrostatic Fields
3
¾ 2= ¾ 0
F x = F 0
1
¾ 1= ¾ 0; W = W 0
0
0 0 .2 0 .4 0 .6 0 .8 1
x p= b
Fig. 2.14–1: Charges σ1,2 /σ0 and force Fx /F0 in dependence on the filling ratio xp /b
with ε/ε0 = 5 ; F0 = W0 /b
With this result and the relation F = −grad W the solution for the x–directed force
Fx acting on the dielectric becomes
ε
W0 /b −1
∂W ε0
Fx = − = 2 > 0 ; F0 = W0 /b.
∂xp xp ε
1+ −1
b ε0
The drawing below shows the problem with a given potential V0 (x) on the boundary
y = 0 in dependence on the x-coordinate, whereas all other boundaries own the poten-
tial V = 0.
Calculate the potential inside the closed boundary.
2.15 2D-Problem with Homogeneous Boundary Conditions 51
y
V = 0
d
V = 0
V = 0
V = 0
0
0 b a x
V 0(x )
It is useful to split the computation area by the plane y = c, because of the given
potential in dependence on the x-coordinate and the homogeneous boundary values
V = 0 on x = (0; a) and x = (0; b). The potentials V1,2 in the domains 0 ≤ y ≤ c
and c ≤ y ≤ d satisfy Laplace’s equation ΔV1,2 = 0. The functions which satisfy
the homogeneous boundary conditions are trigonometric functions sin(p1,2,n x) with
p1n = nπ/a and p2n = nπ/b. For the y-coordinate the hyperbolic functions are chosen,
so that in the area c ≤ y ≤ d the potential V2 (x, y) at y = d vanishes and in the
area 0 ≤ y ≤ c the potential V1 (x, y) is a combination of two linear independent parts,
where the first has zeros at y = 0 and the second at y = c for b < x < a.
∞
sinh(p1n y) sinh(p1n (y − c))
V1 (x, y) = An + Bn sin(p1n x) ; 0≤y≤c
n=1
sinh(p1n c) sinh(p1n (−c))
∞
sinh(p2n (y − d))
V2 (x, y) = Cn sin(p2n x) ; c≤y≤d
n=1
sinh(p2n (c − d))
a
2
Bn = V0 (x) sin(p1n x) dx .
a
0
52 2 Electrostatic Fields
∞
Bn
⇒ An coth(p1n c) − p1n sin(p1n x) =
n=1
sinh(p1n c)
∞
= Cn p2n sin(p2n x) coth(p2n (c − d)) .
n=1
∞
2 Bn
Cm = p1n αmn An coth(p1n c) − .
mπ coth(p2m (c − d)) n=1 sinh(p1n c)
The potential and the field are determined, when both systems of equations are solved.
0
0 b a
x
Fig. 2.15–1: Equipotential lines and electric lines of force for b/a = 2/3, c/a = 1/3,
d/a = 2/3, and V0 = 1 [V]
The plane z = 0 divides the space into two half-spaces z > 0 of permittivity ε1 and
z < 0 of permittivity ε2 . The point rq1 in z > 0 marks the position of a point charge
∗
Q1 and the point rq2 = rq1 − 2ez (rq1ez ) = rq1 in z < 0 marks a second point charge
Q2 .
Q 1
z P
r q 1
" e
1 r p z
" 2
r q 2
Q 2
If only the point charge Q1 in z > 0 is considered, then the method of images for
dielectric half-spaces yields the potential
⎧
⎪
⎪
1
−
k
z ≥ 0; k = (ε2 − ε1 )/(ε2 + ε1 )
Q1 ⎨ |rp − rq1 | |rp − rq1
∗ |;
V (rp ) =
4πε1 ⎪
⎪ 1−k ∗
⎩ ; z ≤ 0; rq1,2 = rq1,2 − 2ez (ez rq1,2 ).
|rp − rq1 |
Accordingly in the presence of the second point charge Q2 in the half-space z < 0 the
total solution for the potential is
Q1 1 k Q2 1+k
V (rp ) = − ∗ | + 4πε | ; z≥0
4πε1 |rp − rq1 | |rp − rq1 2 rp −
rq2 |
Q1 1−k Q2 1 k
V (rp ) = + + ∗ | ; z ≤ 0.
4πε1 |rp − rq1 | 4πε2 |rp − rq2 | |rp − rq2
We can proof the result for the case ε1 = ε2 = ε, thus k = 0, and get the expected
result
Q1 1 Q2 1
V (rp ) = + .
4πε |rp − rq1 | 4πε |rp − rq2 |
With ε2 ε1 , thus k → 1, we get
Q1 1 1
V (rp ) = − ∗ | ; z > 0; V = 0; z < 0.
4πε1 |rp − rq1 | |rp − rq1
This result agrees with the method of images for the conducting half-space. In the
same way with ε1 ε2 and k → −1 we get
Q2 1 1
V = 0 ; z > 0 ; V (rp ) = − ∗ | ; z < 0.
4πε2 |rp − rq2 | |rp − rq2
2.16 Method of Images for Dielectric Half-Spaces 55
z
z q 1
-2
-2 0 2
x = z q 1
Fig. 2.16–1: Equipotential lines and electric lines of force for Q1 /Q2 = −2,
ε1 /ε2 = 1/3, and zq2 = −zq1
The forces F1,2 = ez F1,2 acting on the charges Q1,2 are
ε1 ε2
−kQ1 + (1 + k)Q2 −kQ2 − (1 − k)Q1
Q1 ε2 Q2 ε1
F1 = ; F2 = .
4πε1 (z1 − z2 )2 4πε2 (z1 − z2 )2
56 2 Electrostatic Fields
z
z q 1
-2
-2 0 2
x = z q 1
Fig. 2.16–2: Equipotential lines and electric lines of force for Q1 /Q2 = 2,
ε1 /ε2 = 1/3, and zq2 = −zq1
The cylinder of radius = b has the potential V1 . Within this cylinder a second
cylinder of radius = a < b has the potentials
Find the electric field in the whole space. The permittivity ε is constant.
2.17 Concentric Cylinders with Given Potentials 57
"
"
V 1
V 0
%
"
'
a
b
-V 0
B0 + C0 = 0 ;
π
An π = (Bn + Cn )π = 2V0 sin(nϕ)dϕ =
0
2V0 4V0
= [1 − cos(nπ)] = ; n = 2k + 1 ; k = 0, 1, 2, . . .
n n
58 2 Electrostatic Fields
y = a
-1
-2
-2 -1 0 1 2
x = a
Fig. 2.17–1: Equipotential lines and electric lines of force for b/a = 2 and V1 /V0 = 0.5
4V0
⇒ An = = Bn + Cn ; n = 2k + 1 .
nπ
On the cylinder = b applies
∞
ln(b)
B0 + C0 + [Bn (b/a)n + Cn (a/b)n ] sin(nϕ) = V1
ln(a) n=1
and hence
4V0
Bn = ; Cn = −Bn (b/a)2n .
nπ[1 − (b/a)2n ]
is
The electric field E
with n = 2k + 1.
The ansatz with the zeros x0r of the Bessel functions satisfies the boundary condition
VE ( = a, z) = 0 and the exponential function describes the decay of the potential away
from the charge. It is useful to express the ring charge in terms of a surface charge
60 2 Electrostatic Fields
¸ q
c
V = 0
" 0
b
a %
"
σ( ) in the plane z = c, where the charge density is zero except for points on the ring
( = b; z = c). With the charge σ( ) the boundary condition for the z-component of
in z = c leads to
the electric flux density D
ez D z>c − D z<c = σ( )
z→c z→c
∂V ∂V
σ( )
⇔ − z>c = ε0
∂z z<c ∂z
z→c z→c
∞
2 1
⇒ x0r Ar J0 (x0r /a) = σ( ) .
a r=1 ε0
The Bessel functions J0 (x0r /a) belong to an orthogonal system of functions, thus
the multiplication with an orthogonal function J0 (x0s /a) and the weight function
results after integration from = 0 to = a in
a
2 1 1
x0s As a2 J12 (x0s ) = σ( ) J0 (x0s /a) d ,
a 2 ε0
0
where the orthogonality relation for Bessel functions has been applied. For the remain-
ing integral the integration interval restricts to the location of the ring charge = b ± δ
2.18 Force on a Ring Charge inside a Conducting Cylinder 61
and we get
b+δ
1
Ar = σ( ) J0 (x0r /a) d
ε0 x0r aJ12 (x0r )
b−δ
b+δ
J0 (x0r b/a)b
= σ( )d with δ→0
ε0 ax0r J12 (x0r )
b−δ
& '( )
λq
When the space z < 0 has the permittivity ε, then the variation in the field results
from the difference between the exciting potential and the resulting potential V ( , z).
∞
J0 (x0r b/a)
V ( , z) = VE ( , z) + V0 Br J0 (x0r /a) exp(−x0r |z|/a)
r=1
x0r J12 (x0r )
The additional potential in this approach uses the solution for the exciting potential,
but applies different constants Br . The total approach already satisfies all boundary
conditions for the potential and also for the electric flux density in z = c. Finally
evaluating the continuity relation of the flux density in the boundary layer of different
permittivities z = 0 leads to the constants Br .
∂V ∂V
ε0 = ε
∂z z>0 ∂z z<0
z→0 z→0
∞
J0 (x0r b/a)
V ( , z) = VE ( , z) − kV0 J0 (x0r /a)·
x J 2 (x0r )
r=1 0r 1
⎧
⎨ exp(−x z/a) ; z≥0
0r
· exp(−x0r c/a)
⎩ exp(x0r z/a) ; z≤0
62 2 Electrostatic Fields
0 b a
Fig. 2.18–1: Equipotential lines and electric lines of force for k = 0.9 and a = 2b = 4c
⎧
⎨ V ( , z) − kV ( , −z) ; z≥0
E E
⇒ V ( , z) =
⎩ (1 − k) VE ( , z) ; z≤0
This result agrees with the method of images for dielectric half spaces.
with Vs = V ( , z) − VE ( , z)
2.19 Geometry with Circular Symmetry and Given Potentials on Parallel Planes 63
∞
J 2 (x0r b/a)
0
⇒ F = −2 π b λq k V0 /a exp(−2x0r c/a)
r=1
J12 (x0r )
∞
J 2 (x0r b/a)
0
= −k Q E0 exp(−2x0r c/a)
r=1
J12 (x0r )
∞
exp(−2x0r c/a)
F = −k Q E0 .
b→0
Q = const r=1
J12 (x0r )
If in addition the limit case c → 0 is considered, we get an divergent series and thus
an unbounded force. Within a boundary layer of different permittivities only surface
charges can exist.
a + V 0
V = 0
c %
"
-a -V 0
Because the potential on = c depends on the z-coordinate, the solution of the differ-
ential equation ΔVs = 0 is
∞
I0 (pn )
Vs ( , z) = An sin(pn z) .
n=1
I0 (pn c)
sin(pn a) = 0 ; pn = nπ/a ; n = 1, 2, 3, . . .
and at = c
∞
−V0 z/a = An sin(pn z) .
n=1
2V0
An = (−1)n
pn a
and finally the resulting potential is
∞
z n
(−1) I0 (pn )
V ( , z) = V0 +2 sin(pn z) .
a p a I0 (pn c)
n=1 n
Within the space of permittivity ε0 the cylinder ≤ a is filled with material of per-
mittivity ε. On its surface = a a charge distribution depending on the ϕ-coordinate
is given. ⎧
⎨ σ |ϕ| < ϕ0
0
σ(ϕ) =
⎩ −σ0 ϕ0 < ϕ < 2π − ϕ0
" 0
% ' = ' 0
"
'
¾ 0
a
' = -' 0
-¾ 0
This approach already satisfies the continuity relation of the potential on the cylindrical
surface = a.
The remaining constants follow from the boundary condition for the normal component
on the cylindrical surface.
D of the electric flux density D
∂V ∂V
⇔ ε − ε0 = σ(ϕ)
∂ <a ∂ >a
→a →a
∞
C0 nCn
⇒ −ε0 + (ε + ε0 ) cos(nϕ) = σ(ϕ)
a ln a n=1 a
2π
a ln a σ0 a ln a ϕ0
C0 = − σ(ϕ) dϕ ; C0 = 1−2 .
2πε0 ε0 π
0
2π
nCn σ0 ϕ0 2π−ϕ0
π(ε + ε0 ) = σ(ϕ) cos(nϕ)dϕ = sin(nϕ)|−ϕ 0
− sin(nϕ)| ϕ0
a n
0
4σ0 a sin(nϕ0 )
⇒ Cn =
π(ε + ε0 ) n2
66 2 Electrostatic Fields
y = a
-1
-1 0 1 x = a
Fig. 2.20–1: Electric lines of force and equipotential lines for ϕ0 = π/4 and ε/ε0 = 3/2
⎧ ⎫
∞
n ⎨ ( /a)n−1 ⎬ <a
Eϕ = Cn sin(nϕ) ; .
a ⎩ (a/ )n+1 ⎭ >a
n=1
2.21 Potential and Field of Dipole Layers 67
The sphere r = a carries an electric dipole layer of density m = ±m r/r. In the area
0 ≤ ϑ < ϑ0 the sign of the radially directed dipoles is positive, whereas in the remaining
area ϑ0 < ϑ ≤ π the sign is negative. The absolute value m and the permittivity ε are
constant.
Calculate the potential and the field at points on the z-axis.
r m
r
# 0
"
a r q
q R
# 1
r m
r
It relates to the surface area on the sphere with radius 1 and origin p, which results
from the projection of the boundary contour of the surface aD on the unit sphere.
68 2 Electrostatic Fields
This result shows, that the solid angle depends only on the boundary of the surface
aD and on the orientation of its normal n to the observation point p. Therefore the
potential of the dipole layer m r/r at r = a and 0 ≤ ϑ < ϑ0 on the z-axis in the interval
−∞ < z < a can also be calculated by means of the solid angle of a spherical calotte
with radius R, the same boundary contour like aD , and with an apex angle 2ϑ1 . In
this way the solid angle in z < a is
−z + a cos ϑ0
Ω = − n r/r3 daD = 2π(1 − cos ϑ1 ) = 2π 1 −
R
aD
2 2
R2 = (−z + a cos ϑ0 ) + (a sin ϑ0 ) = z 2 + a2 − 2az cos ϑ0
and the potential becomes
m
V (z) = − [1 + Z0 /R] ; z<a; Z0 = z − a cos ϑ0 .
2ε
For points z > a the orientation of the surface normal to the observation point changes,
thus the solid angel of the dipole layer with coordinates r = a and 0 ≤ ϑ < ϑ0 becomes
The dipole layer on the spherical calotte r = a and ϑ0 < ϑ ≤ π with density −r/r m
analogously provides the potentials
m
V = − [1 + Z0 /R] = −V0 [1 + Z0 /R] /2 ; z < −a
2ε
and
m
V = [1 − Z0 /R] = V0 [1 − Z0 /R] /2 ; z > −a .
2ε
The superposition of the potentials gives the total potential
⎧
⎪
⎪
m ⎨ m/ε ; z>a
V (z) = − Z0 /R + 0 ; −a < z < a
ε ⎪
⎪
⎩
−m/ε ; z < −a .
2.21 Potential and Field of Dipole Layers 69
8
V 0 1
# 0
0 .7 5
0 .5
0 .2 5
# 0
0
-0 .2 5
# 0
-0 .5
# 0 = { n ¼ = 1 2 ; n = 1 , 2 , 3 , 4 , 5 , 6 }
-3 -2 -1 0 1 2 3
z = a
E z
E 0
# 0 = { n ¼ = 1 2 ; n = 1 , 2 , 3 , 4 , 5 , 6 }
3
# 0
0
-3 -2 -1 0 1 2 3
z = a
Due to the rotational symmetry the electric field is z-directed and is derived by
∂V m ∂ Z0 V0
E = −grad V ; Ez (z) = − = ; E0 = ; z = ±a
∂z ε ∂z R a
2
m 1 2 1 a Z0
Ez = − (z − a cos ϑ0 ) = E0 1− ; z = ±a .
ε R R3 R R
The area 0 ≤ ϑ < ϑ0 of the spherical conducting surface r = a has the potential V0 ,
whereas the remaining part ϑ0 < ϑ ≤ π, isolated from the first one, has the Potential
−V0 . The permittivity ε is constant.
Calculate the potential in the whole space and the charge density on the subarea
(r = a; 0 ≤ ϑ < ϑ0 ).
V 0
# r
# 0
"
-V 0
Because of the rotational symmetry with ∂V /∂ϕ = 0 the solution functions of Laplace’s
equation ΔV (r, ϑ, ϕ) = 0 using spherical coordinates (r, ϑ, ϕ) are
⎧ ⎫
∞ ⎨ (r/a)n ⎬ r≤a
V (r, ϑ) = An Pn (cos ϑ) ; .
⎩ (a/r)n+1 ⎭ r≥a
n=0
Since the rotation axis is part of the computation domain the functions Qn (cos ϑ) can
be dropped.
2.22 Sphere with Given Potential 71
With the substitution u = cos ϑ, u0 = cos ϑ0 , and the orthogonality relation for the
Legendre polynomials Pn (u) it follows
⎡ ⎤
∞
1 1 u0
An Pn (u)Pk (u)du = V0 ⎣ Pk (u)du − Pk (u)du⎦
n=0 −1 u0 −1
⎧
2 ⎨ −2u0 ; k=0
Ak = V0 1 u0
2k + 1 ⎩ Pk+1 − Pk−1 − Pk+1 − Pk−1 ; k>0
2k + 1 u0 2k + 1 −1
A0 = −u0 V0
"
An = V0 /2 Pn+1 (1) − Pn−1 (1) + Pn+1 (−1) − Pn−1 (−1)
#
− Pn+1 (u0 ) + Pn−1 (u0 ) − Pn+1 (u0 ) + Pn−1 (u0 ) ; n>0
The inner and outer charge densities on the surface (r = a, u0 < u < 1) follow from
the normal component of the electric flux density.
∂V
σi = −er D =ε
∂V
σa = er D >a = −ε
<a
∂r <a ; ∂r >a
→a →a
→a →a
∞ ∞
ε εV0
σi = An nPn (u) = n [Pn−1 (u0 ) − Pn+1 (u0 )] Pn (u)
a n=0 a n=0
∞ ∞
ε εV0
σa = An (n + 1)Pn (u) = (n + 1) [Pn−1 (u0 ) − Pn+1 (u0 )] Pn (u)
a n=0 a n=0
The formal calculation of the total charge leads to a non-convergent series, because of
the singularity of the electric field at the discontinuity of the potential.
72 2 Electrostatic Fields
z = a
-1
-2 -1 0 1 2 x = a
Fig. 2.22–1: Equipotential lines and electric lines of force for ϑ0 = π/6
"
V 0
-a a x
∞
= −grad V = 2V0 /π
E sin(pa) [ex sin(px) + ey cos(px) sign (y)] exp(−p|y|) dp .
0
74 2 Electrostatic Fields
1 0
¾
¾ 0
8
0
-0 .8 -0 .4 0 0 .4 0 .8
x = a
Fig. 2.23–1: Charge density for the interval |x| < 0.9 a
In this special case the expressions for the integrals are known and it follows
V0
V (x, y) = [arctan((a − x)/|y|) + arctan((a + x)/|y|)] .
π
From this result the charge density σ(x) can easily be deduced.
⎛ ⎞
∂V ∂V V0 1 1
σ(x) = ε ⎝ − ⎠ = 2ε +
∂y y→0
y<0
∂y y→0
y>0
π a−x a+x
x=±a x=±a
2σ0 2εV0
= ; σ0 =
1 − (x/a)2 πa
The integral for the determination of the total charge per unit length l
a
Q/l = σ(x)dx
−a
" 0
a %
"
If the surface charge is located in a plane z = c > 0, then the calculation of the
potential could be done by means of the method of images. In this case an additional
image charge −kσ with k = (ε − ε0 )/(ε + ε0 ) in the plane z = −c must be considered
to determine the potential in z ≥ 0, whereas the permittivity is ε0 in both half-spaces.
The potential for points z < 0 results from a surface charge (1 − k)σ in z = c. For the
present problem the limit case c → 0 holds, thus a surface charge (1 − k)σ in z = 0 is
the substitute in space of homogeneous permittivity ε0 .
Hence the potential on the z-axis is
2
2π a z +a2
1 (1 − k)σ (1 − k)Q d dϕ (1 − k)Q dζ
V (z) = daS = = √
4πε0 r 4π 2 ε0 a2 2
z + 2 2πε0 a2 2 ζ
aS ϕ=0 =0 ζ=z 2
(1 − k)Q 2
2 − |z| .
V (z) = z + a
2πε0 a2
The electric field results from
= −grad V = ez Ez + e E ;
E = 0, z) = ez Ez = −ez ∂V
E(
∂z
76 2 Electrostatic Fields
and becomes
⎡ ⎧ ⎫⎤
∂V (1 − k)Q ⎣ z ⎨ −1 ⎬ z>0
Ez = − = − √ + ⎦;
∂z 2πε0 a2 z 2 + a2 ⎩ +1 ⎭ z<0
⎡ ⎧ ⎫⎤
σ ⎣ z ⎨ −1 ⎬ z>0
= − √ + ⎦; .
ε + ε0 z 2 + a2 ⎩ +1 ⎭ z<0
Q c
" V 0
"
¸ q
%
% = b
Two conducting plates at surfaces (y = ±a ; x < 0) and with potentials ±V0 build a
parallel-plate capacitor.
y
a
"
x
-a
For the calculation of the field one can make use of a conformal transformation in form
of a Schwarz-Christoffel-Mapping. As the plane y = 0 is a symmetry plane with the
potential V = 0, it is sufficient to map the half-space y ≥ 0 with complex coordinates
z = x + jy to the half-space u ≥ 0 with coordinates w = u + jv. Furthermore the
plates of the capacitor are assumed to be of infinite length in negative x-direction.
The drawing below shows the position of three corresponding points in the z- and
w-plane.
j y j v
3
3
1
ja
2 1 2
x u 1 = -1 u 2 = 0 u
If one mapping point in the w-plane is at infinity, then the general mapping function
for a polygon with n-edges and interior angles αλ is
n−1
6
z(w) = C 1 (w − uλ )−αλ dw + C 2 .
λ=1
2.26 Boundary Field of a Parallel-Plate Capacitor 79
where the two arbitrary real numbers u1 and u2 have been set to u1 = −1 and u2 = 0.
For the determination of the constants C 1 and C 2 the correspondence of the first point
yields
z(w = −1) = ja = C 1 [−1 + ln(−1)] + C 2 = C 1 [jπ − 1] + C 2 .
0
1 + exp(jϕ)
−ja = C 1 exp(jϕ) jdϕ = C 1 (−jπ)
exp(jϕ) →0
ϕ=π
a
C 1 = C 2 = a/π ; z(w) = [1 + w + ln(w)]
π
With this, the mapping function is known. The contour in the z-plane is mapped onto
the real axis in the w-plane with the potential V = V0 in u < 0 and the potential V = 0
in u > 0. These potentials excite a field with rotational symmetry in the w-plane and
thus with equipotential surfaces arg(w) = const . This potential is that of a virtual
line source λq in w = 0, but multiplied with j to swap the real and imaginary part of
the complex potential
λq λq ϕ
P e (w) = −j ln(w) with V (u, v) = Re{P e (w)} = .
2πε ε 2π
Now with V (u > 0, v = 0) = 0 and V (u < 0, v = 0) = V0 it follows
V0
P e (w) = −j ln(w) .
π
The inverse transformation of the complex potential P e (w(z)) = P e (z) has to be done
numerically. Like in many other cases the analytic expression of the inverse function
of w(z) is not known.
80 2 Electrostatic Fields
y = a
0
-1 0 1 x = a 2
The complex electric field E(z) = Ex (x, y) + j Ey (x, y) is derived from the potential by
∗ ∗ ∗
dP e dP e dw dP e (w) 1
E(z) = − = − = −
dz dw dz dw (dz/dw)∗
∗
V0 1 1 V0 1
= − −j = −j = Ex + j Ey .
π w a/π(1 + 1/w∗ ) a 1 + w∗ (z)
3. Stationary Current Distributions
The radial component of the electric field E on the surface of a cylinder ≤ a with
conductivity κ is given by
⎧
⎨ cos(kϕ) ; |ϕ| < π/3
E ( = a) = E0 .
⎩ 0 ; π/3 ≤ ϕ ≤ 5π/3
For which values of k is this problem well-defined? Calculate the current distribution
inside the cylinder ≤ a and the dissipated power per unit length.
' = ¼ = 3
%
E % (' )
'
k
k = 0 ' = -¼ = 3
r ) = κE
The stationary current density J( = −κ grad V can be described by a potential
V (r) with ΔV = 0. On the surface of the conducting cylinder the normal derivative
∂V /∂n of the potential is given. As the cylinder is free of sources, the normal derivative
must fulfill
2π 2π π/3
∂V ∂V
adϕ = 0 = adϕ = −E0 cos(kϕ) a dϕ = 0.
∂n =a ∂ =a
ϕ=0 0 ϕ=−π/3
π/3
−mπ Am /a = E0 cos(3pϕ) cos(mϕ)dϕ
−π/3
⎧
⎪
aE0 ⎨ π/3 ; n = 3p
An = − .
nπ ⎩⎪ sin((3p − n)π/3) sin((3p + n)π/3)
3p − n + 3p + n ; n = 3p
The constant A0 remains unknown for boundary value problems of this kind.
3.2 Current Distribution around a Hollow Sphere 83
Calculate the potential and the electric field inside and outside the sphere. What is
the charge density on the surface of the sphere?
∂V0
J0 = −ez J0 = −ez κ E0 = −κ grad V0 = −ez κ
∂z
is described by the potential V0 = E0 z = E0 r cos ϑ.
V = V0 + Vs with Δ Vs = 0
is then given by
⎧
⎪ 1
∞
⎪
⎨ V0 + An (a/r)
n+1
Pn (cos ϑ) ; r ≥ a
n=0
V = 1
∞ .
⎪
⎪ n
⎩ Bn (r/a) Pn (cos ϑ) ; r≤a
n=0
84 3 Stationary Current Distributions
J 0 = -e zJ 0
k = 0 , " 0
As the potential of the initial field only includes the Legendre polynomial P1 = cos ϑ
r
V0 = E0 a P1 (cos ϑ),
a
we can leave out all other polynomials in our solution, thus:
⎧
⎨ V0 + A (a/r)2 cos ϑ ; r ≥ a
V (r, ϑ) = .
⎩
B r/a cos ϑ ; r≤a
At the boundary to the non-conducting sphere the current streamlines are tangential.
Thus the normal derivative of the potential vanishes at r = a.
∂V a E0
= (E0 − 2A/a) cos ϑ = 0 → A = ,
∂r r>a,r→a 2
z = a
-1
-2
-2 -1 0 1 x = a 2
= −grad V | ∂ 3 3
E r<a = −ez E0 (r cos ϑ) = −ez E0 .
∂z 2 2
Thus in r < a the electric field remains homogeneous.
The charge density at r = a results from the normal derivative of the potential V |r<a
∂V ∂V 3
σ = −ε0 = ε = ε0 E0 cos ϑ .
∂n r=a ∂r r=a
0
2
Finally the current density J = −κ grad V in r > a is
" #
J = −κ E0 (er cos ϑ − eϑ sin ϑ) − (a/r)3 (er cos ϑ + 1/2eϑ sin ϑ)
" #
= −J0 ez − (a/r)3 (er cos ϑ + 1/2eϑ sin ϑ) .
86 3 Stationary Current Distributions
Calculate the current distribution, the power loss Pv /l, and the resistance
R/l = Pv /(lI 2 ).
I = l
y
k ® ¥
2 s
d
c a x
2 s
k ® ¥
I = l
On the boundary of the cylinder the normal derivation of the potential ∂V /∂n = 0
vanishes, except for the range of the conducting electrodes, where the potential is
constant because of the perfect conductivity. Hence the problem is a mixed boundary
value problem. However, as 2s b we can assume, that the normal derivative of
the potential at the ends of the electrodes is constant. The normal component of the
current density on the boundary y = b becomes
I/l ∂V ∂V ∂V I/l
Jn = = −κ = −κ ; = − ; |x − c| < s
2s ∂n ∂y ∂y y=b 2κs
3.3 Current Distribution inside a Rectangular Cylinder 87
and on y = 0
I/l ∂V ∂V ∂V I/l
Jn = − = −κ = κ ; = − ; |x − d| < s .
2s ∂n ∂y ∂y y=0 2κs
The boundary conditions for the normal derivative of the potential on x = 0 and x = b
are homogeneous. This is not the case for the planes y = 0 and y = b. As with
= 0 ; div J = 0, and J = κE
rot E the potential V is a solution of ΔV = 0, a practical
approach for the current problem is
∞
cos(nπx/a) " #
V (x, y) = A0 + B0 y + An cosh(nπy/a) + Bn cosh(nπ(y − b)/a) .
n=1
sinh(nπb/a)
The homogeneous boundary conditions ∂V /∂x|x=0,a = 0 are already satisfied and for
the normal derivative ∂V /∂y in y = b it follows
⎧
∞ ⎨ −I/(2κsl) ; |x − c| < s
∂V
= B + cos(nπx/a) (nπ/a) A =
∂y y=b
0 n
⎩ 0 ; |x − c| ≥ s.
n=1
Making use of the orthogonality of the trigonometric functions results after multipli-
cation with cos(mπx/a) and integration in
c+s
I/l I/l
B0 a = − 2s ; Am (mπ/a) a/2 = − cos(mπx/a)dx
2κs 2κs
c−s
b a b a
Pv /l = JEdx
dy = κ (grad V )2 dx dy .
y=0 x=0 y=0 x=0
88 3 Stationary Current Distributions
A straight forward evaluation of this surface integrals is possible, but with Green’s
second identity the surface integral is transformed into a simpler integral over the
contour C of the conducting cylinder.
c+s
d+s
∂V ∂V ∂V
Pv /l = κ V ds = κV (x = c , y = b) dx + κV (x = d, y = 0) dx
∂n ∂n ∂n
C c−s d−s
c+s d+s
∂V ∂V
Pv /l = V (x = c, y = b) κ dx +V (x = d, y = 0) −κ dx
∂y ∂y
c−s d−s
& '( ) & '( )
−I/l I/l
∞
cos(nπd/a)
− An cosh(nπb/a) + Bn .
sinh(nπb/a)
In the limit s → 0 the supply current is impressed on a line. With lim sin(x)/x = 1
x→0
the constants in the series for the potential and current density are
2I/l 2I/l
An |s→0 = − cos(nπc/a) ; Bn |s→0 = cos(nπd/a).
κnπ κnπ
However, the series for the power loss and the resistance are divergent, i.e. with s → 0
the resistance increases unrestricted.
This problem resembles problem 3.3, but the rectangular cylinder has been replaced
by a circular cylinder with conductivity κ and radius a. The input current supply lies
in the half plane ϕ = ϕ0 and the output current supply in ϕ = 0.
Calculate the current distribution and the resistance per axial length. In the limit case
c → 0 it is possible to calculate the current distribution with the method of images.
Calculate the potential by means of equivalent line sources.
3.4 Current Distribution inside a Circular Cylinder 89
k ® ¥
k = 0 2 c
I = l
%
'
' k ® ¥
0
I = l
k 2 c
The potential on the surface of the cylinder is constant in the region of the perfect
conducting electrodes. Furthermore because of 2c a we can assume that the normal
component of the current density is either constant at the end of the electrodes. Thus
the constant normal derivative of the potential is
∂V ∂V
− = − = En (ϕ)|=a = E (ϕ)|=a =
∂n =a ∂ =a
⎧
⎪
⎪
I c
− <ϕ<
c
⎪
⎪ ;
⎪
⎪ 2cκl a a
⎨
= I c c
⎪ − ; ϕ0 − < ϕ < ϕ0 +
⎪
⎪ 2cκl a a
⎪
⎪
⎪
⎩ 0 ; otherwise
2π
∂V
with adϕ = 0 .
∂n =a
0
On account of the symmetry plane ϕ = ϕ0 /2 the approach for the potential V ( , ϕ)
reads
∞ n
∞
2π 2π
n
− Cn sin n(ϕ − ϕ0 /2) sin k(ϕ − ϕ0 /2) dϕ = E (ϕ)|=a sin k (ϕ−ϕ0 /2) dϕ
n=1
a
0 0
& '( )
⎧
⎪
⎪
⎨ π ; n=k
=
⎪
⎪
⎩ 0 ; n = k
0+ a
c
c/a ϕ
n I
−I
c/a
I "
#
= sin n (ϕ − ϕ0 /2) − sin n (ϕ + ϕ0 /2) dϕ
2cκl
−c/a
c/a
−I 2I sin (n c/a)
= sin(n ϕ0 /2) cos (n ϕ) dϕ = − sin (n ϕ0 /2) .
cκl cκl n
−c/a
V ( , ϕ) = sin n (ϕ − ϕ0 /2) .
κπl n=1 n n c/a a
This series is divergent. The potentials for the points ( = a;ϕ = ϕ0 ) and ( = a; ϕ = 0)
and the resistance R/l increase unrestricted for c/a → 0.
This result can also be derived by means of the method of images, in analogy to the
method for line charges in front of a dielectric cylinder. For line currents ±I on the
surface of a conducting cylinder one gets image currents at the same positions, thus in
total line currents ±2I at positions ( = a; ϕ = ϕ0 ) and ( = a; ϕ = 0) in a space of
homogeneous conductivity κ. Now the resulting potential V ( , ϕ) is the real part of
the complex potential P (z).
I z − z1
P (z) = − ln ; z1 = a ; z2 = a exp(jϕ0 )
πκl z − z2
I
V ( , ϕ) = Re {P (z)} = − ln( 1 / 2 ) ; i = |z − zi | ; i = 1, 2
πκl
The complex electric field is
I
E(z) = −(dP /dz)∗ = [1/(z − z1 ) − 1/(z − z2 )]∗
πκl
and the lines of force are defined by the electric stream function
Ψs = Im {P (z)} = const .
Although the resulting solution from the method of images is of different form it is
identical with the preceding result. For the proof we look at a series expansion of the
logarithm of the distance |z − z0 | = 0 .
b
a %
k
I
following approach
V1 ( , z) = −E0 z + Am1 J0 (m1 ) exp(−m1 z) ; z≥0
m1
a 2
V2 ( , z) = −E0 z+ Bm2 J0 (m2 ) exp(m2 z) ; z ≤ 0.
b m2
The radial component of the current density vanishes on the boundary of the cylinder.
This leads for = a and = b with J0 = −J1 to
∂V1 ∂V2
= 0 ⇒ J (m a) = 0 ; = 0 ⇒ J1 (m2 b) = 0
∂ =a ∂ =b
1 1
x1r x1r
J1 (x1r ) = 0 ; r = 1, 2, 3, . . . ⇒ m1r = ; m2r =
a b
∞
V1 ( , z) = −E0 z + Ar J0 (x1r /a) exp(−x1r z/a) ; z≥0
r=1
a 2 ∞
V2 ( , z) = −E0 z+ Br J0 (x1r /b) exp(x1r z/b) ; z ≤ 0.
b r=1
a 2 b ∞
x1r
b
E0 J0 (x1s ) d − Br J0 (x1r ) J0 (x1s ) d =
b b r=1
b b b
0 0
& '( ) & '( )
I1 I2
a a ∞
x1r
= E0 J0 (x1s /b) d + Ar J0 (x1r /a) J0 (x1s /b) d
r=1
a
0 0
& '( ) & '( )
I3 I4
94 3 Stationary Current Distributions
z = a
-1
-2
-2 -1 0 1 x = a 2
with
1 2 2 b
I1 = 0 ; I2 = b J0 (x1s ) ; I3 = a J1 (x1s a/b)
2 x1s
x1s a/b J1 (x1s a/b) J0 (x1r )
I4 = − .
(x1r /a)2 − (x1s /b)2
Finally the boundary condition for the tangential component of the electric field, which
3.6 Current Distribution around a Conducting Sphere 95
Now the multiplication with J0 (x1s /a) and integration from = 0 to = a results
in
∞ a
1 2 2
As a J0 (x1s ) = Br J0 (x1r /b) J0 (x1s /a) d
2 r=1 0
& '( )
x1r a/b J1 (x1r a/b) J0 (x1s )
(x1r /b)2 − (x1s /a)2
and leads to a second system of equations.
∞
2 x1r a/b J1 (x1r a/b)
As = Br
J0 (x1s ) r=1 (x1r a/b)2 − x21s
The electric field and thus the current distribution is determined after solving both
systems of equations numerically.
The initial current density JE = −ez JE in the homogeneous space of conductivity κ2
is perturbed by a conducting sphere of radius a and conductivity κ1 .
Calculate the perturbed current distribution.
The initial current density in the homogeneous space of conductivity κ2 is
E.
JE = −ez JE = − [er cos ϑ − eϑ sin ϑ] JE = κ2 E
z
J E
# r
k 2
k a
1
The additional functions that describe the perturbation by the inserted conducting
sphere must be solutions of Laplace’s equation ΔV (r, ϑ) = 0. In spherical coordinates
these are the powers of the distance (rn , r−(n+1) ) and the Legendre’s polynomials
Pn (cos ϑ) and Qn (cos ϑ). Since the description of the initial field restricts to functions
of order n = 1, also the perturbed field is described by these functions and the approach
for the resulting potential reads as follows
⎧ ⎫
JE a ⎨ r/a + C1 (a/r)2 ⎬ r≥a
V (r, ϑ) = P1 (cos ϑ) ; .
κ2 ⎩ C2 r/a ⎭ r≤a
The unknown constants C1 and C2 result from the boundary conditions in the surface
r = a.
1 + C1 = C2 .
A second equation follows from the continuity of the normal component of the current
density
∂V ∂V
κ2 = κ1 ⇒ 1 − 2 C1 = κ1 /κ2 C2 .
∂r r>a ∂r r<a
r→a r→a
3.6 Current Distribution around a Conducting Sphere 97
ϑ
I(r, ϑ) = a = 2π
Jd er J r2 sin ϑ dϑ = const .
a ϑ =0
The result is
⎧ ⎫
ϑ ⎨ 3 ⎬
JE a 1 − 2 C1 (a/r)
I(r, ϑ) = −2πκ1,2 r2 cos ϑ sin ϑ dϑ
κ2 a ⎩ C2 ⎭
ϑ =0
⎧ ⎫
⎨ 1 − 2 C (a/r)3 ⎬ " #
1
I(r, ϑ)/I0 = cos 2ϑ − 1 = const
⎩ κ1 /κ2 C2 ⎭
z = a
-1
-2
-2 -1 0 1 x = a 2
In the plane z = 0 lies a circular conductor loop with radius a and with its center in
the point of origin. For the current supply radial line currents have been connected at
the points ( = a ; ϕ = ∓ϕ0 ), starting at infinity and carrying the constant current I.
I
%
'
' 0
e z
x
-' 0
a I
In general the magnetic field of a line current with contour C(rq ) is given by
rp ) = I
H( dsq × r/r3 ; r = rp − rq .
4π
C
Hence line currents with radial direction don’t contribute to the field at the point of
origin rp = 0, because the vectors r = rp − rq = −rq and dsq have the same direction
and therefore the cross product vanishes. The result of the integral over the contour
= 0
H
If the current I would split-up equally onto the wire elements, then the contributions
to the field are ez I/(4a)(ϕ0 /π) and −ez I/(4a)(1 − ϕ0 /π), and in total the magnetic
field in (z = 0; = 0) would be
= ez I/(4a)(2ϕ0 /π − 1).
H
Consider a circular ring with outer radius a, inner radius b > a, conductivity κ and a
very small thickness s a. The ring is open between −α < ϕ < α and a constant
current I is impressed at its ends, which are perfectly conducting.
Calculate the magnetic field in the center of the ring.
The current distribution and the excited magnetic field can be determined indepen-
dently in case of currents constant in time. As the thickness s is small and the electrodes
for the current supply are perfectly conducting, it is possible to describe the current
distribution as a current sheet K( ) = eϕ K( ). It is
2 2 2
ds = −
E grad V ds = − dV = V1 − V2 =
1 1 1
I
K( ) = .
ln(b/a)
4.2 Magnetic Field of a Current Sheet 101
r %
K (% ) q
® '
-®
a
b
I
k ;s
For the calculation of the field in the center rp = 0 we first calculate the contribution
of a differential current K( ) d on a circle of radius .
2π−α
= 1 −e π−α d
dH eϕ K( ) d × dϕ = ez K( )
4π 3 2π
ϕ=α
b b
= ez π − α
H K( )
d
= ez
π−α I d
2π 2π ln(b/a) 2
=a =a
= ez π − α
H
I
(1 − a/b) .
2πa ln(b/a)
In this problem the energy of the magnetic field and the inductance of conductors with
circular symmetry should be examined.
¹ %
¹ 0 K = e zK
' 0
With this result the energy in the magnetic field per unit length Wm /l would lead to
the integral
2π ∞ 2 c
1 = 1 a μ 2 d
Wm /l = Bda
H μK02 2
d dϕ = K0 2πa lim ,
2 2 2 c→∞
a ϕ=0 =a a
which is not convergent in the limit c → ∞. Obviously the field is non-regular and thus
the problem is not well-defined. (Because of the endless inductance it is not possible
to impress the current.)
4.3 Energy and Inductance of Conductors with Circular Symmetry 103
%
K = -e zK 0 a = b
' K = e zK 0
a
¹ 0 ¹ ¹ 0
= −ez K0 a/b at
Now the problem is modified by adding a second current sheet K
= b > a and the permeability in > b is μ0 .
The solution for the magnetic field is
⎧
⎪
⎪
⎨ 0 ; <a
= eϕ K0
H a/ ; a< <b
⎪
⎪
⎩
0 ; >b
and the energy of the magnetic field per unit length l becomes
2π a 2
1 a b 1
Wm /l = μ K02 d dϕ = μK02 a2 π ln = L/l I 2 .
2 a 2
ϕ=0 =0
The problem is regular and with the current I = 2aπ K0 the inductance per unit length
L/l is given by
μ b
L/l = ln .
2π a
This result is independent of μ0 , because there is no magnetic field in the surrounding
medium, i.e. the total current is zero.
= ez A( ):
It is possible to describe the magnetic field by means of a vector potential A
H = −eϕ 1 dA ;
= 1 rot A A( ) = −μK0
a
d = −μK0 a ln .
μ μ d c
c
104 4 Magnetic Field of Stationary Currents
In this case we can proof the calculation of the magnetic energy by the relation
1 1 1
Wm /l =
AJda = Kds
A = L/l I 2
2 2 2
a s
2π
1 μ b
= [A(a)K(a)a − A(b)K(b)b] dϕ = (2πaK0 )2 ln
2 4π a
0
and get the same result as before. Thus if the total current in planes z = const
vanishes, then the problem is regular and the relation for the inductance is valid.
Now we modify the problem a second time. The current sheet on the inner surface
= ez K0 cos ϕ. Again the
vanishes and the outer surface carries the current sheet K
total current vanishes:
2π
I = K( )bdϕ = 0 .
0
Vector potential and field follow from
= ez A( , ϕ) ;
A = b
ΔA( , ϕ) = 0 ;
= 1 rot A
H = e H ( , ϕ) + eϕ Hϕ ( , ϕ) = 1 e 1 ∂A − eϕ ∂A .
μ μ ∂ϕ ∂
Here the general solution function of Laplace’s equation for the vector potential is
∞ n n
a
A( , ϕ) = an + bn [cn cos(nϕ) + dn sin(nϕ)] .
n=1
a
where the constants an , bn , cn , and dn have been redefined. The solutions with n = 0
can be dropped as the total current vanishes.
For the determination of the constants firstly the continuity relation for the normal
component of the magnetic flux density, which is in this case identical to the continuity
of the vector potential, is evaluated with the result
cn = an + b n
n a n
b
dn = an + bn .
a b
4.3 Energy and Inductance of Conductors with Circular Symmetry 105
An additional equation results from the continuity of the tangential component of the
magnetic field in = a:
1 ∂A 1 ∂A
Hϕ | a<<b = Hϕ | <a ; =
→a →a μ0 ∂ <a μ ∂ a<<b
→a →a
μ0
⇒ cn = [an − bn ] = an + bn .
μ
Finally on = b applies
1 ∂A 1 ∂A
Hϕ | >b − Hϕ | a<<b = − + = K0 cos ϕ
→b →b μ0 ∂ >b μ ∂ a<<b
→b →b
∞
1 n 1n b
n−1 a n+1
dn + an − bn cos (n ϕ) = K0 cos ϕ
n=1
μ0 b μa a b
a
μ0 b
⇒ n = 1; d1 + a1 − b1 = μ0 K0 b .
μ a b
The solution of the system of equations is
μ0 b μ0 a
a1 1 + + b1 1 − = μ0 K0 b
μ a μ b
μ0 K0 b μ0 K0 a
a1 = ; a1 |μ→∞ =
(1 + μ0 /μ)b/a − k(1 − μ0 /μ)a/b 1 − (a/b)2
μ − μ0 −μ0 K0 a
b1 = −a1 k ; k = ; b1 |μ→∞ =
μ + μ0 1 − (a/b)2
c1 = a1 (1 − k) ; c1 |μ→∞ = 0; d1 = a1 [b/a − k a/b] ; d1 |μ→∞ = μ0 K0 b .
With this results the energy of the magnetic field calculates to
2π ∞ 2π
1 2 1 bdϕ =
Wm /l = μ H ( , ϕ) d dϕ = K(b) A(b)
2 2
ϕ=0 =0 ϕ=0
2π
1 b a π b a
= bK0 a1 −k cos2 ϕdϕ = b K 0 a1 −k .
2 a b 2 a b
0
The direct calculation of the integral of the square of the magnetic field over the cross
section is more complicated. In the limit μ → ∞ the energy is
π
Wm /l|μ→∞ = μ0 K02 b2 .
2
Although the total current vanishes in this problem, and thus it is a regular problem,
it is not possible to define a reasonable inductance, because the current sheet is a
continuous function of ϕ.
106 4 Magnetic Field of Stationary Currents
¹ 0 % ' 0
¹ I
'
¹ 0
a I
b -' 0
As the conductors of the excitation currents are thin wires, it is possible to apply
a scalar magnetic potential Vm . The line currents at positions ( = c; ϕ = ±ϕ0 )
can be equivalently replaced by a magnetic dipole layer m
m = e μ0 I on the surface
( = c; −ϕ0 < ϕ < ϕ0 ). The field equations are
= 0;
rot H = μH
B ; = −grad Vm ;
H ΔVm = 0 ; = c .
At first the exciting field H e = −grad Vme of the line currents and thus that of the
equivalent dipole layer in the homogeneous space of permeability μ0 is determined.
1 ∂ ∂Vme 1 ∂ 2 Vme
Vme = Vme ( , ϕ) ; ΔVme = + 2 = 0
∂ ∂ ∂ϕ2
Vme ( , ϕ) = R( )Φ(ϕ)
Due to the symmetry of the field with respect to ϕ = 0 and ϕ = π the approach for
the potential Vme is given by
⎧ ⎫ ⎧ ⎫
⎨ E ⎬ ∞ ⎨
En1 ( /c)n ⎬ <c
01
Vme ( , ϕ) = + cos(nϕ) ;
⎩ E02 + D02 ln / ln c ⎭ ⎩ En2 (c/ ) ⎭
n
> c.
n=1
The normal component e Be is continuous when passing the dipole layer at = c.
This results in
∂Vme ∂Vme
= ; En1 = −En2 = −En ; D02 = 0 .
∂ >c ∂ <c
→c →c
Now making use of the orthogonality relation for the trigonometric functions by mul-
tiplication with cos(kϕ) and integration from ϕ = 0 to ϕ = 2π leads to
ϕ0
2 I Iϕ0
2En π = I cos(nϕ)dϕ = I sin(nϕ0 ) ; En = sin(nϕ0 ) ; E0 = .
n nπ π
−ϕ0
The coefficient E0 = E02 −E01 in the Fourier series expansion of the potential difference
leads to a source field. Hence this part does not contribute to the magnetic field and
can be neglected.
Therewith the exciting potential is
⎧ ⎫
∞ ⎨ n ⎬
I −( /c) sin(nϕ0 )
Vme = cos(nϕ).
π n=1 ⎩ (c/ )n ⎭ n
A suitable approach for the resulting potential after the insertion of the permeable
hollow cylinder is
⎧⎧ ⎫ ⎫
⎪ ⎨ −( /c)n ⎬ ⎪
⎪
⎪ ⎪
⎪
⎪ + An ( /a) ⎪
n ⎪
⎪ <a
I ∞
sin(nϕ0 ) ⎨ ⎩ (c/ )n ⎭ ⎬
Vm = cos(nϕ) ;
π n=1 n ⎪
⎪ n n ⎪
⎪ a< <b
⎪ Bn ( /a) + Cn (a/ )
⎪ ⎪
⎪
⎪
⎩ ⎪
⎭
n > b.
Dn (b/ )
108 4 Magnetic Field of Stationary Currents
Another two equations follow from the continuity of the normal component of the
magnetic flux density:
∂Vm ∂Vm ∂Vm ∂Vm
μ0 = μ ; μ = μ
∂ <a ∂ a<<b ∂ a<<b ∂ >b
0
→a →a →b →b
and thus we get the method of images for a highly permeable cylinder.
4.4 Shielding of the Magnetic Field of a Parallel-Wire Line 109
ϕ) = ez A( , ϕ) ;
A( , ΔA( , ϕ) = 0 ; B = e 1 ∂A − eϕ ∂A
= rot A
∂ϕ ∂
e ( , ϕ) = ez Ae ( , ϕ) in the homogeneous
For the calculation of the exciting potential A
space of permittivity μ0 we replace the currents equivalently by a z-directed current
sheet K(ϕ) that vanishes except for ϕ = ±ϕ0 . With respect to the symmetry the
solutions of Laplace’s equation ΔAe ( , ϕ) = 0 are
⎧ ⎫
∞ ⎨ ( /c)n ⎬ ≤c
Ae ( , ϕ) = An sin(nϕ) ;
⎩ (c/ )n ⎭ ≥ c.
n=1
∞
2n/c An sin(nϕ) = μ0 K(ϕ) .
n=1
Making use of the orthogonality relation for trigonometric functions leads to
2π π
2μ0 I
2n/cAn π = μ0 K(ϕ) sin(nϕ)dϕ = 2μ0 /c sin(nϕ0 ) K(ϕ)cdϕ = sin(nϕ0 )
c
0 0
& '( )
I
μ0 I sin(nϕ0 )
⇒ An = .
π n
Now the ansatz for the resulting potential after inserting the permeable cylinder is
⎧⎧ ⎫ ⎫
⎪ ⎨ ( /c)n ⎬ ⎪
⎪
⎪ ⎪
⎪
⎪ + Bn ( /a)n ⎪⎪
⎪ <a
∞
μ0 I sin(nϕ0 ) ⎨ ⎩ (c/ ) n ⎭ ⎬
A( , ϕ) = sin(nϕ) ;
π n=1 n ⎪
⎪ Cn ( /a)n + Dn (a/ )n ⎪ ⎪ a< <b
⎪
⎪ ⎪
⎪
⎪
⎩ ⎪
⎭
n > b.
En (b/ )
The continuity of the normal component of the flux density and accordingly of the
vector potential on the surfaces = a and = b yields
and from the continuity of the tangential component of the magnetic field it follows
Cn = −2k(c/b)n(a/b)n /M ; Dn = 2(c/a)n /M
M = (1 + μ0 /μ)(1 − k 2 (a/b)2n ) ; k = (μ − μ0 )/(μ + μ0 ) .
1 .5
y = a
0 .5
-0 .5
-1
-1 .5
-1 .5 -1 -0 .5 0 0 .5 1 1 .5
x = a
Fig. 4.4–1: Magnetic lines of force for μ/μ0 = 10, b/a = 1.5, c/a = 0.5, and ϕ0 = π/3
⎧⎧ ⎫ ⎫
⎪
⎪ ⎨ ( /c)n ⎬ ⎪
⎪
∞ ⎪
⎨ n⎪
⎬
μ0 I sin(nϕ0 ) ⎩ + k(c/a)n
( /a)
n ⎭
A( , ϕ)|b→∞ = (c/ ) sin(nϕ)
π n=1 n ⎪
⎪ ⎪
⎪
⎪
⎩ 2 (c/ )n ⎪
⎭
1 + μ0 /μ
for < a and > a. This expression has the form
⎧
⎨ A ( , ϕ) + k A (a2 / , ϕ) ; < a
e e
A( , ϕ)|b→∞ =
⎩ (1 + k)Ae ( , ϕ) ; > a.
It is the analog expression for vector potentials of the method of images for a permeable
cylinder.
112 4 Magnetic Field of Stationary Currents
b
a %
k , ¹
I 0
Electric and magnetic fields can be determined independently if all fields are constant in
time. The current distribution of the present problem has been analyzed in chapter 3.
Outside the conducting material the ϕ-directed magnetic field is equal to the field of
a line current I0 on the z-axis. Inside the cylinder the magnetic field satisfies the
differential equation
2
= rot rot [eϕ H( , z)] = eϕ − ∂ H − ∂
rot rot H
1 ∂
( H) = 0
∂z 2 ∂ ∂
∂2H 1 ∂H H ∂2H
⇒ 2
+ − 2+ = 0; H( , z) = R( ) Z(z) .
∂ ∂ ∂z 2
The solutions of this differential equation are the ordinary Bessel functions Rm ( ) =
= {J1 (m ); N1 (m )} and the exponential functions Zm (z) = exp(±mz) with m = 0.
In case of m = 0 the solution functions are powers of .
4.5 Magnetic Field and Stationary Current Flow in a Cylinder 113
Making use of the orthogonality relation for Bessel functions results after multiplication
with J1 (x1s /b) and integration from = 0 to = b in
∞
a 1
S4 + Bs b2 J02 (x1s ) = S1 + Ar S2 + S3
b 2 r=1
with
a
ab
S1 = J1 (x1s /b) d = J2 (x1s a/b)
a x1s
0
a
x1r J1 (x1s a/b) J0 (x1r )
S2 = J1 (x1r /a) J1 (x1s /b) d =
(x1s /b)2 − (x1r /a)2
0
b
a ab
S3 = J1 (x1s /b) d = [J0 (x1s a/b) − J0 (x1s )]
x1s
a
b
b2
S4 = J1 (x1s /b) d = J2 (x1s ) .
b x1s
0
where the previously defined integral values have been used again.
Therewith all constants are defined by a system of equations and hence the magnetic
field is known.
The current density is calculated by
= −e ∂H + ez 1 ∂ ( H)
J = rot H
∂z ∂
and finally
∞
J1 ( , z) = H0 /a 2ez + x1r Ar ez J0 (x1r /a) + e J1 (x1r /a) exp(−x1r z/a)
r=1
∞
J2 ( , z) = H0 a/b2 2ez − x1r Br e J1 (x1r /b) − ez J0 (x1r /b) exp(x1r z/b) .
r=1
For z → ±∞ the current densities are homogeneous:
# 0
b
a r
#
¹ 0 ¹
It is practical to use the scalar magnetic potential Vm to calculate the magnetic field
H = −grad Vm . At first the exciting potential Vme of the conductor loop in the
homogeneous space of permeability μ0 is determined. For this purpose we consider the
magnetic dipole layer m m = r/rμ0 I = r/r mm on the spherical cap (r = b, ϑ < ϑ0 ),
that replaces the conductor loop equivalently.
The scalar magnetic potential
1 r
Vme = m
m da
4πμ0 r3
a
gives
∂Vme ∂Vme
μ0 = μ0
∂r r>b ∂r r<b
r→b r→b
n
⇒ −(n + 1) An = n Bn ; An = −
Bn .
n+1
And from the boundary condition for the potential at the dipole layer we get
⎧
∞ ⎨ I ; 0≤ϑ<ϑ
0
(An − Bn ) Pn (cos ϑ) = .
⎩ 0 ; ϑ0 < ϑ ≤ π
n=0
The unknown constants are derived by means of the orthogonality of the Legendre
polynomials Pn . With u = cos ϑ and u0 = cos ϑ0 we get
∞
+1 1
(An − Bn ) · Pn (u)Pk (u)du =I Pk (u)du
n=0 −1 u0
& '( )
$
2/(2k + 1) ; n = k
=
0 ; n = k
1
2 Pn+1 − Pn−1
(An − Bn ) = I ; Pn (1) = 1
2n + 1 2n + 1 u0
I 2n + 1 2n + 1
An − Bn = [Pn−1 (u0 ) − Pn+1 (u0 )] = − Bn = An
2 n+1 n
It should be noted, that the complementary magnetic dipole layer −r/rμ0 I on the
spherical segment (r = b; ϑ0 < ϑ ≤ π) would lead to the same result.
After insertion of the permeable sphere the approach for the resulting potential Vm
with ΔVm = 0 becomes
⎧ ⎫ ⎧ ⎫
⎨ V ⎬ ∞ ⎨
Cn (a/r)n+1 ⎬ r>a
me
Vm = + Pn (cos ϑ) ; .
⎩ 0 ⎭ ⎩ Dn (r/a)n ⎭ r<a
n=0
4.6 Force on a Conductor Loop in Front of a Permeable Sphere 117
For points r < a the exciting potential has the same dependencies on the coordinates
and ϑ as the added potential, thus it is omitted here.
The tangential component of the magnetic field and accordingly the potential is con-
tinuous in r = a:
I a n
−(n + 1) Kn + Cn = Dn .
2 b
A second condition is the continuity of the magnetic flux density in r = a:
a n−1 1
1 n
μ0 Bn n − (n + 1) Cn = μ Dn .
b b a a
(n + 1)(2n + 1)μ0 /μ I a n
Dn = − Kn .
n + (n + 1)μ0 /μ 2 b
In the limit μ → μ0 again only the exciting potential remains. The limit μ → ∞ leads
to
I a n
Cn |μ→∞ = (n + 1) Kn ; Dn |μ→∞ = 0
2 b
I
∞ a n a n+1
Vm |μ→∞ = Vme + (n + 1) Kn Pn (cos ϑ) ; r≥a
2 n=1 b r
a
This relation is the method of images for magnetic fields at highly permeable spheres.
Unlike the method of images for permeable half-spaces and permeable cylinders, in
case of the sphere we can write this law only for μ → ∞.
For ϑ0 = π/2, b a, and I/b = const the exciting magnetic field acting on the sphere
is almost homogeneous and z-directed
Thus in the approach for the resulting potential we just need terms with n = 1.
118 4 Magnetic Field of Stationary Currents
This results in
(n + 1) a n+2
∞
∂
Bsr (b, ϑ0 ) = −μ0 (Vm − Vme )| r=b = μ0 Cn Pn (cos ϑ0 )
∂r ϑ=ϑ0
n=0
a b
μ0 ∂ μ0
∞ a n+1
Bsϑ (b, ϑ0 ) = − (Vm − Vme )| r=b = Cn Pn (cos ϑ0 ) sin ϑ0
r ∂ϑ ϑ=ϑ0 b n=0 b
∞
a n+1 " #
F = −2π μ0 I sin2 ϑ0 Cn (n + 1) Pn (cos ϑ0 ) + cos ϑ0 Pn (cos ϑ0 ) .
n=0
b
Calculate the resulting magnetic field H and the screening factor ηs = |H e |/|H|
for
points inside the cylinder < a. In addition, give the limits μ = μ0 and μ → ∞.
4.7 Shielding of a Homogeneous Magnetic Field by a Hollow Cylinder 119
H e
'
a
x
¹ ¹
b 0
The homogeneous field might emerge from two homogeneous current sheets in planes
y = const with current flow in opposite z-direction. Thus it is possible to describe the
= ez A( , ϕ).
magnetic field by a z-directed vector potential A
e = ez Ae ( , ϕ) of the exciting field the magnetic flux density
From the vector potential A
Be = rot Ae = μ0 He = ex μ0 He is calculated to
∂Ae eϕ ∂Ae e ∂Ae ∂Ae
Be = grad Ae × ez = e + × ez = − eϕ
∂ ∂ϕ ∂ϕ ∂
= μ0 He = μ0 He ex = μ0 He [e cos ϕ − eϕ sin ϕ]
1 ∂Ae ∂Ae
= μ0 He cos ϕ ; = μ0 He sin ϕ .
∂ϕ ∂
The vector potential satisfies the differential equation ΔAe ( , ϕ) = 0 and the general
solution is
∞
Ae ( , ϕ) = (a0 + b0 ln ) (c0 + d0 ϕ) + (an n + bn −n ) (cn cos(nϕ) + dn sin(nϕ)) .
n=1
Because of the given exciting field only the solutions with n = 1 are needed and
furthermore b1 = 0, c1 = 0, and d1 = 1.
⇒ Ae ( , ϕ) = a1 sin ϕ = μ0 He sin ϕ
Accordingly the ansatz for the resulting vector potential A( , ϕ) in presence of the
cylinder is
⎧ ⎫
⎪
⎪ a1 /a ⎪
⎪ ≤a
⎨ ⎬
A( , ϕ) = μ0 He b a2 /a + b2 a/ sin(ϕ) ; a≤ ≤b
⎪
⎪ ⎪
⎪
⎩ ⎭
/b + b3 b/ ≥ b.
120 4 Magnetic Field of Stationary Currents
y = a
-1
-2
-2 -1 0 1 2 x = a
Fig. 4.7–1: Magnetic lines of force for μ/μ0 = 10 and b/a = 1.5
The evaluation of the boundary conditions leads to the unknown constants. The conti-
nuity of the normal component of the magnetic flux density is identical to the continuity
of the vector potential:
a 2 a
a1 = a2 + b 2 ; a2 + b 2 = (1 + b3 ) .
b b
Another two equations follow from the continuity of the tangential component Hϕ of
the magnetic field:
μ0 a 2 μ a
a1 = (a2 − b2 ) ; a2 − b 2 = (1 − b3 ) .
μ b μ0 b
Now there are four equations for four unknowns. With k = (μ − μ0 )/(μ + μ0 ) and
4.7 Shielding of a Homogeneous Magnetic Field by a Hollow Cylinder 121
´ s
4
1 0
a = b = { 0 .1 ; 2 /3 ; 0 .8 ; 0 .9 ; 0 .9 5 ; 0 .9 9 }
3
1 0
a = b
2
1 0
1
1 0
0
1 0
1 0 0
1 0 1 1 0 2
1 0 3
1 0 4 ¹ = ¹ 0
−2a/b
a1 |μ→∞ = 0 ; b2 |μ→∞ = ; a2 |μ→∞ = −b2 |μ→∞ ; b3 |μ→∞ = 1 .
1 − (a/b)2
The magnetic field vanishes inside the cylinder < b and on the cylinder = b the
field has only a radial component H .
122 4 Magnetic Field of Stationary Currents
With a finite permeability the magnetic field inside the cylinder < a is constant.
a
1 e ∂A ∂A 1 a1
H = − eϕ = He b e cos ϕ − eϕ sin ϕ
μ0 ∂ϕ ∂ a a
a1 b
= He b [e cos ϕ − eϕ sin ϕ] = ex He a1
a a
= ex He 2(1 − k)
H .
(1 + μ0 /μ)[1 − k 2 (a/b)2 ]
e |/|H|
Hence the screening factor ηs = |H is
a 2
1 (μ + μ0 )2
ηs = 1 − k2 .
4 μμ0 b
Consider two thin conductor loops of infinite extension in axial direction at positions
(r11 , r12 ) and (r21 , r22 ). The half-space y > 0, that contains the conductor loops, is of
permeability μ0 and the half-space y ≤ 0 is of permeability μ.
Calculate the mutual inductance of the two conductor loops per unit length.
r 2 2
r 1 1
y
r 2 1 r 1 2
¹ 0
A current I in the loop at positions (r11 ; r12 ) induces a magnetic flux Ψm per unit
length l through the second loop.
Ψm /l = B da/l = ds/l = A(r21 ) − A(r22 ) ;
A B = rot A
a C
4.8 Mutual Inductance of Plane Conductor Loops 123
y
m a
a a '
m a
a r 1 2
r r 2 2
2 1
¹ 0
a
r 1 1
x
2 a ¹
ª m = (lI )
L 0 = l
m
2
-2
-4
0 ¼ = 2 ¼ 3 ¼ = 2 ' 2 ¼
Fig. 4.8–1: Design (top) and normalized flux Ψm /(lI) (bottom) in dependence on the
angle ϕ
Parameters: μ/μ0 = 100, L0 /l = μ/(4π), m = 0.5; 0.7; 0.9
124 4 Magnetic Field of Stationary Currents
The vector potential A( r ) = ez A(r ) can be calculated from the currents ±I of the
∗
exciting loop at positions (r11 ; r12 ) and the image currents kI at r11 = r11 − 2ey (ey r11 )
∗
and −kI at r12 = r12 − 2ey (ey r12 ) in the homogeneous space of permeability μ0 .
∗
μ0 I |r − r11 | |r − r11 | μ − μ0
A = − ln + k ln ∗ | ; k =
2π |r − r12 | |r − r12 μ + μ0
The image currents describe the influence of the half-space z < 0 with permeability μ.
Ψm /l = A(r21 ) − A(r22 )
∗ ∗
μI |r21 − r11 | |r22 − r12 | |r21 − r11 | |r22 − r12 |
= − ln + k ln ∗ | | ∗ |
2π |r21 − r12 | |r22 − r11 | |r21 − r12 r22 − r11
The plane z = 0 contains a conductor loop consisting of two parallel wires with distance
2a and infinite extension in axial direction. A second loop with contour C is positioned
as shown in the figure. The conductor loops should be magnetically decoupled.
What is the condition for the geometry parameters? Find the parameter d if b = 2a
and c = 3a.
The conductor loops are decoupled if the magnetic flux, that is excited by the first loop
and runs through the second loop, vanishes.
We assume, that the parallel conductors carry the currents ±I and thus own the vector
potential
-a a b c d x
-h
k 2, ¹ 2
%
'
k 1, ¹ 1
i 0 a
e le c tric a l c o n n e c tio n
d2 A 1 dA
= ez A( ) ;
A = rot A
B ; ΔA = + = α2 A ; α2 = jωκμ
d 2 d
Inside the conducting material the solution of the differential equation requires modified
Bessel functions.
⎧
⎪ I (α )
⎨ C1 0 1 ; <a
I0 (α1 a)
A( ) =
⎪
⎩ C 2 I0 (α2 ) + C 3 K0 (α2 ) ; a < < b
I0 (α2 a) K0 (α2 a)
2j
α21 = jωκ1 μ1 ; α22 = jωκ2 μ2 =
δ2
= 1 rot [ez A] = −eϕ 1 ∂A = eϕ H( ) ;
H = −jω A
E ; J = κE
μ μ ∂
With the boundary conditions
1
|H( /b)|
|H 0 |
0.8
0.6 δ/b
0.4
0.2 κ1 /κ2 = 3
3 κ1 /κ2 = 3
|J( /b)|
|J 0 | δ/b = {0.25, 0.333, 0.5, 100}
2.5
δ/b
2
1.5
0.5
δ/b
0
0 0.2 0.4 0.6 0.8 1 /b 1.2
Fig. 5.1–1: Magnitude of the magnetic field (top) and current density (bottom) for
different normalized skin depths δ/b and b/a = 2
H 0 = i0 /(2πb), J 0 = i0 /(πb2 )
128 5 Quasi Stationary Fields – Eddy Currents
0.05
Im[i1 /i0 ] κ1 /κ2 = {0.5, 1, 2, 3, 5}
0
δ/b = 0.1 δ/b = 2
-0.05
-0.1
-0.15
-0.2
-0.25 κ1 /κ2
-0.3
-0.35
-0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
Re[i1 /i0 ]
0.35
Im[i2 /i0 ]
0.3
κ1 /κ2
0.25
0.2
0.15
0.1
0.05
δ/b = 2 δ/b = 0.1
0
κ1 /κ2 = {0.5, 1, 2, 3, 5}
-0.05
0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
Re[i2 /i0 ]
Fig. 5.1–2: Parametric plot of the normalized currents i1 /i0 (top) and i2 /i0 (bottom)
in dependence on the normalized skin depth δ/b for different conductivity
ratios κ1 /κ2
5.2 Rotating Conductor Loop 129
μ2 i0 I1 (α2 b) K1 (α2 b)
− = C2 − C3
2πα2 b I0 (α2 a) K0 (α2 a)
The constants follow from a straightforward solution of the equations.
Outside the conductor the field is given by
i0
H( ) = ; ≥ b.
2π
Finally the currents i1 and i2 in the conductors < a and a < < b are
2πa ∂A α1 I1 (α1 a)
i1 = 2πa H(a) = − = −2πa C 1 ; i2 = i0 − i1 .
μ1 ∂ =a μ1 I0 (α1 a)
A thin conductor loop with rectangular contour C and dimensions 2a and b rotates
ω (t) = ez ω0 t/T around its axis (x = c < a; y = 0). In x ≤ 0
with angular velocity
the loop is exposed to a homogeneous magnetic field H = ey H. The whole space is of
permeability μ.
What is the voltage u(t) measured by a highly resistive voltmeter in the time interval,
when the loop rotates by an angle of 180◦ from its initial position in the plane x = c
at the time t = 0?
The induced voltage results from a contour integration of the vector v × B, where
ω × r is the velocity of the conductor element at position r.
v =
u(t) = ds = μ (
v × B ds = μ [r (ωH)
ω × r) × H −ω (rH)]
ds
& '( )
C C C =0
0 0
= −μ ez ω0 t/T (r(ey H)) (ez dz) = −μ H ω0 t/T r ey dz = −μH ab ω0 t/T cos(ϕ(t))
−b −b
130 5 Quasi Stationary Fields – Eddy Currents
y
'
a
! (t)
z
c x
a
H = e yH
V
H = e yH
b ! (t) V
y x = c x
t
t2
ϕ(t) = ω(t ) dt = ω0
2T
0
1/2
ω0 t21 2T
ϕ(t1 ) = arcsin (c/a) = ; t1 = arcsin (c/a)
2 T ω0
1/2
ω0 t22 2T
ϕ(t2 ) = π − arcsin (c/a) = ; t2 = (π − arcsin (c/a))
2 T ω0
5.3 Force Caused by an Induced Current Distribution inside a Conducting Sphere 131
The result for t < 2π T /ω0 is
⎧
⎪
⎪
⎨ 0 ; 0 < t < t1
2
u(t) = −Ψm ω0 t/T cos(ω0 t /(2T )) ; t1 < t < t2
⎪
⎪
⎩
0 ; t2 < t < 2πT /ω0 .
A thin conductor loop that carries the current i = i0 cos(ωt + ϕ) = Re{i0 exp(jωt)} is
positioned at (r = b; ϑ = ϑ0 ) in front of a sphere with radius a < b and conductivity
κ. The permeability μ is constant.
Calculate the induced current distribution and determine the force acting on the con-
ductor loop. Use the result to analyze the limit cases, where the exciting field is
homogeneous and where the source is a magnetic dipole. Furthermore calculate the
limit of a high-frequency stimulation (ω → ∞).
i 0
b
# 0
# r
k , ¹
a
k = 0 , ¹
132 5 Quasi Stationary Fields – Eddy Currents
= rot A
; = eϕ A(r, ϑ) ; ∂A
B A = 0.
∂ϕ
In r < a it must satisfy the skin effect equation
∂
rot rot (eϕ A) = −κμ eϕ A ; rot rot (eϕ A) = −jωκμeϕA
∂t
and in r > a, r = b because of κ = 0 the equation
With spherical coordinates (r, ϑ, ϕ) the skin effect equation for A reads
∂2A ∂A ∂ 1 ∂
r2 + 2r + (sin ϑ A) = jωκμA · r2 = (αr)2 A ,
∂r2 ∂r ∂ϑ sin ϑ ∂ϑ
has been applied two times on the vector A = eϕ A. We try to find a solution of the
differential equation by means of the product:
A = R(r) · Θ(ϑ) .
This leads to
2
1 2 d R dR 1 d 1 d
r + 2r + (Θ sin ϑ) −(αr)2 = 0 .
R dr2 dr Θ dϑ sin ϑ dϑ
& '( )
−n(n + 1)
2 dΘ m2
1−u + n(n + 1) − Θ = 0.
du du 1 − u2
The solutions are the spherical harmonics Pnm (u) and Qm
n (u), that restrict to m = 1
in this case. Furthermore the functions of second kind Q1n (u) must be omitted as
they increase unbounded for points u = ±1, i.e. ϑ = 0, π, which are part of the
computational domain.
For the function R(r) the following differential equation remains.
d2 R dR " #
r2 2
+ 2r − n(n + 1) + (αr)2 R = 0 ; α2 = jωκμ = 2j/δ 2
dr dr
Its solution are the modified spherical Bessel functions
π π
Rn (αr) = In+ 21 (αr) and Tn (αr) = K 1 (αr) .
2αr 2αr n+ 2
Here only the functions Rn (αr) have to be regarded, because the functions Tn (αr) are
unbounded for r → 0.
Thus, all solution functions of the skin effect equation for a ϕ-directed vector potential
inside a conducting sphere are known.
In r > a outside the sphere α2 = 0 holds. The differential equation for the function
Θ(ϑ) applies unmodified and the differential equation for R becomes
d2 R dR
r2 + 2r − n(n + 1)R = 0 .
dr2 dr
Solutions of this equation are the functions rn and r−(n+1) .
Exciting Potential in the Homogeneous Space
A first step to the solution of the current problem is the determination of the exciting
e = eϕ Ae (r, ϑ) of the conductor loop at position (r = b; ϑ = ϑ0 ),
vector potential A
when the surrounding space is homogeneous and of permeability μ.
For this purpose we replace the loop by a current sheet K(ϑ) on the sphere r = b, that
is ϕ-directed and depends on the coordinate ϑ. In the surrounding medium outside of
the surface r = b the vector potential must satisfy
thus with the considerations of the previous section the solution is:
⎧ ⎫
∞ ⎨ B n+1 ⎬
r>b
n (b/r)
Ae = Pn1 (cos ϑ) ; .
⎩ ⎭
n=1 C n (r/b)n r<b
134 5 Quasi Stationary Fields – Eddy Currents
Bn = C n .
A second condition follows from the boundary condition for the tangential component
= 1/μB
of the magnetic field H = 1/μ rot A .
e e e
eϑ rot A
− A = μK(ϑ)
r<b
e e
r>b r→b
With the previously given expression for the curl-operator in spherical coordinates it
follows
1 ∂ "
#
r Ae |r<b − Ae |r>b r→b = μK(ϑ)
r ∂r
and with the expressions for Ae
n + 1 r n
∞ n+1
n b
Bn + Pn1 (cos ϑ) = μ K (ϑ) .
n=1
r r r b
r=b
the result is
π
(n + 1)! 2
(2n + 1)B n = μ Pn1 (cos ϑ0 ) sin ϑ0 K (ϑ)b d ϑ .
(n − 1)! 2n + 1
0
& '( )
i0
The current sheet K(ϑ) vanishes except for ϑ = ϑ0 and thus the integral on the right
yields the current i0 . Hence the constants are determined.
This approach already satisfies all boundary conditions in r = b. The evaluation of the
boundary conditions in r = a leads to the new constants Dn and E n .
From the continuity of the normal component of the magnetic flux density in r = a
follows the continuity of the vector potential.
n
(a/b) n+1
Dn Rn (αa) = + E n (b/a)
n(n + 1)
A second equation follows from the continuity of the tangential component of the
magnetic field
∂ ∂
(r A) − (r A)
= 0
∂r r<a ∂r r>a r→a
with the result
n+1
1 a n b
Dn [Rn (αa) + αa Rn (αa)] = Dn [(n + 1) Rn + αa Rn+1 ] = −n En .
n b a
Here it is
αr Rn (αr) = [n Rn (αr) + αr Rn+1 (αr)]
and the recurrence relation reads
Now, multiplying the first equation by n and adding the second equation yields
2n + 1 a n
Dn [(2n + 1) Rn (αa) + αa Rn+1 (αa)] =
& '( ) n(n + 1) b
αa Rn−1 (αa)
In this problem the current density J can be derived by the time derivative of the
vector potential.
J = −jωκA = eϕ J
2n + 1 Rn (αr) a n 1
∞
1 i0
J(r, ϑ) = − αa sin ϑ0 Pn (cos ϑ0 )Pn1 (cos ϑ)
2 a2 n=1
n(n + 1) Rn−1 (αa) b
In the limit ω → 0 it follows with the asymptotic expansion for spherical Bessel func-
tions
(αr)n
Rn (αr) |α→0 ≈
1 · 3 · 5 · · · (2n + 1)
r n
Rn (αr) (αr)n
1
= = .
αa Rn−1 (αa) α→0 (2n + 1)(αa)n α→0 2n + 1 a
μi0 sin ϑ0
∞
1 r n
A|ω→0 = Pn1 (cos ϑ0 ) Pn1 (cos ϑ) = Ae .
2 n=1
n(n + 1) b
As expected only the exciting potential remains. The differential equation for the
magnetic lines of force at the time t = 0 is
Because of the rotational symmetry the closed paths C are concentric circles around
the z-axis and we obtain the equation
With
1
Pn1 (cos ϑ0 ) = sin ϑ0 Pn (cos ϑ0 ) ; Pn (1) = n(n + 1)
2
it follows
M αa
∞ a n R (αr)
P 1 (cos ϑ) .
n
J(r, ϑ) = − (2n + 1)
4πμ (ab)2 n=1 b Rn−1 (αa) n
The flux density BS (t) is solely excited by the induced eddy currents, which follow
with
from the vector potential A S
) exp(jωt)}.
BS = Re{e rot (A S
F (t) = −2πb sin ϑ0 i(t) [BSr (t) sin ϑ + BSϑ (t) cos ϑ]|r=b;ϑ=ϑ0
1 ∂ cos ϑ ∂
= −2πb sin ϑ0 i(t) (sin ϑAS ) − (rAS )
r ∂ϑ r ∂r r=b;ϑ=ϑ0
∂AS ∂AS
= −2π sin ϑ0 i(t) sin ϑ − cos ϑb
∂ϑ ∂r r=b;ϑ=ϑ0
∞
μ sin ϑ0
AS (t) = A(t) − Ae (t) = Re {i0 E n exp(jωt)} (b/r)n+1 Pn1 (cos ϑ)Pn1 (cos ϑ0 )
2 n=1
u = cos ϑ ; u0 = cos ϑ0
1
∞
2
F = μπ(1 − u20 ) Re |i0 | E ∗n + i20 E n exp(2jωt) ·
2 n=1
1
2 dPn
· (1 − u0 ) 1
− (n + 1)u0 Pn (u0 ) Pn1 (u0 ) = F + F∼
du0
& '( )
1
−nPn+1 (u0 )
With the recurrence relation for the spherical harmonics the result reads
∞
1 2
F = − μ |i0 | π(1 − u20 ) nRe {E n } Pn+11
(u0 )Pn1 (u0 )
2 n=1
∞ 2
1 2 i0
F∼ = − μ |i0 | π(1 − u20 ) nRe E
2 n
exp(2jωt) 1
Pn+1 (u0 )Pn1 (u0 ) .
2 n=1
|i 0 |
The force is zero if the conductor loop is located in the plane ϑ0 = π/2, because
1
Pn+1 (0) Pn1 (0) = 0.
5.3 Force Caused by an Induced Current Distribution inside a Conducting Sphere 139
z = a
-1
-2 -1 0 1 x = a 2
z = a
-1
-2 -1 0 1 x = a 2
Fig. 5.3–1: Magnetic lines of force for b/a = 1.5, ϑ0 = π/6, with
δ/a = 1 (top) and δ/a = 0.2 (bottom)
Skin depth δ = 2/(ωκμ)
140 5 Quasi Stationary Fields – Eddy Currents
F
# = ¼ = 6
F 0
0
-1
b = a = { 1 .2 ; 1 .5 ; 2 }
1 0
b = a
-2
1 0
1 0 -3
0 0 .2 0 .4 0 .6 0 .8 1
± = a
Fig. 5.3–2: Average value of the force on the conductor loop as a function of the nor-
malized skin depth δ/a (F0 = μ0 |i0 |2 )
The inner conductor of radius a and the outer conductor of radii b > a and c > b of a
coaxial cable have the conductivity κ. The permeability μ is constant.
Calculate the frequency-dependent impedance per unit length Z/l = R/l + jωL/l and
analyze the limits of high and low frequencies.
The complex amplitude of the axially directed vector potential A( ) = ez A( ) satisfies
the differential equation
2j
ΔA = α2 A ; α2 = jωκμ = .
δ2
The solution functions are the modified Bessel functions I0 (α ) and K0 (α ). As the
magnetic field follows from the derivative in radial direction
If the coaxial cable carries the current i0 (positive z-directed in the inner conductor),
5.4 Impedance of a Coaxial Cable 141
k , ¹
k = 0 , ¹
%
k , ¹ '
a
i 0
i 0
b
Inside the conducting material the electric field is given by J = κE = ez κE( ).
Now applying the complex Poynting theorem with S = 1/2 E ∗ and integration
×H
over the conductor surfaces aL leads to
P ve + 2jω W m = − S daL = 1 |i |2 R + 1 jωμ |H|2 dvL .
2 0 2
aL vL
Here P ve is the time-averaged power loss and W m is the time-averaged energy of the
magnetic field in the conductor of volume vL . Hence it is
2 daL
R + jωLi = − 2 S
|i0 |
aL
2πa b
R/l + jωLi /l = E(a) H ∗ (a) − ∗
E(b) H (b)
|i0 |2 a
1 α I0 (αa) a I0 (αb) K1 (αc) + K0 (αb) I1 (αc)
= − .
2πa κ I1 (αa) b I1 (αb) K1 (αc) − K1 (αb) I1 (αc)
1
I0 (z) ≈ 1 ; I1 (z) ≈ z/2 ; K0 (z) ≈ − ln z ; K1 (z) ≈
z
hold an the resistance per unit length is
1 a2
R/l|ω→0 = 1+ = R0 /l .
κπa2 c2 − b 2
For high frequencies (|z| 1) one can make use of the asymptotic expansions of the
modified Bessel functions
1 π
In (z) ≈ exp(z) ; Kn (z) ≈ exp(−z)
2πz 2z
R = R 0
1
! L i= R 0
0
0 2 4 6 8 a = ± 1 0
Fig. 5.4–1: Real- and imaginary part of the impedance as a function of the normalized
reciprocal skin depth a/δ for b/a = 1.5
The frequency-dependent outer self-inductance La /l per unit length follows from the
time-averaged energy W ma /l of the field in the non-conducting area. The integral of
the energy density over the cross-section aq gives
1 ∗ daq =
1/4La/l|i0 |2 = W ma /l = Re HB
4
aq
b b
1 1 d
= πμ |H( )| d = πμ |H 0 |2 a2
2
2 2
a a
μ
⇒ La /l = ln(b/a) .
2π
Find the induced current distribution and calculate the time-averaged power loss in a
section 2a on the x-axis per unit length in z-direction. Analyze the limits ω → 0 and
ω → ∞.
y
x
¹ ,k
The constants C 1 and C 2 follow from the boundary conditions in y = 0. In this case
is equivalent to the
the continuity of the y-component of the magnetic flux density B
continuity of the vector potential:
exp(−πc/a) + C 1 = C 2 .
y = a
0 .5
-0 .5
-1
-1 0 x = a 1
y = a
0 .5
-0 .5
-1 0 x = a 1
Fig. 5.5–1: Magnetic lines of force for c/a = 0.2, μ/μ0 = 10, and δ/a = 1 (top) and
δ/a = 0.2 (bottom)
5.5 Induced Current Distribution in the Conducting Half-Space 147
k = (μ − μ0 )/(μ + μ0 ) .
Here we can recognize the method of images for vector potentials of stationary magnetic
fields at a permeable half-space, where the exciting vector potential, like in the present
case, has no y-component:
⎧
⎨ A e (x, y) + k A
e (x, −y) ; y ≥ 0
y)
A(x, = ; e = 0
ey A
ϕ=0
⎩
ω→0 (1 + k) Ae (x, y) ; y≤0
with A e (x, y) = A (x, y) .
e ϕ=0
At high frequencies it is
C 1 |ω→∞ = − exp(−πc/a)
and the vector potential becomes
y)
A(x, e (x, −y) ;
= Ae (x, y) − A y>0; e = 0.
ey A
ϕ=0
ω→∞
This is the method of images for a conducting half-space with high-frequency excitation
in y > 0, which is different to the method for stationary fields at a high-permeable half-
space:
y) ϕ=0 = A
A(x, e (x, y) + A
e (x, −y) ; y > 0 ; e = 0 .
ey A
ω→0
μ→∞
With a high-frequency excitation the field vanishes in y ≤ 0 and the plane y = 0 carries
a current sheet
0 0
K(x) = J(x, y)dy = −κjω A(x, y)dy
−∞ −∞
x
jωκμ0 aK 0 1
= − cos π C2
2π a (π/a)2 + α2 ω→∞
x 1 α2
μ0 aK 0
= − cos π C2
2π a μ (π/a) + α ω→∞
2 2
x 1 2πμ/μ
μ0 aK 0 0
= − cos π exp(−πc/a) α
2π a μ αa
ω→∞
1 ∂Ae
= −2 = 2 ey × H ez .
μ0 ∂y y=0 e
y=0
Hence the high-frequency current sheet can be calculated directly from the exciting
magnetic field.
148 5 Quasi Stationary Fields – Eddy Currents
The time-averaged power loss P v /l in a section 2a on the x-axis per unit length in
z-direction is
2a
1 ∗
×H
P v /l = − Re E ey dx
2 y=0
x=0
2a
1
= − Re {ez E(x, y = 0) × ex H ∗ (x, y = 0)} ey dx
2
x=0
⎧ 2a ⎫
1 ⎨ ⎬
= − Re E(x, 0) H ∗ (x, 0) dx
2 ⎩ ⎭
x=0
⎧ 2a ⎫
⎨ ∗ ⎬
1
1 ∂A(x, y < 0)
= − Re −jωA(x, y = 0) dx
2 ⎩ μ ∂y y→0 ⎭
x=0
⎧ 2a %
1 ⎨ μ20 a2 |K 0 |2 x |C |2
∗
2
= Re jω cos2 π (π/a)2 + α2 dx
2 ⎩ 4π 2 a μ
x=0
2
|K | 2
μ0
P v /l = − 02 |C 2 |2 Re (αa)2 π 2 + (αa)2 .
8π κ μ
= ds × 1/μ rot [ez A(x, y)] = 1/μ ds × (grad A × ez )
ds × H
⎡ ⎤
⎢ ⎥
= 1/μ ⎣grad A (ds ez ) −ez (ds grad A)⎦ = 0 ,
& '( ) & '( )
=0 dA
thus at the time t = 0 it is Re {A(x, y)} = const .
z = a
¹ 1, k 1
v K (x , t= 0 ) = e y K 0 c o s(¼ x = a )
0 x = a
¹ 2, k 2
∂A ∂A ∂ξ ∂A ∂2A ∂2A ∂A
= = −v ; 2
+ 2
+ κμv = 0
∂t ∂ξ ∂t ∂ξ ∂ξ ∂z ∂ξ
and for the given excitation the following approach holds
⎧⎧ ⎫ ⎫
⎨⎨ C exp(−γ z) ⎬ ⎬ z>0
1 1
A1,2 (ξ, z) = Re exp(j π ξ/a) ;
⎩⎩ C exp(γ2 z) ⎭ ⎭ z<0
2
with
0
2 2
− (π/a) + γ1,2 + jκ1,2 μ1,2 vπ/a = 0 ; γ1,2 = (π/a)2 − jκ1,2 μ1,2 vπ/a =
κ1,2 μ1,2 va
= π/a 1 − jλ1,2 ; .
λ1,2 =
π
The constants C1 and C2 follow from the evaluation of the boundary conditions in
z = 0. For the tangential component of the magnetic field it follows
" #
−Hξ (ξ, z < 0) + Hξ (ξ, z > 0) z→0 = K(ξ) = K0 cos(πξ/a) .
With
1 ∂A
Hξ = −
μ ∂z
one gets
7 8
Re (C 1 γ1 /μ1 + C 2 γ2 /μ2 ) exp(jπξ/a) = K0 cos(πξ/a) .
150 5 Quasi Stationary Fields – Eddy Currents
The normal component of the magnetic flux density is continuous in z = 0 and therefore
in this case also the vector potential. Hence it is C 1 = C 2 = C.
K0
C =
γ1 /μ1 + γ2 /μ2
7 8
A1,2 (x − vt, z) = A1,2 (x, z, t) = Re C exp(∓γ1,2 z) exp(jπ(x − vt)/a)
The fields are
∂
E1,2 (x, z, t) = − A1,2 (x, z, t) =
∂t
7 8
= Re j Cπ v/a exp(∓ γ1,2 z) exp(j π ξ/a) ; ξ = x − vt
and
1,2 (x, z, t) = 1 ∂A1,2 ∂A1,2
H ez − ex =
μ1,2 ∂x ∂z
1 7
8
= Re C ±ex γ1,2 + ez j π/a exp(∓γ1,2 z) exp(jπξ/a) .
μ1,2
In the limit of a high velocity the induced currents take the form of current sheets next
to the plane z = 0. For the lower half-space applies
0 0
K2 (ξ) = J2 (x, z, t)|v→∞ dz = κ E2 (x, z, t)|v→∞ dz
−∞ −∞
0 0
∂A2 ∂A2 ∂ξ
K2 (ξ) = −κ dz = −κ dz
∂t v→∞ ∂ξ ∂t v→∞
−∞ −∞
⎧ ⎫
⎨ κ vπ K0
0 ⎬
exp(γ2 z)dz
2
K2 (ξ) = Re j exp(jπξ/a)
⎩ a γ1 /μ1 + γ2 /μ2 ⎭
−∞
$ % v→∞
jvπ K0 exp(jπξ/a)
= Re =
a −j κ1 /κ2 μ2 /μ1 (vπ/a)2 − jvπ/a
v→∞
K0
K2 (ξ) = − cos(πξ/a) .
1+ κ1 /κ2 μ2 /μ1
The analog result for the upper half-space is
∞
∂A1 ∂ξ K0
K1 (ξ) = −κ dz = − cos(πξ/a) .
∂ξ ∂t v→∞ 1 + κ2 /κ1 μ1 /μ2
0
5.6 Induced Current Distribution by a Moving Conductor 151
0
z = a
-0 .2
-0 .4
-0 .6
-0 .8
-0 .8 -0 .4 0 0 .4 » = a 0 .8
0
z = a
-0 .2
-0 .4
-0 .6
-0 .8 -0 .4 0 0 .4 » = a 0 .8
Fig. 5.6–1: Magnetic lines of force for λ2 = 0.5 (top) and λ2 = 3 (bottom) with μ1 → ∞
Finally the equation for the magnetic lines of force at the time t = 0 is
7 8
A1,2 (ξ, z) = Re C exp(∓γ1,2 z + jπξ/a) = const
$ %
μ1 K 0 a exp (∓πz/a 1 − jλ1,2 ) exp (jπξ/a)
A1,2 (ξ, z) = Re √ √ = const
π 1 − jλ1 + μ1 /μ2 1 − jλ2
! H e
¹ 0
%
k , ¹ '
x
a
For the exciting rotating field in the homogeneous non-conducting space we can write
The solution of the problem requires to calculate the effect of only one component
ex H0 cos(ωt) = ex Re {H 0 exp(ωt)} of the exciting field. The effect of the second com-
ponent results from changing the phase and applying a transformation of coordinates.
The exiting vector potential of the subproblem is given by
= ez A ( , ϕ) ;
A = rot A
B = μ0 H
= ex μ0 H ,
e e e e e 0
with
= e 1 ∂Ae − eϕ ∂Ae ;
B = (e cos ϕ − eϕ sin ϕ)H
H
e e 0
∂ϕ ∂
and thus
Ae ( , ϕ) = μ0 H 0 sin ϕ .
It is a solution of the differential equation ΔAe ( , ϕ) = Δ(R( )Φ(ϕ)) = 0 with the gen-
eral solution functions Rn ( ) = { n , −n } and Φn (ϕ) = {cos(nϕ) ; sin(nϕ)}. Obviously
the solution for the exciting field is restricted to the terms with n = 1.
5.7 Conducting Cylinder Exposed to a Rotating Magnetic Field 153
∂ 2 A 1 ∂A 1 ∂2A
ΔA = 2
+ + 2 = α2 A ; α2 = jωκμ = 2j/δ 2 .
∂ ∂ ∂ϕ2
With the given excitation and because of the singularity of the modified Bessel function
Kn at = 0 the approach for the resulting vector potential reads with n = 1
⎧ ⎫
⎨ c I1 (α ) ⎪
⎪ ⎬ ≤a
A( , ϕ) = μ0 H 0 a I1 (αa) sin ϕ ; .
⎩ /a + b a/ ⎪
⎪ ⎭ ≥a
c = 1 + b.
A second equation results from the continuity of the tangential component of the
magnetic field e × H
1 ∂ 1 ∂
A ( < a, ϕ)|→a − A ( > a, ϕ)|→a = 0
μ ∂ μ0 ∂
that leads to
μ0 αaI1 (αa) 2I1 (αa)
c = 1 − b; c = .
μ I1 (αa) I1 (αa)(1 − μ0 /μ) + μ0 /μ αaI0 (αa)
= −H 0 (1 − b) sin ϕ|ω→∞ = −2H 0 sin ϕ = 2eϕ H e =a
ω→∞
K(ϕ) = −2H 0 sin ϕ = 2H eϕ =a
ω→∞
The time-averaged power loss in the cylinder per unit length is given by
⎧ 2π ⎫
1 ⎨ ∗
⎬
P v /l = − Re
E(a) ×H (a)e adϕ =
2 ⎩ ⎭
0
⎧ 2π ⎫ ⎧ 2π ⎫
1 ⎨ ⎬ 1 ⎨ ⎬
= Re E(a)H ∗ϕ (a)adϕ = Re jωA(a)H ∗0 (1 − b∗ ) sin ϕadϕ
2 ⎩ ⎭ 2 ⎩ ⎭
0 0
1 7 8
P v /l = a2 π μ0 |H 0 |2 Re jω(2c − |c|2 ) = −μ0 πa2 |H 0 |2 ω Im {c} .
2
The calculation of the total current distribution induced by the rotating field requires
the superposition of the field described by the second part of the initial z-directed
vector potential Ae = −μ0 H 0 cos ϕ .
This is done by the substitution of sin ϕ with sin ϕ − j(− cos ϕ) = j exp(−jϕ)
in the expression for the current density.
= ez A of the rotating field applies in analogy
For the vector potential AD D
⎧ ⎫
⎨ c I1 (α ) ⎬
AD ( , ϕ) = jμ0 H 0 a I1 (αa) exp(−jϕ)
⎩ ⎭
/a + b a/
and the field follows from the given relations.
= 0 at the time t = 0 is Re {A ( , ϕ)} =
The equation for magnetic lines of force ds × H D
const . Examples of magnetic lines of force are shown in the figures.
5.7 Conducting Cylinder Exposed to a Rotating Magnetic Field 155
y = a
-1
-1 0 1 x = a
y = a
-1
-1 0 1 x = a
Fig. 5.7–1: Magnetic lines of force in a conducting permeable cylinder, that is exposed
to a rotating magnetic field at time t = 0
Parameters: μ/μ0 = 1, a/δ = 3 (top), a/δ = 10 (bottom)
156 5 Quasi Stationary Fields – Eddy Currents
y = a
-1
-1 0 1 x = a
y = a
-1
-1 0 1 x = a
Fig. 5.7–2: Magnetic lines of force in a conducting permeable cylinder, that is exposed
to a rotating magnetic field at time t = 0
Parameters: μ/μ0 = 103 , a/δ = 5 (top), a/δ = 10 (bottom)
5.8 Power Loss and Energy Balance inside a Conducting Sphere 157
Find the induced current distribution in the sphere, when in contrast to problem 5.3
the exciting current i(t) has an arbitrary time-dependence. Additionally calculate the
time-averaged power loss in case of a time-harmonic excitation and analyze the energy
balance when a direct current is turned off.
i(t)
r E
# E
# r
k , ¹
¹ 0 a
In this case the rotationally symmetric field is described by a second order vector
potential
= r W (r, ϑ, t) ; ∂W
W = 0.
∂ϕ
With it, the field is given by
∂A
= rot W
A = −r × grad W ; = rot A
B = rot rot W
; J = −κ .
= κE
∂t
158 5 Quasi Stationary Fields – Eddy Currents
Time-Harmonic Field
W (r, ϑ) = r W ; ΔW E = 0 ; rE
r=
E E
⎧ ⎫
∞ ⎨
E n (r/rE )n ⎬ r < rE
W E (r, ϑ) = Pn (cos ϑ) ;
⎩ D (rE /r)n+1 ⎭ r > rE
n=0 n
= −r × grad W ; = rot rot W
= grad ∂
AE E B E E (rW E ) = μ0 H E
∂r
For the determination of the constants Dn and E n we evaluate the boundary conditions
on the sphere r = rE . The continuity of the normal component of the magnetic flux
density requires
∂2 ∂2
(rW E ) = (rW E ) ; Dn = E n .
∂r2 r<rE ∂r2 r>rE
r→rE r→rE
With
1 1 ∂2
H Eϑ = (rW E )
μ0 r ∂r∂ϑ
and u = cos ϑ it follows
⎧ ⎫
1 1
∞ ⎨ (n + 1)(r/r )n ⎬ r < rE
E
H Eϑ = − E 1 − u2 Pn (u) ;
μ0 r n=0 n ⎩ −n(rE /r)n+1 ⎭ r > rE
and thus it is
1
∞
(2n + 1)E n 1 − u2 Pn (u) = rE K(u) .
μ0 n=0
√
Now the multiplication with 1 − u2 Pk (u) and integration from u = −1 to u = +1
5.8 Power Loss and Energy Balance inside a Conducting Sphere 159
leads to
+1
μ0
En = K(u)rE 1 − u2 Pn (u)du
2n(n + 1)
u=−1
π
μ0
= P (uE )(1 − u2E ) K(ϑ)rE dϑ ; uE = cos ϑE
2n(n + 1) n
ϑ=0
& '( )
i0
1 P (uE ) 1 Pn−1 (uE ) − Pn+1 (uE )
En = μ0 i0 (1 − u2E ) n = μ0 i0
2 n(n + 1) 2 2n + 1
In the presence of the conducting sphere the field is described by the second order
potential
W = r W ; ΔW = α2 W
⎧ ⎧ ⎫⎫
⎪
⎪ ⎨ (r/r )n ⎬⎪ ⎪
∞ ⎨⎪ E ⎪
⎬
An (a/r)n+1 + E n r>a
W = ⎩ (rE /r)n+1 ⎭ Pn (cos ϑ) ;
⎪ ⎪
n=0 ⎪
⎪ ⎪
⎪
⎩ B Rn (αr) ⎭ r<a
n
= grad ∂ (rW ) − r α2 W
B
∂r
160 5 Quasi Stationary Fields – Eddy Currents
and because of the continuity of the normal component er B in r = a and thus the
continuity of W it is
An + E n (a/rE )n = B n Rn (αa) .
= B/μ
A second equation results from the continuity of the tangential component of H
μ0
B [(n + 1) Rn (αa) + αa Rn+1 (αa)] = −nAn + (n + 1)E n (a/rE )n .
μ n
This system of equations has the solution
n
2n + 1 a
Bn = En
nRn (αa) + μ0 /μ [(n + 1) Rn (αa) + αa Rn+1 (αa)] rE
n
(2n + 1)Rn (αa) a
An = −1 En .
nRn (αa) + μ0 /μ [(n + 1) Rn (αa) + αa Rn+1 (αa)] rE
1
∞
Hϑ = − B n [(n + 1)Rn (αr) + αrRn+1 (αr)] 1 − u2 Pn (u)
μr n=0
Transient Field
The basis for the calculation of the field in case of an arbitrary time-dependence of
the exciting current is the field, that emerges when a constant current is turned off.
Therefore we need at first knowledge of the stationary field, that is excited by a constant
current I0 in the conductor loop. The stationary field follows from the corresponding
field equations, or alternatively from the limit ω → 0 in the time-periodic case.
n r n
2n + 1 r
B n Rn (αr)|ω→0 = En = Bn(0)
n + μ0 /μ(n + 1) rE a
n
2n + 1 a
Bn(0) = En
n + μ/μ0 (n + 1) rE
n
1 − u2E Pn (uE ) 2n + 1 a
En = μ0 I0 ; An |ω→0 = −1 En
2 n(n + 1) n + μ0 /μ(n + 1) rE
n
(n + 1)(1 − μ0 /μ) a
An |ω→0 = En = A(0)
n
n + (n + 1)μ0 /μ rE
After turning off the current I0 at the time t = 0 the second order potential
A = r WA (r, ϑ, t) is given by
W
⎧ ⎫
∞ ∞ ⎨ ⎬
ans jn (λns r) κμ
WA = Pn (cos ϑ) exp(−t/τns ) ; τns = 2 .
⎩ bns jn (λns a)(a/r)n+1 ⎭ λ ns
n=0 s=1
A = rot (rWA ) ;
A A = rot rot (rWA ) = grad ∂ (rWA ) − κμ ∂ WA .
B
∂r ∂t
162 5 Quasi Stationary Fields – Eddy Currents
For the determination of the constants ans and bns we evaluate the boundary conditions
in r = a and furthermore at the time t = 0 the field matches the initial stationary
solution. The continuity of the normal component of the magnetic flux density requires
2
∂ 2 ∂WA
∂2
r 2 (rWA ) − κμr = r 2 (rWA ) .
∂r ∂t r<a ∂r r>a
r→a r→a
With the recurrence relations and the derivatives for spherical Bessel functions this
results in
ans = bns .
Now with this result the second continuity relation for the tangential component of the
magnetic field results in the eigenvalue equation
and thus leads to the eigenvalues λns = xns /a. Hence the potential is
⎧ ⎫
∞
∞ ⎨ j (x r/a) ⎬ κμa2
n ns
WA = ans Pn (u) exp(−t/τns ) ; τns = 2 .
⎩ jn (xns )(a/r)n+1 ⎭ xns
n=0 s=1
Finally the constant ans follows from the initial stationary field at the time t = 0.
∞
∞ ∞
ans jn (xns r/a) Pn (u) = Bn(0) (r/a)n Pn (u)
n=0 s=1 n=0
∞
Bn(0) (r/a)n = ans jn (xns r/a) ; u = cos ϑ
s=1
Due to the orthogonality relation for the the spherical Bessel functions the multiplica-
tion with jn (xnp r/a)w(r) and w(r) = 2λnp r2 /π and integration results in
2jn+1 (xns )
ans = B (0) .
xns [jn2 (xns ) − jn−1 (xns )jn+1 (xns )] n
With it, the transient field for t > 0 is determined. In case of an arbitrary time-
dependence of the current i(t) for t > 0 and i(t) = 0 for t < 0 the potential WA (r, t) is
given, if the exponential function exp(−t/τns ) in WA is replaced by the function
⎡ ⎤
t
1 ⎣ di(t − τ )
fns (t) = i(0) [1 − exp(−t/τns )] + [1 − exp(−τ /τns )] dτ ⎦ .
I0 dt
0
5.8 Power Loss and Energy Balance inside a Conducting Sphere 163
Energy Balance
When the initially constant current is turned off the energy balance requires
Here Wm0 is the energy of the magnetic field inside the sphere before the turn-off,
WmA is the energy transferred trough the surface r = a after the turn-off, and Wms is
the dissipated energy inside the sphere for t > 0. It follows
∞
1 2 A (r, t) daK dt = Wm0 − WmA =
A (r, t) × H
μ HA (r, t < 0)dv − E
2
vK t=0 aK
∞
= A (r, t) JA (r, t)dvdt = Wms .
E
t=0 vK
With
1 ∂ ∂WA ∂WA
HA × EA = grad (rWA ) − κμr × r × grad = −S
μ ∂r ∂t ∂t
r r ∂ ∂WA eϑ ∂
S = − grad s (rWA ) grad s ; grad s =
r μ ∂r ∂t r ∂ϑ
1 ∂ ∂ 1 ∂ ∂WA
= − (rWA )
μ ∂ϑ ∂r r ∂ϑ ∂t
164 5 Quasi Stationary Fields – Eddy Currents
z = a
-1
-2 -1 0 1 x = a 2
Fig. 5.8–1: Magnetic lines of force before the turn-off with μ = 100μ0 ,
κ = 2 · 106 [S/m], rE /a = 2, and ϑE = π/6
it follows:
∞ 2π +1
∞
∞
1
WmA = − ans [xns jn−1 (xns ) − n jn (xns )] ·
μ n=0 s=1
t=0 0 −1
∞ ∞
1
2
· − 1 − u Pn (u) exp(−t/τns ) akr jk (xkr ) − 1 − u2 Pk (u) ·
a
k=0 r=1
−x2kr
· exp(−t/τkr ) · a2 (−du) dϕdt
κμa2
∞
∞
∞ ∞
2πa
= − ans (xns jn−1 (xns ) − njn (xns )) anr jn (xnr ) ·
μ
t=0 n=0 s=1 r=1
+1
−x2nr
· (1 − u2 )Pn2 (u)du exp(−t/τns ) exp(−t/τnr ) dt
κμa2
−1
5.8 Power Loss and Energy Balance inside a Conducting Sphere 165
z = a
-1
-2 -1 0 1 x = a 2
1
z = a
-1
-1 0 1 x = a
Fig. 5.8–2: Magnetic lines of force at the time t/τ01 = 10−3 (top) and t/τ01 = 10−2
(bottom) after the turn-off
Parameters see Fig. 5.8–1.
∞ ∞ ∞
2πa 2n(n + 1)
WmA = ans anr jn (xnr )·
μ n=0 2n + 1 s=1 r=1
∞
x2
· [xns jn−1 (xns ) − njn (xns )] · nr2 exp(−t/τns ) exp(−t/τnr ) dt
κμa
t=0
166 5 Quasi Stationary Fields – Eddy Currents
∞ ∞ ∞
4πa n(n + 1) [xns jn−1 (xns ) − njn (xns )]
WmA = ans anr jn (xnr ) x2nr .
μ n=0 2n + 1 s=1 r=1 x2ns + x2nr
∞ 2
∂ 2 WA
Wms = κ dvdt
∂ϑ∂t
t=0 vK
∞
+1 ∞ a
∞
2 r exp(−2t/τns )
= 2πκ (1 − u )Pn2 (u) du a2ns jn2 (xns )r2 dr 2
dt
n=0u=−1 s=1r=0
a τns
t=0
∞ ∞
πa n(n + 1) 2 2 " 2 #
Wms = ans xns jn (xns ) − jn−1 (xns )jn+1 (xns ) .
μ n=0 2n + 1 s=1
The proof of the validity of the energy balance is done by a numerical calculation of
the expressions above.
k = 0 , ¹
%
k , ¹ ' K = e z K 0
a
K = -e z K 0
Solution functions are the modified Bessel functions In (α ) and Kn (α ) in < a, the
powers of the distance to the axis n and −n in > a, and the trigonometric functions
sin(nϕ) and cos(nϕ).
At first the exciting vector potential A ( , ϕ) = ez A ( , ϕ) of the current sheet at
e e
= b in the homogeneous space of permeability μ will be determined.
⎧ ⎫
∞ ⎨ ( /b)n ⎬ ≤b
ΔAe = 0 ; = b ; Ae ( , ϕ) = an sin(nϕ) ;
⎩ (b/ )n ⎭ ≥b
n=1
The continuity of the normal component of the magnetic flux density B is identical to
the continuity of the vector potential itself and is already satisfied. Here the solution
with n = 0 is omitted because of
2π
K(ϕ) dϕ = 0 .
0
= rot A
The magnetic flux density B is given by
e e
" #
= rot [ez A ] = grad A × ez = e 1 ∂Ae − eϕ ∂Ae = μ e H + eϕ H
B e e e e eϕ .
∂ϕ ∂
For the determination of an the boundary condition for the tangential component of
the magnetic field in = b is evaluated.
⎧
⎨ K
0 ; 0 < ϕ< π
H eϕ >b − H eϕ <b = K(ϕ) =
→b →b ⎩ −K ; π < ϕ < 2π 0
168 5 Quasi Stationary Fields – Eddy Currents
∂Ae ∂Ae
− = μK(ϕ)
∂ <b ∂ >b
→b →b
⎧
∞ ⎨ K
0 ; 0 < ϕ< π
2n an /b sin(nϕ) = μ
⎩ −K
n=1 0 ; π < ϕ < 2π
After inserting the conducting cylinder the ansatz for the resulting vector potential
A( , ϕ) is
⎧ ⎫
⎪
⎪ bn In (α )/In (αa) ⎪
⎪ ≤a
∞ ⎨⎪ ⎧ ⎫ ⎪
⎬
2μK 0 b sin(nϕ)
A( , ϕ) = ⎨ ( /b)n ⎬ ; a≤ ≤b
π ⎪ ⎪ n2
k=0 ⎪ ⎪
n
⎪
⎩ ⎩ (b/ )n ⎭ + dn (a/ ) ⎪
⎭
≥b.
in = a
Again the continuity of the normal component of the magnetic flux density B
is identical to the continuity of the vector potential.
bn = (a/b)n + dn
The time-averaged power loss per unit length l in the cylinder follows from the inte-
gration of the Poynting vector over the surface of the cylinder.
2π 2π ∗
P v /l = −
1 ∗
×H
Re E e ds =
1
Re
jωA ∂A
adϕ
2 =a 2 μ ∂ =a
ϕ=0 ϕ=0
2
∞ ∗
ωa 2μb|K 0 | In (αa) 1
= Re j|bn |2 α π
2μ π In (αa) n4
k=0
5.9 Induced Current Distribution in a Conducting Cylinder 169
3
y = a
2
-1
-2
-3
-3 -2 -1 0 1 2 3
x = a
3
y = a
2
-1
-2
-3
-3 -2 -1 0 1 2 3
x = a
∞ ∗
2 I (αa)
P v /l = − ωμ (|K 0 | b)2 |bn |2 /n4 Im αa n ; n = 2k + 1
π In (αa)
k=0
= 0 at the time t = 0 reads
Finally the equation for the magnetic lines of force ds × H
k = 0 , ¹
b
a %
k , ¹
i 0
∂H
+ κμ
rot rot H = 0;
rot rot H( , z) + α2 H( ,
z) = 0
∂t
2
∂H 1 ∂ ∂ H ∂ 1 ∂
rot −e + ez ( H) = eϕ − 2 − ( H) = −α2eϕ H
∂z ∂ ∂z ∂ ∂
∂2H 1 ∂H 1 ∂2H
⇒ 2
+ − 2H + = α2 H .
∂ ∂ ∂z 2
The ansatz
H( , z) = R( ) Z(z)
leads to the equation
1 d2 R 1 dR R 1 d2 Z
2
+ − 2
+ = α2 = jωκμ = 2j/δ 2
R d d Z dz 2
that separates to
d2 R 1 dR R 2
2
+ − 2 + m2 R = 0 ; δ =
d d ωκμ
and
d2 Z m
− (α2 + m2 )Z m = 0.
dz 2
Solutions are the Bessel functions of order one
Rm ( ) = C 1 J1 (m ) + C 2 N1 (m )
with H 0 = i0 /(2πa). The preceding terms with modified Bessel functions describe the
homogeneous field far away from the discontinuity in z = 0.
J1 (m1 a) = 0 ⇒ m1 = x1r /a ; r = 1, 2, 3, . . .
J1 (m2 b) = 0 ⇒ m2 = x1r /b ; r = 1, 2, 3, . . . .
⎧ ∞
∞ ⎪
⎨ I1 (α )
a I1 (α ) a + C J1 (x1r /a) ; ≤ a
+ C 2r J1 (x1r /b) = I1 (αa) r=1 1r
b I1 (αb) b r=1 ⎪
⎩
a/ ; a≤ ≤b
∞
I1 (α ) x1r J1 (x1r /a)
= −2 2 + (αa)2
I1 (αa) x
r=1 1r
J0 (x1r )
5.10 Cylinder with Stepped Down Diameter 173
it is
∞
a 2x1r 1
C 2r − 2 J1 (x1r /b) =
b r=1 x1r + (αb)2 J0 (x1r )
⎧ ∞
⎪
⎪ 2x1r 1
⎨ C 1r − 2 J1 (x1r /a) ; < a
= x1r + (αa)2 J0 (x1r )
⎪
⎪
r=1
⎩ a/ ; a < < b.
2 ∞
a 2x1s b 2 2x1r
C 2s − 2 J (x1s ) = C 1r − 2 ·
b (x1s + (αb)2 )J0 2 0 r=1
(x1r + (αa)2 )J0 (x1r )
a b
· J1 (x1r /a)J1 (x1s /b) d + a/ J1 (x1s /b) d .
0 a
& '( ) & '( )
x1r J1 (x1s a/b)J0 (x1r ) ab
[J0 (x1s a/b) − J0 (x1s )]
(x1s /b)2 − (x1r /a)2 x1s
A second system results from the continuity of the tangential component J of the
current density
∂H ∂H ∂H
J = − ; = ; <a
∂z ∂z z<0 ∂z z>0
z→0 z→0
∞
∞
a
− C 1r J1 (x1r ) α2 + (x1r /a)2 = C 2r J1 (x1r ) α2 + (x1r /b)2 .
r=1
a b r=1 b
174 5 Quasi Stationary Fields – Eddy Currents
z = a
-1
-2
-2 -1 0 1 x = a 2
= a/b 2
C 2r α + (x1r /b)2 J1 (x1r /b)J1 (x1s /a) d .
r=1 0
& '( )
x1s J1 (x1r a/b)J0 (x1s )
(x1r /b)2 − (x1s /a)2
The combination of both systems leads to a linear system of equations for the deter-
mination of the unknown constants. Exemplary results for the streamlines at different
times are presented in the figures.
5.10 Cylinder with Stepped Down Diameter 175
z = a
-1
-2
-2 -1 0 1 x = a 2
2
z = a
-1
-2
-2 -1 0 1 x = a 2
Fig. 5.10–2: Streamlines at the time ωt = π/4 (top) and ωt = π/2 (bottom) with
b/a = 2 and b/δ = 3
176 5 Quasi Stationary Fields – Eddy Currents
Find the magnetic field inside the groove and the current distribution in the conductors.
Furthermore calculate the ratio η = i2 /i1 of the currents in conductor (1) and (2) in
dependence on the frequency ω.
y
e le c tric a l c o n n e c tio n
h
i 0
k 1, ¹ k 2, ¹ k = 0 , ¹
0 a b x
¹ ® ¥
H(x) = ey H(x)
d2 H
= −jωκμH
rot rot H ; = α21,2 H
dx2
with
1+j
α21 = jωκ1 μ ; α22 = jωκ2 μ = 2j/δ 2 ⇒ α2 = .
δ
The current density is given by
dH(x)
J(x) = rot H(x) = ez J(x) = ez
dx
5.11 Frequency-Dependent Current Distribution in Conductors 177
C2 + C3 = 1 .
;
J = rot H 1/κ1 dH/dx| x<a = 1/κ2 dH/dx| x>a
x→a x→a
i1 = H(a) · h ; i1 /i0 = C 1
i2 = H(b) · h − i1 ; i2 /i0 = C 2 + C 3 − C 1
1 1
i1 /i0 = ; i2 /i0 = 1 − .
cosh(α2 (b − a))(1 + N ) cosh(α2 (b − a))(1 + N )
Hence the current ratio η is
i2
η = = (N + 1) cosh(α2 (b − a)) − 1 .
i1
178 5 Quasi Stationary Fields – Eddy Currents
|H ( x = b ) | 1
|H 0 | k 1/k 2 = 3
0 .8 ± = b = { 0 .2 5 ; 0 .3 3 3 ; 0 .5 ; 1 0 0 }
0 .6
0 .4
± = b
0 .2
0
0 0 .2 0 .4 0 .6 0 .8 x = b 1
|J ( x = b ) |
|J 0 | 5 k 1/k 2 = 3
± = b = { 0 .2 5 ; 0 .3 3 3 ; 0 .5 ; 1 0 0 }
4
2 ± = b
± = b
1
0
0 0 .2 0 .4 0 .6 0 .8 x = b 1
Fig. 5.11–1: Magnitude of the magnetic field (top) and current density (bottom) for
different normalized skin depths δ/b and b/a = 2 (J 0 = H 0 /b)
5.11 Frequency-Dependent Current Distribution in Conductors 179
0 ± = b = 2 k 1/k 2 = { 1 / 2 ; 1 ; 2 ; 3 ; 5 }
I m [i1= i0]
± = b = 0 .1
-0 .1 k 1/k 2
-0 .2
-0 .3
-0 .4
0 0 .2 0 .4 0 .6 0 .8
R e [i1= i0]
0 .4
I m [i2= i0]
0 .3
0 .2
k 1/k 2
0 .1
± = b = 2
k 1/k = { 1 / 2 ; 1 ; 2 ; 3 ; 5 } ± = b = 0 .1
0 2
0 .2 0 .4 0 .6 0 .8 1
R e [i2= i0]
Fig. 5.11–2: Parametric plot of the normalized currents i1 /i0 (top) and i2 /i0 (bottom)
in dependence on the normalized skin depth δ/b for different conductivity
ratios κ1 /κ2
180 5 Quasi Stationary Fields – Eddy Currents
Calculate the current distribution inside the conductors. The influence of the junctions
in z = 0 and z = l on the field is negligible and the permeability μ inside the groove is
constant.
R E
z = l
y
¹ ® ¥ k , ¹ k = 0 , ¹ k , ¹ k = 0 , ¹
0 a b c x
z = 0
U R G
0
We suppose that the conductor x < a carries a current with complex amplitude i in
positive z-direction and the conductor b < x < c carries the reverse current in negative
direction. Because of the boundary conditions in y = 0 and y = h the magnetic
field is y-directed H(x) = ey H(x) and inside the conductors it satisfies the differential
equation
d2 H
= α2 H with α2 = jωκμ .
dx2
is given by ⎧
⎪ sinh(αx)
⎪
⎨ ; 0≤x≤a
i sinh(αa)
H(x) = .
h⎪⎪ sinh(α(x − c))
⎩ ; b≤x≤c
sinh(α(b − c))
In the space a ≤ x ≤ b between the conductors the magnetic field is homogeneous
H(x) = ey i/h.
Finally for the determination of the amplitude i the law of induction has to be analyzed.
∂ a
Eds = − Bd
∂t
C a
A suitable contour for the integration runs through the inner boundary x = a, b of the
conductors. This results in
dH(x)
E(x) = 1/κ rot [ey H(x)] = 1/κ ez = ez E(x)
dx
⎧
⎪ cosh(αx)
⎪
⎨ ; 0≤x≤a
i sinh(αa)
E(x) = αa .
κah ⎪
⎪ cosh(α(x − c))
⎩ ; b≤x≤c
sinh(α(b − c))
It follows
il i
αa [coth(αa) − coth(α(b − c))] + i (RE + RG ) + jωμ(b − a) l = U 0
κah h
and for the complex amplitude of the current applies
−1
i = U 0 [RE + RG + jωμ(b − a)l/h + α/κ [coth(αa) − coth(α(b − c))] l/h]
LA = μ(b − a)l/h
0
I m [ i ]= A
-0 .1
l= a = { 1 0 2; 1 0 3; 1 0 4; 1 0 5}
-0 .2
l= a
-0 .3
± = a
-0 .4
0 0 .2 0 .4 0 .6 0 .8 1
R e [ i ]= A
0
Im [ i ]= i 0
-0 .1
l= a = { 1 0 ; 1 0 4; 1 0 6}
-0 .2
l= a
-0 .3
± = a
-0 .4
-0 .5
0 0 .2 0 .4 0 .6 0 .8 1
R e [ i ]= i 0
b = a = { 1 ; 1 .2 ; 1 ,5 ; 2 ; 3 ; 5 ; 1 0 }
4 0
I m [ i ]= A ±
2 0
b = a
0
-2 0
-4 0
0 2 0 4 0 6 0 8 0
R e [ i ]= A
Fig. 5.12–2: See caption of Fig. 5.12–1, but RE = 1/(jωC) and C = 10−5 [F]
with a frequency-dependent resistance R(ω) and the reactance defines the inner self-
inductance Li (ω).
In the limit ω = 0 the current is continuous and described by
U0
i = ,
RE + RG + l(1/a + 1/(c − b))/(κh)
where we have assumed that the complex amplitudes RE , RG , and U 0 switch over to
the real valued amplitudes.
R R E z = l
k = 0 , ¹ S k = 0 , ¹
¹ ® ¥
h
k , ¹ k , ¹ k , ¹ k , ¹
k = 0 , ¹
0
x
-d -c -b -a a b c d ¹ ® ¥
z = c
R G U 0
Calculate the currents in the conductors, when the influence of the junctions in z = 0
and z = l on the field is negligible.
We suppose that the conductors a ≤ |x| ≤ b and c ≤ |x| ≤ d carry the currents i2
and i1 . Because of the high-permeable boundary and the geometry of the conductors
=
the magnetic field has only a y-component and depends just on the x-coordinate H
= ey H(x). Inside the conductors it satisfies the differential equation
d2 H 2j 1+j
= α2 H ; α2 = jωκμ = ⇒ α= .
dx2 δ2 δ
As the magnetic field vanishes in |x| > d and is homogeneous between the conductors
the ansatz for H(x) is
⎧
⎪ i1 sinh(α(|x| − d))
⎪
⎪ ; c ≤ |x| ≤ d
⎪
⎪ h sinh(α(c − d))
⎪
⎪
⎪
⎪ i1
⎨ ; b ≤ |x| ≤ c
H(x) = h
⎪
⎪ (i1 + i2 ) sinh(α(|x| − b)) i1 sinh(α(|x| − a))
⎪
⎪ + ; a ≤ |x| ≤ b
⎪
⎪ h sinh(α(a − b)) h sinh(α(b − a))
⎪
⎪
⎪
⎩ (i1 + i2 ) ; |x| ≤ a .
h
This approach already satisfies the continuity relations of the magnetic field.
The homogeneous field between the conductors is defined by the currents i1 and i2 ,
which are unknown so far. For the calculation of the currents we evaluate the law of
5.13 Magnetically Coupled System of Conductors 185
induction
s = −∂
Ed a;
Bd E = ez 1 dH = ez E(x)
= 1 rot H
∂t κ κ dx
C a
along the boundaries of the conductors. It is convenient to integrate along the inner
contour of both conductor pairs. The integration over x = ±c leads to
1 dH α (i1 + i2 ) cosh(α(a − b)) − i1
E(a) = =
κ dx a<x<b κh sinh(α(a − b))
x→a
b
1 − cosh(α(a − b)) cosh(α(b − a)) − 1
H(x)dx = (i1 + i2 )/h + i1 /h
α sinh(α(a − b)) α sinh(α(b − a))
a
2i1 + i2 1 − cosh(α(a − b))
=
hα sinh(α(a − b))
we can identify two equations for the determination of the complex amplitudes i1 and
i2 .
αl 2i1 2(i1 + i2 )
− 2i1 coth(α(c − d)) + jωμ (c − b) l + al +
hκ h h
2l 1 − cosh(α(a − b))
+ (2i1 + i2 ) + i1 (RG + RS ) = U 0
hα sinh(α(a − b))
|i1 |
i0
0 .1 5
0 .1 2 5
l= h = { 1 0 ; 1 0 2; 1 0 3}
0 .1
0 .0 7 5
l= h
0 .0 5
0 .0 2 5
0
0 0 .2 0 .4 0 .6 0 .8 ± = h
|i2 |
i0
l= h = { 1 0 ; 1 0 2; 1 0 3; 1 0 4; 1 0 5; 1 0 6}
0 .2 5
0 .2
0 .1 5
l= h
0 .1
0 .0 5
0
0 1 2 3 4 ± = h
Fig. 5.13–1: Magnitude of the normalized primary current i1 /i0 (top) and secondary
current i2 /i0 (bottom) in dependence on the normalized skin depth δ/h
with 2a/h = 1, RG = RS = 0, and RE = 1/(jωC) with C = 10−5 [F]
(i0 = i1 |ω→0 )
5.13 Magnetically Coupled System of Conductors 187
l= h = { 1 0 2; 1 0 3; 1 0 4; 1 0 5; 1 0 6}
Im [ i2 ]= i 0
0 .1
± = h
0 .0 5
-0 .0 5
-0 .1
l= h
-0 .1 5
-0 .2 5 -0 .2 -0 .1 5 -0 .1 -0 .0 5 0 R e [ i2 ]= i 0
|´ |
1 0
l= h = { 1 0 ; 1 0 2; 1 0 3; 1 0 4; 1 0 5; 1 0 6}
8
6
± = h
4
0
0 0 .5 1 1 .5 ± = h 2
Fig. 5.13–2: Parametric plot of the secondary current i2 /i0 (top) and magnitude of the
ratio |η| = |i2 /i1 | (bottom) in dependence on the normalized skin depth
δ/h
Parameters see Fig. 5.13–1
188 5 Quasi Stationary Fields – Eddy Currents
on planes x = ±b with b > a in z-direction. The space |x| > a has the permeability μ0 .
Find the magnetic field and the current density inside the conducting slab in case of
an arbitrary time-dependency of the exciting current sheets. The dissipated power in
the conductor for t ≥ 0 has to be determined, when direct currents are turned off at
the time t = 0. Compare the dissipated energy with the energy of the magnetic field
stored in the conducting material before the turn-off. Finally analyze the case of a
harmonic time-dependency.
The exciting field H E (x, t) of the current sheets without presence of the slab is homo-
geneous in |x| < b.
⎧
⎨ e H (t) = e K(t) ; |x| < b
y E y
H E (x, t) =
⎩ 0 ; |x| > b
K (t) -K (t)
-b -a a b x
¹ 0 k , ¹ ¹ 0
when the field in a < |x| < b has already vanished, but inside the slab remains un-
changed. With cos(pa) = 0 and pn a = (2n + 1)π/2 with n = 0, 1, 2, . . . it follows
⎧
∞
⎨ −μK ; |x| < a
0
pn cn cos(pn x) = .
⎩ 0 ; x = ±a
n=0
Now because of the orthogonality of the trigonometric functions the constants cn are
given by
a
sin(pn x) −2μK0 2μK0
pn cn a = −μK0 = sin(pn a) = − (−1)n
pn −a pn pn
2μK0
⇒ cn = (−1)n+1 .
p2n a
Hence the magnetic field is
∞
(−1)n κμ
H(x, t) = 2K0 cos(pn x) exp(−t/τn ) ; τn =
n=0
pn a p2n
∞
2K0
⇒ J(x, t) = (−1)n+1 sin(pn x) exp(−t/τn ) .
a n=0
The dissipated power per unit lengths ly and lz in y- and z-direction after turning off
the currents is given by
a a
Pv (t) 1
= J(x, t)E(x, t)dx = J 2 (x, t)dx
ly lz κ
−a −a
∞ a ∞
Wv 4K02 2
= sin (pn x)dx exp(−2t/τn )dt
ly lz κa2 n=0
−a t=0
2 ∞ ∞ ∞ ∞
4K0 exp(−2t/τn ) 2K02 μ 2μK02 −2
= a = = p .
κa2 n=0 −2/τn 0 a n=0 p2n a n=0 n
5.14 Induced Current Distribution in a Conducting Slab 191
Before the turn-off the energy of the magnetic field stored in |x| < a per unit lengths
is
a ∞ a
Wm 1 μ 2 1
= Bdx
H = 4K0 2
cos2 (pn x)dx
ly lz 2 2 n=0
(p n a)
−a −a
∞
2μK02 Wv
= p−2
n =
a n=0
ly lz
This result is also obtained from the differential equation for the complex amplitude.
d2 H cosh(αx)
H(x) = ey H(x) ; = α2 H ; H(x) = C
dx2 cosh(αa)
192 5 Quasi Stationary Fields – Eddy Currents
H 1
K 0
t= ¿ 0
0 .8
0 .6
0 .4
0 .2
0
-1 -0 .5 0 0 .5 x = a 1
3
J
K 0 = a
2
0
t= ¿ 0
-1
-2
-3
-1 -0 .5 0 0 .5 x = a 1
Fig. 5.14–1: Magnetic field H/K0 and current density J/(K0 /a) at different times t/τ0
after turning off the current sheets
t/τ0 = 2·10−4 ; 5·10−4 ; 10−3 ; 2·10−3 ; 5·10−3 ; 10−2 ; 2·10−2 ; 5·10−2 ; 10−1 ; 2·
10−1 ; 5 · 10−1 ; 1; 2; 5; 10
5.14 Induced Current Distribution in a Conducting Slab 193
cosh(αx)
H(x) = K0 .
cosh(αa)
∞
cosh(αx)
= bn cos(pn x)
cosh(αa) n=0
leads to
a
cosh(αx) 2pn a
bn a = cos(pn x)dx ; bn = (−1)n
cosh(αa) (αa)2 + (pn a)2
−a
and thus ∞
(−1)n (pn a)
H(x) = 2K0 cos(pn x)
n=0
(αa)2 + (pn a)2
agrees with the previous result.
6. Electromagnetic Waves
Two ideal transmission lines with inductance L1 and capacitance C1 in z > 0 and L2
and C2 in z < 0, per unit length respectively, that are of infinite length in one direction
on the z-axis, are interconnected by a switch S in z = 0 at the time t = 0. For times
t < 0 the voltage on the transmission line in z > 0 is U0 and on the transmission line
in z < 0 the voltage is zero.
Calculate the amplitudes of the partial waves the voltages u1,2 (z, t) and the currents
i1,2 (z, t) in z > 0 and z < 0 after the switch has been turned on. Furthermore check
the invariance of the total energy for times t > 0.
After the interconnection at the time t = 0 the partial waves u1 (z − v1 t) ; i1 (z − v1 t)
and u2 (z + v2 t) ; i2 (z + v2 t) propagate in positive and negative direction with phase
velocities v1,2 = (L1,2 C1,2
)−1/2 starting in z = 0. For the amplitudes of the partial
waves the relations
0 0
u1 /i1 = Z1 = L1 /C1 ; −u2 /i2 = Z2 = L2 /C2
with ⎧
⎨ 0 ; x≤0
σ(x) =
⎩ 1 ; x > 0.
u (z , t)
U 0
v 1
Z 2
U 0
v Z 1 + Z 2
2
z
v 2 t v 1 t
-i(z , t)
U 0
v 2 Z + Z
1 2 v 1
The energy We (t) + Wm (t) of the electromagnetic field stored in the section
−v2 t < z < v1 t must be equal to the electrostatic energy, that has been stored in
0 < z < v1 t before the interconnection. The energy We (t) of the electric field for t > 0
is
and this is also the energy, that has been stored in the electrostatic field in 0 < z < v1 t
before the interconnection.
−1/2
We (t < 0)|0<z<v1 t = 1/2 C1 U02 (L1 C1 ) t = 1/2 U02 /Z1 t
1 1
H = H TE + HTM = rot rot F T E + rot A
TM
jωμε μ
1
E = E T E + E T M = −1/ε rot F T E + rot rot A TM
jωμε
1
E TM = rot rot A TM =
jωμε
1 ∂ 2 AT M ∂ 2 AT M M
2
= ex + ey + ez β 2 − βzTmn AT M
jωμε ∂x∂z ∂y∂z
6.2 Excitation of Hybrid Waves in a Rectangular Waveguide 197
¹ , "
I 0 c o s(¼ y = b )
a = 2 a x
1 1 ∂F T E ∂F T E
E T E = − rot F T E = ey − ex
ε ε ∂x ∂y
TE = 1
H rot rot F T E =
jωμε
1 ∂2F T E ∂ 2F T E 2
T E
2
= ex + ey + ez β − βzmn F TE
jωμε ∂x∂z ∂y∂z
β 2 = ω 2 με ; βz2mn = β 2 − βmn
2
Δ AT M + β 2 AT M = 0 ; Δ F T E + β2 F T E = 0 .
At the perfect conducting boundary with contour C and surface normal n the tangential
component of the electric field E and the normal component of the magnetic field
= μH
B vanish. This requirement leads to the following conditions for the potentials.
∂F T E
A T M |C = 0 ; = 0
∂n C
The excited waves propagate in positive and negative z-direction starting from the
point of excitement in z = 0. With respect to the boundary conditions the solution of
198 6 Electromagnetic Waves
As the eigenvalues of the rectangular waveguide are degenerated applies for the prop-
agation constants
" #1/2
βzTmn
E
= βzTmn
M
= βzmn = β 2 − βmn
2
; 2
βmn = (m π/a)2 + (n π/b)2 .
= (β 2 − β 2 ) F /(jωε)
ez B zmn TE
is continuous in z = 0. Hence it is
A(1) (2)
mn = Amn = Amn .
C (1) (2)
mn = −C mn = C mn .
The continuity of the tangential component of the electric field E in z = 0 gives the
same result. Finally the last boundary condition in z = 0 is
ez × H z>0 − H
z<0 = K(x, y) = ey K(x, y)
z→0 z→0
H y z>0 − H y z<0 = 0 ; H x | z>0 − H x | z<0 = K(x, y)
z→0 z→0 z→0 z→0
6.2 Excitation of Hybrid Waves in a Rectangular Waveguide 199
1 ∂AT M 1 ∂2F T E
Hy = − +
μ ∂x jωμε ∂y∂z
1 ∂AT M 1 ∂2F T E
Hx = + .
μ ∂y jωμε ∂x∂z
2 mπ 2nπjβzmn
cos(mπ x/a) sin(nπ y/b) − C mn + Amn = 0
m,n
μ a jωμεb
βzmn n a
⇒ C mn = A
ωε m b mn
and with it the evaluation of the x-component results in
1 nπ 2mπjβzmn
sin(mπ x/a) cos(nπ y/b) 2C mn + Amn = K(x, y)
m,n
μ b jωμεa
2 βzmn mπ " #
Amn 1 + (n/m)2 (a/b)2 sin(mπx/a) cos(nπ y/b) =
μ m,n ωε a
Thus because of the given exciting current only solutions with n = 1 exist and with
the orthogonality of the trigonometric functions we obtain after multiplication with
sin(kπ x/a) and integration from x = 0 to x = a
⎧
2 βzm1 mπ " # ⎨ 0 ; m = 2k
Am1 1 + (a/b)2 /m2 a/2 = I0 sin(mπ/2) = I0
μ ωε a ⎩ (−1)k ; m = 2k + 1
with k = 0, 1, 2, . . ..
μI0 (−1)k ωε 1
A2k+1,1 = ; A2k,1 = 0
π 1 + (a/b)2 /m2 βz2k+1,1 2k + 1
K(y, t) = ex K0 cos(πy/a) cos(ωt) = ex Re {K0 cos(πy/a) exp(jωt)}
in (z = b ; 0 < y < a). Inside the waveguide (0 < y < a ; z > 0) the permittivity is ε
and the permeability is μ.
K
¹ , " ¹ , "
b z
The field inside the waveguide is divided into transverse electric (TE-) and transverse
magnetic (TM-)Waves with respect to the z-direction and is thus described by the
vector potentials
2
A TM =
ez AT M ; ΔA T M + β AT M = 0
F T E = ez F T E ; ΔF T E + β 2 F T E = 0
1 1 1 ∂AT M ∂AT M
H TM = rot (A TM ) = grad A TM ×
e z = −
e y +
e x
μ μ μ ∂x ∂y
1 1 ∂ 2 AT M ∂ 2 AT M 2 T M2
E TM = rot rot A TM =
e x +
e y +
e z (β − β z )A TM
jωμε jωμε ∂x∂z ∂y∂z
6.3 Excitation of Transverse Electric Waves in a Parallel-Plate Waveguide 201
TE 1 1 ∂F T E ∂F T E
E = − rot F T E = ey − ex
ε ε ∂x ∂y
1 1 ∂2F T E ∂2F T E 2 T E2
H TE = rot rot F T E = ex + ey + ez (β − βz )F T E .
jωμε jωμε ∂x∂z ∂y∂z
As the field is independent of the x-coordinate, the transverse magnetic field has only
a x-component and thus is not excited by a current sheet in x-direction. Hence the
resulting field is purely transverse electric with respect to the z-direction.
At the perfect conducting boundary the tangential component of the electric field
and the normal component of the magnetic flux density vanish. These conditions are
satisfied, when the normal derivative of the potential F T E vanishes at the boundaries
y = 0, y = a:
∂F T E ∂F T E
= = 0,
∂y y=0 ∂y y=a
and in z = 0 the potential itself vanishes:
F T E |z=0 = 0 .
If only the first condition is met, the solution of the differential equation for the vector
potential F T E reads
∞
F TE = TE
C n cos(nπy/a) exp(∓jβzn z) ; TE
βzn = β 2 − (nπ/a)2 = βzn .
n=1
The continuity of the z-component of the magnetic flux density in z = b requires the
continuity of the potential F T E , and hence C (1) (2)
n = C n = C n.
This constant follows from the boundary condition for the tangential component of the
magnetic field at the position of the current sheet.
1 ∂ 2
F ∂ 2
F
H y z<b − H y z>b = K0 cos(πy/a) = T E − T E =
z→b z→b jωμε ∂y∂z z<b ∂y∂z z>b
z→b z→b
202 6 Electromagnetic Waves
0 b z
y
0 b z
Fig. 6.3–1: Magnetic lines of force below the cut-off (ω = 1011 [s−1 ], top) and for an
excitation of one propagating mode
(ω = 2 · 1011 [s−1 ], bottom)
Parameters: a = 0.01 [m], b/a = 0.8, ε = ε0 , and μ = μ0
6.3 Excitation of Transverse Electric Waves in a Parallel-Plate Waveguide 203
0 b z
y
0
b z
Fig. 6.3–2: Magnetic lines of force when two (ω = 4 · 1011 [s−1 ], top) and three (ω =
6 · 1011 [s−1 ], bottom) propagating modes are excited, Parameters: see Fig.
6.3–1
204 6 Electromagnetic Waves
∞
1 nπ exp(jβzn b)
= − C βzn sin(nπy/a) = K0 cos(πy/a)
jωμε n=1 n a sin(βzn b)
Now, due to the orthogonality of the trigonometric functions the multiplication with
sin(kπy/a) and integration from y = 0 to y = a leads to
a
1 kπ exp(jβzk b) a
− C k βzk = K0 · cos(πy/a) sin(kπy/a)dy
jωμε a sin(βzk b) 2
0
& '( )
= 0 for k = 2m + 1 and
= a/π2k/(k 2 − 1) for k = 2m
and thus
4 ωμεa sin(βzn b) exp(−jβzn b)
C n = −jK0 ; n = 2m; m = 1, 2, 3, . . . .
π 2 βzn n2 − 1
So only modes with even modal numbers n = 2m are excited. The fundamental TE-
mode with n = 1 is not excited by the current sheet.
A coaxial cable with perfect conductors, inner radius a, and outer radius b > a is
of infinite extension in the negative z-direction and is short circuited in z = c. The
permittivity is ε0 in −∞ < z < 0 and ε in 0 ≤ z < c. The permeability μ is constant.
In z < 0 a TEM-wave with the electric field
= e E a/ exp(−jβ0 z) ; √
E e 0 β0 = ω με0
is excited.
Calculate the resulting field and find the voltage u(z) and the current i(z) on the cable.
Choose a value c such that the voltage u(z) at position z = −λ0 = −2π/β0 vanishes.
The magnetic field of the excited wave is
H e = eϕ E 0 /Z0 a/ exp(−jβ0 z) ;
e = 1/Z0ez × E Z0 = μ/ε0 .
An ansatz for the resulting electric field, that accounts for the reflections at z = 0 and
z = c, is
⎧
⎨ exp(−jβ0 z) + A exp(jβ0 z) ; z ≤ 0
z) = e E a
E( , 0
⎩ B exp(−jβz) + C exp(jβz) ; 0 ≤ z ≤ c
6.4 Coaxial Cable with Inhomogeneous Dielectric 205
b
a
c z
√
with waves, that propagate in ±z-direction with phase constants β0 and β = ω με.
= ±1/Z ez × E
For each partial wave holds the fundamental relation H and thus it is
⎧
E a ⎨ exp(−jβ0 z) − A exp(jβ0 z) ; z≤0
H( , z) = eϕ 0
Z0 ⎩ Z0 /Z(B exp(−jβz) − C exp(jβz)) ; 0 ≤ z ≤ c
with Z = μ/ε for 0 < z < c.
The unknown constants follow from the boundary conditions
ez × E = 0; E( , z)| z<0 = E( , z)| z>0 ; H( , z)| z<0 = H( , z)| z>0 .
z=c z→0 z→0 z→0 z→0
r − exp(−2jβc) 2Z
A = ; B = (1 − r exp(−2jβc))−1
1 − r exp(−2jβc) Z0 + Z
−2Z exp(−2jβc) Z − Z0
C = ; r = .
Z0 + Z 1 − r exp(−2jβc) Z + Z0
This is the solution for the total reflection of incident waves at a perfect conducting
wall in z = 0. With ε = ε0 it is r = 0 and
Now this describes the total reflection at the plane z = c in a coaxial cable with
homogeneous permittivity ε0 .
The voltage and the current on the cable are
b b
u(z) = z) e d =
E( , E( , z) d =
a a
⎧
⎨ exp(−jβ z) + A exp(jβ z) ; z ≤ 0
0 0
= u0 ln(b/a)
⎩ B exp(−jβz) + C exp(jβz) ; 0 ≤ z ≤ c
2π
i(z) =
H( , z)eϕ dϕ
0
⎧
⎨ exp(−jβ z) − A exp(jβ z) ; z≤0
0 0
= 2πi0
⎩ Z0 /Z(B exp(−jβz) − C exp(jβz)) ; 0≤z≤c
A = −1 = (r − 1)/(1 − r) .
A F T E = ez F T E
TM =
ez AT M ;
2 2
ΔAT M + β1,2 AT M = 0 ; ΔF T E + β1,2 F T E = 0,
6.5 Cylindrical Waveguide Resonator with Inhomogeneous Permittivity 207
" 2, ¹
a %
" 1,¹
-c
T M = rot A
B T M = −ez × grad AT M = μ H
TM
1
TM = 2 (1,2)2
E grad s (∂AT M /∂z) + ez (β1,2 − βT M )AT M
jωμε1,2
E T E = −1/ε1,2 rot F T E = 1/ε1,2
ez × grad F T E ; grad s = grad − ez ∂/∂z
B 1
2 (1,2)2
H TE =
TE
= grad s (∂F T E /∂z) + ez (β1,2 − βT E )F T E .
μ jωμε1,2
and accordingly the normal
At the cylinder wall = a the tangential component of E
component of B vanish. Therefore it is
AT M |=a = 0 ; ∂F T E /∂ |=a = 0
and in the domain (1) for −c < z < 0 and (2) for 0 < z < d the solution may be
written in the following form.
(1,2) (1,2) (1,2)
AT M ( , ϕ, z) = C T M UT Mnr ( ) φn (ϕ) ZT Mnr (z)
(1,2) (1,2) (1,2)
F T E ( , ϕ, z) = C T E UT Enr ( ) φn (ϕ) ZT Enr (z)
0 0 r = 1, 2, . . .
(1,2) 2 − (x /a)2 ; (1,2) 2 − (x /a)2 ;
βT Mnr = β1,2 nr βT Enr = β1,2 nr
n = 0, 1, 2, . . .
and
Also at the top and bottom wall z = d and z = −c the tangential component of E
accordingly the normal component of B vanish and
∂AT M /∂z|z=−c,d = 0 ; F T E |z=−c,d = 0
holds. The solution functions, that describe waves propagating in positive and negative
z-direction, have to be superposed with respect to these boundary conditions. It follows
(1) (1) (2) (2)
ZT Mnr (z) = cos(βT Mnr (z + c)) ; ZT Mnr (z) = cos(βT Mnr (z − d))
(1) (1) (2) (2)
ZT Enr (z) = sin(βT Enr (z + c)) ; ZT Enr (z) = sin(βT Enr (z − d)) .
In z = 0 the tangential components of E and H
are continuous.
(1) (1) (2) (2)
C T M cos(βT Mnr c) = C T M cos(βT Mnr d)
Therewith the eigenvalue equations for the TE- and TM-Modes follow from the division
of the related equations.
(1) (1) (2) (2)
βT Mnr tan(βT Mnr c) = −ε1 /ε2 βT Mnr tan(βT Mnr d)
(1)2 (2)2
βT Mnr − βT Mnr = β12 − β22 = ωnr
2
μ(ε1 − ε2 )
(1)
(1) βT Enr (2)
tan(βT Enr c) = − (2)
tan(βT Enr d)
βT Enr
(1)2 (2)2
βT Enr − βT Enr = β12 − β22 = ωnr
2
μ(ε1 − ε2 )
The determination of the resonance frequencies ωnr requires a numerical solution of
the transcendental eigenvalue equations.
6.6 Guided Waves in a Parallel-Plate Waveguide with Layered Permittivity 209
a b y
The guided waves are either transverse electric (TE) modes or transverse magnetic
(TM) modes. Therefore the field is described by the vector potentials
2
A TM =
ez AT M (y, z) ; ΔAT M + β1,2 AT M = 0
2 2
F T E = ez F T E (y, z) ; ΔF T E + β1,2 F TE = 0 ; β1,2 = ω 2 με1,2
with B T M = rot AT M and E T E = −1/ε rot F T E .
∂A 1 ∂ 2 AT M 2 2
B ex T M ;
TM = E TM = e y +
e z (β 1,2 − β TM )A TM
∂y jωμε1,2 ∂y ∂z
T E = − 1 ex ∂F T E ; 1 ∂2F T E 2
E TE =
H ey + ez (β1,2 − βT2 E )F T E
ε ∂y jωμε1,2 ∂y ∂z
210 6 Electromagnetic Waves
√ "√ #
ε1 /ε2 c2 − u2 tanh c2 − u2 (b/a − 1)
4 0
2 0
0
u ta n u
-2 0
-4 0
0 ¼ / 2 ¼ 3 ¼ / 2 2 ¼ u
1
¯ T M
¯ 1
0 .8
0 .6
0 1 0 2 0 3 0 ¯ 1a
√ "√ #
− c2 − u2 coth c2 − u2 (b/a − 1)
2 0
1 0
u c o t u
0
-1 0
-2 0
0 ¼ / 2 ¼ 3 ¼ / 2 2 ¼ u
1
¯ T E
¯ 1
0 .8
0 .6
1 0 2 0 3 0 ¯ 1a
a
2
a
2
0 a b
y
An approach for the solution of the differential equation for the TM-modes is
∂ 2 AT M 2 ∂AT M
+ (β1,2 − βT2 M ) AT M = 0 ; = 0; 2
β1,2 = ω 2 με1,2
∂y 2 ∂x
⎧ ⎫
⎨ C sin(β 2 − β 2 y) ⎬ 0<y<a
1 1 TM
AT M = exp(∓jβT M z); .
⎩ C sin( β 2 − β 2 (y − b)) ⎭ a<y<b
2 2 TM
6.6 Guided Waves in a Parallel-Plate Waveguide with Layered Permittivity 213
a
2
a
2
0 a b
y
This approach already satisfies the first boundary condition and the two other boundary
conditions require
∂AT M ∂AT M β12 − βT2 M β22 − βT2 M
= ; AT M = AT M y>a .
∂y y<a ∂y y>a ε1 y<a
y→a
ε2 y→a
y→a y→a
a
2
a
2
0 a b
y
a
2
a
2
0 a b
y
ε2 C 3 β12 − βT2 E sin( β22 − βT2 E (a − b))
=
ε1 C 4 β22 − βT2 E sin( β12 − βT2 E a)
216 6 Electromagnetic Waves
and
ε2 C 3 β12 − βT2 E cos( β22 − βT2 E (a − b))
= ,
ε1 C 4 β22 − βT2 E cos( β12 − βT2 E a)
and leads to the eigenvalue equation
0 0 0 0
β12 − βT2 E cot( β12 − βT2 E a) = β22 − βT2 E cot( β22 − βT2 E (a − b)) ,
out of which the phase constants βT E of the TE-Modes follow. Again, with the sub-
stitutions c2 = a2 (β12 − β22 ) and u2 = a2 (β12 − (βT E )2 ) one gets
u cot(u) = − c2 − u2 coth c2 − u2 (b/a − 1) .
Four Hertzian dipoles with moments ez Is (s λ) are located at positions (±a; 0; 0)
and (0; 0; ±b) in the homogeneous space of permittivity ε and permeability μ.
Calculate the far-field of the group and find the zeros of the radiation pattern in the
plane x = 0 for b = 2λ.
z
P
I s
¹ ," r p
b
I s I s
x
-a a
I s
-b
The magnetic far-field H at the observation point rp of a single dipole ez Is located
0
at the point of origin is
0 = jβ Is ez × rp exp(−jβrp ) ; √
H β = ω με = 2π/λ .
4π rp rp
6.7 Group of Hertzian Dipoles 217
If the same dipole is located at position rq , then the excited far field is
jβ rp exp(−jβrp ) exp(jβrq rp /rp ) .
H = Is ez × exp(jβrq rp /rp ) = H 0
4π rp rp
With rq = ±ex a, rq = ±ez b, and
rp /rp = ex sin ϑ cos ϕ + ey sin ϑ sin ϕ + ez cos ϑ = er
E = −2Z er × H
= −Z er × H f (ϑ, ϕ)
0
jβ exp(−jβrp )
= −Z Is · er × (ez × er ) · f (ϑ, ϕ)
2π & '( ) rp
= −eϑ sin ϑ
z
3 0 ° 6 0 °
1
0 .5
0 1 2 3 x
0 .5
1
1 5 0 ° 1 2 0 °
z
3 0 ° 6 0 °
1 .5
0 .5
0 1 3
2 y
0 .5
1 .5
1 5 0 ° 1 2 0 °
Fig. 6.7–1: Vertical radiation pattern Φ(ϑ, ϕ)/Φ0 for ϕ = 0 (top) and
ϕ = π/2 (bottom) with βb = 4π, βa = 2π, and Φ0 = Z/2(|I|s/λ)2
6.8 Linear Antenna in Front of a Conducting Plane 219
y
1 2 0 ° 6 0 °
1 5 0 ° 3 0 °
2
4 3 2 1 1 2 3 4 x
2
2 1 0 ° 3 3 0 °
2 4 0 ° 3 0 0 °
Fig. 6.7–2: Horizontal radiation pattern Φ(ϑ, ϕ)/Φ0 for ϑ = π/2 and βa = 2π
A linear antenna of length 0 < z < 3/4λ is positioned in front of a perfect conducting
plane in z = 0 and carries the current
√
i(z, t) = Re {i0 cos(2π z/λ) exp(jωt)} ; 2π/λ = β = ω με .
Calculate the far-field and the radiation pattern in z > 0, where the permittivity is ε
and the permeability is μ.
At the perfect conducting boundary the tangential component ez × E| z=0 vanishes
and it is thus possible to apply the method of images to describe the field in z > 0.
For this purpose the given current i(z, t) on the z-axis is extended to the total length
−3/4λ < z < 3/4λ and the antenna is assumed to be in the homogeneous space with
material properties μ, ε.
220 6 Electromagnetic Waves
3 ¸
4
i0 c o s(2 ¼ z = ¸ )
r
#
¹ ,"
k ® ¥
According to the specified current in −λ/4 < z < λ/4 (λ/2−dipole) this leads the
.
following magnetic far-field H D
λ/4
(rp ) = j β i
H ez ×
rp exp(−jβrp )
cos(βzq ) exp(jβzq cos ϑ) dzq
D
4π 0 rp rp
−λ/4
The resulting field H G of the antenna with the current i(z, t) in −3/4λ < z < 3/4λ
is the superposition of the fields of three λ/2-dipoles with midpoints at zq = 0, ±λ/2,
taking into account the correct phase of the currents.
6.8 Linear Antenna in Front of a Conducting Plane 221
3 0 °
1 .5
1
6 0 °
0 .5
0 0 .5 1 1 .5 9 0 °
0 .5
1 2 0 °
1
1 .5
1 5 0 °
= eϕ H [1 − exp(jβλ/2 cos ϑ) − exp(−jβλ/2 cos ϑ)] = eϕ H [1 − 2 cos(π cos ϑ)]
H G D D
is
Finally the time-averaged Poynting vector S
− = 1 Re E
S ∗ = 1 Z|H
×H |2 er = er S−
G G G
2 2
222 6 Electromagnetic Waves
¹ ," z
2 a 2 b
x
2 c
(rp ) = 2H
H cos(βb sin ϑ cos ϕ) .
2 0
where r0 is an arbitrary constant and H E describes the specific field of the source.
Thus, related to the far-field the system of antennas behaves like a point source in the
origin. The radiation pattern is given by
1 1
Φ0 (ϑ, ϕ) = Z|H 0 (rp )|2 rp2 = Z|H rp /rp )|2 r02 .
E (
2 2
Find the magnetic far-field and the radiation pattern when four systems are positioned
equidistantly on the x-axis. The distance between two adjacent systems is 2a.
where rqi points to the position of the i-th antenna. Thus for the present problem we
get
rp ) = H
H( (rp ) exp(jβex arp /rp ) + exp(−jβex arp /rp )+
0
+ exp(jβex 3arp /rp ) + exp(−jβex 3arp /rp ) .
where
2
ΦG4 (ϑ, ϕ) = 2 cos(βac) + 2 cos(β3ac)
An alternative description of the magnetic field follows from the approach of prob-
lem 6.9. In the first step only two antennas are combined
(rp ) = H
H (rp ) exp(jβac) + exp(−jβac) = H (rp ) 2 cos(βac) .
1 0 0
Finally the new combination of two combined antennas with distance 4a leads to the
far-field
rp ) = H
H( 1 (rp ) 2 cos(2βac)
1 (rp ) exp(j2βac) + exp(−j2βac) = H
where
ΦG2 (ϑ, ϕ, d) = 4 cos2 (1/2βdc)
x
a
A rectangular waveguide with perfectly conducting boundaries and cross sectional di-
mensions a and b < a is of infinite length in z-direction. In z < 0 a dielectric material
with permittivity ε1 and in z > 0 material with permittivity ε2 < ε1 is inserted. The
permeability μ is constant.
Find the function, that specifies the dependence on the z-coordinate of the Modes
in z > 0, when in z < 0 the fundamental mode is excited, propagating in positive
z-direction at the frequency
√
ω = 1/2 (ωc1 + ωc2 ) = 1/2 ωc1(1 + ε1 /ε2 ) ; ωc1,2 = π/(a με1,2 ).
6.11 Waveguide with Sections of Different Dielectrics 225
(1)
The phase constants βzmn of the modes in z < 0 are
0 0
βz(1)
mn
= β12 − ((mπ/a)2 + (nπ/b)2 ) = β12 − βmn
2 ; β12 = ω 2 με1 .
For the fundamental mode holds m = 1 and n = 0 because of a > b and thus the
cut-off frequency ωc1 is
√
βz(1)
10
= 0; β12 = β10
2
= (π/a)2 = ωc1
2
με1 ; ωc1 = π/(a με1 ) .
Hence the excitation of the waves in z < 0 occurs at a frequency ω greater than the
cut-off limit ωc1 , as long as ε1 > ε2 .
2
βz(2)
mn
= 2
1/4 ωc1 2 .
1 + ε1 /ε2 με2 − βmn
(2) √
For the fundamental mode βz10 applies with β10 = π/a = ωc1 με1
0
βz(2)
10
= 2 (1 +
1/4 ωc1 2 με
ε1 /ε2 )2 με2 − ωc1 1
= ωc1 με1 1/4(1 + ε1 /ε2 )2 ε2 /ε1 − 1
0
βz(2)
10
= π/a 1/4 ε2/ε1 + 1/2 ε2 /ε1 − 3/4 .
Now the radicand is negative because of ε2 /ε1 < 1 and thus in z > 0 the fundamental
mode and all other modes decay exponentially with increasing distance z. For the
fundamental mode holds
0
exp −π/a 3/4 − 1/4 ε2/ε1 − 1/2 ε2 /ε1 z .
Calculate the field in z ≤ d and the current sheet in z = d. Use your solutions to
analyze the limits κ = 0, μ = μ0 , and ε = ε0 .
¹ 0," 0 k ,¹ ," k ® ¥
d z
E 0 e x p (-j ¯ 0z )
= ex E(z) in z ≤ d is
The approach for the complex amplitude of the electric field E
⎧
⎨ E exp(−jβ z) + E exp(jβ z) ; z ≤ 0
0 0 02 0
E(z) =
⎩ E exp(−jγz) + E exp(jγz) ; 0 ≤ z ≤ d
11 12
0 0
γr = β 1/2( 1 + (κ/(ωε))2 + 1) ; γi = β 1/2( 1 + (κ/(ωε))2 − 1) .
For the complex amplitude of the magnetic flux density B = μH applies B = ±1/ω γ ×
where the leading sign denotes waves propagating in positive or negative z-direction.
E,
6.12 Reflection of a Plane Wave at a Conducting Half-Plane 227
= ey H(z) is
Consequently, the complex amplitude H
⎧
⎨ 1/Z (E exp(−jβ z) − E exp(jβ z)) ; z ≤ 0
0 0 0 02 0
H(z) =
⎩ γ/(ωμ)(E exp(−jγz) − E exp(jγz)) ; 0 ≤ z ≤ d
11 12
with Z0 = μ0 /ε0 .
The unknown complex amplitudes E 02 , E 11 , and E 12 follow from the boundary condi-
tions in z = 0 and z = d.
ez × E =0;
ez × E =
e z ×
E ;
e z ×
H =
e z ×
H
z<0 z>0 z<0 z>0
z=d z→0 z→0 z→0 z→0
1+r
E 02 = r E0 ; E 11 = E0
1 − exp(−2jγd)
exp(−2jγd)
E 12 = −E0 (1 + r) ; r = (1 − q)/(1 + q)
1 − exp(−2jγd)
The current sheet K
= −ez × H = ex K is
z=d
1+r γ
K = γ/(ωμ)(E 11 exp(−jγd) − E 12 exp(jγd)) = −j E0 .
sin(γd) ωμ
With κ = 0 it is γi = 0 and γ = β and furthermore
Z0 cos(βd)
q = ; Z = μ/ε .
jZ sin(βd)
If, in addition, the limits μ = μ0 and ε = ε0 are built, then with Z = Z0 and β = β0
the reflection factor is r = − exp(−2jβ0 d). The complex amplitudes are
E (z , t)
! t = 0
-4 d -3 d -2 d -d 0 d
H (z , t)
! t = 0
-4 d -3 d -2 @ -d 0 d
Fig. 6.12–1: Graphs of the electric and magnetic field as a function of the z-coordinate
at the time t = 0. Additionally the partial waves propagating in different
directions are plotted separately.
μ = μ0 ; ε = ε0 ; κ = 2 [S/m]; d/|γ| = 2.42 · 10−5
6.12 Reflection of a Plane Wave at a Conducting Half-Plane 229
E (z , t)
! t = ¼ = 2
-4 d -3 d -2 d -d 0 d
H (z , t)
! t = ¼ = 2
-4 d -3 d -2 d -d 0 d
Fig. 6.12–2: Graphs of the electric and magnetic field as a function of the z-coordinate
at the time ωt = π/2. Additionally the partial waves propagating in
different directions are plotted separately. Parameters: see Fig. 6.12–1
230 6 Electromagnetic Waves
E (z , t)
-4 d -3 d -2 d -d 0 d
H (z , t)
-4 d -3 d -2 d -d 0 d
Fig. 6.12–3: Graphs of the electric and magnetic field as a function of the z-coordinate
at various times t. Parameters: see Fig. 6.12–1
6.13 Guided Waves in a Dielectric Slab Waveguide 231
A dielectric slab a < x < b with permittivity ε1 is positioned beside a perfectly con-
ducting wall x ≤ 0. Outside of the slab the permittivity is ε0 and the permeability μ is
constant. The slab is capable to guide both transverse magnetic (TM-) and transverse
electric (TE-)waves in the z-direction.
Find the eigenvalue equations for the determination of the phase constants of the
modes, that are independent of the y-coordinate.
a b x
TM-Waves (-Modes)
T M = ez AT M (x, z) with
These modes can be derived from the vector potential A
B ey μH T M (x, z) = −ey ∂AT M /∂x .
T M = rot [
ez AT M (x, z)] = μH TM =
Obviously the magnetic field has only one component, that satisfies the differential
equation
2 ∂2H T M ∂2H T M 2 2
ΔH T M + β0,1 HTM = 2
+ + β0,1 H T M = 0; β0,1 = ω 2 με0,1 .
∂x ∂z 2
An approach for H T M , that describes wave propagation in ±z-direction with the phase
constant βzT M , is
H T M (x, z) = H 0 X 1,2 (x) exp(∓jβzT M z)
232 6 Electromagnetic Waves
d2 X 1 2
+ (β02 − βzT M ) X 1 = 0 ; 0 < x < a; x>b
dx 2 & '( )
−α2
d2 X 2 2
+ (β12 − βzT M ) X 2 = 0 ; a < x < b.
dx2 & '( )
βx2 > 0
The phase constants βzT M of each mode are of the same value in any region of the
half-space x > 0 and lie within an interval defined by the material parameters ε0 and
ε1 .
2 2
β0 < βzT M < β1 ; α2 = βzT M − β02 > 0 ; βx2 = β12 − βzT M > 0
Hence the complex amplitude H T M (x, z) takes the form
⎧ ⎫
⎪
⎪ C 1 cosh(αx)/ cosh(αa) ⎪ ⎪ 0≤x≤a
⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪
⎪
⎨ C sin(β (x − a))+ ⎪
⎬
2 x
H T M (x, z) = H 0 exp(∓jβzT M z) ; a ≤ x ≤ b
⎪
⎪ +C cos(β (x − a)) ⎪
⎪
⎪
⎪ 3 x ⎪
⎪
⎪
⎪ ⎪
⎪
⎪
⎩ C exp(−α(x − b)) ⎪
⎭
4 x≥b
1 1 ∂H T M ∂H T M
⇒ ET M = rot H T M = −ex + ez .
jωε jωε ∂z ∂x
The field must satisfy the boundary conditions
ez E TM = 0; H T M | 0≤x≤a = H T M | a≤x≤b
x=0 x→a x→a
H T M | a≤x≤b = H T M | x≥b ; ez E
T M 0≤x≤a =
e z E T M a≤x≤b
x→b x→b
x→a x→a
ez E T M a≤x≤b = ez E T M x≥b .
x→b x→b
The first one is already satisfied by the given approach and the other ones lead to the
following equations.
"
# √
c2 −v 2
√
tan v ab − 1 · εε01 √c2v−v2 − ε1
ε0 v tanh c2 − v 2
1 5
-1 5
0 ¼ = 4 ¼ = 2 3 ¼ = 4 ¼ 5 ¼ = 4 v
T M
¯ z
¯ 1
0 .9
0 .8
0 .7
0 1 0 2 0 3 0 4 0 ¯ 1a
Fig. 6.13–1: Determination of the eigenvalues at ω = 15·1010 [s−1 ] and phase constants
as a function of the frequency for TM-modes with ε1 /ε0 = 2 and b/a = 2
v 2 = a2 (β12 − (βzT M )2 ); c2 = a2 (β12 − β02 )
234 6 Electromagnetic Waves
-a
0 a b x
with
0 0
(αa)2 + (βx a)2 = (β1 a)2 − (β0 a)2 ; βzT M = β12 − βx2 = α2 + β02 .
-a
0 a b x
TE-Waves (-Modes)
The vector potential of the field is F T E = ez F T E (x, z) with
E T E = −1/ε rot F T E = ey E T E (x, z) = ey 1/ε∂F T E /∂x.
Now the electric field has only a y-component, that satisfies
2 ∂2ET E ∂ 2E T E 2
ΔE T E + β0,1 ET E = + + β0,1 ET E = 0 .
∂x2 ∂z 2
236 6 Electromagnetic Waves
"
# √ 2 2 √
tan v ab − 1 · − c v−v + √ v
c2 −v 2
tanh c2 − v 2
2
-1
0 ¼ = 8 ¼ = 4 v
¯ T E
z
¯ 1
0 .9
0 .8
0 .7
0 1 0 2 0 3 0 4 0 ¯ 1a
Fig. 6.13–4: Determination of the eigenvalues at ω = 3 · 1010 [s−1 ] and phase constants
as a function of the frequency for TE-modes with ε1 /ε0 = 2 and b/a = 2
v 2 = a2 (β12 − (βzT E )2 ); c2 = a2 (β12 − β02 )
6.13 Guided Waves in a Dielectric Slab Waveguide 237
-a
0 a b x
Thus again the approach for the field is E T E = E 0 X 1,2 (x) exp(∓jβzT E z) where the
functions X 1,2 (x) and the constants have to be redefined.
⎧ ⎫
⎪
⎪ C 1 sinh(αx)/ sinh(αa) ⎪ ⎪ 0≤x≤a
⎪
⎪ ⎪
⎪
⎪
⎪ ⎪
⎪
⎪
⎨ C sin(β (x − a))+ ⎪
⎬
2 x
E T E (x, z) = E 0 exp(∓jβzT E z) ; a ≤ x ≤ b
⎪
⎪ +C cos(β (x − a)) ⎪
⎪
⎪
⎪ 3 x ⎪
⎪
⎪
⎪ ⎪
⎪
⎪
⎩ C exp(−α(x − b)) ⎪
⎭
4 x≥b
2 2
α2 = βzT E − β02 > 0 ; βx2 = β12 − βzT E > 0
238 6 Electromagnetic Waves
-a
0 a b x
1 1 ∂E T E ∂E T E
H T E (x, z) = − rot (ey E T E ) = ex − ez
jωμ jωμ ∂z ∂x
The boundary conditions are
E T E | 0<x<a = E T E | a<x<b ;
x→a
ez H
x→a
T E 0<x<a =
e z H T E a<x<b
x→a x→a
E T E | a<x<b = E T E | x>b ; T E
ez H =
e z
T E
H
x→b x→b a<x<b x>b
x→b x→b
C 1 α coth(αa) = βx C 2
1 + tanh(αa)
= −α/βx tan(βx (b − a))
1 − (βx /α)2 tanh(αa)
with
0 0
(αa)2 + (βx a)2 = (β1 a)2 − (β0 a)2 ; βzT E = β12 − βx2 = α2 + β02
for the determination of the finite number of solutions for βx , βzT E , and α.
√ 2
b c − v2 v
2 2
1 + tanh c − v = tan v −1 · − +√ 2 2
tanh c − v .
a v c2 − v 2
A dielectric slab waveguide consists of a core in |x| < a with permittivity ε2 and a
cladding in a < |x| < b with permittivity ε1 < ε2 . The surrounding space has the
permittivity ε0 and the permeability μ is constant. In the waveguide both transverse
magnetic (TM-)waves and transverse electric (TE-)waves can propagate in z-direction,
while there is no dependence on the y-coordinate.
Find the eigenvalue equations for the guided modes whose magnetic (TM) or electric
(TE) field is either an even or odd function of the x-coordinate.
T M = ez AT M (x, z)
A and F T E = ez F T E (x, z).
As in the present case the electric field (TE-waves) or the magnetic field (TM-waves)
has only a y-component, it is useful to solve the field equations directly.
240 6 Electromagnetic Waves
a b x
TM-Waves
The wave propagation with the unknown phase constant βzT M is described by
T M = ey H T M (x, z) ; 2
H ΔH T M + β0,1,2 HT M = 0
∂2H T M ∂2H T M 2
+ + β0,1,2 HTM = 0
∂x2 ∂z 2
H TM
T M (x, z) = H 0 X 0,1,2 (x) exp(∓jβz z)
d2 X 0,1,2 2 2
2
+ (β0,1,2 − βzT M )X 0,1,2 = 0 ; β0,1,2 = ω 2 με0,1,2
dx2
1 1 ∂H T M ∂H T M
E TM = rot H TM = −
e x +
e z .
jωε jωε ∂z ∂x
In case of even modes the solution functions satisfy X(x) = X(−x) and the solution
reads as follows
√
"√
#
c2 − v 2 f 2 − v 2 + ε0
ε1 · c2 − v 2 · tanh c2 − v 2 ab − 1
3 0
2 0
1 0
-1 0
0 ¼ = 2 ¼ v
"√
# √
tan v ε1
ε2 ·v f 2 − v 2 · tanh c2 − v 2 ab − 1 + ε0
ε2 · v c2 − v 2
Fig. 6.14–1: Solutions of the eigenvalue equation for the determination of the phase
constants βzT M of the even TM-modes with
c = 4.622, ε2 /ε0 = 5, ε1 /ε0 = 2, b/a = 1.5, and a = 0.1 [cm]
with βzT M = β12 − βx1
2 = β22 − βx2
2 = α2 + β02 .
The other case of odd modes requires X(x) = −X(−x) and leads to
Here the constants C 1,2,3,4 are redefined in contrast to the first approach for even
modes.
The tangential components of the electric and magnetic field are continuous in x = a
and x = b. This requires
¯ T M
z
¯ 2
0 .9
0 .8
0 .7
0 1 0 2 0 3 0 4 0 5 0 ¯ 2a
Fig. 6.14–2: Phase constants of the first five odd TM-modes as a function of
√
β2 a = ω με2 a with ε1 /ε0 = 2, ε2 /ε0 = 5, b/a = 1.5, and a = 0.1 [cm]
2 a
-a
-2 a
-b -a 0 a b
x
Fig. 6.14–3: Electric lines of force of the second even TM-Mode with
ε1 /ε0 = 2, ε2 /ε0 = 10, b/a = 1.5, a = 0.1 [cm], and ω = 4.37 ·1011 [s−1 ]
After a similar calculation we obtain the eigenvalue equation for the odd TM-modes.
The equations are
2 a
-a
-2 a
-b -a 0 a b
x
Fig. 6.14–4: Electric lines of first of the first odd TM-mode with
ε1 /ε0 = 2, ε2 /ε0 = 10, b/a = 1.5, a = 0.1 [cm], and ω = 2.99 ·1011 [s−1 ]
2 a
-a
-2 a
-b -a 0 a b
x
Fig. 6.14–5: Magnetic lines of force of the second even TE-mode with
ε1 /ε0 = 2, ε2 /ε0 = 10, b/a = 1.5, a = 0.1 [cm], and ω = 3.58 ·1011 [s−1 ]
TE-Waves
The analog approach for the electric field of the TE-modes is
T E (x, z) = E 0 X 0,1,2 (x) exp(∓jβzT E z)
E
The functions X 0,1,2 remain unchanged but with new constants C 1,2,3,4 .
After evaluating the boundary conditions we finally get the eigenvalue equations
α/βx1 − βx2 /βx1 tan(βx2 a)
tan(βx1 (b − a)) = ; (even TE-modes)
1 + α/βx1 βx2 /βx1 tan(βx2 a)
246 6 Electromagnetic Waves
2 a
-a
-2 a
-b -a 0 a b
x
Fig. 6.14–6: Magnetic lines of force of the first odd TE-Mode with
ε1 /ε0 = 2, ε2 /ε0 = 5, b/a = 1.5, a = 0.1 [cm], and ω = 3.35 ·1011 [s−1 ]
A dielectric cylinder has the radius a and is of infinite length in z-direction. Its per-
meability and permittivity are μ and ε, whereas the surrounding medium is the free
space with material parameters μ0 and ε0 .
1. Calculate the electromagnetic field excited by a plane wave, when the phase
√
vector β0 = −ex β0 ; β0 = ω μ0 ε0 is directed perpendicular to the axis of the
cylinder. In the first case a) the wave is z–polarized and in the second case b)
y–polarized.
2. Find the electromagnetic field excited by a line source i(t) = Re {i0 exp(jωt)}
with distance c > a parallel to the axis of the cylinder.
¯ 0
¹ 0," 0
¹ ," ' i 0
c x
a
(2)
The diffracted wave is described by Hankel functions Hn (β0 ) and inside the cylinder
√
the field is described by Bessel functions Jn (β ) with β = ω με. Hence the resulting
electric field is
⎧ ⎫
∞
⎨ J (β ) + a H (2) (β ) ⎬ ≥a
n 0 n n 0
= ez E0
E (2−δ0n ) exp(jnπ/2) cos(nϕ) ; ,
⎩ b Jn (β ) ⎭ ≤a
n=0 n
and the continuity of the tangential component H ϕ of the magnetic field yields
1 ∂E ∂E
rot [ez E( , ϕ)] = grad E × ez = e − eϕ = −jωμ H
∂ϕ ∂
1 1 ∂E 1 ∂E
H = − ; Hϕ =
jωμ ∂ϕ jωμ ∂
1 ∂E 1 ∂E
=
μ ∂ <a μ0 ∂ >a
→a →a
β0
β
⇒ Jn (β0 a) + an Hn(2) (β0 a) = [bn Jn (βa)] .
μ0 μ
Thus, there are two equations for the determination of the constants an and bn .
−1 Z0
an = M n Jn (β0 a)Jn (βa) − Jn (β0 a)Jn (βa)
Z
bn = M −1 n Hn(2) (β0 a)Jn (β0 a) − Hn(2) (β0 a)Jn (β0 a)
Z0 (2)
with Mn = Hn(2) (β0 a)Jn (βa) − Hn (β0 a)Jn (βa)
Z
μ μ0
Z = ; Z0 =
ε ε0
6.15 Diffraction by a Dielectric Cylinder 249
In the special case of a perfect conducting cylinder the field in < a vanishes.
E z |=a = 0 ; H =a = 0
Jn (β0 a)
Jn (β0 a) + an Hn(2) (β0 a) = 0 ; an = − (2)
Hn (β0 a)
= ez K . It follows
Now the surface of the cylinder = a carries a current sheet K E E
H ϕ >a − H ϕ <a = K E
→a →a
1 ∂E
K E = H ϕ >a = =
→a jωμ0 ∂ >a
→a
∞
jβ0 Jn (β0 a)
= − E0 (2 − δ0n ) exp(jnπ/2) Jn (β0 a) − (2)
Hn(2) (β0 a) cos(nϕ)
ωμ0 n=0 Hn (β0 a)
Finally with
n (2) n (2)
Jn (z) = −Jn+1 (z) + Jn (z) ; Hn(2) (z) = −Hn+1 (z) + H (z)
z z n
Hn(2) (z) = Jn (z) − jNn (z) ;
2
Jn+1 (z)Nn (z) − Jn (z) Nn+1 (z) = (Wronskian determinant)
πz
∞
2E0 (2 − δ0n ) exp(jnπ/2)
⇒ KE = (2)
cos(nϕ) .
πβ0 aZ0 n=0 Hn (β0 a)
β ε0 Z
Jn (β0 a) + cn Hn(2) (β0 a) = d J (βa) = d J (βa) .
β0 ε n n Z0 n n
Hence the constants are
Z (2)
Mn = Hn(2) (β0 a) Jn (βa) − H (β0 a) Jn (βa)
Z0 n
Z
cn = M −1
n J (β
n 0 a)J
n (βa) − J n (βa)J
n 0(β a)
Z0
dn = M −1
n Hn(2) (β0 a)Jn (β0 a) − Hn(2) (β0 a)Jn (β0 a) .
K = eϕ K H ; H| <a − H| >a = K H
H →a →a
KH = − H|=a
∞
= H0 (2 − δ0n ) exp(jnπ/2) ·
n=0
Jn (β0 a)
· (2)
Hn(2) (β0 a)
− Jn (β0 a) cos(nϕ)
Hn (β0 a)
∞
2H0 (2 − δ0n ) exp(jnπ/2)
KH = j (2)
cos(nϕ) .
πβ0 a n=0 Hn (β0 a)
The vector potential of a line source in the homogeneous space with μ0 , ε0 is:
= ez A ( , ϕ) ;
A0 0 ΔA0 + β02 A0 = 0 ; = c
6.15 Diffraction by a Dielectric Cylinder 251
y = a
-2
-4
-4 -2 0 2 4 x = a
Fig. 6.15–1: Diffraction of a plane, z-polarized wave: Magnetic lines of force for
λ0 /λ = 1.1 and λ/a = 2 (λ0 = 2 π/β0 ; λ = 2π/β)
⎧ ⎫
∞
⎨ J (β c) Hn(2) (β ) ⎬ ≥c
n 0 0
A0 ( , ϕ) = an cos(n(ϕ − ϕ0 )) ; .
⎩ Hn(2) (β0 c) Jn (β0 ) ⎭ ≤c
n=0
This approach already satisfies the continuity relation of the electric field and the
= rot A
continuity of the radial component of B at = c.
0 0
252 6 Electromagnetic Waves
y = a
-2
-4
-4 -2 0 2 4 x = a
1 ∂A0 ∂A0
H 0ϕ >c − H 0ϕ <c = K(ϕ) = −
→c →c μ0 ∂ <c ∂ >c
→c →c
β0 (2)
∞
K(ϕ) = an Hn (β0 c)Jn (β0 c) − Jn (β0 c) Hn(2) (β0 c) cos(n(ϕ − ϕ0 )) .
μ0 n=0
6.15 Diffraction by a Dielectric Cylinder 253
2 0
y = a
1 0
-1 0
-2 0
-2 0 -1 0 0 1 0 x = a 2 0
an (1 + δ0n ) μ0
= β0 an π(1 + δ0n )j [Jn+1 Nn − Jn Nn+1 ] = 2j = i
c c 0
254 6 Electromagnetic Waves
y = a
-2
-4
-4 -2 0 2 4 x = a
μ0 i0
⇒ an = −j .
2(1 + δ0n )
y = a
-1
-2
-2 -1 0 1 x = a 2
⎧ ⎫
⎨ β2 ⎬ > a ; = c
= ez A( , ϕ) ; 0
A ΔA + A = 0;
⎩ β2 ⎭ <a
256 6 Electromagnetic Waves
y = a
-2
-4
-4 -2 0 2 x = a 4
Fig. 6.15–6: Diffraction of the cylindrical wave excited by a line source at a dielectric
cylinder: Magnetic lines of force for λ0 /λ = 3, λ/a = 1, and c/a = 2
∞
μ0 1
A( , ϕ) = −j i0 cos(n(ϕ − ϕ0 ))·
2 n=0 1 + δ0n
⎧ ⎫
⎪
⎪ b J (β ) ; <a ⎪
⎪
⎪
⎨ n n ⎧ ⎫ ⎪
⎬
· ⎨ J (β c)H (2) (β ); >c ⎬ .
(2) n 0 0
⎪
⎪
n
; >a ⎪
⎪
⎪ cn Hn (β0 ) + ⎩ (2) ⎪
⎩ Hn (β0 c)Jn (β0 ); <c ⎭ ⎭
2 0
y = a
1 0
-1 0
-2 0
-2 0 -1 0 0 1 0 x = a 2 0
At first the continuity of the normal component of the magnetic flux density leads to
and a second equation results from the continuity of the tangential component of the
magnetic field:
β0 μ
bn Jn (βa) = cn Hn(2) (β0 a) + Hn(2) (β0 c) Jn (β0 a) .
β μ0
258 6 Electromagnetic Waves
3
K
K 0
2
-1
-¼ -¼ = 2 0 ¼ = 2 ' ¼
Fig. 6.15–8: Current sheet K/K0 on the surface of a conducting cylinder exposed to
the field of a line source at ( = c, ϕ0 = 0) with
K0 = −|i0 |/(πa), λ0 /λ = 1.1, and λ/a = 1/2
Hn (β0 c)
(2)
bn = Jn (β0 a)Hn(2) (β0 a) − Jn (β0 a) Hn(2) (β0 a) .
Mn
Applying the recurrence relations and the Wronskian-determinants they become
(2)
2 Hn (β0 c)
bn = −j
π β0 aM n
with
Z0
M n = Jn (βa)Hn(2) (β0 a) − J (βa) Hn(2) (β0 a) .
Z n
For the constants cn it follows
(2)
Hn (β0 c) Z0
cn = − Jn (βa) Jn (β0 a) − Jn (βa) Jn (β0 a) .
Mn Z
In the limit μ = μ0 ; ε = ε0 ; β = β0 ; Z = Z0 it is
2
M n | μ=μ0 = −j
ε=ε0 πβ0 a
and thus, just for the purpose of validation, we again get the exciting field of the line
source in free space.
If the cylinder is a perfect conductor, then the tangential component of the electric
field vanishes at = a and this leads to
∞
i0
(2)
1 Hn (β0 c) (2)
K = j β0 − Jn (β0 a) (2) Hn (β0 a)+
2 n=0 1 + δ0n Hn (β0 a)
+ Hn(2) (β0 c)Jn (β0 a) · cos(n(ϕ − ϕ0 ))
and with the recurrence relations and the Wronskian-determinants this expression sim-
plifies to
∞
i
(2)
1 Hn (β0 c)
K = − 0 cos(n(ϕ − ϕ0 )) .
πa n=0 1 + δ0n Hn(2) (β0 a)
Appendix
Vector Calculus
B
A( × C)
A
= C( × B)
= B(
C × A)
× (B
A × C)
A
= B( C)
− C(
A
B)
× B)
(A (C × D)
C)
= (A (B D)
− (A D)
(B C)
Metric factors: h1 , h2 , h3
V = V (u1 , u2 , u3 ) ;
A = e1 A1 (u1 , u2 , u3 ) + e2 A2 (u1 , u2 , u3 ) + e3 A3 (u1 , u2 , u3 )
1 ∂V 1 ∂V 1 ∂V
grad V = e1 + e2 + e3
h1 ∂u1 h2 ∂u2 h3 ∂u3
1 ∂ ∂ ∂
div A = (h2 h3 A1 ) + (h3 h1 A2 ) + (h1 h2 A3 )
h1 h2 h3 ∂u1 ∂u2 ∂u3
e /(h h ) e /(h h ) e /(h h )
1 2 3 2 1 3 3 1 2
= ∂/∂u1
rot A ∂/∂u2 ∂/∂u3
h1 A1 h2 A2 h3 A3
div grad V = ΔV =
1 ∂ h2 h3 ∂V ∂ h1 h3 ∂V ∂ h1 h2 ∂V
= + +
h1 h2 h3 ∂u1 h1 ∂u1 ∂u2 h2 ∂u2 ∂u3 h3 ∂u3
Cartesian coordinates:
u1 = x ; u2 = y ; u3 = z ; h1 = h2 = h3 = 1
∂V ∂V ∂V
grad V = ex + ey + ez
∂x ∂y ∂z
∂ 2V ∂2V ∂2V
ΔV = 2
+ 2
+
∂x ∂y ∂z 2
Circular-cylinder coordinates:
u1 = ; u2 = ϕ ; u3 = z ; h1 = 1 ; h2 = ; h3 = 1
∂V eϕ ∂V ∂V
grad V = e + + ez
∂ ∂ϕ ∂z
= 1 ∂ ( A ) + 1 ∂Aϕ + ∂Az
div A
∂ ∂ϕ ∂z
= e 1 ∂Az ∂Aϕ ∂A ∂Az ez ∂ ∂A
rot A − + eϕ − + ( Aϕ ) −
∂ϕ ∂z ∂z ∂ ∂ ∂ϕ
Spherical coordinates:
u1 = r ; u2 = ϑ ; u3 = ϕ ; h1 = 1 ; h2 = r ; h3 = r sin ϑ
r = er (ϑ, ϕ)r ; dr = er dr + eϑ rdϑ + eϕ r sin ϑdϕ
dv = r2 sin ϑ dr dϑ dϕ
∂V eϑ ∂V eϕ ∂V
grad V = er + +
∂r r ∂ϑ r sin ϑ ∂ϕ
= 1 ∂ 2 1 ∂ 1 ∂Aϕ
div A (r Ar ) + (Aϑ sin ϑ) +
r2 ∂r r sin ϑ ∂ϑ r sin ϑ ∂ϕ
= er ∂ ∂Aϑ eϑ 1 ∂Ar ∂
rot A (Aϕ sin ϑ) − + − (rAϕ ) +
r sin ϑ ∂ϑ ∂ϕ r sin ϑ ∂ϕ ∂r
eϕ ∂ ∂Ar
+ (rAϑ ) −
r ∂r ∂ϑ
∂2V 2 ∂V 1 ∂ ∂V 1 ∂2V
ΔV = + + 2 sin ϑ + 2 2
∂r2 r ∂r r sin ϑ ∂ϑ ∂ϑ r sin ϑ ∂ϕ2
dU r
grad U (r) =
dr r
· grad ) r
(A
= A
= U div A
div (U A) +A grad U
div (A × B)
rot A
= B −A rot B
dU
div [r U (r)] = 3 U (r) + r
dr
div grad U = Δ U
div rot A = 0
div (e × r) = 0 ; e constant vector
Appendix 263
rot (U A) + grad U × A
= U rot A
× B)
rot (A grad ) A
= (B − (A
grad ) B
+A
div B
−B
div A
rot (r U (r)) = 0
rot (e × r) = 2e ; e constant vector
rot grad U = 0
rot rot A − ΔA
= grad div A
rot rot rot (e U ) = e × grad (ΔU )
rot rot rot (r U ) = r × grad (ΔU )
Theorems
Stokes theorem
da
rot A = ds
A
a C
da × grad U = U ds
a C
Gauss theorem
dv
div A = da
A
v a
grad U dv = U da
v a
dv
rot A =
da × A
v a
2 ∂U
(grad U ) dv = U da for V = U and Δ U = 0
∂n
v a
264 Appendix
Stratton’s theorems
rot B
rot A −A
rot rot B
dv = × rot B
A da
v a
A rot rot B
−B
rot rot A
dv = × rot A
B −A
× rot B
da
v a
Bibliography
[1] Spiegel, Murray R.: “Theory and Problems of Vector Analysis”, Schaum’s Outline
Series, Mc Graw Hill Book Company, 1974
[3] Nascar, Syed A.: “2000 Solved Problems in Electromagnetics”, Schaum’s Solved
Problems Series, Mc Graw Hill Book Company, 1992
[4] Jackson, John D.: “Classical Electrodynamics”, John Wiley & Sons, 1998
[5] Jones, Douglas S.: “The Theory of Electromagnetism”, Pergamon Press, 1964
[6] Smythe, William R.: “Static and dynamic electricity”, Mc Graw Hill Book Com-
pany, 1968
[7] Stratton, Julius A.: “Electromagnetic Theory”, John Wiley & Sons, 2007
[9] Lehner, Günther: “Electromagnetic Field Theory for Engineers and Physicists”,
Springer Verlag, Berlin Heidelberg New York, 2010
[10] Morse, Philip, M., Feshbach, Herman: “Methods of Theoretical Physics, Part I,
Part II”, Mc Graw-Hill Book Company, Boston, 1953
[12] Moon, P., Spencer: “Field Theory Handbook”, Springer Verlag, Berlin Heidelberg
New York, 1971
of the magnetic flux density, 108 exciting, 25, 33, 40, 59, 61, 95, 107,
Normal derivative 108, 158
of the potential, 81, 84, 86, 89 given, 37, 51, 73, 79
given on a cylinder, 41, 56
Orthogonal series expansion, see Expan- given on a plane in free space, 72
sion given on a sphere, 70
of hyperbolic functions, 193
of a hemispherical charge distribu-
Orthogonality relation
tion, 44
for Bessel functions, 113, 173
of a space charge, 27
for Legendre polynomials, 71, 159
on a cylindrical surface, 89
for spherical Bessel functions, 162
on an axis, 44
for spherical harmonics, 134
resulting, 33, 35, 41, 61, 96, 97, 107,
for trigonometric functions, 57, 64,
110, 117
82, 90, 168, 190, 253
scalar magnetic, 5, 106, 115
Parallel-plate capacitor, 48, 78 Potential coefficients, 4, 14
Parallel-plate waveguide, 200, 209 Potential jump, 68, 74
Parallel-wire line, 106 Potential mirroring, 34
Partial waves, 194 Power
Permeability, 1 complex, 7
finite, 122 Power loss, 83, 86, 87, 157, 190
high, 117, 121 time-averaged, 8, 142, 144, 152, 157,
Permittivity, 1 160, 166
Phase constant, 11, 205, 209, 225, 231 Poynting theorem, 7
Phase vector, 247 Poynting vector, 9, 163, 168, 217
Phase velocity, 8, 11, 194 complex, 8, 142
Plane Primary current, 187
complex, 78
conducting, 20, 219 Radiation pattern, 9, 216
Point charge, 3, 29, 53, 76 Recurrence relation
Polarization for Bessel functions, 258
orthogonal, 249 for spherical Bessel functions, 135,
Potential, 2 162
circular symmetry, 12 for spherical harmonics, 138
complex electric, 4, 79, 91 Reflection, see Wave
complex magnetic, 6 Resistance, 86, 88, 90, 99, 143
electrodynamic, 7 frequency-dependent, see Impedance
electrostatic, 2, 12 Resonance frequency, 206
Index 271
Wave
electromagnetic, 8
excited, 197, 204
incident, 205, 249
plane, 247
reflected, 204, 226
transverse electric, 200, 209, 239
transverse magnetic, 200, 209, 239
Wave impedance, 8, 217
Waveguide
cylindrical, 10
dielectric slab, 231, 239
excitation, 225
rectangular, 196
Waveguide resonator
cylindrical, 206
Wavelength, 8
Weight function, 60