Computational fluid dynamics simulation to evaluate aortic coarctation gradient with contrast-enhanced CT

Coarctation of aorta (CoA) is a narrowing of the aorta leading to a pressure gradient (ΔP) across the coarctation, increased afterload and reduced peripheral perfusion pressures. Indication to invasive treatment is based on values of maximal (systolic) trans-coarctation ΔP. A computational fluid dynamic (CFD) approach is herein presented for the non-invasive haemodynamic assessment of ΔP across CoA. Patient-specific CFD simulations were created from contrast-enhanced computed tomography (CT) and appropriate flow boundary conditions. Computed ΔP was validated with invasive intravascular trans-CoA pressure measurements. Haemodynamic indices, including pressure loss coefficient (PLc), time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI), were also quantified. CFD-estimated ΔP values were comparable to the invasive ones. Moreover, the aorta proximal to CoA was exposed to altered TAWSS and OSI suggesting hypertension. PLc was found as a further geometric marker of CoA severity. Finally, CFD-estimated ΔP confirmed a significant reduction after percutaneous balloon dilatation and stenting of the CoA in one patient (e.g. from ΔP～52 mmHg to ΔP～3 mmHg). The validation of the ΔP computations with catheterisation measurements suggests that CFD simulation, based on CT-derived anatomical data, is a useful tool to readily quantify CoA severity.     

Numerical simulation of haemodynamics of the descending aorta in the non-diabetic and diabetic rabbits

Diabetes mellitus (DM) is a predisposing risk factor leading to macrovascular diseases. Changes in haemodynamics of the diabetic aortas remain largely unclear and relevant computational analyses are lacking in the literature. Ten adult rabbits (1.6–2.2 kg) were collected and the type I diabetic rabbit model was induced by injection of alloxan. A total of five control and five diabetic rabbit aortas were considered for subsequent numerical simulation. The CT scanning was performed to reconstruct three-dimensional model of the individual rabbit descending aorta. The flow velocity waveforms were measured by ultrasound machine and were set to be the inlet boundary conditions. The reconstructed aortas were then imported into ANSYS to perform mesh generation and computational analysis. Results showed that the distributions of haemodynamic indicators time-averaged wall shear stress (TAWSS), oscillating shear index (OSI) and transverse wall shear stress (transWSS) in the non-diabetic rabbit aortas were similar to those in the diabetic rabbit aortas. However, the mean values of TAWSS and transWSS in the non-diabetic rabbit aortas were significantly higher than those values in the diabetic rabbit aortas (TAWSS: p = 0.04; transWSS: p = 0.02). The back of right renal artery tended to have high OSI in both the non-diabetic and the diabetic rabbit aortas. Notably, the regions with high OSI tended to have intense disturbed flow and low TAWSS in the most diabetic rabbit aortas. The results suggest that diabetes leads to changes in haemodynamic parameters in the rabbit aortas. In particular, the lower TAWSS and the higher OSI within the diabetic aortas may further contribute to aortic wall remodeling.     

Numerical simulation of two-phase non-Newtonian blood flow with fluid-structure interaction in aortic dissection

The behavior of blood cells and vessel compliance significantly influence hemodynamic parameters, which are closely related to the development of aortic dissection. Here the two-phase non-Newtonian model and the fluid-structure interaction (FSI) method are coupled to simulate blood flow in a patient-specific dissected aorta. Moreover, three-element Windkessel model is applied to reproduce physiological pressure waves. Important hemodynamic indicators, such as the spatial distribution of red blood cells (RBCs) and vessel wall displacement, which greatly influence the hemodynamic characteristics are analyzed. Results show that the proximal false lumen near the entry tear appears to be a vortex zone with a relatively lower volume fraction of RBCs, a low time-averaged wall shear stress (TAWSS) and a high oscillatory shear index (OSI), providing a suitable physical environment for the formation of atherosclerosis. The highest TAWSS is located in the narrow area of the distal true lumen which might cause further dilation. TAWSS distributions in the FSI model and the rigid wall model show similar trend, while there is a significant difference for the OSI distributions. We suggest that an integrated model is essential to simulate blood flow in a more realistic physiological environment with the ultimate aim of guiding clinical treatment.     

Safety and efficacy of stenting for aortic arch hypoplasia in patients with coarctation of the aorta

Netherlands Heart Journal - Despite a successful repair procedure for coarctation of the aorta (CoA), up to two-thirds of patients remain hypertensive. CoA is often seen in combination...     

Hemodynamic analysis of a novel stent graft design with slit perforations in thoracic aortic aneurysm

Thoracic endovascular aortic repair (TEVAR) has been introduced as a less invasive approach to the treatment of thoracic aortic aneurysm (TAA). However, the effectiveness of TEVAR in the treatment of TAA is often limited due to the complex anatomy of aortic arch. Flow preservation at the three supra-aortic branches further increases the overall technical difficulty. This study proposes a novel stent graft design with slit perforations that can positively alter the hemodynamics at the aortic arch while maintaining blood flow to supra-aortic branches. We carried out a computational fluid dynamic (CFD) analysis to evaluate flow characteristics near stented aortic arch in simplified TAA models, followed by in-vitro experiments using particle image velocimetry (PIV) in a mock circulatory loop. The hemodynamics result was studied in terms of time-averaged wall shear stress (TAWSS), oscillating shear index (OSI), and endothelial cell action potential (ECAP). The results showed that the stent graft with slit perforations can reduce the disturbed flow region considerably. Furthermore, the effect of the slits on flow preservation to the supra-aortic branches was simulated and compared with experimental results. The effectiveness of the stent graft with slit perforations in preserving flow to the branches was demonstrated by both simulated and experimental results. Low TAWSS and elevated ECAP were observed in the aortic arch aneurysm after the placement of the stent graft with slits, implying the potential of thrombus formation in the aneurysm. On the other hand, the effects of the stent grafts with full-slit design and half-slit design on the shear stress did not differ significantly. The present analysis indicated that not only could the stent graft with slit perforations shield the aneurysm from rupture, but also it resulted in a favorable environment for thrombus that can contribute to the shrinkage of the aneurysm.     

Abdominal aortic hemodynamics in young healthy adults at rest and during lower limb exercise: quantification using image-based computer modeling

Localization of atherosclerotic lesions in the abdominal aorta has been previously correlated to areas of adverse hemodynamic conditions, such as flow recirculation, low mean wall shear stress, and high temporal oscillations in shear. Along with its many systemic benefits, exercise is also proposed to have local benefits in the vasculature via the alteration of these regional flow patterns. In this work, subject-specific models of the human abdominal aorta were constructed from magnetic resonance angiograms of five young, healthy subjects, and computer simulations were performed under resting and exercise (50% increase in resting heart rate) pulsatile flow conditions. Velocity fields and spatial variations in mean wall shear stress (WSS) and oscillatory shear index (OSI) are presented. When averaged over all subjects, WSS increased from 4.8 +/- 0.6 to 31.6 +/- 5.7 dyn/cm2 and OSI decreased from 0.22 +/- 0.03 to 0.03 +/- 0.02 in the infrarenal aorta between rest and exercise. WSS significantly increased, whereas OSI decreased between rest and exercise at the supraceliac, infrarenal, and suprabifurcation levels, and significant differences in WSS were found between anterior and posterior sections. These results support the hypothesis that exercise provides localized benefits to the cardiovascular system through acute mechanical stimuli that trigger longer-term biological processes leading to protection against the development or progression of atherosclerosis.     

Diagnostic Value of Transthoracic Echocardiography in Patients with Coarctation of Aorta: The Chinese Experience in 53 Patients Studied between 2008 and 2012 in One Major Medical Center

Although aortography is well known as the “gold standard” for the diagnosis of coarctation of aorta (CoA), the method is invasive, expensive and not readily accepted by some patients. Ultrasound diagnosis for CoA is non-invasive, inexpensive, readily accepted by every patient, and can be repeated as frequently as necessary. The purpose of this presentation is to evaluate the applicability of transthoracic echocardiography for the diagnosis of CoA. The echocardiographic appearances of 53 patients with CoA who had undergone surgery during a 5-year period from January 2008 to October 2012 were analyzed retrospectively, and the results were compared with findings at surgery. Fifty-three patients with CoA include six with isolated CoA and 47 of CoA associated with other cardiac anomalies. Of the 53 operated patients, 48 were correctly diagnosed preoperatively by echocardiography, while two were misdiagnosed as interrupted aortic arch and the diagnosis were missed in three other patients. Thus the diagnostic accuracy rate was 90.6%, and the misdiagnosis rate was 9.4%. Preoperative echocardiographic evaluation offers very satisfactory anatomic assessment in most patients with CoA. It makes preoperative angiography unnecessary. Thus transthoracic echocardiography should be the first-line method for the diagnosis of coarctation of the aorta.     

Realistic non-Newtonian viscosity modelling highlights hemodynamic differences between intracranial aneurysms with and without surface blebs

Most computational fluid dynamic (CFD) simulations of aneurysm hemodynamics assume constant (Newtonian) viscosity, even though blood demonstrates shear-thinning (non-Newtonian) behavior. We sought to evaluate the effect of this simplifying assumption on hemodynamic forces within cerebral aneurysms, especially in regions of low wall shear stress, which are associated with rupture. CFD analysis was performed for both viscosity models using 3D rotational angiography volumes obtained for 26 sidewall aneurysms (12 with blebs, 12 ruptured), and parametric models incorporating blebs at different locations (inflow/outflow zone). Mean and lowest 5% values of time averaged wall shear stress (TAWSS) computed over the dome were compared using Wilcoxon rank-sum test. Newtonian modeling not only resulted in higher aneurysmal TAWSS, specifically in areas of low flow and blebs, but also showed no difference between aneurysms with or without blebs. In contrast, for non-Newtonian analysis, bleb-bearing aneurysms showed significantly lower 5% TAWSS compared to those without (p=0.005), despite no significant difference in mean dome TAWSS (p=0.32). Non-Newtonian modeling also accentuated the differences in dome TAWSS between ruptured and unruptured aneurysms (p<0.001). Parametric models further confirmed that realistic non-Newtonian viscosity resulted in lower bleb TAWSS and higher focal viscosity, especially when located in the outflow zone. The results show that adopting shear-thinning non-Newtonian blood viscosity in CFD simulations of intracranial aneurysms uncovered hemodynamic differences induced by bleb presence on aneurysmal surfaces, and significantly improved discriminant statistics used in risk stratification. These findings underline the possible implications of using a realistic model of blood viscosity in predictive computational hemodynamics.     

An unexpected finding late after repair of coarctation of the aorta

We describe a late complication in a 75-year-old man 50 years after repair of a coarctation of the aorta (CoA). Two years after an aortic valve replacement, mitral valve repair and radiofrequency MAZE the patient presented with dyspnoea and right-sided heart failure, based on a large pseudoaneurysm of the descending aorta, compressing the main bronchus and possibly temporarily the pulmonary arterial system. After sealing the aneurysm with an endovascular stent the patient recovered uneventfully. Recommendations are made for follow-up in patients after repair of CoA. (Neth Heart J 2008;16:260-3.)     

Toward translating near-infrared spectroscopy oxygen saturation data for the non-invasive prediction of spatial and temporal hemodynamics during exercise

Image-based computational fluid dynamics (CFD) studies conducted at rest have shown that atherosclerotic plaque in the thoracic aorta (TA) correlates with adverse wall shear stress (WSS), but there is a paucity of such data under elevated flow conditions. We developed a pedaling exercise protocol to obtain phase contrast magnetic resonance imaging (PC-MRI) blood flow measurements in the TA and brachiocephalic arteries during three-tiered supine pedaling at 130, 150, and 170 % of resting heart rate (HR), and relate these measurements to non-invasive tissue oxygen saturation \((\hbox {StO}_{2})\) acquired by near-infrared spectroscopy (NIRS) while conducting the same protocol. Local quantification of WSS indices by CFD revealed low time-averaged WSS on the outer curvature of the ascending aorta and the inner curvature of the descending aorta (dAo) that progressively increased with exercise, but that remained low on the anterior surface of brachiocephalic arteries. High oscillatory WSS observed on the inner curvature of the aorta persisted during exercise as well. Results suggest locally continuous exposure to potentially deleterious indices of WSS despite benefits of exercise. Linear relationships between flow distributions and tissue oxygen extraction calculated from \(\hbox {StO}_{2}\) were found between the left common carotid versus cerebral tissue \((r^{2}=0.96)\) and the dAo versus leg tissue \((r^{2}=0.87)\). A resulting six-step procedure is presented to use NIRS data as a surrogate for exercise PC-MRI when setting boundary conditions for future CFD studies of the TA under simulated exercise conditions. Relationships and ensemble-averaged PC-MRI inflow waveforms are provided in an online repository for this purpose.     

Fluid Dynamics of Coarctation of the Aorta and Effect of Bicuspid Aortic Valve

Up to 80% of patients with coarctation of the aorta (COA) have a bicuspid aortic valve (BAV). Patients with COA and BAV have elevated risks of aortic complications despite successful surgical repair. The development of such complications involves the interplay between the mechanical forces applied on the artery and the biological processes occurring at the cellular level. The focus of this study is on hemodynamic modifications induced in the aorta in the presence of a COA and a BAV. For this purpose, numerical investigations and magnetic resonance imaging measurements were conducted with different configurations: (1) normal: normal aorta and normal aortic valve; (2) isolated COA: aorta with COA (75% reduction by area) and normal aortic valve; (3) complex COA: aorta with the same severity of COA (75% reduction by area) and BAV. The results show that the coexistence of COA and BAV significantly alters blood flow in the aorta with a significant increase in the maximal velocity, secondary flow, pressure loss, time-averaged wall shear stress and oscillatory shear index downstream of the COA. These findings can contribute to a better understanding of why patients with complex COA have adverse outcome even following a successful surgery.     

Blood pressure during exercise in healthy children

During maximal dynamic exercise the blood pressure (BP) was measured in 497 healthy 9- to 18-year-old children. Systolic BP increased more in the postpubertal groups than in the prepubertal ones. It was also higher in the boys than in the girls of the same age. This was due to a higher work load in boys than girls. Twenty-two subjects had a systolic BP of 200 mmHg or more during the exercise. Only 2 had a resting systolic BP exceeding the mean by 2 standard deviations or more. Three postpubertal boys reached a systolic BP of 240 mmHg at heart rate 170. None had an elevated resting BP. It may be concluded that it cannot be predicted on the basis of the resting BP whether or not an individual is going to have an excessive increase in systolic BP during exercise. The increase in systolic BP to dangerous levels, e.g. 240 mmHg or more, during exercise can only be excluded by means of an individual exercise test.     

Limitations of Doppler echocardiography for the post-operative evaluation of aortic coarctation.

Doppler blood flow measurements and derived pressure differences, through the Bernoulli equation, are used in the diagnosis of aortic coarctation, a congenital stenosis distal to the left subclavian artery. Doppler velocities remain elevated at the coarctation site after successful repair of coarctation, leading to high Doppler derived pressure differences without significant arm-leg pressure differences. We studied this apparent contradiction of two diagnostic methods, in vivo using patient and control data, and in vitro using a hydraulic model. Clinical and echocardiographic data from 31 patients, aged 13.0 +/- 4.0, 10.5 +/- 4.7 yr after coarctectomy by end-to-end anastomosis, and 18 age-matched healthy subjects were reviewed. Doppler peak velocities at the aortic isthmus were elevated in patients (2.2 +/- 0.4 vs. 1.2 +/- 0.2m/s, P < 0.001), corresponding to significant Doppler differences (20 +/- 7 mmHg), however, without significant arm-leg pressure differences. In all patients, a mild anatomic stenosis could still be observed. Local stiffness was increased. The hypothesis that the less distensible surgical scar in post-coarctectomy patients leads to a significant dynamic obstruction in systole was validated in a latex model of the aorta. Rigid rings (0.5-1.5 cm), matching the unloaded aortic diameter, were mounted around the aorta. Under loading conditions, Doppler peak velocities increased by 40 +/-7%, yielding Doppler differences of 21 +/- 3 mmHg, without a significant pressure drop. An alternative expression to calculate pressure differences, using both velocity and geometric information, was validated in the model. In conclusion, post-operatively, Doppler velocities remain elevated due to a mild anatomical and significant dynamic narrowing, but the specific geometry, resembling a tubular hypoplasia rather than an abrupt stenosis, permits an almost complete pressure recovery explaining the occurrence of Doppler differences in disagreement with the negligible arm-leg pressure difference.     

Evaluation of 4D flow MRI-based non-invasive pressure assessment in aortic coarctations

Severity of aortic coarctation (CoA) is currently assessed by estimating trans-coarctation pressure drops through cardiac catheterization or echocardiography. In principle, more detailed information could be obtained non-invasively based on space- and time-resolved magnetic resonance imaging (4D flow) data. Yet the limitations of this imaging technique require testing the accuracy of 4D flow-derived hemodynamic quantities against other methodologies.With the objective of assessing the feasibility and accuracy of this non-invasive method to support the clinical diagnosis of CoA, we developed an algorithm (4DF-FEPPE) to obtain relative pressure distributions from 4D flow data by solving the Poisson pressure equation. 4DF-FEPPE was tested against results from a patient-specific fluid-structure interaction (FSI) simulation, whose patient-specific boundary conditions were prescribed based on 4D flow data. Since numerical simulations provide noise-free pressure fields on fine spatial and temporal scales, our analysis allowed to assess the uncertainties related to 4D flow noise and limited resolution.4DF-FEPPE and FSI results were compared on a series of cross-sections along the aorta. Bland-Altman analysis revealed very good agreement between the two methodologies in terms of instantaneous data at peak systole, end-diastole and time-averaged values: biases (means of differences) were +0.4 mmHg, −1.1 mmHg and +0.6 mmHg, respectively. Limits of agreement (2 SD) were ±0.978 mmHg, ±1.06 mmHg and ±1.97 mmHg, respectively. Peak-to-peak and maximum trans-coarctation pressure drops obtained with 4DF-FEPPE differed from FSI results by 0.75 mmHg and −1.34 mmHg respectively. The present study considers important validation aspects of non-invasive pressure difference estimation based on 4D flow MRI, showing the potential of this technology to be more broadly applied to the clinical practice.     

On the relative importance of rheology for image-based CFD models of the carotid bifurcation

BACKGROUND: Patient-specific computational fluid dynamics (CFD) models derived from medical images often require simplifying assumptions to render the simulations conceptually or computationally tractable. In this study, we investigated the sensitivity of image-based CFD models of the carotid bifurcation to assumptions regarding the blood rheology. METHOD OF APPROACH: CFD simulations of three different patient-specific models were carried out assuming: a reference high-shear Newtonian viscosity, two different non-Newtonian (shear-thinning) rheology models, and Newtonian viscosities based on characteristic shear rates or, equivalently, assumed hematocrits. Sensitivity of wall shear stress (WSS) and oscillatory shear index (OSI) were contextualized with respect to the reproducibility of the reconstructed geometry, and to assumptions regarding the inlet boundary conditions. RESULTS: Sensitivity of WSS to the various rheological assumptions was roughly 1.0 dyn/cm(2) or 8%, nearly seven times less than that due to geometric uncertainty (6.7 dyn/cm(2) or 47%), and on the order of that due to inlet boundary condition assumptions. Similar trends were observed regarding OSI sensitivity. Rescaling the Newtonian viscosity based on time-averaged inlet shear rate served to approximate reasonably, if overestimate slightly, non-Newtonian behavior. CONCLUSIONS: For image-based CFD simulations of the normal carotid bifurcation, the assumption of constant viscosity at a nominal hematocrit is reasonable in light of currently available levels of geometric precision, thus serving to obviate the need to acquire patient-specific rheological data.     

Blood pressure reference intervals for healthy adult chimpanzees (Pan troglodytes)

Background Cardiovascular disease is the primary cause of morbidity and mortality among captive chimpanzees. But there are no clinical definitions of normotension or hypertension in chimpanzees. Methods We analyzed 1 year of blood pressure (BP) data from a population of 261 healthy captive adult chimpanzees using a consistent set of criteria to ascertain health. Results Systolic BP varied by body weight. Diastolic BP varied by age. Median normotension was 126/63 mmHg, with an upper limit of 147/84 mmHg. We defined categories of pre‐hypertension (148/85—153/88 mmHg) and hypertension (≥154/89 mmHg). The prevalence of elevated BP was 15%. The relative risk of mortality was 2.60, compared to normotensive animals. Conclusions We used contemporary methods from human laboratory medicine to define reliable reference intervals for chimpanzee BP. Results allow accurate diagnosis of hypertension and pre‐hypertension, and demonstrate an effect of elevated BP on mortality.     

Familial coarctation of the aorta in three generations.

Four members in three generations of a family were affected by a coarctation of the aorta (CoA) which was mild or severe, either isolated or in association with other cardiac defects. This family suggests that a rare form of CoA could be the result of an autosomal dominant mutation with incomplete penetrance and variable expressivity rather than polygenic inheritance.     

The effect of inlet waveforms on computational hemodynamics of patient-specific intracranial aneurysms

Due to the lack of patient-specific inlet flow waveform measurements, most computational fluid dynamics (CFD) simulations of intracranial aneurysms usually employ waveforms that are not patient-specific as inlet boundary conditions for the computational model. The current study examined how this assumption affects the predicted hemodynamics in patient-specific aneurysm geometries. We examined wall shear stress (WSS) and oscillatory shear index (OSI), the two most widely studied hemodynamic quantities that have been shown to predict aneurysm rupture, as well as maximal WSS (MWSS), energy loss (EL) and pressure loss coefficient (PLc). Sixteen pulsatile CFD simulations were carried out on four typical saccular aneurysms using 4 different waveforms and an identical inflow rate as inlet boundary conditions. Our results demonstrated that under the same mean inflow rate, different waveforms produced almost identical WSS distributions and WSS magnitudes, similar OSI distributions but drastically different OSI magnitudes. The OSI magnitude is correlated with the pulsatility index of the waveform. Furthermore, there is a linear relationship between aneurysm-averaged OSI values calculated from one waveform and those calculated from another waveform. In addition, different waveforms produced similar MWSS, EL and PLc in each aneurysm. In conclusion, inlet waveform has minimal effects on WSS, OSI distribution, MWSS, EL and PLc and a strong effect on OSI magnitude, but aneurysm-averaged OSI from different waveforms has a strong linear correlation with each other across different aneurysms, indicating that for the same aneurysm cohort, different waveforms can consistently stratify (rank) OSI of aneurysms.     

Superficial Femoral Artery Endothelial Responses to a Short-Term Altered Shear Rate Intervention in Healthy Men

In animal and in-vitro models, increased oscillatory shear stress characterized by increased retrograde shear-rate (SR) is associated with acutely decreased endothelial cell function. While previous research suggests a possible detrimental role of elevated retrograde SR on endothelial-function in the brachial artery in humans, little research has been conducted examining arteries in the leg. Examinations of altered shear pattern in the superficial femoral artery (SFA) are important, as this vessel is both prone to atherosclerosis and leg exercise is a common form of activity in humans. Seven healthy men participated; bilateral endothelial-function was assessed via flow-mediated-dilation (FMD) before and after 30-minute unilateral inflations of a thigh blood pressure cuff to either 75 mmHg or 100 mmHg on two separate visits. Inflation of the cuff induced increases in maximum anterograde (p<0.05), maximum retrograde (p<0.01), and oscillatory shear index (OSI) (p<0.001) in the cuffed leg at both inflation pressures. At 100 mmHg the increases in SR were larger in the retrograde than the anterograde direction evidenced by a decrease in mean SR (p<0.01). There was an acute decrease in relative FMD in the cuffed leg alone following inflation to both pressures. These results indicate that in the SFA, altered SR profiles incorporating increased retrograde and OSI influence the attenuation in FMD after a 30-minute unilateral thigh-cuff inflation intervention. Novel information highlighting the importance of OSI calculations and assessments of flow profiles add to current body of knowledge regarding the influence of changes in SR patterns on FMD. Findings from the current study may provide additional insight when designing strategies to combat impaired vascular function in the lower extremity where blood vessels are more prone to atherosclerosis in comparison to the upper extremity.