A nonlinear biphasic model of flow-controlled infusions in brain: Mass transport analyses

A biphasic nonlinear mathematical model is proposed for the mass transport that occurs during constant flow-rate infusions into brain tissue. The model takes into account geometric and material nonlinearities and a hydraulic conductivity dependent upon strain. The biphasic and convective–diffusive transport equations were implemented in a custom-written code assuming spherical symmetry and using an updated Lagrangian finite element algorithm. Results of the model indicate that the inclusion of these nonlinearities produced modest changes in the interstitial concentration but important variations in drug penetration and bulk concentration. Increased penetration of the drug but smaller bulk concentrations were obtained at smaller strains caused by combination of parameters such as increased Young’s modulus and initial hydraulic conductivity. This indicates that simulations of constant flow-rate infusions under the assumption of infinitesimal deformations or rigidity of the tissue may yield lower bulk concentrations near the infusion cavity and over-predictions of the penetration of the infused agent. The analyses also showed that decrease in the infusion flow rate of a fixed amount of drug results in increased penetration of the infused agent. From the clinical point-of-view, this may promote a safer infusion that delivers the therapeutic range over the desired volume while avoiding damage to the tissue by minimizing deformation and strain.     

Implications of Transvascular Fluid Exchange in Nonlinear,Biphasic Analyses of Flow-Controlled Infusion in Brain

A nonlinear, coupled biphasic-mass transport model that includes transvascular fluid exchange is proposed for flow-controlled infusions in brain tissue. The model accounts for geometric and material nonlinearities, a hydraulic conductivity dependent on deformation, and transvascular fluid exchange according to Starling’s law. The governing equations were implemented in a custom-written code assuming spherical symmetry and using an updated Lagrangian finite-element algorithm. Results of the model indicate that, using normal physiological values of vascular permeability, transvascular fluid exchange has negligible effects on tissue deformation, fluid pressure, and transport of the infused agent. As vascular permeability may be increased artificially through methods such as administering nitric oxide, a parametric study was conducted to determine how increased vascular permeability affects flow-controlled infusion. Increased vascular permeability reduced both tissue deformation and fluid pressure, possibly reducing damage to tissue adjacent to the infusion catheter. Furthermore, the loss of fluid to the vasculature resulted in a significantly increased interstitial fluid concentration but a modestly increased tissue concentration. From a clinical point of view, this increase in concentration could be beneficial if limited to levels below which toxicity would not occur. However, the modestly increased tissue concentration may make the increase in interstitial fluid concentration difficult to assess in vivo using co-infused radiolabeled agents.     

Flow-induced deformation from pressurized cavities in absorbing porous tissues

The behaviour of a cavity during an injection of fluid into biological tissue is considered. High cavity pressure drives fluid into the neighbouring tissue where it is absorbed by capillaries and lymphatics. The tissue is modelled as a nonlinear deformable porous medium with the injected fluid absorbed by the tissue at a rate proportional to the local pressure. A model with a spherical cavity in an infinite medium is used to find the pressure and displacement of the tissue as a function of time and radial distance. Analytical and numerical solutions for a step change in cavity pressure show that the flow induces a radial compression in the medium together with an annular expansion, the net result being an overall expansion of the medium. Thus any flow induced deformation of the material will aid in the absorption of fluid.     

Sensitivity analysis of permeability parameters of bovine nucleus pulposus obtained through inverse fitting of the nonlinear biphasic equation: effect of sampling strategy

Permeability controls the fluid flow into and out of soft tissue, and plays an important role in maintaining the health status of such tissue. Accurate determination of the parameters that define permeability is important for the interpretation of models that incorporate such processes. This paper describes the determination of strain-dependent permeability parameters from the nonlinear biphasic equation from experimental data of different sampling frequencies using the Nelder-Mead simplex method. The ability of this method to determine the global optimum was assessed by constructing the whole manifold arising from possible parameter combinations. Many parameter combinations yielded similar fits with the Nelder-Mead algorithm able to identify the global maximum within the resolution of the manifold. Furthermore, the sampling strategy affected the optimum values of the permeability parameters. Therefore, permeability parameter estimations arising from inverse methods should be utilised with the knowledge that they come with large confidence intervals.     

A simple method to estimate the exponential material parameters of heart valve tissue based on analogy between uniaxial tension and micropipette aspiration

Micropipette aspiration (MA) has been widely used to measure the biomechanical properties of cells and biomaterials. To estimate material parameters from MA experimental data, analytical half-space models and inverse finite element (FE) analyses are typically used. The half-space model is easy to implement but cannot account for nonlinear material properties and complex geometrical boundary conditions that are inherent to MA. Inverse FE approaches can account for geometrical and material nonlinearities, but their implementation is resource-intensive and not widely available. Here, by making analogy between an analytical uniaxial tension model and a FE model of MA, we proposed an easily implementable and accurate method to estimate the material parameters of tissues tested by MA. We first adopted a strain invariant-based isotropic exponential constitutive model and implemented it in both the analytical uniaxial tension model and the FE model. The two models were fit to experimental data generated by MA of porcine aortic valve tissue (45 spots on four leaflets) to estimate material parameters. We found no significant differences between the effective moduli estimated by the two models ( $p > 0.39$ ), with the effective moduli estimated by the uniaxial tension model correlating significantly with those estimated by the FE model ( $p < 0.001; R^{2}= 0.96$ ) with a linear regression slope that was not different than unity ( $p = 0.38$ ). Thus, the analytical uniaxial tension model, which avoids solving resource-intensive numerical problems, is as accurate as the FE model in estimating the effective modulus of valve tissue tested by MA.     

Strain energy descriptions of biological swelling. I: Single fluid compartment models

Strain energy functions are derived from biphasic soft tissue models in order to describe large-deformation, large-swelling, elastic behavior of nonlinear materials. The resulting analysis leads to calculations of stress-extension relations and tissue fluid pressure. Also explored are the elastic stability of the biphasic tissue models and the manner in which tissue pressure is altered by material deformation.     

Mathematical simulation of the body fluid regulating system in dog

A mathematical model of body fluid volume and osmolality regulation was developed which incorporated the major nonlinearities of fluid assimilation, exchange, distribution and excretion. The non-linear differential equations define compartmental material balances for water, urea, sodium, protein and antidiuretic hormone (ADH). The parameters of these equations were calculated using analytical solutions and available steady-state experimental data. The model was used to simulate the renal response to five input forcings: (1) intraesophageal water infusion; (2) water ingestion; (3) intravenous ADH injection; (4) intravenous water infusion; and (5) intermittent water loading. The model yielded continuous simulation curves which agreed reasonably well with the available transient and steady-state experimental data. The model predicted that stimulating volume receptors via changes in left atrial pressure accounts for only 15–20% of changes in ADH secretion rate, whereas stimulation of the osmotic receptors via changes in plasma osmolality accounts for the remaining 80–85% of changes. Thus, it appears that regulation of ADH secretion is largely dependent upon plasma osmolality during forcings which do not appreciably alter the cardiovascular blood volume.     

The validation of a generalized Hooke's law for coronary arteries

The exponential form of constitutive model is widely used in biomechanical studies of blood vessels. There are two main issues, however, with this model: 1) the curve fits of experimental data are not always satisfactory, and 2) the material parameters may be oversensitive. A new type of strain measure in a generalized Hooke's law for blood vessels was recently proposed by our group to address these issues. The new model has one nonlinear parameter and six linear parameters. In this study, the stress-strain equation is validated by fitting the model to experimental data of porcine coronary arteries. Material constants of left anterior descending artery and right coronary artery for the Hooke's law were computed with a separable nonlinear least-squares method with an excellent goodness of fit. A parameter sensitivity analysis shows that the stability of material constants is improved compared with the exponential model and a biphasic model. A boundary value problem was solved to demonstrate that the model prediction can match the measured arterial deformation under experimental loading conditions. The validated constitutive relation will serve as a basis for the solution of various boundary value problems of cardiovascular biomechanics.     

On the infusion of a therapeutic agent into a solid tumor modeled as a poroelastic medium

The direct infusion of an agent into a solid tumor, modeled as a spherical poroelastic material with anisotropic dependence of the tumor hydraulic conductivity upon the tissue deformation, is treated both by solving the coupled fluid/elastic equations, and by expressing the solution as an asymptotic expansion in terms of a small parameter, ratio between the driving pressure force in the fluid system, and the elastic properties of the solid. Results at order one match almost perfectly the solutions of the full system over a large range of infusion pressures. Comparison with experimental results is acceptable after the hydraulic conductivity of the medium is properly calibrated. Given the uncertain estimates of some model constants, the order zero solution of the expansion, for which fluid and porous matrix are decoupled, yields acceptable values and trends for all the physical fields of interest, rendering the coupled analysis (in the limit of small displacements) of little use. When the deformation of the tissue becomes large nonlinear elasticity theory must be resorted to.     

Comparison of biophysical stimuli for mechano-regulation of tissue differentiation during fracture healing

Most long-bone fractures heal through indirect or secondary fracture healing, a complex process in which endochondral ossification is an essential part and bone is regenerated by tissue differentiation. This process is sensitive to the mechanical environment, and several authors have proposed mechano-regulation algorithms to describe it using strain, pore pressure and/or interstitial fluid velocity as biofeedback variables. The aim of this study was to compare various mechano-regulation algorithms' abilities to describe normal fracture healing in one computational model. Additionally, we hypothesized that tissue differentiation during normal fracture healing could be equally well regulated by the individual mechanical stimuli, e.g. deviatoric strain, pore pressure or fluid velocity. A biphasic finite element model of an ovine tibia with a 3mm fracture gap and callus was used to simulate the course of tissue differentiation during normal fracture healing. The load applied was regulated in a biofeedback loop, where the load magnitude was determined by the interfragmentary movement in the fracture gap. All the previously published mechano-regulation algorithms studied, simulated the course of normal fracture healing correctly. They predicted (1) intramembranous bone formation along the periosteum and callus tip, (2) endochondral ossification within the external callus and cortical gap, and (3) creeping substitution of bone towards the gap from the initial lateral osseous bridge. Some differences between the effects of the algorithms were seen, but they were not significant. None of the volumetric components, i.e. pore pressure or fluid velocity, alone were able to correctly predict spatial or temporal tissue distribution during fracture healing. However, simulation as a function of only deviatoric strain accurately predicted the course of normal fracture healing. This suggests that the deviatoric component may be the most significant mechanical parameter to guide tissue differentiation during indirect fracture healing.     

Wear particle diffusion and tissue differentiation in TKA implant fibrous interfaces

In the context of mechanical loosening, we studied the hypothesis that wear-particle migration in the fibrous membrane under tibial plateaus after total knee arthroplasty can be explained by the pumping effects of the interstitial fluid in the tissue. Further, as a secondary objective we investigated the possibility that interface-tissue differentiation is influenced by interstitial fluid flow and strain, as mechanical effects of interface motions. For comparative reasons, we analyzed a previously published simplified two-dimensional finite-element model, this time assuming biphasic tissue properties. We wanted to determine hydrostatic pressure and flow velocities in the fluid phase, in addition to stresses and strains, for time-dependent loading of the plateau. We found that fluid flow in the interface was extremely slow, except in the periphery. Hence, loosening due to particle-induced bone resorption appears improbable. The results, however, do support the idea that particles migrate with fluid flow, when such flow occurs. Where fibrous tissue tends to be prominent in reality, the fluid is repeatedly extruded and reabsorbed in the model. Where these values are low, fibrocartilage is commonly found. When material properties were varied to subsequently represent fibrocartilage and two stages of mineralization, the strains and fluid velocities is reduced. Fluid pressure, however, did not change. Our results refute the hypothesis that wear particles are pumped through the interface tissue below a TKA but support the hypothesis that interface tissue type and loosening processes are influenced by mechanical tissue variables such as tissue strain and interstitial fluid velocity.     

A comparison between mechano-electrochemical and biphasic swelling theories for soft hydrated tissues

Biological tissues like intervertebral discs and articular cartilage primarily consist of interstitial fluid, collagen fibrils and negatively charged proteoglycans. Due to the fixed charges of the proteoglycans, the total ion concentration inside the tissue is higher than in the surrounding synovial fluid (cation concentration is higher and the anion concentration is lower). This excess of ion particles leads to an osmotic pressure difference, which causes swelling of the tissue. In the last decade several mechano-electrochemical models, which include this mechanism, have been developed. As these models are complex and computationally expensive, it is only possible to analyze geometrically relatively small problems. Furthermore, there is still no commercial finite element tool that includes such a mechano-electrochemical theory. Lanir (Biorheology, 24, pp. 173-187, 1987) hypothesized that electrolyte flux in articular cartilage can be neglected in mechanical studies. Lanir's hypothesis implies that the swelling behavior of cartilage is only determined by deformation of the solid and by fluid flow. Hence, the response could be described by adding a deformation-dependent pressure term to the standard biphasic equations. Based on this theory we developed a biphasic swelling model. The goal of the study was to test Lanir's hypothesis for a range of material properties. We compared the deformation behavior predicted by the biphasic swelling model and a full mechano-electrochemical model for confined compression and 1D swelling. It was shown that, depending on the material properties, the biphasic swelling model behaves largely the same as the mechano-electrochemical model, with regard to stresses and strains in the tissue following either mechanical or chemical perturbations. Hence, the biphasic swelling model could be an alternative for the more complex mechano-electrochemical model, in those cases where the ion flux itself is not the subject of the study. We propose thumbrules to estimate the correlation between the two models for specific problems.     

Embryonic brain enlargement requires cerebrospinal fluid pressure

The brain of the chick embryo begins to enlarge abruptly on the second day of incubation. Shortly thereafter, major flexures and torsions of the brain occur, and many bulges and furrows appear. The onset of enlargement coincides with closure of the spinal canal which makes the neural tube a closed compartment filled with cerebrospinal fluid. We propose that cerebrospinal fluid pressure is a necessary driving force for normal brain enlargement. We have experimentally tested this hypothesis by intubating brains of chick embryos and comparing brain cavity and tissue volumes in normal and intubated embryos. The increase in cavity volume is greatly reduced, whereas brain tissue continues to grow at a reduced rate and folds into the ventricles.     

Biphasic indentation of articular cartilage--I. Theoretical analysis

A mathematical solution has been obtained for the indentation creep and stress-relaxation behavior of articular cartilage where the tissue is modeled as a layer of linear KLM biphasic material of thickness h bonded to an impervious, rigid bony substrate. The circular (radius = a), plane-ended indenter is assumed to be rigid, porous, free-draining, and frictionless. Double Laplace and Hankel transform techniques were used to solve the partial differential equations. The transformed equations and boundary conditions yielded an integral equation of the Fredholm type which was analyzed asymptotically and solved numerically. Our asymptotic analyses showed that the linear KLM biphasic material behaves like an incompressible (v = 0.5) single-phase elastic solid at t = 0+; the instantaneous response of the material is governed by the shear modulus (mu s) of the solid matrix. The linear KLM biphasic material behaves like a compressible elastic solid with material properties defined by those of the solid matrix, i.e. (lambda s, mu s) or (mu s, v s) as t----infinity. The transient viscoelastic creep and stress-relaxation behavior, 0 less than t less than infinity, of this material is controlled by the frictional drag (which is inversely proportional to the permeability k) associated with the flow of the interstitial fluid through the porous-permeable solid matrix. For given values of the Poisson's ratio of the solid matrix v s and the aspect ratio a/h, where a is the radius of the indenter and h is the thickness of the layer, the creep behavior with respect to the dimensionless time H Akt/a2 is completely controlled by the load parameter P/2 mu sa2 and the stress relaxation behavior is completely controlled by the rate of compression parameter R0 = kH A/V0h where H A = lambda s + 2 mu s and the equilibrium strain u0/h. This mathematical solution may now be used to describe an indentation experiment on articular cartilage to determine the intrinsic material properties of the tissue, i.e. permeability k, and the elastic coefficients of the solid phase (lambda s, mu s) or (mu s, v s).     

A microstructurally based continuum model of cartilage viscoelasticity and permeability incorporating measured statistical fiber orientations

The remarkable mechanical properties of cartilage derive from an interplay of isotropically distributed, densely packed and negatively charged proteoglycans; a highly anisotropic and inhomogeneously oriented fiber network of collagens; and an interstitial electrolytic fluid. We propose a new 3D finite strain constitutive model capable of simultaneously addressing both solid (reinforcement) and fluid (permeability) dependence of the tissue’s mechanical response on the patient-specific collagen fiber network. To represent fiber reinforcement, we integrate the strain energies of single collagen fibers—weighted by an orientation distribution function (ODF) defined over a unit sphere—over the distributed fiber orientations in 3D. We define the anisotropic intrinsic permeability of the tissue with a structure tensor based again on the integration of the local ODF over all spatial fiber orientations. By design, our modeling formulation accepts structural data on patient-specific collagen fiber networks as determined via diffusion tensor MRI. We implement our new model in 3D large strain finite elements and study the distributions of interstitial fluid pressure, fluid pressure load support and shear stress within a cartilage sample under indentation. Results show that the fiber network dramatically increases interstitial fluid pressure and focuses it near the surface. Inhomogeneity in the tissue’s composition also increases fluid pressure and reduces shear stress in the solid. Finally, a biphasic neo-Hookean material model, as is available in commercial finite element codes, does not capture important features of the intra-tissue response, e.g., distributions of interstitial fluid pressure and principal shear stress.     

Squeeze-film lubrication of the human ankle joint with synovial fluid filtrated by articular cartilage with the superficial zone worn out

A squeeze-film lubrication model of the human ankle joint in standing that takes into account the fluid transport across the articular surface is presented. Articular cartilage is a biphasic mixture of the ideal interstitial fluid and an elastic permeable isotropic homogeneous intrinsically incompressible matrix. The simple homogeneous model for articular cartilage models the case of early osteoarthritis, when the intact superficial zone of the normal articular cartilage, much stiffer in tension than the bulk material, has been already disrupted or worn out. The calculations indicate for this case that in normal approach motion the lubricating fluid film is quickly depleted and turned into a synovial gel film that is supposed to serve as a boundary lubricant if sliding motion follows