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,
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.
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 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
|
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.
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 Sda. With these field-theoretical generalizations, the power flowing into a volume V, enclosed by the surface S must be given by
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.
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. Because V is arbitrary, the integrand must be zero and a differential statement of energy conservation follows.
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.
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).
Faraday's law, including the effects of magnetization, is (11.0.3).
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
The rate of work on the particle is If the particle density is N and only one species of charged particles exists, then the rate of work per unit volume is 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,
Thus, subtracting (2) dot multiplied by H from (1) dot multiplied by E one obtains In writing the first and third terms on the right, we have exploited the relation u du = d( 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 becomeswith and 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),A second example offers the opportunity to apply the integral version of (8) to a simple MQS system.while the magnetic field that follows from the integral form of Ampère's law is (3.3.10).
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).
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.
where 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. 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), Comparison of this expression with the electric energy storage found in (15) shows that the EQS approximation is valid provided that For a sinusoidal excitation of frequency , this gives 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.
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),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].
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 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. where 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.
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).
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. 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.
Consider a stationary material described by the constitutive lawswhere 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 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
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
In recognition of the balance between collision forces and electrical forces, the forces of (4) are replaced by |q+|E and -|q-|E, respectively.
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 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.\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,First, within the circular cylindrical conductor [region (b) in Fig. 11.3.1], the electric field was found to be uniform, (7.5.7),
while in the surrounding free space region, it was [from (7.5.11)] and in the disk-shaped conductor [from (7.5.9)] 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, 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.
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 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 where i is the total current through the cylinder, or from an evaluation of the right-hand side of the integral conservation law.
and the power input term on the left in the integral conservation law, (11.1.1), can be expressed as Next, the vector identity is used to write the right-hand side of (19) as 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. In this way, we have determined an alternative expression for S, valid only in the electroquasistatic approximation.
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 isBy 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)] and it follows that the power flux density is simply 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. 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?
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). 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 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,
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.
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.
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 volumeThe 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 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
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
where 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
where 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.
Further, w = we, the stored electric energy. Thus, one concludes from (6) that 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 To complete this integral, each of the terminal voltages must be a known function of the associated charges.
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.
Then, (6) leads to To complete this integral, we require the terminal currents as functions of the terminal flux linkages.
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.
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).
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 Note that this expression takes the form where In terms of the terminal variables, where q = Cv, the total energy follows from an evaluation of (9).
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 relationsIn this case, (12) becomes 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 and regroup terms where the coenergy is defined as 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).
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 Substitution for the flux linkages from (18) and evaluation of the integrals then gives 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.
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
Evaluation of the energy stored in a unity-coupled transformer, where the inductances take the form of (9.7.20), gives 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 formHere, 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
and then regroup terms so that the expression becomes where 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
It follows from (31) that the energy density is 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).
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
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.
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.
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 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 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,
This can be shown by using the identity 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
According to (10.3.9), the complex amplitude of the response, the magnetic field inside the shell, is where m = o a. The complex amplitude of the surface current density circulating in the shell follows from (10.3.8).
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.
It follows from the application of the time average theorem, (6), that the total time average dissipation is 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).
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. To evaluate this expression, (10) is used to determine E Evaluation of (14) then gives 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.
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 averageIn view of the time average theorem, (6), this becomes 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. Here \hat is the complex permittivity with real and imaginary parts ' and -", respectively. Evaluation of (18) using this constitutive law gives Thus, " represents the electrical dissipation associated with the polarization process. 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, From this definition, it follows from Eulers formula that Given the complex amplitude of the electric field, D is 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
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
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 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 The time average power dissipation density follows from this expression and (20).
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
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) thatThis 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 The time average power dissipation, (29), then is the sum of these two contributions divided by T.
Thus, the area within the hysteresis loop is the energy dissipated in one cycle.
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. Thenwhere 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 and 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.
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
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 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.
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 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
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.
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. Such a path
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
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?
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 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 Elimination of Eb between these equations gives 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 and so it follows from (18) that the electrical terminal relation can be written in terms of a capacitance C. 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).
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.
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 .With the letters A through D used to designate states of the system, we consider the following energy conversion cycle. Is electrical energy converted to mechanical form, or vice versa?
\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. 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).
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.
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 .
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
In view of the expression for C( ), (20), this expression can be written as 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.
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.
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).
It follows from (15) that there is a force tending to draw the dielectric into the region between the electrodes.
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.
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.
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.
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.
This expression can be solved for (r).
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.
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.
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. 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
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
and substitute into (2) to obtain where the coenergy function, seen before in Sec. 11.4, is defined as 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 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, 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 thatAlthough 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.
Application of (8) shows that energy conservation requires the equality of the mutual inductances. 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.
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 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).
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.
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.
A second relation between Ha and Hc follows from requiring that the net flux out of the disk must be zero.
Using this last expression to replace Hc in (14) results in The magnetic flux linking every turn in the coil is o Hc c2. Thus, the total flux linked by the coil is
Finally, the force follows from an evaluation of (13) (specialized to the single terminal pair system of this example).
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?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.
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).
The force per unit length in the direction then follows from an evaluation of (13) (again adapted to the single terminal pair situation).
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 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).
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.
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 becomeswhere is defined as the torque.
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 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 densitieswhere 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 formwhere 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).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 thatUnder 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 A trigonometric identity
1\over 2 (sin (x + y) + sin (x - y)) makes the implications of this result more apparent.
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, the time average torque is 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.
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.
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.
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.
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.
Force on an Electric Dipole
The force on a stationary electric charge is given by the Lorentz law with v = 0.
A dipole is the limit of two charges of equal magnitude and opposite sign spaced a distance d apart, in the limit 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. 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
The first and last terms cancel. In more compact notation, this expression is therefore 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 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 fieldwhich is familiar from Example 4.1.1. It follows from (5) that the force is 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.
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.
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
(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.
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 where stands for the differential change of o E. The change of energy is where V includes all of space. The electroquasistatic E field is the negative gradient of the potential Introducing this into (10), one has 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
from Gauss' law, where u is the change of unpaired charge. Thus, from (12), 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 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 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
and therefore 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 iswhere, for the uniformly magnetized rod, and the magnetic charge 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), 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), where m is the magnetic dipole moment. 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 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
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,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 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.
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 Directly below the magnet, H has only a z component. Thus, the dipole moment follows from (27) as Evaluation of (22) therefore gives 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.
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
- P. Penfield, Jr., and H. A. Haus, Electrodynamics of Moving Media, MIT Press, Cambridge, Mass. (1967).
- R. M. Fano, L. J. Chu, and R. B. Adler, Electromagnetic Fields, Energy, and Forces, John Wiley and Sons, New York (1960).
- D. B. H. Tellegen, "Magnetic-dipole models," Am. J. Phys. Vol. 30 (Sept. 1962), pp. 650-652.
- H. A. Haus and P. Penfield, Jr., "Force on a current loop," Phys. Lett., Vol. 26A (March 1968), pp. 412-413.
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.
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.
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.
Substitution of (2) and (3) into this expression gives the Lorentz force density where u is the unpaired charge density (7.1.6) and J is the current density.
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.
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
This is often called the Kelvin polarization force density.
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?For the liquid dielectric, the polarization constitutive law is taken as linear [(6.4.2) and (6.4.4)] 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 By using a vector identity
A A= ( x A) x A + \nabla (A A) and invoking the EQS approximation where x E = 0, this expression is written asA second identity
11 ( ) = + \nablaconverts this expression into one that will now prove useful in picturing the distribution of force density.
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
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.
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)
we have derived here and the Korteweg-Helmholtz force density (KH) for incompressible media 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 + (A A) to write (13) asThe MQS form of Ampère's law makes it possible to substitute x H for J in this expression, which then becomes The second term in this expression is then expanded using a second vector identity
13 ( ) = \nabla +This expression differs from (14) by the last term, which indeed takes the form where Now consider Newton's force law for an elemental volume of material. Using the Korteweg-Helmholtz force density, (14), it takes the form
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
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 densityand 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 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
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 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
Integration of (25) over the volume of the block therefore gives 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.
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).
The differential form of the conservation of energy statement is implied by the above.
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 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).
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 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.
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).
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).
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.
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.
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 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,
while for MQS systems, 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]
or from the coenergy functions [(11.7.7) and the analogous expression for electric systems].
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.
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 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].
| |||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Introduction
| |||||||||||||||
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
| ||||||||||||||
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.
| ||||||||||||||
11.2.2 | The 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.
| ||||||||||||||
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).
| ||||||||||||||
11.3.2 | Coaxial 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.
| ||||||||||||||
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.
| ||||||||||||||
11.3.4 | In 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).
| ||||||||||||||
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 = Li2 to show that L is as given by
(8.5.20).
| ||||||||||||||
11.4.2 | In 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 =
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.4 | The 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.5 | The cross-section of a motor or generator is shown in Fig.
11.7.7.
the total coenergy wm'.
| ||||||||||||||
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
| ||||||||||||||
11.4.7 | Consider 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
| ||||||||||||||
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
| ||||||||||||||
11.5.2 | In 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.
| ||||||||||||||
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.
| ||||||||||||||
11.5.4 | In 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.
| ||||||||||||||
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
| ||||||||||||||
11.5.6 | In 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
| ||||||||||||||
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.8 | The conducting block shown in Fig. 10.7.2 has a length d in the z
direction.
| ||||||||||||||
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.10 | The 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).
| ||||||||||||||
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.
| ||||||||||||||
11.6.2 | In 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.
| ||||||||||||||
11.6.4 | The 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 .
| ||||||||||||||
11.6.5* | Show that the vertical force on the nonlinear dielectric material
of Prob. 11.4.6 is
| ||||||||||||||
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
| ||||||||||||||
11.7.2 | Determine the force f(i, ) acting in the x direction on the
plunger of the magnetic circuit shown in Fig. P9.7.6.
| ||||||||||||||
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
| ||||||||||||||
11.7.4 | The 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.
| ||||||||||||||
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/.
| ||||||||||||||
11.7.6 | A "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
| ||||||||||||||
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).
| ||||||||||||||
11.8.2 | The 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.
| ||||||||||||||
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
| ||||||||||||||
11.8.4 | A 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.
| ||||||||||||||
Macroscopic Force Densities | |||||||||||||||
11.9.1 | In 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.)
| ||||||||||||||
11.9.3 | In 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.)
|
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
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 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 densityand 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 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
Integration of this force density over the volume amounts to a multiplication by the cross-sectional area ab and integration on x.
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 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
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
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.
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.
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.
Substitution of (2) and (3) into this expression gives the Lorentz force density.
where u is the unpaired charge density (7.1.6) and J is the current density.
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.
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
This is often called the Kelvin polarization force density.
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?For the liquid dielectric, the polarization constitutive law is taken as linear (6.4.2) 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 By using a vector identity \footnote* and invoking the EQS approximation where x E = 0, this expression is written asA second vector identity \footnote** converts this expression into one that will now prove useful in picturing the distribution of force density.
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
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.
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
we have derived here and the Korteweg-Helmholtz force density for incompressible media 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 + A A to write (13) asThe MQS form of Ampère's law makes it possible to substitute J for x H in this expression, which then becomes The second term in this expression is then expanded using a second vector identity \footnote* x (\Psi A) = \Psi \nabla A + A \cdot This expression differs from (14) by the last term, which indeed takes the form where 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.
With the objective of converting this volume integral to one over the enclosing surface, we write this expression in the equivalent form It follows from Gauss' Theorem that the volume integral on V is equivalent to one over the enclosing surface S.
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 densityand 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 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
Integration of this force density over the volume amounts to a multiplication by the cross-sectional area ab and integration on x.
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 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
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!
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.
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.
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.
Substitution of (2) and (3) into this expression gives the Lorentz force density.
where u is the unpaired charge density (7.1.6) and J is the current density.
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.
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.
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
This is often called the Kelvin polarization force density.
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?For the liquid dielectric, the polarization constitutive law is taken as linear (6.4.2) and (6.4.4) 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 By using a vector identity
A A= ( x A) x A + \nabla (A A)and invoking the EQS approximation where x E = 0, this expression is written asA 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.
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
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.
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
we have derived here and the Korteweg-Helmholtz force density for incompressible media 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
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 = (\nabla x A) x A + A
A to write (13) asThe MQS form of Ampère's law makes it possible to substitute J for x H in this expression, which then becomes The second term in this expression is then expanded using a second vector identity
) = \Psi \nabla A + A \cdot \PsiThis expression differs from (14) by the last term, which indeed takes the form where 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.
With the objective of converting this volume integral to one over the enclosing surface, we write this expression in the equivalent form It follows from Gauss' Theorem that the volume integral on V is equivalent to one over the enclosing surface S.
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 permeabilityis 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 densityand 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 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
Integration of this force density over the volume amounts to a multiplication by the cross-sectional area ab and integration on x.
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 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
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!
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)Implied by this is the differential form of the conservation of energy statement. 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 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)
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
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.
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) 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) 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.
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.
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 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,
while for MQS systems 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)
or from the coenergy functions ((11.7.7) and the analogous expression for electric systems).
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 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].