Theory Manual Version 3.4
$\newcommand{\lyxlock}{}$
Subsection 3.3.3: Discretization Up Chapter 3: The Nonlinear FE Method Subsection 3.4.1: Linearization along

## 3.4 Weak Formulation for Multiphasic Materials

The virtual work integral for a mixture of intrinsically incompressible constituents combines the balance of momentum for the mixture, the balance of mass for the mixture, and the balance of mass for each of the solutes. In addition, for charged mixtures, the condition of (2.7.1-8) may be enforced as a penalty constraint on each solute mass balance equation: where is the virtual velocity of the solid, is the virtual effective fluid pressure, and is the virtual molar energy of solute . Here, represents the mixture domain in the spatial frame and is an elemental volume in . Applying the divergence theorem, may be split into internal and external contributions to the virtual work, , where and In these expressions, , is the boundary of , and is an elemental area on . In this finite element formulation, , and are used as nodal variables, and essential boundary conditions may be prescribed on these variables. Natural boundary conditions are prescribed to the mixture traction, , normal fluid flux, , and normal solute flux, , where is the outward unit normal to . To solve the system for nodal values of , and , it is necessary to linearize these equations, as shown for example in Sections 3.3.1↑-3.3.2↑ for biphasic-solute materials. If the mixture is charged, it is also necessary to solve for the electric potential by solving the algebraic relation of the electroneutrality condition in (2.7.1-4), which may be rewritten as In the special case of a triphasic mixture, where solutes consist of two counter-ions (, this equation may be solved in closed form to produce Only the positive root is valid in the argument of the logarithm function.
Subsection 3.3.3: Discretization Up Chapter 3: The Nonlinear FE Method Subsection 3.4.1: Linearization along