Electric Field for a cylinder without using Gauss's law?

[gsingh2011]

I recently learned how to calculate the field from a cylinder (inside and outside the cylinder) using Gauss's law. I was wondering how I would be able to derive the same formula without using Gauss's law (just for practice). My idea is that you would need to integrate the electric field from cylindrical shells of a radius 0 to R (where R is the radius of the cylinder), and to find the electric field of a cylindrical shell, you would need to integrate the field from rings along the length of the cylinder. The field from a ring is easy to calculate along its axis, but how would you calculate it in a less symmetric situation? Or is there a better way of doing this?

 

[jtbell]

to find the electric field of a cylindrical shell, you would need to integrate the field from rings along the length of the cylinder.

How about integrating the fields from a collection of line charges arranged as a cylindrical shell?

I haven't tried either approach myself, but I do know that the field of a ring at off-axis points is difficult to solve exactly. A line charge is easy, and its field is radial (no longitudinal or azimuthal components) which makes it easier (than for a ring) to integrate a bunch of them together.

 

[vanhees71]

The direct solution with the electrostatic Maxwell equations is usually simpler than using the integral form or the solution with the Green's function.

In your case you just write down the Poisson equation (Heaviside-Lorentz units),

\Delta \Phi=-\rho

in cylinder coordinates, assuming cylinder symmetry, i.e., you assume that \Phi and \rho are both functions of r alone. Then you have an ordinary differential equation to solve,

\frac{1}{r}[r \Phi'(r)]'=-\rho(r).

For you problem of a cylindrical shell, you have

\rho(r)=\frac{\lambda}{2 \pi r} \delta(r-a),

where a is the radius of the cylinder and \lambda the charge per unit length on the surface.

Now except for r=a your charge density is 0. So we first solve the homogeneous equation

[r \Phi'(r)]'=0 \; \Rightarrow \; r \Phi'(r)=A \; \Rightarrow \; \Phi(r)=A \ln(r/a)+B.

Here, I've written the integration constant such that the argument of the logarithm becomes dimensionless as it must be.

The above solution is valid for r<a and r>a, but the integration constants have to be chosen differently in the two regions such as to solve the equation including the charge density. For r<a we must have A=0, and then we can choose also B=0. So we have

\Phi(r)=0 \quad \text{for} \quad r<a.

To find the solution for \quad r>a we must make sure that the potential is continuous at r=a, leading to B=0. To find A, we integrate the equation,

[r \Phi'(r)]'=\frac{\lambda}{2 \pi} \delta(r-a)

over an infinitesimally small interval around r=a. Using the solution for r<a this gives

A=\frac{\lambda}{2 \pi},

i.e., the solution reads

\Phi(r)=\frac{\lambda}{2 \pi} \ln(r/a) \quad \text{for} \quad r>a.

The electric field is

\vec{E}=\vec{e}_r \frac{\lambda}{2 \pi r} \Theta(r-a),

where \Theta is the Heaviside-unitstep function.

 

[Delta2]

I think OP is interested in the case where the charge density is not zero inside the cylinder.

 

[vanhees71]

Ah, I've misunderstood the question since he wrote about cylindrical shells.

For a given cylinder-symmetric charge distribution it's also easily evaluated in the way of my previous posting. Let's take the most simple case of a homogeneous charge density inside the cylinder. In cylindrical coordinate it reads

\rho(r)=\rho_0 \Theta(a-r).

Then you can simply integrate:

\frac{1}{r} [r \Phi'(r)]'=-\rho(r).

For r<a we have

[r \Phi'(r)]'=-r \rho_0, \; \Rightarrow \; \Phi'(r)=-\rho_0 \frac{r}{2}+\frac{A}{r} \; \Rightarrow \; \Phi(r)=-\rho_0 \frac{r^2}{4} + A \ln(r/a)+B.

Since the potential should not be singular at r=0 we must have A=0.

For r>a we can use the previous solution

\Phi(r)=A' \ln(r/a)+B', \quad \Phi'(r)=\frac{A'}{r}.

To find the integration constants, we have to fulfill the boundary conditions that \Phi and \Phi' are continuous at r=a. This leads to

-\rho_0 \frac{a^2}{4}+B=B', \quad -\rho_0 \frac{a}{2}=\frac{A'}{a}.

Setting

B=\rho_0 \frac{a^2}{4}

leads to B'=0, and we finally have

\Phi(r)=\frac{\rho_0}{4}(a^2-r^2) \Theta(a-r)-\frac{\rho_0 a^2}{2} \ln(r/a) \Theta(r-a).

The electric field reads

\vec{E}=-\vec{\nabla} \Phi=\vec{e}_r \frac{\rho_0}{2} \left[r \Theta(a-r) + \frac{a^2}{r} \Theta(r-a) \right ].