Laplace equation

Laplace Equations

1. Special Techniques for Calculating Potential
3.1 Laplace’s Equation
3.2 The Method of Images
3.3 Separation of Variables
3.4 Multipole Expansion

2. 3.1 Laplace’s Equation
3.1.1 Introduction
3.1.2 Laplace’s Equation in One Dimension
3.1.3 Laplace’s Equation in Two Dimensions
3.1.4 Laplace’s Equation in Three Dimensions
3.1.5 Boundary Conditions and Uniqueness
Theorems
3.1.6 Conducts and the Second Uniqueness
Theorem

3. 3.1.1 Introduction
The primary task of electrostatics is to study the interaction
(force) of a given stationary charges.
since
this integrals can be difficult (unless there is symmetry)

we usually calculate
This integral is often too tough to handle analytically.
2
0
ˆ
1
( )
4
F q E
test
R
E d
R
r t
pe
=
= ò
v v
v
ˆˆ ˆ
x y z
E E x E y E z= + + +´´´
v
0
1 1
( )
4
V d
R
E V
r t
pe
= ò
=-Ñ
v
\
Q

4. In differential form
Eq.(2.21):Poisson’s eq.
Laplace’s eq
or
The solutions of Laplace’s eq are called harmonic function.
3.1.1
2
0
V
r
e
Ñ =-
2
0VÑ =
2
2 2 2
2 2 2
1 1
0 ( )
0
i
k i i i
k
T
V A
h q h q
V V V
x y z
¶ ¶
=Ñ =
¶ ¶Õ
¶ ¶ ¶
+ + =
¶ ¶ ¶
 to solve a differential eq. we need boundary conditions.
 In case of ρ = 0, Poisson’s eq. reduces to

5. 3.1.2 Laplace’s Equation in One Dimension
2
2
0
d V
V mx b
dx
= Þ = +
m, b are to be determined by B.C.s
e.q.,
1. V(x) is the average of V(x + R) and V(x - R), for any R:
2. Laplace’s equation tolerates no local maxima or minima.
( 1) 4 1
5
( 5) 0 5
V x m
V x
V x b
= = =-ì ì
Þ Þ =- +í í
= = =î î
V=4 V=0
1 5
x
[ ]
1
2
( ) ( ) ( )V x V x R V x R= + + -

6. Method of relaxation:
3.1.2 (2)
 A numerical method to solve Laplace equation.
 Starting V at the boundary and guess V on a grid of interior points.
Reassign each point with the average of its nearest neighbors.
Repeat this process till they converge.

7. 3.1.2 (3)
The example of relaxation
nV(x=0)V(x=1)V(x=2)V(x=3)V(x=4)
0 4 0 0 0 0
1 4 2 0 0 0
2 4 2 1 0 0
3 4 1 0
4 4 0
5 4 0
6 4 0
7 4 0
8 4 0
9 4 0
104 0
4 3 2 1 0

8. 3.1.3 Laplace’s Equation in Two Dimensions
A partial differential eq. :

There is no general solution. The solution will be given in 3.3.
We discuss certain general properties for now.
1.The value of V at a point (x, y) is the average of those around
the point.
2. V has no local maxima or minima; all
extreme occur at the boundaries.
3. The method of relaxation can be applied.
2 2
2 2
0
V V
x y
¶ ¶
+ =
¶ ¶
1
( , )
2 circle
V x y Vdl
Rp
=òÑ

9. 3.1.4 Laplace’s Equation in Three Dimensions
1. The value of V at point P is the average value of V
over a spherical surface of radius R centered at P:
The same for a collection of q by the superposition principle.
2. As a consequence, V can have no local maxima or minima;
the extreme values of V must occur at the boundaries.
3. The method of relaxation can be applied.
2
1
( )
4
sphere
V p Vda
Rp
=òÑ
2 2 2
0
1
, 2 cos
4
q
V r R rR q
pe
= = + -
r
r
1
2 2 2
2
2
0
1
[ 2 cos ] sin
4
4
ave
q
V r R rR R d d
R
q q q j
pe
p
-
= + - ò
2 2
0
1
2 cos
04 2
q
r R rR
rR
p
q
pe
= + -
( )( )
at the center of the sphere
0 0
1 1
4 2 4
q q
r R r R V
rR rpe pe
é ù= + - - = =
ë û

10. 3.1.5 Boundary Conditions and Uniqueness Theorems
in 1D,
·
one end
V
a
·
the other end
V
b
V is uniquely determined by its value at the boundary.
First uniqueness theorem :
the solution to Laplace’s equation in some region is uniquely
determined, if the value of V is specified on all their surfaces;
the outer boundary could be at infinity, where V is ordinarily taken
to be zero.
V mx b= +

11. 3.1.5
Proof: Suppose V
1
, V
2
are solutions
at boundary at boundary.
V
3
= 0 everywhere
hence V
1
= V
2
everywhere
2
1
0VÑ =
2
2
0VÑ =
1 2
?V V=
3 1 2
V V V= -
2 2 2
3 1 2
0V V VÑ =Ñ -Ñ =
1 2 3
0V V V= Þ =
\

12. Proof.
at boundary.
Corollary : The potential in some region is uniquely determined if
(a) the charge density throughout the region, and
(b) the value of V on all boundaries, are specified.
3.1.5
The first uniqueness theorem applies to regions with charge.
2
1
0
V
r
e
Ñ =-
2
2
0
V
r
e
Ñ =-
3 1 2
V V V= -
2 2 2
3 1 2
0 0
0V V V
r r
e e
Ñ =Ñ -Ñ =- + =
3 1 2
3
0
0, . .,
1 2
V V V
V ie V V
= - =
\ = =

13. 3.1.6 Conductors and the Second Uniqueness Theorem
 Second uniqueness
theorem:
Proof:
Suppose both and are satisfied.
and
0
1
1
E
e
rѴ =
v
0
1
2
E
e
rѴ =
v
In a region containing conductors and filled
with a specified charge density ρ,
the electric field is uniquely determined if the
total charge on each conductor is given.
(The region as a whole can be bounded by
another conductor, or else unbounded.)
2
E
v
1
E
v
0 0
1 1
1 1 2 2
,
ith conducting ith conducting
surface surface
i i
E da Q E da Q
e e
´ = ´ =ò ò
v vv v
Ñ Ñ

0 0
1 1
1 2
,
outer bourndary outer bourndary
tot tot
E da Q E da Q
e e
´ = ´ =ò ò
v vv v
Ñ Ñ

14. for V
3
is a constant over each conducting surface
define
3.1.6
3 1 2
E E E= -
v v v
0 0
3 1 2 1 2
( ) 0E E E E E
r r
e e
Ѵ =Ѵ - =Ѵ -Ѵ = - =
rv v v v
3 1 2
0E da E da E da´ = ´ - ´ =ò ò ò
v v vv v v
Ñ Ñ Ñ
2
3 3 3 3 3 3 3
( ) ( ) ( )V E V E E V EѴ = Ѵ + Ñ =-
v v v
2
3 3 3 3 3
( )
volume surface
E d V E d V E dat t=- Ѵ =- ´ò ò ò
v v v
Ñ
3 3
0
surface
V E da=- ´ =ò
vv
Ñ
3 1 2
0 E E E\ = =
v v v
i. e.,

15. 3.2 The Method of Images
3.2.1 The Classical Image Problem
3.2.2 The Induced Surface Charge
3.2.3 Force and Energy
3.2.4 Other Image Problems

16. 3.2.1 The Classical Image Problem
What is V (z>0) ?
( )
2 2 2 2
1. 0 0
2. 0
V z
V for x y z d
= =
® + + >>
B.C.
The first uniqueness theorem guarantees that
there is only one solution.
If we can get one any means, that is the only answer.

17. Trick :
Only care
z > 0
z < 0
is not of
concern
3.2.1
in original
problem
for
V(z=0) = const = 0
( )0 0E z< =
v
0z³
( ) ( )
2 22 2 2 2
0
1
( , , )
4
q q
V x y z
x y z d x y z d
pe
é ù
-
ê ú
= +
ê ú
+ + - + + +
ê úë û
ˆ
0
0
0
z
x y
E E z
E E
at z
or E
=-
= =æ ö
ç ¸
=è ø
=
P
v

18. 3.2.2 The Induced Surface Charge
Eq.(2.49)
0
ˆˆ
V
E z V z
z
s
e
¶
= =-Ñ =-
¶
v
0
0z
V
z
s e
=
¶
=-
¶
0 0
z z
A
E A E
s s
e e
´
- = =-
0
0
Q
E da atz
e
´ = @ò
vv
QÑ
( ) ( )
2 22 2 2 2
0
1
4
z
q q
z
x y z d x y z d
s
p
=
ì ü
- ¶ï ï
= -í ý
¶
ï ï+ + - + + +
î þ

19. total induced surface charge
3.2.2
Q=-q
3 3
2 2 2 2 2 2
2 2
0
1 1
2( ) 2( )
2 2
4
( ) ( )
z
z d z d
q
x y z d x y z d
p
=
ì ü
- ´ - - ´ +ï ï
-ï ï
= -í ý
ï ï
é ù é ù+ + - + + +
ï ïë û ë ûî þ
3/2 3/2
2 2 2 2 2
2 2
qd qd
x y d r d
s
p p
- -
= =
é ù é ù+ + +
ë û ë û
0
1/2
2 2
1
qd q
r d
¥-
=- =-
é ù+
ë û
2
3/20 0
2 2
2
qd
Q da rdrd
r d
p
s q
p
¥ -
= =
é ù+
ë û
ò ò ò

20. 3.2.3 Force and Energy
The charge q is attracted toward the plane.
The force of attraction is
With 2 point charges and no conducting plane, the energy is
Eq.(2.36)
( )
( ) ( )
2
2 2
0 0
1 1
ˆˆ
4 4
2
q q q
F z z
dd d
pe pe
-
= =-
é ù- -
ë û
v
()
() ( )
( )
( )
2
1
2
0 0
2
0
1
2
1 1 1

2 4 4
1

4 2
i i i
i
W qV p
q q
q q
d d
d d
q
d
pe pe
pe
=
=
ì ü é ù
ï ï é ù ê ú
= ´ - + - ´í ý
ê ú ê ú+ë û
ï ï - -
ê úë ûî þ
=-
å

21. 3.2.3 (2)
or
For point charge q and the conducting plane at z = 0 the energy
is half of the energy given at above, because the field exist
only at z ≥ 0 ,and is zero at z < 0 ; that is
( )
2
0
1
4 4
q
W
dpe
=-
()
( )
2
2
0
2 2
0 0
1
4 2
1 1

4 4 4 4
d d
d
q
W F dl dz
z
q q
z d
pe
pe pe
¥ ¥
¥
= ´ =
æ ö
= - =-ç ¸
ç ¸
è ø
ò ò
vv
ˆ( )dl dlz=-
v
Q

22. 3.2.4 Other Image Problems
Stationary distribution of charge
1
q
2
q
3
q
1
q-
2
q-
3
q-

23. 3.2.4 (2)
Conducting sphere of radius R
Image charge
at
R
q q
a
¢=-
2
R
b
a
=
( )
0
0
0
1
,
4
1

ˆˆˆˆ4
1

4
ˆˆˆˆ
q q
V r
q q
rr as rr bs
q q
a r
r r s b r s
r b
q
pe
pe
pe
¢æ ö
= +
ç ¸
è ø
æ ö ¢
= +ç ¸
- -
è ø
æ ö
ç ¸
¢
= +ç ¸
ç ¸
- -
ç ¸
è ø
¢rr
1
2 2
22 cosr a ra qé ù= + -
ë û
r
1
2 2
22 cosr b rb qé ù= + -
ë û¢r

24. Force
[Note : how about conducting circular cylinder?]
3.2.4 (3)
( )
0
1
, 0
4
ˆˆˆˆ
q q
V R
a R
Rr s b r s
R b
q
pe
æ ö
ç ¸
¢
= + =ç ¸
ç ¸
- -
ç ¸
è ø
q q
R b
¢
=-
b R
q q q
R a
¢=- =-
a R
R b
=
2
R
b
a
=
( ) ( )
2
2 2
2 20 0
1 1
4 4
qq q Ra
F
a b
a R
pe pe
¢
= =-
-
-

25. 3.3 Separation of Variables
3.3.0 Fourier series and Fourier transform
3.3.1 Cartesian Coordinate
3.3.2 Spherical Coordinate

26. ˆˆˆ, ,x y zare unique because
3.3.0 Fourier series and Fourier transform
Basic set of unit vectors in a certain coordinate can express

any vector uniquely in the space represented by the coordinate.
e.g.
Completeness: a set of function ()xf
n
() ()xfCxf
n
n
nå
¥
=
=
1
()xf
is complete.
if for any function
Orthogonal: a set of functions is orthogonal
if ()()0=ò
dxxfxf
m
b
a
n
mn¹
=const
for
for mn=
in 3D. Cartesian Coordinate. zVyVxViVV
zyx
N
i
i
ˆˆˆˆ
1
++==å
=

zyx
VVV,, are orthogonal.
jiˆˆ×=
i j=
i j¹0
1

27. 3.3.0 (2)
A complete and orthogonal set of functions forms
a basic set of functions.
e.g.
NnkÎ,
( ) ()sin sinkx kx- =-
( ) ()cos coskx kx- =
odd
{
nkif
nkif
dxnxkx
¹
=
-ò
=
0
coscos
p
p
p
0cossinò
-
=
p
p
dxnxkx
{
nkif
nkif
dxnxkx
¹
=
-ò
=
0
sinsin
p
p
p
even
sin(nx) is a basic set of functions for any odd function.
cos(nx) is a basic set of functions for any even function.
sin(nx) and cos(nx) are a basic set of functions for any functions.

28. ÷
ø
ö
ç
è
æ
for any()xf ()
() ()
2
xfxf
xg
--
= ()
() ()
2
xfxf
xh
-+
=odd
() () ()xhxgxf +=
even
3.3.0 (3)
,
­
odd
­
even
[sinkx][coskx]
Fourier transform and Fourier series
() ( )
0
(*) sin cos
n n
n
f x A nx B nx
¥
=
= +å
()()
1
sin
n
A f x nx dx
p
pp-
=ò
()()
1
cos
n
B f x nx dx
p
pp-
=ò
() 0
2
1
00
==ò
-
AdxxfB
p
pp
}
¥××××=,2,1n
;

29. 3.3.0 (4)
Proof
() ()
()
* sin
sin
kx dx
f x kxdx
p
p
p
p
-
-
´
Þ ´
ò
ò
( ) ()
0
sin cos sin
n n
n
A nx B nx kx dx
p
p
¥
-
=
= +åò
()()
0
sin sin
n
n
A nx kx dx
p
p
¥
-
=
=åò
() ()
0
1
sin
k
A f x kx dx
p
p
\ = ò
if 0
k
A kp= ¹

30. 0
cos cos
n
n
B nx kx dx
p
p
¥
-
=
=åò
() ()
() ()
* cos
cos
kx dx
f x kx dx
p
p
p
p
-
-
´
Þ
ò
ò
( ) ()
0
sin cos cos
n n
n
A nx B nx kx dx
p
p
¥
-
=
= +åò
()0
1
2
B f x dx
p
pp-
=ò
{=
() ()
1
cos
n
B f x nx dx
p
pp-
\ = ò
3.3.0 (5)
B
0
2p for k = 0
B
k
p for k = 1, 2, …

31. 3.3.1 Cartesian Coordinate
Use the method of separation of variables to solve the Laplace’s eq.
Example 3
0)(
)()0(
0)(
0)0(
0
®¥®
==
==
==
xV
yVxV
yV
yV
p
Find the potential inside this “slot”?
Laplace’s eq.
2 2
2 2
0
V V
x y
¶ ¶
+ =
¶ ¶
set )()(),( yYxXyxV =
2 2
2 2
0
X Y
Y X
x y
¶ ¶
+ =
¶ ¶

32. 0
11
)(
2
2
)(
2
2
=
¶
¶
+
¶
¶

onlyydepedentyg
onlyxdepedentxf
y
Y
Yx
X
X
f and g are constant
1 2
0C C+ =
12
2
1
)( C
dx
Xd
X
xf ==
2
2
2
1
( )
d Y
g y C
Ydy
= =
2
21 kCC =-=set
so Xk
dx
Xd
2
2
2
=
2
2
2
d Y
k Y
dy
=-
kyDkyCyYBeAexX
kxkx
cossin)(,)( +=+=
-
=),(yxV )(
kxkx
BeAe
-
+ )cossin( kyDkyC +
3.3.1 (2)
1 2
, areconstant
( ) ( ) 0
C C
f x g y+ =

33. 3.3.1 (3)
B.C. (iv)( ) 0 0, 0V x A k®¥ ® Þ = >
B.C. (i)( 0) 0 0V y D= = Þ =
=),(yxV
kx
e
-
)cossin( kyDkyC +
( , ) sin
kx
V x y Ce ky
-
=
B.C. (ii)( ) 0 sin 0 1,2,3V y k kp p= = Þ = = L
B.C. (iii) )()0(
0
yVxV ==
=)(
0
yV
1
sin
kx
k
k
C e ky
¥
-
=
å A fourier series for odd function
00
2
( )sin
k
C V y kydy
p
p
=ò
The principle of superposition =),(yxV
1
sin
kx
k
Ce ky
¥
-
=
å

34. 3.3.1 (4)
For antconstVyV ==
00
)(
0
0
0
4
0
2
sin
0
2
(1 cos )
k
V
k
V
C kydy
if k even
V
k
k if k odd
p
p
p
p
p
=
=ì
ï
= - = í
=
ï
î
ò
sin10 0
sinh
1,3,5,
4 2 1
( , ) sin tan ( )
yk
x
k
V V
V x y e ky
kp p
- -
=
= =å
L

35. 3.4 Multipole Expansion
3.4.1 Approximate Potentials at Large distances
3.4.2 The Monopole and Dipole Terms
3.4.3 Origin of Coordinates
in Multipole Expansions
3.4.4 The Electric Field of a Dipole

36. 1 1
2
1 cos
24
0
( ) cos
( )
d
r
qd
r
V p
q
pe
q- @
@
r r
+ -
3.4.1 Approximate Potentials at Large distances
Example 3.10 A dipole (Fig 3.27). Find the approximate potential
at points far from the dipole.
Þ
1
4
0
( ) ( )
q q
V p
pe
+ -
= -
r r
2 2 2
2
( ) cos
d
r rd q
±
= + mr
2
2
2
4
(1 cos )
d d
r
r
r q= ± +
2
(1 cos )
r d
d
r
r q
>>
@ m
1
1 1 1
2
2
(1 cos ) (1 cos )
d d
r r r r
q q
-
±
@ @ ±m
r

37. 3.4.1
Example 3.10 For an arbitrary localized charge distribution.
Find a systematic expansion of the potential?
0
1 1
( )
4
V p d r t
pe
=ò
r
2 2 2
2 2
2 cos
[1 ( ) 2( )cos ]
r r
r r
r r rr
r
q
q
¢ ¢
¢ ¢= + -
= + -
r
1 ( )( 2cos )
r r
r r
r q
¢ ¢
= +Î Î= -r

38. 3.4.1(2)
1 1 3 5 2 2 3 3
2 8 16
[1 ( )( 2cos ) ( ) ( 2cos ) ( ) ( 2cos ) ]
r r r r r r
r r r r r r r
q q q
¢ ¢ ¢ ¢ ¢ ¢
= - - + - - - + L
1 3 1 5 3 2 2 3 3
2 2 2 2
[1 ( )cos ( ) ( cos ) ( ) ( cos cos ) ]
r r r
r r r r
q q q q
¢ ¢ ¢
= + + - + - + L
1
0
( ) (cos )
rn
n
r r
n
P dq r t
¥
¢
=
¢=å
Monopole term Dipole term Quadrupole term
Multipole expansion
1
1 1 1 1 3 5 2 32
2 8 16
(1 ) (1 )
r r
-
Î Î Î= +Î = - + - + L
r
1
( 1)4
0
0
1 1 1 1 3 1 2
2 34 2 2
0
1
( ) ( ) (cos )
[ cos ( ) ( cos ) ]
n
n
n
n
r
r r
V r r P d
r
d r d r d
pe
pe
q r t
r t qr t q r t
¥
+
=
¢=
¢ ¢= + + - +
å ò
ò ò ò
L

39. 3.4.2 The Monopole and Dipole Terms
monopole
dipole
dominates if r >> 1
dipole moment (vector)
A physical dipole is consist of a pair of equal and
opposite charge,
1
4
0
( )
Q
mon
r
V p
pe
=
1 1
24
0
1 1
24
0
( ) cos
ˆ
dip
r
r
P
V p r d
r r d
pe
pe
qr t
r t
¢=
¢= ´
ò
ò
v
v
14243
ˆcosr r rq¢ ¢= ´
v
ˆ1
24
0
1
( )
p r
dip
r
n
i i
i
V p
p r d qr
pe
r t
´
=
=
¢ ¢= = åò
v
v v v
q±
( )p qr qr q r r qd
+ - + -
¢ ¢ ¢ ¢= - = - =
v
v v v v v

40. 3.4.3 Origin of Coordinates in Multipole Expansions
if
( )
p r d
r d d
r d d d
p dQ
r t
r t
r t r t
¢=
¢= -
¢= -
= -
ò
ò
ò ò
vv
v
v
v
v
0Q p p= =
v

41. 3.4.4 The Electric Field of a Dipole
A pure dipole
ˆ cos
2 2
4 4
0 0
2 cos
3
4
0
1 sin
3
4
0
1
sin
3
4
0
( , )
0
ˆ
ˆ( , ) (2cos sin )
P r P
dip
r r
V P
r
r
r
V P
r
r
V
r
P
dip
r
V r
E
E
E
E r r
q
pe pe
q
pe
q
q
q
pe
j
q j
pe
q
q q qq
´
¶
-
¶
¶
-
¶
¶
-
¶
= =
= =
= =
= =
= +
v
v