Prev Section 2.13: Hybrid Biphasic Material Up Section 2.13: Hybrid Biphasic Material Subsection 2.13.2: BFSI Continuous Variables Next
2.13.1 BFSI Governing Equations
We model the BFSI domain as an unconstrained mixture of an isothermal, compressible, and viscous fluid and an isothermal, deformable porous solid whose skeleton material is intrinsically incompressible. As usual in mixture theory, each constituent ( for solid and fluid, respectively) has its own apparent density , which represents the ratio of the elemental mass of constituent in the elemental mixture volume . This apparent density is related to the true density (mass of per volume of ) via , where is the volume fraction of in the mixture, satisfying . Since the solid skeleton is intrinsically incompressible, is constant. The boundaries of a biphasic domain are defined on the porous solid matrix. The deformation gradient of the solid is denoted by and its determinant, , represents the ratio of the mixture elemental volumes in the current ( ) and reference ( ) configurations. Thus, since the solid skeleton is intrinsically incompressible, purely represents the compressibility of the pore volume as fluid enters or leaves the pore space, or as the compressible fluid within the pores changes in volume. The axiom of mass balance for the solid may be integrated in closed form to produce , where is the solid apparent density in the mixture reference configuration. Therefore, we may define the referential solid volume fraction as . In the finite element analysis, and are both specified by the user, then , , and are evaluated from the above relations, given a solution for (and thus, ). The solution for is obtained by solving for the nodal solid displacement vector , where represents the gradient operator in the material frame.
For the compressible fluid phase of a hybrid biphasic material, the true density varies with the intrinsic fluid volumetric strain, or dilatation, , where is the volume ratio of the fluid in its current and reference configurations, is the element mass of fluid and is the elemental volume of that fluid in the pores of the mixture, in the current configuration. It follows from this definition that where is the elemental volume of this same fluid (having the same elemental mass ) in the fluid's reference configuration. Thus, is a constant representing the true fluid density in its reference configuration, which is specified by the user. Accordingly, in the limit when the fluid is idealized to be intrinsically incompressible. This definition of is consistent with that for a pure fluid, as shown in our earlier formulation of computational fluid dynamics . As shown above, may be evaluated as , which now expands into an expression involving solid and fluid volume ratios, The limiting case when the fluid within the solid matrix pores has been completely squeezed out corresponds to (or ). Based on the above relation, this is equivalent to having the entire mixture volume reduce to the solid volume in the current configuration, or equivalently . In this limiting case, the mixture becomes an intrinsically incompressible solid and the finite element formulation presented in this study no longer applies.
We assume that (or equivalently, ) and at the start of a finite element analysis. The fluid dilatation is included as a nodal degree of freedom, implying that it is continuous across finite element boundaries. This assumption is verified below, when we review the jump conditions on axioms of mass, momentum and energy balance across interfaces.
Since the axiom of mass balance for the solid has already been solved in closed form, we only need to be concerned with the axiom of mass balance for the fluid, or alternatively, that of the mixture, which takes the form where is the solid velocity (the material time derivative of the solid displacement ), and is the volumetric flux of the fluid relative to the solid , with representing the fluid velocity. Here, is the material time derivative operator in the spatial frame, following the fluid motion. Since the finite element mesh is defined on the porous solid matrix of a biphasic mixture, and since the fluid flows relative to the solid, material time derivatives need to follow the solid motion in the finite element implementation. A similar scheme was used in our implementations of solute transport within deformable porous domains (Sections 2.7↑), to account for the motion of solutes relative to the porous solid matrix, as well as in our FSI formulation (Section 2.12↑) [11, 52, 53]. Thus, we substitute the following identity, into the axiom of mixture mass balance in (2.13.1-3) to produce the final form In this expression, the dot operator in represents the material time derivative following the solid motion. In the solid material frame, this material time derivative reduces to the partial time derivative. Accordingly, we may now write . In the BFSI implementation, the fluid volumetric flux relative to the solid is added to the list of nodal DOFs, giving us the complete set . This choice implies that is continuous across finite element boundaries, as justified further below when we review jump conditions across interfaces.
Based on the constitutive assumptions of our hybrid biphasic formulation , the momentum balance equations for the fluid and solid constituents reduce to and where is the acceleration, is the body force per mass acting on constituent , is the fluid pressure, is the apparent fluid viscous stress, and is the apparent solid elastic stress. These stress tensors are called apparent because their associated traction vectors represent a force acting on constituent per elemental mixture area. Here, is the hydraulic permeability tensor which regulates frictional drag between the fluid and solid constituents; setting to implies that this frictional drag is non-existent. Since the fluid is compressible, its pressure must be given by a function of state. In the isothermal framework used here, this function only depends on , and the form adopted in FEBio is where is the fluid bulk modulus (a user-specified material property). In the limit when inertia and body forces are neglected ( ), the fluid momentum balance (2.13.1-7) produces the classical Darcy-Brinkman relation . If the fluid viscous stress is also neglected ( ), we recover Darcy's law, .
Since the finite element implementation requires all material time derivatives to follow the motion of the solid, we recognize that and we evaluate the fluid acceleration as . However, since is not a nodal degree of freedom, we use eq.(2.13.1-4) to substitute into this expression. It follows that the fluid velocity gradient may be evaluated as where and . Now, the fluid acceleration takes the form