Prev Subsubsection 3.12.2.3: Mixture Normal Traction Up Subsection 3.12.2: Surface Loads Subsubsection 3.12.2.5: Solute Flux Next
3.12.2.4 Fluid Flux
In a biphasic mixture of intrinsically incompressible solid and fluid constituents, the formulation adopted in FEBio implies that the normal component of the relative fluid flux is a natural boundary condition. If this boundary condition is not explicitly prescribed, the code automatically assumes that it is equal to zero. Therefore, biphasic boundaries are impermeable by default. (To implement a free-draining boundary, the fluid pressure nodal degrees of freedom should be set to zero.)
The flux of fluid relative to the solid matrix is given by the vector . Since viscosity is not explicitly modeled in a biphasic material, the tangential component of on a boundary surface may not be prescribed. Only the normal component of the relative fluid flux, , represents a natural boundary condition. To prescribe a value for on a surface, use:
<surface_load type="fluidflux" surface="surf1"> <flux lc="1">1.0</flux> <linear>0</linear> <mixture>1</mixture> </fluidflux>
The flux parameter defines the flux that will be applied to the surface. The optional parameter lc defines a loadcurve for the normal flux evolution. If omitted a constant fluid flux is applied.
When linear is set to zero (default) it means that the flux matches the prescribed value even if the surface on which it is applied changes in area as it deforms. Therefore, the net volumetric flow rate across the surface changes with changes in area. This type of boundary condition is useful if the fluid flux is known in the current configuration.
When linear is set to non-zero it means that the prescribed flux produces a volumetric flow rate based on the undeformed surface area in the reference configuration. Therefore, the flux in the current configuration does not match the prescribed value. This type of boundary condition is useful if the net volumetric flow rate across the surface is known. For example: Let be the known volumetric flow rate, let be the surface area in the reference configuration (a constant). Using “linear” means that the user prescribes / as the flux boundary condition. (However, regardless of the type, the fluid flux saved in the output file has a normal component equal to / , where area in current configuration.)
Prescribing on a free surface works only if the nodal displacements of the corresponding faces are also prescribed. If the nodal displacements are not known a priori, the proper boundary condition calls for prescribing the normal component of the mixture velocity, . To prescribe the value of on a surface, use
<surface_load type="fluidflux" surface="surf1"> <flux lc="1">1.0</flux> <linear>0</linear> <mixture>1</mixture> </surface_load>
For example, this boundary condition may be used when modeling a permeation problem through a biphasic material, when the upstream fluid velocity is prescribed, . If the upstream face is free, the companion boundary condition would be to let on that face as well.