Hemodynamics of human carotid artery bifurcations: computational studies with models reconstructed from magnetic resonance imaging of normal subjects

Reconstructing blood flow in data-poor regimes: a vasculature network kernel for Gaussian process regression

Blood flow reconstruction in the vasculature is important for many clinical applications. However, in clinical settings, the available data are often quite limited. For instance, transcranial Doppler ultrasound is a non-invasive clinical tool that is commonly used in clinical settings to measure blood velocity waveforms at several locations. This amount of data is grossly insufficient for training machine learning surrogate models, such as deep neural networks or Gaussian process regression. In this work, we propose a Gaussian process regression approach based on empirical kernels constructed by data generated from physics-based simulations-enabling near-real-time reconstruction of blood flow in data-poor regimes. We introduce a novel methodology to reconstruct the kernel within the vascular network. The proposed kernel encodes both spatiotemporal and vessel-to-vessel correlations, thus enabling blood flow reconstruction in vessels that lack direct measurements. We demonstrate that any prediction made with the proposed kernel satisfies the conservation of mass principle. The kernel is constructed by running stochastic one-dimensional blood flow simulations, where the stochasticity captures the epistemic uncertainties, such as lack of knowledge about boundary conditions and uncertainties in vasculature geometries. We demonstrate the performance of the model on three test cases, namely, a simple Y-shaped bifurcation, abdominal aorta and the circle of Willis in the brain.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

Fluid structure interaction analysis to stratify the behavior of different atheromatous carotid plaques

BACKGROUND

In asymptomatic carotid stenosis (ACS), different plaque types, i.e. lipidic (LP), fibrous (FP), and calcific (CP), could have different hemodynamic and structural behaviors.

METHODS

Different carotid plaques, reconstructed from medical imaging of ACS >70%, were analyzed by computing fluid structure interaction (FSI), modeling the spatial distribution of wall shear stresses (WSS), plaque displacements (D), von Mises stresses (VMS), and absorbed elastic energy (AEE) together with their maximum-in-space values at the systole (WSS<inf>syst</inf>, D<inf>syst</inf>, VMS<inf>syst</inf> and AEE<inf>syst</inf>).

RESULTS

WSS resulted significantly higher in CP, whereas D and VMS showed the highest values for LP. Regarding AEE<inf>syst</inf> stored by the plaques, LP absorbed in average 2320 J/m³, FP 408 J/m³ (470%) and CP 99 J/m³ (2240%), (P<0.01, P<0.01, and P<0.01, respectively).

CONCLUSIONS

Depending upon their nature, plaques store different deformations and inner distributions of forces, thus potentially influencing their vulnerability.

Atherogenic potential of microgravity hemodynamics in the carotid bifurcation: a numerical investigation

Long-duration spaceflight poses multiple hazards to human health, including physiological changes associated with microgravity. The hemodynamic adaptations occurring upon entry into weightlessness have been associated with retrograde stagnant flow conditions and thromboembolic events in the venous vasculature but the impact of microgravity on cerebral arterial hemodynamics and function remains poorly understood. The objective of this study was to quantify the effects of microgravity on hemodynamics and wall shear stress (WSS) characteristics in 16 carotid bifurcation geometries reconstructed from ultrasonography images using computational fluid dynamics modeling. Microgravity resulted in a significant 21% increase in flow stasis index, a 22-23% decrease in WSS magnitude and a 16-26% increase in relative residence time in all bifurcation branches, while preserving WSS unidirectionality. In two anatomies, however, microgravity not only promoted flow stasis but also subjected the convex region of the external carotid arterial wall to a moderate increase in WSS bidirectionality, which contrasted with the population average trend. This study suggests that long-term exposure to microgravity has the potential to subject the vasculature to atheroprone hemodynamics and this effect is modulated by subject-specific anatomical features. The exploration of the biological impact of those microgravity-induced WSS aberrations is needed to better define the risk posed by long spaceflights on cardiovascular health.

Characterizing Intracranial Hemodynamics in Sickle Cell Anemia: Impact of Patient-Specific Viscosity

PURPOSE

Pediatric and adult patients with sickle cell anemia (SCA) are at increased risk of stroke and cerebrovascular accident. In the general adult population, there is a relationship between arterial hemodynamics and pathology; however, this relationship in SCA patients remains to be elucidated. The aim of this work was to characterize circle of Willis hemodynamics in patients with SCA and quantify the impact of viscosity choice on pathophysiologically-relevant hemodynamics measures.

METHODS

Based on measured vascular geometries, time-varying flow rates, and blood parameters, detailed patient-specific simulations of the circle of Willis were conducted for SCA patients (n = 6). Simulations quantified the impact of patient-specific and standard blood viscosities on wall shear stress (WSS).

RESULTS

These results demonstrated that use of a standard blood viscosity introduces large errors into the estimation of pathophysiologically-relevant hemodynamic parameters. Standard viscosity models overpredicted peak WSS by 55% and 49% for steady and pulsatile flow, respectively. Moreover, these results demonstrated non-uniform, spatial patterns of positive and negative WSS errors related to viscosity, and standard viscosity simulations overpredicted the time-averaged WSS by 32% (standard deviation = 7.1%). Finally, differences in shear rate demonstrated that the viscosity choice alters the simulated near-wall flow field, impacting hemodynamics measures.

CONCLUSIONS

This work presents simulations of circle of Willis arterial flow in SCA patients and demonstrates the importance and feasibility of using a patient-specific viscosity in these simulations. Accurately characterizing cerebrovascular hemodynamics in SCA populations has potential for elucidating the pathophysiology of large-vessel occlusion, aneurysms, and tissue damage in these patients.

Experimental Validation of Enhanced Magnetic Resonance Imaging (EMRI) Using Particle Image Velocimetry (PIV)

Flow-sensitive four-dimensional Cardiovascular Magnetic Resonance Imaging (4D Flow CMR) has increasingly been utilised to characterise patients' blood flow, in association with patiens' state of health and disease, even though spatial and temporal resolutions still constitute a limit. Computational fluid dynamics (CFD) is a powerful tool that could expand these information and, if integrated with experimentally-obtained velocity fields, would enable to derive a large variety of the flow descriptors of interest. However, the accuracy of the flow parameters is highly influenced by the quality of the input data such as the anatomical model and boundary conditions typically derived from medical images including 4D Flow CMR. We previously proposed a novel approach in which 4D Flow CMR and CFD velocity fields are integrated to obtain an Enhanced 4D Flow CMR (EMRI), allowing to overcome the spatial-resolution limitation of 4D Flow CMR, and enable an accurate quantification of flow. In this paper, the proposed approach is validated in a U bend channel, an idealised model of the human aortic arch. The flow patterns were studied with 4D Flow CMR, CFD and EMRI, and compared with high resolution 2D PIV experiments obtained in pulsatile conditions. The main strengths and limitations of 4D Flow CMR and CFD were illustrated by exploiting the accuracy of PIV by comparing against PIV velocity fields. EMRI flow patterns showed a better qualitative and quantitative agreement with PIV results than the other techniques. EMRI enables to overcome the experimental limitations of MRI-based velocity measurements and the modelling simplifications of CFD, allowing an accurate prediction of complex flow patterns observed experimentally, while satisfying mass and momentum balance equations.

Influence of near-wall PIV data on recirculation hemodynamics in a patient-specific moderate stenosis: Experimental-numerical comparison

BACKGROUND

Recirculation zones within the blood vessels are known to influence the initiation and progression of atherosclerotic lesions. Quantification of recirculation parameters with accuracy remains subjective due to uncertainties in measurement of velocity and derived wall shear stress (WSS).

OBJECTIVE

The primary aim is to determine recirculation height and length from PIV experiments while validating with two different numerical methods: finite-element (FE) and -volume (FV). Secondary aim is to analyze how FE and FV compare within themselves.

METHODS

PIV measurements were performed to obtain velocity profiles at eight cross sections downstream of stenosis at flow rate of 200 ml/min. WSS was obtained by linear/quadratic interpolation of experimental velocity measurements close to wall.

RESULTS

Recirculation length obtained from PIV technique was 1.47 cm and was within 2.2% of previously reported in-vitro measurements. Derived recirculation length from PIV agreed within 6.8% and 8.2% of the FE and FV calculations, respectively. For lower shear rate, linear interpolation with five data points results in least error. For higher shear rate either higher order (quadratic) interpolation with five data points or lower order (linear) with lesser (three) data points leads to better results.

CONCLUSION

Accuracy of the recirculation parameters is dependent on number of near wall PIV data points and the type of interpolation algorithm used.

The structural properties of carotid arteries in carotid artery diseases - a retrospective computed tomography angiography study

Purpose

Aim of the study was to find answers to the following questions: What haemodynamic changes may occur in patients with stenotic, aneurysmal, dissection of the carotid artery and its branches? How do these changes differ in patients with normal and carotid disease?

Material and methods

In order to achieve this aim, the cranio-cervical CT angiography images of patients who were referred to our clinic for any reason and received the diagnosis of carotid stenosis, carotid dissection, and extra or intracranial aneurysm were reviewed retrospectively.

Results

Significant differences were detected in the carotid arteries of normal patients and those with aneurysm and dissection. When normal and aneurismal patients were compared, right and left ICA diameters (p = 0.000, p = 0.002, respectively), total ICA diameters (p = 0.000), carotid left Ø diameters (p = 0.026), right and left total Ø diameters (p = 0.024), and Murray’s and our cosine values of Ø angles (p = 0.001 and p = 0.022, respectively) were found to be different. Also, in a comparison made between normal patients and patients with dissection, right CCA (p = 0.000), ICA (p = 0.001), ECA (p = 0.004) diameters, total CCA (p = 0.001), ICA (p = 0.009), and ECA (p = 0.003) diameters were also found to be different.

Conclusions

This study showed that the presence of aneurysm plays an important role in the remodelling of the carotid arteries. Also, it is understood that Murray’s laws are still valid for the detection of structural deterioration in carotid artery diseases. Hence, it is believed that these data can be used in artificial intelligence studies.

Investigation of the hemodynamic flow conditions and blood-induced stresses inside an abdominal aortic aneurysm by means of a SPH numerical model

The estimation of blood flow-induced loads occurring on the artery wall is affected by uncertainties hidden in the complex interaction of the pulsatile flow, the mechanical parameters of the artery, and the external support conditions. To circumvent these difficulties, a specific tool is developed by combining the aorta displacements measured by an electrocardiogram-gated-computed tomography angiography, with the blood velocity field computed by a smoothed particle hydrodynamics (SPH) numerical model. In the present work, the SPH model has been specifically adapted to the solution of the 3D Navier-Stokes equations inside a domain with boundaries of prescribed motion. Images of the abdominal aorta aneurysm (AAA) of a 44-year-old female patient were acquired during a stabilized cardiac cycle by electrocardiogram-gated-computed tomography angiography. The in vivo kinematic field inside the pulsating arterial wall was estimated by using recent technology, which makes it possible to follow the shape of the arterial wall during a cardiac cycle. We compare the flow conditions and the blood-induced loads, computed by the numerical model under the assumption of a moving arterial wall, with the corresponding results obtained assuming three rigid wall geometries of the vessel during the cardiac cycle. Significant differences were found for the wall shear stress distribution.

Enhancing Magnetic Resonance Imaging With Computational Fluid Dynamics

The analysis of the blood flow in the great thoracic arteries does provide valuable information about the cardiac function and can diagnose the potential development of vascular diseases. Flow-sensitive four-dimensional flow cardiovascular magnetic resonance imaging (4D flow CMR) is often used to characterize patients' blood flow in the clinical environment. Nevertheless, limited spatial and temporal resolution hinders a detailed assessment of the hemodynamics. Computational fluid dynamics (CFD) could expand this information and, integrated with experimental velocity field, enable to derive the pressure maps. However, the limited resolution of the 4D flow CMR and the simplifications of CFD modeling compromise the accuracy of the computed flow parameters. In this article, a novel approach is proposed, where 4D flow CMR and CFD velocity fields are integrated synergistically to obtain an enhanced MR imaging (EMRI). The approach was first tested on a two-dimensional (2D) portion of a pipe, to understand the behavior of the parameters of the model in this novel framework, and afterwards in vivo, to apply it to the analysis of blood flow in a patient-specific human aorta. The outcomes of EMRI are assessed by comparing the computed velocities with the experimental one. The results demonstrate that EMRI preserves flow structures while correcting for experimental noise. Therefore, it can provide better insights into the hemodynamics of cardiovascular problems, overcoming the limitations of MRI and CFD, even when considering a small region of interest. EMRI confirmed its potential to provide more accurate noninvasive estimation of major cardiovascular risk predictors (e.g., flow patterns, endothelial shear stress) and become a novel diagnostic tool.

Ultrasound vector projectile imaging for detection of altered carotid bifurcation hemodynamics during reductions in cardiac output

PURPOSE

Complex blood flow is commonly observed in the carotid bifurcation, although the factors that regulate these patterns beyond arterial geometry are unknown. The emergence of high-frame-rate ultrasound vector flow imaging allows for noninvasive, time-resolved analysis of complex hemodynamic behavior in humans, and it can potentially help researchers understand which physiological stressors can alter carotid bifurcation hemodynamics in vivo. Here, we seek to pursue the first use of vector projectile imaging (VPI), a dynamic form of vector flow imaging, to analyze the regulation of carotid bifurcation hemodynamics during experimental reductions in cardiac output induced via a physiological stressor called lower body negative pressure (LBNP).

METHODS

Seven healthy adults (age: 27 ± 4 yr, 4 men) underwent LBNP at -45 mmHg to simulate a postural hemodynamic response in a controlled environment. Using a research-grade, high-frame-rate ultrasound platform, vector flow estimation in each subject's right carotid bifurcation was performed through a multi-angle plane wave imaging (two transmission angles of 10° and -10°) formulation, and VPI cineloops were generated at a frame rate of 750 fps. Vector concentration was quantified by the resultant blood velocity vector angles within a region of interest; lower concentration indicated greater flow dispersion. Discrete concentration values during peak and late systole were compared across different segments of the carotid artery bifurcation before, and during, LBNP.

RESULTS

Vector projectile imaging revealed that external and internal carotid arteries exhibited regional hemodynamic changes during LBNP, which acted to reduce both the subject's cardiac output (Δ - 1.2 ± 0.5 L/min, -19%; P < 0.01) and peak carotid blood velocity (Δ - 6.30 ± 8.27 cm/s, -7%; P = 0.05). In these carotid artery branches, the vector concentration time trace before and during LBNP were observed to be different. The impact of LBNP on flow complexity in the two carotid artery branches showed variations between subjects.

CONCLUSIONS

Using VPI, intuitive visualization of complex hemodynamic changes can be obtained in healthy humans subjected to LBNP. This imaging tool has potential for further applications in vascular physiology to identify and quantify complex hemodynamic features in humans during different physiological stressor tests that regulate hemodynamics.

A Re-Engineered Software Interface and Workflow for the Open-Source SimVascular Cardiovascular Modeling Package

Patient-specific simulation plays an important role in cardiovascular disease research, diagnosis, surgical planning and medical device design, as well as education in cardiovascular biomechanics. simvascular is an open-source software package encompassing an entire cardiovascular modeling and simulation pipeline from image segmentation, three-dimensional (3D) solid modeling, and mesh generation, to patient-specific simulation and analysis. SimVascular is widely used for cardiovascular basic science and clinical research as well as education, following increased adoption by users and development of a GATEWAY web portal to facilitate educational access. Initial efforts of the project focused on replacing commercial packages with open-source alternatives and adding increased functionality for multiscale modeling, fluid-structure interaction (FSI), and solid modeling operations. In this paper, we introduce a major SimVascular (SV) release that includes a new graphical user interface (GUI) designed to improve user experience. Additional improvements include enhanced data/project management, interactive tools to facilitate user interaction, new boundary condition (BC) functionality, plug-in mechanism to increase modularity, a new 3D segmentation tool, and new computer-aided design (CAD)-based solid modeling capabilities. Here, we focus on major changes to the software platform and outline features added in this new release. We also briefly describe our recent experiences using SimVascular in the classroom for bioengineering education.

Assessment of boundary conditions for CFD simulation in human carotid artery

Computational fluid dynamics (CFD) is an increasingly used method for investigation of hemodynamic parameters and their alterations under pathological conditions, which are important indicators for diagnosis of cardiovascular disease. In hemodynamic simulation models, the employment of appropriate boundary conditions (BCs) determines the computational accuracy of the CFD simulation in comparison with pressure and velocity measurements. In this study, we have first assessed the influence of inlet boundary conditions on hemodynamic CFD simulations. We selected two typical patients suspected of carotid artery disease, with mild stenosis and severe stenosis. Both patients underwent digital subtraction angiography (DSA), magnetic resonance angiography, and the invasive pressure guide wire measured pressure profile. We have performed computational experiments to (1) study the hemodynamic simulation outcomes of distributions of wall shear stress, pressure, pressure gradient and (2) determine the differences in hemodynamic performances caused by inlet BCs derived from DSA and Womersley analytical solution. Our study has found that the difference is related to the severity of the stenosis; the greater the stenosis, the more the difference ensues. Further, in our study, the two typical subjects with invasively measured pressure profile and thirty subjects with ultrasound Doppler velocimeter (UDV) measurement served as the criteria to evaluate the hemodynamic outcomes of wall shear stress, pressure, pressure gradient and velocity due to different outlet BCs based on the Windkessel model, structured-tree model, and fully developed flow model. According to the pressure profiles, the fully developed model appeared to have more fluctuations compared with the other two models. The Windkessel model had more singularities before convergence. The three outlet BCs models also showed good correlation with the UDV measurement, while the Windkessel model appeared to be slightly better ([Formula: see text]). The structured-tree model was seen to have the best performance in terms of available computational cost and accuracy. The results of our numerical simulation and the good correlation with the computed pressure and velocity with their measurements have highlighted the effectiveness of CFD simulation in patient-specific human carotid artery with suspected stenosis.

Impact of Patient-Specific Inflow Velocity Profile on Hemodynamics of the Thoracic Aorta

Computational fluid dynamics (CFD) provides a noninvasive method to functionally assess aortic hemodynamics. The thoracic aorta has an anatomically complex inlet comprising of the aortic valve and root, which is highly prone to different morphologies and pathologies. We investigated the effect of using patient-specific (PS) inflow velocity profiles compared to idealized profiles based on the patient's flow waveform. A healthy 31 yo with a normally functioning tricuspid aortic valve (subject A), and a 52 yo with a bicuspid aortic valve (BAV), aortic valvular stenosis, and dilated ascending aorta (subject B) were studied. Subjects underwent MR angiography to image and reconstruct three-dimensional (3D) geometric models of the thoracic aorta. Flow-magnetic resonance imaging (MRI) was acquired above the aortic valve and used to extract the patient-specific velocity profiles. Subject B's eccentric asymmetrical inflow profile led to highly complex velocity patterns, which were not replicated by the idealized velocity profiles. Despite having identical flow rates, the idealized inflow profiles displayed significantly different peak and radial velocities. Subject A's results showed some similarity between PS and parabolic inflow profiles; however, other parameters such as Flowasymmetry were significantly different. Idealized inflow velocity profiles significantly alter velocity patterns and produce inaccurate hemodynamic assessments in the thoracic aorta. The complex structure of the aortic valve and its predisposition to pathological change means the inflow into the thoracic aorta can be highly variable. CFD analysis of the thoracic aorta needs to utilize fully PS inflow boundary conditions in order to produce truly meaningful results.

Mechanical Characterization of the Vessel Wall by Data Assimilation of Intravascular Ultrasound Studies

Atherosclerotic plaque rupture and erosion are the most important mechanisms underlying the sudden plaque growth, responsible for acute coronary syndromes and even fatal cardiac events. Advances in the understanding of the culprit plaque structure and composition are already reported in the literature, however, there is still much work to be done toward in-vivo plaque visualization and mechanical characterization to assess plaque stability, patient risk, diagnosis and treatment prognosis. In this work, a methodology for the mechanical characterization of the vessel wall plaque and tissues is proposed based on the combination of intravascular ultrasound (IVUS) imaging processing, data assimilation and continuum mechanics models within a high performance computing (HPC) environment. Initially, the IVUS study is gated to obtain volumes of image sequences corresponding to the vessel of interest at different cardiac phases. These sequences are registered against the sequence of the end-diastolic phase to remove transversal and longitudinal rigid motions prescribed by the moving environment due to the heartbeat. Then, optical flow between the image sequences is computed to obtain the displacement fields of the vessel (each associated to a certain pressure level). The obtained displacement fields are regarded as observations within a data assimilation paradigm, which aims to estimate the material parameters of the tissues within the vessel wall. Specifically, a reduced order unscented Kalman filter is employed, endowed with a forward operator which amounts to address the solution of a hyperelastic solid mechanics model in the finite strain regime taking into account the axially stretched state of the vessel, as well as the effect of internal and external forces acting on the arterial wall. Due to the computational burden, a HPC approach is mandatory. Hence, the data assimilation and computational solid mechanics computations are parallelized at three levels: (i) a Kalman filter level; (ii) a cardiac phase level; and (iii) a mesh partitioning level. To illustrate the capabilities of this novel methodology toward the in-vivo analysis of patient-specific vessel constituents, mechanical material parameters are estimated using in-silico and in-vivo data retrieved from IVUS studies. Limitations and potentials of this approach are exposed and discussed.

Does the Ventrica Magnetic Vascular Positioner (MVP®) for Coronary Artery Bypass Grafting Significantly alter Local Fluid Dynamics? a Numeric Study

Objective

Automatic devices have been recently introduced to make the anastomosis procedure quick and efficient when creating a coronary bypass on the beating heart. However, the implantation of these devices could modify the graft configuration, consistently affecting the hemodynamics usually found in the traditional anastomosis. As local fluid dynamics could play a significant role in the onset of vessel wall pathologies, in this article a computational approach was designed to investigate flow patterns in the presence of the Ventrica magnetic vascular positioner (Ventrica MVP®) device.

Methods

A model of standard hand-sewn anastomosis and of automated magnetic anastomosis were constructed, and the finite volume method was used to simulate in silico realistic graft hemodynamics. Synthetic analytical descriptors - i.e., time-averaged wall shear stress (TAWSS), oscillating shear index (OSI) and helical flow index (HFI) - were calculated and compared for quantitative assessment of the anastomosis geometry hemodynamic performance.

Results

In this case study, the same most critical region was identified for the 2 models as the one with the lowest TAWSS and the highest OSI (TAWSS=0.229, OSI=0.255 for the hand-sewn anastomosis; TAWSS=0.297, OSI=0.171 for the Ventrica MVP®). However, the shape of the Ventrica MVP® does not induce more critical wall shear stresses, oscillating flow and damped helicity in the graft fluid dynamics, as compared with conventional anastomosis.

Conclusions

We found that the use of the Ventrica MVP® for the case study under investigation was not associated with more critical fluid dynamics than with conventional hand-sewn anastomosis. Thereby, the device could facilitate beating heart and minimally invasive coronary artery bypass grafting without increasing local hemodynamic-related risks of failure. (Int J Artif Organs 2007; 30: 628–39)

Merging computational fluid dynamics and 4D Flow MRI using proper orthogonal decomposition and ridge regression

Time resolved phase-contrast magnetic resonance imaging 4D-PCMR (also called 4D Flow MRI) data while capable of non-invasively measuring blood velocities, can be affected by acquisition noise, flow artifacts, and resolution limits. In this paper, we present a novel method for merging 4D Flow MRI with computational fluid dynamics (CFD) to address these limitations and to reconstruct de-noised, divergence-free high-resolution flow-fields. Proper orthogonal decomposition (POD) is used to construct the orthonormal basis of the local sampling of the space of all possible solutions to the flow equations both at the low-resolution level of the 4D Flow MRI grid and the high-level resolution of the CFD mesh. Low-resolution, de-noised flow is obtained by projecting in vivo 4D Flow MRI data onto the low-resolution basis vectors. Ridge regression is then used to reconstruct high-resolution de-noised divergence-free solution. The effects of 4D Flow MRI grid resolution, and noise levels on the resulting velocity fields are further investigated. A numerical phantom of the flow through a cerebral aneurysm was used to compare the results obtained using the POD method with those obtained with the state-of-the-art de-noising methods. At the 4D Flow MRI grid resolution, the POD method was shown to preserve the small flow structures better than the other methods, while eliminating noise. Furthermore, the method was shown to successfully reconstruct details at the CFD mesh resolution not discernible at the 4D Flow MRI grid resolution. This method will improve the accuracy of the clinically relevant flow-derived parameters, such as pressure gradients and wall shear stresses, computed from in vivo 4D Flow MRI data.

Investigation of Ultrasound-Measured Flow Velocity, Flow Rate and Wall Shear Rate in Radial and Ulnar Arteries Using Simulation

Parameters of blood flow measured by ultrasound in radial and ulnar arteries, such as flow velocity, flow rate and wall shear rate, are widely used in clinical practice and clinical research. Investigation of these measurements is useful for evaluating accuracy and providing knowledge of error sources. A method for simulating the spectral Doppler ultrasound measurement process was developed with computational fluid dynamics providing flow-field data. Specific scanning factors were adjusted to investigate their influence on estimation of the maximum velocity waveform, and flow rate and wall shear rate were derived using the Womersley equation. The overestimation in maximum velocity increases greatly (peak systolic from about 10% to 30%, time-averaged from about 30% to 50%) when the beam-vessel angle is changed from 30° to 70°. The Womersley equation was able to estimate flow rate in both arteries with less than 3% error, but performed better in the radial artery (2.3% overestimation) than the ulnar artery (15.4% underestimation) in estimating wall shear rate. It is concluded that measurements of flow parameters in the radial and ulnar arteries with clinical ultrasound scanners are prone to clinically significant errors.

Computational medical imaging and hemodynamics framework for functional analysis and assessment of cardiovascular structures

Cardiac dysfunction constitutes common cardiovascular health issues in the society, and has been an investigation topic of strong focus by researchers in the medical imaging community. Diagnostic modalities based on echocardiography, magnetic resonance imaging, chest radiography and computed tomography are common techniques that provide cardiovascular structural information to diagnose heart defects. However, functional information of cardiovascular flow, which can in fact be used to support the diagnosis of many cardiovascular diseases with a myriad of hemodynamics performance indicators, remains unexplored to its full potential. Some of these indicators constitute important cardiac functional parameters affecting the cardiovascular abnormalities. With the advancement of computer technology that facilitates high speed computational fluid dynamics, the realization of a support diagnostic platform of hemodynamics quantification and analysis can be achieved. This article reviews the state-of-the-art medical imaging and high fidelity multi-physics computational analyses that together enable reconstruction of cardiovascular structures and hemodynamic flow patterns within them, such as of the left ventricle (LV) and carotid bifurcations. The combined medical imaging and hemodynamic analysis enables us to study the mechanisms of cardiovascular disease-causing dysfunctions, such as how (1) cardiomyopathy causes left ventricular remodeling and loss of contractility leading to heart failure, and (2) modeling of LV construction and simulation of intra-LV hemodynamics can enable us to determine the optimum procedure of surgical ventriculation to restore its contractility and health This combined medical imaging and hemodynamics framework can potentially extend medical knowledge of cardiovascular defects and associated hemodynamic behavior and their surgical restoration, by means of an integrated medical image diagnostics and hemodynamic performance analysis framework.

Gap-free segmentation of vascular networks with automatic image processing pipeline

Current image processing techniques capture large vessels reliably but often fail to preserve connectivity in bifurcations and small vessels. Imaging artifacts and noise can create gaps and discontinuity of intensity that hinders segmentation of vascular trees. However, topological analysis of vascular trees require proper connectivity without gaps, loops or dangling segments. Proper tree connectivity is also important for high quality rendering of surface meshes for scientific visualization or 3D printing. We present a fully automated vessel enhancement pipeline with automated parameter settings for vessel enhancement of tree-like structures from customary imaging sources, including 3D rotational angiography, magnetic resonance angiography, magnetic resonance venography, and computed tomography angiography. The output of the filter pipeline is a vessel-enhanced image which is ideal for generating anatomical consistent network representations of the cerebral angioarchitecture for further topological or statistical analysis. The filter pipeline combined with computational modeling can potentially improve computer-aided diagnosis of cerebrovascular diseases by delivering biometrics and anatomy of the vasculature. It may serve as the first step in fully automatic epidemiological analysis of large clinical datasets. The automatic analysis would enable rigorous statistical comparison of biometrics in subject-specific vascular trees. The robust and accurate image segmentation using a validated filter pipeline would also eliminate operator dependency that has been observed in manual segmentation. Moreover, manual segmentation is time prohibitive given that vascular trees have more than thousands of segments and bifurcations so that interactive segmentation consumes excessive human resources. Subject-specific trees are a first step toward patient-specific hemodynamic simulations for assessing treatment outcomes.

Automatic recognition of subject-specific cerebrovascular trees

PURPOSE

An image filter designed for reconstructing cerebrovascular trees from MR images is described. Current imaging techniques capture major cerebral vessels reliably, but often fail to detect small vessels, whose contrast is suppressed due to limited resolution, slow blood flow rate, and distortions around bifurcations or nonvascular structures. An incomplete view of angioarchitecture limits the information available to physicians.

METHODS

A novel Hessian-based filter for contrast-enhancement in MR angiography and venography for blood vessel reconstruction without introducing dangling segments is presented. We quantify filter performance with receiver-operating-characteristic and dice-similarity-coefficient analysis. Total extracted vascular length, number-of-segments, volume, surface-to-distance, and positional error are calculated for validation.

RESULTS

Reconstruction of cerebrovascular trees from MR images of six volunteers show that the new filter renders more complete representations of subject-specific cerebrovascular networks. Validation with phantom models shows the filter correctly detects blood vessels across all length scales without failing at bifurcations or distorting diameters.

CONCLUSION

The novel filter can potentially improve the diagnosis of cerebrovascular diseases by delivering metrics and anatomy of the vasculature. It also facilitates the automated analysis of large datasets by computing biometrics free of operator subjectivity. The high quality reconstruction enables computational mesh generation for subject-specific hemodynamic simulations. Magn Reson Med 77:398-410, 2017. © 2016 Wiley Periodicals, Inc.

Functional assessment of thoracic aortic aneurysms - the future of risk prediction?

INTRODUCTION

Treatment guidelines for the thoracic aorta concentrate on size, yet acute aortic dissection or rupture can occur when aortic size is below intervention criteria. Functional imaging and computational techniques are a means of assessing haemodynamic parameters involved in aortic pathology.

SOURCES OF DATA

Original articles, reviews, international guidelines.

AREAS OF AGREEMENT

Computational fluid dynamics and 4D flow MRI allow non-invasive assessment of blood flow parameters and aortic wall biomechanics.

AREAS OF CONTROVERSY

Aortic valve morphology (particularly bicuspid aortic valve) is associated with aneurysm of the ascending aorta, although the exact mechanism of aneurysm formation is not yet established.

GROWING POINTS

Haemodynamic assessment of the thoracic aorta has highlighted parameters which are linked with both clinical outcome and protein changes in the aortic wall. Wall shear stress, flow displacement and helicity are elevated in patients with bicuspid aortic valve, particularly at locations of aneurysm formation.

AREAS TIMELY FOR DEVELOPING RESEARCH

With further validation, functional assessment of the aorta may help identify patients at risk of aortic complications, and introduce new haemodynamic indices into management guidelines.

Patient-specific structural effects on hemodynamics in the ischemic lower limb artery

Lower limb peripheral artery disease is a prevalent chronic non-communicable disease without obvious symptoms. However, the effect of ischemic lower limb peripheral arteries on hemodynamics remains unclear. In this study, we investigated the variation of the hemodynamics caused by patient-specific structural artery characteristics. Computational fluid dynamic simulations were performed on seven lower limb (including superficial femoral, deep femoral and popliteal) artery models that were reconstructed from magnetic resonance imaging. We found that increased wall shear stress (WSS) was mainly caused by the increasing severity of stenosis, bending, and branching. Our results showed that the increase in the WSS value at a stenosis at the bifurcation was 2.7 Pa. In contrast, the isolated stenosis and branch caused a WSS increase of 0.7 Pa and 0.5 Pa, respectively. The WSS in the narrow popliteal artery was more sensitive to a reduction in radius. Our results also demonstrate that the distribution of the velocity and pressure gradient are highly structurally related. At last, Ultrasound Doppler velocimeter measured result was presented as a validation. In conclusion, the distribution of hemodynamics may serve as a supplement for clinical decision-making to prevent the occurrence of a morbid or mortal ischemic event.

SimVascular: An Open Source Pipeline for Cardiovascular Simulation

Patient-specific cardiovascular simulation has become a paradigm in cardiovascular research and is emerging as a powerful tool in basic, translational and clinical research. In this paper we discuss the recent development of a fully open-source SimVascular software package, which provides a complete pipeline from medical image data segmentation to patient-specific blood flow simulation and analysis. This package serves as a research tool for cardiovascular modeling and simulation, and has contributed to numerous advances in personalized medicine, surgical planning and medical device design. The SimVascular software has recently been refactored and expanded to enhance functionality, usability, efficiency and accuracy of image-based patient-specific modeling tools. Moreover, SimVascular previously required several licensed components that hindered new user adoption and code management and our recent developments have replaced these commercial components to create a fully open source pipeline. These developments foster advances in cardiovascular modeling research, increased collaboration, standardization of methods, and a growing developer community.

Modeling the Effect of Red Blood Cells Deformability on Blood Flow Conditions in Human Carotid Artery Bifurcation

The purpose of this work is to predict the effect of impaired red blood cells (RBCs) deformability on blood flow conditions in human carotid artery bifurcation. First, a blood viscosity model is developed that predicts the steady-state blood viscosity as a function of shear rate, plasma viscosity, and mechanical (and geometrical) properties of RBC's. Viscosity model is developed by modifying the well-known Krieger and Dougherty equation for monodisperse suspensions by using the dimensional analysis approach. With the approach, we manage to account for the microscopic properties of RBC's, such as their deformability, in the macroscopic behavior of blood via blood viscosity. In the second part of the paper, the deduced viscosity model is used to numerically predict blood flow conditions in human carotid artery bifurcation. Simulations are performed for different values of RBC's deformability and analyzed by investigating parameters, such as the temporal mean wall shear stress (WSS), oscillatory shear index (OSI), and mean temporal gradient of WSS. The analyses show that the decrease of RBC's deformability decrease the regions of low WSS (i.e., sites known to be prevalent at atherosclerosis-prone regions); increase, in average, the value of WSS along the artery; and decrease the areas of high OSI. These observations provide an insight into the influence of blood's microscopic properties, such as the deformability of RBC's, on hemodynamics in larger arteries and their influence on parameters that are known to play a role in the initiation and progression of atherosclerosis.

Asymmetric pulsation of rat carotid artery bifurcation in three-dimension observed by ultrasound imaging

The arterial structure cyclically fluctuates in three-dimensions (3-D) caused by pulsatile blood flow. The evaluation of arterial wall motion and hemodynamics contributes to early diagnosis of carotid atherosclerosis. Ultrasound is one of the most appropriate imaging modalities to evaluate arterial wall motion in real time. Although many previous studies have discussed the mechanical properties of the carotid artery bifurcation (CAB) from the two-dimensional (2-D) view, the spatio-temporal variation of carotid artery geometry in 3-D has not yet been investigated in detail. In this study, the 3-D data set of CAB from rats was acquired using a high spatio-temporal resolution ultrasound imaging system with a 40 MHz probe using mechanical sector scanning. A total of 31 slices of cross-section images were stored and a spoke scan algorithm was implemented to radially scan the lumen area in polar coordinates based on a pre-tracked seed point. The boundary of the arterial lumen was segmented using intensity-threshold-based boundary detection and fitted by polynomial regression. Two operators, who were trained with the same protocol to minimize inter- and intra-operator variability, manually segmented the lumen boundary on systolic and diastolic phase from the gray-scale images. Finally, the 3-D lumen geometries of CAB during one cardiac cycle were constructed based on the segmented lumen boundaries. From this constructed 3-D geometry, we observed that the CAB geometry favorably expanded to the anterior/posterior direction, parallel to the sagittal plane; and the manually segmented geometry also confirmed the asymmetrical change in bifurcation geometry. This is the first study on visualization and quantification on the asymmetrical variation of the CAB geometry of a rat in 3-D during a whole cardiac cycle. This finding may be useful in understanding hemodynamic etiology of various cardiovascular diseases such as arterial stenosis and its complications, and also provides reference information for numerical simulation studies on arterial wall motion.

Carotid Artery Dynamics during Head Movements: A Reason for Concern with Regard to Carotid Stenting?

Purpose:

To evaluate carotid artery mobility patterns during head movements following carotid angioplasty/stenting (CAS).

Methods:

In 7 patients (all men; mean age 69 years, range 65–76) who had undergone unilateral CAS, 3D time-of-flight magnetic resonance angiography was performed, visualizing both carotid arteries in 5 different head positions (neutral, turned left and right, and bent forward and backward). Maximum intensity projection reconstructions were obtained to measure angulation at the proximal and distal stent junction. Configuration changes of the stented section of the carotid artery and the unstented contralateral artery were judged. Secondly, transverse sections at the level of the carotid bifurcation and at the skull base were used to calculate torsion shear in the common and internal carotid arteries (CCA, ICA) during turned left and right head position. Results were expressed as median (range).

Results:

In neutral head position, maximal angulation at the distal stent junction was 34.3° (32.3°–55.6°). With the head bent forward, this angulation changed to 47.6° (42.6°–85.2°, p=0.028) and when bent backward to 26.5° (25.0°–48.7°, p=0.027). In all patients, configuration changes of the stented sections were absent. The contralateral unstented side showed diffuse configuration changes without specific angulation at one location. With the head turned left and right, the CCA on the stented side was subjected to 28.6° (13.6°–53.7°) and 24.9° (2.0°–50.6°) of torsion shear, respectively. Torsion of the ICA was subsequently 18.1° (12.7°–40.5°) and 15.2° (2.9°–69.4°).

Conclusions:

Following carotid stenting, sharp ICA angulations that are aggravated by forward bending of the head occur at the distal stent junction. The stented section of the carotid artery shows complete lack of flexibility despite highly flexible features of the stents ex vivo. Both the CCA and ICA are subjected to considerable torsion shear with the head turned left and right. This shear is not accommodated by the current stent designs.

Investigation of blood flow in the external carotid artery and its branches with a new 0D peripheral model

BACKGROUND

Patient-specific modelling in clinical studies requires a realistic simulation to be performed within a reasonable computational time. The aim of this study was to develop simple but realistic outflow boundary conditions for patient-specific blood flow simulation which can be used to clarify the distribution of the anticancer agent in intra-arterial chemotherapy for oral cancer.

METHODS

In this study, the boundary conditions are expressed as a zero dimension (0D) resistance model of the peripheral vessel network based on the fractal characteristics of branching arteries combined with knowledge of the circulatory system and the energy minimization principle. This resistance model was applied to four patient-specific blood flow simulations at the region where the common carotid artery bifurcates into the internal and external carotid arteries.

RESULTS

Results of these simulations with the proposed boundary conditions were compared with the results of ultrasound measurements for the same patients. The pressure was found to be within the physiological range. The difference in velocity in the superficial temporal artery results in an error of 5.21 ± 0.78 % between the numerical results and the measurement data.

CONCLUSIONS

The proposed outflow boundary conditions, therefore, constitute a simple resistance-based model and can be used for performing accurate simulations with commercial fluid dynamics software.

Coronary Computed Tomography Angiography Derived Fractional Flow Reserve and Plaque Stress

Fractional flow reserve (FFR) measured during invasive coronary angiography is an independent prognosticator in patients with coronary artery disease and the gold standard for decision making in coronary revascularization. The integration of computational fluid dynamics and quantitative anatomic and physiologic modeling now enables simulation of patient-specific hemodynamic parameters including blood velocity, pressure, pressure gradients, and FFR from standard acquired coronary computed tomography (CT) datasets. In this review article, we describe the potential impact on clinical practice and the science behind noninvasive coronary computed tomography (CT) angiography derived fractional flow reserve (FFR_CT) as well as future applications of this technology in treatment planning and quantifying forces on atherosclerotic plaques.

Computational fluid dynamics simulations of blood flow regularized by 3D phase contrast MRI

Background

Phase contrast magnetic resonance imaging (PC-MRI) is used clinically for quantitative assessment of cardiovascular flow and function, as it is capable of providing directly-measured 3D velocity maps. Alternatively, vascular flow can be estimated from model-based computation fluid dynamics (CFD) calculations. CFD provides arbitrarily high resolution, but its accuracy hinges on model assumptions, while velocity fields measured with PC-MRI generally do not satisfy the equations of fluid dynamics, provide limited resolution, and suffer from partial volume effects. The purpose of this study is to develop a proof-of-concept numerical procedure for constructing a simulated flow field that is influenced by both direct PC-MRI measurements and a fluid physics model, thereby taking advantage of both the accuracy of PC-MRI and the high spatial resolution of CFD. The use of the proposed approach in regularizing 3D flow fields is evaluated.

Methods

The proposed algorithm incorporates both a Newtonian fluid physics model and a linear PC-MRI signal model. The model equations are solved numerically using a modified CFD algorithm. The numerical solution corresponds to the optimal solution of a generalized Tikhonov regularization, which provides a flow field that satisfies the flow physics equations, while being close enough to the measured PC-MRI velocity profile. The feasibility of the proposed approach is demonstrated on data from the carotid bifurcation of one healthy volunteer, and also from a pulsatile carotid flow phantom.

Results

The proposed solver produces flow fields that are in better agreement with direct PC-MRI measurements than CFD alone, and converges faster, while closely satisfying the fluid dynamics equations. For the implementation that provided the best results, the signal-to-error ratio (with respect to the PC-MRI measurements) in the phantom experiment was 6.56 dB higher than that of conventional CFD; in the in vivo experiment, it was 2.15 dB higher.

Conclusions

The proposed approach allows partial or complete measurements to be incorporated into a modified CFD solver, for improving the accuracy of the resulting flow fields estimates. This can be used for reducing scan time, increasing the spatial resolution, and/or denoising the PC-MRI measurements.

Cardiovascular magnetic resonance phase contrast imaging

Cardiovascular magnetic resonance (CMR) phase contrast imaging has undergone a wide range of changes with the development and availability of improved calibration procedures, visualization tools, and analysis methods. This article provides a comprehensive review of the current state-of-the-art in CMR phase contrast imaging methodology, clinical applications including summaries of past clinical performance, and emerging research and clinical applications that utilize today's latest technology.

Determination of vessel cross-sectional area by thresholding in Radon space

The cross-sectional area of a blood vessel determines its resistance, and thus is a regulator of local blood flow. However, the cross-sections of penetrating vessels in the cortex can be non-circular, and dilation and constriction can change the shape of the vessels. We show that observed vessel shape changes can introduce large errors in flux calculations when using a single diameter measurement. Because of these shape changes, typical diameter measurement approaches, such as the full-width at half-maximum (FWHM) that depend on a single diameter axis will generate erroneous results, especially when calculating flux. Here, we present an automated method--thresholding in Radon space (TiRS)--for determining the cross-sectional area of a convex object, such as a penetrating vessel observed with two-photon laser scanning microscopy (2PLSM). The thresholded image is transformed back to image space and contiguous pixels are segmented. The TiRS method is analogous to taking the FWHM across multiple axes and is more robust to noise and shape changes than FWHM and thresholding methods. We demonstrate the superior precision of the TiRS method with in vivo 2PLSM measurements of vessel diameter.

Turbulence intensity measurements using particle image velocimetry in diseased carotid artery models: effect of stenosis severity, plaque eccentricity, and ulceration

Clinical decision-making for the treatment of patients with diseased carotid artery is mainly based on the severity of the stenosis. However, stenosis severity alone is not a sensitive indicator, and other local factors for the assessment of stroke risk are required. Flow disturbance is of particular interest due to its proven association with increased thromboembolic activities. The objective of this study was to investigate the level of turbulence intensity (TI) with regards to certain geometrical features of the plaque - namely stenosis severity, eccentricity, and ulceration. A family of eight carotid-artery bifurcation models was examined using particle image velocimetry. Results showed a marked difference in turbulence intensity among these models; increasing degree of stenosis severity resulted in increased turbulence intensity, going from 0.12 m/s for mild stenosis to 0.37 m/s for severe stenosis (with concentric geometry). Moreover, independent of stenosis severity, eccentricity led to further elevations in turbulence intensity, increasing TI by 0.05-0.10 m/s over the counterpart concentric plaque. The presence of ulceration (in a 50% eccentric plaque) produced a larger portion of moderate turbulence intensity (~0.10 m/s) compared to the non-ulcerated model, more proximal to the bifurcation apex in the post-stenotic recirculation zone. The effect of plaque eccentricity and ulceration in enhancing the downstream turbulence has potential clinical implications for a more sensitive assessment of stroke risk beyond stenosis severity alone.

Effects of carotid artery stenting on arterial geometry

BACKGROUND

The role of carotid artery stenting (CAS) for the treatment of carotid artery disease continues to evolve, despite higher stroke and restenosis risks for CAS compared with conventional open endarterectomy. Understanding the effects of CAS on arterial geometry, which strongly influence hemodynamics and wall mechanics, can assist in better stratifying the inherent risk of CAS to individual patients.

STUDY DESIGN

Fifteen consecutive patients undergoing CAS had pre- and post-stenting CT angiograms. These images were used to reconstruct the 3-dimensional geometries of the bilateral carotid arteries from their origin to the skull base. Quantitative assessment of the carotid bifurcation angle, cross-sectional area, tortuosity and artery length, were compared pre- and post-stenting. Plaque volume and calcification were also measured. Mathematical models were devised to determine the mechanisms of CAS-induced geometric changes, and their mechanical and hemodynamic significances.

RESULTS

Major and moderate changes in arterial tortuosity and elongation were seen in 5 (33%) patients. Characteristics most associated with the development of CAS-induced geometric changes were stenoses located in the internal carotid artery distal to the carotid bulb, circumferential distribution of plaque, and plaque calcification. Modeling did not demonstrate substantial alterations in wall shear stress due to geometric changes, but did show considerable increases in arterial wall axial stress.

CONCLUSIONS

Carotid artery stenting can produce geometric changes to the artery that promote favorable conditions for complications and recurrent disease. Patients with circumferential, highly calcified plaques that are located relatively distal in the internal carotid artery are most likely to have post-stenting geometric changes.

Understanding the role of hemodynamics in the initiation, progression, rupture, and treatment outcome of cerebral aneurysm from medical image-based computational studies

About a decade ago, the first image-based computational hemodynamic studies of cerebral aneurysms were presented. Their potential for clinical applications was the result of a right combination of medical image processing, vascular reconstruction, and grid generation techniques used to reconstruct personalized domains for computational fluid and solid dynamics solvers and data analysis and visualization techniques. A considerable number of studies have captivated the attention of clinicians, neurosurgeons, and neuroradiologists, who realized the ability of those tools to help in understanding the role played by hemodynamics in the natural history and management of intracranial aneurysms. This paper intends to summarize the most relevant results in the field reported during the last years.

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.

Three-dimensional geometry of the human carotid artery

Accurate characterization of carotid artery geometry is vital to our understanding of the pathogenesis of atherosclerosis. Three-dimensional computer reconstructions based on medical imaging are now ubiquitous; however, mean carotid artery geometry has not yet been comprehensively characterized. The goal of this work was to build and study such geometry based on data from 16 male patients with severe carotid artery disease. Results of computerized tomography angiography were used to analyze the cross-sectional images implementing a semiautomated segmentation algorithm. Extracted data were used to reconstruct the mean three-dimensional geometry and to determine average values and variability of bifurcation and planarity angles, diameters and cross-sectional areas. Contrary to simplified carotid geometry typically depicted and used, our mean artery was tortuous exhibiting nonplanarity and complex curvature and torsion variations. The bifurcation angle was 36 deg ± 11 deg if measured using arterial centerlines and 15 deg ± 14 deg if measured between the walls of the carotid bifurcation branches. The average planarity angle was 11 deg ± 10 deg. Both bifurcation and planarity angles were substantially smaller than values reported in most studies. Cross sections were elliptical, with an average ratio of semimajor to semiminor axes of 1.2. The cross-sectional area increased twofold in the bulb compared to the proximal common, but then decreased 1.5-fold for the combined area of distal internal and external carotid artery. Inter-patient variability was substantial, especially in the bulb region; however, some common geometrical features were observed in most patients. Obtained quantitative data on the mean carotid artery geometry and its variability among patients with severe carotid artery disease can be used by biomedical engineers and biomechanics vascular modelers in their studies of carotid pathophysiology, and by endovascular device and materials manufacturers interested in the mean geometrical features of the artery to target the broad patient population.

Vascular remodeling in autogenous arterio-venous fistulas by MRI and CFD

Hemodynamic parameters play an important role in regulating vascular remodeling in arterio-venous fistula (AVF) maturation. Investigating the changes in hemodynamic parameters during AVF maturation is expected to improve our understanding of fistula failure, but very little data on actual temporal changes in human AVFs is available. The present study aimed to assess the feasibility of using a noncontrast-enhanced MRI protocol combined with CFD modeling to relate hemodynamic changes to vascular remodeling following native AVF placement. MR angiography (MRA) and MR velocimetry (MRV) data was acquired peri-operatively, 1 month, and 3 months later in three patients. Vascular geometries were obtained by segmentation of the MRA images. Pulsatile flow simulations were performed in the patient specific vascular geometries with time-dependent boundary conditions prescribed from MRV measurements. A principal result of the study is the description of WSS changes over time in the same patients. The disturbed flow observed in the venous segments resulted in a variability of the WSS distribution and could be responsible for the non-uniform remodeling of the vessel. The artery did not show regions of disturbed flow upstream from the anastomosis, which would be consistent with the uniform remodeling. MRI use demonstrated the ability to provide a comprehensive evaluation of clinically relevant information for the investigation of upper extremity AVFs. 3D geometry from MRA in combination with MRV provides the opportunity to perform detailed CFD analysis of local hemodynamics in order to determine flow descriptors affecting fistula maturation.

THE INFLUENCE OF WALL COMPLIANCE ON FLOW PATTERN IN A CURVED ARTERY EXPOSED TO A DYNAMIC PHYSIOLOGICAL ENVIRONMENT: AN ELASTIC WALL MODEL VERSUS A RIGID WALL MODEL

Plenty of well-established medical research works have shown that many vascular diseases such as stenosis and atherosclerosis are prone to appear in curved arteries. In this paper, we investigated the influence of wall compliance on flow pattern in curved arteries exposed to dynamic physiological environments in order to understand the hemodynamic mechanism and provide a basis for clinical research in related areas. Representative curved arteries with elastic and rigid walls are constructed by computers. The fluid-structure interaction (FSI) effect is considered in our calculations. Physiological velocity profile is assigned as the inlet boundary condition. No-slip boundary condition is applied at the blood-wall interface. Our results show that the maximum axial velocity in the rigid wall model is larger than that in the elastic wall model. Wall compliance also has a remarkable effect on backflow patterns. Significant differences in pressure distribution are found between the elastic and rigid wall models. Blood strain rate distribution patterns in the two models were also compared. It was interesting to discover that in the straight part of the artery, the flexible wall made the counter-rotating vortices induced by the curvature gradually disappear along a downstream direction. However, for the flow feature in the rigid wall model, strong vortices existed throughout the entire straight part of the artery. It revealed that the increment of wall rigidity results in a reduction in wall movement capacity, thus affecting the physiological function of the arterial wall, making it incapable of effectively regulating the flow pattern inside the artery. The current work indicates that the influence of wall compliance on flow pattern in curved artery is significant.

Alterations in Pulse Pressure Affect Artery Function

Pulse pressure changes in response to cardiovascular diseases and interventions, but its effect on vascular wall structure and function is poorly understood. We examined the effect of increased or decreased pulse pressure on artery function, cellular function, and extracellular matrix remodeling. Porcine carotid arteries were cultured under non-pulsatile (100 mmHg), pulsatile (70-130 mmHg), or hyper-pulsatile pressure (50-150 mmHg) for 1 to 3 days. Vasomotor response, wall permeability, cell proliferation, apoptosis, extracellular matrix remodeling, and proteins involved in atherogenesis were examined. Our results showed that hyper-pulsatile pressure decreased the artery response to sodium nitroprusside, basal tone, and wall permeability after three days. Non-pulsatile pressure increased cell proliferation. Neither hyper-pulsatile nor non-pulsatile pressure caused a change in the extracellular matrix or in the expression of matrix metalloproteinase-2 (MMP-2), MMP-9, caveolin-1, or α-actin. Hyper-pulsatile pressure increased monocyte chemotactic protein-1 gene expression. Taken together, these changes indicate that pulse pressure has a limited effect on the artery immediately after its application. Specifically an increase in pulse pressure alters the artery tone and wall permeability while a decrease in pulse pressure alters cell proliferation. Overall these results provide insight into how the artery initially responds to changes in pulse pressure.

Model reduction techniques for fast blood flow simulation in parametrized geometries

In this paper, we propose a new model reduction technique aimed at real-time blood flow simulations on a given family of geometrical shapes of arterial vessels. Our approach is based on the combination of a low-dimensional shape parametrization of the computational domain and the reduced basis method to solve the associated parametrized flow equations. We propose a preliminary analysis carried on a set of arterial vessel geometries, described by means of a radial basis functions parametrization. In order to account for patient-specific arterial configurations, we reconstruct the latter by solving a suitable parameter identification problem. Real-time simulation of blood flows are thus performed on each reconstructed parametrized geometry, by means of the reduced basis method. We focus on a family of parametrized carotid artery bifurcations, by modelling blood flows using Navier-Stokes equations and measuring distributed outputs such as viscous energy dissipation or vorticity. The latter are indexes that might be correlated with the assessment of pathological risks. The approach advocated here can be applied to a broad variety of (different) flow problems related with geometry/shape variation, for instance related with shape sensitivity analysis, parametric exploration and shape design.