The Electromagnetic Stress Tensor
Within a vacuum or other medium, forces between charged bodies can be computed on the assumption that a fictitious state of stress exists within the field. Historically the consistency of this approach led 19th-century physicists to postulate the existence of the ether, a ubiquitous medium through which electromagnetic forces propagate. While these ideas have been superseded by the development of special relativity, the use of an electromagnetic stress tensor (also known as the Maxwell stress tensor) remains an accurate and convenient technique to compute electromagnetic forces.
Maxwell stresses exists in all dielectric materials, including air and even free space. However, the force magnitude is rather small, and it usually can cause significant deformations only at small dimension.
The electromagnetic stress tensor in a vacuum (in the absence of magnetic fields) is given by (Ref. 1):
(5-2)
Where E is the electric field, I is the identity tensor, ε0 is the permittivity of free space, and
Within an isotropic linear dielectric under small deformations, the electric displacement vector is related to the electric field as (Ref. 2):
where the electric susceptibility χ can be a function of the mechanical strain in the material
where the small strain tensor is given by
The corresponding electromagnetic tress tensor takes the following form:
(5-3)
where the relative permittivity of the material with the absence of deformation is introduced as . It can be written equivalently as
where the two terms represent, respectively, the contributions from: the underlying free space and the material polarization.
The Electromechanical Forces coupling feature applies the material and underlying free space contributions to the stress as a body load. In weak form: .
Large deformations
For finite deformations, the expressions for the electromagnetic stress and material polarization can be derived using the following thermodynamic potential called electric enthalpy:
where the subscript m indicates that the vector components must be taken on the material frame, and the right Cauchy-Green deformation tensor is
with , and . The mechanical energy function depends on the solid model used.
The total second Piola-Kirchhoff stress tensor is given by
and the electric displacement can be calculated as
The electromagnetic stress tensor can be used to compute the forces acting on a dielectric body. From Equation 5-1 the balance of forces at the surface of a dielectric body (material 1) in vacuum or air (material 2) implies:
where the total stress tensors in media 1 and 2 are, respectively:
where σm is a mechanical component of the total stress and the ambient pressure p2 in to surrounding air or other medium (if present). Using Equation 5-2:
where the electric field E is computed in material 2.
The Electromechanical Forces coupling feature applies the above defined traction as a boundary load on the boundaries which are external for the coupling feature selection but internal for the corresponding Electrostatics interface.
COMSOL Multiphysics does not explicitly include the ambient pressure definition on the Electromechanical Forces coupling feature. However, it is possible to add an additional surface force to the corresponding Solid Mechanics interface if the pressure is known or computed by another physics interface.
About different stress forms
There exists a long-lasting controversy in scientific literature about the definition of electromagnetic forces acting on solids. An extensive review can be found in Ref. 3. Many classic textbooks (for example, Ref. 1 and Ref. 2) operate with the so-called Minkowski stress tensor that is usually written as:
The corresponding electromagnetic body force can be written as
which sometimes is referred to as the Korteweg-Helmholtz force. For homogeneous materials without deformation, one has ∇χ = 0. Hence, in the absence of free charges (ρe = 0), the body force becomes zero. Thus, the whole electromechanical load on the solid is due to the Maxwell stress jumps at the boundaries between domains with different material properties.
There exists two most often used alternatives. The first one is called Einstein-Laub stress tensor:
This form is widely accepted in modern electroelasticity and material science, see Ref. 4. The corresponding body force can be written as:
which is also called Kelvin force. Note that it is nonzero as soon as there are electric field variation and polarization within the material:
The other alternative is Chu stress tensor:
(5-4)
and the corresponding electromagnetic body force is given by
Note that in contrast to the above equation, the expression for the stress tensor itself (Equation 5-4) seems to ignore the material polarization properties. However, if the frame difference is important because of the deformations, one has
(5-5)
The electric enthalpy can be written as
Then,
and the electromechanical second Piola-Kirchhoff stress tensor is given by
(5-6)
where σEM is given by Equation 5-4, which shows the consistency of the theory. The key assumption behind the derivation of Equation 5-6 is that the material components of the electric susceptibility χm are independent of deformations.
COMSOL Multiphysics provides a choice of the electromagnetic stress model, which is available in the Electromechanical Forces coupling feature. Any of the above presented three options can be selected; the default choice is the Minkowski stress tensor.