†These authors contributed equally.
Academic Editor: Peter A. McCullough
Backgrounds: Coronary cavitation is supposed to be generated by both
concentric and eccentric coronary artery stenosis which propagates downstream the
vessel, creating microbubbles which exploded when the fluid pressure was lower
than the vapor pressure at a local thermodynamic state. Objective: To assess,
using numerical and computational fluid dynamic analysis (CFD), the potential of
cavitation to both induce damage to coronary artery endothelium and to promote
atherosclerotic plaque progression. Methods: We retrospectively reviewed
the data 12 consecutive patients evaluated between 1st January 2013 and 1st
January 2014 with an isolated hemodynamically significant Left Main (LM) disease.
The patient specific geometries have been reconstructed. Bubble velocity has been
calculated in accordance with Newton’s second law. Both the forces arising from
the bubbles’ interaction with the continuous phase and impact with the
endothelium have been evaluated. The impact of turbulence on the motion of
bubbles have been modelled with a dispersion model. Results: Among the
12 patients retrospectively analysed [8 males, mean age 68.2
Atherosclerosis is currently considered as an inflammatory and/or
immunomodulatory reaction due to the presence of oxidized low-density
lipoproteins (LDL) within the arterial wall concomitantly in the setting of
traditional cardiovascular risk factors [1, 2, 3, 4]. Nowadays, the triggers of coronary
lesion progression and vulnerability from an asymptomatic fibro-atheromatous to a
“vulnerable plaque” are not yet fully understood . As we previously
demonstrated, a rheological phenomenon called cavitation is generated by both
concentric and eccentric coronary artery stenosis (
We retrospectively analysed the data of 12 consecutive patients evaluated at
Rovigo General Hospital between 1st January 2013 and 1st January 2014 reporting
an isolated hemodynamically significant LM body/distal disease (
The CCTA protocol has been already described . Briefly, CCTA was acquired using a 64-slice multi-detector computed tomography scanner (64-detector row Lightspeed VCT scanner, GE Healthcare, Milwaukee, WI, USA). Electrocardiographic (ECG) gating was used for all scans. The scanning protocol was adjusted for patient weight and HR. Blood pressure was measured with an arm-band sphygmomanometer just before the examination. Each examination was independently interpreted by a trained CT radiologist using a combination of axial images, 3D volume-rendered images, multiplanar reformations, and maximum intensity projections. CA has been performed using a 6F radial approach whenever possible and bilateral angiographic 5F catheter. Angiographic data were independently reviewed by two 20-year experience angiographic operator (G.R. and L.R.). Any discordances were resolved via collegial discussion.
As previously performed  the computational domain and patient specific geometries have been reconstructed using the open-source software package Vascular Modelling Toolkit (VMTK, https://www.vmtk.org), identifying the vessel centreline and then splitting the branches of the bifurcation. Therefore, external wall surfaces were reconstructed using circumferential cross sections perpendicular to the centrelines using the CAD software Rhinoceros 4.0 Evaluation (McNeel & Associates, Indianapolis, IN, USA) . Atherosclerotic plaques have been modelled considering the distance between each node and the centreline of the external wall surface. Coronary artery stenosis (DS%) and minimal lumen area have been calculated considering the reconstructed coronary artery model. Specifically, the severity of the stenosis has been computed as 100% minus the percentage of minimal lumen area to the reference lumen area. Subsequently, the patient-specific geometries have been discretized using Ansys ICEM CFD (ANSYS Inc., Canonsburg, PA, USA) into linear tetrahedral volume elements which had a maximum size of 0.2 mm at the model surface. Meshing parameters included a default global maximal wall-face edge sizes of 1.0 and 0.15 mm, for the proximal length of both LAD and LCX (defined as the first fertile of the vessels lengths) and carina, respectively. Moreover, a mesh independence analysis was performed before the CFD simulation using the region adaptation function in FLUENT to refine the mesh until the modification in mesh size did not produce results with significant differences. A prism layer consisting of four layers with a height ratio of 1.2 was employed to capture the flow close to the lumen . The mesh sizes used into these simulations were 280.090 elements and 45.375 nodes, respectively.
Blood has been modelled using a laminar non-Newtonian and incompressible fluid using the Navier–Stokes Eqn. 1
and the continuity Eqn. 2
where v is the 3D velocity vector, P the pressure,
Vapor has been modelled as discrete vapour bubbles whose trajectory is determined using a Lagrangian frame of reference . Cavitation has been assumed to stary by cavitation nuclei which subsequently grow into bubbles undergoing different physical processes, which have been determined in a stochastic Monte-Carlo approximation . The size of the nuclei has been sampled from a probability density function. When the pressure of blood fallen below its vapor pressure, the volume under tension was identified and the most probable locations for bubble nuclei formation were randomly calculated by sampling from a distribution function which considers the bubble volume fraction, volume and the non-dimensional tension of each Eulerian cell .
Bubble velocity has been calculated in accordance with Newton’s second law. The forces arising from the bubbles’ interaction with the continuous phase and impact with the endothelium have been also evaluated. The effect of turbulence on the motion of bubbles have been modelled with a dispersion model [20, 21].
During the study period, 15 patients were retrospectively identified. Of them,
three were excluded due the presence of concomitant significant CAD. As a result,
12 consecutive patients [9 males, mean age 74.1
|N = 12|
|Mean or (%)|
|Age (years)||74.1 |
|Gender (Male) %||9 (75)|
|Overweight, %||7 (58.3)|
|Obese, %||1 (8.3)|
|SBP (mmHg)||140.5 |
|DBP (mmHg)||77.3 |
|HR (bpm)||82.2 |
|TC (mg/dL)||222.4 |
|TG (mg/dL)||168.3 |
|HDL (mg/dL)||32.3 |
|Diabetes, %||6 (50)|
|HT, %||11 (91.6)|
|Lesion location on angiography|
|-Ostial only||0 (0)|
|BMI, body mass index; SBP, Systolic blood pressure; DBP, Diastolic blood pressure; HR, Heart rate; TC, Total cholesterol; TG, triglycerides; HDL, high density lipoprotein; LDL, Low density lipoprotein; IGF, Impaired fasting glucose; HT, arterial hypertension.|
Cavitation profile into a left main stenosis. (A) The simulation illustrates the vapor fraction iso-surfaces and the scattered bubble plots as predicted by the Langrangian model. Notably, most cavitation bubbles that form at the inlet of left main bifurcation do not collapse immediately but they are transported towards the vessel determining several interactions with the endothelium. (B) Graphical representation of the bubble radius modification and related pressure transmitted to the endothelium if the bubble collapse happened near to this last one. The magnification in boxes (B1–B4) evidences the dynamic modification of bubbles. The re-entry jet causing the collapse is evidenced with a red arrow.
As expectable from the general cavitation theory, due to the pressure gradient
generated by the stenosis as well as the relative turbulent flow, bubbles firstly
increase their size (Fig. 1, Panels B1–B2) until their collapse above the vapor
pressure level (Fig. 1, Panel B3). This event was promoted by the formation of a
re-entry jet which penetrates each bubble (Fig. 1, Panel B4, red arrow)
generating a shock wave and, if near the endothelium, transferring the related
energy observed as an instantaneous pressure-peaks on the coronary surface. The
bubble collapse happened quickly, after 2.1
The present study evaluates the relation between flow dynamics and cavitation phenomenon in promoting atherosclerotic plaque progression into coronary arteries. The results evidenced that shocks and instantaneous pressure-peaks, determined by the collapses of vapor bubbles, resulted in a measurable dynamic loads on endothelium which has the potential of causing damage or erosion. In this regard, CFD represents an important technique for providing non-invasive detailed blood flow data.
Obviously, the complexity of the mathematical analysis inherent to our demonstration has been simplified in order to be of practical value for the non-bioengineers professionals, maintaining the scientific quality of the analysis.
From a physical point of view, cavitation occurs when the pressure (P), in an area of flowing fluid, decreases below the level of vapor pressure (Pv). More precisely, cavitation bubbles form immediately below the (Pv) level, reaching the highest velocity at the lowest pressure level (Pmin), which is defined as vena contracta point. Subsequently, bubbles increase their size with the pressure increasing, leading to bubble collapse above the Pv level. This event represents the result of two different aspects, which form the essential substrates of cavitation: the former is characterized by a large pressure fluctuation (dp/dt), while the latter is due to the turbulent flow . Conversely, the cavitation phenomenon in a hyperelastic or inhomogeneous soft solid occurs when in a neo-Hookean solid (which is an hyperleastic material model that can be used for predicting the non-linear stress-strain behaviour of materials undergoing deformation) under sufficient loading conditions with large triaxial forces cavities appear . From a translational point of view, the cavitation within a fluid (e.g., blood) needs the presence of a significant vessel stenosis to be generated while the cavitation of the hyperelastic material (e.g., the atherosclerotic plaque) needs sufficient triaxialities to be generated.
The hypothesis that the energy released by microbubbles collapse on the plaque surface be able to produce a significant destabilization of the plaque, promoting cardiovascular events, is similar to the findings reported some years before, regarding the occurrence of haemolysis and generation of valvular lesions, due to the cavitation, observed in patients with MHVs . The cavitation phenomenon, alone, is not able to explain the growth as well as the structural and biochemical modifications of the coronary atherosclerotic plaque but may act as a concurrent rather a causal factor in the precipitation of the atherosclerotic disease. Indeed, the cavitation is observed in bloodstream, in accordance to the fluid dynamic laws, only when a significant abrupt stenosis [6, 25], even concentric or eccentric, determines the sudden decrease of P below the Pv . Similarly, when the cavitation is analysed into the plaque, which has been compared to a soft inhomogeneous solid, requires a specific loading conditions to become significant and able to determine the “failure” of the plaque . Intriguingly, some studies have evidenced that the cavitation can be used to induce a localized and controlled atherosclerotic process in the abdominal aorta .
Doubtless, due the relative few number of patients enrolled and considering the peculiar location of the stenotic lesion as well as the distance to the carina, and the bifurcation angle, which may significantly vary among patients, our preliminary data must be cautiously considered since their integration into may have led to some intrinsic bias.
The cavitation phenomenon probably produces a specific rheologic effect which contribute to plaque progression. In other words, the lesion induced by long-lasting and continuous explosion of bubbles colliding with the endothelium may transfer enough energy to the vessel wall to generate an endothelial injury which triggered the progression of the plaque by the well-known cascade.
Our study suffers from several intrinsic limitations which firstly include its retrospective single center small patients sample nature. Secondly, the application of a phenomena that usually produces erosion on different surface than the endothelial wall of the coronary tree. Thirdly, we considered a diastolic steady flow for the analysis: in coronary vascular tree the flow is mainly diastolic with a small systolic compound and we cannot avoid a mild model distortion induced by this latter flow component. Fourthly the plaque has not been characterized in its components, fibrous, fibrous -fatty, necrotic core, calcified, because of the huge number of different model combinations which would make the analysis inconsistent. Finally, the absence of a histological confirmation of the damage which obviously cannot be made in vivo.
Our study demonstrated by means of numerical and CFD analysis that in coronary model of LM bifurcation with a significant disease, the cavitation phenomena is effectively generated and instantaneous pressure-peaks due to collapses of vapor bubbles resulted in dynamic loads on endothelium which reached a measurable mean pressure of 3.9 MPa, which may have the potential to provoke endothelial damage. This demonstration offers the basis for further experimental studies about the significance of cavitation in coronary plaque progression.
GR—Conceptualization; write the manuscript; reviewed the manuscript; data analysis; supervision; MZ—Conceptualization; write the manuscript; reviewed the manuscript; data analysis; CB—Critical revision and final supervision; TN—Conceptualization; write the manuscript; reviewed the manuscript; data analysis.
Ethical approval was waived due to the retrospective nature of the study. All partecipants provide the consent to partecipate to the study.
This research received no external funding.
The authors declare no conflict of interest.