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.
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.
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 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.
FIGURE 10.1
Sphere in uniform external field.
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
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
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
The assumed solution is convenient as the boundary lies at r = constant.
Combining Equations 10.3 and 10.4
Dividing by M(r)N(θ),
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:
where:
p is a positive constant
Another solution is obtained when the separation constant is zero.
Hence, the zero separation constant solution
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
where:
M(r) is non-zero
Therefore,
Integrating and incorporating constants of integration
Next the second term of Equation 10.7 is
where:
N(θ) is non-zero
Therefore,
Integrating and incorporating constant of integration
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,
Then from Equations 10.4, 10.8 and 10.10, the zero separation constant solution can be obtained as
where:
C1 = A20C10 and C2 = A20C20
Determination of the equal and opposite separation constant solution:
The first term of Equation 10.6 is
Putting M(r) = Crn in Equation 10.12, we get
Hence,
The second term of Equation 10.6 is
Putting N(θ) = Bcosθ in Equation 10.14,
Hence, from Equation 10.13 n = +1, −2.
From Equations 10.4 and 10.15,
where:
C3 = C′B and C4 = C′B
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.
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.
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 (r ≫ a) is given by
If the sphere is charged with a finite amount of charge Q, then
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
On the conductor surface, that is, for r = a,
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
The r and θ components of the electric field intensity could be obtained as follows:
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:
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.
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
and
The potential at large distance r (r ≫ a) from the sphere
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
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
At r = a, both Equations 10.29 and 10.30 should yield the same electric potential. Therefore,
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,
From Equations 10.31 and 10.32
and
Therefore, the complete solutions for potential functions inside and outside the dielectric sphere are given by
and
Noting that rcos θ = z, potential function inside the dielectric sphere can be written as
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
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.
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 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.
FIGURE 10.3
Cylinder in uniform external field.
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
Separating the terms of the LHS into functions of only one variable by multiplying r2 with Equation 10.38, it may be written that
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
The assumed solution is convenient as the boundary lies at r = constant.
Combining Equations 10.39 and 10.40,
Dividing by M(r)N(θ),
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
where:
M(r) is non-zero
Therefore,
Integrating and incorporating constants of integration
Next the second term of Equation 10.41 is
where:
N(θ) is non-zero
Therefore,
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,
Determination of the equal and opposite separation constant solution:
The first term of Equation 10.41 is
where:
p is a positive constant
Substituting M(r) = Crn in Equation 10.45, it may be obtained that
Hence,
Again, the second term of Equation 10.41 is
Substituting N(θ) = eaθ in Equation 10.47, it may be obtained that
Hence,
Equations 10.46 and 10.48 lead to
where:
C2 = C′B and C3 = C″B
From Equations 10.44 and 10.49, the complete solution for potential function at all values of r and θ can be obtained as
The potential at large distance r (r ≫ a) from the cylinder is given by
Matching Equations 10.50 and 10.51, and α = 0.
Hence, the complete solution as given by Equation 10.50 reduces to
Comparing Equations 10.51 and 10.52, C1 = V0 and C3 = E0.
Therfore,
On the conductor surface, that is, for r = a,
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
The r and θ components of the electric field intensity could be obtained as follows:
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:
Potential function valid for the region within the cylinder having dielectric of permittivity εi is
and the potential function valid for the region outside the cylinder having dielectric of permittivity εe is
The potential at large distance r (r ≫ a) from the cylinder
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
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
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,
From the boundary condition of normal component of electric flux density at r = a,
From Equations 10.63 and 10.64
Therefore, the complete solutions for potential functions inside and outside the dielectric cylinder are given by
As rcos θ = x, potential function inside the dielectric cylinder can be written as
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
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.
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
1) d;
2) c;
3) a;
4) b;
5) c;
6) c;
7) b;
8) d