11
Energy,
Power Flow,
and Forces

11.0 Introduction

One way to decide whether a system is electroquasistatic or magnetoquasistatic is to consider the relative magnitudes of the electric and magnetic energy storages. The subject of this chapter therefore makes a natural transition from the quasistatic laws to the complete set of electrodynamic laws. In the order introduced in Chaps. 1 and 2, but now including polarization and magnetization,

1 For polarized and magnetized media at rest. these are Gauss' law [(6.2.1) and (6.2.3)]
equation GIF #11.1
Ampère's law (6.2.11),
equation GIF #11.2
Faraday's law (9.2.7),
equation GIF #11.3
and the magnetic flux continuity law (9.2.2).

equation GIF #11.4
Circuit theory describes the excitation of a two-terminal element in terms of the voltage v applied between the terminals and the current i into and out of the respective terminals. The power supplied through the terminal pair is vi. One objective in this chapter is to extend the concept of power flow in such a way that power is thought to flow throughout space, and is not associated only with current flow into and out of terminals. The basis for this extension is the laws of electrodynamics, (1)-(4).

Even if a system can be represented by a circuit, the need for the generalization of the circuit-theoretical power flow concept is apparent if we try to understand how electrical energy is transferred within, rather than between, circuit elements. The limitations of the circuit viewpoint would be crucial to testimony of an expert witness in litigation concerning the authority of the Federal Power Commission \footnote2Now the Federal Energy Regulatory Commission. to regulate power flowing between states. If the view is taken that passage of current across a border is a prerequisite for power flow, either of the devices shown in Fig. 11.0.1 might be installed at the border to "launder" the power. In the first, the state line passes through the air gap between capacitor plates, while in the second, it separates the primary from the secondary in a transformer.

\footnote3To be practical, the capacitor would be constructed with an enormous number of interspersed plates, so that in order to keep the state line in the air gap, a gerrymandered border would be required. Contemplation of the construction of a practical transformer, as described in Sec. 9.7, reveals that the state line would be even more difficult to explain in the MQS case. In each case, the current never leaves the state where it is generated. Yet in the examples shown, power generated in one state can surely be consumed in another, and a meaningful discussion of how this takes place must be based on a broadened view of power flow.

floating figure GIF #1
From the circuit-theoretical viewpoint, energy storage and rate of energy dissipation are assigned to circuit elements as a whole. Power flowing through a terminal pair is expressed as the product of a potential difference v between the terminals and the current i in one terminal and out of the other. Thus, the terminal voltage v and current i do provide a meaningful description of power flow into a surface S that encloses the circuit shown in Fig. 11.0.2. The surface S does not pass "inside" one of the elements.

floating figure GIF #2

Power Flow in a Circuit

For the circuit of Fig. 11.0.2, Kirchhoff's laws combine with the terminal relations for the capacitor, inductor, and resistor to give
equation GIF #11.5
equation GIF #11.6
Motivated by the objective to obtain a statement involving vi, we multiply the first of these laws by the terminal voltage v. To eliminate the term viL on the right, we also multiply the second equation by iL. Thus, with the addition of the two relations, we obtain
equation GIF #11.7
Because L and C are assumed to be constant, we can use the relation udu = d(fraction GIF #22 u2) to rewrite this expression as
equation GIF #11.8
where
equation GIF #11.9
With its origins solely in the circuit laws, (8) can be regarded as giving no more information than inherent in the original laws. However, it gives insights into the circuit dynamics that are harbingers of what can be expected from the more general statement to be derived in Sec. 11.1. These come from considering some extremes.

resistor is absent (G = 0), w is constant. Thus, the energy w is conserved in this limiting case. The solution to the circuit laws must lead to the conclusion that the sum of the electric energy fraction GIF #23 Cv2 and the magnetic energy 1\over 2 LiL2 is constant.

\bullet If the terminals are open (i = 0), and if the
\bullet Again, with G = 0, but now with a current supplied to the terminals, (8) becomes
equation GIF #11.10
Because the right-hand side is a perfect time derivative, the expression can be integrated to give
equation GIF #11.11
Regardless of the details of how the currents and voltage vary with time, the time integral of the power vi is solely a function of the initial and final total energies w. Thus, if w were zero to begin with and vi were positive, at some later time t, the total energy would be the positive value given by (10). To remove the total energy from the inductor and capacitor, vi must be reversed in sign until the integration has reduced w to zero. Because the process is reversible, we say that the energy w is stored in the capacitor and inductor.

resistor is present, (8) shows that the stored energy w must decrease with time. Because Gv2 is positive, this process is not reversible and we therefore say that the energy is dissipated in the resistor. In circuit theory terms, (8) is an example of an energy conservation theorem. According to this theorem, electrical energy is not conserved. Rather, of the electrical energy supplied to the circuit at the rate vi, part is stored in the capacitor and inductor and indeed conserved, and part is dissipated in the resistor. The energy supplied to the resistor is not conserved in electrical form. This energy is dissipated in heat and becomes a new kind of energy, thermal energy. Just as the circuit laws can be combined to describe the flow of power between the circuit elements, so Maxwell's equations are the basis for a field-theoretical view of power flow. The reasoning that casts the circuit laws into a power flow statement parallels that used in the next section to obtain the more general field-theoretical law, so it is worthwhile to review how the circuit laws are combined to obtain a statement describing power flow.

Overview

The energy conservation theorem derived in the next two sections will also not be a conservation theorem in the sense that electrical energy is conserved. Rather, in addition to accounting for the storage of energy, it will include conversion of energy into other forms as well. Indeed, one of the main reasons for our interest in power flow is the insight it gives into other subsystems of the physical world [e.g. the thermodynamic, chemical, or mechanical subsystems]. This will be evident from the topics of subsequent sections.

The conservation of energy statement assumes as many special forms as there are different constitutive laws. This is one reason for pausing with Sec. 11.1 to summarize the integral and differential forms of the conservation law, regardless of the particular application. We shall reference these expressions throughout the chapter. The derivation of Poynting's theorem, in the first part of Sec. 11.2, is motivated by the form of the general conservation theorem. As subsequent sections evolve, we shall also make continued reference to this law in its general form.

By specializing the materials to Ohmic conductors with linear polarization and magnetization constitutive laws, it is possible to make a clear identification of the origins of electrical energy storage and dissipation in media. Such systems are considered in Sec. 11.3, where the flow of power from source to "sinks" of thermal dissipation is illustrated. Processes of energy storage and dissipation are developed in greater depth in Secs. 11.4 and 11.5.

Through Sec. 11.5, the assumption is that materials are at rest. In Secs. 11.6 and 11.7, the power input is studied in the presence of motion of materials. These sections illustrate how the energy conservation law is used to determine electric and magnetic forces on macroscopic media. The discussion in these sections is confined to a determination of total forces. Consistent with the field theory point of view is the concept of a distributed force per unit volume, a force density. Rigorous derivations of macroscopic force densities are based on energy arguments paralleling those of Secs. 11.6 and 11.7.

In Sec. 11.8, we shall look at microscopic models of force density distributions that provide a picture of the origin of these distributions. Finally, Sec. 11.9 is an introduction to the macroscopic force densities needed to put electromechanical coupling on a continuum basis.

11.1 Integral and Differential Conservation Statements

The circuit with theoretical conservation theorem (11.0.8) equates the power flowing into the circuit to the rate of change of the energy stored and the rate of energy dissipation. In a field, theoretical generalization, the energy must be imagined distributed through space with an energy density W (joules/m3), and the power is dissipated at a local rate of dissipation per unit volume Pd (watts/m3). The power flows with a density S (watts/m2), a vector, so that the power crossing a surface Sa is given by \intSa S
da. With these field-theoretical generalizations, the power flowing into a volume V, enclosed by the surface S must be given by

boxed equation GIF #11.1
where the minus sign takes care of the fact that the term on the left is the power flowing into the volume.

According to the right-hand side of this equation, this input power is equal to the rate of increase of the total energy stored plus the power dissipation. The total energy is expressed as an integral over the volume of an energy density, W. Similarly, the total power dissipation is the integral over the volume of a power dissipation density Pd.

floating figure GIF #3
The volume is taken as being fixed, so the time derivative can be taken inside the volume integration on the right in (1). With the use of Gauss' theorem, the surface integral on the left is then converted to one over the volume and the term transferred to the right-hand side.
equation GIF #11.12
Because V is arbitrary, the integrand must be zero and a differential statement of energy conservation follows.

boxed equation GIF #11.2
With an appropriate definition of S, W and Pd, (1) and (3) could describe the flow, storage, and dissipation not only of electromagnetic energy, but of thermal, elastic, or fluid mechanical energy as well. In the next section we will use Maxwell's equations to determine these variables for an electromagnetic system.

11.2 Poynting's Theorem

The objective in this section is to derive a statement of energy conservation from Maxwell's equations in the form identified in Sec. 11.1. The conservation theorem includes the effects of both displacement current and of magnetic induction. The EQS and MQS limits, respectively, can be taken by neglecting those terms having their origins in the magnetic induction o (H + M)/ t on the one hand, and in the displacement current density (o E + P)/ t on the other.

Ampère's law, including the effects of polarization, is (11.0.2).

equation GIF #11.13
Faraday's law, including the effects of magnetization, is (11.0.3).

equation GIF #11.14
These field-theoretical laws play a role analogous to that of the circuit equations in the introductory section. What we do next is also analogous. For the circuit case, we form expressions that are quadratic in the dependent variables. Several considerations guide the following manipulations. One aim is to derive an expression involving power dissipation or conversion densities and time rates of change of energy storages. The power per unit volume imparted to the current density of unpaired charge follows directly from the Lorentz force law (at least in free space). The force on a particle of charge q is

equation GIF #11.15
The rate of work on the particle is
equation GIF #11.16
If the particle density is N and only one species of charged particles exists, then the rate of work per unit volume is
equation GIF #11.17
Thus, one must anticipate that an energy conservation law that applies to free space must contain the term Ju E. In order to obtain this term, one should dot multiply (1) by E.

A second consideration that motivates the form of the energy conservation law is the aim to obtain a perfect divergence of density of power flow. Dot multiplication of (1) by E generates ( x H) E. This term is made into a perfect divergence if one adds to it -( x E) H, i.e., if one subtracts (2) dot multiplied by H.

Indeed,

equation GIF #11.18
Thus, subtracting (2) dot multiplied by H from (1) dot multiplied by E one obtains
boxed equation GIF #11.3
In writing the first and third terms on the right, we have exploited the relation u du = d(fraction GIF #24 u2). These two terms now take the form of the energy storage term in the power theorem, (11.1.3). The desire to obtain expressions taking this form is a third consideration contributing to the choice of ways in which (1) and (2) were combined. We could have seen at the outset that dotting E with (1) and subtracting (2) after it had been dotted with H would result in terms on the right taking the desired form of "perfect" time derivatives.

In the electroquasistatic limit, the magnetic induction terms on the right in Faraday's law, (2), are neglected. It follows from the steps leading to (7) that in the EQS approximation, the third and fourth terms on the right of (7) are negligible. Similarly, in the magnetoquasistatic limit, the displacement current, the last two terms on the right in Ampère's law, (1), is neglected. This implies that for MQS systems, the first two terms on the right in (7) are negligible.

Systems Composed of Perfect Conductors and Free Space

Quasistatic examples in this category are the EQS systems of Chaps. 4 and 5 and the MQS systems of Chap. 8, where perfect conductors are surrounded by free space. Whether quasistatic or electrodynamic, in these configurations, P = 0, M = 0; and where there is a current density Ju, the perfect conductivity insures that E = 0. Thus, the second and last two terms on the right in (7) are zero. For perfect conductors surrounded by free space, the differential form of the power theorem becomes
boxed equation GIF #11.4
with
boxed equation GIF #11.5
and
boxed equation GIF #11.6
where S is the Poynting vector and W is the sum of the electric and magnetic energy densities. The electric and magnetic fields are confined to the free space regions. Thus, power flow and energy storage pictured in terms of these variables occur entirely in the free space regions.

Limiting cases governed by the EQS and MQS laws, respectively, are distinguished by having predominantly electric and magnetic energy densities. The following simple examples illustrate the application of the power theorem to two simple quasistatic situations. Applications of the theorem to electrodynamic systems will be taken up in Chap. 12.

Example 11.2.1. Plane Parallel Capacitor

The plane parallel capacitor of Fig. 11.2.1 is familiar from Example 3.3.1. The circular electrodes are perfectly conducting, while the region between the electrodes is free space. The system is driven by a voltage source distributed around the edges of the electrodes. Between the electrodes, the electric field is simply the voltage divided by the plate spacing (3.3.6),
equation GIF #11.19
while the magnetic field that follows from the integral form of Ampère's law is (3.3.10).

equation GIF #11.20
floating figure GIF #4
Consider the application of the integral version of (8) to the surface S enclosing the region between the electrodes in Fig. 11.2.1. First we determine the power flowing into the volume through this surface by evaluating the left-hand side of (8). The density of power flow follows from (11) and (12).

equation GIF #11.21
The top and bottom surfaces have normals perpendicular to this vector, so the only contribution comes from the surface at r = b. Because S is constant on that surface, the integration amounts to a multiplication.

equation GIF #11.22
where
equation GIF #11.23
Here the expression has been written as the rate of change of the energy stored in the capacitor. With E again given by (11), we double-check the expression for the time rate of change of energy storage.
equation GIF #11.24
From the field viewpoint, power flows into the volume through the surface at r = b and is stored in the form of electrical energy in the volume between the plates. In the quasistatic approximation used to evaluate the electric field, the magnetic energy storage is neglected at the outset because it is small compared to the electric energy storage. As a check on the implications of this approximation, consider the total magnetic energy storage. From (12),
equation GIF #11.25
Comparison of this expression with the electric energy storage found in (15) shows that the EQS approximation is valid provided that
equation GIF #11.26
For a sinusoidal excitation of frequency , this gives
equation GIF #11.27
where c is the free space velocity of light (3.1.16). The result is familiar from Example 3.3.1. The requirement that the propagation time b/c of an electromagnetic wave be short compared to a period 1/ is equivalent to the requirement that the magnetic energy storage be negligible compared to the electric energy storage.
A second example offers the opportunity to apply the integral version of (8) to a simple MQS system.

Example 11.2.2. Long Solenoidal Inductor

The perfectly conducting one-turn solenoid of Fig. 11.2.2 is familiar from Example 10.1.2. In terms of the terminal current i = Kd, the magnetic field intensity inside is (10.1.14),
equation GIF #11.28
while the electric field is the sum of the particular and conservative homogeneous parts [(10.1.15) for the particular part and Eh for the conservative part].

equation GIF #11.29
floating figure GIF #5
Consider how the power flow through the surface S of the volume enclosed by the coil is accounted for by the time rate of change of the energy stored. The Poynting flux implied by (19) and (20) is
equation GIF #11.30
This Poynting vector has no component normal to the top and bottom surfaces of the volume. On the surface at r = a, the first term in brackets is constant, so the integration on S amounts to a multiplication by the area. Because Eh is irrotational, the integral of Eh ds = E h rd around a contour at r = a must be zero. For this reason, there is no net contribution of Eh to the surface integral.
equation GIF #11.31
where
equation GIF #11.32
Here the result shows that the power flow is accounted for by the rate of change of the stored magnetic energy. Evaluation of the right hand side of (8), ignoring the electric energy storage, indeed gives the same result.

equation GIF #11.33
The validity of the quasistatic approximation is examined by comparing the magnetic energy storage to the neglected electric energy storage. Because we are only interested in an order of magnitude comparison and we know that the homogeneous solution is proportional to the particular solution (10.1.21), the latter can be approximated by the first term in (20).

equation GIF #11.34
We conclude that the MQS approximation is valid provided that the angular frequency is small compared to the time required for an electromagnetic wave to propagate the radius a of the solenoid and that this is equivalent to having an electric energy storage that is negligible compared to the magnetic energy storage.
equation GIF #11.35
A note of caution is in order. If the gap between the "sheet" terminals is made very small, the electric energy storage of the homogeneous part of the E field can become large. If it becomes comparable to the magnetic energy storage, the structure approaches the condition of resonance of the circuit consisting of the gap capacitance and solenoid inductance. In this limit, the MQS approximation breaks down. In practice, the electric energy stored in the gap would be dominated by that in the connecting plates, and the resonance could be described as the coupling of MQS and EQS systems as in Example 3.4.1.

In the following sections, we use (7) to study the storage and dissipation of energy in macroscopic media.

11.3 Ohmic Conductors With Linear Polarization andMagnetization

Consider a stationary material described by the constitutive laws
equation GIF #11.36
equation GIF #11.37
equation GIF #11.38
where the susceptibilities e and m, and hence the permittivity and permeability and , as well as the conductivity , are all independent of time. Expressed in terms of these constitutive laws for P and M, the polarization and magnetization terms in (11.2.7) become
equation GIF #11.39
equation GIF #11.40
Because these terms now appear in (11.2.7) as perfect time derivatives, it is clear that in a material having "linear" constitutive laws, energy is stored in the polarization and magnetization processes.

With the substitution of these terms into (11.2.7) and Ohm's law for Ju, a conservation law is obtained in the form discussed in Sec. 11.1. For an electrically and magnetically linear material that obeys Ohm's law, the integral and differential conservation laws are (11.1.1) and (11.1.3), respectively, with

boxed equation GIF #11.7
boxed equation GIF #11.8
boxed equation GIF #11.9
The power flux density S and the energy density W appear as in the free space conservation theorem of Sec. 11.2. The energy storage in the polarization and magnetization is included by simply replacing the free space permittivity and permeability by and , respectively.

The term Pd is always positive and seems to represent a rate of power loss from the electromagnetic system. That Pd indeed represents power converted to thermal form is motivated by considering the origins of the Ohmic conduction law. In terms of the bipolar conduction model introduced in Sec. 7.1, positive and negative carriers, respectively, experience the forces f+ and f-. These forces are balanced by collisions with the surrounding particles, and hence the work done by the field in forcing the migration of the particles is converted into thermal energy. If the velocity of the families of particles are, respectively, v+ and v-, and the number densities N+ and N-, respectively, then the rate of work performed on the carriers (per unit volume) is

equation GIF #11.41
In recognition of the balance between collision forces and electrical forces, the forces of (4) are replaced by |q+|E and -|q-|E, respectively.

equation GIF #11.42
If, in turn, the velocities are written as the products of the respective mobilities and the macroscopic electric field, (7.1.3), it follows that
equation GIF #11.43
where the definition of the conductivity (7.1.7) has been used.

The power dissipation density Pd = E E (watts/m3) represents a rate of energy loss from the electromagnetic system to the thermal system.

Example 11.3.1. The Poynting Vector of a Stationary Current Distribution

In Example 7.5.2, we studied the electric fields in and around a circular cylindrical conductor fed by a battery in parallel with a disk-shaped conductor. Here we determine the Poynting vector field and explore its spatial relationship to the dissipation density.

First, within the circular cylindrical conductor [region (b) in Fig. 11.3.1], the electric field was found to be uniform, (7.5.7),

equation GIF #11.44
while in the surrounding free space region, it was [from (7.5.11)]
equation GIF #11.45
and in the disk-shaped conductor [from (7.5.9)]
equation GIF #11.46
By symmetry, the magnetic field intensity is directed. The component of H is most easily evaluated from the integral form of Ampère's law. The current density in the circular conductor follows from (7) as Jo = v/L. Then,
equation GIF #11.47
equation GIF #11.48
The magnetic field distribution in the disk conductor is also deduced from Ampère's law. In this region, it is easiest to evaluate the r component of Ampère's differential law with the current density Jc = Ec, with Ec given by (9). Integration of this partial differential equation on z then gives a linear function of z plus an "integration constant" that is a function of r. The latter is determined by the requirement that H be continuous at z = -L.

equation GIF #11.49
It follows from these last four equations that the Poynting vector inside the circular cylindrical conductor, in the surrounding space, and in the disk-shaped electrode is
equation GIF #11.50
equation GIF #11.51
equation GIF #11.52
floating figure GIF #6
This distribution of S is sketched in Fig. 11.3.1. Wherever there is a dissipation density, there must be a negative divergence of S. Thus, in the conductors, the S lines terminate in the volume. In the free space region (a), S is solenoidal. Even with the fields perfectly stationary in time, the power is seen to flow through the open space to be absorbed in the volume where the dissipation takes place. The integral of the Poynting vector over the surface surrounding the inner conductor gives what we would expect either from the circuit point of view
equation GIF #11.53
where i is the total current through the cylinder, or from an evaluation of the right-hand side of the integral conservation law.

equation GIF #11.54
\sectnonumAn Alternative Conservation Theorem for Electroquasistatic Systems In describing electroquasistatic systems, it is inconvenient to require that the magnetic field intensity be evaluated. We consider now an alternative conservation theorem that is specialized to EQS systems. We will find an alternative expression for S that does not involve H. In the process of finding an alternative distribution of S, we illustrate the danger of ascribing meaning to S evaluated at a point, rather than integrated over a closed surface. In the EQS approximation, E is irrotational. Thus,
equation GIF #11.55
and the power input term on the left in the integral conservation law, (11.1.1), can be expressed as
equation GIF #11.56
Next, the vector identity
equation GIF #11.57
is used to write the right-hand side of (19) as
equation GIF #11.58
The first integral on the right is zero because the curl of a vector is divergence free and a field with no divergence has zero flux through a closed surface. Ampère's law can be used to eliminate curl H from the second.
equation GIF #11.59
In this way, we have determined an alternative expression for S, valid only in the electroquasistatic approximation.

boxed equation GIF #11.10
The density of power flow, expressed by (23) as the product of a potential and total current density consisting of the sum of the conduction and displacement current densities, has a form similar to that used in circuit theory. The power flux density of (23) is convenient in describing EQS systems, where the effects of magnetic induction are not significant. To be consistent with the EQS approximation, the conservation law must be used with the magnetic energy density neglected.

Example 11.3.2. Alternative EQS Power Flux Density for Stationary Current Distribution

To contrast the alternative EQS power flow density with the Poynting flux density, consider again the coaxial resistor configuration of Example 11.3.1. Because the fields are stationary, the EQS power flux density is
equation GIF #11.60
By contrast with the Poynting flux density, this vector field is zero in the free space region. In the circular cylindrical conductor, the potential and current density are [(7.5.6) and (7.5.7)]
equation GIF #11.61
and it follows that the power flux density is simply
equation GIF #11.62
There is a similar, radially directed flux density in the disk-shaped resistor. The alternative distribution of S, shown in Fig. 11.3.2, is clearly very different from that shown in Fig. 11.3.1 for the Poynting flux density.
floating figure GIF #7

Poynting Power Density Related to Circuit Power Input

Suppose that the surface S described by the conservation theorem encloses a system that is accessed through terminal pairs, as shown in Fig. 11.3.3. Under what circumstances is the integral of S da over S equivalent to summing the voltage-current product of the terminals of the wires connected to the system?

floating figure GIF #8
Two attributes of the fields on the surface S enclosing the system are required. First, the contribution of the magnetic induction to E must be negligible on S. If this is so, then regardless of what is inside S (for example, both EQS and MQS systems), on the surface S, the electric field can be taken as irrotational. It follows that in taking the integral over a closed surface of the Poynting power density, we can just as well use (23).
equation GIF #11.63
By contrast with the EQS systems treated in deriving this expression, it now holds only on the surface S, not necessarily on surfaces inside the volume enclosed by S. Second, on the surface S, the contribution of the displacement current must be negligible. This is equivalent to requiring that S is chosen parallel to the displacement flux density. In this case, the total power into the system reduces to
equation GIF #11.64
The integrand has value only where the surface S intersects a wire. If taken as perfectly conducting (but nevertheless in a region where B/ t is zero and hence E is irrotational), the wires have potentials that are uniform over their cross-sections. Thus, in (28), is equal to the voltage of the terminal. In integrating the current density over the cross-section of the wire, note that da is directed out of the surface, while a positive terminal current is directed into the surface. Thus,

equation GIF #11.65
and the input power expressed by (28) is equivalent to what would be expected from circuit theory.

Poynting Flux and Electromagnetic Radiation

Power cannot be supplied to or lost by a quasistatic system of finite extent through a surface at infinity. Such a power supply or loss requires radiation, and electromagnetic waves are neglected when either the magnetic induction or the displacement current density are neglected. To prove this statement, consider an EQS system of finite net charge. Its electric field intensity decays like 1/r2 at infinity, where r is the distance to a far-off point from some origin chosen within the system. At a great distance, the currents appear equivalent to current loop sources. Hence, the magnetic field intensity has the 1/r3 decay typical of a magnetic dipole. It follows that the Poynting vector decays at least as fast as 1/r5, so that the flux of E x H integrated over the "sphere" at infinity of area 4 r2 gives zero contribution. Because it is only that part of E x H resulting from electromagnetic radiation that contributes at infinity, Poynting's theorem is shown in Sec. 12.5 to be a powerful tool for dealing with antennae.

11.4 Energy Storage

In the conservation theorem, (11.2.7), we have identified the terms E P/ t and H \cdot o M / t as the rate of energy supplied per unit volume to the polarization and magnetization of the material. For a linear isotropic material, we found that these terms can be written as derivatives of energy density functions. In this section, we seek a more general description of energy storage. First, nonlinear materials are considered from the field viewpoint. Then, for those systems that can be described in terms of electrical terminal pairs, energy storage is formulated in terms of terminal variables. We will find the results of this section directly applicable to finding electric and magnetic forces in Secs. 11.6 and 11.7.

floating figure GIF #9

Energy Densities

Consider a material in which E and D \equiv (o E + P) are collinear. With E and D representing the magnitudes of these vectors, this material is presumed to be described by a constitutive law in which E is a single-valued function of D, such as that sketched in Fig. 11.4.1a. In the case of a linear constitutive law, the curve is a straight line with a slope equal to the permittivity . Consider a material in which E and P are collinear (isotropic material). Then, of course, E and D \equiv o E + P are collinear as well. One may graph the magnitude of D versus E and obtain a complete characterization of the material. Now the power per unit volume imparted to the polarization is E P/ t. If one adds to it the rate of energy supply to the field per unit volume (the free space part) E oE/ t, one obtains for the power per unit volume
boxed equation GIF #11.11
The power supplied to the unit volume can now be written as the time derivative of a function of D, We(D). Indeed, if we define the area above the graph in Fig. 11.4.1 as We, then
boxed equation GIF #11.12
Thus, E ( D/ t) is the derivative of the function We(D). This function is the energy stored per unit volume, because the energy supplied per unit volume expressed by the integral

equation GIF #11.66
is a function of the final value D of the displacement flux, and we assumed that the fields E and D were zero at t = -. Here, D represents the differential of D, usually denoted by dD. We will use rather than d to avoid confusion between differentials used in carrying out volume, surface and line integrals and the differential used here, which implies an integration in a "state space" having the "dimension" D.

Similar arguments show that if B \equiv o (H + M) and H are collinear, and if H is a single-valued function of B, then

boxed equation GIF #11.13
where
boxed equation GIF #11.14
With (1) and (4) replacing the first four terms on the right in the energy theorem of (11.2.7), it is clear that the energy density W = We + Wm. The electric and magnetic energy densities have the geometric interpretations as areas on the graphs representing the constitutive laws in Fig. 11.4.1.

Energy Storage in Terms of Terminal Variables

It was shown in Sec. 11.3 that the power input to a system could be represented by the sum of the vi products for each of the terminal pairs, (11.3.29), provided certain conditions were met in the neighborhoods of the terminals. The description of energy storage in a loss-free system in terms of terminal variables will be found useful in determining electric and magnetic forces. With the assumption that all of the power input to a system is accounted for by a time rate of change of the energy stored, the energy conservation statement for a system becomes

equation GIF #11.67
where
equation GIF #11.68
and the integral is carried over the volume of the system. If the system is electroquasistatic, conservation of charge requires that the terminal current be the time rate of change of the charge on the electrode to which the positive terminal is attached.

boxed equation GIF #11.15
Further, w = we, the stored electric energy. Thus, one concludes from (6) that
equation GIF #11.69
The second expression states that with the addition of an incremental amount of charge dqi to an electrode having the voltage vi goes an incremental change in the stored energy we. Integration on the charges then gives the total energy
boxed equation GIF #11.16
floating figure GIF #10
To complete this integral, each of the terminal voltages must be a known function of the associated charges.

boxed equation GIF #11.17
Integration is then carried out along any path in the state space (q1 \ldots qn) that begins at the origin and ends with the desired charges on the electrodes (and hence the desired terminal voltages). For a single terminal pair, the energy can be pictured as the area shown in Fig. 11.4.2a.

If the system is magnetoquasistatic, the conservation law for a lossless system that can be described by terminal relations again takes the form of (6). However, rather than expressing the currents as derivatives of electrode charges, the voltages are derivatives of the fluxes linked by the respective terminal pairs.

boxed equation GIF #11.18
Then, (6) leads to
boxed equation GIF #11.19
To complete this integral, we require the terminal currents as functions of the terminal flux linkages.

boxed equation GIF #11.20
For a single terminal pair system, wm is portrayed in Fig. 11.4.2b.

The most general way to compute the total energy stored in a system is to integrate the energy densities given by (3) and (5) over the volumes of the respective systems. If systems can be described in terms of terminal relations and are loss free, (9) and (12) must lead to the same answers. Note that (D, E) and (q, v) are the field and circuit variables in the EQS systems, while (B, H) and ( ,i) have corresponding roles in MQS systems.

Example 11.4.1. An Electrically Linear System

A dielectric slab of permittivity partially fills the region between plane parallel perfectly conducting electrodes, as shown in Fig. 11.4.3. With the fringing field ignored, we find the total energy stored by two methods. First, the energy density is integrated over the volume. Then, the terminal relation is used to evaluate the total energy.

floating figure GIF #11
An exact solution for the electric field well between the electrodes is simply E =ix (v/a). Note that this field satisfies the boundary conditions at the interface between the dielectric slab and the free space region above and at the electrodes. We assume that a \ll b and therefore neglect the fringing fields.

The energy density in the linear dielectric, where D = E, follows from evaluation of (3).

equation GIF #11.70
equation GIF #11.71
In the free space region, the same result applies with o. Integration of these energy densities over the regions in which they apply amounts to a multiplication by the respective volumes. Thus, the total energy is
equation GIF #11.72
Note that this expression takes the form
equation GIF #11.73
where
equation GIF #11.74
In terms of the terminal variables, where q = Cv, the total energy follows from an evaluation of (9).

equation GIF #11.75
Once the integration has been carried out, the last expression is written by again using the relation q = Cv. Note that the volume integration of the energy density and the integration in terms of the terminal variables give the same result.

The next example considers an MQS system with two terminal pairs and thus illustrates the integration called for in evaluating the energy from the terminal relations. Also, the energy stored in coupled inductors is often of practical interest.

Example 11.4.2. Coupled Coils; Transformers

An example of a two terminal pair lossless MQS system is a pair of coupled coils having the terminal relations
equation GIF #11.76
In this case, (12) becomes
equation GIF #11.77
To evaluate this expression, we need to substitute for the currents written in terms of the flux linkages. This requires the inversion of (18). For linear systems, this is easily done, but not for nonlinear systems. To avoid inversion, we rewrite the right-hand side of (19), which becomes
equation GIF #11.78
and regroup terms
boxed equation GIF #11.21
where the coenergy is defined as
boxed equation GIF #11.22
Equation (21) can be integrated when the flux linkages are expressed in terms of the currents, and that is the form in which the terminal relations are given by (18). Once the coenergy wm' has been found, wm follows from (22).

floating figure GIF #12
The integration of (19) is a line integral in a state space (i1, i2). If energy is conserved, we must be able to carry out this integration along any path that begins with the currents turned off and ends with the currents at the desired values. In the path represented by Fig. 11.4.4, the current i1 is turned up first while holding the current i2 to zero. Then, with i1 held fixed at its final value, the current i2 is raised from zero to its final value. For this path, the integration of (22) becomes
equation GIF #11.79
Substitution for the flux linkages from (18) and evaluation of the integrals then gives
equation GIF #11.80
If the integration is carried out along a path where the roles of i1 and i2 are reversed, the expression obtained is (24) with L12 L21. To make the energy stored independent of path, the mutual inductances must be equal.

equation GIF #11.81
This relation, which we found to hold for the transformer of Example 9.7.4, is required if energy is to be conserved. The energy is now evaluated by substituting this expression and the flux linkages expressed using (18) into (22) solved for wm. It follows that

equation GIF #11.82
Evaluation of the energy stored in a unity-coupled transformer, where the inductances take the form of (9.7.20), gives
equation GIF #11.83
Operating under "ideal" conditions [in the sense that i2/i1 = -N1/N2, (9.7.13)], the transformer does not store energy, wm = 0. Thus, according to the power theorem in the form of (6), under ideal operating conditions, the power input at one terminal pair instantaneously appears as a power output at the second terminal pair.

Examples have so far involved linear polarization and magnetization constitutive laws. In the following, the EQS energy storage in a material having a nonlinear polarization constitutive law is determined.

Example 11.4.3. Energy Storage in Electrically Nonlinear Material

To represent the tendency of the polarization to saturate as the electric field is raised, a constitutive law might take the form
equation GIF #11.84
Here, 1 and 2 are parameters descriptive of the specific material, and D is collinear with E. This constitutive law is portrayed graphically in Fig. 11.4.5.

Because D is given as a function of E that is not easily solved for E as a function of D, the computation of the electric energy density using (3) is inconvenient. However, we can observe that

equation GIF #11.85
and then regroup terms so that the expression becomes
boxed equation GIF #11.23
where
boxed equation GIF #11.24
Integration now leads to the coenergy density We', but the energy density We can then be found using (31) and the constitutive law.

Specifically, evaluation of (30) using (28) gives the coenergy density

equation GIF #11.86
It follows from (31) that the energy density is
equation GIF #11.87
floating figure GIF #13
A graphical representation of the energy and coenergy functions is given in Fig. 11.4.5. The area "under the curve" with D as the integration variable is We, (3), and the area under the curve with E as the integration variable is We', (31).

11.5 Electromagnetic Dissipation

The heat generated by electromagnetic fields is often the controlling feature of an engineering design. Semiconductors inevitably produce heat, and the distribution and magnitude of the heat source is an important consideration whether the application is to computers or power conversion. Often, the generation of heat poses a fundamental limitation on the performance of equipment. Examples where the generation of heat is desirable include the heating coil of an electric stove and the microwave irradiation of food in a microwave oven. Ohmic conduction is the primary cause of heat generation in metals, but it also operates in semiconductors, electrolytes, and (at low frequencies) in semi-insulating liquids and solids. The mechanism responsible for this type of heating was discussed in Sec. 11.3. The dissipation density associated with Ohmic conduction is E E.

An Ohmic current can be imposed by making electrical contact with the material, as for the heating element in a stove. If the material is a good conductor, such currents can also be induced by magnetic induction (without electrical contact). The currents induced by time-varying magnetic fields in Chap. 10 are an example. Induction heating is an MQS process and often used in processing metals. Currents induced in transformer cores by the time-varying magnetic flux are an example of undesirable heating. In this context, the associated losses (which are minimized by laminating the core) are said to be due to eddy currents.

Ohmic heating can also be induced by "capacitive" coupling. In the EQS examples of Sec. 7.9, dielectric heating is caused by the currents associated with the accumulation of unpaired charges.

Whether due to magnetic induction or capacitive coupling, the generation of heat is described by the dissipation density Pd = E E identified in Sec. 11.3. However, the polarization and magnetization terms in the conservation theorem, (11.2.7), can also be responsible for energy dissipation. This occurs when the (electric or magnetic) dipoles do not align instantaneously with the fields. The polarization and magnetization constitutive laws differ from the laws postulated in Sec. 11.3.

As an example suggesting how the polarization term in (11.2.7) can represent dissipation, picture the artificial dielectric of Demonstration 6.6.1 (the ping-pong ball dielectric) but with spheres that are highly resistive rather than perfectly conducting. The accumulation of charge on the poles of the spheres in response to the application of an electric field is described by a rate, rather than a magnitude, that is proportional to the field. Thus, we would expect P/ t rather than P to be proportional to E. With a coefficient representing the properties and geometry of the spheres, the polarization constitutive law would then take the form

equation GIF #11.88
If this law is used to express the polarization term in the conservation law, the second term on the right in (11.2.7), a positive definite quantity results.

equation GIF #11.89
As might be expected from the physical origins of the constitutive law, the polarization term now represents dissipation rather than energy storage.

When materials are placed in electric fields having frequencies so high that conduction effects are negligible, losses due to the polarization of dipoles become the dominant heating mechanism. The artificial diamagnetic material considered in Demonstration 9.5.1 suggests how analogous losses are associated with the dynamic magnetization of a material. If the spherical particles comprising the artificial diamagnetic material have a finite conductivity, the induced dipole moments are not in phase with an applied sinusoidal field. What amounts to Ohmic dissipation on the particle scale is accounted for on the macroscopic scale by a modified constitutive law of magnetization.

The most common losses due to magnetization are encountered in ferromagnetic materials. Hysteresis losses occur because of the coercion required to obtain alignment of ferromagnetic domains. We will end this section with the relationship between the hysteresis curve of Fig. 9.4.6 and the dissipation density.

Energy Conservation for Temporally Periodic Systems

Many practical situations involve fields that vary with time in a periodic fashion. The sinusoidal steady state is the most common example. If the energy conservation law (11.0.8) is integrated over one period T, the energy storage term makes no contribution.

equation GIF #11.90
As a result, the time average of the conservation law states that the time average of the input power goes into the time average of the dissipation. The time average of the integral form of the conservation law, (11.1.1), becomes
equation GIF #11.91
This expression, which assumes that the dynamics are periodic but not necessarily sinusoidal, gives us two ways to compute the total energy dissipation. Either we can use the right-hand side and integrate the power dissipation density over the volume, or we can use the left-hand side and integrate the time average of S da over the surface enclosing the volume. Consider the sinusoidal steady state as a particular case. If P and M are related to E and H by linear differential equations, an approach can be taken that is familiar from circuit theory. The phase and amplitude of each field at a given location are represented by a complex amplitude. For example, the electric and magnetic field intensities are written as
equation GIF #11.92
A complex vector \hatE (r) has three complex scalar components \hat Ex (r), \hat Ey (r), and \hat Ez(r). The meaning of each is the same as the meaning of a complex voltage in circuit theory: e.g., the magnitude of \hat Ex(r), |\hat Ex(r)|, gives the peak amplitude of the x component of the electric field varying cosinusoidally with time, and the phase of \hat Ex(r) gives the phase advance of the cosine time function.

In determining the time averages of products of quantities that are in the sinusoidal steady state, it is helpful to make use of the time average theorem. With * designating the complex conjugate,

boxed equation GIF #11.25
This can be shown by using the identity
equation GIF #11.93

Induction Heating

In this case, the heating is represented by Ohmic conduction and Pd given by (11.3.3c). The examples from Chaps. 7 and 10 involving conductors of finite conductivity offer the opportunity to apply this relation to the evaluation of the right-hand side of (4). If the same total time average power is calculated using the left-hand side of this expression, it may seem that Ohm's law is not required. However, remember that this law is also reflected in the field quantities used to calculate S.

Example 11.5.1. Induction Heating of the Thin Shell

The thin conducting shell of Fig. 11.5.1, in a field Ho(t) applied collinear with its axis, was described in Example 10.3.1. Here the applied field is in the sinusoidal steady state

equation GIF #11.94
floating figure GIF #14
According to (10.3.9), the complex amplitude of the response, the magnetic field inside the shell, is
equation GIF #11.95
where m = fraction GIF #25 o a.

The complex amplitude of the surface current density circulating in the shell follows from (10.3.8).

equation GIF #11.96
Because the current density is uniform over the radial cross-section of the shell, the dissipation density can be written in terms of the surface current density K = E.

equation GIF #11.97
It follows from the application of the time average theorem, (6), that the total time average dissipation is
equation GIF #11.98
where l is the shell length. To complete the derivation based on an integration of the density over the volume of the conductor, this expression can be evaluated using (10).

equation GIF #11.99
The same result is found by evaluating the time average of the Poynting flux density integrated over a surface that is just outside the shell at r = a. To see this, we again use the time average theorem, (6), and recognize that the surface integral amounts to a multiplication by the surface area of the shell.
equation GIF #11.100
To evaluate this expression, (10) is used to determine E
equation GIF #11.101
Evaluation of (14) then gives
equation GIF #11.102
which is the same result as found by integrating the dissipation density over the volume, (13).

The dependence of the time average power dissipation on the normalized frequency is shown in Fig. 11.5.2. At very low frequencies, the induced current is not large enough to have an appreciable effect on the imposed field. Thus, the electric field is proportional to the time rate of change of the applied field, and because the dissipation is proportional to the square of E, the power dissipation increases as the square of . At high frequencies, the induced current can be no more than that required to shield the imposed field from the region inside the shell. As a result, the dissipation reaches an asymptotic limit.

floating figure GIF #15
Which of the two approaches is best for finding the total power dissipation? The answer depends on what field information is available. Certainly, the notion that the total heat generated can be found by integrating over a surface that is completely outside the heated material is a fundamental consequence of Poynting's theorem.

Dielectric Heating

In the sinusoidal steady state, we can identify the power dissipation density associated with polarization by finding the time average
equation GIF #11.103
In view of the time average theorem, (6), this becomes
equation GIF #11.104
If the polarization P does not follow the electric field E instantaneously, yet the material is still linear and isotropic, the complex vector \hat P can be related to \hatE by a complex susceptibility. Or, instead, the complex displacement flux density vector \hatD is related to \hatE by a complex dielectric constant.
equation GIF #11.105
Here \hat is the complex permittivity with real and imaginary parts ' and -", respectively. Evaluation of (18) using this constitutive law gives
equation GIF #11.106
Thus, " represents the electrical dissipation associated with the polarization process.
floating figure GIF #16
In the literature, the loss tangent \tan is often used to represent dissipation. It is the tangent of the phase angle of the complex dielectric constant defined in terms ' and " in Fig. 11.5.3. Thus,
equation GIF #11.107
From this definition, it follows from Eulers formula that
equation GIF #11.108
Given the complex amplitude of the electric field, D is
equation GIF #11.109
If the electric field is Eo cos ( t), then D is |\hat | Eo cos ( t - ). The electric displacement lags the electric field by the phase angle .

In terms of the loss tangent defined by (21), the time average electrical dissipation density of (20) becomes

equation GIF #11.110
Usually the loss tangent and ' are measured. In the following example, we compute the complex permittivity from a model of the polarizable medium and find the electrical dissipation on a macroscopic basis. In this special case we have the option of finding the time average loss by considering each of the dipoles on a microscopic basis. This is not generally possible, because the interactions among dipoles that are neglected in this example are usually too complicated for an analytic treatment.

Example 11.5.2. An Artificial Lossy Dielectric

By putting together examples considered in Chaps. 6 and 7, we can illustrate the origins of the complex permittivity. The artificial dielectric of Example 6.6.1 and Demonstration 6.6.1 had "molecules" consisting of perfectly conducting spheres. As a result, the polarization was pictured as instantaneously in step with the applied field. We consider now the result of having spheres that have finite conductivity.

The response of a single sphere having a finite conductivity and permittivity surrounded by free space is a special case of Example 7.9.3. The response to a sinusoidal drive is summarized by (7.9.36), where we set a = 0, a = o, b = , and b = . All that is required from this solution for the potential is the moment of a dipole that would give rise to the same exterior field as does the sphere. Comparison of the potential of a dipole, (4.4.10), to that given by (7.9.36a) shows that the complex amplitude of the moment is

equation GIF #11.111
where e \equiv (2o + )/. If mutual interactions between dipoles are ignored, the polarization density P is this moment of a single dipole multiplied by the number of dipoles per unit volume, N. For a cubic array with a distance s between the dipoles (the centers of the spheres), N = 1/s3. Thus, the complex amplitude of the electric displacement is
equation GIF #11.112
Combining this result with the moment given by (25) yields the desired constitutive law in the form \hat D = \hat \hat E, where the complex permittivity is
equation GIF #11.113
The time average power dissipation density follows from this expression and (20).

equation GIF #11.114
The dependence of the power dissipation on frequency has the same form as for the induction heating example, Fig. 11.5.2. At low frequencies, the surface charges induced at the north and south poles of each sphere are completely determined by the external field. Thus, the current density within the sphere that makes possible the accumulation of these surface charges is proportional to the time rate of change of the applied field. At low frequencies, the dissipation is proportional to the square of the volume current and hence to the square of the time rate of change of the applied field. As a result, at low frequencies, the dissipation density increases with the square of the frequency. As the frequency is raised, less surface charge is induced on the spheres. Although the amount of charge induced is inversely proportional to the frequency, there is a compensating effect because the volume currents are responsible for the dissipation, and these are proportional to the time rate of change of the charge. Thus, the dissipation density reaches a saturation value as the frequency becomes very high.

One tool used to form a picture of atomic, molecular, and domain physics is dielectric spectroscopy. Using this approach, the frequency dependence of the complex permittivity is used to gain insight into the microscopic structure.

Magnetization, like polarization, can also be the source of dissipation. The time average dissipation density due to magnetization follows by taking the time averge of the third and fourth terms on the right in the basic power theorem, (11.2.7). Combined, these terms give

equation GIF #11.115
For small-signal applications, this source of dissipation is dealt with by introducing a complex permeability \hat such that \hat B = \hat \hat H. The role of the complex permeability is similar to that of the complex permittivity. The artificial diamagnetic material of Example 9.5.2 and Demonstration 9.5.1 can be used to exemplify the concept. Instead of perfectly conducting spheres that give rise to a magnetic moment instantaneously induced antiparallel to the applied field, spherical shells of finite conductivity would be used. The dipole moment induced in the individual spherical shells would be deduced following the same approach as in Sec. 10.4. The resulting dipole moment would not be in phase with an applied sinusoidally varying magnetic field. The derivation of an equivalent complex permeability would follow from the same line of reasoning as used in the previous example.

Hysteresis Losses

Under periodic conditions in magnetizable solids, B and H are related by the hysteresis curve described in Sec. 9.4 and illustrated again in Fig. 11.5.4. What time average power dissipation is implied by the hysteresis? As before, B and H are collinear. However, neither is now a single-valued function of the other. Evaluation of (29) is accomplished by breaking the cycle into two parts, each involving a single-valued relationship between B and H. The first is the upswing "trajectory" from A C in Fig. 11.5.4. Over this half-cycle, which takes B from BA to BC, the trajectory is H+(B). With B taken as BA when t = 0, it follows from (11.4.4) and (11.4.5) that
equation GIF #11.116
This is the area under the curve of H versus B between A and C in Fig. 11.5.4, traversed on the "upswing." A similar evaluation for the "downswing," where the trajectory is H-(B), gives
equation GIF #11.117
The time average power dissipation, (29), then is the sum of these two contributions divided by T.

equation GIF #11.118
Thus, the area within the hysteresis loop is the energy dissipated in one cycle.

floating figure GIF #17

11.6 Electrical Forces on Macroscopic Media

Electrical forces on macroscopic materials have their origins in the forces exerted on the microscopic particles of which the materials are composed. Macroscopic fields have been used to describe conduction, polarization, and magnetization. In Chaps. 6, 7 and 9, polarization, current, and magnetization densities, respectively, were related to the macroscopic field variables through constitutive laws. Typically, the parameters in these laws are determined from measurements. Thus, the experimentally determined relations make it unnecessary to take detailed account of how the microscopic fields are averaged. Because the definition of the average is already implicit in our macroscopic formulation of Maxwell's equations, we must now take care that our use of macroscopic field quantities for representing electromagnetic forces is self-consistent. The force on a macroscopic volume element V of a material is the sum of the forces on the charged particles and magnetic dipoles constituting the material. Consider the simple case in which no magnetic dipoles are present. Then
equation GIF #11.119
where the summation is over all the charges within V at their respective positions. Now, the fields E(ri) and H(ri) are the microscopic fields that vary greatly from point ri to point rj in the material. The macroscopic fields E(r) and H(r) are averaged (smoothed) versions of these fields, whose sources are the averaged charge densities
equation GIF #11.120
and
equation GIF #11.121
where the velocity vi of the microscopic particles should be distinguished from that of the macroscopic material in which they are embedded or through which they move. The average of a product is not equal to the product of the averages. Thus, one could not find the force density F = f/ V from the expression E + J x o H, as the product of the averaged charge density and averaged electric field plus averaged current density times averaged magnetic flux density. Other methods have to be used to determine the force. One of the most useful is the energy method. Given the constitutive law for the material, which represents the interrelationship between macroscopic field variables, conservation of energy provides a way of deducing the self-consistent force acting on the material.

In this and the next section, we illustrate how total forces can be determined using conservation of energy as a premise. In this section, the EQS systems considered have only one mechanical degree of freedom and only one electrical terminal pair. In the next section, MQS systems are considered and the approach is broadened to a somewhat more general class of systems. A parallel approach determines the force density rather than the total force. After expanding on microscopic forces in Sec. 11.8, we shall review macroscopic force densities in Sec. 11.9.

floating figure GIF #18
Typical of the electroquasistatic problems considered in this section is the pair of metallic electrodes shown in Fig. 11.6.1. With the application of a voltage, unpaired charges of opposite polarity are induced on the electrode surfaces. The electrical state of the system is specified by giving the geometry and the potential difference v between the electrodes. Here we picture one electrode as movable, with its position denoted by . The two terminal pair system of Fig. 11.6.1b is useful to include mechanical effects via an additional terminal pair. If we think of the net unpaired charge q on the electrode as an electrical terminal variable complementing v, then the force of electrical origin f complements the mechanical displacement .

Given the electrical terminal relation v = v(q, ), we now use an energy conservation principle to determine the force f = f(q, ) that acts to increase the displacement . The electrical terminal relation can either be regarded as a measured function or be predicted using the macroscopic field laws and constitutive laws for the materials within the "box." It is now assumed that there is no conversion of electrical energy to thermal form within the box of Fig. 11.6.1b. Mechanisms for conversion of energy to heat are modeled by elements outside the box. For example, the finite conductivity of any dielectric is taken into account by a resistance external to the system. Thus, the electrical power input to what is defined as the "box," the electroquasistatic subsystem, must either result in a change in the electrical energy stored or mechanical power expended as the force f acts on the mechanical system. The integral form of the power conservation theorem, (11.1.1), is generalized to include the rate of work by the force f

equation GIF #11.122
In Sec. 11.4, we represented the quasistatic net electrical power input on the left in this expression in circuit theory terms. With the total energy we defined as the integral of the energy density over the entire volume of the system, (4) becomes
equation GIF #11.123
where the electrical power input is the product vi = vdq/dt. Multiplication of (4) by dt converts a statement of power flow to one of energy conservation.

boxed equation GIF #11.26
If an increment of charge dq is placed on an electrode at potential v, an increment of energy vdq is added to the system that produces a change in the total stored energy dwe, an increment of work fd done on an external mechanical system, or some combination of both. Here, f(q, ) is the as yet unknown force. Solved for dwe, this energy conservation statement is
equation GIF #11.124
This expression describes what might be termed a quasistatic electrical and mechanical subsystem. The state of this subsystem is specified by prescribing the geometry () and the charge on the electrode, for then the voltage of the electrode follows from the terminal relation v(q, ). The state of the subsystem is fully determined by the variables (q, ), which are therefore regarded as independent variables. In terms of the two terminal pairs shown in Fig. 11.6.1b, one of each pair of terminal variables has been chosen as an independent variable.

The incremental change in we (q, ) associated with incremental changes of dq and d in the independent variables is

equation GIF #11.125
Because q and can be independently specified, (7) and (8) must hold for any combination of dq and d. For example, they must hold if the position of the electrode is held fixed so that d = 0 and the charge is changed by the incremental amount dq. They must also describe the change in energy resulting from making an incremental displacement d of the electrode under open circuit conditions, where dq = 0. Indeed, (7) and (8) hold if q and are changed by arbitrary incremental amounts, and so it follows that the coefficients of dq in (7) and (8) must be equal to each other, as must the coefficients of d.

boxed equation GIF #11.27
Given the total energy, written in terms of the independent variables (q, ), the second of these relations provides the desired force. Integration of the energy density over the volume of the system is one way to determine we. Another is to integrate (7) along a line in the state space (q, ) designed so that the integral can be carried out without having to know f.
equation GIF #11.126
floating figure GIF #19
Such a path

4 Note the analogy with the line integral \int (Ex dx + Ey dy) of a two-dimensional conservation field that results in the potential (x, y). is shown in Fig. 11.6.2, where it is assumed that the force of electrical origin f is zero if the charge q is zero. Thus, in integrating along the contour q = 0 from A B, dq = 0 and f = 0, so there is no contribution. The remainder of the integral, from B C, is carried out with fixed, so d = 0, and (10) reduces to
boxed equation GIF #11.28
We have accounted for the energy required to place the subsystem in the state (q, ). In physical terms, the mathematical steps represent first assembling the subsystem mechanically with no electrical excitation. Because there is no force acting on the electrode as it is put in place, no work is involved. Then, with its location fixed, the electrode is charged by means of an electrical source.

Suppose that the subsystem is electrically linear, so that either as a result of mathematical modeling or of measurements on the actual system, the electrical terminal relation takes the form

equation GIF #11.127
Then, with this relation used to evaluate (11), it follows that the energy is
equation GIF #11.128
Finally, the desired force of electrical origin follows from substituting this expression into (9b).
equation GIF #11.129
Note that with a similar substitution into (9a), the terminal relation of (12) is obtained. Once the partial derivative with respect to has been taken while holding the proper independent variable (q) fixed, the force can be written in terms of variables other than the independent ones. Thus, with the use of the terminal relation, (12), the force is written in terms of the terminal voltage v as

boxed equation GIF #11.29
The following example gives the opportunity to apply this result to a specific configuration.

Example 11.6.1. Force on a Capacitor Plate

The region between the plane parallel electrodes shown in Fig. 11.6.3 is filled by a layer of dielectric having permittivity and thickness b and an air gap . The total distance between electrodes, b + , is small compared to the linear dimensions of the plates, so fringing fields will be ignored. Thus, the electric fields Ea and Eb in the air gap and in the dielectric, respectively, are uniform. What force on the upper electrode results from applying the voltage v between the electrodes?

floating figure GIF #20
First we determine the charge q on the upper electrode. To this end, the integral of E from the upper electrode to the lower one must be equal to the applied voltage, so
equation GIF #11.130
Further, there is presumably no unpaired surface charge at the interface between the dielectric layer and the air gap. Thus, Gauss' continuity condition requires that
equation GIF #11.131
Elimination of Eb between these equations gives
equation GIF #11.132
In terms of this electric field at the surface of the upper electrode, Gauss' continuity condition shows that the total charge on the upper electrode is
equation GIF #11.133
and so it follows from (18) that the electrical terminal relation can be written in terms of a capacitance C.
equation GIF #11.134
Because the dielectric is described by a linear constitutive law, we have obtained an electrical terminal relation where v is a linear function of q.

The force acting on the upper electrode follows from a substitution of (20) into (15).

equation GIF #11.135
By definition, if f is positive, it acts in the direction of . Here we find that regardless of the polarity of the applied voltage, f is negative. This is to be expected, because charges of one polarity on the upper electrode are attracted toward those of opposite polarity on the lower electrode.

In describing energy conversion, a minus sign can be extremely important. For example, vdq is the incremental energy into the electroquasistatic subsystem, while f d is the energy leaving that subsystem as the force of electrical origin acts on the external mechanical system. Thus, if f is positive, it acts on the mechanical system in such a direction as to increase the associated displacement.

Rotating motors and generators are examples where the conversion of energy between electrical and mechanical form is a cyclic process. In these cases, the subsystem returns to its original state once each cycle. The energy converted per cycle is determined by integrating the energy conservation law, (6), around the closed path in the state space representing this process.

equation GIF #11.136
Because the energy stored in the system must return to its original value, there is no net contribution of the energy storage term in (22). For a cyclic process, the net electrical energy input per cycle must be equal to the net mechanical power output per cycle.

equation GIF #11.137
The following demonstration is primarily intended to give further insight into the implications of the conservation of energy principle for a cyclic process.

Demonstration 11.6.1. An Energy Conversion Cycle

The experiment shown in Fig. 11.6.4 is based on the plane parallel capacitor configuration analyzed in Example 11.6.1. The lower electrode, aluminum foil mounted on a table top, is covered by a thin sheet of plastic. The upper one is also foil, but taped to an insulating sheet which is attached to a frame. This electrode can then be manually raised and lowered to effectively control the displacement .
floating figure GIF #21
With the letters A through D used to designate states of the system, we consider the following energy conversion cycle.
\bullet If the terminals are again open (i = 0) but the
\bullet A B. With v = 0, the upper electrode rests on the plastic sheet. A voltage Vo is applied.
\bullet B C. With the voltage source removed so that the upper electrode is electrically isolated, it is raised to the position = L.
\bullet C D. The upper electrode is shorted, so that its voltage returns to zero.

\bullet D A. The upper electrode is returned to its original position at = 0.
Is electrical energy converted to mechanical form, or vice versa?
floating figure GIF #22
The process in carrying out the closed integrals on the left-hand and right-hand sides of (23) as the cycle is carried out can be pictured in the (q, v) and (f, ) planes, respectively, as shown in Fig. 11.6.5. From A B, q = C(0)v, where C( ) is given by (20). Thus, the trajectory in the (q, v) plane is a straight line ending at the voltage v = Vo of the source. Because the upper electrode has remained at its original position, the trajectory in the (f, ) plane is along the = 0 axis. The force on the electrode caused by raising its voltage to Vo follows from (21).

equation GIF #11.138
Now, from B C, the voltage source is removed so that as the upper electrode is raised to = L, its charge is conserved. This means that the trajectory in the (q, v) plane is one where q = constant = Vo C(0). The voltage reached by the upper electrode can be found by requiring that q be conserved.

equation GIF #11.139
In the experiment, the thickness b of the dielectric sheet is a fraction of a millimeter, while the final elevation = L might be 20 cm. (If the displacement is larger than this, the fringing field comes into play and the expression for the capacitance is no longer valid.) Thus, as the sheet is raised, an original voltage of 500 V is easily amplified to 10 - 20 kV. This is readily observed by means of an electrostatic voltmeter attached to the upper electrode, as shown in Fig. 11.6.4.

To determine the trajectory B C in the (f, ) plane, observe from (14) or (24) that as a function of q, the force is independent of .

equation GIF #11.140
The trajectory B C in the (f, ) plane is therefore one of constant fB given by (24). In general, f(q, ) is not independent of , but in plane parallel geometry, it is. The system is now returned to its original state in two steps. First, from C D, the upper electrode remains at = L and is shorted to ground. In the (q, v) plane, the state returns to the origin along the straight line given by q = C(L)v, (20). In the (f, ) plane, the force drops to zero with = L. Second, from D A, the upper electrode is returned to its original position. The values of (q, v) remain zero, while the trajectory in the (f, ) plane is f = 0.

The experiment is simple enough so that we can use physical reasoning to decide the direction of energy conversion. Although the force of gravity is likely to exceed the electrical force of attraction between the electrodes, as far as the electrical subsystem is concerned, the upper electrode is raised against a downward electrical force. Because the charge is removed before it is lowered, there is no electrical force on the electrode as it is lowered. Thus, net work is done on the EQS subsystem. The right-hand side of (23), the net work done by the subsystem on the external mechanical system, is thus negative. Evaluation of one or the other of the two sides of the energy conversion law, (23), provides two other ways to determine the direction of energy conversion. Consider first the electrical input energy. The integral has contributions from A B (where the source is used to charge the upper electrode) and from C D (where the electrode is discharged). The areas under the respective triangles representing the integral of vdq are

equation GIF #11.141
In view of the expression for C( ), (20), this expression can be written as
equation GIF #11.142
This expression is clearly negative, indicating that the net electrical energy flow is out of the electrical terminal pair. This is consistent with having a net mechanical energy input to the system.

The net mechanical output energy per cycle expressed by the right-hand side of (23) should be equal to (28). To see that this is so, we recognize that the integral consists of two possible contributions, from B C and D A. During the latter, f = 0, so the magnitude of the integral is simply the area of the rectangle in Fig. 11.6.5b.

equation GIF #11.143
As required by the conservation law, this mechanical energy output is negative and is equal to the net electrical input energy given by (28). With the sequence of electrical and mechanical terminal constraints described above, there is a net conversion of energy from mechanical to electrical form. The system acts as an electrical generator with energy provided by whoever raises and lowers the upper electrode. With the voltage applied when the electrode is at its largest displacement, and the electrode grounded before it is raised, the energy flow is from the electrical voltage source to the mechanical system. In this case, the system acts as a motor. We will have the opportunity to exemplify a practical motor in the next section. Most motors are MQS rather than EQS. However, a practical EQS device that is designed to convert energy from mechanical to electrical form is the capacitor microphone, a version of which was described in Example 6.3.3.

Electrical forces have their origins in forces on unpaired charges and on dipoles. The force on the upper capacitor plate of Example 11.6.1 is due to the unpaired charges. The equal and opposite force on the combination of electrode and dielectric is in part due to unpaired charges and in part to dipoles induced in the dielectric. The following exemplifies a total force that is entirely due to polarization.

Example 11.6.2. Force on a Dielectric Material

We return to the configuration of Example 11.4.1, where a dielectric slab extends a distance into the region between plane parallel electrodes, as shown in Fig. 11.6.6. The capacitance was found in Example 11.4.1 to be (11.4.16).

equation GIF #11.144
It follows from (15) that there is a force tending to draw the dielectric into the region between the electrodes.

equation GIF #11.145
This force results because dipoles induced by the fringing field experience forces in the direction that are passed on to the material in which they are embedded. Even though the force is due to the fringing fields, the net force does not depend on the details of that field. This is evident from our energy arguments because, at least as long as the upper edge of the slab is well within the region between the plates and the lower edge never reaches the vicinity of the electrodes, the energy storage in the fringing fields does not change when the slab is moved.

floating figure GIF #23
Further discussion of the force density responsible for the force on the dielectric will be given in Sec. 11.9. Its physical reality is demonstrated next.

Demonstration 11.6.2. Force on a Liquid Dielectric

In the experiment shown in Fig. 11.6.7, capacitor plates are dipped into a dish full of dielectric liquid. Thus, with the application of the voltage, it rises against gravity. To demonstrate the relationship between the voltage and the force, the spacing between the electrodes is a slowly varying function of radial position. With r denoting the radial distance from an axis where an extension of the electrodes would join, and the angle between the electrodes, the spacing at a distance r is r. Thus, in (31), the spacing between electrodes a r and the force per unit radial distance tending to push the liquid upward is a function of the radial position.

floating figure GIF #24
equation GIF #11.146
This force must raise a column of liquid having a height and width r. With the mass density of the liquid defined as , the total mass per radial distance raised by this force is therefore r . Force equilibrium is therefore represented by setting the force per unit radial length equal to this mass multiplied by the gravitational acceleration g.

equation GIF #11.147
This expression can be solved for (r).

equation GIF #11.148
In the experiment, \footnote5See film Electric Fields and Moving Media from series by National Committee for Electrical Engineering Films, Education Development Center, 39 Chapel St., Newton, Mass. 02160. the electrodes are constructed from tin-oxide coated glass. They are then both conducting and transparent. As a result, the height to which the liquid rises can be seen to obey (34), both as to its magnitude and its radial dependence on r. To obtain an appreciable rise of the liquid without exceeding the field strength for electrical breakdown between the electrodes, the atmosphere over the liquid must be pressurized. Also, to avoid effects of unpaired charges injected at high field strengths by the electrodes, the applied voltage is alternating and, because the force is proportional to the square of the applied field, the height of rise is proportional to the rms value of the voltage.

11.7 Macroscopic Magnetic Forces

In this section, an energy principle is applied to the determination of net forces in MQS systems. With Sec. 11.6 as background, it is appropriate to include the case of multiple terminal pairs. As in Sec. 11.4, the coenergy is again found to be a convenient alternative to the energy.

floating figure GIF #25
The MQS system is shown schematically in Fig. 11.7.1. It has two electrical terminal pairs and one mechanical degree of freedom. The magnetoquasistatic subsystem now described by an energy principle excludes electrical dissipation and all aspects of the mechanical system, mechanical energy storage and dissipation. The energy principle then states that the input of electrical power through the electrical terminal pairs either goes into a rate of change of the stored magnetic energy or into a rate of change of the work done on the external mechanical world.
equation GIF #11.149
As in Sec. 11.6, our starting point in finding the force is a postulated principle of energy conservation. Because the system is presumably MQS, in accordance with (11.3.29), the left-hand side represents the net flux of power into the system. With the addition of the last term and the inherent assumption that there is no electrical dissipation in the subsystem being described, (1) is more than the recasting of Poynting's theorem.

In an MQS system, the voltages are the time rates of change of the flux linkages. With these derivatives substituted into (1) and the expression multiplied by dt, it becomes

boxed equation GIF #11.30
This energy principle states that the increments of electrical energy put into the MQS subsystem (as increments of flux d1 and d2 through the terminals multiplied by their currents i1 and i2, respectively) either go into the total energy, which is increased by the amount dwm, or into work on the external mechanical system, subject to the force f and experiencing a displacement d.

With the energy principle written as in (2), the flux linkages are the independent variables. We saw in Example 11.4.2 that it is inconvenient to specify the flux linkages as functions of the currents. With the objective of casting the currents as the independent variables, we now recognize that

equation GIF #11.150
and substitute into (2) to obtain
equation GIF #11.151
where the coenergy function, seen before in Sec. 11.4, is defined as
equation GIF #11.152
We picture the MQS subsystem as having flux linkages 1 and 2, a force f and a total energy wm that are specified once the currents i1 and i2 and the displacement are stipulated. According to (4), the coenergy is a function of the independent variables i1, i2, and , wm' = wm' (i1, i2, ), and the change in wm' can also be written as
equation GIF #11.153
Because the currents and displacement are independent variables, (4) and (6) can hold only if the coefficients of like terms on the right are equal. Thus,
boxed equation GIF #11.31
The last of these three expressions is the key to finding the force f.

Reciprocity Condition

Before we find the coenergy and hence f, consider the implication of the first two expressions in (7) for the electrical terminal relations. Taking the derivative of 1 with respect to i2, and of 2 with respect to i1, shows that
boxed equation GIF #11.32
Although this reciprocity condition must reflect conservation of energy for any lossless system, magnetically linear or not, consider its implications for a system described by the linear terminal relations.

equation GIF #11.154
Application of (8) shows that energy conservation requires the equality of the mutual inductances.
equation GIF #11.155
This relation has been derived in Example 11.4.2 from a related but different point of view.

Finding the Coenergy

To find wm', we integrate (4) along a path in the state space (i1, i2, ) arranged so that the integral can be carried out without having to know f.

boxed equation GIF #11.33
Thus, the first leg of the line integral is carried out on with the currents equal to zero. Provided that f = 0 in the absence of these currents, this means that the integral of fd makes no contribution. The payoff from our formulation in terms of the coenergy rather than the energy comes in being able to carry out the remaining integration using terminal relations in which the flux linkages are expressed in terms of the currents. For the linear terminal relations of (9), this line integration was illustrated in Example 11.4.2, where it was found that
equation GIF #11.156

Evaluation of the Force

In general, the inductances in this expression are functions of . Thus, the force f follows from substituting this expression into (7c).

boxed equation GIF #11.34
Of course, this expression applies to systems having a single electrical terminal pair as a special case where i2 = 0. This generalization of the energy method to multiple electrical terminal pair systems suggests how systems with two or more mechanical degrees of freedom are treated.

Example 11.7.1. Driver for a Matrix Printer

A transducer that is similar to those used to drive an impact printer is shown in Fig. 11.7.2. The device, which is symmetric about the axis, might be one of seven used to drive wires in a high-speed matrix printer. The objective is to transduce a current i that drives the N-turn coil into a longitudinal displacement of the permeable disk at the top. This disk is attached to one end of a wire, the other end of which is used to impact the ribbon against paper, imprinting a dot.

floating figure GIF #26
The objective here is to determine the force f acting on the plunger at the top. For simplicity, we make a highly idealized model in which the magnetizable material surrounding the coil and filling its core, as well as that of the movable disk, is regarded as perfectly permeable. Moreover, the air gap spacing is small compared to the radial dimension a, so the magnetic field intensity is approximated as uniform in the air gap. The wire is so fine that the magnetizable material removed to provide clearance for the wire can be disregarded.

With Ha and Hc defined as shown by the inset to Fig. 11.7.2, Ampère's integral law is applied to a contour passing upward through the center of the core, across the air gap at the center, radially outward in the disk, and then downward across the air gap and through the outer part of the stator to encircle the winding in the infinitely permeable material.

equation GIF #11.157
A second relation between Ha and Hc follows from requiring that the net flux out of the disk must be zero.

equation GIF #11.158
Using this last expression to replace Hc in (14) results in
equation GIF #11.159
The magnetic flux linking every turn in the coil is o Hc c2. Thus, the total flux linked by the coil is

equation GIF #11.160
Finally, the force follows from an evaluation of (13) (specialized to the single terminal pair system of this example).

equation GIF #11.161
As might have been expected by one who has observed magnetizable materials pulled into a magnetic field, the force is negative. Given the definition of in Fig. 11.7.2, the application of a current will tend to close the air gap. To write a dot, the current is applied. To provide for a return of the plunger to its original position when the current is removed, a spring is inserted in the air gap.

In the magnetic transducer of the previous example, the force on the driver disk is due to magnetization. The next example illustrates the force associated with the current density.

Example 11.7.2. Force on a Wire over a Perfectly Conducting Plane

The cross-section of a perfectly conducting wire with its center a distance above a perfectly conducting ground plane is shown in Fig. 11.7.3. The configuration is familiar from Demonstration 8.6.1. The current carried by the wire is returned in the ground plane. The distribution of this current on the surfaces of the wire and ground plane is consistent with the requirement that there be no flux density normal to the perfectly conducting surfaces. What is the force per unit length f acting on the wire?

floating figure GIF #27
The inductance per unit length is half of that for a pair of conductors having the center-to-center spacing 2. Thus, it is half of that given by (8.6.12).

equation GIF #11.162
The force per unit length in the direction then follows from an evaluation of (13) (again adapted to the single terminal pair situation).

equation GIF #11.163
The dependence of this force on the elevation above the ground plane is shown in Fig. 11.7.4. In the limit where the elevation is large compared to the radius of the conductor, (20) becomes
equation GIF #11.164
floating figure GIF #28
In Sec. 11.8, we will identify the force density acting on materials carrying a current density J as being J x o H. Note that the upward force predicted by (20) is indeed consistent with the direction of this force density.

The force on a thin wire, (21), can be derived from this force density by recognizing that the contribution of the self-field of the wire to the total force per unit length is zero. Thus, the force per unit length can be computed using for B the flux density caused by the image current a distance 2 away. The flux density due to this image current has a magnitude that follows from Ampère's integral law as o i/2 (2 ). This field is essentially uniform over the cross-section of the wire, so the integral of the force density J x B over the cross-section of the wire amounts to an integration of the current density J over the cross-section. The latter is the total current i, and so we are led to a force per unit length of magnitude o i2 /2 (2 ), which is in agreement with (21).

The following is a demonstration of the force on current-carrying conductors exemplified previously. It also provides a dramatic demonstration of the existence of induced currents.

Demonstration 11.7.1. Steady State Magnetic Levitation

In the experiment shown in Fig. 11.7.5, the current-carrying wire of the previous example has been wound into a pancake shaped coil that is driven by about 20 amps of 60 Hz current. The conductor beneath is an aluminum sheet of 1.3 cm thickness. Even at 60 Hz, this conductor tends to act as a perfect conductor. This follows from an evaluation of the product of the angular frequency and the time constant m estimated in Sec. 10.2. From (10.2.17), m = o a \approx 20, where the average radius is a = 9 cm, is the sheet thickness and the sheet conductivity is given by Table 7.1.1.

floating figure GIF #29
The time average force, of the type described in Example 11.7.2, is sufficiently large to levitate the coil. As the current is increased, its height above the aluminum sheet increases, as would be expected from the dependence of the force on the height for a single wire, Fig. 11.7.4.

The Torque of Electrical Origin

In some of the most important transducers, the mechanical response takes the form of a rotation rather than a translation. The shaft shown in Fig. 11.7.6a might be attached to the rotor of a motor or generator. A force f acting through a lever arm of length r that rotates the shaft through an incremental angle d causes a displacement d = rd. Thus, the incremental work done on the mechanical system fd becomes
equation GIF #11.165
where is defined as the torque.

floating figure GIF #30
If the two terminal pair MQS system of Fig. 11.7.1 had a rotational rather than a displacement degree of freedom, the representation would be the same as has been outlined, except that f and \rightarrow . The mechanical terminal pair is now represented as in Fig. 11.7.6b. The torque follows from (13) as
equation GIF #11.166
Among the types of magnetic rotating motors and generators that could be used to exemplify the torque of (23), we now choose a synchronous machine. Although other types of motors are more common, it is a near certainty that if these words are being read with the aid of electrical illumination, the electricity used is being generated by means of a synchronous generator.

Example 11.7.3. A Synchronous Machine

The cross-section of a stator and rotor modeling a rotating machine is shown in Fig. 11.7.7. The rotor consists of a highly permeable circular cylindrical material mounted on a shaft so that it can undergo a rotation measured by the angle . Surrounding this rotor is a stator, composed of a highly permeable material. In slots, on the inner surface of the stator and on the outer surface of the rotor, respectively, are windings with sinusoidally varying turn densities. These windings, driven by the currents i1 and i2, respectively, give rise to current distributions that might be modeled by surface current densities
floating figure GIF #31
equation GIF #11.167
equation GIF #11.168
where Ns and Nr are constants descriptive of the windings. Thus, the current distribution shown on the stator in Fig. 11.7.7 is fixed and gives rise to a magnetic field having the fixed vertical axis shown in the figure. The rotor coil gives rise to a similar field except that its axis is at the angle of the rotor. The rotor magnetic axis, also shown in Fig. 11.7.7, therefore rotates with the rotor.

Electrical Terminal Relations

With the rotor and stator materials taken as infinitely permeable, the air gap fields are determined by using (24) and (25) to write boundary conditions on the tangential H and then solving Laplace's equation for the air gap magnetic fields (Secs. 9.6 and 9.7). The flux linked by the respective coils is then of the form
equation GIF #11.169
where the self-inductances Ls and Lr and peak mutual inductance M are constants. The dependence of the inductance matrix on the angle of the rotor, , can be reasoned physically. The rotor is modeled as a smooth circular cylinder, so in the absence of a rotor current i2, there can be no effect of the rotor angle on the flux linked by the stator winding. Hence, the stator self-inductance is independent of . Similar reasoning shows that the rotor self-inductance must be independent of rotor angle . The dependence of the mutual inductance is plausible because the flux 1 linked by the stator, due to the current in the rotor, must peak when the magnetic axes of the coils are aligned ( = 0) and must be zero when they are perpendicular ( = 90 degrees).

Torque Evaluation

The magnetic torque on the rotor follows directly from using (26) to evaluate (23).
equation GIF #11.170
This torque depends on the currents and in such a way that the magnetic axis of the rotor tends to align with that of the stator. With = 0, the axes are aligned and there is no torque. If is slightly positive and the currents are both positive, the torque is negative. This is as would be expected with the magnetic axis of the stator vertical and that of the rotor in the first quadrant (as in Fig. 11.7.7).

Synchronous Operation

In the synchronous mode of operation, the stator current is constrained to be sinusoidal while that on the rotor is a constant. To avoid having to describe the mechanical system, we will assume that the shaft is attached to a mechanical load that makes the angular velocity \Omega constant. Thus, the electrical and mechanical terminals are constrained so that
equation GIF #11.171
equation GIF #11.172
equation GIF #11.173
Under what circumstances can we derive a time average torque on the shaft, and hence a net conversion of energy with each rotation? With the constraints of (28), the torque follows from (27) as
equation GIF #11.174
A trigonometric identity

1\over 2 (sin (x + y) + sin (x - y)) makes the implications of this result more apparent.

equation GIF #11.175
There is no time average value of either of these sinusoidal functions of time unless one or the other of the frequencies, ( + \Omega ) and ( - \Omega ), is zero. For example, with the rotation frequency equal to that of the excitation,
equation GIF #11.176
the time average torque is
equation GIF #11.177
The dependence of the time average torque on the phase angle is shown in Fig. 11.7.8a. With between 0 and 180 degrees, there is a positive time average torque acting on the external mechanical system in the direction of rotation \Omega. In this range, the machine acts as a motor to convert energy from electrical to mechanical form. In the range of from 180 degrees to 360 degrees, energy is converted from mechanical to electrical form and operation is as a generator.

floating figure GIF #32

Stator Field Analyzed into Traveling Waves

From (30), it is clear that a time-average torque results from either a forward (\Omega = ) or a backward (\Omega = - ) rotation. This suggests that the field produced by the stator winding is the superposition of fields having magnetic axes rotating in the clockwise and counterclockwise directions. Formally, this can be seen by rewriting the stator surface current density, (24), using the electrical and mechanical constraints of (28). With the use once again of the double-angle trigonometric identity, the distribution of surface current density is separated into two parts.

equation GIF #11.178
The sinusoidal excitation produces a standing-wave surface current with nodes at = 0 and 180 degrees. This is the first distribution in Fig. 11.7.9. Analyzed as it is on the right in (33), and pictured in Fig. 11.7.9, it is the sum of two countertraveling waves. The magnetic axis of the wave traveling to the right is at = t.

floating figure GIF #33
We now have the following picture of the synchronous operation found to give rise to the time average torque. The field of the stator is composed of rotating parts, one with a magnetic axis that rotates in a clockwise direction at angular velocity , and the other rotating in the opposite direction. The frequency condition of (31) therefore represents a synchronous condition in which the "forward" component of the stator field and the magnetic axis of the rotor rotate at the same angular velocity. In view of the definition of given in (28), if is positive, the rotor magnetic axis lags the stator axis by the angle , as shown in Fig. 11.7.8. When the machine operates as a motor, the forward component of the stator magnetic field "pulls" the rotor along. When the device operates as a generator, is negative and the rotor magnetic axis leads that of the forward component of the stator field. For generator operation, the rotor magnetic axis "pulls" the forward component of the stator field.

11.8 Forces on Microscopic Electric and Magnetic Dipoles

The energy principle was used in the preceding sections to derive the macroscopic forces on polarizable and magnetizable materials. The same principle can also be applied to derive the force distributions, the force densities. For this purpose, one needs more than a purely electromagnetic description of the system. In order to develop the simple model for the force density distribution, we need the expression for the force on an electric dipole for polarizable media, and on a magnetic dipole for magnetizable media. The force on an electric dipole will be derived simply from the Lorentz force law. We have not stated a corresponding force law for magnetic charges. Even though these are not found in nature as isolated charges but only as dipoles, it is nevertheless convenient to state such a law. This will be done by showing how the electric force law follows from the energy principle. By analogy a corresponding law on magnetic charges will be derived from which the force on a magnetic dipole will follow.

floating figure GIF #34

Force on an Electric Dipole

The force on a stationary electric charge is given by the Lorentz law with v = 0.

equation GIF #11.179
A dipole is the limit of two charges of equal magnitude and opposite sign spaced a distance d apart, in the limit
equation GIF #11.180
with p being finite. Charges q of opposite polarity, separated by the vector distance d, are shown in Fig. 11.8.1. The total force on the dipole is the sum of the forces on the individual charges.
equation GIF #11.181
Unless the electric field at the location r + d of the positive charge differs from that at the location r of the negative charge, the separate contributions cancel.

In order to develop an expression for the force on the dipole in the limit where the spacing d of the charges is small compared to distances over which the field varies appreciably, (2) is written in Cartesian coordinates and the field at the positive charge expanded about the position of the negative charge. Thus, the x component is

equation GIF #11.182
The first and last terms cancel. In more compact notation, this expression is therefore
equation GIF #11.183
where we have identified the dipole moment p \equiv qd. The other force components follow in a similar fashion, with y and z playing the role of x. The three components are then summarized in the vector expression
boxed equation GIF #11.35
The derivation provides an explanation of how p E is evaluated in Cartesian coordinates. The i-th component of (5) is obtained by dotting p with the gradient of the i-th component of E.

Illustration. Force on a Dipole

Suppose that a dipole finds itself in the field
equation GIF #11.184
which is familiar from Example 4.1.1. It follows from (5) that the force is
equation GIF #11.185
According to this expression, the y-directed dipole on the x axis in Fig. 11.8.2 experiences a force in the x direction. The y-directed force is zero because Ey is the same at the respective locations of the charges. The x-directed force exists because Ex goes from being positive just above the x axis to negative just below. Thus, the x-directed contributions to the force of each of the charges is in the same direction.

floating figure GIF #35

Force on Electric Charge Derived from Energy Principle

The force on an electric charge is stated in the Lorentz law. This law is also an ingredient in Poynting's theorem, and in the identification of energy and power flow. Indeed, E Ju was recognized from the Lorentz law as the power density imparted to the current density of unpaired charge. The energy principle can be used to derive the force law on a microscopic charge "in reverse". This seems to be the hard way to obtain the Lorentz law of force on a stationary charge. Yet we go through the derivation for three reasons.

shown to be consistent with the Lorentz force on a stationary charge.

(a) The derivation of force from the EQS energy principle is
(b) The derivation shows that the field can be produced by permanently polarized material objects, and yet the energy principle can be employed in a straightforward manner.

(c) The same principle can be applied to derive the microscopic MQS force on a magnetic charge.
Let us consider an EQS field produced by charge distributions and permanent polarizations Pp in free space as sketched in Fig. 11.8.3a. By analogy, we will then have found the force on a magnetic charge in the field of a permanent magnet, Fig. 11.8.3b. The Poynting theorem identifies the rate of energy imparted to the polarization per unit volume as

equation GIF #11.186
floating figure GIF #36
Because Pp is a permanent polarization, Pp / t = 0, and the permanent polarization does not contribute to the change in energy associated with introducing a point charge. Hence, as charge is brought into the vicinity of the permanent polarization, the change of energy density is
equation GIF #11.187
where stands for the differential change of o E. The change of energy is
equation GIF #11.188
where V includes all of space. The electroquasistatic E field is the negative gradient of the potential
equation GIF #11.189
Introducing this into (10), one has
equation GIF #11.190
where we have "integrated by parts," using an identity.

\footnote7 ( ) A = \nabla \cdot (A ) - (\nabla A) The first integral can be written as an integral over the surface enclosing the volume V. Since V is all of space, the surface is at infinity. Because E vanishes at infinity at least as fast as 1/r3 (1/r2 for E, 1/r for , where r is the distance from the origin of a coordinate system mounted within the electroquasistatic structure), the surface integral vanishes. Now

equation GIF #11.191
from Gauss' law, where u is the change of unpaired charge. Thus, from (12),
equation GIF #11.192
so the change of energy is equal to the charge increment u dv introduced at r times the potential at r, summed over all the charges. Suppose that one introduces only a small test charge q, so that u dv = q at point r. Then
equation GIF #11.193
The change of energy is the potential at the point at which the charge is introduced times the charge. This form of the energy interprets the potential of an EQS field as the work to be done in bringing a charge from infinity to the point of interest. If the charge is introduced at r + r, then the change in total energy associated with introducing that charge is
equation GIF #11.194
Introduction of a charge q at r, subsequent removal of the charge, and introduction of the charge at r + r is equivalent to the displacement of the charge from r to r + r. If there is a net energy decrease, then work must have been done by the force f exerted by the field on the charge. The work done by the field on the charge is

equation GIF #11.195
and therefore
equation GIF #11.196
Thus, the Lorentz law for a stationary charge is implied by the EQS laws. Before we attack the problem of force on a magnetic charge, we explore some features of the electroquasistatic case. In (17), q is a small test charge. Electric test charges are available as electrons. But suppose that in analogy with the magnetic case, no free electric charge was available. Then one could still produce a test charge by the following artifice. One could polarize a very long-thin rod of cross-section a, with a uniform polarization density P along the axis of the rod (Sec. 6.1). At one end of the rod, there would be a polarization charge q = Pa, at the other end there would be a charge of equal magnitude and opposite sign. If the rod were of very long length, while the end with positive charge could be used as the "test charge," the end of opposite charge would be outside the field and experience no force. Here the charge representing the polarization of the rod has been treated as unpaired. We are now ready to derive the force on a magnetic charge.

Force on a Magnetic Charge and Magnetic Dipole

The attraction of a magnetizable particle to a magnet is the result of the force exerted by a magnetic field on a magnetic dipole. Even in this case, because the particle is macroscopic, the force is actually the sum of forces acting on the microscopic atomic constituents of the material. As pointed out in Secs. 9.0 and 9.4, the magnetization characteristics of macroscopic media such as the iron particle relate back to the magnetic moment of molecules, atoms, and even individual electrons. Given that a particle has a magnetic moment m as defined in Sec. 8.2, what is the force on the particle in a magnetic field intensity H? The particle can be comprised of a macroscopic material such as a piece of iron. However, to distinguish between forces on macroscopic media and microscopic particles, we should consider here that the force is on an elementary particle, such as an atom or electron. We have shown how one derives the force on an electric charge in an electric field from the energy principle. The electric field could have been produced by permanently polarized dielectric bodies. In analogy, one could produce a magnetic field by permanently magnetized magnetic bodies. In the EQS case, the test charge could have been produced by a long, uniformly and permanently polarized cylindrical rod. In the magnetic case, an "isolated" magnetic charge could be produced by a long, uniformly and permanently magnetized rod of cross-sectional area a. If the magnetization density is M, then the analogy is
equation GIF #11.197
where, for the uniformly magnetized rod, and the magnetic charge
equation GIF #11.198
is located at one end of the rod, the charge -qm at the other end of the rod (Example 9.3.1). The force on a magnetic charge is thus, in analogy with (18),
boxed equation GIF #11.36
which is the extension of the Lorentz force law for a stationary electric charge to the magnetic case. Of course, the force on a dipole is, in analogy to (5) (see Fig. 11.8.4),
boxed equation GIF #11.37
where m is the magnetic dipole moment.
floating figure GIF #37
We have seen in Example 8.3.2 that a magnetic dipole of moment m can be made up of a circulating current loop with magnitude m = ia, where i is the current and a the area of the loop. Thus, the force on a current loop could also be evaluated from the Lorentz law for electric currents as
equation GIF #11.199
with i the total current in the loop. Use of vector identities indeed yields (22) in the MQS case. Thus, this could be an alternate way of deriving the force on a magnetic dipole. We prefer to derive the law independently via a Lorentz force law for stationary magnetic charges, because an important dispute on the validity of the magnetic dipole model rested on the correct interpretation of the force law[1-3]. While the details of the dispute are beyond the scope of this textbook, some of the issues raised are fundamental and may be of interest to the reader who wants to explore how macroscopic formulations of the electrodynamics of moving media based on magnetization represented by magnetic charge (Chap. 9) or by circulating currents are reconciled.

The analogy between the polarization and magnetization was emphasized by Prof. L. J. Chu[2], who taught the introductory electrical engineering course in electromagnetism at MIT in the fifties. He derived the force law for moving magnetic charges, of which (21) is the special case for a stationary charge. His approach was soon criticized by Tellegen[3], who pointed out that the accepted model of magnetization is that of current loops being the cause of magnetization. While this in itself would not render the magnetic charge model invalid, Tellegen pointed out that the force computed from (23) in a dynamic field does not lead to (22), but to

equation GIF #11.200
Because the force is different depending upon whether one uses the magnetic charge model or the circulating current model for the magnetic dipole, so his reasoning went, and because the circulating current model is the physically correct one, the magnetic charge model is incorrect. The issue was finally settled[4] when it was shown that the force (24) as computed by Tellegen was incorrect. Equation (23) assumes that i could be described as constant around the current loop and pulled out from under the integral. However, in a time-varying electric field, the charges induced in the loop cause a current whose contribution precisely cancels the second term in (24). Thus, both models lead to the same force on a magnetic dipole and it is legitimate to use either model. The magnetic model has the advantage that a stationary dipole contains no "moving parts," while the current model does contain moving charges. Hence the circulating current formalism is by necessity more complicated and more likely to lead to error.

\sectnonumComparison of Coulomb's Force on an Electron to the Force on its Magnetic Dipole


Why is it possible to accurately describe the motions of an electron in vacuum by the Lorentz force law without including the magnetic force associated with its dipole moment? The answer is that the magnetic dipole effect on the electron is relatively small. To obtain an estimate of the magnitude of the magnetic dipole effect, we compare the forces produced by a typical (but large) electric field achievable without electrical breakdown in air on the charge e of the electron, and by a typical (but large) magnetic field gradient acting on the magnetic dipole moment of the electron. Taking for E the value 106 V/m, with e \approx 1.6 x 10-19 coulomb,
equation GIF #11.201
A B of 1 tesla (10,000 gauss) is a typical large flux density produced by an iron core electromagnet. Let us assume that a flux density variation of this order can be produced over a distance of 1 cm, which is, in practice, a rather high gradient. Yet taking this value and a moment of one Bohr magneton (9.0.1), we obtain from (22) for the force on the electron
equation GIF #11.202
Note that the electric force associated with the net charge is much greater than the magnetic one due to the magnetic dipole moment. Because of the large ratio fe/fm for fields of realistic magnitudes, experiments designed to detect magnetic dipole effects on fundamental particles did not utilize particles having a net charge, but rather used neutral atoms (most notably, the Stern-Gerlach experiment

\footnote8W. Gerlach and O. Stern, "Uber die Richtungsquantelung im Magnetfeld," Ann. d. Physik, 4th series, Vol. 74, (1924), pp. 673-699.). Indeed, a stray electric field on the order of 10-6 V/m would deflect an electron as strongly as a magnetic field gradient of the very large magnitude assumed in calculating (26).

The small magnetic dipole moment of the electron can become very important in solid matter because macroscopic solids are largely neutral. Hence, the forces exerted upon the positive and negative charges within matter by an applied electric field more or less cancel. In such a case, the forces on the electronic magnetic dipoles in an applied magnetic field can dominate and give rise to the significant macroscopic force observed when an iron filing is picked up by a magnet.

Example 11.8.1. Magnetization Force on a Macroscopic Particle

Suppose that we wanted to know the force exerted on an iron particle by a magnet. Could the microscopic force, (22), be used? The energy method derivation shows that, provided the particle is surrounded by free space, the answer is yes. The particle is taken as being spherical, with radius R, as shown in Fig. 11.8.5. It is assumed to have such a large magnetizability that its permeability can be taken as infinite. Further, the radius R is much smaller than other dimensions of interest, especially those characterizing variations in the applied field in the neighborhood of the particle.

floating figure GIF #38
Because the particle is small compared to dimensions over which the field varies significantly, we can compute its moment by approximating the local field as uniform. Thus, the magnetic potential is determined by solving Laplace's equation in the region around the particle subject to the conditions that H be the uniform field Ho at "infinity" and \Psi be constant on the surface of the particle. The calculation is fully analogous to that for the electric potential surrounding a perfectly conducting sphere in a uniform electric field. In the electric analog, the dipole moment was found to be (6.6.5), p = 4 o R3 E. Therefore, it follows from the analogy provided by (19) that the magnetic dipole moment at the particle location is
equation GIF #11.203
Directly below the magnet, H has only a z component. Thus, the dipole moment follows from (27) as
equation GIF #11.204
Evaluation of (22) therefore gives
equation GIF #11.205
where Hz and its derivative are evaluated at the location of the particle.

A typical axial distribution of Hz is shown in Fig. 11.8.6 together with two pictures aimed at gaining insight into the origins of the magnetic dipole force. In the first, the dipole is again depicted as a pair of magnetic monopoles, induced to form a moment collinear with the H. Because the field is more intense at the north pole of the particle than at the south pole, there is then a net force.

floating figure GIF #39
Alternatively, suppose that the dipole is actually a circulating current, so that the force is given by (23). Even though the energy argument makes it clear that the force is again given by (22), the physical picture is different. Because H is solenoidal, an intensity that increases with z implies that the field just off axis has a component that is directed radially inward. It is this radial component of the flux density crossed with the current density that results in an upward force on each segment of the loop.

R E F E R E N C E S

  1. P. Penfield, Jr., and H. A. Haus, Electrodynamics of Moving Media, MIT Press, Cambridge, Mass. (1967).

  2. R. M. Fano, L. J. Chu, and R. B. Adler, Electromagnetic Fields, Energy, and Forces, John Wiley and Sons, New York (1960).

  3. D. B. H. Tellegen, "Magnetic-dipole models," Am. J. Phys. Vol. 30 (Sept. 1962), pp. 650-652.

  4. H. A. Haus and P. Penfield, Jr., "Force on a current loop," Phys. Lett., Vol. 26A (March 1968), pp. 412-413.

11.9 Macroscopic Force Densities

A macroscopic force density F(r ) is the force per unit volume acting on a medium in the neighborhood of r. Fundamentally, the electromagnetic force density is the result of forces acting on those microscopic particles embedded in the material that are charged, or that have electric or magnetic dipole moments. The forces acting on these individual particles are passed along through interparticle forces to the macroscopic material as a whole. In the limit where that volume becomes small, the force density can then be regarded as the sum of the microscopic forces over a volume element V.

equation GIF #11.206
Of course, the linear dimensions of V are large compared to the microscopic scale.

Strictly, the forces in this sum should be evaluated using the microscopic fields. However, we can gain insight concerning the form taken by the force density by using the macroscopic fields in this evaluation. This is the basis for the following discussions of the force densities associated with unpaired charges and with conduction currents (the Lorentz force density) and with the polarization and magnetization of media (the Kelvin force density). To be certain that the usage of macroscopic fields in describing the force densities is consistent with that implicit in the constitutive laws already introduced to describe conduction, polarization, and magnetization, the electromagnetic force densities should be derived using energy arguments. These derivations are extensions of those of Secs. 11.6 and 11.7 for forces. We end this section with a discussion of the results of such derivations and of circumstances under which they will predict the same total forces or even material deformations as those derived here.

The Lorentz Force Density

Without restricting the generality of the resulting force density, suppose that the electrical force on a material is due to two species of charged particles. One has N+ particles per unit volume, each with a charge q+, while the other has density N- and a charge equal to -q-. With v denoting the velocity of the macroscopic material and v\pm representing the respective velocities of the carriers relative to that material, the Lorentz force law gives the force on the individual particles.

equation GIF #11.207
equation GIF #11.208
Note that q- is a positive number.

In typical solids and fluids, the charged particles are either bonded to the material or migrate relative to the material, suffering many collisions with the neutral material during times of interest. In either case, the inertia of the particles is inconsequential, so that on the average, the forces on the individual particles are passed along to the macroscopic material. In either situation, the force density on the material is the sum of (2) and (3), respectively, multiplied by the charged particle densities.

equation GIF #11.209
Substitution of (2) and (3) into this expression gives the Lorentz force density
boxed equation GIF #11.38
where u is the unpaired charge density (7.1.6) and J is the current density.

equation GIF #11.210
Because the material is in motion, with velocity v, the current density J has not only the contribution familiar from Sec. 7.1 (7.1.4) due to the migration of the carriers relative to the material, but one due to the net charge carried by the moving material as well.

In EQS systems, the first term in (5) usually outweighs the second, while in MQS systems (where the unpaired charge density is negligible), the second term tends to dominate.

The derivation and Fig. 11.9.1 suggest why the electric term is proportional to the net charge density. In a given region, the force density resulting from the positively charged particles tends to be canceled by that due to the negatively charged particles, and the net force density is therefore proportional to the difference in absolute magnitudes of the charge densities. We exploited this fact in Chap. 7 to let electrically induced material motions evidence the distribution of the unpaired charge density. For example, in Demonstration 7.5.1, the unpaired charge density was restricted to an interface, and as a result, the motion of the fluid was suppressed by constraining the interface. A more recent example is the force on the upper electrode in the capacitor transducer of Example 11.6.1. Here again the force density is confined to a thin region on the surface of the conducting electrode.

floating figure GIF #40
floating figure GIF #41
The magnetic term in (6), pictured in Fig. 11.9.2 as acting on a current-carrying wire, is also familiar. This force density was responsible for throwing the metal disk into the air in the experiment described in Sec. 10.2. The force responsible for the levitation of the pancake coil in Demonstration 11.7.1 was also the net effect of the Lorentz force density, acting either over the volume of the coil conductors or over that of the conducting sheet below. In MQS systems, where the contribution of the "convection" current uv is negligible, the current density is typically due to conduction. Note that this means that the velocity of the charge carriers is determined by the electric field they experience in the conductor, and not simply by the motion of the conductor. The current density J in a moving conductor is generally not in the direction of motion.

\footnote9Indeed, it is fortunate that the carriers do not have the same velocity as the material, for if they did, it would not be possible to use the magnetic Lorentz force density for electromechanical energy conversion. If we recognize that the rate at which a force f does work on a particle that moves at the velocity v is v f, then it follows from the Lorentz force law, (1.1.1), that the rate of doing work on individual particles through the agent of the magnetic field is v (v x o H). The cross-product is perpendicular to v, so this rate of doing work must be zero.


The Kelvin Polarization Force Density

If microscopic particles carrying a net charge were the only contributors to a macroscopic force density, it would not be possible to explain the forces on polarized materials that are free of unpaired charge. Example 11.6.2 and Demonstration 11.6.2 highlighted the polarization force. The experiment was carried out in such a way that the dielectric material did not support unpaired charge, so the force is not explained by the Lorentz force density.

In EQS cases where u = 0, the macroscopic force density is the result of forces on the microscopic particles with dipole moments. The resulting force density is fundamentally different from that due to unpaired charges; the forces p E on the individual microscopic particles are passed along by interparticle forces to the medium as a whole. A comparison of Fig. 11.9.3 with Fig. 11.9.1 emphasizes this point. For a single species of particle, the force density is the force on a single dipole multiplied by the number of dipoles per unit volume Np. By definition, the polarization density P = Np p, so it follows that the force density due to polarization is

boxed equation GIF #11.39
This is often called the Kelvin polarization force density.

floating figure GIF #42

Example 11.9.1. Force on a Dielectric Material

In Fig. 11.9.4, the cross-section of a pair of electrodes that are dipped into a liquid dielectric is shown. The picture might be of a cross-section from the experiment of Demonstration 11.6.2. With the application of a potential difference to the electrodes, the dielectric rises between the electrodes. According to (7), what is the distribution of force density causing this rise?
floating figure GIF #43
For the liquid dielectric, the polarization constitutive law is taken as linear [(6.4.2) and (6.4.4)]
equation GIF #11.211
so that with the understanding that is a function of position (uniform in the liquid, o in the gas, and taking a step at the interface), the force density of (7) becomes
equation GIF #11.212
By using a vector identity

A A
= ( x A) x A + fraction GIF #26 \nabla (A A) and invoking the EQS approximation where x E = 0, this expression is written as
equation GIF #11.213
A second identity

11 ( ) = + \nabla
converts this expression into one that will now prove useful in picturing the distribution of force density.

equation GIF #11.214
Provided that the interface is well removed from the fringing fields at the top and bottom edges of the electrodes, the electric field is uniform not only in the dielectric and gas above and below the interface between the electrodes, but through the interface as well. Thus, throughout the region between the electrodes, there is no gradient of E, and hence, according to (7), no Kelvin force density. The Kelvin force density is therefore confined to the fringing field region where the fluid surrounds the lower edges of the electrodes. In this region, is uniform, so the force density reduces to the first term in (11). Expressed by this term, the direction and magnitude of the force density is determined by the gradient of the scalar E E. Thus, where E is varying in the fringing field, it is directed generally upward and into the region of greater field intensity, as suggested by Fig. 11.9.4. The force on the dipole shown by the inset lends further credence to the dipolar origins of the force density.

Although there is no physical basis for doing so, it might seem reasonable to take the force density caused by polarization as being p E. After all, it is the polarization charge density p that was used in Chap. 6 to represent the effect of the media on the macroscopic electric field intensity E. The experiment of Demonstration 11.6.2, pictured in Fig. 11.6.7, makes it clear that this force density is not correct. With the interface well removed from the fringing fields, there is no polarization charge density anywhere in the liquid, either at the interface or in the fringing field. If pE were the correct force density, it would be zero throughout the fluid volume except at the interfaces with the conducting electrodes. There, the forces are perpendicular to the surface of the electrodes. Such a force distribution could not cause the fluid to rise.

The Kelvin Magnetization Force Density

Forces caused by magnetization are probably the most commonly experienced electromagnetic forces. They account for the attraction between a magnet and a piece of iron. In Example 11.7.1, this force density acts on the disk of magnetizable material.

Given that the magnetizable material is made up of microscopic dipoles, each experiencing a force of the nature of (11.8.22), and that the magnetization density M is the number of these per unit volume multiplied by m, it follows from the arguments of the preceding section that the force density due to magnetization is

boxed equation GIF #11.40
This is sometimes called the Kelvin magnetization force density.

Example 11.9.2. Force Density in a Magnetized Fluid

With the dielectric liquid replaced by a ferrofluid having a uniform permeability , and the electrodes replaced by the pole faces of an electromagnet, the physical configuration shown in Fig. 11.9.4 becomes the one of Fig. 11.9.5, illustrating the magnetization force density. In such fluids[1], the magnetization results from an essentially permanent suspension of magnetized particles. Each particle comprises a magnetic dipole and passes its force on to the liquid medium in which it is suspended. Provided that the magnetization obeys a linear law, the discussion of the distribution of force density given in Example 11.9.1 applies equally well here.

floating figure GIF #44

Alternative Force Densities

We now return to comments made at the beginning of this section. The fields used to express the Lorentz and Kelvin force densities are macroscopic. To assure consistency between the averages implied by these force densities and those already inherent in the constitutive laws, an energy principle can be used. The approach is a continuum version of that exemplified for lumped parameter systems in Secs. 11.7 and 11.8. In the lumped parameter systems, electrical terminal relations were used to determine a total energy, and energy conservation was used to determine the force. In the continuum system[2], the electrical constitutive law is used to find an energy density, and energy conservation used, in turn, to find a force density. This energy method, like the one exemplified in Secs. 11.7 and 11.8 for lumped parameter systems, describes systems that are loss free. In making practical use of the result, it is assumed that it will be applicable even if there are losses. A more general method, which invokes a principle of virtual power[3], allows for dissipation but requires more empirical information than the polarization or magnetization constitutive law as a starting point.

Force densities derived from more rigorous arguments than given here can have very different distributions from the superposition of the Lorentz and Kelvin force densities. We would expect that the arguments break down when the microscopic particles become so densely packed that the field experienced by one is significantly altered by its nearest neighbor. But surely the difference between the magnetic force density of Lorentz and Kelvin (LK)

equation GIF #11.215
we have derived here and the Korteweg-Helmholtz force density (KH) for incompressible media
equation GIF #11.216
cited in the literature[2] is not due to interactions between microscopic particles. This latter force density is often obtained for an incompressible material from energy arguments. [Note that with - and H H, respectively, playing the roles of dL/d and i2, the magnetization term in (14) takes a form found for the force on a magnetizable material in Sec. 11.7.] In Example 11.9.2 (where J = 0), we found the force density of (13) to be confined to the fringing field. By contrast, (14) gives no force density in the fringing region (where is uniform), but rather puts it all at the interface. According to this latter equation, through the agent of a surface force density (a force density that is a spatial impulse at the interface), the field pulls upward on the interface.

The question may then be asked whether, and how, the two force density expressions can be reconciled. The answer is that if o M = ( - o)H, they predict the same motion for any volume-conserving material deformations such as those of an incompressible fluid. We shall demonstrate this for the case of a liquid, such as shown in Fig. 11.9.5, but allowing for the action of a current J as well. As the first step in the derivation, we shall show that (13) and (14) differ by the gradient of a scalar, (r). To see this, use a vector identity


A A
= (\nabla x A) x A + fraction GIF #27 (A A) to write (13) as
equation GIF #11.217
The MQS form of Ampère's law makes it possible to substitute x H for J in this expression, which then becomes
equation GIF #11.218
The second term in this expression is then expanded using a second vector identity

13 ( ) = \nabla +
equation GIF #11.219
This expression differs from (14) by the last term, which indeed takes the form where
equation GIF #11.220
Now consider Newton's force law for an elemental volume of material. Using the Korteweg-Helmholtz force density, (14), it takes the form

equation GIF #11.221
where p is the internal fluid pressure and Fm is the sum of all other mechanical contributions to the force density. Alternatively, using (13) written as (17) as the force density, this same law is represented by

equation GIF #11.222
For an incompressible material, none of the other laws needed to describe the continuum (such as mass conservation) involve the pressure.


14 For example, for a compressible fluid, the pressure depends on mass density and temperature, so the pressure does appear in the physical laws. Indeed, in the constitutive law relating these values, the pressure has a well-defined value. However, in an incompressible fluid, the constitutive law relating the pressure to mass density and temperature is not relevant to the prediction of material motion. Thus, if (19) is used, p appears only in that equation and if (20) is used, p' \equiv p - appears only in that expression. This means that p and p' play identical roles in predicting the deformation. In an incompressible material, it is the role of the pressure to adjust itself so that only volume conserving deformations are allowed.

\footnote15Like the "perfectly permeable material" of magnetic circuits, in which B remains finite as H goes to zero, the "perfectly incompressible" material is one in which the pressure remains finite even as the material becomes infinitely "stiff" to all but those deformations that conserve volume. The two formulations would differ in what one would call the pressure, but would result in the same material deformation and velocity. An example is the height of rise of the fluid between the parallel plates in Fig. 11.9.5.

Included in the class of incompressible deformations are rigid body motions. If used self-consistently, force densities that differ by the gradient of a "" will predict the same motions of rigid bodies. Thus, the net force on a body surrounded by free space will be the same whether found using the Lorentz-Kelvin or the Korteweg-Helmholtz force density. The following example illustrates this concept.

Example 11.9.3. Magnetic Force on a Magnetizable Current-Carrying Material

A block of conducting material having permeability is shown in Fig. 11.9.6 sandwiched between perfectly conducting plates. A current source, distributed over the left edges of these electrodes, drives a constant surface current density K in the +x direction along the left edge of the lower electrode. This current passes through the block in the y direction as a current density
floating figure GIF #45
equation GIF #11.223
and is returned to the source in the -x direction at the left edge of the upper electrode. The thickness a of the block is small compared to its other two dimensions, so the magnetic field between the electrodes is z directed and dependent only on x. From Ampère's law it follows that
equation GIF #11.224
in the conducting block.

The alternative force densities, (13) and (14), have very different distributions in the block. Yet we must find that the net force on the block, found by integrating each over its volume, is the same. To see that this is so, consider first the sum of the Lorentz and Kelvin force densities, (13).

There is no x component of the magnetic field intensity, so for this particular configuration, the magnetization term makes no contribution to (13). Evaluation of the first term using (19) and (20) then gives

equation GIF #11.225
Integration of this force density over the volume amounts to a multiplication by the cross-sectional area ad, and integration on x. Thus, the net force predicted by using the force density of Lorentz and Kelvin is
equation GIF #11.226
Now, the Korteweg-Helmholtz force density given by (14) is evaluated. The permeability is uniform throughout the interior of the block, so the magnetization term is again zero there. However, is a step function at the ends of the block, where x = -b and x = 0. Thus, is an impulse there and we must take care to include the contributions from the surface regions in our integration. Evaluation of the x component of (14) using (21) and

equation GIF #11.227
Integration of (25) over the volume of the block therefore gives
equation GIF #11.228
Note that Hz is constant through the interface at x = -b. Thus, the integration of the last term can be carried out. Simplification of this expression gives the same total force as found before, (24).

The distributions of the force densities given by (13) and (14) are generally different, even very different. It is therefore natural to ask which of the two is the "right" one. In general, until the "other" force densities acting on the medium in question are specified, this question cannot be answered. Here, where a discussion of continuum mechanics is beyond our purview, we have identified a class of mechanical deformations (namely, those that are volume conserving or "incompressible"), where these force densities are equally valid. In fact, any other force density differing from these by a term having the form would also be valid. The combined Lorentz and Kelvin force densities have the advantage of a satisfying physical interpretation. However, the derivation has the weakness of making an ad hoc use of the macroscopic fields. Force densities resulting from an energy argument have the advantage of dealing rigorously with the macroscopic fields.

R E F E R E N C E S

  1. R. E. Rosensweig, "Magnetic Fluids," Scientific American, (Oct. 1982), pp. 136-145.

  2. J. R. Melcher, Continuum Electromechanics, MIT Press, Cambridge, Mass. (1981), chap. 3.
  3. P. Penfield and H. A. Haus, Electrodynamics of Moving Media, MIT Press, Cambridge, Mass. (1967).

11.10 Summary

Far reaching as they are, the laws summarized by Maxwell's equations are directly applicable to the description of only one of many physical subsystems of scientific and engineering interest. Like those before it, this chapter has been concerned with the electromagnetic subsystem. However, by casting the electromagnetic laws into statements of power flow, we have come to recognize how the electromagnetic subsystem couples to the thermodynamic subsystem through the power dissipation density and to the mechanical subsystem through forces and force densities of electromagnetic origin.

The basis for a self-consistent macroscopic description of any continuum subsystem is a power flow statement having the forms identified in Sec. 11.1. Describing the energy and power flow in and into a volume V enclosed by a surface S, the integral conservation of energy statement takes the form (11.1.1).

equation GIF #11.229
The differential form of the conservation of energy statement is implied by the above.

equation GIF #11.230
Poynting's theorem, the subject of Sec. 11.2, is obtained starting from the laws of Faraday and Ampère to obtain an expression of the form of (2). For materials that are Ohmic (J = E) and that are linearly polarizable and magnetizable (D = E and B = H), the power flux density S (or Poynting's vector), energy density W, and power dissipation density Pd were shown in Sec. 11.3 to be
equation GIF #11.231
equation GIF #11.232
equation GIF #11.233
Of course, taking the free space limit where and assume their free space values and = 0 gives the free space conservation statement discussed in Sec. 11.2.

In Sec. 11.3, we found that in EQS systems, an alternative to Poynting's vector is (11.3.24).

equation GIF #11.234
This expression is of practical importance, because it can be evaluated without determining H, which is generally not of interest in EQS systems. An important application of the integral form of the energy conservation statement is to lumped parameter systems. In these cases, the surface S of (1) encloses a system that is connected to the outside world through terminals. It is then convenient to describe the power flow in terms of the terminal variables. It was shown in Sec. 11.3 (11.3.29), that the net power into the system represented by the left-hand side of (1) becomes
equation GIF #11.235
provided that the magnetic induction and the electric displacement current through the surface S are negligible. This set the stage for the application of the integral form of the energy conservation theorem to lumped parameter systems.

In Sec. 11.4, attention focused on the energy storage term, the first terms on the right in (1) and (2). The energy density concept was broadened to include materials having constitutive laws relating the flux densities to the field intensities that were single valued and collinear. With E, D, H, and B representing the field magnitudes, the energy density was found to be the sum of electric and magnetic energy densities.

equation GIF #11.236
Integrated over the volume V of a system, this function leads to the total energy w. For quasistatic lumped parameter systems, the total electric or magnetic energy is often conveniently found following a different route. First the terminal relations are determined and then the total energy is found by adding up the increments of energy put into the system as it is energized. In the case of an n terminal pair EQS system, where the relation between terminal voltage vi and associated charge qi is vi (q1, q2, \ldots qn ), the increment of energy is vi dqi, and the total electric energy is (11.4.9).

equation GIF #11.237
The line integration in an n-dimensional space representing the n independent qi's was illustrated by Example 11.4.2. Similarly, for an n terminal pair MQS system where the current ii is related to the flux linkage i by ii = ii (1, 2, \ldots n), the total energy is (11.4.12).

equation GIF #11.238
Note the analogy between these expressions for the total energy of EQS and MQS lumped parameter systems and the electric and magnetic energy densities, respectively, of (8). The transition from the field picture afforded by the energy densities to the lumped parameter characterization is made by E v, D q and by H \rightarrow i, B \rightarrow .

Especially in using the energy to evaluate forces of electrical origin, we found it convenient to define coenergy density functions.

equation GIF #11.239
It followed that these functions were natural when it was desirable to use E and H as the independent variables rather than D and B.

equation GIF #11.240
The total coenergy functions for lumped parameter EQS and MQS systems could be found either by integrating these densities over the volume or by again viewing the system in terms of its terminal variables. With the total coenergy functions defined by
equation GIF #11.241
it followed that the coenergy functions could be determined from the terminal relations by again carrying out line integrations, but this time with the voltages and currents as the independent variables. For EQS systems,

equation GIF #11.242
while for MQS systems,
equation GIF #11.243
Again, note the analogy to the respective terms in (12).

The remaining sections of the chapter developed some of the possible implications of the "dissipation" term in the energy conservation statement, the last terms in (1) and (2). In Sec. 11.5, coupling to a thermal subsystem was discussed. In this section, the disparity between the power input and the rate of increase of the energy stored was accounted for by heating. In addition to Ohmic heating, caused by collisions between the migrating carriers and the neutral media, we considered losses associated with the dynamic polarization and magnetization of materials.

In Secs. 11.6-11.9, we considered coupling to a mechanical subsystem as a second mechanism by which energy could be extracted from (or put into) the electromagnetic subsystem. With the displacement of an object denoted by , we used an energy conservation postulate to infer the total electric or magnetic force acting on the object from the energy functions [(11.6.9), and its magnetic analog]

equation GIF #11.244
or from the coenergy functions [(11.7.7) and the analogous expression for electric systems].

equation GIF #11.245
In Sec. 11.8, where the Lorentz force on a particle was generalized to account for electric and magnetic dipole moments, one objective was a microscopic picture that would lend physical insight into the forces on polarized and magnetized materials. The Lorentz force was generalized to include the force on stationary electric and magnetic dipoles, respectively.

equation GIF #11.246
The total macroscopic forces resulting from microscopic forces had already been encountered in the previous two sections. The force density describes the interaction between a volume element of the electromagnetic subsystem and a mechanical continuum. The force density inferred by averaging over the forces identified in Sec. 11.8 as acting on microscopic particles was
equation GIF #11.247
A more rigorous approach to finding the force density could be based on a generalization of the energy method introduced in Secs. 11.6 and 11.7. As background for further pursuit of this subject, we have illustrated the importance of including the mechanical continuum with which the force density acts. Before there can be a meaningful answer to the question, "Which force density is correct?" the other force densities acting on the material must be specified. As an illustration, we found that very different electric or magnetic force densities would result in the same deformations of an incompressible material and in the same net force on an object surrounded by free space[1,2].

R E F E R E N C E S

  1. P. Penfield, Jr., and H. A. Haus, Electrodynamics of Moving Media, MIT Press, Cambridge, Mass. (1967).

  2. J. R. Melcher, Continuum Electromechanics, MIT Press, Cambridge, Mass. (1981), chap. 3.

A capacitor C, an inductor L, and a resistor R are in series, driven by the voltage v(t) and carrying the current i(t). With vc defined as the voltage across the capacitor, show that vi = dw/dt + i2 R where w = fraction GIF #28 C vc2 + 1\over 2 Li2. Argue that w is the energy stored in the inductor and capacitor, while i2 R is the power dissipated in the resistor.

P R O B L E M S

 

Introduction
\problem1*

 

Integral and Differential Conservation Statements

11.1.1*Consider a system in which the fields are y and/or z directed and independent of y and z. Then S = Sx (x, t)ix, W = W(x, t), and Pd = Pd(x, t).

and located between x = x1 and x = x2, (1) becomes

equation GIF #11.248
(a) Show that for a volume having area A in any y - z plane
(b) Take the limit where x1 - x2 = x 0 and show that the one-dimensional form of (3) results.
(c) Based on (a), argue that Sx is the power flux density in the x direction.
 

Poynting's Theorem

11.2.1*The perfectly conducting plane parallel electrodes of Fig. 13.1.1 are driven at the left by a voltage source Vd (t) and are "open circuit" at the right, as shown in Fig. 13.1.4. The system is EQS.

y/a2)Vd dVd/dt.

(a) Show that the power flux density is S = iy (-o
(b) Using S, show that the power input is d(fraction GIF #29 CVd2)/dt, where C = o b w/a.
(c) Evaluate the right-hand side of (11.1.1) to show that if the magnetic energy storage is neglected, the same result is obtained.
(d) Show that the magnetic energy storage is indeed negligible if b/c is much shorter than times of interest.

11.2.2The perfectly conducting plane parallel electrodes of Fig. 13.1.1 are driven at the left by a current source Id(t), as shown in Fig. 13.1.3. The system is MQS.

(a) Determine S.
(b) From S, find the input power.
(c) Evaluate the right-hand side of (11.1.1) for a volume enclosing the region between the electrodes, and show that if the electric energy storage is neglected, it is indeed equal to the left-hand side.
(d) Under what conditions is the electric energy storage negligible?

 

Ohmic Conductors with Linear Polarization and Magnetization

11.3.1*In Example 7.3.2, a three-dimensional dipole current source drives circulating currents through a uniformly conducting material. This source is so slowly varying with time that time rates of change have a negligible effect. Consider first the power flow as pictured in terms of the Poynting flux density, (3).

equation GIF #11.249
(a) Show that
(b) Show that
equation GIF #11.250
(c) Using these results, show that (11.1.3) is indeed satisfied.

(d) Now, using the alternative EQS power theorem, evaluate S as given by (23) and again show that (11.1.3) is satisfied.
(e) Observe that the latter evaluation is much simpler to carry out and that the latter power flux density is easier to picture.
11.3.2Coaxial perfectly conducting circular cylindrical electrodes make contact with a uniformly conducting material of conductivity in the annulus b < r < a, as shown in Fig. P11.3.2. The length l is large compared to a. A voltage source v drives the system at the left, while the electrodes are "open" at the right. Assume that v(t) is so slowly varying that the voltage can be regarded as independent of z.

floating figure GIF #46
in the annulus.
(a) Determine E, , and H
(b) Evaluate the Poynting power flux density S [as given by (3)] in the annulus.
(c) Use S to evaluate the total power dissipation by integration over the surface enclosing the annulus.
(d) Show that the same result is obtained by integrating Pd over the volume.
(e) Evaluate S as given by (23), and use that distribution of the power flux density to determine the total power dissipation.
(f) Make sketches of the alternative distributions of S.
(g) Show that the input power is vi, where i is the total current from the voltage source.

11.3.3*A pair of perfectly conducting circular plates having a spacing d form parallel electrodes in a system having cylindrical symmetry about the z axis and the cross-section shown by Fig. P11.3.3. The central region between the plates is filled out to the radius b by a uniformly conducting material having conductivity and uniform permittivity , while the surrounding region, where b < r < a, is free space. A distributed voltage source v(t) constrains the potential difference between the outer edges of the electrodes. Assume that the system is EQS.
equation GIF #11.251
(a) Show that the Poynting power flux density is
(b) Integrate this flux density over a surface enclosing the region between the plates, and show that it is equal to the sum of the rate of change of electric energy storage and the power dissipation.
(c) Now show that the alternative power flux density given by (23) is
equation GIF #11.252
(d) Carry out part (b) using this distribution of S, and show that the result is the same.
(e) Show that the power input is equal to vi, where i is the total current from the voltage source.

floating figure GIF #47
11.3.4In Example 7.5.1, the steady current distribution in and around a conducting circular cylindrical rod immersed in a conducting material was determined. Assume that Eo is so slowly varying that it can be regarded as static. flux density S, as given by (3).
(a) Determine the distribution of Poynting power
(b) Determine the alternative S given by (23).
(c) Find the power dissipation density Pd in and around the rod.
(d) Show that the differential energy conservation law [(11.1.3) with W/ t = 0] is satisfied at each point in and around the rod using either of these distributions of S.

 

Energy Storage

11.4.1*In Example 8.5.1, the inductance L of a spherically shaped coil was found by "adding up" the flux linkages of the individual windings. Taking an alternative approach to finding L, use the fields found in that example to determine the total energy storage, wm. Then use the fact that wm = fraction GIF #30 Li2 to show that L is as given by (8.5.20).

11.4.2In Prob. 9.6.3, a coil has turns at the interface between a magnetizable material and a circular cylindrical core of free space, as shown in Fig. P9.6.3. Assume that the system has a length l in the z direction and determine the total energy, wm. (Assume that the rotatable coil carries no current.) Use the fact that wm = fraction GIF #31 L i2 to find L.

11.4.3*In Example 8.6.4, the fields of a coil distributed throughout a volume were found. Using these fields to evaluate the total energy storage, show that the inductance is as given by (8.6.35).
11.4.4The magnetic circuit described in Prob. 9.7.5 and shown in Fig. P9.7.5 has two electrical excitations. Determine the total magnetic coenergy, w'm (i1, i2, x).

11.4.5The cross-section of a motor or generator is shown in Fig. 11.7.7.

the total coenergy wm'.
(a) Determine the magnetic coenergy density Wm', and hence
(b) By writing wm' in the form of (11.4.24), determine L11, L12, and L22.

11.4.6*The material in the system of Fig. 11.4.3 has the constitutive law of (28). Show that the total coenergy is
equation GIF #11.253
11.4.7Consider the system shown in Fig. P9.5.1 but with a = o and the region where B = b H now filled with a material having the constitutive law
equation GIF #11.254
(a) Determine B and H in each region.
(b) Find the coenergy density in each region and hence the total coenergy wm' as a function of the driving current i.

 

Electromagnetic Dissipation

11.5.1*In Example 7.9.2, the Maxwell capacitor has an area A (perpendicular to x), and the terminals are driven by a source v = Re [\hat v \exp (j t)]. The sinusoidal steady state has been established. Show that the time average power dissipation in the lossy dielectrics is
equation GIF #11.255
11.5.2In Example 7.9.3, the potential is found in the EQS approximation in and around a lossy dielectric sphere embedded in a lossy dielectric and stressed by a uniform field having a sinusoidal dependence on time (7.9.36). region.
(a) Find the time average power dissipation density in each
(b) What is the total time average power dissipated in the sphere?
11.5.3*Plane parallel perfectly conducting plates having the spacing d are shorted by a perfectly conducting sheet in the plane x = 0, as shown in Fig. P11.5.3. A sheet having thickness and conductivity is in the plane x = -b and makes contact with the perfectly conducting plates above and below. At their left edges, in the plane x = -(a + b), a source of surface current density, K(t), is connected to the plates. The regions to left and right of the resistive sheet are free space, and w is large compared to a, b, and d.
figure GIF #1
stored as defined on the right in (11.1.1), are
equation GIF #11.256
(a) Show that the total power dissipation and magnetic energy
(b) Show that the integral on the left in (11.1.1) over the surface indicated by the dashed line in the figure gives the same result as found in part (a).

11.5.4In Example 10.4.1, the applied field is Ho (t) = Hm cos ( t) and sinsuoidal steady state conditions prevail. Determine the time average power dissipation in the conducting sheet.
figure GIF #2
11.5.5*The cross-section of an N-turn circular solenoid having radius a is shown in Fig. P11.5.5. It surrounds a thin cylindrical shell of square cross-section, with length b on a side. This shell has thickness and conductivity , and is filled by a material having permeability . Both the shell and the solenoid have a length d perpendicular to the paper that is large compared to a. current i1 = io cos t and the sinusoidal steady state has been established, integrate the time average power dissipation density over the volume of the shell to show that the total time-average power dissipation is
equation GIF #11.257
(a) Given that the terminals of the solenoid are driven by the
(b) In the sinusoidal steady state, the time average Poynting flux through a surface enclosing the shell goes into the time average dissipation. Use this fact to obtain (a).

11.5.6In describing the response of macroscopic media to fields in the sinusoidal steady state, it is convenient to use complex constitutive laws. The complex permittivity is introduced by (19). Here we introduce and illustrate the complex permeability. Suppose that field quantities take the form
equation GIF #11.258
density, the MQS laws require that
equation GIF #11.259
equation GIF #11.260
equation GIF #11.261
(a) Show that in a region where there is no macroscopic current
(c) Given that the spherical shell of Prob. 10.4.3 comprises each element in the cubic array of Fig. P11.5.6, each sphere with spacing s such that s \gg R, what is the complex permeability defined such that \hat B = \hat \hat H?
(d) A macroscopic material composed of this array of spheres is placed in the one-turn solenoid of rectangular cross-section shown in Fig. P11.5.6. This configuration is long enough in the z direction so that fringing fields can be ignored. At their left edges, the perfectly conducting plates composing the top and bottom of the solenoid are driven by a distributed current source, K(t). With the fringing fields in the neighborhood of the left end ignored, the resulting fields take the form H = Hz (x, t) iz and E = Ey (x, t)iy. Use an evaluation of the Poynting flux to determine the total time average power dissipated in the length l, width d, and height a of the material.
floating figure GIF #48
11.5.7*In the limit where the skin depth is small compared to the length b, the magnetic field distribution in the conductor of Fig. 10.7.2 is given by (10.7.15). Show that (per unit y - z area) the time average power dissipation associated with the current flowing in the "skin" region is |Ks|2/2 watts/m2.

11.5.8The conducting block shown in Fig. 10.7.2 has a length d in the z direction.
(a) Determine the total time average power dissipation.
(b) Show that in the case \ll b this expression reduces to that obtained in Prob. 11.5.7, while in the limit \gg b, the result is i2R where R is the dc resistance of the slab and i is the total current.

11.5.9*The toroid of Fig. 9.4.1 is filled with an insulating material having the magnetization constitutive law of Prob. 9.4.3. Show that from the terminals of the N1-turn coil, the circuit is equivalent to one having an inductance L = o N12 w2 /8R in series with a resistance Rm = o N12 w2/8R.

11.5.10The toroid of Fig. 9.4.1 is filled by a material having the magnetization characteristic shown in Fig. P11.5.10. A sinusoidal current is supplied with a particular amplitude, i = (2Hc 2 R/N1) cos ( t).

figure GIF #3
(a) Draw a dimensioned plot of B(t).
(b) Find the terminal voltage v(t) and also make a dimensioned plot.
(c) Compute the time average power input, defined as
equation GIF #11.262
where T = 2 /.
(d) Show that the result of part (c) can also be found by recognizing that, during one cycle, there is an energy/unit volume dissipated which is equal to the area enclosed by the B-H characteristic.

 

Electrical Forces on Macroscopic Media

11.6.1*A pair of perfectly conducting plates, the upper one fixed and the lower one free to move with the horizontal displacement , have a fixed spacing a as shown in Fig. P11.6.1. Show that the force of electrical origin acting on the lower electrode in the direction is f = - o v2 d/2a.

figure GIF #4
11.6.2In Example 4.6.3, the capacitance per unit length of the pair of parallel circular cylindrical conductors shown in Fig. 4.6.6 was found. Determine the force per unit length acting on the right cylinder in the x direction.

11.6.3*The electric transducer shown in cross-section by Fig. P11.6.3 has cylindrical symmetry about the center line. A coaxial pair of perfectly conducting electrodes having length l are excited at the left end by a voltage source v(t). A perfectly insulating dielectric material having permittivity is free to slide in and out of the annular region between electrodes.
figure GIF #5
dielectric material in the axial direction is f = v2 ( - o )/ln (a/b).
(a) Show that the force of electric origin acting on the
(b) Show that if the electrical terminals are constrained by the circuit shown, R is very small and the plunger suffers the displacement (t) the output voltage is vo = -2 RV( - o )(d /dt)/ln (a/b).

11.6.4The electrometer movement shown in Fig. P11.6.4 consists of concentric, perfectly conducting tubes, the inner one free to move in the axial direction. force of electrical origin acting in the direction of .
(a) Ignore the fringing field and determine the
(b) For the energy conversion cycle of Demonstration 11.6.1, but for this transducer, make dimensioned plots of the cycle in the (q, v) and (f, ) planes (analogous to those of Fig. 11.6.5).
(c) By calculating both, show that the electrical energy input in one cycle is equal to the work done on the external mechanical system.
floating figure GIF #49
11.6.5*Show that the vertical force on the nonlinear dielectric material of Prob. 11.4.6 is
equation GIF #11.263
 

Macroscopic Magnetic Fields

11.7.1*Show that the force acting in the x direction on the movable element of Prob. 9.7.5 (Note Prob. 11.4.4.) is
equation GIF #11.264
11.7.2Determine the force f(i, ) acting in the x direction on the plunger of the magnetic circuit shown in Fig. P9.7.6.

floating figure GIF #50
11.7.3*The magnetic transducer shown in Fig. P11.7.3 consists of a magnetic circuit in which the lower element is free to move in the x and y directions. From the energy principle, ignoring fringing fields, show that the force on this element is
equation GIF #11.265
11.7.4The magnetic circuit shown in cross-section by Fig. P11.7.4 has cylindrical symmetry. A plunger of permeability having outer and inner radii a and b can suffer a displacement into the annular gap of a magnetic circuit otherwise made of infinitely permeable material. The coil has N turns. Assume that the left end of the plunger is well within the magnetic circuit, so that fringing fields can be ignored, and determine the force f(i, ) acting to displace the plunger in the direction.

floating figure GIF #51
11.7.5*The "variable reluctance" motor shown in cross-section in Fig. P11.7.5 consists of an infinitely permeable yoke and an infinitely permeable rotor element forming a magnetic circuit with two air gaps of length \ll R. The system has depth d \gg \Delta into the paper. Assume that 0 < < , as shown, and show that the torque caused by passing a current i through the two N-turn coils is = - o Rd N2 i2/.

floating figure GIF #52
11.7.6A "two-phase" synchronous machine is constructed having a cross-section like that shown in Fig. 11.7.7, except that there is an additional winding on the stator. This is identical to the one shown except that it is rotated 90 degrees in the clockwise direction. The current in the stator winding shown in Fig. 11.7.7 is denoted by ia, while that in the additional winding is ib. Thus, the magnetic axes of ia and ib, respectively, are upward and to the right. With Ls, Lr, and M given constants, the inductance matrix is
equation GIF #11.266
(a) Determine the coenergy wm' (ia, ib, ).
(b) Find the torque on the rotor, (ia, ib, ).

(c) With ia = I cos ( t) and ib = I sin ( t), where I and are given constants, argue that the magnetic axis produced by the stator rotates with the angular velocity .

(d) Using these current constraints together with ir = Ir and = \Omega t - , where Ir, and \Omega are constants, show that under synchronous conditions (where = \Omega), the torque is = MI Ir sin ( ).

 

Forces on Microscopic Electric and Magnetic Dipoles

11.8.1*In a uniform electric field E, a perfectly conducting particle having radius R has a dipole moment p = 4 o R3 E. Provided that R is short compared to distances over which the field varies, this gives a good approximation to p, even where the field is not uniform. Such a particle is shown at the location x = X, y = Y in Fig. P11.8.1, where it is subject to the field produced by a periodic potential = Vo cos ( x) imposed in the plane y = 0.

Vo cos ( x) exp (- y).
(a) Show that the potential imposed in the region 0 < y is
(b) Show that, provided that the particle has no net charge, the force on the particle is
equation GIF #11.267

figure GIF #6
11.8.2The perfectly conducting particle described in Prob. 11.8.1, carrying no net charge but polarized by the imposed electric field, is subjected to the field of a charge Q located at the origin of a spherical coordinate system. In terms of its location R relative to the charged particle at the origin, determine the force on the particle.

figure GIF #7
11.8.3*In Fig. P11.8.3, permanent magnets in the lower half-space are represented by the magnetization density M = Mo cos ( x)iy, where Mo and are given positive constants.

the upper half-space is

equation GIF #11.268
(a) Show that the resulting magnetic potential in
(b) A small infinitely permeable particle having the radius R is located at x = X, y = Y. Show that the magnetization force on the particle is as given by (a) of Prob. 11.8.1, with Vo (Mo/2 ) and o \rightarrow o.

11.8.4A small "infinitely permeable" particle of radius R is a distance Z above an infinitely permeable plane, as shown in Fig. P11.8.4. A uniform field H = Ho iz is imposed. Assume that R \ll Z, and use (27) to approximate the dipole moment induced in the particle. The effect of the infinitely permeable plane on the field induced by this dipole is equivalent to that of a second image dipole located at z = -Z. Thus, there is a force of attraction between the magnetized particle and the infinite plane that is equivalent to that attracting the dipole to its image. Determine the force in the z direction on the particle.
floating figure GIF #53
 

Macroscopic Force Densities

11.9.1In Prob. 11.7.2, the total force on a magnetizable plunger is found (Fig. P9.7.6). Find this same force by integrating the force density, (14), over the volume of the plunger.
11.9.2*In Example 10.3.1, the transient current induced by applying a magnetic field intensity Ho to a conducting shell is determined.

acting on the shell Tr = o K(Ho + Hi)/2. (Note that the thin-shell model implies that H varies in an essentially linear fashion with R inside the shell.)
(a) Show that there is a radial magnetic force per unit area
(b) Specifically, show that
equation GIF #11.269

11.9.3In Example 10.4.1, the transient current induced in a conducting shell by the application of a transverse magnetic field is found. Suppose that the magnetizable core is absent.

shell is Tr = o K(Ho + Hi)/2. (Note that according to the thin-shell model, H has an essentially linear dependence on r within the shell.)
(a) Show that the radial force per unit area acting on the
(b) Determine Tr ( ,t) and relate the result to Demonstration 10.4.1.

*********-here Hold on to this material for now (8-1-86). It may be used in Chapter 15.

The third objective of this chapter has been to see the quasistatic approximations from the perspective of electrodynamics. This began in Sec. 12.2 with the consideration of the turn-on transient of an electric dipole. It continued in Sec. 12.5 with a study of the sinusoidal steady state standing waves on the parallel plate transmission line. In the low frequency limit

@eq[s=0.2,n=15] an "open-circuit" termination resulted in a capacitor while with the "short-circuit" termination the system took on the characteristics of an inductor. These limiting cases were found in Sec. 12.6 by respectively making the EQS and MQS approximations at the outset. That section concluded with a formal procedure for making the quasistatic approximations in systems composed of perfect conductors and perfect insulators. The EQS and MQS approximations were identified as the lowest order fields in a time-rate-parameter expansion.

In materials of finite conductivity, additional processes that depend on the time-rate-of-change contribute to the distribution of the fields. Because fields are not only influenced by material properties, length scales and time scales but by the topology as well, there is no simple procedure for identifying quasistatic systems or subsystems. However, for systems having all relevant dimensions of the same order, typically of length @s[L], it was shown in Sec. 12.8 that a rough idea of the relevant physical processes that could occur on a time scale @g[t] could be obtained by considering the subsystems position in the @s[L]-@g[t] plane. A summary of regimes for electromagnetic subsystems having one characteristic length @s[L] and composed of linear materials is shown in Fig. 12.9.1. The length scale has been normalized to the matching length @s[L]@+[@m[8]] defined by Eq. 13 while the characteristic time @g[t] has been normalized to the charge relaxation time /@g[s]. Thus, the lines demarking the regimes are

@eq[s=0.6,n=15] For a simple system having one characteristic length to be quasistatic we must have @g[t]@-[em]@s[ @s[L] is large compared to @s[L]@+[@m[8]] and EQS if @s[L] is small compared to @s[L]@+[@m[8]].

-here But surely the difference between the magnetic force density of Lorentz and Kelvin

equation GIF #11.270
we have derived here and the Korteweg-Helmholtz force density for incompressible media
equation GIF #11.271
is not due to interactions between microscopic particles. This latter force density is often obtained for an incompressible material energy arguments. (Note that, with - and H H respectively playing the roles of dL/d and i2, the magnetization term in (14) takes a form found for the force on a magnetizable material in Sec. 11.7.) In Example 11.9.2 (where J = 0), we found the force density of (13) to be confined to the fringing field. By contrast, (14) gives no force density in the fringing region (where is uniform) but rather puts it all at the interface. According to this latter equation, through the agent of a surface force density (a force density that is a spatial impulse at the interface) the field pulls upward on the interface.

Even though the force densities of (13) and (14) have very different distributions, they predict the same height of rise of the liquid! This is because the liquid deformations being considered are essentially incompressible, in the sense that, (with v denoting the liquid velocity) v = 0 \footnote*Solids can also deform in an essentially incompressible fashion. An example is the low frequency motion of jello or muscle.. In the force equation representing an incompressible material, there will always be another force density taking the form p, where the pressure p assumes whatever distribution it must to insure that the deformations are incompressible. As a result, contributions to the force density that take the form will have no effect on the incompressible deformations. The contribution of is simply balanced by that due to p.

For a magnetically linear material (o M = ( = \mo )H) the force densities of (13) and (14 do indeed only differ by a term taking the form . To see this, use a vector identity


A A = (\nabla x A) x A + fraction GIF #32 A A to write (13) as
equation GIF #11.272
The MQS form of Ampère's law makes it possible to substitute J for x H in this expression, which then becomes
equation GIF #11.273
The second term in this expression is then expanded using a second vector identity \footnote* x (\Psi A) = \Psi \nabla A + A \cdot
equation GIF #11.274
This expression differs from (14) by the last term, which indeed takes the form where
equation GIF #11.275
Thus, although they have very different distributions, the force densities of (13) and (14) will predict the same incompressible deformations of a material. Will the integration of these two force densities over the volume of an object result in two different total forces? Provided that the object is surrounded by media having the properties of free space, the answer is no. This is because the integral over the volume of the force density by which the two differ is zero. To see this, consider the x component (say) of that integral.

equation GIF #11.276
With the objective of converting this volume integral to one over the enclosing surface, we write this expression in the equivalent form
equation GIF #11.277
It follows from Gauss' Theorem that the volume integral on V is equivalent to one over the enclosing surface S.

equation GIF #11.278
With the surface S taken as enclosing an object surrounded by free space, = o on S and is zero everywhere on S. We conclude that integration of either (13) or (14) over the volume will give the same total force. In summary, if an object is surrounded by free space, integration over its volume of two force densities that differ by the gradient of a scalar that is zero in free space will result in the same net force.

Example 11.10.3. Magnetic Force on a Magnetizable Current Carrying Material

A block of conducting material having permeability is shown in Fig. 11.9.6 sandwiched between perfectly conducting plates. A current source, distributed over the left edges of these electrodes, drives a constant surface current density K in the +x direction along the left edge of the lower electrode. This current passes through the block in the y direction as a current density
floating figure GIF #54
equation GIF #11.279
and is returned to the source in the -x direction at the left edge of the upper electrode. The thickness a of the block is small compared to its other two dimensions, so the magnetic field between the electrodes is z directed and dependent only on x. From Ampère's law it follows that
equation GIF #11.280
in the conducting block.

The alternative force densities, (13) and (14), have very different distributions in the block. Yet, we must find that the net force on the block found by integrating each over its volume is the same. To see that this is so consider first the sum of the Lorentz and Kelvin force densities, (13).

There is no x component of the magnetic field intensity, so for this particular configuration, the magnetization term makes no contribution to (13). Evaluation of the first term using (22) and (23) then gives

equation GIF #11.281
Integration of this force density over the volume amounts to a multiplication by the cross-sectional area ab and integration on x.

equation GIF #11.282
Now, the force density given by (14) is evaluated. The permeability is uniform throughout the interior of the block, so the magnetization term is again zero there. However, is a step function at the ends of the block, where x = -b and x = 0. Thus, is an impulse there and we must take care to include the contributions from the surface regions in our integration. Evaluation of the x component of (14) using (22) and (23) gives
equation GIF #11.283
Note that Hz is K at x = -b and is zero at x = 0. Integration of (26) over the volume of the block therefore gives

equation GIF #11.284
Note that Hz is constant through the interface at x = -b. So, the integration of the last term can be carried out. Simplification of this expression gives the same total force as found before, (25).

The distributions of the force densities given by (13) and (14) are generally different, even very different. It is natural to therefore ask which of the two is the "right" one. In general, until the "other" force densities acting on the medium in question are specified, this question cannot be answered. Here, where a discussion of continuum mechanics is beyond our purvue, we have identified a class of mechanical continua (namely incompressible materials), where these force densities are equally valid. In fact so would any other force density differing from these by a term having the form . The combined Lorentz and Kelvin force densities have the advantage of a satisfying physical interpretation. However, the derivation has the weakness of making an "ad hoc" use of the macroscopic fields. Force densities resulting from an energy argument have the advantage of dealing rigorously with the macroscopic fields. The form illustrated by (14) also has the advantage of concentrating the magnetization force density at interfaces. For example, it is then clear that the height to which the liquid rises in the experiment of Fig. 11.9.4 depends on the field intensity at the interface and not on the details of the fringing field!

-here

11.11 Macroscopic Force Densities

A macroscopic force density f(r ) is the force per unit volume acting on a medium in the neighborhood of r. Fundamentally, the electromagnetic force density is the result of forces acting on those microscopic particles embedded in the material that are charged, or that have electric or magnetic dipole moments. The understanding is that the forces acting on these individual particles is passed along through interparticle forces to the macroscopic material as a whole. In the limit where that volume becomes small, the force density can then be regarded as the sum of the microscopic forces over a volume element V.

equation GIF #11.285
Of course, the linear dimensions of V are not so small as the microscopic scale.

Strictly, the forces in this sum should be evaluated using the microscopic fields. However, we can gain insight concerning the form taken by the force density by using the macroscopic fields in this evaluation. This is the basis for the following discussions of the force densities associated with unpaired charges and with conduction currents (the Lorentz force density) and with the polarization and magnetization of media (the Kelvin Force Density). To be certain that the usage of macroscopic fields in describing the force densities is consistent with that implicit in the constitutive laws already introduced to describe conduction, polarization and magnetization, the electromagnetic force densities should be derived using energy arguments. These derivations are extensions of those of Secs. 11.6 and 11.7 for forces. We end this section with a discussion of the results of such derivations and of circumstances under which they will predict the same total forces or even material deformations as those derived here.

The Lorentz Force Density

Without restricting the generality of the resulting force density, suppose that the electrical force on a material is due to two species of charged particles. One has N+ particles per unit volume, each with a charge q+, while the other has density N- and a charge equal to -q-. With v denoting the velocity of the macroscopic material and v\pm representing the respective velocities of the carriers relative to that material, the Lorentz force law gives the force on the individual particles.

equation GIF #11.286
equation GIF #11.287
Note that q- is a positive number.

In typical solids and fluids, the charged particles are either bonded to the material or migrate relative to the material, suffering many collisions with the neutral material during times of interest. In either case, the inertia of the particles is inconsequential so that on the average the forces on the individual particles is passed along to the macroscopic material. In either case, the force density on the material is the sum of (2) and (3) respectively multiplied by the charged particle densities.

equation GIF #11.288
Substitution of (2) and (3) into this expression gives the Lorentz force density.

equation GIF #11.289
where u is the unpaired charge density (7.1.6) and J is the current density.

equation GIF #11.290
Because the material is in motion, with velocity v, the current density J has not only the contribution familiar from Sec. 7.1 (7.1.4) due to the migration of the carriers relative to the material, but one due to the net charge carried by the moving material as well.

In EQS systems, the first term in (5) usually outweighs the second, while in MQS systems (where the unpaired charge density is negligible) the second term dominates.

The derivation and Fig. 11.9.1 suggest why the electric term is proportional to the net charge density. In a given region, the force density resulting from the positively charged particles tends to be canceled by that due to the negatively charged particles and the net force density is therefore proportional to the difference in absolute magnitudes of the charge densities. We exploited this fact in Chap. 7 to let electrically induced material motions evidence the distribution of the unpaired charge density. For example, in Demonstration 7.5.1, the unpaired charge density was restricted to an interface and as a result the motion of the fluid was suppressed by constraining the interface. A more recent example is the force on the upper electrode in the capacitor transducer of Example 11.6.1. Here again, the force density is confined to a thin region on the surface of the conducting electrode.

floating figure GIF #55
The magnetic term in (6), pictured in Fig. 11.9.2 as acting on a current carrying wire, is also familiar. It was this force density that was responsible for throwing the metal disk into the air in the experiment described in Sec. 10.2. The force responsible for the levitation of the pan-cake coil in Demonstration 11.7.1 was also the net effect of the Lorentz force density, either acting over the volume of the coil conductors, or over that of the conducting sheet below. In MQS systems, where the contribution of the "convection" current uv is negligible, the current density is typically due to conduction. Note that this means that the velocity of the charge carriers is determined by the electric field they experience in the conductor, and not simply by the motion of the conductor. The current density J in a moving conductor is generally not in the direction of motion.

\footnote*Indeed, it is fortunate that the carriers do not have the same velocity as the material, for if they did it would not be possible to use the magnetic Lorentz force density for electromechanical energy conversion. If we recognize that the rate at which a force f does work on a particle that moves at the velocity v is v f then it follows from the Lorentz force, (1.1.1), that the rate of doing work on individual particles through the agent of the magnetic field is v (v x o H). The cross-product is perpendicular to v, so this rate of doing work must be zero. How the law of induction and Ohm's law can be used to determine the current density in a moving conductor is considered in Appendix A4.

The Kelvin Polarization Force Density

If microscopic particles carrying a net charge were the only contributors to a macroscopic force density, it would not be possible to explain the forces on polarized materials that are free of unpaired charge. Example 11.6.2 and Demonstration 11.6.2 highlighted the polarization force. The experiment was carried out in such a way that the dielectric material did not support unpaired charge, so the force is not explained by the Lorentz force density.

In these EQS cases where u = 0, the macroscopic force density is the result of forces on microscopic particles that have dipole moments. The resulting force density is fundamentally different from that due to unpaired charges, the forces p E on the individual microscopic particles are passed along by interparticle forces to the medium as a whole. A comparison of Fig. 11.9.3 to Fig. 11.9.1 emphasizes this point. For a single species of particles, the force density is the force on a single dipole multiplied by the number of dipoles per unit volume Np. By definition, the polarization density P = Np p, so it follows that the force density due to polarization is

equation GIF #11.291
This is often called the Kelvin polarization force density.

floating figure GIF #56

Example 11.11.1. Force on a Dielectric Material

In Fig. 11.9.4, the cross-section of a pair of electrodes that are dipped into a liquid dielectric is shown. The picture might be of a cross-section from the experiment of Demonstration 11.6.2. With the application of a potential difference to the electrodes, the dielectric rises between the electrodes. According to (7), what is the distribution of force density causing this rise?
floating figure GIF #57
For the liquid dielectric, the polarization constitutive law is taken as linear (6.4.2)
equation GIF #11.292
so that with the understanding that is a function of position (uniform in the liquid, o in the gas and taking a step at the interface) the force density of (7) becomes
equation GIF #11.293
By using a vector identity \footnote*
equation GIF #11.294

and invoking the EQS approximation where x E = 0, this expression is written as
equation GIF #11.295
A second vector identity \footnote**
equation GIF #11.296
converts this expression into one that will now prove useful in picturing the distribution of force density.

equation GIF #11.297
Provided that the interface is well removed from the fringing fields at the top and bottom edges of the electrodes, the electric field is uniform not only in the dielectric and gas above and below the interface between the electrodes, but through the interface as well. Thus, throughout the region between the electrodes there is no gradient of E, and hence, according to (7), no Kelvin force density. The Kelvin force density is therefore confined to the fringing field region where the fluid surrounds the lower edges of the electrodes. In this region, is uniform, so the force density reduces to the first term in (11). Expressed by this term, the direction and magnitude of the force density is determined by the gradient of the scalar E E. Thus, where E is varying in the fringing field, it is directed generally upward and into the region of greater field intensity as suggested by Fig. 11.9.4. The force on the dipole shown by the inset lends further credence to the dipolar origins of the force density.

Although there is no physical basis for doing so, it might seem reasonable to take the force density caused by polarization as being p E. After all, it is the polarization charge density p that was used in Chap. 6 to represent the effect of the media on the macroscopic electric field intensity E. The experiment of Demonstration 11.6.2, pictured in Fig. 11.9.4, is classic because it makes it clear that this force density is not correct. With the interface well removed from the fringing fields, there is no polarization charge density anywhere in the liquid, either at the interface or in the fringing field! If pE were the correct force density, it would be zero throughout the fluid volume.

The Kelvin Magnetization Force Density

Probably forces caused by magnetization are the most commonly experienced of electromagnetic forces. They account for the attraction between a magnet and a piece of iron. In Example 11.7.1, it is this force density that acts on the disk of magnetizable material.

Given that the magnetizable material is made up of microscopic dipoles, each experiencing a force of the nature of (11.8.6), and that the magnetization density M is the number of these per unit volume multiplied by m, it follows that the force density due to magnetization is

equation GIF #11.298
This is sometimes called the Kelvin magnetization force density.

Example 11.11.2. Force Density in a Magnetized Fluid

With the dielectric liquid replaced by a ferrofluid having a uniform permeability and the electrodes replaced by the pole faces of an electromagnet, the physical configuration shown in Fig. 11.9.4 becomes the one of Fig. 11.9.5, illustrating the magnetization force density. In such fluids

\footnote*Rosensweig, R.E., Magnetic Fluids, Scientific American, Oct. 1982, pp136-145., the magnetization results from an essentially permanent suspension of magnetized particles. Each particle comprises a magnetic dipole and passes its force on to the liquid medium in which it is suspended. Provided that the magnetization obeys a linear law, the discussion of the distribution of force density given in Example 11.9.1 applies equally well here.

floating figure GIF #58

Alternative Force Densities

We now return to comments used to introduce this section. The fields used to express the Lorentz and Kelvin force densities are macroscopic. To assure consistency between the averages implied by these force densities and those already inherent to the constitutive laws, an energy principle can be used. The approach is a continuum version of that exemplified for lumped-parameter systems in Secs. 11.7 and 11.8. In the lumped parameter systems, electrical terminal relations were used to determine a total energy and energy conservation was used to determine the force in turn. In the continuum system,

\footnote*Melcher, J.R., Continuum Electromechanics, M.I.T. Press, 1981, Chap. 3. the electrical constitutive law is used to find an energy density and energy conservation used in turn to find a force density. This energy method, like the one exemplified in Secs. 11.7 and 11.8 for lumped parameter systems, describes systems that are loss-free. In making practical use of the result it is assumed that it will be applicable even if there are losses. A more general method, which invokes a principle of virtual power, \footnote**Penfield, P. and Haus, H.A., Electrodynamics of Moving Media, M.I.T. Press, Cambridge, Mass., 1967 allows for dissipation but requires more empirical information than the polarization or magnetization constitutive law as a starting point.

Force densities derived from more rigorous arguments then given here can have very different distributions from the superposition of the Lorentz and Kelvin force densities. We expect that the arguments used here prove inadequate when the microscopic particles become so densely packed that the field experienced by one is significantly altered by its nearest neighbor. But surely the difference between the magnetic force density of Lorentz and Kelvin

equation GIF #11.299
we have derived here and the Korteweg-Helmholtz force density for incompressible media
equation GIF #11.300
is not due to interactions between microscopic particles. This latter force density is often obtained for an incompressible material energy arguments. (Note that, with - and H H respectively playing the roles of dL/d and i2, the magnetization term in (14) takes a form found for the force on a magnetizable material in Sec. 11.7.) In Example 11.9.2 (where J = 0), we found the force density of (13) to be confined to the fringing field. By contrast, (14) gives no force density in the fringing region (where is uniform) but rather puts it all at the interface. According to this latter equation, through the agent of a surface force density (a force density that is a spatial impulse at the interface) the field pulls upward on the interface.

Even though the force densities of (13) and (14) have very different distributions, they predict the same height of rise of the liquid! This is because the liquid deformations being considered are essentially incompressible, in the sense that, (with v denoting the liquid velocity) v = 0 \footnote*Solids can also deform in an essentially incompressible fashion. An example is the low frequency motion of jello or muscle.. In the force equation representing an incompressible material, there will always be another force density taking the form p, where the pressure p assumes whatever distribution it must to insure that the deformations are incompressible. As a result, contributions to the force density that take the form will have no effect on the incompressible deformations. The contribution of is simply balanced by that due to p.

For a magnetically linear material (o M = ( = \mo )H) the force densities of (13) and (14 do indeed only differ by a term taking the form . To see this, use a vector identity


A A = (\nabla x A) x A + fraction GIF #33 A A to write (13) as
equation GIF #11.301
The MQS form of Ampère's law makes it possible to substitute J for x H in this expression, which then becomes
equation GIF #11.302
The second term in this expression is then expanded using a second vector identity \footnote* x (\Psi A) = \Psi \nabla A + A \cdot
equation GIF #11.303
This expression differs from (14) by the last term, which indeed takes the form where
equation GIF #11.304
Thus, although they have very different distributions, the force densities of (13) and (14) will predict the same incompressible deformations of a material. Will the integration of these two force densities over the volume of an object result in two different total forces? Provided that the object is surrounded by media having the properties of free space, the answer is no. This is because the integral over the volume of the force density by which the two differ is zero. To see this, consider the x component (say) of that integral.

equation GIF #11.305
With the objective of converting this volume integral to one over the enclosing surface, we write this expression in the equivalent form
equation GIF #11.306
It follows from Gauss' Theorem that the volume integral on V is equivalent to one over the enclosing surface S.

equation GIF #11.307
With the surface S taken as enclosing an object surrounded by free space, = o on S and is zero everywhere on S. We conclude that integration of either (13) or (14) over the volume will give the same total force. In summary, if an object is surrounded by free space, integration over its volume of two force densities that differ by the gradient of a scalar that is zero in free space will result in the same net force.

Example 11.11.3. Magnetic Force on a Magnetizable Current Carrying Material

A block of conducting material having permeability is shown in Fig. 11.9.6 sandwiched between perfectly conducting plates. A current source, distributed over the left edges of these electrodes, drives a constant surface current density K in the +x direction along the left edge of the lower electrode. This current passes through the block in the y direction as a current density
floating figure GIF #59
equation GIF #11.308
and is returned to the source in the -x direction at the left edge of the upper electrode. The thickness a of the block is small compared to its other two dimensions, so the magnetic field between the electrodes is z directed and dependent only on x. From Ampère's law it follows that
equation GIF #11.309
in the conducting block.

The alternative force densities, (13) and (14), have very different distributions in the block. Yet, we must find that the net force on the block found by integrating each over its volume is the same. To see that this is so consider first the sum of the Lorentz and Kelvin force densities, (13).

There is no x component of the magnetic field intensity, so for this particular configuration, the magnetization term makes no contribution to (13). Evaluation of the first term using (22) and (23) then gives

equation GIF #11.310
Integration of this force density over the volume amounts to a multiplication by the cross-sectional area ab and integration on x.

equation GIF #11.311
Now, the force density given by (14) is evaluated. The permeability is uniform throughout the interior of the block, so the magnetization term is again zero there. However, is a step function at the ends of the block, where x = -b and x = 0. Thus, is an impulse there and we must take care to include the contributions from the surface regions in our integration. Evaluation of the x component of (14) using (22) and (23) gives
equation GIF #11.312
Note that Hz is K at x = -b and is zero at x = 0. Integration of (26) over the volume of the block therefore gives

equation GIF #11.313
Note that Hz is constant through the interface at x = -b. So, the integration of the last term can be carried out. Simplification of this expression gives the same total force as found before, (25).

The distributions of the force densities given by (13) and (14) are generally different, even very different. It is natural to therefore ask which of the two is the "right" one. In general, until the "other" force densities acting on the medium in question are specified, this question cannot be answered. Here, where a discussion of continuum mechanics is beyond our purvue, we have identified a class of mechanical continua (namely incompressible materials), where these force densities are equally valid. In fact so would any other force density differing from these by a term having the form . The combined Lorentz and Kelvin force densities have the advantage of a satisfying physical interpretation. However, the derivation has the weakness of making an "ad hoc" use of the macroscopic fields. Force densities resulting from an energy argument have the advantage of dealing rigorously with the macroscopic fields. The form illustrated by (14) also has the advantage of concentrating the magnetization force density at interfaces. For example, it is then clear that the height to which the liquid rises in the experiment of Fig. 11.9.4 depends on the field intensity at the interface and not on the details of the fringing field!

11.12 Macroscopic Force Densities

A macroscopic force density F(r ) is the force per unit volume acting on a medium in the neighborhood of r. Fundamentally, the electromagnetic force density is the result of forces acting on those microscopic particles embedded in the material that are charged, or that have electric or magnetic dipole moments. The forces acting on these individual particles are passed along through interparticle forces to the macroscopic material as a whole. In the limit where that volume becomes small, the force density can then be regarded as the sum of the microscopic forces over a volume element V.

equation GIF #11.314
Of course, the linear dimensions of V are large compared to the microscopic scale.

Strictly, the forces in this sum should be evaluated using the microscopic fields. However, we can gain insight concerning the form taken by the force density by using the macroscopic fields in this evaluation. This is the basis for the following discussions of the force densities associated with unpaired charges and with conduction currents (the Lorentz force density) and with the polarization and magnetization of media (the Kelvin Force Density). To be certain that the usage of macroscopic fields in describing the force densities is consistent with that implicit in the constitutive laws already introduced to describe conduction, polarization and magnetization, the electromagnetic force densities should be derived using energy arguments. These derivations are extensions of those of Secs. 11.6 and 11.7 for forces. We end this section with a discussion of the results of such derivations and of circumstances under which they will predict the same total forces or even material deformations as those derived here.

The Lorentz Force Density

Without restricting the generality of the resulting force density, suppose that the electrical force on a material is due to two species of charged particles. One has N+ particles per unit volume, each with a charge q+, while the other has density N- and a charge equal to -q-. With v denoting the velocity of the macroscopic material and v\pm representing the respective velocities of the carriers relative to that material, the Lorentz force law gives the force on the individual particles.

equation GIF #11.315
equation GIF #11.316
Note that q- is a positive number.

In typical solids and fluids, the charged particles are either bonded to the material or migrate relative to the material, suffering many collisions with the neutral material during times of interest. In either case, the inertia of the particles is inconsequential so that on the average the forces on the individual particles is passed along to the macroscopic material. In either case, the force density on the material is the sum of (2) and (3) respectively multiplied by the charged particle densities.

equation GIF #11.317
Substitution of (2) and (3) into this expression gives the Lorentz force density.

boxed equation GIF #11.41
where u is the unpaired charge density (7.1.6) and J is the current density.

equation GIF #11.318
Because the material is in motion, with velocity v, the current density J has not only the contribution familiar from Sec. 7.1 (7.1.4) due to the migration of the carriers relative to the material, but one due to the net charge carried by the moving material as well.

In EQS systems, the first term in (5) usually outweighs the second, while in MQS systems (where the unpaired charge density is negligible) the second term tends to dominate.

The derivation and Fig. 11.9.1 suggest why the electric term is proportional to the net charge density. In a given region, the force density resulting from the positively charged particles tends to be canceled by that due to the negatively charged particles and the net force density is therefore proportional to the difference in absolute magnitudes of the charge densities. We exploited this fact in Chap. 7 to let electrically induced material motions evidence the distribution of the unpaired charge density. For example, in Demonstration 7.5.1, the unpaired charge density was restricted to an interface and as a result the motion of the fluid was suppressed by constraining the interface. A more recent example is the force on the upper electrode in the capacitor transducer of Example 11.6.1. Here again, the force density is confined to a thin region on the surface of the conducting electrode.

floating figure GIF #60
The magnetic term in (6), pictured in Fig. 11.9.2 as acting on a current carrying wire, is also familiar. It was this force density that was responsible for throwing the metal disk into the air in the experiment described in Sec. 10.2. The force responsible for the levitation of the pan-cake coil in Demonstration 11.7.1 was also the net effect of the Lorentz force density, either acting over the volume of the coil conductors, or over that of the conducting sheet below. In MQS systems, where the contribution of the "convection" current uv is negligible, the current density is typically due to conduction. Note that this means that the velocity of the charge carriers is determined by the electric field they experience in the conductor, and not simply by the motion of the conductor. The current density J in a moving conductor is generally not in the direction of motion.


9 Indeed, it is fortunate that the carriers do not have the same velocity as the material, for if they did it would not be possible to use the magnetic Lorentz force density for electromechanical energy conversion. If we recognize that the rate at which a force f does work on a particle that moves at the velocity v is v f then it follows from the Lorentz force, (1.1.1), that the rate of doing work on individual particles through the agent of the magnetic field is v (v x o H). The cross-product is perpendicular to v, so this rate of doing work must be zero.

The Kelvin Polarization Force Density

If microscopic particles carrying a net charge were the only contributors to a macroscopic force density, it would not be possible to explain the forces on polarized materials that are free of unpaired charge. Example 11.6.2 and Demonstration 11.6.2 highlighted the polarization force. The experiment was carried out in such a way that the dielectric material did not support unpaired charge, so the force is not explained by the Lorentz force density.

In these EQS cases where u = 0, the macroscopic force density is the result of forces on microscopic particles that have dipole moments. The resulting force density is fundamentally different from that due to unpaired charges, the forces p E on the individual microscopic particles are passed along by interparticle forces to the medium as a whole. A comparison of Fig. 11.9.3 to Fig. 11.9.1 emphasizes this point. For a single species of particles, the force density is the force on a single dipole multiplied by the number of dipoles per unit volume Np. By definition, the polarization density P = Np p, so it follows that the force density due to polarization is

boxed equation GIF #11.42
This is often called the Kelvin polarization force density.

floating figure GIF #61

Example 11.12.1. Force on a Dielectric Material

In Fig. 11.9.4, the cross-section of a pair of electrodes that are dipped into a liquid dielectric is shown. The picture might be of a cross-section from the experiment of Demonstration 11.6.2. With the application of a potential difference to the electrodes, the dielectric rises between the electrodes. According to (7), what is the distribution of force density causing this rise?
floating figure GIF #62
For the liquid dielectric, the polarization constitutive law is taken as linear (6.4.2) and (6.4.4)
equation GIF #11.319
so that with the understanding that is a function of position (uniform in the liquid, o in the gas and taking a step at the interface) the force density of (7) becomes
equation GIF #11.320
By using a vector identity

A A
= ( x A) x A + fraction GIF #34 \nabla (A A)
and invoking the EQS approximation where x E = 0, this expression is written as
equation GIF #11.321
A second vector identity

A ) = A \Psi + \Psi \nabla \cdot A converts this expression into one that will now prove useful in picturing the distribution of force density.

equation GIF #11.322
Provided that the interface is well removed from the fringing fields at the top and bottom edges of the electrodes, the electric field is uniform not only in the dielectric and gas above and below the interface between the electrodes, but through the interface as well. Thus, throughout the region between the electrodes there is no gradient of E, and hence, according to (7), no Kelvin force density. The Kelvin force density is therefore confined to the fringing field region where the fluid surrounds the lower edges of the electrodes. In this region, is uniform, so the force density reduces to the first term in (11). Expressed by this term, the direction and magnitude of the force density is determined by the gradient of the scalar E E. Thus, where E is varying in the fringing field, it is directed generally upward and into the region of greater field intensity as suggested by Fig. 11.9.4. The force on the dipole shown by the inset lends further credence to the dipolar origins of the force density.

Although there is no physical basis for doing so, it might seem reasonable to take the force density caused by polarization as being p E. After all, it is the polarization charge density p that was used in Chap. 6 to represent the effect of the media on the macroscopic electric field intensity E. The experiment of Demonstration 11.6.2, pictured in Fig. 11.9.7, is classic because it makes it clear that this force density is not correct. With the interface well removed from the fringing fields, there is no polarization charge density anywhere in the liquid, either at the interfaces or in the fringing field! If pE were the correct force density, it would be zero throughout the fluid volume except at the interfaces with the conducting electrodes. Such a force distribution could not cause the fluid to rise.

The Kelvin Magnetization Force Density

Probably forces caused by magnetization are the most commonly experienced of electromagnetic forces. They account for the attraction between a magnet and a piece of iron. In Example 11.7.1, it is this force density that acts on the disk of magnetizable material.

Given that the magnetizable material is made up of microscopic dipoles, each experiencing a force of the nature of (11.8.28), and that the magnetization density M is the number of these per unit volume multiplied by m, it follows that the force density due to magnetization is

equation GIF #11.323
This is sometimes called the Kelvin magnetization force density.

Example 11.12.2. Force Density in a Magnetized Fluid

With the dielectric liquid replaced by a ferrofluid having a uniform permeability and the electrodes replaced by the pole faces of an electromagnet, the physical configuration shown in Fig. 11.9.4 becomes the one of Fig. 11.9.5, illustrating the magnetization force density. In such fluids[1], the magnetization results from an essentially permanent suspension of magnetized particles. Each particle comprises a magnetic dipole and passes its force on to the liquid medium in which it is suspended. Provided that the magnetization obeys a linear law, the discussion of the distribution of force density given in Example 11.9.1 applies equally well here.

floating figure GIF #63

Alternative Force Densities

We now return to comments used to introduce this section. The fields used to express the Lorentz and Kelvin force densities are macroscopic. To assure consistency between the averages implied by these force densities and those already inherent to the constitutive laws, an energy principle can be used. The approach is a continuum version of that exemplified for lumped-parameter systems in Secs. 11.7 and 11.8. In the lumped parameter systems, electrical terminal relations were used to determine a total energy and energy conservation was used to determine the force in turn. In the continuum system[2], the constitutive laws are used to find an energy density and energy conservation used in turn to find a force density. This energy method, like the one exemplified in Secs. 11.7 and 11.8 for lumped parameter systems, describes systems that are loss-free. In making practical use of the result it is assumed that it will be applicable even if there are losses. A more general method, which invokes a principle of virtual power[3], allows for dissipation but requires more empirical information than the polarization or magnetization constitutive law as a starting point.

Force densities derived from more rigorous arguments then given here can have very different distributions from the superposition of the Lorentz and Kelvin force densities. We expect that the arguments used here prove inadequate when the microscopic particles become so densely packed that the field experienced by one is significantly altered by its nearest neighbor. But surely the difference between the magnetic force density of Lorentz and Kelvin

equation GIF #11.324
we have derived here and the Korteweg-Helmholtz force density for incompressible media
equation GIF #11.325
is not due to interactions between microscopic particles. This latter force density is often obtained for an incompressible material by energy arguments. (Note that, with - and H H respectively playing the roles of dL/d and i2, the magnetization term in (14) takes a form found for the force on a magnetizable material in Sec. 11.7.) In Example 11.9.2 (where J = 0), we found the force density of (13) to be confined to the fringing field. By contrast, (14) gives no force density in the fringing region (where is uniform) but rather puts it all at the interface. According to this latter equation, through the agent of a surface force density (a force density that is a spatial impulse at the interface) the field pulls upward on the interface.

Even though the force densities of (13) and (14) have very different distributions, they predict the same height of rise of the liquid! This is because the liquid deformations being considered are essentially volume conserving (incompressible), in the sense that, (with denoting the liquid displacement) = 0


12 Solids can also deform in an essentially incompressible fashion. An example is the low frequency motion of jello or muscle.. In the force equation representing an incompressible material, there will always be another force density taking the form p, where the pressure p assumes whatever distribution it must to insure that the deformations are incompressible. As a result, contributions to the force density that take the form will have no effect on the incompressible deformations. The contribution of is simply balanced by that due to p.

For a magnetically linear material (o M = ( - \mo )H) the force densities of (13) and (14 do indeed only differ by a term taking the form . To see this, use a vector identity


13 A
A = (\nabla x A) x A + fraction GIF #35 A
A to write (13) as
equation GIF #11.326
The MQS form of Ampère's law makes it possible to substitute
J for x H in this expression, which then becomes
equation GIF #11.327
The second term in this expression is then expanded using a second vector identity

14 x (\Psi A
) = \Psi \nabla A + A \cdot \Psi
equation GIF #11.328
This expression differs from (14) by the last term, which indeed takes the form
where
equation GIF #11.329
Thus, although they have very different distributions, the force densities of (13) and (14) will predict the same incompressible deformations of a material. Will the integration of these two force densities over the volume of an object result in the same total force? Provided that the object is surrounded by media having the properties of free space, the answer is yes. This is because the integral over the volume of the force density by which the two differ is zero. To see this, consider the
x component (say) of that integral.

equation GIF #11.330
With the objective of converting this volume integral to one over the enclosing surface, we write this expression in the equivalent form
equation GIF #11.331
It follows from Gauss' Theorem that the volume integral on
V is equivalent to one over the enclosing surface S.

equation GIF #11.332
With the surface
S taken as enclosing an object surrounded by free space, = o on S and is zero everywhere on S. We conclude that integration of either (13) or (14) over the volume will give the same total force. In summary, if an object is surrounded by free space, integration over its volume of two force densities that differ by the gradient of a scalar that is zero in free space will result in the same net force.

Example 11.12.3. Magnetic Force on a Magnetizable Current Carrying Material

A block of conducting material having permeability is shown in Fig. 11.9.6 sandwiched between perfectly conducting plates. A current source, distributed over the left edges of these electrodes, drives a constant surface current density K in the +x direction along the left edge of the lower electrode. This current passes through the block in the y direction as a current density
floating figure GIF #64
equation GIF #11.333
and is returned to the source in the
-x direction at the left edge of the upper electrode. The thickness a of the block is small compared to its other two dimensions, so the magnetic field between the electrodes is z directed and dependent only on x. From Ampère's law it follows that
equation GIF #11.334
in the conducting block.

The alternative force densities, (13) and (14), have very different distributions in the block. Yet, we must find that the net force on the block found by integrating each over its volume is the same. To see that this is so consider first the sum of the Lorentz and Kelvin force densities, (13).

There is no x component of the magnetic field intensity, so for this particular configuration, the magnetization term makes no contribution to (13). Evaluation of the first term using (22) and (23) then gives

equation GIF #11.335
Integration of this force density over the volume amounts to a multiplication by the cross-sectional area ab and integration on x.

equation GIF #11.336
Now, the force density given by (14) is evaluated. The permeability
is uniform throughout the interior of the block, so the magnetization term is again zero there. However, is a step function at the ends of the block, where x = -b and x = 0. Thus, is an impulse there and we must take care to include the contributions from the surface regions in our integration. Evaluation of the x component of (14) using (22) and (23) gives
equation GIF #11.337
Note that
Hz is K at x = -b and is zero at x = 0. Integration of (26) over the volume of the block therefore gives

equation GIF #11.338
Note that
Hz is constant through the interface at x = -b. So, the integration of the last term can be carried out. Simplification of this expression gives the same total force as found before, (25).

The distributions of the force densities given by (13) and (14) are generally different, even very different. It is natural to therefore ask which of the two is the "right" one. In general, until the "other" force densities acting on the medium in question are specified, this question cannot be answered. Here, where a discussion of continuum mechanics is beyond our purvue, we have identified a class of mechanical deformations (namely those that are volume conserving or "incompressible"), where these force densities are equally valid. In fact so would any other force density differing from these by a term having the form . The combined Lorentz and Kelvin force densities have the advantage of a satisfying physical interpretation. However, the derivation has the weakness of making an "ad hoc" use of the macroscopic fields. Force densities resulting from an energy argument have the advantage of dealing rigorously with the macroscopic fields. The form illustrated by (14) also has the advantage of concentrating the magnetization force density at interfaces. For example, it is then clear that the height to which the liquid rises in the experiment of Fig. 11.9.4 depends on the field intensity at the interface and not on the details of the fringing field!

R E F E R E N C E S

  1. R. E. Rosensweig, "Magnetic fluids," Scientific American, 136-145, Oct. 1982.

  2. J. R. Melcher, Continuum Electromechanics, MIT Press, 1981, Chap. 3.
  3. P. Penfield and H. A. Haus, Electrodynamics of Moving Media, MIT Press, Cambridge, Mass., 1967.

11.13 Summary

The basis for a self-consistent macroscopic description of any continuum subsystem is a power-flow statement having the forms identified in Sec. 11.1. Describing a volume V enclosed by a surface S, the integral conservation of energy statements takes the form (11.1.1)
equation GIF #11.339
Implied by this is the differential form of the conservation of energy statement.
equation GIF #11.340
Poynting's Theorem, the subject of Sec. 11.2, is obtained starting from the laws of Faraday and Ampère to obtain an expression of the form of (2). For materials that are Ohmic
(J = E) and that are linearly polarizable and magnetizable (D = E and B = H) the power flux density S (or Poynting's vector), energy density W and power dissipation density Pd were shown in Sec. 11.3 to be
equation GIF #11.341
equation GIF #11.342
equation GIF #11.343
Of course, taking the free space limit where
and assume their free-space values and = 0 gives the free-space conservation statement discussed in Sec. 11.2.

In Sec. 11.3, we found that in EQS systems an alternative to Poynting's vector is (11.3.24)

equation GIF #11.344
This expression is of practical importance because it can be evaluated without determining
H, which is generally not of interest in EQS systems. Perhaps a more important reason for introducing the alternative EQS forms of S is that it emphasizes that S is not unique. In Sec. 11.3 we saw that Poynting's flux density and (6) assumed entirely different spatial distributions. Yet, evaluation of the integral of S da over a closed surface (the left hand side of (1)) and of the divergence of S at any point must give the same result using either expression. The integral over a closed surface of S da or the divergence of S is the physically significant quantity and not the value of S at a given point.

An important application of the integral form of the energy conservation statement is to lumped parameter systems. In these cases, the surface S of (1) encloses a system that is connected to the outside world through terminals. It is then convenient to describe the power flow in terms of the terminal variables. It was shown in Sec. 11.3 (11.3.29) that the net power into the system represented by the left hand side of (1) becomes

equation GIF #11.345
provided that the magnetic induction is negligible on the surface S, and the electric displacement current is parallel to the surface S. This set the stage for the application of the integral form of the energy conservation theorem to lumped parameter systems.

In Sec. 11.4, attention focused on the energy storage term, the first terms on the right in (1) and (2). The energy density concept was broadened to include materials having constitutive laws relating the flux densities to the field intensities that were single valued and collinear. With E, D, H and B representing the field magnitudes, the energy density was found to be the sum of electric and magnetic energy densities.

equation GIF #11.346
Integrated over the volume V of a system, this function leads to the total energy w. For quasistatic lumped parameter systems, the total electric or magnetic energy is often conveniently found following a different route. First, the terminal relations are determined and then the total energy is found by adding up the increments of energy put into the system as it is energized. In the case of an n terminal pair EQS system, where the relation between terminal voltage vi and associated charge qi is vi (q1, q2, \ldots qn ), the increment of energy is vi dqi and the total electric energy is (11.4.9)
equation GIF #11.347
The line integration in an
n-dimensional space representing the n independent qi's was illustrated by Example 11.4.2. Similarly, for an n terminal pair MQS system where the current ii is related to the flux linkage i by ii = ii (1, 2, \ldots n), the total energy is (11.4.12)
equation GIF #11.348
Note the analogy between these expressions for the total energy of EQS and MQS lumped parameter systems and the respective electric and magnetic energy densities of (8). The transition from the field picture afforded by the energy densities to the lumped parameter characterization is made by
E v, D \rightarrow q and by H i, B \rightarrow .

Especially in using the energy to evaluate forces of electrical origin, we found it convenient to define coenergy density functions.

equation GIF #11.349
It followed that these functions were natural when it was desirable to use
E and H as the independent variables rather than D and B.

equation GIF #11.350
The total coenergy functions for lumped parameter EQS and MQS systems could be found either by integrating these densities over the volume or by again viewing the system in terms of its terminal variables. With the total coenergy functions defined by
equation GIF #11.351
it followed that the coenergy functions could be determined from the terminal relations by again carrying out line integrations, but this time with the voltages and currents as the independent variables. For EQS systems,

equation GIF #11.352
while for MQS systems
equation GIF #11.353
Again, note the analogy to the respective terms in (12).

The remaining sections of the chapter developed some of the possible implications of the "dissipation" term in the energy conservation statement, the last terms in (1) and (2). In Sec. 11.5, coupling to a thermal subsystem was the point. In this section, the disparity between the power input and the rate of increase of the energy stored was accounted for by heating. In addition to Ohmic heating, caused by collisions between the migrating carriers and the neutral media, we considered losses associated with the dynamic polarization and magnetization of materials.

In Secs. 11.6-11.9 we considered coupling to a mechanical subsystem as a second mechanism by which energy gain or energy loss from the electromagnetic system could occur. With the displacement of an object denoted by , we used an energy conservation postulate to infer the total electric or magnetic force acting on the object from the energy functions ((11.6.9), and its magnetic analog)

equation GIF #11.354
or from the coenergy functions ((11.7.7) and the analogous expression for electric systems).

equation GIF #11.355
In Sec. 11.8, where the Lorentz force on a particle was generalized to account for electric and magnetic dipole moments, one objective was a microscopic picture that would lend physical insight to the forces on polarized and magnetized materials. The Lorentz force was generalized to indicate the force on stationary electric and magnetic dipoles, respectively
equation GIF #11.356
The total macroscopic forces resulting from microscopic forces had already been encountered in the previous two sections. The force density describes the interaction between a volume element of the electromagnetic subsystem and a mechanical continuum. A rigorous approach to finding the force density is based on a generalization of the energy method introduced in Secs. 11.6 and 11.7.

Such an approach has to include information on the effect of thermodynamic variables (such as density and entropy) on the constitutive laws[1,2].

R E F E R E N C E S

  1. P. Penfield, Jr., and H. A. Haus, Electrodynamics of Moving Media, MIT Press, Cambridge, MA, 1967.

  2. J. R. Melcher, Continuum Electromechanics, MIT Press, Cambridge, MA, 1981, Chap. 3.