Outflow boundary conditions for 3D simulations of non-periodic blood flow and pressure fields in deformable arteries

The simulation of blood flow and pressure in arteries requires outflow boundary conditions that incorporate models of downstream domains. We previously described a coupled multidomain method to couple analytical models of the downstream domains with 3D numerical models of the upstream vasculature. This prior work either included pure resistance boundary conditions or impedance boundary conditions based on assumed periodicity of the solution. However, flow and pressure in arteries are not necessarily periodic in time due to heart rate variability, respiration, complex transitional flow or acute physiological changes. We present herein an approach for prescribing lumped parameter outflow boundary conditions that accommodate transient phenomena. We have applied this method to compute haemodynamic quantities in different physiologically relevant cardiovascular models, including patient-specific examples, to study non-periodic flow phenomena often observed in normal subjects and in patients with acquired or congenital cardiovascular disease. The relevance of using boundary conditions that accommodate transient phenomena compared with boundary conditions that assume periodicity of the solution is discussed

Pulsatile non-newtonian blood flow in image-based models of carotid bifurcation

Present hemodynamical study is motivated by the ongoing clinical research at the University Hospital in Pilsen, Czech Republic. On the basis of provided CT scans, several carotid artery models were reconstructed and used for this numerical study of pulsatile blood flow. The blood is modelled as a shear-dependent incompressible fluid, motion of which is described by the non-linear system of Navier-Stokes equations coupled with the Carreau-Yasuda model. The mathematical model is solved using in-house software based on the principle of the SIMPLE algorithm and cell-centred finite volume method (FVM) formulated for hybrid unstructured tetrahedral grids. The discussion of obtained numerical results is performed with special emphasis placed on the analysis of velocity field and distribution of main hemodynamic factors such as cycle-averaged WSS and oscillatory shear index (OSI) in areas prone to atherosclerosis

Computational Simulations for Aortic Coarctation: Representative Results From a Sampling of Patients

Treatments for coarctation of the aorta (CoA) can alleviate blood pressure (BP) gradients(D), but long-term morbidity still exists that can be explained by altered indices of hemodynamics and biomechanics. We introduce a technique to increase our understanding of these indices for CoA under resting and nonresting conditions, quantify their contribution to morbidity, and evaluate treatment options. Patient-specific computational fluid dynamics (CFD) models were created from imaging and BP data for one normal and four CoA patients (moderate native CoA: D12 mmHg, severe native CoA: D25 mmHg and postoperative end-to-end and end-to-side patients: D0 mmHg). Simulations incorporated vessel deformation, downstream vascular resistance and compliance. Indices including cyclic strain, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) were quantified. Simulations replicated resting BP and blood flow data. BP during simulated exercise for the normal patient matched reported values. Greatest exercise-induced increases in systolic BP and mean and peak DBP occurred for the moderate native CoA patient (SBP: 115 to 154 mmHg; mean and peak DBP: 31 and 73 mmHg). Cyclic strain was elevated proximal to the coarctation for native CoA patients, but reduced throughout the aorta after treatment. A greater percentage of vessels was exposed to subnormal TAWSS or elevated OSI for CoA patients. Local patterns of these indices reported to correlate with atherosclerosis in normal patients were accentuated by CoA. These results apply CFD to a range of CoA patients for the first time and provide the foundation for future progress in this area

Mathematical modeling of local perfusion in large distensible microvascular networks

Microvessels -blood vessels with diameter less than 200 microns- form large, intricate networks organized into arterioles, capillaries and venules. In these networks, the distribution of flow and pressure drop is a highly interlaced function of single vessel resistances and mutual vessel interactions. In this paper we propose a mathematical and computational model to study the behavior of microcirculatory networks subjected to different conditions. The network geometry is composed of a graph of connected straight cylinders, each one representing a vessel. The blood flow and pressure drop across the single vessel, further split into smaller elements, are related through a generalized Ohm's law featuring a conductivity parameter, function of the vessel cross section area and geometry, which undergo deformations under pressure loads. The membrane theory is used to describe the deformation of vessel lumina, tailored to the structure of thick-walled arterioles and thin-walled venules. In addition, since venules can possibly experience negative transmural pressures, a buckling model is also included to represent vessel collapse. The complete model including arterioles, capillaries and venules represents a nonlinear system of PDEs, which is approached numerically by finite element discretization and linearization techniques. We use the model to simulate flow in the microcirculation of the human eye retina, a terminal system with a single inlet and outlet. After a phase of validation against experimental measurements, we simulate the network response to different interstitial pressure values. Such a study is carried out both for global and localized variations of the interstitial pressure. In both cases, significant redistributions of the blood flow in the network arise, highlighting the importance of considering the single vessel behavior along with its position and connectivity in the network

A structural approach including the behavior of collagen cross-links to model patient-specific human carotid arteries

The final publication is available at Springer via http://dx.doi.org/10.1007/s10439-014-0995-7The objective of this work is to develop a remodeling model for biological matter coupling two different processes in a 3D framework: reorientation of the preferential direction of a given fibered structure and reorientation of the fibrils or filaments that make up such a structure. This work uses the microsphere-based approach to take into account the micro mechanics involved in biological fibered structures regarding both their passive behavior and the reorientation of their micro constituents. Moreover, the macro behavior of the material as a whole is obtained by means of homogenizing the underlying micro response. We associate the orientation space of the integration directions to the physical space of micro-fibrils. To approximate the directional distribution of the fibrils within each fiber bundle, a Bingham probability orientation density function is introduced into the Helmholtz energy function. With all these assumptions, the problem is studied from an energetic point of view, describing the dissipation inherent to remodeling processes, and the evolution equations for both reorientations (change in preferential direction of the network and change in shape of the fibril distribution) re obtained. The model is included in a finite element code which allows computing different geometries and boundary value problems. This results in a complete methodology for characterizing the reorientation evolution of different fibered biological structures, such as cells. Our results show remodeling of fibered structures in two different scales, presenting a qualitatively good agreement with experimental findings in cell mechanics. Hierarchical structures align in the direction of the maximum principal direction of the considered stimulus and narrow in the perpendicular direction. The dissipation rates follows predictable trends although there are no experimental findings to date for comparison. The incorporation of metabolic processes and an insight into cell-oriented mechano-sensing processes can help to overcome the limitations involved.Peer ReviewedPostprint (author's final draft

In silico modeling of oxygen-enhanced MRI of specific ventilation.

Specific ventilation imaging (SVI) proposes that using oxygen-enhanced 1H MRI to capture signal change as subjects alternatively breathe room air and 100% O2 provides an estimate of specific ventilation distribution in the lung. How well this technique measures SV and the effect of currently adopted approaches of the technique on resulting SV measurement is open for further exploration. We investigated (1) How well does imaging a single sagittal lung slice represent whole lung SV? (2) What is the influence of pulmonary venous blood on the measured MRI signal and resultant SVI measure? and (3) How does inclusion of misaligned images affect SVI measurement? In this study, we utilized two patient-based in silico models of ventilation, perfusion, and gas exchange to address these questions for normal healthy lungs. Simulation results from the two healthy young subjects show that imaging a single slice is generally representative of whole lung SV distribution, with a calculated SV gradient within 90% of that calculated for whole lung distributions. Contribution of O2 from the venous circulation results in overestimation of SV at a regional level where major pulmonary veins cross the imaging plane, resulting in a 10% increase in SV gradient for the imaging slice. A worst-case scenario simulation of image misalignment increased the SV gradient by 11.4% for the imaged slice