**Fractional
order Legendre functions and electric fields (why are lightning-rods
pointy)?**

This is a bit of mathematical amusement which came out of a study of field concentration and wave scattering by sharp edges and points. It is well known that electric fields are strongly concentrated in the vicinity of sharp points. In fact, near perfectly conducting points, the electric field will approach infinity. What is the nature of this singularity?

To answer this question, let us assume that we can model the pointy structure as an infinite circular cone with the point at the origin (0,0,0) with angle alpha.

In this situation, we shall solve the Laplace Equation for the scalar electric potential V

subjected to the boundary conditions

At this point, we shall focus on azimuthally symmetric solutions (where V does not depend on phi, the azimuthal angle). This allows us to write Laplace's equation in spherical coordinates as

This partial differential equation (PDE) is separable by assuming

where
R(r) and P() are functions of
their respective variables *only!*
This allows us to write the PDE as two ordinary differential
equations

where *l(l+1)* is the
so-called separation constant.

It is easy to show that the differential equation in r has the solutions

where *A _{l}*

We can simplify the appearance
of the differential equation for *P* by making the substitution

The second DE now looks like

This is the famous Legendre
differential equation and it is an *eigenvalue* equation in *
*,
where *l *is the eigenvalue. If there are no boundary
conditions, but we constrain the solution to regular (non-singular)
solutions in
,
we get the familiar Legendre functions of integer order (the Legendre
polynomials), viz.

The eigenvalues *l* in this
case are 0, 1, 2, 3, 4, .....

In the case of the pointy
circular cone, however, things are not quite so simple. The
eigenvalues do not necessarily take integer values and the polynomial
series is not necessarily finite. Perhaps the most straightforward
way to see what these functions look like is to use some type of
numerical method to generate *P* given the boundary condition

We constructed a simple solution method based on a 1-dimensional finite element method to generate the solutions to the eigenvalue problem. By assuming the cone angle alpha to be 180 degrees, we get the flat ground plane problem (in effect, the cone becomes a perfect electric conducting plane). We expect the solutions to be all the odd Legendre polynomials of integer order. In fact, using 100 cubic Hermite polynomials to approximate the Legendre function we get the eigenvalues

Order |
Computed Eigenvalue |
---|---|

1 |
1.0000000 |

2 |
3.0000000 |

3 |
5.0000000 |

4 |
7.0000000 |

5 |
9.0000000 |

6 |
11.000000 |

7 |
13.000000 |

8 |
15.000000 |

9 |
17.000001 |

10 |
19.000002 |

We see that the expected eigenvalues are computed to 7 decimal places quite easily. The first five functions look like

By comparing with the analytic
expressions with the computed solutions, we are satisfied the
numerical method is working properly. We see that there are no
singular fields in this case, because the eigenvalues *l* are
all greater than 1.

Now, let's move on to something
more interesting. Let us assume that we have a cone whose angle
is
150 degrees (i.e. the cone point spans an angle of 60 degrees). The
first 10 eigenvalues *l *are given by

Eigenvalue Order |
Computed Eigenvalue |
---|---|

1 |
0.3461839 |

2 |
1.568297 |

3 |
2.777734 |

4 |
3.982932 |

5 |
5.186204 |

6 |
6.388442 |

7 |
7.590066 |

8 |
8.791295 |

9 |
9.992257 |

10 |
11.19303 |

If we consider the first eigenvalue, we see that the solution for the potential near the tip of the cone can be written as

where
refers
to the *fractional* order
Legendre function which describes the angular behaviour of the
potential near the cone tip. It looks like

By recognising that the electric field is or

Numerically speaking, we have

We
can see that the electric field goes to infinity as we get closer to
the cone tip (*r*=0) .

Now,
we'll try a very pointy cone, where alpha is 350 degrees. The
computed eigenvalue that gives rise to singular fields is *l *=
0.1581469. The fractional-order Legendre function looks like P_0.158
in the plot below.

It
is interesting to see that for the sharper cone point, the slope of
the Legendre function drops off much more rapidly. This makes sense
physically, because the this cone is thinner at any given* r*
than the previous one. The fields will vary much more rapidly. In
fact, we have

*
*

The
singularity is significantly stronger in this case. The exponent is
-0.8418531 vesus -0.6538161 in the previous “less pointy”
cone case... Hence, if we want the enhance the possibility of
causing a discharge (as in the case of the lightning rod), we want a
*high *electric field. This means, the pointier we make the
lightningrod, the higher the field and the more likely we can cause
charge to leak away, avoiding a possibly damaging lightning strike!

**References:**

S. Ramo, J. R. Whinnery and T. D. Van Duzer,

*Fields and Waves in Communication Electronics,*Wiley, 1984.J. Van Bladel,

*Singular Electromagnetic Fields and Sources,*IEEE Press, 1995.