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

Computational fluid dynamicaccuracy in mimicking changes in blood hemodynamics in patients with acute type IIIb aortic dissection treated with TEVAR

Background: We aimed to verify the accuracy of the Computational Fluid Dynamics (CFD) algorithm for blood flow reconstruction for type IIIb aortic dissection (TBAD) before and after thoracic endovascular aortic repair (TEVAR). Methods: We made 3D models of the aorta and its branches using pre- and post-operative CT data from five patients treated for TBAD. The CFD technique was used to quantify the displacement forces acting on the aortic wall in the areas of endograft, mass flow rate/velocity and wall shear stress (WSS). Calculated results were verified with ultrasonography (USG-Doppler) data. Results: CFD results indicated that the TEVAR procedure caused a 7-fold improvement in overall blood flow through the aorta (p = 0.0001), which is in line with USG-Doppler data. A comparison of CFD results and USG-Doppler data indicated no significant change in blood flow through the analysed arteries. CFD also showed a significant increase in flow rate for thoracic trunk and renal arteries, which was in accordance with USG-Doppler data (accuracy 90% and 99.9%). Moreover, we observed a significant decrease in WSS values within the whole aorta after TEVAR compared to pre-TEVAR (1.34 ± 0.20 Pa vs. 3.80 ± 0.59 Pa, respectively, p = 0.0001). This decrease was shown by a significant reduction in WSS and WSS contours in the thoracic aorta (from 3.10 ± 0.27 Pa to 1.34 ± 0.11Pa, p = 0.043) and renal arteries (from 4.40 ± 0.25 Pa to 1.50 ± 0.22 Pa p = 0.043). Conclusions: Post-operative remodelling of the aorta after TEVAR for TBAD improved hemodynamic patterns reflected by flow, velocity and WSS with an accuracy of 99%

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

Estimation of wall shear stress using 4D flow cardiovascular MRI and computational fluid dynamics

Electronic version of an article published as Journal of mechanics in medicine and biology, 0, 1750046 (2016), 16 pages. DOI:10.1142/S0219519417500464 © World Scientific Publishing CompanyIn the last few years, wall shear stress (WSS) has arisen as a new diagnostic indicator in patients with arterial disease. There is a substantial evidence that the WSS plays a significant role, together with hemodynamic indicators, in initiation and progression of the vascular diseases. Estimation of WSS values, therefore, may be of clinical significance and the methods employed for its measurement are crucial for clinical community. Recently, four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has been widely used in a number of applications for visualization and quantification of blood flow, and although the sensitivity to blood flow measurement has increased, it is not yet able to provide an accurate three-dimensional (3D) WSS distribution. The aim of this work is to evaluate the aortic blood flow features and the associated WSS by the combination of 4D flow cardiovascular magnetic resonance (4D CMR) and computational fluid dynamics technique. In particular, in this work, we used the 4D CMR to obtain the spatial domain and the boundary conditions needed to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. Similar WSS distributions were found for cases simulated. A sensitivity analysis was done to check the accuracy of the method. 4D CMR begins to be a reliable tool to estimate the WSS within the entire thoracic aorta using computational fluid dynamics. The combination of both techniques may provide the ideal tool to help tackle these and other problems related to wall shear estimation.Peer ReviewedPostprint (author's final draft

The influence of anesthesia and fluid-structure interaction on simulated shear stress patterns in the carotid bifurcation of mice

Background: Low and oscillatory wall shear stresses (WSS) near aortic bifurcations have been linked to the onset of atherosclerosis. In previous work, we calculated detailed WSS patterns in the carotid bifurcation of mice using a Fluid-structure interaction (FSI) approach. We subsequently fed the animals a high-fat diet and linked the results of the FSI simulations to those of atherosclerotic plaque location on a within-subject basis. However, these simulations were based on boundary conditions measured under anesthesia, while active mice might experience different hemodynamics. Moreover, the FSI technique for mouse-specific simulations is both time- and labor-intensive, and might be replaced by simpler and easier Computational Fluid Dynamics (CFD) simulations. The goal of the current work was (i) to compare WSS patterns based on anesthesia conditions to those representing active resting and exercising conditions; and (ii) to compare WSS patterns based on FSI simulations to those based on steady-state and transient CFD simulations. Methods: For each of the 3 computational techniques (steady state CFD, transient CFD, FSI) we performed 5 simulations: 1 for anesthesia, 2 for conscious resting conditions and 2 more for conscious active conditions. The inflow, pressure and heart rate were scaled according to representative in vivo measurements obtained from literature. Results: When normalized by the maximal shear stress value, shear stress patterns were similar for the 3 computational techniques. For all activity levels, steady state CFD led to an overestimation of WSS values, while FSI simulations yielded a clear increase in WSS reversal at the outer side of the sinus of the external carotid artery that was not visible in transient CFD-simulations. Furthermore, the FSI simulations in the highest locomotor activity state showed a flow recirculation zone in the external carotid artery that was not present under anesthesia. This recirculation went hand in hand with locally increased WSS reversal. Conclusions: Our data show that FSI simulations are not necessary to obtain normalized WSS patterns, but indispensable to assess the oscillatory behavior of the WSS in mice. Flow recirculation and WSS reversal at the external carotid artery may occur during high locomotor activity while they are not present under anesthesia. These phenomena might thus influence plaque formation to a larger extent than what was previously assumed. (C) 2016 Elsevier Ltd. All rights reserved

Computational fluid dynamic study on effect of Carreau-Yasuda and Newtonian blood viscosity models on hemodynamic parameters

Pulsatile blood flow through the human carotid artery is studied using Computational Fluid Dynamics (CFD) in order to investigate the effect of blood rheology on the hemodynamic parameters. The carotid artery model used is segmented and reconstructed from the Magnetic Resonance Images (MRI) of a specific patient. The results of a non-Newtonian (Carreau-Yasuda) model and a Newtonian model are studied and compared. The results are represented for each peak systole where it is observed that there is significant variation in the spatial parameters between the two models considered in the study. Comparison of local shear stress magnitude in different branches namely Common Carotid Artery (CCA), Internal Carotid Artery (ICA) and External Carotid Artery (ECA) show that the shear thinning property of blood influences the Wall Shear Stress (WSS) variation. This is observed in branches where there is reduction in diameter and where the diameter reduces due to plaque deposition and also in the region where there is flow recirculation like carotid sinus

Pulsatile spiral blood flow through arterial stenosis

Pulsatile spiral blood flow in a modelled three-dimensional arterial stenosis, with a 75% cross-sectional area reduction, is investigated by using numerical fluid dynamics. Two-equation k-ω model is used for the simulation of the transitional flow with Reynolds numbers 500 and 1000. It is found that the spiral component increases the static pressure in the vessel during the deceleration phase of the flow pulse. In addition, the spiral component reduces the turbulence intensity and wall shear stress found in the post-stenosis region of the vessel in the early stages of the flow pulse. Hence, the findings agree with the results of Stonebridge et al. (2004). In addition, the results of the effects of a spiral component on time-varying flow are presented and discussed along with the relevant pathological issues

Multiscale Fluid-Structure Interaction Models Development and Applications to the 3D Elements of a Human Cardiovascular System

Cardiovascular diseases (CVD) are the number one cause of death of humans in the United States and worldwide. Accurate, non-invasive, and cheaper diagnosis methods have always been on demand as cardiovascular monitoring increase in prevalence. The primary causes of the various forms of these CVDs are atherosclerosis and aneurysms in the blood vessels. Current noninvasive methods (i.e., statistical/medical) permit fairly accurate detection of the disease once clinical symptoms are suggestive of the existence of hemodynamic disorders. Therefore, the recent surge of hemodynamics models facilitated the prediction of cardiovascular conditions. The hemodynamic modeling of a human circulatory system involves varying levels of complexity which must be accounted for and resolved. Pulse-wave propagation effects and high aspect-ratio segments of the vasculature are represented using a quasi-one-dimensional (1D), non-steady, averaged over the cross-section models. However, these reduced 1D models do not account for the blood flow patterns (recirculation zones), vessel wall shear stresses and quantification of repetitive mechanical stresses which helps to predict a vessel life. Even a whole three-dimensional (3D) modeling of the vasculature is computationally intensive and do not fit the timeline of practical use. Thus the intertwining of a quasi 1D global vasculature model with a specific/risk-prone 3D local vessel ones is imperative. This research forms part of a multiphysics project that aims to improve the detailed understanding of the hemodynamics by investigating a computational model of fluid-structure interaction (FSI) of in vivo blood flow. First idealized computational a 3D FSI artery model is configured and executed in ANSYS Workbench, forming an implicit coupling of the blood flow and vessel walls. Then the thesis focuses on an approach developed to employ commercial tools rather than in-house mathematical models in achieving multiscale simulations. A robust algorithm is constructed to combine stabilization techniques to simultaneously overcome the added-mass effect in 3D FSI simulation and mathematical difficulties such as the assignment of boundary conditions at the interface between the 3D-1D coupling. Applications can be of numerical examples evaluating the change of hemodynamic parameters and diagnosis of an abdominal aneurysm, deep vein thrombosis, and bifurcation areas

Multiscale Modeling of Hemodynamics in Human Vessel Network and Its Applications in Cerebral Aneurysms

Three-dimensional (3D) simulation of patient-specific morphological models has been widely used to provide the hemodynamic information of individual patients, such as wall shear stress (WSS), oscillatory shear index (OSI), and flow patterns, etc. Since patient-specific morphological segment was only restricted locally, boundary conditions (BCs) are required to implement the CFD simulation. Direct measurements of the flow and pressure waveforms were often required as input BCs for 3D CFD simulations of patient-specific models. However, as the morphology develops, the feedback from this topological deformation may lead to BCs being altered, and hence without this feedback, the flow characteristics of the morphology are only computed locally. A one-dimensional (1D) numerical model containing the entire human vessel network has been proposed to compute the global hemodynamics. In the meantime, experimental studies of blood flow in the patient-specific modeling of the circle of Willies (CoW) was conducted. The flow and pressure waveforms were quantified to validate the accuracy of the pure 1D model. This 1D model will be coupled with a 3D morphological model to account for the effects of the altered BCs. The proposed 1D-3D multi-scale modeling approach investigates how the global hemodynamic changes can be induced by the local morphological effects, and in consequence, may further result in altering of BCs to interfere with the solution of the 3D simulation. Validation of the proposed multi-scale model has also been made by comparing the solution of the flow rate and pressure waveforms with the experimental data and 3D numerical simulations reported in the literature. Moreover, the multi-scale model is extended to study a patient-specific cerebral aneurysm and a stenosis model. The proposed multi-scale model can be used as an alternative to current approaches to study intracranial vascular diseases such as an aneurysm, stenosis, and combined cases