A Computational Pipeline to Investigate Longitudinal Blood Flow Changes in the Circle of Willis of Patients with Stable and Growing Aneurysms.

Changes in cerebral blood flow are often associated with the initiation and development of different life-threatening medical conditions including aneurysm rupture and ischemic stroke. Nevertheless, it is not fully clear how haemodynamic changes in time across the Circle of Willis (CoW) are related with intracranial aneurysm (IA) growth. In this work, we introduced a novel reduced-order modelling strategy for the systematic quantification of longitudinal blood flow changes across the whole CoW in patients with stable and unstable/growing aneurysm. Magnetic Resonance Angiography (MRA) images were converted into one-dimensional (1-D) vessel networks through a semi-automated procedure, with a level of geometric reconstruction accuracy controlled by user-dependent parameters. The proposed pipeline was used to systematically analyse longitudinal haemodynamic changes in seven different clinical cases. Our preliminary simulation results indicate that growing aneurysms are not necessarily associated with significant changes in mean flow over time. A concise sensitivity analysis also shed light on which modelling aspects need to be further characterized to have reliable patient-specific predictions. This study poses the basis for investigating how time-dependent changes in the vasculature affect the haemodynamics across the whole CoW in patients with stable and growing aneurysms.

A new model for evaluating pressure-induced vascular tone in small cerebral arteries.

The capacity of small cerebral arteries (SCAs) to adapt to pressure fluctuations has a fundamental physiological role and appears to be relevant in different pathological conditions. Here, we present a new computational model for quantifying the link, and its contributors, between luminal pressure and vascular tone generation in SCAs. This is assembled by combining a chemical sub-model, representing pressure-induced smooth muscle cell (SMC) signalling, with a mechanical sub-model for the tone generation and its transduction at tissue level. The devised model can accurately reproduce the impact of luminal pressure on different cytoplasmic components involved in myogenic signalling, both in the control case and when combined with some specific pharmacological interventions. Furthermore, the model is also able to capture and predict experimentally recorded pressure-outer diameter relationships obtained for vessels under control conditions, both in a Ca 2 + -free bath and under drug inhibition. The modularity of the proposed framework allows the integration of new components for the study of a broad range of processes involved in the vascular function.

On the Ca2＋ elevation in vascular endothelial cells due to inositol trisphosphate-sensitive store receptors activation: A data-driven modeling approach.

Agonist-induced Ca2+ signaling is essential for the regulation of many vital functions in endothelial cells (ECs). A broad range of stimuli elevate the cytosolic Ca2＋ concentration by promoting a pathway mediated by inositol 1,4,5 trisphosphate (IP3) which causes Ca2＋ release from intracellular stores. Despite its importance, there are very few studies focusing on the quantification of such dynamics in the vascular endothelium. Here, by using data from isolated ECs, we established a minimalistic modeling framework able to quantitatively capture the main features (averaged over a cell population) of the cytosolic Ca2＋ response to different IP3 stimulation levels. A suitable description of Ca2＋-regulatory function of inositol 1,4,5 trisphosphate receptors (IP3Rs) and corresponding parameter space are identified by comparing the different model variants against experimental mean population data. The same approach is used to numerically assess the relevance of cytosolic Ca2＋ buffering, as well as Ca2＋ store IP3-sensitivity in the overall cell dynamics. The variability in the dynamics' features observed across the population can be explained (at least in part) through variation of certain model parameters (such as buffering capacity or Ca2＋ store sensitivity to IP3). The results, in terms of experimental fitting and validation, support the proposed minimalistic model as a reference framework for the quantification of the EC Ca2＋ dynamics induced by IP3Rs activation.

Modeling Reactive Hyperemia to Better Understand and Assess Microvascular Function: A Review of Techniques.

Reactive hyperemia is a well-established technique for the non-invasive evaluation of the peripheral microcirculatory function, measured as the magnitude of limb re-perfusion after a brief period of ischemia. Despite widespread adoption by researchers and clinicians alike, many uncertainties remain surrounding interpretation, compounded by patient-specific confounding factors (such as blood pressure or the metabolic rate of the ischemic limb). Mathematical modeling can accelerate our understanding of the physiology underlying the reactive hyperemia response and guide in the estimation of quantities which are difficult to measure experimentally. In this work, we aim to provide a comprehensive guide for mathematical modeling techniques that can be used for describing the key phenomena involved in the reactive hyperemia response, alongside their limitations and advantages. The reported methodologies can be used for investigating specific reactive hyperemia aspects alone, or can be combined into a computational framework to be used in (pre-)clinical settings.

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.

On the poro-elastic models for microvascular blood flow resistance: An in vitro validation.

Nowadays, adequate and accurate representation of the microvascular flow resistance constitutes one of the major challenges in computational haemodynamic studies. In this work, a theoretical, porous media framework, ultimately designed for representing downstream resistance, is presented and compared against an in vitro experimental results. The resistor consists of a poro-elastic tube, with either a constant or variable porosity profile in space. The underlying physics, characterizing the fluid flow through the porous media, is analysed by considering flow variables at different network locations. Backward reflections, originated in the reservoir of the in vitro model, are accounted for through a reflection coefficient imposed as an outflow network condition. The simulation results are in good agreement with the measurements for both the homogenous and heterogeneous porosity conditions. In addition, the comparison allows identification of the range of values representing experimental reservoir reflection coefficients. The pressure drops across the heterogeneous porous media increases with respect to the simpler configuration, whilst flow remains almost unchanged. The effect of some fluid network features, such as tube Young's modulus and fluid viscosity, on the theoretical results is also elucidated, providing a reference for the invitro and insilico simulation of different microvascular conditions.

A novel numerical modelling approach for keratoplasty eye procedure.

Objective of the work is to investigate stress and deformation that conrneal tissue and donor graft undergo during endothelial keratoplasty. In order to attach the donor graft to the cornea, different air bubble pressure profiles acting on the graft are considered. This study is carried out by employing a three-dimensional nonlinear finite element methodology, combined with a contact algorithm. The ocular tissues are treated as isotropic, hyper-elastic and nearly-incompressible materials. The contact algorithm, based on the penalty-based node-to-surface approach, is used to model the donor graft-corneal interface region. First, the proposed computational methodology is tested against benchmark data for bending of the plates over a cylinder. Then, the influence of geometrical and material parameters of the graft on the corneal contact-structural response is investigated. The results are presented in terms of Von Mises stress intensity, displacement and mean contact force. Results clearly indicate that the air bubble pressure plays a key role in the corneal stress and strain, as well as graft stiffness and thickness.

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.

Influence of ageing on human body blood flow and heat transfer: A detailed computational modelling study.

Ageing plays a fundamental role in arterial blood transport and heat transfer within a human body. The aim of this work is to provide a comprehensive methodology, based on biomechanical considerations, for modelling arterial flow and energy exchange mechanisms in the body accounting for age-induced changes. The study outlines a framework for age-related modifications within several interlinked subsystems, which include arterial stiffening, heart contractility variations, tissue volume and property changes, and thermoregulatory system deterioration. Some of the proposed age-dependent governing equations are directly extrapolated from experimental data sets. The computational framework is demonstrated through numerical experiments, which show the impact of such age-related changes on arterial blood pressure, local temperature distribution, and global body thermal response. The proposed numerical experiments show that the age-related changes in arterial convection do not significantly affect the tissue temperature distribution. Results also highlight age-related effects on the sweating mechanism, which lead to a significant reduction in heat dissipation and a subsequent rise in skin and core temperatures.

A multiscale active structural model of the arterial wall accounting for smooth muscle dynamics.

Arterial wall dynamics arise from the synergy of passive mechano-elastic properties of the vascular tissue and the active contractile behaviour of smooth muscle cells (SMCs) that form the media layer of vessels. We have developed a computational framework that incorporates both these components to account for vascular responses to mechanical and pharmacological stimuli. To validate the proposed framework and demonstrate its potential for testing hypotheses on the pathogenesis of vascular disease, we have employed a number of pharmacological probes that modulate the arterial wall contractile machinery by selectively inhibiting a range of intracellular signalling pathways. Experimental probes used on ring segments from the rabbit central ear artery are: phenylephrine, a selective α1-adrenergic receptor agonist that induces vasoconstriction; cyclopiazonic acid (CPA), a specific inhibitor of sarcoplasmic/endoplasmic reticulum Ca2+-ATPase; and ryanodine, a diterpenoid that modulates Ca2+ release from the sarcoplasmic reticulum. These interventions were able to delineate the role of membrane versus intracellular signalling, previously identified as main factors in smooth muscle contraction and the generation of vessel tone. Each SMC was modelled by a system of nonlinear differential equations that account for intracellular ionic signalling, and in particular Ca2+ dynamics. Cytosolic Ca2+ concentrations formed the catalytic input to a cross-bridge kinetics model. Contractile output from these cellular components forms the input to the finite-element model of the arterial rings under isometric conditions that reproduces the experimental conditions. The model does not account for the role of the endothelium, as the nitric oxide production was suppressed by the action of L-NAME, and also due to the absence of shear stress on the arterial ring, as the experimental set-up did not involve flow. Simulations generated by the integrated model closely matched experimental observations qualitatively, as well as quantitatively within a range of physiological parametric values. The model also illustrated how increased intercellular coupling led to smooth muscle coordination and the genesis of vascular tone.

Modelling accidental hypothermia effects on a human body under different pathophysiological conditions.

Accidental exposure to cold water environment is one of the most challenging situations in which hypothermia occurs. In the present work, we aim to characterise the energy balance of a human body subjected to such extreme environmental conditions. This study is carried out using a recently developed computational model and by setting boundary conditions needed to simulate the effect of cold surrounding environment. A major finding is the capacity of the body core regions to maintain their temperature high for a substantial amount of time, even under the most extreme environmental conditions. We also considered two disease states that highlight the spectrum of possible pathologies implicated in thermal regulation of the human body. These states are (i) cardiomyopathy, which affects the operating capacity of the heart, and (ii) malnutrition, which directly impairs the body's ability to regulate heat exchange with the environment. We have found that cardiomyopathy has little influence on the thermal balance of the human body, whereas malnutrition has a profound negative effect on the thermal balance and leads to dramatic reduction in core temperature.

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.