Thrombotic risk stratification using computational modeling in patients with coronary artery aneurysms following Kawasaki disease

Kawasaki disease (KD) is the leading cause of acquired heart disease in children and can result in life-threatening coronary artery aneurysms in up to 25 % of patients. These aneurysms put patients at risk of thrombus formation, myocardial infarction, and sudden death. Clinicians must therefore decide which patients should be treated with anticoagulant medication, and/or surgical or percutaneous intervention. Current recommendations regarding initiation of anticoagulant therapy are based on anatomy alone with historical data suggesting that patients with aneurysms \(\ge \) 8 mm are at greatest risk of thrombosis. Given the multitude of variables that influence thrombus formation, we postulated that hemodynamic data derived from patient-specific simulations would more accurately predict risk of thrombosis than maximum diameter alone. Patient-specific blood flow simulations were performed on five KD patients with aneurysms and one KD patient with normal coronary arteries. Key hemodynamic and geometric parameters, including wall shear stress, particle residence time, and shape indices, were extracted from the models and simulations and compared with clinical outcomes. Preliminary fluid structure interaction simulations with radial expansion were performed, revealing modest differences in wall shear stress compared to the rigid wall case. Simulations provide compelling evidence that hemodynamic parameters may be a more accurate predictor of thrombotic risk than aneurysm diameter alone and motivate the need for follow-up studies with a larger cohort. These results suggest that a clinical index incorporating hemodynamic information be used in the future to select patients for anticoagulant therapy.     

Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were 'treated' with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were ‘treated’ with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

Mechanical instability of normal and aneurysmal arteries

Tortuous arteries associated with aneurysms have been observed in aged patients with atherosclerosis and hypertension. However, the underlying mechanism is poorly understood. The objective of this study was to determine the effect of aneurysms on arterial buckling instability and the effect of buckling on aneurysm wall stress. We investigated the mechanical buckling and post-buckling behavior of normal and aneurysmal carotid arteries and aorta’s using computational simulations and experimental measurements to elucidate the interrelationship between artery buckling and aneurysms. Buckling tests were done in porcine carotid arteries with small aneurysms created using elastase treatment. Parametric studies were done for model aneurysms with orthotropic nonlinear elastic walls using finite element simulations. Our results demonstrated that arteries buckled at a critical buckling pressure and the post-buckling deflection increased nonlinearly with increasing pressure. The presence of an aneurysm can reduce the critical buckling pressure of arteries, although the effect depends on the aneurysm’s dimensions. Buckled aneurysms demonstrated a higher peak wall stress compared to unbuckled aneurysms under the same lumen pressure. We conclude that aneurysmal arteries are vulnerable to mechanical buckling and mechanical buckling could lead to high stresses in the aneurysm wall. Buckling could be a possible mechanism for the development of tortuous aneurysmal arteries such as in the Loeys–Dietz syndrome.     

Fluid-structure interaction modeling of aneurysmal arteries under steady-state and pulsatile blood flow: a stability analysis

Tortuous aneurysmal arteries are often associated with a higher risk of rupture but the mechanism remains unclear. The goal of this study was to analyze the buckling and post-buckling behaviors of aneurysmal arteries under pulsatile flow. To accomplish this goal, we analyzed the buckling behavior of model carotid and abdominal aorta with aneurysms by utilizing fluid-structure interaction (FSI) method with realistic waveforms boundary conditions. FSI simulations were done under steady-state and pulsatile flow for normal (1.5) and reduced (1.3) axial stretch ratios to investigate the influence of aneurysm, pulsatile lumen pressure and axial tension on stability. Our results indicated that aneurysmal artery buckled at the critical buckling pressure and its deflection nonlinearly increased with increasing lumen pressure. Buckling elevates the peak stress (up to 118%). The maximum aneurysm wall stress at pulsatile FSI flow was (29%) higher than under static pressure at the peak lumen pressure of 130 mmHg. Buckling results show an increase in lumen shear stress at the inner side of the maximum deflection. Vortex flow was dramatically enlarged with increasing lumen pressure and artery diameter. Aneurysmal arteries are more susceptible than normal arteries to mechanical instability which causes high stresses in the aneurysm wall that could lead to aneurysm rupture.     

The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions

In the abdominal segment of the human aorta under a patient's average resting conditions, pulsatile blood flow exhibits complex laminar patterns with secondary flows induced by adjacent branches and irregular vessel geometries. The flow dynamics becomes more complex when there is a pathological condition that causes changes in the normal structural composition of the vessel wall, for example, in the presence of an aneurysm. This work examines the hemodynamics of pulsatile blood flow in hypothetical three-dimensional models of abdominal aortic aneurysms (AAAs). Numerical predictions of blood flow patterns and hemodynamic stresses in AAAs are performed in single-aneurysm, asymmetric, rigid wall models using the finite element method. We characterize pulsatile flow dynamics in AAAs for average resting conditions by means of identifying regions of disturbed flow and quantifying the disturbance by evaluating flow-induced stresses at the aneurysm wall, specifically wall pressure and wall shear stress. Physiologically realistic abdominal aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 < or = Rem < or = 300, corresponding to a range of peak Reynolds numbers 262.5 < or = Repeak < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is depicted by a sequence of four different flow phases in one period of the cardiac pulse. Peak wall shear stress and peak wall pressure are reported as a function of the time-average Reynolds number and aneurysm asymmetry. The effect of asymmetry in hypothetically shaped AAAs is to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole.     

Numerical simulation of saccular aneurysm hemodynamics: influence of morphology on rupture risk

The governing equations for pulsatile fluid flow were solved in their finite volume formulation in order to simulate blood flow in a variety of three-dimensional aneurysm geometries. The influence of geometric factors on flow patterns and fluid mechanical forces was studied with the goal of identifying the risk of aneurysm rupture. Aneurysm morphology was characterized by quantitative shape indices reflecting the three dimensionality of the vasculature derived from clinical studies. Recirculation zones and secondary flows were observed in aneurysms and arteries. Regions of extreme and alternating shear stress were observed and identified as sites for potential aneurysm rupture. The ellipticity of an aneurysm was observed to be strongly correlated with wall shear stress at the aneurysm fundus, while its non-sphericity, volume, and degree of undulation were more weakly correlated.     

Can aspect ratio be used to categorize intra-aneurysmal hemodynamics?-A study of elastase induced aneurysms in rabbit

Clinical studies suggest that aneurysm aspect ratio (AR) is an important indicator of rupture likelihood. The importance of AR is hypothesized to arise from its influence on intra-aneurysmal hemodynamics. It has been conjectured that slower flow in high AR sacs leads to a cascade of biological activities that weaken the aneurysm wall (Ujiie et al.,1999). However, the connection between AR, hemodynamics and wall weakening has never been proven. Animal models of saccular aneurysms provide a venue for evaluating this conjecture. The focus of this work was to evaluate whether a commonly used elastase induced aneurysm model in rabbits is suitable for a study of this kind from a hemodynamic perspective. In particular, to assess whether hemodynamic factors in low and high AR sacs are statistically different. To achieve this objective, saccular aneurysms were created in 51 rabbits and pulsatile computational fluid dynamics (CFD) studies were performed using rabbit specific inflows. Distinct hemodynamics were found in the low AR (AR<1.8, n=25), and high AR (AR>2.2, n=18) models. A single, stable recirculation zone was present in all low AR aneurysms, whereas a second, transient recirculation zone was also found in the superior aspect of the aneurysm dome for all high AR cases. Aneurysms with AR between 1.8 and 2.2 displayed transitional flow patterns. Differences in values and distributions of hemodynamic parameters were found between low and high AR cases including time averaged wall shear stress, oscillatory shear index, relative residence time and non-dimensional inflow rate. This work lays the foundation for future studies of the dependence of growth and remodeling on AR in the rabbit model and provides a motivation for further studies of the coupling between AR and hemodynamics in human aneurysms.     

Blood flow patterns in the proximal human coronary arteries: relationship to atherosclerotic plaque occurrence

Atherosclerotic plaques in human coronary arteries are focal manifestations of systemic disease, and biomechanical factors have been hypothesized to contribute to plaque genesis and localization. We developed a computational fluid dynamics (CFD) model of the ascending aorta and proximal sections of the right and left coronary arteries of a normal human subject using computed tomography (CT) and magnetic resonance imaging (MRI) and determined the pulsatile flow field. Results demonstrate that flow patterns in the ascending aorta contribute to a pro-atherosclerotic flow environment, specifically through localization of low and oscillatory wall shear stress in the neighborhood of coronary orifices. Furthermore, these patterns differ in their spatial distribution between right and left coronary arteries. Entrance effects of aortic flow diminish within two vessel diameters. We examined relationships between spatial distributions of wall shear stress and reports of plaque occurrence in the literature. Results indicate low wall shear stress is co-located with increased incidence of lesions, and higher wall shear stresses are associated with lesion-resistant areas. This investigation does not consider plaque progression or advanced lesions, inasmuch as the CFD model was developed from a normal individual and the clinical data used for comparisons were obtained from autopsy specimens of subjects who died from non-cardiovascular causes. The data reported are consistent with the hypothesis that low wall shear stress is associated with the localization of atherosclerotic lesions, and the results demonstrate the importance of aortic flow on flow patterns in the proximal segments of the coronary arteries.     

Computational modeling of blood component transport related to coronary artery thrombosis in Kawasaki disease

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case—mechanical or biochemical—could help define improved strategies for thrombosis prevention tailored for each patient.     

Impact of stents and flow diverters on hemodynamics in idealized aneurysm models

Cerebral aneurysms constitute a major medical challenge as treatment options are limited and often associated with high risks. Statistically, up to 3% of patients with a brain aneurysm may suffer from bleeding for each year of life. Eight percent of all strokes are caused by ruptured aneurysms. In order to prevent this rupture, endovascular stenting using so called flow diverters is increasingly being regarded as an alternative to the established coil occlusion method in minimally invasive treatment. Covering the neck of an aneurysm with a flow diverter has the potential to alter the hemodynamics in such a way as to induce thrombosis within the aneurysm sac, stopping its further growth, preventing its rupture and possibly leading to complete resorption. In the present study the influence of different flow diverters is quantified considering idealized patient configurations, with a spherical sidewall aneurysm placed on either a straight or a curved parent vessel. All important hemodynamic parameters (exchange flow rate, velocity, and wall shear stress) are determined in a quantitative and accurate manner using computational fluid dynamics when varying the key geometrical properties of the aneurysm. All simulations are carried out using an incompressible, Newtonian fluid with steady conditions. As a whole, 72 different cases have been considered in this systematic study. In this manner, it becomes possible to compare the efficiency of different stents and flow diverters as a function of wire density and thickness. The results show that the intra-aneurysmal flow velocity, wall shear stress, mean velocity, and vortex topology can be considerably modified thanks to insertion of a suitable implant. Intra-aneurysmal residence time is found to increase rapidly with decreasing stent porosity. Of the three different implants considered in this study, the one with the highest wire density shows the highest increase of intra-aneurysmal residence time for both the straight and the curved parent vessels. The best hemodynamic modifications are always obtained for a small aneurysm diameter.     

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.     

Non-Newtonian blood flow in human right coronary arteries: steady state simulations

This study looks at blood flow through four different right coronary arteries, which have been reconstructed from bi-plane angiograms. Five non-Newtonian blood models, as well as the usual Newtonian model of blood viscosity, are used to study the wall shear stress in each of these arteries at a particular point in the cardiac cycle. It was found that in the case of steady flow in a given artery, the pattern of wall shear stress is consistent across all models. The magnitude of wall shear stress, however, is influenced by the model used and correlates with graphs of shear stress versus strain for each model. For mid-range velocities of around 0.2 m s(-1) the models are virtually indistinguishable. Local and global non-Newtonian importance factors are introduced, in an attempt to quantify the types of flows where non-Newtonian behaviour is significant. It is concluded that, while the Newtonian model of blood viscosity is a good approximation in regions of mid-range to high shear, it is advisable to use the Generalised Power Law model (which tends to the Newtonian model in those shear ranges in any case) in order to achieve better approximation of wall shear stress at low shear.     

Is arterial wall-strain stiffening an additional process responsible for atherosclerosis in coronary bifurcations?: an in vivo study based on dynamic CT and MRI

Coronary bifurcations represent specific regions of the arterial tree that are susceptible to atherosclerotic lesions. While the effects of vessel compliance, curvature, pulsatile blood flow, and cardiac motion on coronary endothelial shear stress have been widely explored, the effects of myocardial contraction on arterial wall stress/strain (WS/S) and vessel stiffness distributions remain unclear. Local increase of vessel stiffness resulting from wall-strain stiffening phenomenon (a local process due to the nonlinear mechanical properties of the arterial wall) may be critical in the development of atherosclerotic lesions. Therefore, the aim of this study was to quantify WS/S and stiffness in coronary bifurcations and to investigate correlations with plaque sites. Anatomic coronary geometry and cardiac motion were generated based on both computed tomography and MRI examinations of eight patients with minimal coronary disease. Computational structural analyses using the finite element method were subsequently performed, and spatial luminal arterial wall stretch (LW(Stretch)) and stiffness (LW(Stiff)) distributions in the left main coronary bifurcations were calculated. Our results show that all plaque sites were concomitantly subject to high LW(Stretch) and high LW(Stiff), with mean amplitudes of 34.7 ± 1.6% and 442.4 ± 113.0 kPa, respectively. The mean LW(Stiff) amplitude was found slightly greater at the plaque sites on the left main coronary artery (mean value: 482.2 ± 88.1 kPa) compared with those computed on the left anterior descending and left circumflex coronary arteries (416.3 ± 61.5 and 428.7 ± 181.8 kPa, respectively). These findings suggest that local wall stiffness plays a role in the initiation of atherosclerotic lesions.     

Investigation of hemodynamics in the development of dissecting aneurysm within patient-specific dissecting aneurismal aortas using computational fluid dynamics (CFD) simulations

Aortic dissecting aneurysm is one of the most catastrophic cardiovascular emergencies that carries high mortality. It was pointed out from clinical observations that the aneurysm development is likely to be related to the hemodynamics condition of the dissected aorta. In order to gain more insight on the formation and progression of dissecting aneurysm, hemodynamic parameters including flow pattern, velocity distribution, aortic wall pressure and shear stress, which are difficult to measure in vivo, are evaluated using numerical simulations. Pulsatile blood flow in patient-specific dissecting aneurismal aortas before and after the formation of lumenal aneurysm (pre-aneurysm and post-aneurysm) is investigated by computational fluid dynamics (CFD) simulations. Realistic time-dependent boundary conditions are prescribed at various arteries of the complete aorta models. This study suggests the helical development of false lumen around true lumen may be related to the helical nature of hemodynamic flow in aorta. Narrowing of the aorta is responsible for the massive recirculation in the poststenosis region in the lumenal aneurysm development. High pressure difference of 0.21 kPa between true and false lumens in the pre-aneurismal aorta infers the possible lumenal aneurysm site in the descending aorta. It is also found that relatively high time-averaged wall shear stress (in the range of 4-8 kPa) may be associated with tear initiation and propagation. CFD modeling assists in medical planning by providing blood flow patterns, wall pressure and wall shear stress. This helps to understand various phenomena in the development of dissecting aneurysm.     

Aortic valve leaflet wall shear stress characterization revisited: impact of coronary flow

Computational characterizations of aortic valve hemodynamics have typically discarded the effects of coronary flow. The objective of this study was to complement our previous fluid–structure interaction aortic valve model with a physiologic coronary circulation model to quantify the impact of coronary flow on aortic sinus hemodynamics and leaflet wall shear stress (WSS). Coronary flow suppressed vortex development in the two coronary sinuses and altered WSS magnitude and directionality on the three leaflets, with the most substantial differences occurring in the belly and tip regions.     

Wall shear over high degree stenoses pertinent to atherothrombosis

Atherothrombosis can induce acute myocardial infarction and stroke by progressive stenosis of a blood vessel lumen to full occlusion. Since thrombus formation and embolization may be shear-dependent, we quantify the magnitude of shear rates in idealized severely stenotic coronary arteries (≥75% by diameter) using computational fluid dynamics to characterize the shear environment that may exist during atherothrombosis. Maximum shear rates in severe short stenoses were found to exceed 250,000 s^?1 (9500 dynes/cm²) and can reach a peak value of 425,000 s^?1 for a 98% stenosis. These high shear rates exceed typical shear used for in vitro blood flow experiments by an order of magnitude, indicating the need to examine thrombosis at very high shear rates. Pulsatility and stenosis eccentricity were found to have minor effects on the maximum wall shear rates in severe stenoses. In contrast, increases in the stenosis length reduced the maximum shear to 107,000 s^?1 (98% stenosis), while surface roughness could increase focal wall shear rates to a value reaching 610,000 s^?1 (90% stenosis). The “shear histories” of circulating platelets in these stenoses are far below reported activation thresholds. Platelets may be required to form bonds in 5 μs and resist shear forces reaching 8000 pN per platelet. Arterial thrombosis occurs in the face of pathological high shear stress, creating rapid and strong bonds without prior activation of circulating platelets.     

Hemodynamics of the mouse abdominal aortic aneurysm

The abdominal aortic aneurysm (AAA) is a significant cause of death and disability in the Western world and is the subject of many clinical and pathological studies. One of the most commonly used surrogates of the human AAA is the angiotensin II (Ang II) induced model used in mice. Despite the widespread use of this model, there is a lack of knowledge concerning its hemodynamics; this study was motivated by the desire to understand the fluid dynamic environment of the mouse AAA. Numerical simulations were performed using three subject-specific mouse models in flow conditions typical of the mouse. The numerical results from one model showed a shed vortex that correlated with measurements observed in vivo by Doppler ultrasound. The other models had smaller aneurysmal volumes and did not show vortex shedding, although a recirculation zone was formed in the aneurysm, in which a vortex could be observed, that elongated and remained attached to the wall throughout the systolic portion of the cardiac cycle. To link the hemodynamics with aneurysm progression, the remodeling that occurred between week one and week two of the Ang II infusion was quantified and compared with the hemodynamic wall parameters. The strongest correlation was found between the remodeled distance and the oscillatory shear index, which had a correlation coefficient greater than 0.7 for all three models. These results demonstrate that the hemodynamics of the mouse AAA are driven by a strong shear layer, which causes the formation of a recirculation zone in the aneurysm cavity during the systolic portion of the cardiac waveform. The recirculation zone results in areas of quiescent flow, which are correlated with the locations of the aneurysm remodeling.