Assuming that a charge density Ο*c* exists everywhere in the domain, the net electrical potential energy is given by Ο*c* Ο*e* , where Ο*e* is the electric potential. The total electrical energy of the system at steady state is given by
Display Formula

$\Xi \xa8e=\beta \x88\xab\Xi \xa9\Omicron \x81c\Omicron \x86ed\Xi \xa9$

(4)

The rate of change of this quantity, assuming equilibrium is reached for electrical potential much faster than the rest of the system, is obtained by taking the material time derivative [

11] of Eq.

4Display Formula$\Xi \xa8\Beta \xb7e=\beta \x88\xab\Xi \xa9\beta \x88\x82(\Omicron \x81c\Omicron \x86e)\beta \x88\x82td\Xi \xa9+\beta \x88\xab\Xi \x93ext(\Omicron \x81c\Omicron \x86e)vnd\Xi \x93+\beta \x88\xab\Xi \x93int[[(\Omicron \x81c\Omicron \x86e)vn]]d\Xi \x93$

(5)

where [[]] is the jump in the quantities across the interface defined as [[

*a*]]β=β

*a*+ β+β

*a*β , referring to the quantity on either side of the interface. Assuming that the material velocity is significantly slower than the time taken to reach steady state, we can assume that the system is always at steady electrical state. Hence, the spontaneous change term,

$\beta \x88\x82(\Omicron \x81c\Omicron \x86e)\beta \x88\x82t\beta \x86\x920$. Also, if the external boundaries of the system are fixed,

$\beta \x88\xab\Xi \x93ext(\Omicron \x81c\Omicron \x86e)vnd\Xi \xa9=0$. Further, for a materialβvoid system, one side of the jump term can be ignored. Hence, the free energy rate due to electrical work can be approximated in electromigration-induced void evolution problems as

Display Formula$\Xi \xa8\Beta \xb7e=\beta \x88\xab\Xi \x93int(\Omicron \x81c\Omicron \x86e)vnd\Xi \x93$

(6)

Similarly, the contribution to free energy from mechanical loads is

Display Formula$\Xi \xa8m=\beta \x88\xab\Xi \xa9\Omicron \x86md\Xi \xa9$

(7)

where Ο

*m* is the strain energy density. To evaluate the mechanical energy contributions, from the balance of mechanical energy on the system, we get

Display Formula$\beta \x88\xab\Xi \x93t\beta \x86\x92\Beta \xb7vd\Xi \x93=\beta \x88\xab\Xi \xa9\Omicron \x83:Dd\Xi \xa9$

(8)

where

*D* is the rate of deformation tensor. Assuming linear elasticity (Οβ=β

*E*: Ι), and therefore the strain energy density

$\Omicron \x86m=(12)\Iota \x9b:E:\Iota \x9b$, we have, after similar considerations for equilibrium of the stresses

Display Formula$\beta \x88\xab\Xi \x93t\beta \x86\x92\Beta \xb7vd\Xi \x93=\beta \x88\xab\Xi \xa9\beta \x88\x82\Omicron \x86m\beta \x88\x82t\beta \x80\x83d\Xi \xa9+\beta \x88\xab\Xi \x93ext\Omicron \x86mvnd\Xi \x93+\beta \x88\xab\Xi \x93int[[\Omicron \x86mvn]]d\Xi \x93$

(9)

Arguing, as before, that the system reaches equilibrium instantly relative to the times involved in mass transport, we have

$\beta \x88\x82\Omicron \x86m\beta \x88\x82t\beta \x86\x920$. Assuming fixed external boundaries, and no loads in the void, we get

Display Formula$\Xi \xa8\Beta \xb7m=\beta \x88\xab\Xi \x93t\beta \x86\x92\Beta \xb7vd\Xi \x93=\beta \x88\xab\Xi \x93int\Omicron \x86mvnd\Xi \x93$

(10)

Combining Eqs.

3,

6,

10, we get

Display Formula$\Xi \xa8\Beta \xb7=\beta \x88\xab\Xi \x93int[\Xi \xb3\Xi \u038a+\Omicron \x81c\Omicron \x86e+\Omicron \x86m]vnd\Xi \x93$

(11)

Following the arguments from Ref. [

12], and introducing surface and bulk mobilities

*Ms * and

*Mb *, we get the following equations for the velocity of the front.