The Problem With National Institute of Standards and Technology Thermodynamics Tables in Continuum Mechanics.

Thermodynamics is a fundamental topic of continuum mechanics and biomechanics, with a wide range of applications to physiological and biological processes. This study addresses two fundamental limitations of current thermodynamic treatments. First, thermodynamics tables distributed online by the U.S. National Institute of Standards and Technology (NIST) report properties of fluids as a function of absolute temperature T and absolute pressure P. These properties include mass density ρ, specific internal energy u, enthalpy h=u+P/ρ, and entropy s. However, formulations of jump conditions across phase boundaries derived from Newton's second law of motion and the first law of thermodynamics employ the gauge pressure p=P-Pr, where Pr is an arbitrarily selected referential absolute pressure. Interchanging p with P is not innocuous as it alters tabulated NIST values for u while keeping h and s unchanged. Using p for functions of state and governing equations solves the problem with using NIST entries for the specific internal energy u in standard thermodynamics tables and analyses of phase transformation in continuum mechanics. Second, constitutive models for the free energy of fluids, such as water and air, are not typically provided in standard thermodynamics treatments. This study proposes a set of constitutive models and validates them against suitably modified NIST data.

Finite Element Implementation of Computational Fluid Dynamics With Reactive Neutral and Charged Solute Transport in FEBio.

The objective of this study was to implement a novel fluid-solutes solver into the open-source finite element software FEBio, that extended available modeling capabilities for biological fluids and fluid-solute mixtures. Using a reactive mixture framework, this solver accommodates diffusion, convection, chemical reactions, electrical charge effects, and external body forces, without requiring stabilization methods that were deemed necessary in previous computational implementations of the convection-diffusion-reaction equation at high Peclet numbers. Verification and validation problems demonstrated the ability of this solver to produce solutions for Peclet numbers as high as 1011, spanning the range of physiological conditions for convection-dominated solute transport. This outcome was facilitated by the use of a formulation that accommodates realistic values for solvent compressibility, and by expressing the solute mass balance such that it properly captured convective transport by the solvent and produced a natural boundary condition of zero diffusive solute flux at outflow boundaries. Since this numerical scheme was not necessarily foolproof, guidelines were included to achieve better outcomes that minimize or eliminate the potential occurrence of numerical artifacts. The fluid-solutes solver presented in this study represents an important and novel advancement in the modeling capabilities for biomechanics and biophysics as it allows modeling of mechanobiological processes via the incorporation of chemical reactions involving neutral or charged solutes within dynamic fluid flow. The incorporation of charged solutes in a reactive framework represents a significant novelty of this solver. This framework also applies to a broader range of nonbiological applications.

A Hybrid Biphasic Mixture Formulation for Modeling Dynamics in Porous Deformable Biological Tissues.

The primary aim of this study is to establish the theoretical foundations for a solid-fluid biphasic mixture domain that can accommodate inertial effects and a viscous interstitial fluid, which can interface with a dynamic viscous fluid domain. Most mixture formulations consist of constituents that are either all intrinsically incompressible or compressible, thereby introducing inherent limitations. In particular, mixtures with intrinsically incompressible constituents can only model wave propagation in the porous solid matrix, whereas those with compressible constituents require internal variables, and related evolution equations, to distinguish the compressibility of the solid and fluid under hydrostatic pressure. In this study, we propose a hybrid framework for a biphasic mixture where the skeleton of the porous solid is intrinsically incompressible but the interstitial fluid is compressible. We define a state variable as a measure of the fluid volumetric strain. Within an isothermal framework, the Clausius-Duhem inequality shows that a function of state arises for the fluid pressure as a function of this strain measure. We derive jump conditions across hybrid biphasic interfaces, which are suitable for modeling hydrated biological tissues. We then illustrate this framework using confined compression and dilatational wave propagation analyses. The governing equations for this hybrid biphasic framework reduce to those of the classical biphasic theory whenever the bulk modulus of the fluid is set to infinity and inertia terms and viscous fluid effects are neglected. The availability of this novel framework facilitates the implementation of finite element solvers for fluid-structure interactions at interfaces between viscous fluids and porous-deformable biphasic domains, which can include fluid exchanges across those interfaces.

A Hybrid Reactive Multiphasic Mixture With a Compressible Fluid Solvent.

Mixture theory is a general framework that has been used to model mixtures of solid, fluid, and solute constituents, leading to significant advances in modeling the mechanics of biological tissues and cells. Though versatile and applicable to a wide range of problems in biomechanics and biophysics, standard multiphasic mixture frameworks incorporate neither dynamics of viscous fluids nor fluid compressibility, both of which facilitate the finite element implementation of computational fluid dynamics solvers. This study formulates governing equations for reactive multiphasic mixtures where the interstitial fluid has a solvent which is viscous and compressible. This hybrid reactive multiphasic framework uses state variables that include the deformation gradient of the porous solid matrix, the volumetric strain and rate of deformation of the solvent, the solute concentrations, and the relative velocities between the various constituents. Unlike standard formulations which employ a Lagrange multiplier to model fluid pressure, this framework requires the formulation of a function of state for the pressure, which depends on solvent volumetric strain and solute concentrations. Under isothermal conditions the formulation shows that the solvent volumetric strain remains continuous across interfaces between hybrid multiphasic domains. Apart from the Lagrange multiplier-state function distinction for the fluid pressure, and the ability to accommodate viscous fluid dynamics, this hybrid multiphasic framework remains fully consistent with standard multiphasic formulations previously employed in biomechanics. With these additional features, the hybrid multiphasic mixture theory makes it possible to address a wider range of problems that are important in biomechanics and mechanobiology.

Finite Element Implementation of Biphasic-Fluid Structure Interactions in febio.

In biomechanics, solid-fluid mixtures have commonly been used to model the response of hydrated biological tissues. In cartilage mechanics, this type of mixture, where the fluid and solid constituents are both assumed to be intrinsically incompressible, is often called a biphasic material. Various physiological processes involve the interaction of a viscous fluid with a porous-hydrated tissue, as encountered in synovial joint lubrication, cardiovascular mechanics, and respiratory mechanics. The objective of this study was to implement a finite element solver in the open-source software febio that models dynamic interactions between a viscous fluid and a biphasic domain, accommodating finite deformations of both domains as well as fluid exchanges between them. For compatibility with our recent implementation of solvers for computational fluid dynamics (CFD) and fluid-structure interactions (FSI), where the fluid is slightly compressible, this study employs a novel hybrid biphasic formulation where the porous skeleton is intrinsically incompressible but the fluid is also slightly compressible. The resulting biphasic-FSI (BFSI) implementation is verified against published analytical and numerical benchmark problems, as well as novel analytical solutions derived for the purposes of this study. An illustration of this BFSI solver is presented for two-dimensional (2D) airflow through a simulated face mask under five cycles of breathing, showing that masks significantly reduce air dispersion compared to the no-mask control analysis. In addition, we model three-dimensional (3D) blood flow in a bifurcated carotid artery assuming porous arterial walls and verify that mass is conserved across all fluid-permeable boundaries. The successful formulation and implementation of this BFSI solver offers enhanced multiphysics modeling capabilities that are accessible via an open-source software platform.

Modeling Pulse Wave Propagation Through a Stenotic Artery With Fluid Structure Interaction: A Validation Study Using Ultrasound Pulse Wave Imaging.

Pulse wave imaging (PWI) is an ultrasound-based method that allows spatiotemporal mapping of the arterial pulse wave propagation, from which the local pulse wave velocity (PWV) can be derived. Recent reports indicate that PWI can help the assessment of atherosclerotic plaque composition and mechanical properties. However, the effect of the atherosclerotic plaque's geometry and mechanics on the arterial wall distension and local PWV remains unclear. In this study, we investigated the accuracy of a finite element (FE) fluid-structure interaction (FSI) approach to predict the velocity of a pulse wave propagating through a stenotic artery with an asymmetrical plaque, as quantified with PWI method. Experiments were designed to compare FE-FSI modeling of the pulse wave propagation through a stenotic artery against PWI obtained with manufactured phantom arteries made of polyvinyl alcohol (PVA) material. FSI-generated spatiotemporal maps were used to estimate PWV at the plaque region and compared it to the experimental results. Velocity of the pulse wave propagation and magnitude of the wall distension were correctly predicted with the FE analysis. In addition, findings indicate that a plaque with a high degree of stenosis (>70%) attenuates the propagation of the pulse pressure wave. Results of this study support the validity of the FE-FSI methods to investigate the effect of arterial wall structural and mechanical properties on the pulse wave propagation. This modeling method can help to guide the optimization of PWI to characterize plaque properties and substantiate clinical findings.

A Formulation for Fluid-Structure Interactions in febio Using Mixture Theory.

Many physiological systems involve strong interactions between fluids and solids, posing a significant challenge when modeling biomechanics. The objective of this study was to implement a fluid-structure interaction (FSI) solver in the free, open-source finite element code FEBio, that combined the existing solid mechanics and rigid body dynamics solver with a recently developed computational fluid dynamics (CFD) solver. A novel Galerkin-based finite element FSI formulation was introduced based on mixture theory, where the FSI domain was described as a mixture of fluid and solid constituents that have distinct motions. The mesh was defined on the solid domain, specialized to have zero mass, negligible stiffness, and zero frictional interactions with the fluid, whereas the fluid was modeled as isothermal and compressible. The mixture framework provided the foundation for evaluating material time derivatives in a material frame for the solid and in a spatial frame for the fluid. Similar to our recently reported CFD solver, our FSI formulation did not require stabilization methods to achieve good convergence, producing a compact set of equations and code implementation. The code was successfully verified against benchmark problems from the FSI literature and an analytical solution for squeeze-film lubrication. It was validated against experimental measurements of the flow rate in a peristaltic pump and illustrated using non-Newtonian blood flow through a bifurcated carotid artery with a thick arterial wall. The successful formulation and implementation of this FSI solver enhance the multiphysics modeling capabilities in febio relevant to the biomechanics and biophysics communities.

Finite Element Framework for Computational Fluid Dynamics in FEBio.

The mechanics of biological fluids is an important topic in biomechanics, often requiring the use of computational tools to analyze problems with realistic geometries and material properties. This study describes the formulation and implementation of a finite element framework for computational fluid dynamics (CFD) in FEBio, a free software designed to meet the computational needs of the biomechanics and biophysics communities. This formulation models nearly incompressible flow with a compressible isothermal formulation that uses a physically realistic value for the fluid bulk modulus. It employs fluid velocity and dilatation as essential variables: The virtual work integral enforces the balance of linear momentum and the kinematic constraint between fluid velocity and dilatation, while fluid density varies with dilatation as prescribed by the axiom of mass balance. Using this approach, equal-order interpolations may be used for both essential variables over each element, contrary to traditional mixed formulations that must explicitly satisfy the inf-sup condition. The formulation accommodates Newtonian and non-Newtonian viscous responses as well as inviscid fluids. The efficiency of numerical solutions is enhanced using Broyden's quasi-Newton method. The results of finite element simulations were verified using well-documented benchmark problems as well as comparisons with other free and commercial codes. These analyses demonstrated that the novel formulation introduced in FEBio could successfully reproduce the results of other codes. The analogy between this CFD formulation and standard finite element formulations for solid mechanics makes it suitable for future extension to fluid-structure interactions (FSIs).

Heterogeneous engineered cartilage growth results from gradients of media-supplemented active TGF-ß and is ameliorated by the alternative supplementation of latent TGF-ß.

Transforming growth factor beta (TGF-ß) has become one of the most widely utilized mediators of engineered cartilage growth. It is typically exogenously supplemented in the culture medium in its active form, with the expectation that it will readily transport into tissue constructs through passive diffusion and influence cellular biosynthesis uniformly. The results of this investigation advance three novel concepts regarding the role of TGF-ß in cartilage tissue engineering that have important implications for tissue development. First, through the experimental and computational analysis of TGF-ß concentration distributions, we demonstrate that, contrary to conventional expectations, media-supplemented exogenous active TGF-ß exhibits a pronounced concentration gradient in tissue constructs, resulting from a combination of high-affinity binding interactions and a high cellular internalization rate. These gradients are sustained throughout the entire culture duration, leading to highly heterogeneous tissue growth; biochemical and histological measurements support that while biochemical content is enhanced up to 4-fold at the construct periphery, enhancements are entirely absent beyond 1 mm from the construct surface. Second, construct-encapsulated chondrocytes continuously secrete large amounts of endogenous TGF-ß in its latent form, a portion of which undergoes cell-mediated activation and enhances biosynthesis uniformly throughout the tissue. Finally, motivated by these prior insights, we demonstrate that the alternative supplementation of additional exogenous latent TGF-ß enhances biosynthesis uniformly throughout tissue constructs, leading to enhanced but homogeneous tissue growth. This novel demonstration suggests that latent TGF-ß supplementation may be utilized as an important tool for the translational engineering of large cartilage constructs that will be required to repair the large osteoarthritic defects observed clinically.

Dynamic mechanical compression of devitalized articular cartilage does not activate latent TGF-ß.

A growing body of research has highlighted the role that mechanical forces play in the activation of latent TGF-ß in biological tissues. In synovial joints, it has recently been demonstrated that the mechanical shearing of synovial fluid, induced during joint motion, rapidly activates a large fraction of its soluble latent TGF-ß content. Based on this observation, the primary hypothesis of the current study is that the mechanical deformation of articular cartilage, induced by dynamic joint motion, can similarly activate the large stores of latent TGF-ß bound to the tissue extracellular matrix (ECM). Here, devitalized deep zone articular cartilage cylindrical explants (n=84) were subjected to continuous dynamic mechanical loading (low strain: ±2% or high strain: ±7.5% at 0.5Hz) for up to 15h or maintained unloaded. TGF-ß activation was measured in these samples over time while accounting for the active TGF-ß that remains bound to the cartilage ECM. Results indicate that TGF-ß1 is present in cartilage at high levels (68.5±20.6ng/mL) and resides predominantly in the latent form (>98% of total). Under dynamic loading, active TGF-ß1 levels did not statistically increase from the initial value nor the corresponding unloaded control values for any test, indicating that physiologic dynamic compression of cartilage is unable to directly activate ECM-bound latent TGF-ß via purely mechanical pathways and leading us to reject the hypothesis of this study. These results suggest that deep zone articular chondrocytes must alternatively obtain access to active TGF-ß through chemical-mediated activation and further suggest that mechanical deformation is unlikely to directly activate the ECM-bound latent TGF-ß of various other tissues, such as muscle, ligament, and tendon.