prev next


 

The Vector Potential and the Vector Poisson Equation

8.1.1A solenoid has radius a, length d, and turns N, as shown in Fig. 8.2.3. The length d is much greater than a, so it can be regarded as being infinite. It is driven by a current i.

continuity law [(8.0.1) and (8.0.2)], as well as the associated continuity conditions [(8.0.3) and (8.0.4)], are satisfied by an interior magnetic field intensity that is uniform and an exterior one that is zero.
(a) Show that Ampère's differential law and the magnetic flux
(b) What is the interior field?
(c) A is continuous at r = a because otherwise the H field would have a singularity. Determine A.
8.1.2*A two-dimensional magnetic quadrupole is composed of four line currents of magnitudes i, two in the positive z direction at x = 0, y = d/2 and two in the negative z direction at x = d/2, y = 0. (With the line charges representing line currents, the cross-section is the same as shown in Fig. P4.4.3.) Show that in the limit where r \gg d, A_z = - (o id2/4 )(r-2) cos 2. (Note that distances must be approximated accurately to order d2.) 8.1.3A two-dimensional coil, shown in cross-section in Fig. P8.1.3, is composed of N turns of length l in the z direction that is much greater than the width w or spacing d. The thickness of the windings in the y direction is much less than w and d. Each turn carries the current i. Determine A.

figure GIF #1
 

The Biot-Savart Superposition Integral

8.2.1*The washer-shaped coil shown in Fig. P8.2.1 has a thickness that is much less than the inner radius b and outer radius a. It supports a current density J = J_o i_. Show that along the z axis,

equation GIF #8.145
floating figure GIF #38
8.2.2*A coil is wound so that the wire forms a spherical shell of radius R with the wire essentially running in the direction. With the wire driven by a current source, the resulting current distribution is a surface current at r = R having the density K = K_o sin i_, where K_o is a given constant. There are no other currents. Show that at the center of the coil, H = (2K_o/3)iz.

8.2.3In the configuration of Prob. 8.2.2, the surface current density is uniformly distributed, so that K = K_o i_, where K_o is again a constant. Find H at the center of the coil. 8.2.4Within a spherical region of radius R, the current density is J = J_o i_, where J_o is a given constant. Outside this region is free space and no other sources of H. Determine H at the origin.

8.2.5*A current i circulates around a loop having the shape of an equilateral triangle having sides of length d, as shown in Fig. P8.2.5. The loop is in the z = 0 plane. Show that along the z axis,

equation GIF #8.146
floating figure GIF #39
8.2.6For the two-dimensional coil of Prob. 8.1.3, use the Biot-Savart superposition integral to find H along the x axis. 8.2.7*Show that A induced at point P by the current stick of Figs. 8.2.5 and 8.2.6 is

equation GIF #8.147
 

The Scalar Magnetic Potential

8.3.1Evaluate the H field on the axis of a circular loop of radius R carrying a current i. Show that your result is consistent with the result of Example 8.3.2 at distances from the loop much greater than R.

8.3.2Determine \Psi for two infinitely long parallel thin wires carrying currents i in opposite directions parallel to the z axis of a Cartesian coordinate system and located along x = a. Show that the lines \Psi = const in the x - y plane are circles. 8.3.3Find the scalar potential on the axis of a stack of circular loops (a coil) of N turns and length l using 8.3.12 for an individual turn, integrating over all the turns. Find H on the axis.

 Conductors 8.4.1*A current loop of radius R is at the center of a conducting spherical shell having radius b. Assume that R \ll b and that i(t) is so rapidly varying that the shell can be taken as perfectly conducting. Show that in spherical coordinates, where R \ll r < b

equation GIF #8.148
8.4.2The two-dimensional magnetic dipole of Example 8.1.2 is at the center of a conducting shell having radius a \gg d. The current i(t) is so rapidly varying that the shell can be regarded as perfectly conducting. What are \Psi and H in the region d \ll r < a? 8.4.3*The cross-section of a two-dimensional system is shown in Fig. P8.4.3. A magnetic flux per unit length so H_o is trapped between perfectly conducting plane parallel plates that extend to infinity to the left and right. At the origin on the lower plate is a perfectly conducting half-cylinder of radius R.

equation GIF #8.149
(a) Show that if s \gg R, then
(b) Show that a plot of H would appear as in the left half of Fig. 8.4.2 turned on its side.

floating figure GIF #40
8.4.4In a three-dimensional version of that shown in Fig. P8.4.3, a perfectly conducting hemispherical bump of radius s \gg R is attached to the lower of two perfectly conducting plane parallel plates. The hemisphere is centered at the origin of a spherical coordinate system such as in Fig. P8.4.3, with . The magnetic field intensity is uniform far from the hemisphere. Determine \Psi and H. 8.4.5*Running from z = - to z = +\infty at (x, y) = (0, -h) is a wire. The wire is parallel to a perfectly conducting plane at y = 0. When t = 0, a current step i = I u_-1(t) is applied in the +z direction to the wire.

equation GIF #8.150
(a) Show that in the region y < 0,
(b) Show that the surface current density at y = 0 is K_z = - ih/ (x2 + h2).

floating figure GIF #41
8.4.6The cross-section of a system that extends to infinity in the z directions is shown in Fig. P8.4.6. Surrounded by free space, a sheet of current has the surface current density K_o iz uniformly distributed between x = b and x = a. The plane x = 0 is perfectly conducting.
(a) Determine \Psi in the region 0 < x.
(b) Find K in the plane x = 0.
 

Piece-Wise Magnetic Fields

8.5.1*The cross-section of a cylindrical winding is shown in Fig. P8.5.1. As projected onto the y = 0 plane, the number of turns per unit length is constant and equal to N/2R. The cylinder can be modeled as infinitely long in the axial direction.

floating figure GIF #42
floating figure GIF #43
equation GIF #8.151
and that therefore
equation GIF #8.152
(a) Given that the winding carries a current i, show that
(b) Show that the inductance per unit length of the winding is L = o N2/8.

8.5.2The cross-section of a rotor, coaxial with a perfectly conducting ``magnetic shield,'' is shown in Fig. P8.5.2. Windings consisting of N turns per unit peripheral length are distributed uniformly at r = b so that at a given instant in time, the surface current distribution is as shown. At r = a, there is the inner surface of a perfect conductor. The system is very long in the z direction. conditions on \Psi at r = b and the boundary condition at r = a?
(a) What are the continuity
(b) Find \Psi, and hence H, in regions (a) and (b) outside and inside the winding, respectively.
(c) With the understanding that the rotor is wound using one wire, so that each turn is in series with the next and a wire carrying the current in the +z direction at returns the current in the -z direction at -, what is the inductance of the rotor coil? Why is it independent of the rotor position _o?
 

Vector Potential

8.6.1*In Example 1.4.1, the magnetic field intensity is determined to be that given by (1.4.7). Define A_z to be zero at the origin.

of r = R, A_z must be continuous there.
(a) Show that if H_ is to be finite in the neighborhood
(b) Show that A is given by
equation GIF #8.153
(c) The loop designated by C' in Fig. 1.4.2 has a length l in the z direction, an inner leg at r = 0, and an outer leg at r = a > R. Use A to show that the flux linked is
equation GIF #8.154
8.6.2For the configuration of Prob. 1.4.2, define A_z as being zero at the origin.
(a) Determine A_z in the regions r < b and b < r < a.
(b) Use A to determine the flux linked by a closed rectangular loop having length \l in the z direction and each of its four sides in a plane of constant . Two of the sides are parallel to the z axis, one at radius r = c and the other at r = 0. The other two, respectively, join the ends of these segments, running radially from r = 0 to r = c.

8.6.3*In cylindrical coordinates, o H = \mo [H_r(r, z)ir + H_z(r, z)iz]. That is, the magnetic flux density is axially symmetric and does not have a component.

equation GIF #8.155
(a) Show that
(b) Show that the flux passing between contours at r = a and r = b is
equation GIF #8.156
8.6.4*For the inductive attenuator considered in Example 8.6.3 and Demonstration 8.6.2: without identifying this MQS problem with its EQS counterpart.
(a) derive the vector potential, (20),
(b) Show that the current is as given by (21).
(c) In the limit where b/a \gg 1, show that the response has the dependence on b/a shown in the plot of Fig. 8.6.11.
(d) Show that in the opposite limit, where b/a \ll 1, the total current in the lower plate (21) is consistent with a magnetic field intensity between the upper and lower plates that is uniform (with respect to y) and hence equal to (\Lambda /bo )ix. Note that
equation GIF #8.157
floating figure GIF #44
8.6.5Perfectly conducting electrodes are composed of sheets bent into the shape of \sqcup's, as shown in Fig. P8.6.5. The length of the system in the z direction is very large compared to the length 2a or height d, so the fields can be regarded as two dimensional. The insulating gaps have a width that is small compared to all dimensions. Passing through these gaps is a magnetic flux (per unit length in the z direction) \Lambda (t). One method of solution is suggested by Example 6.6.3. (b) to the right and left, respectively, of the plane x = 0.
(a) Find A in regions (a) and
(b) Sketch H.
8.6.6*The wires comprising the winding shown in cross-section by Fig. P8.6.6 carry current in the -z direction over the range 0 < x < a and return this current over the range -a < x < 0. These windings extend uniformly over the range 0 < y < b. Thus, the current density in the region of interest is J = - in_o sin ( x/a)iz, where i is the current carried by each wire and |n_o sin ( x/a)| is the number of turns per unit area. This region is surrounded by perfectly conducting walls at y = 0 and y = b and at x = -a and x = a. The length l in the z direction is much greater than either a or b.

floating figure GIF #45
equation GIF #8.158
(a) Show that
(b) Show that the inductance of the winding is
equation GIF #8.159
(c) Sketch H.
8.6.7In the configuration of Prob. 8.6.6, the rectangular region is uniformly filled with wires that all carry their current in the z direction. There are n_o of these wires per unit area. The current carried by each wire is returned in the perfectly conducting walls.
(a) Determine A.
(b) Assume that all the wires are connected to the wall by a terminating plate at z = l and that each is driven by a current source i(t) in the plane z = 0. Note that it has been assumed that each of these current sources is the same function of time. What is the voltage v(x, y, t) of these sources?
8.6.8In the configuration of Prob. 8.6.6, the turns are uniformly distributed. Thus, n_o is a constant representing the number of wires per unit area carrying current in the -z direction in the region 0 < x. Assume that the wire carrying current in the -z direction at the location (x, y) returns the current at (-x, y).
(a) Determine A.
(b) Find the inductance L.

Given the distribution of current density throughout all of space, a superposition integral solution is found, much as the solution of Poisson's equation in Chap. 4. Fields associated with current singularities play a role in MQS systems similar to that for charge singularities in EQS systems. The derivation of two- and three-dimensional dipole fields using the vector potential form of the superposition integral therefore exemplifies the application of the integral while giving results that find many applications.

A superposition integral for the field intensity rather than the vector potential proves convenient, especially if the fields are determined numerically from given current distributions. This superposition integral is provided by Biot-Savart's law in Sec. 8.3.

The determination of fields from given current distributions is important for at least two reasons. First, by contrast with EQS systems, where the charges seldom have known distributions, the use of conductors wrapped with insulation (wires) makes feasible the imposition of prescribed current distributions. In the absence of magnetizable materials such as iron, electromagnets are an example. The Biot-Savart law is often used to design electromagnets.

A second reason for seeking solutions to the magnetic version of the superposition integral is its role as the particular solution in boundary value problems. Typical of this class of situations is a time varying magnetic field within some region bounded at least in part by a highly conducting material. For reasons cited in Sec. 8.4, currents are induced in the surfaces of ``perfect'' conductors in such a way that the normal magnetic flux density is zero. The homogeneous solutions needed to satisfy boundary conditions are familiar from Chap. 5, because they are solutions to Laplace's equation. The method of images and boundary value point of view needed to solve such problems is the subject of Sec. 8.5.

In many practical configurations, the currents are confined to relatively thin regions where they can be modeled by surface currents or line currents. In these situations, where H is irrotational in the volumes of interest, it is natural to exploit a scalar potential that is analogous to the electric potential introduced in Chap. 4. This scalar magnetic potential \Psi, which is the theme of Secs. 8.6 and 8.7, is used extensively in Chap. 9 for systems containing magnetizable materials. The vector potential is the natural variable for evaluating the flux passing through a surface. (We will see in Sec. 8.4 that the time-rate-of-change of flux determines the terminal voltage of an inductor.) Integration of (1) over the open surface S of Fig. 8.1.1 gives

floating figure GIF #46
equation GIF #8.160
and it follows from Stokes' theorem that this flux is equal to the line integral of A ds around the contour enclosing the surface.
boxed equation GIF #8.21
Because (1) and (2) represent the magnetic flux density in terms of another vector function (rather than a scalar) the situation is inherently more complicated than for the EQS laws. In certain important cases, however, A has only one component. Then, a vector field is again represented in terms of one scalar function. Two of these cases are identified in the following subsections.

Vector Potential for Two-Dimensional Fields

Suppose that the flux density is parallel to the x - y plane and is independent of z. Then, it can be represented by a vector potential having only a z component.

equation GIF #8.161
Note that the divergence of this A is automatically zero and that in Cartesian coordinates the components of the flux density are given in terms of A_z by (1), which becomes
equation GIF #8.162
If A_z were written in terms of the polar coordinates (r, ), then the polar components of the flux density would follow from (1) using the curl operator in polar coordinates.

floating figure GIF #47
Consider now the evaluation of the net flux of magnetic flux density through a surface S that has length l in the z direction, as shown in Fig. 8.1.2. The points (a) and (b) denote the coordinates of the corners of the contour enclosing S. The contour consists of a pair of parallel straight segments of length l parallel to the z axis, one at the location (a) in the x - y plane and the other at (b), and contours joining (a) and (b) in x - y planes. Contributions to the contour integral, (4), from these latter segments of C are zero, because A is perpendicular to ds. Integration along the z directed segments amounts to multiplication of A_z evaluated at (a) or (b) by the length of the segment. Thus, (3) becomes

boxed equation GIF #8.22
The vector potential at (a) relative to (b) is the net magnetic flux per unit length passing through a surface of unit length in the z direction subtended between the two points and a corresponding pair at unity distance along the z axis. Note that the flux has a sign, relative to the direction of the contour integration, governed by the right-hand rule (Fig. 1.4.1).

Example 8.7.1. Vector Potential of a Line Current

The magnetic field intensity for the line current i of Fig. 8.1.3 was determined in Sec. 1.4 using Ampère's integral law. Thus, the magnetic flux density is given by (1.4.10) as
equation GIF #8.163
floating figure GIF #48
The vector potential is determined by first writing (1) in polar coordinates.

equation GIF #8.164
Because there is no r component of the flux density, it follows that A_z is independent of . Setting the components of (8) and (9) equal and integrating gives
equation GIF #8.165
To illustrate the use of the vector potential to evaluate the net flux through a surface, consider the flux through the coil shown in Fig. 1.7.2a. The contour running parallel to the z axis in the z direction is at a radius d from the wire while that returning the current in the -z direction is at radius R. The length in the z direction is l. Thus, (7) gives a net flux of
equation GIF #8.166
Note that this result is in agreement with (1.7.5), which was obtained by evaluating the specific surface integral.

Visualization of Two-Dimensional Solenoidal Fields

The representation of the flux density in terms of a component of the vector potential not only completes a step in the evaluation of the flux linked by a coil, it also is the basis for visualizing o H. For the two-dimensional fields pictured by A_z in Fig. 8.1.2, think of (b) as held fixed and let the line denoted by (a) trace out a surface of constant A_z. The flux per unit length passing between the line denoted by (b) and any line in this surface of constant A_z is the same. Thus, there is no flux density normal to a surface of constant A_z. That is, lines of constant A_z are the lines of magnetic flux density.

Illustration. Field Lines Due to a Line Current

For the line current of Fig. 8.1.3, the surfaces of constant A_z given by (10) appear in a constant z plane as concentric circles, as shown in Fig. 8.1.4. With arrows added to indicate direction, these represent the flux density. Drawn for equal increments of A_z, the magnetic flux trapped between pairs of lines is the same. Thus, as the distance between lines increases, the flux density decreases.

floating figure GIF #49
\sectnonumVector Potential for Axisymmetric FIelds in Spherical Coordinates If the magnetic flux density is invarient with respect to rotation around the z axis, having components in the r and directions only, the vector potential again has a single component.

equation GIF #8.167
The net flux through the annular surface ``spanned'' over the contour shown in Fig. 8.1.5, having constant outer and inner radii denoted by (a) and (b) respectively, is given by the contributions to (4) of the azimuthal segments, A_ multiplied by the circumferences. The contour is closed by adjacent oppositely directed segments joining points (a) and (b) in a plane of constant . Thus, the contributions to the line integral of (4) from these segments cancel, even if A had components in the direction of ds on these segments. Thus, the net flux through the annulus is simply the axisymmetric stream function \Lambda at (a) relative to that at (b)
5 With A used to represent the velocity distribution of an incompressible fluid, \Lambda (or \Lambda /2) is called Stokes' stream function.
floating figure GIF #50
boxed equation GIF #8.23
where
equation GIF #8.168
Lines of flux density are tangential to the axisymmetric surfaces of constant \Lambda. Just as A_z provides a ready visualization of the flux lines in two-dimensions, \Lambda portrays the axisymmetric flux lines.

Example 8.7.1. Dipole Field of a Current Loop

A three-dimensional magnetic dipole field is observed if a small circular current loop sustains a circulating current i. By ``small'' we mean that the radius R is much less than the distance to the point of observation. The configuration, shown in Fig. 8.2.1, has rotational symmetry with respect to the z axis, so the magnetic field intensity is independent of . The vector potential is therefore expected to take the form of (8.1.12).

floating figure GIF #51
In the evaluation of (6), the source is confined to the wire which is assumed to have a cross-sectional area, a, of dimensions small compared to the radius R. With da' used to denote an incremental area of the wire perpendicular to the current, the incremental source in (6) can be written as
equation GIF #8.169
Because the current density is confined to a wire having a cross-section that is much smaller than the radius R, integration over the cross-section of the wire can be taken as a multiplication of the current density by the area with distance |r - r' | essentially held fixed. Thus, (6) becomes
equation GIF #8.170
where (r2 = x2 + y2 + z2).

equation GIF #8.171
Because r \gg R, we can expand the source-observer distance and drop R2 compared to r2 in this expression. Thus,

6 [(a + b)n \approx an + nban-1]
equation GIF #8.172
and (9) becomes
equation GIF #8.173
Only the third and fifth terms in the integrand contribute to the integral, which is
equation GIF #8.174
Here, the expression has been written so as to recognize that the Cartesian unit vectors combine to become the unit vector in the direction.
equation GIF #8.175
Written in terms of the magnetic dipole moment m, defined as the product of the area enclosed by the loop and the current i, this is the desired dipole vector potential.

The magnetic field intensity is obtained by taking the curl of this expression, as called for by (8.1.1) in spherical coordinates.

equation GIF #8.176
The lines of magnetic field flux density, and hence H, are tangent to the surfaces of constant \Lambda = 2 r sin A_. These are shown in profile by Fig. 8.2.2. From a comparison of the expressions for the electric and magnetic dipole fields, the negative gradient of (4.4.10) and (14), it is clear that the fields take the same form and in fact are identical if qd /o m.

floating figure GIF #52

Two-Dimensional Current and Vector Potential Distributions

A current density that is z directed but independent of z gives rise to a vector potential that is also in the z direction and dependent only on x and y. This is evident from the vector Poisson's equation, which reduces to (4). It also follows from the a two dimensional version of the superposition integral.

Observe that the vector potential for a z-directed uniform line current on the z axis is given by (8.1.10). Positioned at an arbitrary source location denoted by r' and illustrated in Fig. 8.2.3, this line current results in the vector potential

equation GIF #8.177
floating figure GIF #53
Note that, as in the case of the two-dimensional superposition integral for the electric potential, (4.5.20), the position vectors r and r' are two dimensional in the sense that they lie in any x - y plane. The current i is now replaced by the current density at the source coordinate multiplied by the differential area da'. The latter could be dx' , dy', for example. The total vector potential is then the superposition of contributions from all of the current distribution. Thus, the two-dimensional vector potential superposition integral is
boxed equation GIF #8.24
In essence, this represents the z component of the three dimensional superposition integral, (6), with the integration on z' carried out.

Note that the two-dimensional current distribution is automatically solenoidal.

With the identification /o o J_z, the two-dimensional electric-potential superposition integral becomes identical to (6). As the following examples illustrate, for every two-dimensional electric potential problem that has been solved, there is an antidual magnetic configuration that has also been described.

Thus, the z component of the vector Poisson's equation reduces to Laplace's equation.

equation GIF #8.178
Because the fields are two-dimensional, A_z is the only component of A.

In terms of A_z, the magnetic flux density can be written in polar coordinates (8.1.9) or Cartesian coordinates (8.1.6)

equation GIF #8.179
The Cartesian expression is useful for recognizing the form of A_z as y and hence
equation GIF #8.180
It follows that as r ,
equation GIF #8.181
In terms of the vector potential, the condition that there be no flux density normal to the cylinder surface, where r = R, is met if
equation GIF #8.182
These boundary conditions are met by taking a linear combination of the solutions to (4) that have the same dependence as A_z at infinity,
equation GIF #8.183
To match the uniform vector potential at infinity, C = - o H_o. The coefficient D is then adjusted to satisfy (8). It follows that
equation GIF #8.184
Remember, the field lines follow the ``lines'' of constant A_z, which are shown in Fig. 8.4.2. The field lines are excluded by the cylinder.

floating figure GIF #54

Example 8.7.2. Terminal Voltage Computed From Vector Potential

The perfectly conducting cylinder in an imposed time-varying magnetic field intensity H_o (t) from Example 8.4.1 is shown in Fig. 8.4.5. An N turn coil has been added that has legs of length l in the z direction closed at the ends by segments perpendicular to the z axis. Thus, wires run in the +z direction at (r = a, = 0) and return at (r = a, = ). The terminals of this coil are connected to a very high impedance voltmeter, so that what current is carried by the wires induces a field that is negligible compared to that due to H_o (t). What is the voltage measured at the terminals?

floating figure GIF #55
The total flux linked is N times the flux linked by one turn. In view of (8.1.7), is therefore given by simply
equation GIF #8.185
The specific A_z is given by (10), so (16) becomes
equation GIF #8.186
Thus, the terminal voltage follows from (14) as
equation GIF #8.187
The terminal voltage reflects the rate of change of the flux linked by the coil, regardless of the origins of that flux. In this case, it is due to the currents used to excite the electromagnet responsible for H_o and the currents induced in the cylinder to expel the flux from its interior. We consider next cases where the flux is due to current in the coil itself.

Boundary Value Solution By Superposition of Sources

By judicious arrangements of currents, it is possible to obtain solutions to practical boundary value problems. In Chap. 4, superposition of charges was used to solve boundary value problems involving perfect conductors in EQS systems. As we shall now show, a family of solutions for two-dimensional MQS systems can be obtained by superposition of the fields due to current sources.

In two-dimensional configurations, any surface of constant A_z can be replaced by the surface of a perfect conductor. Moreover, in the free space region between conductors, A_z satisfies Laplace's equation. Thus, any two-dimensional configuration from Chaps. 4 and 5 can be replaced by one where the potential lines are field lines. The equipotential surfaces of the EQS perfect conductors become the perfectly conducting surfaces of an MQS system with tangential lines of magnetic field intensity.

Illustration. Field Trapped Between Hyperbolic Perfect Conductors

The two-dimensional potential distribution of Example 4.1.1 suggests the vector potential A_z = _o xy/a2. The lines of magnetic field intensity, which are the surfaces of constant A_z, are shown in Fig. 8.4.6. Here, the surfaces A_z = _o are taken as being the surfaces of perfect conductors. Thus, the current density on the surfaces of these conductors are, given by using (8.1.6) to determine H and in turn (8.4.2) to find K_z. These currents shield the fields from the volume of the perfect conductors. The net flux per unit length passing downward between the upper pair of conductors is (in view of (8.1.7)) simply 2_o.

floating figure GIF #56
This solution is the superposition of the fields of four line currents. Two directed in the +z direction are at infinity in the first and third quadrants while two in the -z direction are in the second and third quadrants.

Example 8.7.4. Field and Inductance of Oppositely Directed Currents in Parallel Perfectly Conducting Cylinders

The cross-section of a pair of parallel perfectly conducting cylinders that extend to in the z direction is shown in Fig. 8.4.7. The conductors have the same geometry as in the EQS case considered in Example 4.6.3. However, they should be regarded as shorted at one end and driven by a current source i at the other. Thus, current in the +z direction in the right conductor is returned in the left conductor. Although the net current in each conductor is given, its distribution on the surface of the conductors is to be determined.

floating figure GIF #57
Example 4.6.3 suggests our strategy. Instead of superimposing the potentials of a pair of lines charges of opposite sign, we superimpose the A_z of oppositely directed line currents. With r_1 and r_2 the distances from the observer coordinate to the source coordinates, defined in Fig. 8.4.7, it follows from the vector potential for a line current given by (8.2.15) that
equation GIF #8.188
With the identification of variables

7 This is one of the few places where there is the possibility of confusing , the magnetic flux, with the line charge density. In this section we add a prime to to denote the line charge density.
equation GIF #8.189
This expression is identical to that for the antidual EQS configuration, (4.6.18). We can conclude that the line currents should be located at a = (l2 - R2)1/2 and that the constant k used in that deduction (4.6.20) is identified using (26).

equation GIF #8.190
Thus, the surfaces of constant A_z are circular cylinders and represent the field lines shown in Fig. 8.4.8.

floating figure GIF #58
The inductance per unit length L is now deduced from (27).

equation GIF #8.191
In the limit where the conductors represent wires that are thin compared to their spacing, the inductance per unit length of (28) is approximated using (4.6.28).

equation GIF #8.192
Once the vector potential has been determined, it is possible to evaluate the distribution of current density on the conductors. Note that the currents tend to concentrate on the inside surfaces of the conductors, where the magnetic field intensity is more intense.

We are one step short of a general relationship between the capacitance per unit length and inductance per unit length of a pair of parallel perfect conductors, regardless of the cross-sectional geometry. With and A_z defined as zero on one of the conductors, evaluated on the other conductor they respectively represent the voltage and the flux linkage per unit length. Thus, with the understanding that and A_z are evaluated on the second conductor, L = A_z/i, (19), and C = ' /, (4.6.5). Here, i and are respectively the line current and line charge density that give rise to the same fields as do those actually on the surfaces of the conductors. These quantitites are related by (26), so we can conclude that regardless of the cross-sectional geometry, the product of the inductance per unit length and the capacitance per unit length is

equation GIF #8.193
where c is the velocity of light (3.1.16).

Note that inductance per unit length of parallel circular conductors given by (28) and the capacitance per unit length for the same conductors under ``open-circuit'' conditions (4.6.27) satisfy the general relation of (30).

Method of Images

In the presence of a planar perfect conductor, the zero normal flux condition can be satisfied by symmetrically mounting source distributions on both sides of the plane. This approach is familiar from Sec. 4.7, where the boundary condition required a plane of symmetry on which the tangential electric field was zero. Here, we require that the field intensity be tangential to the boundary. For two-dimensional configurations, the analogy between the electric potential and A_z makes the image method of Sec. 4.7 directly applicable here. In both cases, the symmetry plane is one of constant potential ( or A_z).

The most obvious example is an infinitely long line current at a distance d/2 from a perfectly conducting plane. If Fig. 4.7.1 were a picture of line charges rather than point charges, this would be the antidual situation. The appropriate image is then an oppositely directed line current located at a distance d/2 to the other side of the perfectly conducting plane. By making a pair of symmetrically located line currents the image for this pair of currents, the boundary condition on yet another plane can be satisfied, the analog to the configuration of Fig. 4.7.3. The following demonstration is intended to emphasize that the perfectly conducting symmetry plane carries a surface current that terminates the field in the region of interest. \begindemo2Surface Currents Induced in Ground Plane by Overhead Conductor The metal cylinder mounted over a metal ground-plane shown in Fig. 8.4.10 is familiar from Demonstration 4.7.1. Rather than being insulated from the ground plane and driven by a voltage source, this cylinder is shorted to the ground plane at one end and driven by a current source at the other. The height l is small compared to the length so that the two-dimensional model describes the field distribution in the mid-region.

floating figure GIF #59
A probe is used to measure the magnetic flux density tangential to the metal ground-plane. The distribution of this field, and hence of the surface current density in the adjacent metal, can be determined by recognizing that the ground plane boundary condition of no normal flux density is met by symmetrically mounting a distribution of oppositely directed currents below the metal sheet. This is just what was done in determining the fields for the pair of cylindrical conductors, Fig. 8.4.7. Thus, (25) is the image solution for the region x \geq 0. In terms of x and y,
equation GIF #8.194
The flux density tangential to the ground-plane is
equation GIF #8.195
Normalized to H_o = i/ a, this distribution is shown as a function of the probe position, Y, in the inset to Fig. 8.4.10. The role of the surface current density implied by this tangential field is demonstrated by the same probe measurement of the magnetic flux density normal to the conducting sheet. Provided that the frequency is high enough that the sheet does indeed behave as a ``perfect conductor'', this flux density is small compared to that tangential to the sheet. This is also true at the surface of the cylindrical conductor.

To appreciate the physical origins of this distribution, a dc current source is used in place of the ac source. Then, the distribution of current in the sheet is dictated by the rules of steady conduction as enunciated in the first half of Chap. 7. If the sheet is long enough compared to its width, the current is uniformly distributed over the sheet and over the cross-section of the cylinder. By contrast with the high-frequency a-c case where the field is terminated by surface currents in the sheet, the magnetic field now extends below the sheet.

The method of images is not restricted to the two-dimensional situations where there is a convenient analogy between and A_z. In the following example, of a three-dimensional field, the symmetry conditions are viewed without the aid of the vector potential.

Example 8.7.6. Current Loop Above a Perfectly Conducting Plane

A current loop with time varying current i is mounted a distance h above a perfectly conducting plane, as in Fig. 8.4.11. Its axis is inclined at an angle with respect to the normal to the plane. What is the net field produced by the current loop and the currents it induces in the plane?

floating figure GIF #60
To satisfy the boundary condition in the plane of the perfectly conducting sheet, an image loop is mounted as shown in Fig. 8.4.12. For each current segment in the actual loop, there is a segment in the image loop giving rise to an oppositely directed vertical component of H. Thus, the net normal flux density in the plane of the perfect conductor is zero.

8.8 The Boundary Value Point of View

The MQS configurations of this section are approached as were the EQS systems of Secs. 5.1 and 5.6. The volume V of interest is bounded by perfect conductors. There can also be a distribution of prescribed current density within V, perhaps by means of windings driven by current sources.

As in Secs. 5.1 and 5.6, the magnetic field intensity within V is divided into particular and homogeneous parts.

equation GIF #8.196
Although it is often easier to find the particular solution by a more direct method, the given current distribution can be used in the vector potential superposition integral to find a particular solution A_p that satisfies Ampère's law and the flux continuity condition at each point within V. Alternatively, H_p can be found directly from the Biot-Savart law, (8.3.7).

equation GIF #8.197
Note that this integration is confined to the volume V, where the distribution of current density is known.

Then, the homogeneous solution satisfies the source-free forms of Ampère's law and the flux continuity law.

equation GIF #8.198
equation GIF #8.199
In this chapter, we have emphasized the vector potential as a way of representing solenoidal fields. In Chap. 5, we used the scalar potential to represent irrotational fields. By design, the homogeneous field is both solenoidal and irrotational, so we now have a choice of representations, at least for the homogeneous part. The vector potential is convenient if the flux is to be found and if the analogies discussed in Sec. 8.4 are to be exploited. However, unless A has only one component, as it fortunately does in many important cases, it only replaces one vector with another in the field description. A scalar potential representation of H preserves the advantage of letting a scalar represent a vector in three dimensional configurations. Although we could now introduce the scalar potential, we will find that it is particularly useful in dealing with magnetization problems. Thus, the magnetic scalar potential will be introduced in Sec. 8.6 and used extensively in Chap. 9 and we proceed here to use A_h to represent the homogeneous solution. The flux continuity law is automatically satisfied by letting
equation GIF #8.200
and substitution into (3) shows that
equation GIF #8.201
which of course is the vector Poisson's equation, (8.2.3), with the current density set equal to zero.

Example 8.8.1. Inductive Attenuator

The cross-section of two conducting electrodes that extend to infinity in the z directions is shown in Fig. 8.5.1. The time-varying current in the +z direction in the electrode at y = b is returned in the -z direction through the \sqcup shaped electrode. This current is so rapidly varying that the electrodes behave as though they were perfectly conducting. The gaps of width insulating the electrodes from each other are small compared to the other dimensions of interest. The magnetic flux (per unit length in the z direction) passing through these gaps in the directions shown is defined as \Lambda (t).

floating figure GIF #61
The magnetic fields are two-dimensional and there are no sources in the region of interest. Thus, B can be represented in terms of A_z, which satisfies the z component of (6).

equation GIF #8.202
The walls are perfectly conducting in the sense that they are modeled as having no normal B (8.4.1). This means that A_z is constant on these walls (8.4.3). We define A_z to be zero on the vertical and bottom walls. Then, A_z must be equal to \Lambda on the upper electrode so that the flux per unit length in the z direction through the gaps is \Lambda. (Note that the z axis is out of the page, so that according to (8.1.7), \Lambda is positive if directed as shown in the figure.)
equation GIF #8.203
With the identification of variables
equation GIF #8.204
the boundary value problem is now formally identical to the EQS capacitive attenuator that was the theme of Sec. 5.5. Thus, it follows from (5.5.9) that
equation GIF #8.205
The lines of magnetic flux density are the lines of constant A_z. Thus, they are the equipotential ``lines'' of Fig. 5.5.3, shown in Fig. 8.5.1 with arrows added to indicate the field direction. Remember, there is a z-directed surface current density that is proportional to the tangential field intensity. For the flux lines shown, K_z is out of the page in the upper electrode and returned into the page on the side walls and (to an extent determined by b relative to a) on the bottom wall as well.

From the cross-sectional view given by Fig. 8.5.1, the provision for the current through the driven plate at the top to recirculate through the side and bottom plates is not shown. The following demonstration emphasises the implied current paths at the ends of the configuration.

Demonstration 8.8.1. Inductive Attenuator

One configuration described by Example 8.5.1 is shown in Fig. 8.5.2. Here, the upper plate is shorted to the adjacent walls at the near end and driven at the far end through a step-down transformer by a 20 kHz oscillator. The driving voltage v(t) at the far end of the upper plate is measured by means of an oscilloscope. The lower plate is shorted to the side-walls at the far end and also connected to these walls at the near end, but in such a way that the induced current i(t) can be measured by means of a current probe.

The walls and upper and lower plates are made from brass or copper. To insure that the resistances of the plate terminations are negligible, they are made from heavy copper wire with the connections soldered. (To make it possible to adjust the spacing b, braided wire is used for the shorts on the lower electrode.) If the length w of the plates in the z direction is large compared to a and b, H within the volume follows from (10). The surface current density K_z in the lower plate then follows from evaluation of the tangential H on its surface. In turn, the total current follows from integration of K_z over the width, a, of the plate.

equation GIF #8.206
With the objective of relating this current to the driving voltage, note that (8.4.14)
equation GIF #8.207
so that with the driving voltage a sinusoid of magnitude V,
equation GIF #8.208
Thus, in terms of the driving voltage, the output current is i_o sin ( t) where it follows from (11) and (13) that
equation GIF #8.209
We have found that the output current, normalized to U, has the dependence on spacing between upper and lower plates shown by the insert to Fig. 8.5.2. With the spacing b small compared to a, almost all of the current through the upper plate is returned in the lower one and the field between is essentially uniform. As the spacing b becomes comparable to the distance a between the side walls, most of the current through the upper electrode is returned in these side walls. Thus, for large b/a, the normalized output current of Fig. 8.5.2 reflects the exponential decay in the -y direction of the field.

floating figure GIF #62
Value is added to this demonstration if it is compared to its EQS antidual, Demonstration 5.5.1. For the EQS configuration, the lower plate was properly constrained to essentially the same potential as the walls by connecting it to these side walls through a resistance (which was then used to measure the induced current). Up to frequencies above 100 Hz in the EQS case, this resistance could be as high as that of the oscilloscope (say 1 M\Omega) and still constrain the lower plate to essentially the same zero potential as the walls. In the MQS case, we did not use a resistance to connect the lower plate to the side walls (and hence provide a means of measuring the output current) because that resistance would have had to have been extremely low, even at 20 kHz, to prevent flux from leaking through the gaps between the lower plate and the side walls. We used the current probe instead. The effects of finite conductivity in MQS systems are the subject of Chap. 10.

In a final example, we exemplify how the particular and homogeneous solutions are combined to satisfy boundary condtions while also illustrating how the inductance of a distributed winding is determined.

Example 8.8.2. Field and Inductance of Distributed Winding Bounded by Perfect Conductor

The cross-section of a distributed winding of radius a is shown in Fig. 8.5.3. It consists of turns carrying current i in the +z direction at a location (r, ) and returning the current at (r, - ) in the -z direction. The density of turns, each carrying the current i in the +z direction for 0 \leq and in the -z direction for < < 2, is

equation GIF #8.210
floating figure GIF #63
The total number of wires N in the left ball of the coil is
equation GIF #8.211
so that the current density is
equation GIF #8.212
The windings are very long in the z direction so that effects of the end-turns are ignored and the fields taken as independent of z.

The coil is bounded at r = a by a perfect conductor. With the following steps we determine the field distribution throughout the winding and finally its inductance.

The vector potential must satisfy (8.2.4) and is z-independent. In polar coordinates:

equation GIF #8.213
First, we look for a particular solution. If it is to take a product form, inspection shows that sin is the appropriate dependence. Substitution of an r dependence rn shows that the equation can be satisfied if n = 2. Thus, we have ``guessed'' a particular solution
equation GIF #8.214
The magnetic flux density normal to the perfectly conducting surface at r = a must be zero, so the total vector potential must be constant there. It follows that, at r = a, the homogeneous solution, A_zh must be the negative of the particular solution, A_zp.

equation GIF #8.215
A linear combination of the two solutions to Laplace's equation that have the same dependence as this condition is
equation GIF #8.216
The coefficient D must be zero so that the solution is finite at the origin. Then, the coefficient C is adjusted to make (21) satisfy the condition of (20). Hence, the sum of the particular and homogeneous solutions is
equation GIF #8.217
A graphical representation of what has been accomplished is given by Fig. 8.5.4, where the surfaces of constant A_z (and hence the lines of field intensity) are shown for the particular, homogeneous and total solutions.

floating figure GIF #64
Each turn of the coil links a different magnetic flux. Thus, to determine the total flux linked by the distribution of turns, it is necessary to carry out an integration. To do this, first observe that the flux linked by the turns with their right legs within the area rd dr in the neighborhood of (r, ) and their left legs within a similar area in the neighborhood of (r,- ) is
equation GIF #8.218
Here, l is the length of the system in the z direction.

The total flux linked by all of the turns is obtained by integrating over all of the turns.

equation GIF #8.219
Substitution for A_z from (22) then gives
equation GIF #8.220
where L will be recognized as the inductance.

Example 8.8.1. Vector Potential of a Line Current

The magnetic field intensity for the line current i of Fig. 8.6.3 was determined in Sec. 1.4 and in Example 8.3.1 using Ampère's integral law. The magnetic flux density is given by (1.4.10) as
equation GIF #8.221
floating figure GIF #65
The vector potential is determined by expressing curl A in polar coordinates.

equation GIF #8.222
Because there is no r component of the flux density, it follows that A_z is independent of . Setting the components of (6) and (7) equal and integrating gives
equation GIF #8.223
To illustrate the use of the vector potential to evaluate the net flux through a surface, consider the flux through the coil shown in Fig. 1.7.2a. The contour running parallel to the z axis in the z direction is at a radius d from the wire while that returning the current in the -z direction is at radius R. The length in the z direction is l. Thus, (5) gives a net flux of
equation GIF #8.224
Note that this result is in agreement with (1.7.5), which was obtained by evaluating the specific surface integral.