Prev Section 2.7: Triphasic and Multiphasic Materials Up Section 2.7: Triphasic and Multiphasic Materials Section 2.8: Constrained Reactive Mixture of Solids Next
2.7.1 Governing Equations
In multiphasic materials the solvent is assumed to be neutral, whereas the solid and solutes may carry charge. The mixture is isothermal and all constituents are considered to be intrinsically incompressible. Since the viscosity of the fluid constituents (solvent and solutes) is considered negligible relative to the frictional interactions among constituents, the stress tensor for the mixture includes only a contribution from the fluid pressure and the stress in the solid, The mechano-chemical potential of the solvent is given by where is the solvent chemical potential in the solvent standard state, is the absolute temperature, is the true density of the solvent (which is invariant since the solvent is assumed intrinsically incompressible), is the fluid pressure, is the corresponding pressure in the standard state, is the universal gas constant, is the non-dimensional osmotic coefficient, and is the solution volume-based concentration of solute . The summation is taken over all solutes in the mixture. The mechano-electrochemical potential of each solute is similarly given by where is the molar mass of the solute, is its activity coefficient, is its solubility, is its charge number, and is its concentration in the solute standard state; is Faraday's constant, is the electrical potential of the mixture, and is the corresponding potential in the standard state.
In these relations, and are functions of state that describe the deviation of the mixture from ideal physico-chemical behavior; represents the fraction of the pore volume which may be occupied by solute . The standard state represents an arbitrary set of reference conditions for the physico-chemical state of each constituent. Therefore, the values of , , , , and , remain invariant over the entire domain of definition of an analysis. Since and appear together as a ratio, they may be combined into a single material function, , called the effective solubility.
In multiphasic mixtures, it is also assumed that electroneutrality is satisfied at every point in the continuum. Therefore, the net electrical charge summed over all constituents must reduce to zero, and no net charge accumulation may occur at any time. Denoting the fixed charge density of the solid by (moles of equivalent charge per solution volume), and recognizing that the solvent is always considered neutral, the electroneutrality condition may be written as This condition represents a constraint on a mixture of charged constituents. If none of the constituents are charged ( and for all , the constraint disappears.
Each constituent of the mixture must satisfy the axiom of mass balance. In the absence of chemical reactions involving constituent , its mass balance equation is where is the apparent density and is the velocity of that constituent. For solutes, the apparent density is related to the concentration according to , where is the volume fraction of the solid. When the solute volume fractions are negligible, it follows that , where is the solvent volume fraction. The molar flux of the solute relative to the solid is given by , where is the solute velocity. Using these relations, the mass balance relation for the solute may be rewritten as where represents the material time derivative in the spatial frame, following the solid; , where is the deformation gradient of the solid. This form of the mass balance for the solute is convenient for a finite element formulation where the mesh is defined on the solid matrix.
The volume flux of solvent relative to the solid is given by , where is the solvent velocity. When solute volume fractions are negligible, the mass balance equation for the mixture reduces to
Finally, the mass balance for the solid may be reduced to the form , which may be integrated to produce the algebraic relation , where is the solid volume fraction in the stress-free reference state of the solid.
Differentiating the electroneutrality condition in (2.7.1-4) using the material time derivative following the solid, and substituting the mass balance relations into the resulting expressions, produces a constraint on the solute fluxes: Recognizing that is the current density in the mixture, with representing Faraday's constant, the relation of (2.7.1-8) reduces to one of the Maxwell's equation, , in the special case when there can be no charge accumulation (electroneutrality).
As described in Section 2.6.2↑, the fluid pressure and solute concentrations are not continuous across boundaries of a mixture, whereas and 's for the solutes do satisfy continuity. Therefore, in a finite element implementation, the following continuous variables are used as nodal degrees of freedom: which represents the effective fluid pressure, and which represents the effective solute concentration. In the last expression, is the partition coefficient of the solute, which is related to the effective solubility and electric potential according to Physically, since is the osmotic (chemical) contribution to the fluid pressure, may be interpreted as that part of the total (mechano-chemical) fluid pressure which does not result from osmotic effects; thus, it is the mechanical contribution to . Similarly, the effective solute concentration represents the true contribution of the molar solute content to its electrochemical potential.
When using these variables instead of mechano-electrochemical potentials, the momentum equations for the solvent and solutes may be inverted to produce the following flux relations: and where is the effective hydraulic permeability of the solution (solvent solutes) in the mixture. The momentum equation for the mixture is