A novel numerical modelling approach for keratoplasty eye procedure

Biomechanics and Modeling in Mechanobiology - Objective of the work is to investigate stress and deformation that conrneal tissue and donor graft undergo during endothelial keratoplasty. In order...     

A novel porous media-based approach to outflow boundary resistances of 1D arterial blood flow models

In this paper we introduce a novel method for prescribing terminal boundary conditions in one-dimensional arterial flow networks. This is carried out by coupling the terminal arterial vessel with a poro-elastic tube, representing the flow resistance offered by microcirculation. The performance of the proposed porous media-based model has been investigated through several different numerical examples. First, we investigate model parameters that have a profound influence on the flow and pressure distributions of the system. The simulation results have been compared against the waveforms generated by three elements (RCR) Windkessel model. The proposed model is also integrated into a realistic arterial tree, and the results obtained have been compared against experimental data at different locations of the network. The accuracy and simplicity of the proposed model demonstrates that it can be an excellent alternative for the existing models.

     

A framework for incorporating 3D hyperelastic vascular wall models in 1D blood flow simulations

We present a novel framework for investigating the role of vascular structure on arterial haemodynamics in large vessels, with a special focus on the human common carotid artery (CCA). The analysis is carried out by adopting a three-dimensional (3D) derived, fibre-reinforced, hyperelastic structural model, which is coupled with an axisymmetric, reduced order model describing blood flow. The vessel transmural pressure and lumen area are related via a Holzapfel–Ogden type of law, and the residual stresses along the thickness and length of the vessel are also accounted for. After a structural characterization of the adopted hyperelastic model, we investigate the link underlying the vascular wall response and blood-flow dynamics by comparing the proposed framework results against a popular tube law. The comparison shows that the behaviour of the model can be captured by the simpler linear surrogate only if a representative value of compliance is applied. Sobol’s multi-variable sensitivity analysis is then carried out in order to identify the extent to which the structural parameters have an impact on the CCA haemodynamics. In this case, the local pulse wave velocity (PWV) is used as index for representing the arterial transmission capacity of blood pressure waveforms. The sensitivity analysis suggests that some geometrical factors, such as the stress-free inner radius and opening angle, play a major role on the system’s haemodynamics. Subsequently, we quantified the differences in haemodynamic variables obtained from different virtual CCAs, tube laws and flow conditions. Although each artery presents a distinct vascular response, the differences obtained across different flow regimes are not significant. As expected, the linear tube law is unable to accurately capture all the haemodynamic features characterizing the current model. The findings from the sensitivity analysis are further confirmed by investigating the axial stretching effect on the CCA fluid dynamics. This factor does not seem to alter the pressure and flow waveforms. On the contrary, it is shown that, for an axially stretched vessel, the vascular wall exhibits an attenuation in absolute distension and an increase in circumferential stress, corroborating the findings of previous studies. This analysis shows that the new model offers a good balance between computational complexity and physics captured, making it an ideal framework for studies aiming to investigate the profound link between vascular mechanobiology and blood flow.

     

An advanced computational bioheat transfer model for a human body with an embedded systemic circulation

In the present work, an elaborate one-dimensional thermofluid model for a human body is presented. By contrast to the existing pure conduction-/perfusion-based models, the proposed methodology couples the arterial fluid dynamics of a human body with a multi-segmental bioheat model of surrounding solid tissues. In the present configuration, arterial flow is included through a network of elastic vessels. More than a dozen solid segments are employed to represent the heat conduction in the surrounding tissues, and each segment is constituted by a multilayered circular cylinder. Such multi-layers allow flexible delineation of the geometry and incorporation of properties of different tissue types. The coupling of solid tissue and fluid models requires subdivision of the arterial circulation into large and small arteries. The heat exchange between tissues and arterial wall occurs by convection in large vessels and by perfusion in small arteries. The core region, including the heart, provides the inlet conditions for the fluid equations. In the proposed model, shivering, sweating, and perfusion changes constitute the basis of the thermoregulatory system. The equations governing flow and heat transfer in the circulatory system are solved using a locally conservative Galerkin approach, and the heat conduction in the surrounding tissues is solved using a standard implicit backward Euler method. To investigate the effectiveness of the proposed model, temperature field evolutions are monitored at different points of the arterial tree and in the surrounding tissue layers. To study the differences due to flow-induced convection effects on thermal balance, the results of the current model are compared against those of the widely used modelling methodologies. The results show that the convection significantly influences the temperature distribution of the solid tissues in the vicinity of the arteries. Thus, the inner convection has a more predominant role in the human body heat balance than previously thought. To demonstrate its capabilities, the proposed new model is used to study different scenarios, including thermoregulation inactivity and variation in surrounding atmospheric conditions.