Fluid-structure interactions (FSI) and poroelasticity are useful continuum mechanics frameworks for modeling a variety of materials and engineering processes. The interaction of these two types of material domains, which we call fluid-porous-structure interactions (FPSI), has not been studied in detail. Currently, it is not possible to computationally model porous-permeable materials that can simultaneously incorporate finite deformation of the porous solid matrix, dynamics, and a viscous interstitial fluid, or any fluid exchange between such a porous domain and an adjoining viscous fluid domain. Such computational tools would be useful for of variety of applications, most notably in biomechanics of soft biological tissues and their interactions with biological fluids. Currently, researchers cannot model the interstitial fluid exchanged between synovial fluid and cartilage, blood and the arterial wall, or cerebrospinal fluid and the brain. Moreover, few computational tools exist that can solve independently for FSI and large deformation porous materials with viscous interstitial fluid, and even fewer can address FPSI.  Among these, none are available as open-source software, thus lacking the benefits of easy access, and the ability to replicate and disseminate models and results.  The objective of this study was to formulate and implement a FPSI solver in the free, open-source finite element code FEBio (febio.org). FEBio currently offers solvers for solid mechanics, computational fluid dynamics, and porous deformable materials \[1,2\]. The FPSI solver presented here is a generalization of the recently developed FSI solver that defined interactions at the interface between a deformable (non-porous) solid domain and a viscous fluid domain \[3\]. Previously, the deformable fluid domain in an FSI analysis was formulated as a specialized solid-fluid mixture where the mesh was defined on the solid material, with zero mass density and negligible stiffness to regularize its deformation. Here, this formulation is generalized to FPSI materials where the viscous fluid flows through a porous solid matrix that has non-zero mass density. In addition, to account for Darcy-Brinkman-like behavior, frictional interactions between the fluid and porous solid are incorporated using a hydraulic permeability tensor. Constitutive relations for the solid include hyperelastic and viscoelastic models, and those for the fluid include Newtonian and non-Newtonian models. The fluid constituent is assumed to be compressible, using its dilatation as a state variable and a suitable function of state to evaluate the pressure. For nearly incompressible responses, the fluid pressure is calculated from the negated product of its bulk modulus with the dilatation. Here, the mass balance equation is replaced with the equivalent kinematic constraint between the fluid dilatation and its velocity. Therefore, like the FSI solver, our FPSI formulation does not require stabilization methods to achieve good convergence. Governing equations for the fluid and solid constituents were rearranged to satisfy the jump conditions related to momentum balance across interfaces automatically, and as a result, two FPSI domains can be interfaced without requiring additional user-specified boundary conditions. A FPSI domain can be easily reduced to a standard poroelastic domain, or a deformable viscous fluid domain, as was verified using a variety of test problems. The finite element code was also successfully verified against various analytical solutions from the literature, such as Couette and Poiseuille-like flow where rigid-impermeable walls are replaced with porous deformable domains. A newly-derived analytical solution for 1D ultrafiltration of a viscous fluid through a poroelastic domain was used to further validate FEBio results. Finally, the FPSI solver was used to examine poro-elastohydrodynamic squeeze-film lubrication, non-Newtonian viscous blood flow within and through a porous deformable arterial wall, and a flow constriction problem where the deformable impermeable walls of a fluid domain can be squeezed shut. The successful formulation and implementation of this FPSI solver offers enhanced multiphysics modeling capabilities that are accessible via an open source software platform.

References: \[1\] Maas, SA et al., J Biomech Eng, 2012.

\[2\] Ateshian, GA et al., J Biomech Eng, 2018.

\[3\] Shim, JJ et al., J Biomech Eng, 2019.


A General Fluid-Porous-Structure Interactions Formulation Implemented in Febio

Division / Sponsors: 
Aerospace, Bioengineering, Applied Mechanics, Design Engineering Division, Advanced Energy Systems, Manufacturing Engineering Division, Noise Control and Acoustics, Nondestructive Evaluation, Diagnosis, and Prognosis (NDPD), Nuclear Engineering Division, Process Industries Division, Fluids Engineering, Heat Transfer, Materials, Microelectromechanical Systems, Safety Engineering and Risk Analysis, Computers & Information in Engineering


Paper Proceedings: https://asmedigitalcollection.asme.org/IMECE/IMECE2020/volume/84478

IMECE2020 - Virtual Conference

Supporting the need among researchers worldwide for additional online and multimedia content, ASME now offers the ASME Conference Video Collection. This platform will collect and provide long-term access to conference presentation videos, encompassing all technology areas within ASME's technical conference and events program. About 30,000-40,000 engineers across all mechanical engineering and related disciplines participate in ASME conferences each year. This includes authors of approximately 6,000 to 8,000 technical papers published by ASME, the majority of whom also submit video presentations for our virtual conferences. The video collection will be available as an add-on to existing product packages or as a separate subscription.