Converted document
$\newcommand{\lyxlock}{}$
Subsection 3.13.2: Biphasic Contact Up Section 3.13: Contact Section Subsection 3.13.4: Rigid Wall Interfaces

### 3.13.3 Biphasic-Solute and Multiphasic Contact

The sliding-biphasic-solute implementation for sliding interfaces can deal with biphasic-solute contact surfaces (including biphasic-solute-on-biphasic-solute, biphasic-solute-on-biphasic, biphasic-solute-on-elastic, biphasic-solute-on-rigid) and the sliding-multiphasic contact interface can similarly deal with multiphasic contact surfaces. These contact interfaces allow for the possibility to track fluid and solute flow across the contact interface [12]. In other words, fluid and solute can flow from one side of the contact interface to the other. To use this feature, the user must define additional contact parameters, namely:
<pressure_penalty>1.0</pressure_penalty>
<concentration_penalty>1.0</concentration_penalty>
<ambient_pressure>0</ambient_pressure>
<ambient_concentration>0</ambient_concentration> (for sliding-biphasic-solute)

or
<ambient_concentration sol="id">0</ambient_concentration> (for sliding-multiphasic)

In the same way that the penalty parameter controls the contact tractions, these penalty parameters control the penalty values that are used to calculate the Lagrange multipliers for the pressure and concentration constraints. If the laugon flag is set, the augmented Lagrangian method is used to enforce the pressure and concentration constraints. And if the auto_penalty flag is defined, an initial guess for the pressure and concentration penalty is calculated automatically using the following formulas: where is the element's area, is the element's volume, is a measure of the fluid permeability which is defined as one third of the trace of the material's initial permeability tensor, and is a measure of the solute diffusivity which is defined as one third of the trace of the material's initial diffusivity tensor.
When either contact surface is biphasic-solute or multiphasic, the surface outside the contact area(s) is automatically set to ambient conditions (equivalent to setting the effective fluid pressure and effective solute concentration to the <ambient_pressure> and <ambient_concentration> values, respectively). Ambient conditions may also be associated with a load curve, for example:
<ambient_pressure lc="2">1.0</ambient_pressure>
<ambient_concentration lc="3">1.0</ambient_concentration>

When performing biphasic-solute-on-rigid or multiphasic-on-rigid contact, a two-pass analysis should not be used; the rigid surface should be the secondary surface.
Subsection 3.13.2: Biphasic Contact Up Section 3.13: Contact Section Subsection 3.13.4: Rigid Wall Interfaces