Electromagnetic Force

At a boundary or interface, Simcenter STAR-CCM+ calculates the electromagnetic force by evaluating the electromechanical stress tensor.

For a linear material, the electromechanical stress tensor is defined as:

Figure 1. EQUATION_DISPLAY
σEM(B)=μ1BB(12μ1BB)I
(4347)

and is also referred to as the Maxwell stress tensor.

For a nonlinear material, the electromechanical stress tensor is defined as:

Figure 2. EQUATION_DISPLAY
σEM(B)=HB(HB0|B|μ1(B)dB)I
(4348)

At an interface between two materials, the electromechanical stress tensor is generally discontinuous. Simcenter STAR-CCM+ calculates the electromagnetic force density at the interface as:

Figure 3. EQUATION_DISPLAY
fEM=(σEM1σEM0)n01
(4349)

where σEM1 and σEM0 represent the electromechanical stress tensor at each side of the interface, and n01 is the surface normal pointing from side 0 to side 1.

In general, a material is considered force-free when it is free from electric current sources, including user-defined sources, eddy currents, and currents induced by excitation coils or permanent magnets. The total electromagnetic force acting on a body that is completely surrounded by a force-free material can be calculated as:

Figure 4. EQUATION_DISPLAY
FEM=Apda
(4350)

where A is the surface enclosing the body and p is the electromagnetic stress vector:

Figure 5. EQUATION_DISPLAY
p=σn
(4351)

where n is a unit vector normal to A (pointing outwards) and σ is the Maxwell stress tensor at surface A.

From Eqn. (4350), the magnetic torque acting on the body can be calculated as:

Figure 6. EQUATION_DISPLAY
τEM=Ar×pda
(4352)

where r is a user-specified position vector.

The finite element implementation also allows for calculation of the electromagnetic force at mesh element nodes. Simcenter STAR-CCM+ computes the nodal forces as:

Figure 7. EQUATION_DISPLAY
FM=MVσEMNMdV
(4353)

where M identifies the element nodes, V is the set of elements containing the node M, and NM represents H1 Lagrangian shape functions [845].