Electrostatic and Magnetostatic Forces That Arise from Electrostatic and Magnetostatic Pressures

Show more

1. Introduction

In the interaction of electromagnetic fields with matter, there appear forces that theory must explain in order to have some control on these forces. Perhaps the simplest of these are the forces that arise from the interaction of electrostatic and magnetostatic fields with polarizable and magnetizable matter.

The electrostatic case may be illustrated with the force exerted on a dielectric slab partially introduced into a charged parallel plate capacitor. In the magnetostatic case, the force that is exerted on a magnetizable bar inside a solenoid, is more familiar and has many applications. It is generally considered that this magnetic effect is analogous to the electrostatic case. We show in this paper that this is not the case, and present a magnetic device really analogous to the electrostatic case.

In the case of a dielectric slab inside a capacitor, the force is usually explained as the action of the non-uniform fringing electrostatic field on the electric dipoles of the dielectric. We have shown elsewhere [1] that this force arises rather from the action of Maxwell’s electrostatic stresses at the dielectric-vacuum interface. We have also shown that the above magnetic effect arises from the magnetostatic stresses, but from the tension part of the stress [2] . This interpretation is based on the original conception of Faraday and Maxwell that the electromagnetic forces are transmitted through stresses in matter and vacuum. This view is expressed clearly by Maxwell [3] :

“If we further admit that every part of a dielectric medium through which electric induction is taking place there is a tension, like that of a rope, in the direction of the lines of force, and a pressure in all directions at right angles to the lines of force, we may account for all the mechanical actions which take place between electrified bodies.”

Indeed, we have shown elsewhere [1] [2] that in the electric case the force has its origin in the compression around the lines of force, while in the magnetic case the force has its origin in the tension.

2. Theory

In formal terms it is shown in texts [4] [5] [6] that Maxwell’s stress tensor can be decomposed in principal axes, with components given by a tension

$t=\frac{1}{2}\u0454{E}^{2}$ (1)

along the lines of force, and an orthogonal compression

$K=-\frac{1}{2}\u0454{E}^{2}$ (2)

around the lines of force where $\u0454$ is the permittivity and $E$ is the electric field. In spite that these facts are discussed in well-known texts, some authors [7] [8] [9] [10] [11] , express surprise at the fact that an electrostatic field can exert forces in a direction orthogonal to the lines of force.

We propose to analyze these forces from the point of view of balance equations that can be derived from Maxwell’s equations and therefore are well founded. These balance equations have the structure

$\nabla \cdot \stackrel{\leftrightarrow}{T}-{\partial}_{t}g=f$ (3)

where $\stackrel{\leftrightarrow}{T}$ is Maxwell’s stress tensor, $g$ is a momentum density related to Poynting’s vector, and $f$ is a force density. In static conditions this balance equation permits to obtain the force as a volume integration of the force density, or as a surface integral of the stress tensor over a surface surrounding the volume of interest. We have done these calculations [1] [2] , for the cases of the capacitor and solenoid and we find that they are not exactly analogous. In the first case the force arises from the compression around the electrostatic lines of force, while in the second case the force arises rather from the tension of the magnetostatic lines of force.

In the case of the capacitor, the force accepted as correct is [12] [13] ,

$F=\frac{1}{2}{\u0454}_{0}{\chi}_{e}{E}^{2}{A}_{0}\stackrel{^}{k}.$ (4)

Here ${A}_{0}$ is the cross section of the dielectric slab, and $E$ is the uniform electric field inside the capacitor, taken as the potential difference divided by the separation of the plates.

If we use the method of a balance equation in the calculation of this force, the force density adequate to solve the problem is

$f=-\frac{1}{2}E\times \left(\nabla \times P\right).$ (5)

In the case of the solenoid we have that the force on a magnetizable bar partially introduced in it is [12] [13] ,

$F=\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}{A}_{0}\stackrel{^}{k}.$ (6)

Here $H$ is taken as the uniform field inside the solenoid. No reference is made to the fringing field (We follow the convention proposed by Purcell and Morin [14] and Griffiths [11] ) of considering $B$ as the magnetic field and $H$ as an auxiliary field.

The analogy with the capacitor is apparent, with the correspondences

${\mu}_{0}\to {\u0454}_{0},{\chi}_{m}\to {\chi}_{e},H\to E,$ (7)

It is worthwhile to note that in this case the correspondence of fields is

$H\to E$

rather than

$B\to E,$

as frequently proposed [12] [15] .

If the analogy is further pursued, and as force density we propose the analogy to “Equation (5)”

${f}_{newM}=-\frac{1}{2}{\mu}_{0}H\times \left(\nabla \times M\right),$ (8)

we do not obtain the known result. In this case the adequate force density is rather [2]

${f}_{newM}=-\frac{1}{2}{\mu}_{0}\left(H\cdot \nabla \right)M,$ (9)

which is obtained from a magnetic energy density.

It is important to note that this force density is part of a balance equation, and is different from the familiar force density

${f}_{M}=-\frac{1}{2}{\mu}_{0}\left(M\cdot \nabla \right)H,$ (10)

which applies in the case of magnetic dipoles immersed in a non-uniform magnetic field.

The force densities “Equations (5) and (9)” are different, but the forces to which they lead are similar. Then, which is the magnetic system analogous to the capacitor, such that the forces arise from similar force densities? One aim of the present paper is to propose such a magnetic system, another is to explain the origin and relation among different force densities.

The action of stresses across the boundary between two dielectrics may be calculated with the Helmholtz force density [6] , which also appears naturally in a balance equation derived from Maxwell’s equations [16] . However, the force that arises from the action of these stresses can also be calculated with the force density “Equation (5)”; indeed, this force density is equivalent to the Helmholtz force density as consequence of the discontinuity in the permittivity. The action of the stresses can also be seen in the sucking of a liquid dielectric by a laser that falls orthogonally to the surface of the liquid, since in this case we have an average electric field parallel to interface.

3. Analogy between Capacitor with Dielectric Slab and Solenoid with Magnetizable Bar?

Let us consider the space between the poles of a permanent horseshoe magnet. The field can be assumed almost uniform in middle space, while the fringing field is non-uniform. We define the x direction as that of the uniform field and a magnetizable bar of cross section ${A}_{0}$ and length L is partially introduced a distance l (in the z direction), in the region of the uniform field Figure 1.

Figure 1. A magnetizable bar of cross section ${A}_{0}$ is partially introduced a distance l between the poles of a permanent magnet. In the region of length ${L}^{\prime}$ the magnetic field is approximately uniform.

The magnetic energy is

$U=\frac{1}{2}{\displaystyle {\int}_{Vol}^{}\text{d}V}B\cdot H=\frac{1}{2}{\displaystyle {\int}_{Vol}^{}\text{d}V}BH,$ (11)

since the magnetic field $B$ and the $H$ field are parallel (we are considering linear media). If the integration volume is ${A}_{0}{L}^{\prime}$ , where ${L}^{\prime}$ is the length over which $B$ and $H$ are approximately uniform, then the energy is

$U=\frac{1}{2}\left[{\left(BH\right)}_{medium}+{\left(BH\right)}_{vacuum}\right]{A}_{0}{L}^{\prime},$ (12)

since $B$ and $H$ are uniform.

The continuity conditions for the fields are

${B}_{vacuum}={\mu}_{0}{H}_{vacuum}$ (13)

and

${B}_{medium}={\mu}_{0}{\mu}_{r}{H}_{medium}.$ (14)

Then, with the constitutive relation,

${\mu}_{r}=1+{\chi}_{m},$ (15)

we can write “Equation (12)” in the form

$U=\frac{1}{2}\left[{\mu}_{0}{\mu}_{r}{H}^{2}{A}_{0}z+{\mu}_{0}{H}^{2}{A}_{0}\left({L}^{\prime}-z\right)\right],$ (16)

which can be rewritten as

$U=\frac{1}{2}{\mu}_{0}{H}^{2}{A}_{0}\left({\chi}_{m}z+{L}^{\prime}\right),$ (17)

leading to the force on the magnetizable bar

$F=-\frac{\partial U}{\partial l}=\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}{A}_{0}\stackrel{^}{k}.$ (18)

This has the same structure that the known force, and if the analogy expressed in “Equation (7)” is used, “Equation (18)” is transformed into the known force exerted on a dielectric slab inside a capacitor.

4. The Force Density

The magnetic energy density of linear magnetizable matter when introduced into a magnetic field is [12] [13] [17]

$u=-\left(\frac{1}{2}{\mu}_{0}M\cdot H\right).$ (19)

Since the field of the permanent magnet can be considered constant, it is equivalent to a constant current system, and the force density can be calculated with the equation

$f={\left(\nabla u\right)}_{l}.$ (20)

Then,

$f=-\nabla \left(\frac{1}{2}{\mu}_{0}M\cdot H\right),$ (21)

which by means of the vector identity

$\nabla \left(a\cdot b\right)=\left(a\cdot \nabla \right)b+\left(b\cdot \nabla \right)a+a\times \left(\nabla \times b\right)+b\times \left(\nabla \times a\right),$ (22)

is transformed into

$f=-\frac{1}{2}{\mu}_{0}\left[\left(M\cdot \nabla \right)H+\left(H\cdot \nabla \right)M+M\times \left(\nabla \times H\right)+H\times \left(\nabla \times M\right)\right].$ (23)

Since there are not free currents, we have

$\nabla \times H=0.$ (24)

Now, the magnetic susceptibility is different from zero only in matter; therefore the constitutive relation is

$M={\chi}_{m}\Theta \left(l-z\right)H,$ (25)

where ${\chi}_{m}$ is constant and $\Theta $ is the Heaviside distribution.

In this case the $H$ field is in the x direction, so we have

$M\left(z\right)={\chi}_{m}\Theta \left(l-z\right){H}_{x}\stackrel{^}{i}.$ (26)

Therefore,

$\left(H\cdot \nabla \right)M=0.$ (27)

Also, since the interface is in the region where $H$ is uniform,

$\left(M\cdot \nabla \right)H=0.$ (28)

In this way “Equation (23)” reduces to

${f}_{newMag}=-\frac{1}{2}{\mu}_{0}\left[H\times \left(\nabla \times M\right)\right],$ (29)

which is precisely the force density proposed by analogy, “Equation (8)”.

In the next section we show that this force density also leads to the known result obtained by conventional methods, that is, from an energy density.

5. Force on a Magnetizable Bar in the Field of a Permanent Magnet Calculated with the Proposed Force Density

With the constitutive relation “Equation (26)” the proposed new force density “Equation (29)” results

${f}_{newMag}=-\frac{1}{2}{H}_{x}\left[{\partial}_{z}{M}_{x}\left(z\right)\right]\stackrel{^}{k},$ (30)

since the field $H$ is in the x direction.

Now, the mathematical result

${\partial}_{z}\Theta \left(l-z\right)=-\delta \left(l-z\right)$ (31)

permits to write “Equation (30)” in the form

${f}_{newMag}=\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}\delta \left(l-z\right)\stackrel{^}{k}.$ (32)

In order to obtain the force on the magnetizable rod it is necessary to integrate the force density over a volume ${L}^{\prime}h$ around the interface, that is,

$F={\displaystyle {\int}_{vol\text{\hspace{0.17em}}around\text{\hspace{0.17em}}interface}^{}{f}_{newM}\text{d}S\text{d}z}.$ (33)

If the cross section of the interface is ${L}^{\prime}h={A}_{0}$ , then with the force density “Equation (32)” the force results

$F={\displaystyle {\int}_{vol\text{\hspace{0.17em}}around\text{\hspace{0.17em}}interface}^{}\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}h\delta \left(l-z\right)\stackrel{^}{k}\text{d}z.}$ (34)

After integration we obtain the force

$F=\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}{A}_{0}\stackrel{^}{k},$ (35)

which is the result previously obtained with the usual procedure using an energy density, “Equation (18)”. This proves the correctness of the proposed force density.

6. Maxwellian Momentum Balance Equation and the Proposed Force Density

We have shown that the force density “Equation (8)” leads to the same result than the method using the energy density “Equation (18)”. However, the force density “Equation (8)” was proposed by analogy with the case of a capacitor with dielectric. Then, it is necessary to give this force density a solid foundation. This can be done with a momentum balance equation derived directly from Maxwell’s equations. One such equation is [16]

$\begin{array}{l}\nabla \cdot \left\{DE+BH-\frac{1}{2}I\left(D\cdot E+B\cdot H\right)\right\}-{\partial}_{t}\left(D\times B\right)\\ =\rho E+J\times B+\frac{1}{2}\left[\left(\nabla E\right)\cdot D-\left(\nabla D\right)\cdot E+\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right].\end{array}$ (36)

Our aim is to show that this balance equation contains the proposed force density. Since we are dealing with a magnetostatic case there are not free currents, the balance equation “Equation (36)” reduces to

$\nabla \cdot \left\{BH-\frac{1}{2}I\left(B\cdot H\right)\right\}=\frac{1}{2}\left[\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right].$ (37)

We then have to show that the right-hand member of “Equation (37)” contains the proposed force density. This can be done as follows.

The identity

$\left(\nabla v\right)\cdot u=u\times \left(\nabla \times v\right)+\left(u\cdot \nabla \right)v$ (38)

and the constitutive relation

$B={\mu}_{0}\left(H+M\right),$ (39)

permits to express the right-hand member of “Equation (37)” as

$\begin{array}{l}\frac{1}{2}\left[\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right]\\ =\frac{1}{2}{\mu}_{0}\left[M\times \left(\nabla \times H\right)+\left(M\cdot \nabla \right)H-H\times \left(\nabla \times M\right)-\left(H\cdot \nabla \right)M\right].\end{array}$ (40)

In order to see that the force density “Equation (8)” is contained in the right-hand member of “Equation (40)” we must note that there are not free currents and the field $H$ is uniform at the interface, so that “Equations. (24), (27), and (28)” hold, and therefore “Equation (40)” reduces to

$\frac{1}{2}\left[\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right]=-\frac{1}{2}H\times \left(\nabla \times M\right),$ (41)

which is the result looked for.

In this way, the balance equation for these particular conditions is

$\nabla \cdot \left\{BH-\frac{1}{2}I\left(B\cdot H\right)\right\}=-\frac{1}{2}H\times \left(\nabla \times M\right).$ (42)

These arguments show that the proposed force density “Equation (8)” has a firm foundation on Maxwell’s equations. Additional support for the proposed force density can be obtained by a surface integration of the left-hand member of the balance equation “Equation (42)”, that is, of the magnetostatic stress tensor. This is done in the following section.

7. Calculation of the Force by a Surface Integration of the Stress Tensor

The force in terms of the magnetostatic stress tensor is given by

$F={\displaystyle {\oint}_{\sigma}\text{d}S\cdot \stackrel{\leftrightarrow}{T}},$ (43)

where

$\stackrel{\leftrightarrow}{T}=\left[BH-\frac{1}{2}\stackrel{\leftrightarrow}{I}\left(B\cdot H\right)\right]$ (44)

and $\sigma $ is a closed surface. The force is the total magnetic force exerted on the matter enclosed in the surface. In this case a convenient surface is one that surrounds the interface.

The constitutive relation

$B={\mu}_{0}{\mu}_{r}H\stackrel{^}{i}$ (45)

permits to write the tensor $\stackrel{\leftrightarrow}{T}$ as

$\stackrel{\leftrightarrow}{T}={\mu}_{0}{\mu}_{r}{H}^{2}\left[\stackrel{^}{i}\stackrel{^}{i}-\frac{1}{2}\stackrel{\leftrightarrow}{I}\right].$ (46)

With this expression “Equation (43)” can be written in the form

$F={\displaystyle {\oint}_{\sigma}\text{d}S\cdot {\mu}_{0}{\mu}_{r}{H}^{2}\left[\stackrel{^}{i}\stackrel{^}{i}-\frac{1}{2}\stackrel{\leftrightarrow}{I}\right]},$ (47)

where surface $\sigma $ can be taken as a parallelepiped with the surfaces parallel to the interface defined by vectors $\pm \stackrel{^}{k}$ , while the surfaces defined by vectors $\pm \stackrel{^}{i}$ and $\pm \stackrel{^}{j}$ do not contribute since they give contributions that cancel each other. Then the surface integral becomes

$F={\displaystyle \int \text{d}S\left(-\stackrel{^}{k}\right)\cdot {\mu}_{0}{\mu}_{r}{H}_{medium}^{2}\left[\stackrel{^}{i}\stackrel{^}{i}-\frac{1}{2}\stackrel{\leftrightarrow}{I}\right]}-{\displaystyle \int \text{d}S\stackrel{^}{k}\cdot {\mu}_{0}{H}_{vacuum}^{2}\left[\stackrel{^}{i}\stackrel{^}{i}-\frac{1}{2}\stackrel{\leftrightarrow}{I}\right]}\text{\hspace{0.17em}},$ (48)

and the force results

$F=\frac{1}{2}{\mu}_{0}{H}^{2}\left({\mu}_{r}-1\right){\displaystyle \int \stackrel{^}{k}\text{d}S}.$ (49)

Now, since

$\int \stackrel{^}{k}\text{d}S}=\stackrel{^}{k}{A}_{0},$ (50)

the force can be expressed as

$F=\frac{1}{2}{\mu}_{0}\left({\mu}_{r}-1\right){H}^{2}{A}_{0}\stackrel{^}{k},$ (51)

which with the constitutive relation “Equation (15)” takes the form

$F=\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}{A}_{0}\stackrel{^}{k},$ (52)

which is the expected result, “Equation (35)”,momentum balance equation give by “Equation (42)”.

8. Force Densities for Static Fields and Their Equivalence

Force densities like “Equations (5), (8) and (9)” are unfamiliar, but they are firmly founded on Maxwell’s equations, as we proceed to show.

We begin with the momentum balance equation derived from Maxwell’s equations, “Equation (36)”. This equation is equivalent to Maxwell’s equations with linear media, from which it is obtained by means of vector and dyadic identities. For static conditions and in absence of free charge and current densities this balance equation can be separated into two independent equations,

$\nabla \cdot \left\{DE-\frac{1}{2}I\left(D\cdot E\right)\right\}=-\frac{1}{2}\left[\left(\nabla E\right)\cdot D-\left(\nabla D\right)\cdot E\right]={f}_{Maxelecmedia},$ (53)

$\nabla \cdot \left\{BH-\frac{1}{2}I\left(B\cdot H\right)\right\}=-\frac{1}{2}\left[\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right]={f}_{Maxmagmedia},$ (54)

where subscripts in the right-hand members indicate the force densities for the electrostatic and magnetostatic cases.

We want now relate the balance “Equations (53) and (54)” with the electric and magnetic Helmholtz force densities, electric and magnetic,

${f}_{Helm\text{\hspace{0.17em}}elect}=-\frac{1}{2}{E}^{2}\nabla \u0454+\frac{1}{2}\nabla \left[{E}^{2}{\rho}_{m}\frac{\partial \u0454}{\partial {\rho}_{m}}\right]$ (55a)

${f}_{Helm\text{\hspace{0.17em}}mag}=-\frac{1}{2}{H}^{2}\nabla \mu +\frac{1}{2}\nabla \left[{H}^{2}{\rho}_{m}\frac{\partial \mu}{\partial {\rho}_{m}}\right].$ (55b)

It can be shown that these force densities are contained in “Equations (53) and (54)”, but for the present analysis it is enough to consider only the first terms on the right, since the dependence of the permittivity and permeability on mass density is for our purposes irrelevant. Then, using the constitutive relation

$D={\u0454}_{0}{\u0454}_{r}E$ (56)

the right-hand side of “Equation (53)” results

$\frac{1}{2}\left[\left(\nabla E\right)\cdot D-\left(\nabla D\right)\cdot E\right]=\frac{1}{2}{E}^{2}\nabla \u0454.$ (57)

Analogously, for the magnetic case, with the constitutive relation

$B={\mu}_{0}{\mu}_{r}H$ (58)

we obtain

$\frac{1}{2}\left[\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right]=\frac{1}{2}{H}^{2}\nabla \mu .$ (59)

The force densities “Equations (57) and (59)” look very different from the force densities “Equations (5), (8) and (9)”, but now it is evident their relation to Helmholtz’s force densities. Let us see how they are related.

The force density “Equation (5)”, which has been used to calculate the force on a dielectric slab inside a charged capacitor, is convenient to establish the equivalence among different force densities.

If the slab is in the z direction and has been introduced a distance l, the constitutive relation

$P={\u0454}_{0}{\chi}_{e}E,$ (60)

can be written as

$P\left(z\right)={\u0454}_{0}{\chi}_{e}\Theta \left(l-z\right)E\stackrel{^}{i},$ (61)

where $\Theta $ is Heaviside’s distribution, while $\stackrel{^}{i}$ is a unit vector in the x-direction, which is the direction of the electrostatic field. Then the force density “Equation (5)” becomes

$f=-\frac{1}{2}E\left({\partial}_{z}{P}_{x}\left(z\right)\right)\stackrel{^}{k},$ (62)

where $\stackrel{^}{k}$ is a unit vector in the z direction. With “Equation (60)” written as $P\left(z\right)={\u0454}_{0}{\chi}_{e}\left(z\right){E}_{x}$ we can express “Equation (62)” in the form

$f=-{\u0454}_{0}\frac{1}{2}{E}^{2}\left({\partial}_{z}{\chi}_{e}\left(z\right)\right)\stackrel{^}{k}.$ (63)

If we consider the relation between the relative permittivity and the susceptibility,

${\u0454}_{r}=1+{\chi}_{e},$ (64)

we can see that the right-hand member of “Equation (57)” equals “Equation (62)”. This validates firmly “Equation (5)”.

In the case of a magnetizable bar inside a solenoid, we can proceed analogously to establish that “Equation (8)” is equivalent to the right-hand member of “Equation (59)”.

If the axis of the solenoid is in the z direction, then “Equation (9)” turns into.

${f}_{newmag}=-\frac{1}{2}{\mu}_{0}H\left({\partial}_{z}M\left(z\right)\right).$ (65)

Since the constitutive relation is

$M={\chi}_{m}\Theta \left(l-z\right)H,$ (66)

“Equation (65)” becomes

${f}_{newmag}=-\frac{1}{2}{\mu}_{0}{H}^{2}\left({\partial}_{z}{\chi}_{m}\left(z\right)\right)\stackrel{^}{k}.$ (67)

Now, the relative permeability can be expressed as

${\mu}_{r}=1+{\chi}_{m},$ (68)

so that “Equation (67)” and the right-hand member of “Equation (59)” are equal, what establishes the validity of “Equation (8)”, which seems unfamiliar, but is equivalent to “Equation (59)”, a part of the Helmholtz force density.

When we calculate force densities from the gradient of energy densities, the vector identity “Equation (22)” leads in a natural way to force densities like “Equations (5), (8) and (9)”. What we have shown is that these force densities are contained in a momentum balance equation derived from Maxwell’s equations. The Faraday-Maxwell conception of a tension along the lines of force, and compressions around these lines, permits to interpret these forces as follows.

In the case of a dielectric slab inside a capacitor, the force arises from the difference of compressions at the interface. In the case of a magnetizable bar inside a solenoid, it is the difference in tensions at the interface, while in the device here proposed as a proper analogy of the slab inside a capacitor it is again the difference in compressions what produces the force.

9. Relation with Force Densities Derived as Gradients of Energy Densities

In order to establish the equivalence of force densities obtained as gradients of energy densities and the force densities obtained from momentum balance equations we use the dyadic identity

$\nabla \left(a\cdot b\right)=\left(\nabla a\right)\cdot b+\left(\nabla b\right)\cdot a.$ (69)

Then the force density “Equation (53)” can be written as

${f}_{Maxelecmedia}=\left(\nabla E\right)\cdot D-\frac{1}{2}\nabla \left(E\cdot D\right).$ (70)

Therefore, we have

$\frac{1}{2}\left[\left(\nabla E\right)\cdot D-\left(\nabla D\right)\cdot E\right]=\left(\nabla E\right)\cdot D-\frac{1}{2}\nabla \left(E\cdot D\right).$ (71)

Since the force on a dielectric slab inside a capacitor is given in terms of the polarization $P$ , we can use “Equations (56) and (64)”, obtaining

$\frac{1}{2}\left[\left(\nabla E\right)\cdot P-\left(\nabla P\right)\cdot E\right]=\left(\nabla E\right)\cdot P-\frac{1}{2}\nabla \left(E\cdot P\right).$ (72)

We have, on the other hand, the dyadic identity

$u\times \left(\nabla \times v\right)=\left(\nabla v\right)\cdot u-\left(u\cdot \nabla \right)v.$ (73)

If we set $v=E$ and $u=P$ , this identity becomes

$P\times \left(\nabla \times E\right)=\left(\nabla E\right)\cdot P-\left(P\cdot \nabla \right)E;$ (74)

but $\left(\nabla \times E\right)=0$ , so that

$\left(\nabla E\right)\cdot P=\left(P\cdot \nabla \right)E.$ (75)

Then, since the polarization $P\left(z\right)$ is in the x direction and the electric field is uniform,

$\left(P\cdot \nabla \right)E={P}_{x}\frac{\partial E}{\partial x}=0,$ (76)

concluding that in the case of the capacitor we have

$\left(\nabla E\right)\cdot P=0,$ (77)

so that “Equation (72)” reduces to

$\frac{1}{2}\left[\left(\nabla E\right)\cdot P-\left(\nabla P\right)\cdot E\right]=-\frac{1}{2}\nabla \left(E\cdot P\right).$ (78)

We have then that for the electric case the force density expressed in the right-hand member of “Equation (53)”, the Helmholtz force density, the force density given as the gradient of an energy density, “Equation (78)”, and the unfamiliar force density “Equation (5)” are all equivalent.

In the case of the solenoid, we proceed analogously applying identity “Equation (69)” to the right-hand member of “Equation (54)”, obtaining

$\frac{1}{2}\left[\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right]=\left(\nabla H\right)\cdot B-\frac{1}{2}\nabla \left(H\cdot B\right).$ (79)

In order to put this result in terms of the magnetization $M$ it is convenient to use the constitutive relation “Equation (39)” in “Equation (79)”, obtaining

$\frac{1}{2}{\mu}_{0}\left[\left(\nabla H\right)\cdot M-\left(\nabla M\right)\cdot H\right]={\mu}_{0}\left(\nabla H\right)\cdot B-\frac{1}{2}{\mu}_{0}\nabla \left(H\cdot B\right).$ (80)

However, the force expressed in “Equation (18)” is obtained assuming, for paramagnetic and diamagnetic media, that the field $H$ is nearly constant [12] so that

${\mu}_{0}\left(\nabla H\right)\cdot M=0.$ (81)

The results here obtained give support to this procedure, showing that an elementary problem may have interesting aspects.

Then “Equation (80)” is transformed into

${f}_{Maxmagmedia}=\frac{1}{2}{\mu}_{0}\left[\left(\nabla H\right)\cdot M-\left(\nabla M\right)\cdot H\right]=\frac{1}{2}{\mu}_{0}\nabla \left(H\cdot M\right),$ (82)

which is valid for the magnetic tension.

10. Summary of Results

In order to appreciate better the results obtained with our approach, and facilitate the comparison of electrostatic and magnetosatic effects, we present here a summary. These results are implications of Maxwell’s equations and the assumed constitutive relations applied to the electrostatic and magnetostaic conditions.

10.1. Magnetic Force Densities

$\nabla \cdot \left\{BH-\frac{1}{2}I\left(B\cdot H\right)\right\}=-\frac{1}{2}\left[\left(\nabla H\right)\cdot B-\left(\nabla B\right)\cdot H\right]={f}_{Maxmagmedia},$ (Am)

If the permeability does not depend of mass density, then

${f}_{Hmag}=-\frac{1}{2}{H}^{2}\nabla \mu $ (Bm)

Solenoid

$H,B,M$ in direction z, $M$ jumps in interface (z)

$\nabla \times H=0$ (C1)

$\nabla \times M=0$ (D1)

$\left(H\cdot \nabla \right)M\ne 0$ (E1)

$\left(M\cdot \nabla \right)H=0$ (F1)

${f}_{newM}=-\frac{1}{2}{\mu}_{0}\left(H\cdot \nabla \right)M,$ (G1)

$F=\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}{A}_{0}\stackrel{^}{k}.$ (H1)

${f}_{Helm\text{\hspace{0.17em}}mag}=-\frac{1}{2}{H}^{2}\nabla \mu =-\frac{1}{2}{H}^{2}\nabla {\mu}_{0}{\chi}_{m}=-\frac{1}{2}\left(H\cdot \nabla \right)M$ (K1)

The force arises from the tension part of the stress tensor.

10.2. Horseshoe Permanent Magnet

$H,B,M$ in direction x, $M$ jumps in interface (x)

$\nabla \times H=0$ (C2)

$\nabla \times M\ne 0$ (D2)

$\left(H\cdot \nabla \right)M=0$ (E2)

$\left(M\cdot \nabla \right)H=0$ (F2)

${f}_{newmag}=-\frac{1}{2}{\mu}_{0}H\times \left(\nabla \times M\right),$ (G2)

$F=\frac{1}{2}{\mu}_{0}{\chi}_{m}{H}^{2}{A}_{0}\stackrel{^}{k},$ (H2)

The force arises from the compression part of the stress tensor, equivalent to a difference of pressures at the interface.

${f}_{Helm\text{\hspace{0.17em}}mag}=-\frac{1}{2}{H}^{2}\nabla \mu =-\frac{1}{2}{H}^{2}\nabla {\mu}_{0}{\chi}_{m}=-\frac{1}{2}{\mu}_{0}H\times \left(\nabla \times M\right),$ (K2)

Both previous cases

$\nabla \times H=0$

$\left(M\cdot \nabla \right)H=0$

10.2.1. Electric Force Density

$\nabla \cdot \left\{DE-\frac{1}{2}I\left(D\cdot E\right)\right\}=-\frac{1}{2}\left[\left(\nabla E\right)\cdot D-\left(\nabla D\right)\cdot E\right]={f}_{Maxelecmedia},$ (Ae)

If the permittivity does not depend of mass density,

${f}_{Helm\text{\hspace{0.17em}}elect}=-\frac{1}{2}{E}^{2}\nabla \u0454$ (Be)

10.2.2. Capacitor

$E,D,P$ in direction z, $P$ jumps in interface (z)

$\nabla \times E=0$ (C3)

$\nabla \times P\ne 0$ (D3)

$\left(E\cdot \nabla \right)P=0$ (E3)

$\left(P\cdot \nabla \right)E=0$ (F3)

$f=-\frac{1}{2}E\times \left(\nabla \times P\right),$ (G3)

$F=\frac{1}{2}{\u0454}_{0}{\chi}_{e}{E}^{2}{A}_{0}\stackrel{^}{k},$ (H3)

${f}_{Helm\text{\hspace{0.17em}}elect}=-\frac{1}{2}{E}^{2}\nabla \u0454=-\frac{1}{2}{E}^{2}\nabla {\u0454}_{0}{\chi}_{e}=-\frac{1}{2}E\times \left(\nabla \times P\right)$ (K3)

The force arises from the compression part of the stress tensor, equivalent to a difference of pressures at the interface.

It is important to have in mind that

Solenoid: $\left(H\cdot \nabla \right)M\ne 0.$

Horseshoe: $\nabla \times M\ne 0.$

Capacitor: $\left(E\cdot \nabla \right)P\ne 0.$

11. Conclusions

We have shown elsewhere [16] [18] that the macroscopic Maxwell equations can be transformed, by means of vector and dyadic identities, into electromagnetic momentum balance equations. These balance equations involve different momentum flux tensors and force densities.

In the present work we use a particular balance equation to show the similarities and differences between the action of electrostatic and magnetostatic fields on linearly polarizable and magnetizable matter.

As a first point we can conclude that, based on the insights provided by this balance equation, the force exerted on a dielectric slab partially introduced into a charged capacitor is not exactly analogous to the force exerted on a magnetizable rod partially introduced into a solenoid. In the first case the force arises from the orthogonal compressions around the lines of force, equivalent to a difference of pressures at the interface, that is, the force is orthogonal to the electric field $E$ , while in the second case, the force arises from the action of the tension along the lines of force, that is, the force is parallel to the magnetic field $B$ .

As a second point, we have proposed a device that represents an exact analogy with the case of a dielectric slab inside a capacitor, in the sense that the force also arises from the compression around the magnetic lines of force, implying a difference of magnetic pressures at the interface. In this case the force is orthogonal to the magnetic field $B$ . In both cases the fringing fields are irrelevant, though in the electric case it is argued that the fringing field is the cause of the force, while in the magnetic case it is explicitly neglected [12] .

Acknowledgements

J. L. J. thanks the hospitality of Facultad de Ciencias, UNAM, during a sabbatical leave from Departamento de Física, UAM-I. We also acknowledge the valuable suggestions of a referee.

NOTES

*On sabatical leave from Departamento de Física, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana, Unidad Iztapalapa, Av. San Rafael Atlixco No. 186 Col. Vicentina.

References

[1] Campos-Flores, I., Jiménez-Ramírez, J.L. and Roa-Neri, J.A.E. (2018) How Does Arise the Force on a Dielestric Slab a Parallel Plate Capacitor? Journal of Electromagnetic Analysis and Applications, 10, 131-142.

https://doi.org/10.4236/jemaa.2018.107010

[2] Campos-Flores, I., Jiménez-Ramírez, J.L. and Roa-Neri J.A.E. (2018) A New Perspective of the Force on a Magnetizable Rod inside a Solenoid. Journal of Electromagnetic Analysis and Applications, 10, 171-183.

https://doi.org/10.4236/jemaa.2018.1010013

[3] Maxwell, J.C. (2005) Elementary Treatise on Electricity. 2nd Edition, Mineola New York, Dover.

[4] Becker, R. (1982) Electromagnetic Field and Interactions. Mineola.

[5] Stratton, J.A. (1941) Electromagnetic Theory. McGraw-Hill, New York.

[6] Panofsky, W.K.H. and Phillips, M. (1962) Classical Electricity and Magnetism. 2nd Edition, Addison-Wesley, Reading.

[7] Margulies, S. (1984) Force on a Dielectric Slab Inserted into a Parallel-Plate Capacitor. American Journal of Physics, 52, 515-518.

https://doi.org/10.1119/1.13861

[8] Utreras-Diaz, C. (1988) Dielectric Slab in a Parallel-Plate Condenser. American Journal of Physics, 56, 700-701.

https://doi.org/10.1119/1.15504

[9] Naini, A. and Green, M. (1997) Fringing Fields in a Parallel-Plate Capacitor. American Journal of Physics, 45, 877-879.

https://doi.org/10.1119/1.11075

[10] Dietz, E.R. (2004) Force on a Dielectric Slab: Fringing Field Approach. American Journal of Physics, 72, 1499-1500.

https://doi.org/10.1119/1.1764563

[11] Griffiths, D.J. (1999) Introduction to Electrodynamics. 3nd Edition, Prentice-Hall.

[12] Wangsness, R.K. (1986) Electromagnetic Theory. 2nd Edition, John Wiley and Sons, New York.

[13] Reitz, J.R., Milford, F.J and W. Christy, R.W. (1993) Foundations of Electromagnetic Theory. 4th Edition, Addison-Wesley, Reading.

[14] Purcell, E.M. and Morin, D.J. (2013) Electricity and Magnetism. 3rd Edition, Cambridge, New York.

[15] Jackson, J.D. (1998) Classical Electrodynamics. 3rd Edition, Wiley and Sons, New York.

[16] Campos, I., Jiménez, J.L. and López-Mariño, M.A. (2012) Electromagnetic Momentum Balance Equation and the Force Density in Material Media. Revista Brasileira de Ensino de Fisica, 34, 2303.

[17] Landau, L.D., Lifshitz, E.M. and Pitaevskii, L.P. (1984) Electrodynamics of Continuous Media. Vol. 8 Course if Theoretical Physics, 2nd Edition, Pergamon, Oxford.

[18] Jiménez, J.L., Campos, I. and, López-Mariño, M.A. (2013) Maxwell Equations in Material Media, Momentum Balance Equations and Force Densities Associated with Them. European Physical Journal Plus, 128, 46.

https://doi.org/10.1140/epjp/i2013-13046-8