# Maxwell stress tensor

Last updated

The Maxwell stress tensor (named after James Clerk Maxwell) is a symmetric second-order tensor used in classical electromagnetism to represent the interaction between electromagnetic forces and mechanical momentum. In simple situations, such as a point charge moving freely in a homogeneous magnetic field, it is easy to calculate the forces on the charge from the Lorentz force law. When the situation becomes more complicated, this ordinary procedure can become impossibly difficult, with equations spanning multiple lines. It is therefore convenient to collect many of these terms in the Maxwell stress tensor, and to use tensor arithmetic to find the answer to the problem at hand.

## Contents

In the relativistic formulation of electromagnetism, the Maxwell's tensor appears as a part of the electromagnetic stress–energy tensor which is the electromagnetic component of the total stress–energy tensor. The latter describes the density and flux of energy and momentum in spacetime.

## Motivation

As outlined below, the electromagnetic force is written in terms of E and B. Using vector calculus and Maxwell's equations, symmetry is sought for in the terms containing E and B, and introducing the Maxwell stress tensor simplifies the result.

Maxwell's equations in SI units in vacuum
(for reference)
NameDifferential form
Gauss's law (in vacuum)${\displaystyle \nabla \cdot \mathbf {E} ={\frac {\rho }{\epsilon _{0}}}}$
Gauss's law for magnetism${\displaystyle \nabla \cdot \mathbf {B} =0}$
${\displaystyle \nabla \times \mathbf {E} =-{\frac {\partial \mathbf {B} }{\partial t}}}$
Ampère's circuital law (in vacuum)
(with Maxwell's correction)
${\displaystyle \nabla \times \mathbf {B} =\mu _{0}\mathbf {J} +\mu _{0}\epsilon _{0}{\frac {\partial \mathbf {E} }{\partial t}}\ }$
1. Starting with the Lorentz force law
${\displaystyle \mathbf {F} =q(\mathbf {E} +\mathbf {v} \times \mathbf {B} )}$

the force per unit volume is

${\displaystyle \mathbf {f} =\rho \mathbf {E} +\mathbf {J} \times \mathbf {B} }$
2. Next, ρ and J can be replaced by the fields E and B, using Gauss's law and Ampère's circuital law:
${\displaystyle \mathbf {f} =\epsilon _{0}\left({\boldsymbol {\nabla }}\cdot \mathbf {E} \right)\mathbf {E} +{\frac {1}{\mu _{0}}}\left({\boldsymbol {\nabla }}\times \mathbf {B} \right)\times \mathbf {B} -\epsilon _{0}{\frac {\partial \mathbf {E} }{\partial t}}\times \mathbf {B} \,}$
3. The time derivative can be rewritten to something that can be interpreted physically, namely the Poynting vector. Using the product rule and Faraday's law of induction gives
${\displaystyle {\frac {\partial }{\partial t}}(\mathbf {E} \times \mathbf {B} )={\frac {\partial \mathbf {E} }{\partial t}}\times \mathbf {B} +\mathbf {E} \times {\frac {\partial \mathbf {B} }{\partial t}}={\frac {\partial \mathbf {E} }{\partial t}}\times \mathbf {B} -\mathbf {E} \times ({\boldsymbol {\nabla }}\times \mathbf {E} )\,}$

and we can now rewrite f as

${\displaystyle \mathbf {f} =\epsilon _{0}\left({\boldsymbol {\nabla }}\cdot \mathbf {E} \right)\mathbf {E} +{\frac {1}{\mu _{0}}}\left({\boldsymbol {\nabla }}\times \mathbf {B} \right)\times \mathbf {B} -\epsilon _{0}{\frac {\partial }{\partial t}}\left(\mathbf {E} \times \mathbf {B} \right)-\epsilon _{0}\mathbf {E} \times ({\boldsymbol {\nabla }}\times \mathbf {E} ),}$

then collecting terms with E and B gives

${\displaystyle \mathbf {f} =\epsilon _{0}\left[({\boldsymbol {\nabla }}\cdot \mathbf {E} )\mathbf {E} -\mathbf {E} \times ({\boldsymbol {\nabla }}\times \mathbf {E} )\right]+{\frac {1}{\mu _{0}}}\left[-\mathbf {B} \times \left({\boldsymbol {\nabla }}\times \mathbf {B} \right)\right]-\epsilon _{0}{\frac {\partial }{\partial t}}\left(\mathbf {E} \times \mathbf {B} \right).}$
4. A term seems to be "missing" from the symmetry in E and B, which can be achieved by inserting (∇ ⋅ B)B because of Gauss' law for magnetism:
${\displaystyle \mathbf {f} =\epsilon _{0}\left[({\boldsymbol {\nabla }}\cdot \mathbf {E} )\mathbf {E} -\mathbf {E} \times ({\boldsymbol {\nabla }}\times \mathbf {E} )\right]+{\frac {1}{\mu _{0}}}\left[({\boldsymbol {\nabla }}\cdot \mathbf {B} )\mathbf {B} -\mathbf {B} \times \left({\boldsymbol {\nabla }}\times \mathbf {B} \right)\right]-\epsilon _{0}{\frac {\partial }{\partial t}}\left(\mathbf {E} \times \mathbf {B} \right).}$

Eliminating the curls (which are fairly complicated to calculate), using the vector calculus identity

${\displaystyle {\frac {1}{2}}{\boldsymbol {\nabla }}(\mathbf {A} \cdot \mathbf {A} )=\mathbf {A} \times ({\boldsymbol {\nabla }}\times \mathbf {A} )+(\mathbf {A} \cdot {\boldsymbol {\nabla }})\mathbf {A} ,}$

${\displaystyle \mathbf {f} =\epsilon _{0}\left[({\boldsymbol {\nabla }}\cdot \mathbf {E} )\mathbf {E} +(\mathbf {E} \cdot {\boldsymbol {\nabla }})\mathbf {E} \right]+{\frac {1}{\mu _{0}}}\left[({\boldsymbol {\nabla }}\cdot \mathbf {B} )\mathbf {B} +(\mathbf {B} \cdot {\boldsymbol {\nabla }})\mathbf {B} \right]-{\frac {1}{2}}{\boldsymbol {\nabla }}\left(\epsilon _{0}E^{2}+{\frac {1}{\mu _{0}}}B^{2}\right)-\epsilon _{0}{\frac {\partial }{\partial t}}\left(\mathbf {E} \times \mathbf {B} \right).}$
5. This expression contains every aspect of electromagnetism and momentum and is relatively easy to compute. It can be written more compactly by introducing the Maxwell stress tensor,
${\displaystyle \sigma _{ij}\equiv \epsilon _{0}\left(E_{i}E_{j}-{\frac {1}{2}}\delta _{ij}E^{2}\right)+{\frac {1}{\mu _{0}}}\left(B_{i}B_{j}-{\frac {1}{2}}\delta _{ij}B^{2}\right).}$

All but the last term of f can be written as the tensor divergence of the Maxwell stress tensor, giving:

${\displaystyle \nabla \cdot {\boldsymbol {\sigma }}=\mathbf {f} +\epsilon _{0}\mu _{0}{\frac {\partial \mathbf {S} }{\partial t}}\,}$,

As in the Poynting's theorem, the second term on the right side of the above equation can be interpreted as the time derivative of the EM field's momentum density, while the first term is the time derivative of the momentum density for the massive particles. In this way, the above equation will be the law of conservation of momentum in classical electrodynamics.

where the Poynting vector has been introduced

${\displaystyle \mathbf {S} ={\frac {1}{\mu _{0}}}\mathbf {E} \times \mathbf {B} .}$

in the above relation for conservation of momentum, ${\displaystyle \nabla \cdot {\boldsymbol {\sigma }}}$ is the momentum flux density and plays a role similar to ${\displaystyle \mathbf {S} }$ in Poynting's theorem.

The above derivation assumes complete knowledge of both ρ and J (both free and bounded charges and currents). For the case of nonlinear materials (such as magnetic iron with a BH-curve), the nonlinear Maxwell stress tensor must be used. [1]

## Equation

In physics, the Maxwell stress tensor is the stress tensor of an electromagnetic field. As derived above in SI units, it is given by:

${\displaystyle \sigma _{ij}=\epsilon _{0}E_{i}E_{j}+{\frac {1}{\mu _{0}}}B_{i}B_{j}-{\frac {1}{2}}\left(\epsilon _{0}E^{2}+{\frac {1}{\mu _{0}}}B^{2}\right)\delta _{ij}}$,

where ε0 is the electric constant and μ0 is the magnetic constant, E is the electric field, B is the magnetic field and δij is Kronecker's delta. In Gaussian cgs unit, it is given by:

${\displaystyle \sigma _{ij}={\frac {1}{4\pi }}\left(E_{i}E_{j}+H_{i}H_{j}-{\frac {1}{2}}\left(E^{2}+H^{2}\right)\delta _{ij}\right)}$,

where H is the magnetizing field.

An alternative way of expressing this tensor is:

${\displaystyle {\overset {\leftrightarrow }{\boldsymbol {\sigma }}}={\frac {1}{4\pi }}\left[\mathbf {E} \otimes \mathbf {E} +\mathbf {H} \otimes \mathbf {H} -{\frac {E^{2}+H^{2}}{2}}\mathbb {I} \right]}$

where ⊗ is the dyadic product, and the last tensor is the unit dyad:

${\displaystyle \mathbb {I} \equiv {\begin{pmatrix}1&0&0\\0&1&0\\0&0&1\end{pmatrix}}=\left(\mathbf {\hat {x}} \otimes \mathbf {\hat {x}} +\mathbf {\hat {y}} \otimes \mathbf {\hat {y}} +\mathbf {\hat {z}} \otimes \mathbf {\hat {z}} \right)}$

The element ij of the Maxwell stress tensor has units of momentum per unit of area per unit time and gives the flux of momentum parallel to the ith axis crossing a surface normal to the jth axis (in the negative direction) per unit of time.

These units can also be seen as units of force per unit of area (negative pressure), and the ij element of the tensor can also be interpreted as the force parallel to the ith axis suffered by a surface normal to the jth axis per unit of area. Indeed, the diagonal elements give the tension (pulling) acting on a differential area element normal to the corresponding axis. Unlike forces due to the pressure of an ideal gas, an area element in the electromagnetic field also feels a force in a direction that is not normal to the element. This shear is given by the off-diagonal elements of the stress tensor.

## Magnetism only

If the field is only magnetic (which is largely true in motors, for instance), some of the terms drop out, and the equation in SI units becomes:

${\displaystyle \sigma _{ij}={\frac {1}{\mu _{0}}}B_{i}B_{j}-{\frac {1}{2\mu _{0}}}B^{2}\delta _{ij}\,.}$

For cylindrical objects, such as the rotor of a motor, this is further simplified to:

${\displaystyle \sigma _{rt}={\frac {1}{\mu _{0}}}B_{r}B_{t}-{\frac {1}{2\mu _{0}}}B^{2}\delta _{rt}\,.}$

where r is the shear in the radial (outward from the cylinder) direction, and t is the shear in the tangential (around the cylinder) direction. It is the tangential force which spins the motor. Br is the flux density in the radial direction, and Bt is the flux density in the tangential direction.

## In electrostatics

In electrostatics the effects of magnetism are not present. In this case the magnetic field vanishes, ${\displaystyle \mathbf {B} =\mathbf {0} }$, and we obtain the electrostatic Maxwell stress tensor. It is given in component form by

${\displaystyle \sigma _{ij}=\varepsilon _{0}E_{i}E_{j}-{\frac {1}{2}}\varepsilon _{0}E^{2}\delta _{ij}}$

and in symbolic form by

${\displaystyle {\boldsymbol {\sigma }}=\varepsilon _{0}\mathbf {E} \otimes \mathbf {E} -{\frac {1}{2}}\varepsilon _{0}(\mathbf {E} \cdot \mathbf {E} )\mathbf {I} }$

where ${\displaystyle \mathbf {I} }$ is the appropriate identity tensor (usually ${\displaystyle 3\times 3}$).

## Eigenvalue

The eigenvalues of the Maxwell stress tensor are given by:[ citation needed ]

${\displaystyle \{\lambda \}=\left\{-\left({\frac {\epsilon _{0}}{2}}E^{2}+{\frac {1}{2\mu _{0}}}B^{2}\right),~\pm {\sqrt {\left({\frac {\epsilon _{0}}{2}}E^{2}-{\frac {1}{2\mu _{0}}}B^{2}\right)^{2}+{\frac {\epsilon _{0}}{\mu _{0}}}\left({\boldsymbol {E}}\cdot {\boldsymbol {B}}\right)^{2}}}\right\}}$

These eigenvalues are obtained by iteratively applying the Matrix Determinant Lemma, in conjunction with the Sherman-Morrison Formula.

Noting that the characteristic equation matrix, ${\displaystyle {\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } }$, can be written as

${\displaystyle {\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } =-\left(\lambda +V\right)\mathbf {\mathbb {I} } +\epsilon _{0}\mathbf {E} \mathbf {E} ^{\textsf {T}}+{\frac {1}{\mu _{0}}}\mathbf {B} \mathbf {B} ^{\textsf {T}}}$

where

${\displaystyle V={\frac {1}{2}}\left(\epsilon _{0}E^{2}+{\frac {1}{\mu _{0}}}B^{2}\right)}$

we set

${\displaystyle \mathbf {U} =-\left(\lambda +V\right)\mathbf {\mathbb {I} } +\epsilon _{0}\mathbf {E} \mathbf {E} ^{\textsf {T}}}$

Applying the Matrix Determinant Lemma once, this gives us

${\displaystyle \det {\left({\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } \right)}=\left(1+{\frac {1}{\mu _{0}}}\mathbf {B} ^{\textsf {T}}\mathbf {U} ^{-1}\mathbf {B} \right)\det {\left(\mathbf {U} \right)}}$

Applying it again yields,

${\displaystyle \det {\left({\overleftrightarrow {\boldsymbol {\sigma }}}-\lambda \mathbf {\mathbb {I} } \right)}=\left(1+{\frac {1}{\mu _{0}}}\mathbf {B} ^{\textsf {T}}\mathbf {U} ^{-1}\mathbf {B} \right)\left(1-{\frac {\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} }{\lambda +V}}\right)\left(-\lambda -V\right)^{3}}$

From the last multiplicand on the RHS, we immediately see that ${\displaystyle \lambda =-V}$ is one of the eigenvalues.

To find the inverse of ${\displaystyle \mathbf {U} }$, we use the Sherman-Morrison formula:

${\displaystyle \mathbf {U} ^{-1}=-\left(\lambda +V\right)^{-1}-{\frac {\epsilon _{0}\mathbf {E} \mathbf {E} ^{\textsf {T}}}{\left(\lambda +V\right)^{2}-\left(\lambda +V\right)\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} }}}$

Factoring out a ${\displaystyle \left(-\lambda -V\right)}$ term in the determinant, we are left with finding the zeros of the rational function:

${\displaystyle \left(-\left(\lambda +V\right)-{\frac {\epsilon _{0}\left(\mathbf {E} \cdot \mathbf {B} \right)^{2}}{\mu _{0}\left(-\left(\lambda +V\right)+\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} \right)}}\right)\left(-\left(\lambda +V\right)+\epsilon _{0}\mathbf {E} ^{\textsf {T}}\mathbf {E} \right)}$

Thus, once we solve

${\displaystyle -\left(\lambda +V\right)\left(-\left(\lambda +V\right)+\epsilon _{0}E^{2}\right)-{\frac {\epsilon _{0}}{\mu _{0}}}\left(\mathbf {E} \cdot \mathbf {B} \right)^{2}=0}$

we obtain the other two eigenvalues.

## Related Research Articles

Acoustic theory is a scientific field that relates to the description of sound waves. It derives from fluid dynamics. See acoustics for the engineering approach.

In physics, the Navier–Stokes equations, named after Claude-Louis Navier and George Gabriel Stokes, describe the motion of viscous fluid substances.

In special relativity, a four-vector is an object with four components, which transform in a specific way under Lorentz transformation. Specifically, a four-vector is an element of a four-dimensional vector space considered as a representation space of the standard representation of the Lorentz group, the (½,½) representation. It differs from a Euclidean vector in how its magnitude is determined. The transformations that preserve this magnitude are the Lorentz transformations, which include spatial rotations and boosts.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.

Least mean squares (LMS) algorithms are a class of adaptive filter used to mimic a desired filter by finding the filter coefficients that relate to producing the least mean square of the error signal. It is a stochastic gradient descent method in that the filter is only adapted based on the error at the current time. It was invented in 1960 by Stanford University professor Bernard Widrow and his first Ph.D. student, Ted Hoff.

In continuum mechanics, the finite strain theory—also called large strain theory, or large deformation theory—deals with deformations in which strains and/or rotations are large enough to invalidate assumptions inherent in infinitesimal strain theory. In this case, the undeformed and deformed configurations of the continuum are significantly different, requiring a clear distinction between them. This is commonly the case with elastomers, plastically-deforming materials and other fluids and biological soft tissue.

The electromagnetic wave equation is a second-order partial differential equation that describes the propagation of electromagnetic waves through a medium or in a vacuum. It is a three-dimensional form of the wave equation. The homogeneous form of the equation, written in terms of either the electric field E or the magnetic field B, takes the form:

The upper-convected Maxwell (UCM) model is a generalisation of the Maxwell material for the case of large deformations using the upper-convected time derivative. The model was proposed by James G. Oldroyd. The concept is named after James Clerk Maxwell.

A hyperelastic or green elastic material is a type of constitutive model for ideally elastic material for which the stress–strain relationship derives from a strain energy density function. The hyperelastic material is a special case of a Cauchy elastic material.

The magnetic tension force is a restoring force that acts to straighten bent magnetic field lines. It equals:

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.

The Liénard–Wiechert potentials describe the classical electromagnetic effect of a moving electric point charge in terms of a vector potential and a scalar potential in the Lorenz gauge. Built directly from Maxwell's equations, these potentials describe the complete, relativistically correct, time-varying electromagnetic field for a point charge in arbitrary motion, but are not corrected for quantum-mechanical effects. Electromagnetic radiation in the form of waves can be obtained from these potentials. These expressions were developed in part by Alfred-Marie Liénard in 1898 and independently by Emil Wiechert in 1900.

The intent of this article is to highlight the important points of the derivation of the Navier–Stokes equations as well as its application and formulation for different families of fluids.

In continuum mechanics, a compatible deformation tensor field in a body is that unique tensor field that is obtained when the body is subjected to a continuous, single-valued, displacement field. Compatibility is the study of the conditions under which such a displacement field can be guaranteed. Compatibility conditions are particular cases of integrability conditions and were first derived for linear elasticity by Barré de Saint-Venant in 1864 and proved rigorously by Beltrami in 1886.

Plasticity theory for rocks is concerned with the response of rocks to loads beyond the elastic limit. Historically, conventional wisdom has it that rock is brittle and fails by fracture while plasticity is identified with ductile materials. In field scale rock masses, structural discontinuities exist in the rock indicating that failure has taken place. Since the rock has not fallen apart, contrary to expectation of brittle behavior, clearly elasticity theory is not the last work.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used for discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

In optics, the Ewald–Oseen extinction theorem, sometimes referred to as just "extinction theorem", is a theorem that underlies the common understanding of scattering. It is named after Paul Peter Ewald and Carl Wilhelm Oseen, who proved the theorem in crystalline and isotropic media, respectively, in 1916 and 1915. Originally, the theorem applied to scattering by an isotropic dielectric objects in free space. The scope of the theorem was greatly extended to encompass a wide variety of bianisotropic media.

In physics, zilch is a conserved quantity of the electromagnetic field.

## References

1. Brauer, John R. (2014-01-13). Magnetic Actuators and Sensors. John Wiley & Sons. ISBN   9781118754979.
• David J. Griffiths, "Introduction to Electrodynamics" pp. 351–352, Benjamin Cummings Inc., 2008
• John David Jackson, "Classical Electrodynamics, 3rd Ed.", John Wiley & Sons, Inc., 1999.
• Richard Becker, "Electromagnetic Fields and Interactions", Dover Publications Inc., 1964.