A nonlinear biphasic model of flow-controlled infusion in brain: fluid transport and tissue deformation analyses

A biphasic nonlinear mathematical model is proposed for the concomitant fluid transport and tissue deformation that occurs during constant flow rate infusions into brain tissue. The model takes into account material and geometrical nonlinearities, a hydraulic conductivity dependent on strain, and nonlinear boundary conditions at the infusion cavity. The biphasic equations were implemented in a custom written code assuming spherical symmetry and using an updated Lagrangian finite element algorithm. Results of the model showed that both, geometric and material nonlinearities play an important role in the physics of infusions, yielding important differences from infinitesimal analyses. Geometrical nonlinearities were mainly due to the significant enlargement of the infusion cavity, while variations of the parameters that describe the degree of nonlinearity of the stress-strain curve yielded significant differences in all distributions. For example, a parameter set showing stiffening under tension yielded maximum values of radial displacement and porosity not localized at the infusion cavity. On the other hand, a parameter set showing softening under tension yielded a slight decrease in the fluid velocity for a three-fold increase in the flow rate, which can be explained by the substantial increase of the infusion cavity, not considered in linear analyses. This study strongly suggests that significant enlargement of the infusion cavity is a real phenomenon during infusions that may produce collateral damage to brain tissue. Our results indicate that more experimental tests have to be undertaken in order to determine material nonlinearities of brain tissue over a range of strains. With better understanding of these nonlinear effects, clinicians may be able to develop protocols that can minimize the damage to surrounding tissue.