21.3. Electric Potential Method

The second approach for the current density is to solve the electric potential equation and calculate the current density using Ohm’s law. In general, the electric field can be expressed as:

(21–19)

where and are the scalar potential and the vector potential, respectively. For a static field and assuming , Ohm’s law given in Equation 21–8 can be written as:

(21–20)

For sufficiently conducting media, the principle of conservation of electric charge gives:

(21–21)

The electric potential equation is therefore given by:

(21–22)

The boundary condition for the electric potential is given by:

(21–23)

for an insulating boundary, where is the unit vector normal to the boundary, and

(21–24)

for a conducting boundary, where is the specified potential at the boundary. The current density can then be calculated from Equation 21–20.

With the knowledge of the induced electric current, the MHD coupling is achieved by introducing additional source terms to the fluid momentum equation and energy equation. For the fluid momentum equation, the additional source term is the Lorentz force given by:

(21–25)

which has units of in the SI system. For the energy equation, the additional source term is the Joule heating rate given by:

(21–26)

which has units of .

For charged particles in an electromagnetic field, the Lorentz force acting on the particle is given by:

(21–27)

where is the particle charge density () and is the particle velocity. The force has units of .

For multiphase flows, assuming that the electric surface current at the interface between phases can be ignored, the electric conductivity for the mixture is given by:

(21–28)

where and are respectively the electric conductivity and volume fraction of phase . is used in solving the induction equations.