10

Sphere or Cylinder in Uniform External Field

ABSTRACT If a conducting or dielectric object is present within the electric field caused by external sources, then the induced charges on such an object adds to the existing external field. If the size of the conducting or dielectric object is small compared to the external field region, then the electric field is modified only in the vicinity of such an object. In order to get the complete solution of the electric field in such configuration, the boundary conditions imposed by the conducting or dielectric object must also be satisfied. If the object is of well-defined geometric shape such as cylinder or sphere, and if the external field is considered to be uniform in space, then the complete solution for the electric field could be obtained analytically with the help of method of separation of variables. These analytical solutions are often used to validate numerical methodology and could also be used to estimate field modification due to small conducting or dielectric pieces present in large high-voltage insulation arrangement.

10.1 Introduction

Conducting and dielectric components are integral parts of any electrical equipment. If the size of the conducting or dielectric object is very small compared to dimensions of the field region where the object is located, then the object contributes to the field only in the domain near the object. In many cases, such objects are present as stray bodies in high-voltage insulation arrangement. As practical examples, one may cite a small piece of conductor or dielectric floating in liquid insulation of large volume in transformers, metallic dust particles floating in gaseous insulation within gas-insulated system (GIS) and so on. It is important to understand how the presence of a conducting or dielectric object modifies the external field in the vicinity of the object, because any enhancement of the electric field intensity due to the conducting or dielectric object may lead to unwanted discharge or in the worst case failure of the insulation system.

If it is assumed that the source charges (in practical arrangement, the electrodes or conductors with specific potentials) that produce the external field is located far away from the object under consideration, then they are unaffected by the presence of the object. Consequently, the field due to the source charges may be considered to be uniform at the location of the object. If the object is such that its shape is defined by well-known mathematical functions, for example, cylinders or spheres, then the complete solution for the electric field due to the source charges located at far away positions and the induced charges on the surface of the object could be obtained by solving Laplace’s equation considering the field region to be free from any volume charge. However, in order to get the complete solution appropriate boundary conditions on the surface of the object, whether it is conducting or dielectric, need to be satisfied. One of the common methods of getting the analytical solution for a cylinder or sphere in a uniform external field is the method of separation of variables as described in this chapter.

10.2 Sphere in Uniform External Field

Consider a spherical object of radius a within uniform external field, as shown in Figure 10.1. Because the boundary is a sphere of r = constant, hence the system is best described in spherical coordinates, as shown in Figure 10.1. The uniform external field is given by E0=E0u^z and the potential at any point due to the external field is given by E0r cosθ = E0z with respect to the centre of the sphere. In order to get the complete solution for the electric field in this system, Laplace’s equation in spherical coordinates, as given in Equation 10.1, needs to be solved.

Images

FIGURE 10.1
Sphere in uniform external field.

1r2r(r2Vr)+1r2sinθθ(sinθVθ)+lr2sin2θ2Vϕ2=0(10.1)

The field system has azimuthal symmetry with respect to the z-axis, that is, the field system does not change with the rotation around the z-axis. Therefore, z-axis is made the polar axis in the spherical coordinate system. Then the field is independent of coordinate φ and Laplace’s equation reduces to

1r2r(r2Vr)+1r2sinθθ(sinθVθ)=0(10.2)

In order to separate the terms of the left-hand side (LHS) of Equation 10.2 into functions of only one variable, Equation 10.2 may be rewritten by multiplying r2 as

r(r2Vr)+1sinθθ(sinθVθ)=0(10.3)

Then LHS of Equation 10.3 is the sum of two terms, which are functions of only one variable each, that is, the first term is function of r only, whereas the second term is function of θ only. The solution to Equation 10.3 can be obtained as the product of two functions of which one is dependent only on r and the other is dependent only on θ.

Let the assumed solution be

V(r,θ)=M(r)N(θ)(10.4)

The assumed solution is convenient as the boundary lies at r = constant.

Combining Equations 10.3 and 10.4

r[r2M(r)N(θ)r]+1sinθθ[sinθM(r)N(θ)θ]=0or,N(θ)r[r2M(r)r]+M(r)1sinθθ[sinθN(θ)θ]=0

Dividing by M(r)N(θ),

1M(r)ddr[r2dM(r)dr]+1N(θ)1sinθddθ[sinθdN(θ)dθ]=0(10.5)

The partial derivatives become total derivatives in Equation 10.5, as each term is dependent on only one coordinate.

The sum of two terms of the LHS of Equation 10.5 could be zero only when the two terms are separately equal to opposite and equal constant terms as given in Equation 10.6.

Equal and opposite separation constant solution:

1M(r)ddr[r2dM(r)dr]=+p and 1N(θ)1sinθddθ[sinθdN(θ)dθ]=-p(10.6)

where:

p is a positive constant

Another solution is obtained when the separation constant is zero.

Hence, the zero separation constant solution

1M(r)ddr[r2dM(r)dr]=0 and 1N(θ)1sinθddθ[sinθdN(θ)dθ]=0(10.7)

Each of the above-mentioned two solutions is to be obtained separately.

Determination of the zero separation constant solution:

The first term of Equation 10.7 is

1M(r)ddr[r2dM(r)dr]=0

where:

M(r) is non-zero

Therefore,

ddr[r2dM(r)dr]=0,or,r2dM(r)dr=C0,or,dM(r)dr=C0r2

Integrating and incorporating constants of integration

M(r)=C10r+C20(10.8)

Next the second term of Equation 10.7 is

1N(θ)1sinθddθ[sinθdN(θ)dθ]=0

where:

N(θ) is non-zero

Therefore,

1sinθddθ[sinθdN(θ)dθ]=0,or,ddθ[sinθdN(θ)dθ]=0or,sinθdN(θ)dθ=A0

Integrating and incorporating constant of integration

N(θ)=A10ln(tanθ2)+A20(10.9)

Equation 10.9 becomes undefined for θ = π. But this is not feasible in the given system as potential must be a continuous function. Therefore, A10 should be zero in Equation 10.9. Therefore,

N(θ)=A20(10.10)

Then from Equations 10.4, 10.8 and 10.10, the zero separation constant solution can be obtained as

V(r,θ)=C1r+C2(10.11)

where:

C1 = A20C10 and C2 = A20C20

Determination of the equal and opposite separation constant solution:

The first term of Equation 10.6 is

1M(r)ddr[r2dM(r)dr]=+por,ddr[r2dM(r)dr]=+pM(r)(10.12)

Putting M(r) = Crn in Equation 10.12, we get

ddr(r2Cnrn-1)=+pCrn,or,Cn(n+1)rn=pCrn,or,n2+n-p=0

Hence,

n=12(-1±1+4p)(10.13)

The second term of Equation 10.6 is

1N(θ)1sinθddθ[sinθdN(θ)dθ]=-p1sinθddθ[sinθdN(θ)dθ]=-pN(θ)(10.14)

Putting N(θ) = Bcosθ in Equation 10.14,

ddθ(-Bsin2θ)=-pBcosθsinθ,or,p=2.

Hence, from Equation 10.13 n = +1, −2.

Therefore,M(r)=Cr+C"r-2 and N(θ)=Bcosθ(10.15)

From Equations 10.4 and 10.15,

V(r,θ)=(C3r+C4r2)cosθ(10.16)

where:

C3 = CB and C4 = CB

The complete solution for potential function is uniquely given as a linear combination of the two solutions given by Equations 10.11 and 10.16.

V(r,θ)=C1r+C2+(C3r+C4r2)cosθ(10.17)

where:

The constants are determined by satisfying the boundary conditions

It is evident from Equation 10.17 that the first term corresponds to a net charge on the sphere and the second term to a finite potential.

10.2.1 Conducting Sphere in Uniform Field

Consider that the sphere is a conducting one and is isolated and uncharged. Further, consider that the potential at the location of the centre of the sphere due to the external field is V0.

Because the perturbing action of the sphere is negligible at a large distance from the sphere, the potential at a large distance from the sphere (ra) is given by

V(r,θ)=V0+E0rcosθ(10.18)

If the sphere is charged with a finite amount of charge Q, then

V(r,θ)=Q4πε0r+V0+E0rcosθ(10.19)

In practical systems, floating metallic particles are usually not charged and hence Equation 10.18 is taken here for further discussion.

Comparing Equations 10.17 and 10.18 for r → ∞, C2 = V0 and C3 = E0. C1 will be zero for uncharged sphere.

Therefore, Equation 10.17 can be rewritten as

V(r,θ)=V0+(E0r+C4r2)cosθ(10.20)

On the conductor surface, that is, for r = a,

V(a,θ)=V0+(E0a+C4a2)cosθ(10.21)

But the conducting sphere surface is an equipotential and hence electric potential is independent of θ on the conductor surface.

Therefore, from Equation 10.21,

  • C4 = −E0a3

Hence, the complete solution for an electric potential in the domain r > a is given by

V(r,θ)=V0+E0(r-a3r2)cosθ(10.22)

The r and θ components of the electric field intensity could be obtained as follows:

Er=-Vr=-E0(1+2a3r3)cosθ(10.23)

Eθ=-1rVθ=E0(1-a3r3)sinθ(10.24)

On the conducting sphere surface, tangential component of the electric field intensity must be zero as it is an equipotential surface. Equation 10.24 shows that for r = a, Eθ is zero, which in turn validates the solution obtained.

Again, on the conducting sphere surface, Er is the normal component of the electric field intensity, which is given by Er|r=a = − 3 E0 cos θ. Thus, the maximum value of the electric field intensity on the surface of the conducting sphere is 3E0, that is, three times the strength of the uniform external field.

This is the reason why metallic dust particles should be avoided at all costs for GISs. Because the presence of metallic dust particles will increase the local electric field intensity three times, which will result into partial discharge within the GIS, that is very detrimental for the GIS operation.

Induced surface charge density on the surface of the conducting sphere may be obtained as follows:

σsε0=Er|r=a=-3E0cosθ,or,σs=-3ε0E0cosθ(10.25)

As stated earlier, the sphere may be charged with an additional charge Q, which is distributed uniformly on the sphere surface and its effect on the field could be found by superposition.

10.2.2 Dielectric Sphere in Uniform Field

In the case of dielectric sphere present in uniform external field, there will be two solutions to potential function, Vi valid for the region within the sphere having dielectric of permittivity εi and Ve valid for the region outside the sphere having dielectric of permittivity εe. Therefore, from Equation 10.17

Vi(r,θ)=C1ir+C2i+(C3ir+C4ir2)cosθ(10.26)

and

Ve(r,θ)=C1er+C2e+(C3er+C4er2)cosθ(10.27)

The potential at large distance r (ra) from the sphere

V(r,θ)=V0+E0rcosθ(10.28)

where:

  • V0 is the potential at the location of the centre of the sphere due to the external field

Comparing Equations 10.27 and 10.28 for r → ∞, C2e = V0 and C3e = E0. C1e will be zero as a dielectric sphere is not considered to have any free charge.

Hence, Equation 10.27 can be rewritten as

Ve(r,θ)=V0+(E0r+C4er2)cosθ(10.29)

Inside the dielectric sphere electric potential must be finite at all the points. Hence, from Equation 10.26 Cli = C4i = 0, C2i = V0. Hence, Equation 10.26 can be rewritten as

Vi(r,θ)=V0+C3ircosθ(10.30)

At r = a, both Equations 10.29 and 10.30 should yield the same electric potential. Therefore,

V0+(E0a+C4ea2)cosθ=V0+C3i a cosθor,E0a+C4ea2=C3ia(10.31)

On the dielectric–dielectric boundary, the normal component of the electric field intensity should be same on both sides of the boundary. For the spherical boundary, r-component of the electric field intensity is the normal component on the boundary. Hence,

-εi(Vir)|r=a=-εe(Ver)|r=aor,εiC3i=εe(E0-2C4ea3)(10.32)

From Equations 10.31 and 10.32

C4e=εe-εi2εe+εia3E0

and

C3i=3εe2εe+εiE0

Therefore, the complete solutions for potential functions inside and outside the dielectric sphere are given by

Vi(r,θ)=V0+3εe2εe+εiE0rcosθ(10.33)

and

Ve(r,θ)=V0+(r+εe-εi2εe+εia3r2)E0cosθ(10.34)

Noting that rcos θ = z, potential function inside the dielectric sphere can be written as

Vi(x,y,z)=V0+3εe2εe+εiE0z(10.35)

Hence, the electric potential within the dielectric sphere varies in only z-direction, that is, the direction of the external field. Electric field intensity within the dielectric sphere will, therefore, have only the z-component, which is given by

Ezi=-Viz=-3εe2εe+εiE0(10.36)

Images

FIGURE 10.2
Electric field in and around dielectric sphere in uniform field.

Equation 10.36 shows that the magnitude of the electric field intensity within the dielectric sphere is constant. Typical field distribution in and around a dielectric sphere within uniform external field is shown in Figure 10.2.

Equation 10.36 also shows that if εi < εe, then |Ezi| > E0. Consider the case of a spherical air bubble trapped within a moulded solid insulation of relative permittivity 4. If the magnitude of the electric field intensity in solid insulation at the location of the air bubble is E0, then the magnitude of the electric field intensity within the air bubble will be 1.33E0. The operating electric field intensity within solid insulation is usually kept at a higher value as the solid insulation has a higher dielectric strength and hence such increase in field intensity within the air bubble often causes partial discharge within the air bubble as the dielectric strength of air is much lower than solid insulation.

10.3 Cylinder in Uniform External Field

Consider a long cylindrical object of radius a within uniform external field, as shown in Figure 10.3. Because the boundary is a circle of r = constant on the x–y plane, the system is best described in cylindrical coordinates, as shown in Figure 10.3. The uniform external field is given by E0=E0i^ and the potential at any point due to the external field is given by E0r cosθ = E0x with respect to the axis of the cylinder. In order to get the complete solution for the electric field in this system, Laplace’s equation in cylindrical coordinates as given in Equation 10.37 needs to be solved.

Images

FIGURE 10.3
Cylinder in uniform external field.

1rr(rVr)+1r22Vθ2+2Vz2=0(10.37)

For this arrangement, the electric field distribution does not vary along the length of the cylinder, that is, along z-coordinate. Hence, Laplace’s equation reduces to

1rr(rVr)+1r22Vθ2=0(10.38)

Separating the terms of the LHS into functions of only one variable by multiplying r2 with Equation 10.38, it may be written that

rr(rVr)+2Vθ2=0(10.39)

The two terms on the LHS of Equation 10.39 are functions of only one variable each, that is, the first term is function of r only, whereas the second term is function of θ only. The solution to Equation 10.39 can be obtained as the product of two functions of which one is dependent only on r and the other is dependent only on θ.

Let the assumed solution be

V(r,θ)=M(r)N(θ)(10.40)

The assumed solution is convenient as the boundary lies at r = constant.

Combining Equations 10.39 and 10.40,

rr[rM(r)N(θ)r]+2M(r)N(θ)θ2=0N(θ){rr[rM(r)r]}+M(r)2M(r)N(θ)θ2=0

Dividing by M(r)N(θ),

rM(r)ddr[rdM(r)dr]+1N(θ)d2N(θ)dθ2=0(10.41)

The partial derivatives become total derivatives in Equation 10.41 as each term is dependent on only one coordinate.

As in the case of sphere in uniform field, zero separation constant solution and equal and opposite separation constant solution are to be obtained separately in this case, too.

Determination of the zero separation constant solution:

The first term of Equation 10.41 is

rM(r)ddr[rdM(r)dr]=0

where:

M(r) is non-zero

Therefore,

ddr[rdM(r)dr]=0,or,rdM(r)dr=C0,or,dM(r)dr=C0r

Integrating and incorporating constants of integration

M(r)=C10lnr+C20(10.42)

Next the second term of Equation 10.41 is

1N(θ)d2N(θ)dθ2=0

where:

N(θ) is non-zero

Therefore,

N(θ)=A10θ+A20(10.43)

But, from Equations 10.42 and 10.43, it can be seen that there is discontinuity of potential at r = 0 and θ = ∞, which are not feasible in the given arrangement as potential must be a continuous function. Hence, C10 = A10 = 0 in Equations 10.42 and 10.43.

Therefore,

V(r,θ)=C20A20=C1(10.44)

Determination of the equal and opposite separation constant solution:

The first term of Equation 10.41 is

rM(r)ddr[rdM(r)dr]=+p

where:

p is a positive constant

or,r2d2M(r)dr2+rdM(r)dr=+pM(r)(10.45)

Substituting M(r) = Crn in Equation 10.45, it may be obtained that

n(n-1)+n=p,or,n=±p

Hence,

M(r)=Crp+C"rp(10.46)

Again, the second term of Equation 10.41 is

1N(θ)d2N(θ)dθ2=-por,d2N(θ)dθ2=-pN(θ)(10.47)

Substituting N(θ) = eaθ in Equation 10.47, it may be obtained that

a2eaθ=-peaθ,or,a=±ip

Hence,

N(θ)=Bcos(pθ+α)(10.48)

Equations 10.46 and 10.48 lead to

V(r,θ)=M(r)N(θ)=(C2rp+C3rp)cos(pθ+α)(10.49)

where:

C2 = CB and C3 = CB

From Equations 10.44 and 10.49, the complete solution for potential function at all values of r and θ can be obtained as

V(r,θ)=C1+(C2rp+C3rp)cos(pθ+α)(10.50)

The potential at large distance r (ra) from the cylinder is given by

V(r,θ)=V0+E0rcosθ(10.51)

Matching Equations 10.50 and 10.51, p=1 and α = 0.

Hence, the complete solution as given by Equation 10.50 reduces to

V(r,θ)=C1+(C2r+C3r)cosθ(10.52)

10.3.1 Conducting Cylinder in Uniform Field

Comparing Equations 10.51 and 10.52, C1 = V0 and C3 = E0.

Therfore,

V(r,θ)=V0+(C2r+E0r)cosθ

On the conductor surface, that is, for r = a,

V(a,θ)=V0+(C2a+E0a)cosθ(10.53)

But the conducting cylinder surface is an equipotential and hence electric potential is independent of θ on the conductor surface.

Therefore, from Equation 10.53,

  • C2 = −E0a2

Hence, the complete solution for electric potential in the domain r > a is given by

V(r,θ)=V0+E0(r-a2r)cosθ(10.54)

The r and θ components of the electric field intensity could be obtained as follows:

Er=-Vr=-E0(1+a2r2)cosθ(10.55)

Eθ=-1rVθ=E0(1-a2r2)sinθ(10.56)

Equation 10.56 shows that for r = a, Eθ is zero, that is, the tangential component of the electric field intensity is zero on the cylindrical conductor surface as it is an equipotential surface.

Again, on the conducting cylinder surface, Er is the normal component of the electric field intensity, which is given by Er|r=n = −2E0 cosθ. Thus, the maximum value of the electric field intensity on the surface of the conducting cylinder is 2E0, that is, twice the magnitude of the uniform external field. Comparing this maximum electric field intensity with the value obtained for conducting sphere in uniform field, it may be seen that the enhancement of field intensity is more if the conducting object is spherical is shape.

Induced surface charge density on the surface of the conducting cylinder may be obtained as follows:

σsε0=Er|r=a=-2E0cosθ,or,σs=-2ε0E0cosθ(10.57)

10.3.2 Dielectric Cylinder in Uniform Field

Potential function valid for the region within the cylinder having dielectric of permittivity εi is

Vi(r,θ)=C1i+(C2ir+C3ir)cosθ(10.58)

and the potential function valid for the region outside the cylinder having dielectric of permittivity εe is

Ve(r,θ)=C1e+(C2er+C3er)cosθ(10.59)

The potential at large distance r (ra) from the cylinder

V(r,θ)=V0+E0rcosθ(10.60)

where:

  • V0 is the potential at the location of the axis of the cylinder due to the external field

Comparing Equations 10.59 and 10.60 for r → ∞, C1e = and C3e = E0

Hence, Equation 10.59 can be rewritten as

Ve(r,θ)=V0+(C2er+E0r)cosθ(10.61)

Inside the dielectric cylinder electric potential must be finite at all the points. Hence, from Equation 10.58, Cli = V0 and C2i = 0. Hence, Equation 10.58 can be rewritten as

Vi(r,θ)=V0+C3ircosθ(10.62)

At any point on the dielectric cylinder surface, that is, for r = a, electric potential as may be obtained from Equations 10.61 and 10.62 must be unique. Hence,

C2ea+E0a=C3ia(10.63)

From the boundary condition of normal component of electric flux density at r = a,

-εi(Vir)|r=a=-εe(Ver)|r=aor,εiC3i=εe(-C2ea2+E0)(10.64)

From Equations 10.63 and 10.64

C2e=εe-εiεe+εia2E0and C3i=2εeεe+εiE0

Therefore, the complete solutions for potential functions inside and outside the dielectric cylinder are given by

Vi(r,θ)=V0+2εeεe+εiE0rcosθ(10.65)

and Ve(r,θ)=V0+(r+εe-εiεe+εia2r)E0cosθ(10.66)

As rcos θ = x, potential function inside the dielectric cylinder can be written as

Vi(x,y)=V0+2εeεe+εiE0x(10.67)

Hence, the electric potential within the dielectric cylinder varies in only x-direction, that is, the direction of the external field. Electric field intensity within the dielectric cylinder will, therefore, have only the x-component, which is given by

Exi=-Vix=-2εeεe+εiE0(10.68)

Similar to the case of dielectric sphere in the uniform field, Equation 10.68 shows that the magnitude of the electric field intensity within the dielectric cylinder is constant. Typical field distribution on the x–y plane in and around a dielectric cylinder within uniform external field will be the same as that shown in Figure 10.2.

As in the case of dielectric sphere in uniform field, for dielectric cylinder in uniform field also |Ezi| > E0 if εi < εe. If a cylindrical air bubble is trapped within a moulded solid insulation of relative permittivity 4, then the magnitude of the electric field intensity within the air bubble will be 1.6E0, where E0 is the magnitude of the electric field intensity in solid insulation at the location of the air bubble. Comparing this result with the corresponding value in the case of dielectric sphere, it may be seen that field enhancement is more if the gas cavity in liquid or solid insulation is cylindrical in shape.

Objective Type Questions

1. A conducting sphere is placed within an external field of the uniform field intensity E0. Then the maximum value of field intensity on the surface of the conducting sphere will be

a. 0

b. E0

c. 2E0

d. 3E0

2. A conducting cylinder is placed within an external field of the uniform field intensity E0. Then the maximum value of field intensity on the surface of the conducting cylinder will be

a. 0

b. E0

c. 2E0

d. 3E0

3. A dielectric sphere is placed within an external field of the uniform field intensity E0. The external field intensity acts in the direction of the x-axis. Then the field intensity within the dielectric sphere will be in the direction of

a. x-axis

b. y-axis

c. z-axis

d. None of the above

4. A dielectric cylinder is placed within an external field of the uniform field intensity E0. The external field intensity acts in the direction of y-axis. Then the field intensity within the dielectric cylinder will be in the direction of

a. x-axis

b. y-axis

c. z-axis

d. None of the above

5. A dielectric sphere is placed within an external field of the uniform field intensity E0. Then the field intensity within the dielectric sphere will be

a. Zero

b. Equal to E0

c. Constant but not equal to E0

d. None of the above

6. A dielectric cylinder is placed within an external field of the uniform field intensity E0. Then the field intensity within the dielectric cylinder will be

a. Zero

b. Equal to E0

c. Constant but not equal to E0

d. None of the above

7. A gas cavity is present in large volume of solid insulation. Then the ratio of field intensity within the gas cavity for spherical to that for cylindrical cavity will be

a. Zero

b. <1

c. 1

d. >1

8. A conducting object is present in an external field of uniform field intensity. Then the ratio of maximum field intensity on the surface of the conducting object for spherical to that for cylindrical object will be

a. Zero

b. <1

c. 1

d. >1

Answers:

1) d;

2) c;

3) a;

4) b;

5) c;

6) c;

7) b;

8) d

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset