Validation of biomechanical assessment of coronary plaque vulnerability based on intravascular optical coherence tomography and digital subtraction angiography

Background

It has been suggested that biomechanical factors may influence plaque development. However, key determinants for assessing plaque vulnerability remain speculative.

Methods

In this study, a two-dimensional (2D) structural mechanical analysis and a three-dimensional (3D) fluid-structure interaction (FSI) analysis were conducted based on intravascular optical coherence tomography (IV-OCT) and digital subtraction angiography (DSA) data sets. In the 2D study, 103 IV-OCT slices were analyzed. An in-depth morpho-mechanic analysis and a weighted least absolute shrinkage and selection operator (LASSO) regression analysis were conducted to identify the crucial features related to plaque vulnerability via the tuning parameter (λ). In the 3D study, the coronary model was reconstructed by fusing the IV-OCT and DSA data, and a FSI analysis was subsequently performed. The relationship between vulnerable plaque and wall shear stress (WSS) was investigated.

Results

The influential factors were selected using the minimum criteria (λ-min) and one-standard error criteria (λ-1se). In addition to the common vulnerable factor of the minimum fibrous cap thickness (FCTmin), four biomechanical factors were selected by λ-min, including the average/maximal displacements and average/maximal stress, and two biomechanical factors were selected by λ-1se, including the average/maximal displacements. Additionally, the positions of the vulnerable plaques were consistent with the sites of high WSS.

Conclusions

Functional indices are crucial for plaque status assessment. An evaluation based on biomechanical simulations might provide insights into risk identification and guide therapeutic decisions.

Patient-specific air puff-induced loading using machine learning

Introduction: The air puff test is a contactless tonometry test used to measure the intraocular pressure and the cornea's biomechanical properties. Limitations that most challenge the accuracy of the estimation of the corneal material and the intraocular pressure are the strong intercorrelation between the intraocular pressure and the corneal parameters, either the material properties that can change from one person to another because of age or the geometry parameters like central corneal thickness. This influence produces inaccuracies in the corneal deformation parameters while extracting the IOP parametric equation, which can be reduced through the consideration of the patient-specific air puff pressure distribution taking into account the changes in corneal parameters. This air puff pressure loading distribution can be determined precisely from the fluid-structure interaction (FSI) coupling between the air puff and the eye model. However, the computational fluid dynamics simulation of the air puff in the coupling algorithm is a time-consuming model that is impractical to use in clinical practice and large parametric studies. Methods: By using a supervised machine learning algorithm, we predict the time-dependent air puff pressure distribution for different corneal parameters via a parametric study of the corneal deformations and the gradient boosting algorithm. Results: The results confirmed that the algorithm gives the time-dependent air puff pressure distribution with an MAE of 0.0258, an RMSE of 0.0673, and an execution time of 93 s, which is then applied to the finite element model of the eye generating the corresponding corneal deformations taking into account the FSI influence. Using corneal deformations, the response parameters can be extracted and used to produce more accurate algorithms of the intraocular pressure and corneal material stress-strain index (SSI). Discussion: Estimating the distribution of air pressure on the cornea is essential to increase the accuracy of intraocular pressure (IOP) measurements, which serve as valuable indicator of corneal disease. We find that the air puff pressure loading is largely influenced by complex changes in corneal parameters unique to each patient case. With our innovative algorithm, we can preserve the same accuracy developed by the CFD-based FSI model, while reducing the computational time from approximately 101000 s (28 h) to 720 s (12 min), which is about 99.2% reduction in time. This huge improvement in computational cost will lead to significant improvement in the parametric equations for IOP and the Stress-Strain Index (SSI).

Hemodynamic performance evaluation of neonatal ECMO double lumen cannula using fluid-structure interaction

Extra corporeal membrane oxygenation (ECMO) is an artificial oxygenation facility, employed in situations of cardio-pulmonary failure. Some diseases i.e., acute respiratory distress syndrome, pulmonary hypertension, corona virus disease (COVID-19) etc. affect oxygenation performance of the lungs thus requiring the need of artificial oxygenation. Critical care teams used ECMO technique during the COVID-19 pandemic to support the heart and lungs of COVID-19 patients who had an acute respiratory or cardiac failure. Double Lumen Cannula (DLC) is one of the most critical components of ECMO as it resides inside the patient and, connects patient with external oxygenation circuit. DLC facilitates delivery and drainage of blood from the patient's body. DLC is characterized by delicate balance of internal and external flows inside a limited space of the right atrium (RA). An optimal performance of the DLC necessitates structural stability under biological and hemodynamic loads, a fact that has been overlooked by previously published studies. In the past, many researchers experimentally and computationally investigated the hemodynamic performance of DLC by employing Eulerian approach, which evaluate instantaneous blood damage without considering blood shear exposure history (qualitative assessment only). The present study is an attempt to address the aforementioned limitations of the previous studies by employing Lagrangian (quantitative assessment) and incorporating the effect of fluid-structure interaction (FSI) to study the hemodynamic performance of neonatal DLC. The study was performed by solving three-dimensional continuity, momentum, and structural mechanics equation(s) by numerical methods for the blood flow through neonatal DLC. A two-way coupled FSI analysis was performed to analyze the effect of DLC structural deformation on its hemodynamic performance. Results show that the return lumen was the most critical section with maximum pressure drop, velocity, shear stresses, and blood damage. Recirculation and residence time of blood in the right atrium (RA) increases with increasing blood flow rates. Considering the structural deformation has led to higher blood damage inside the DLC-atrium system. Maximum Von-Mises stress was present on the side edges of the return lumen that showed direct proportionality with the blood flow rate.

On the Potential Self-Amplification of Aneurysms Due to Tissue Degradation and Blood Flow Revealed From FSI Simulations

Tissue degradation plays a crucial role in the formation and rupture of aneurysms. Using numerical computer simulations, we study the combined effects of blood flow and tissue degradation on intra-aneurysm hemodynamics. Our computational analysis reveals that the degradation-induced changes of the time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI) within the aneurysm dome are inversely correlated. Importantly, their correlation is enhanced in the process of tissue degradation. Regions with a low TAWSS and a high OSI experience still lower TAWSS and higher OSI during degradation. Furthermore, we observed that degradation leads to an increase of the endothelial cell activation potential index, in particular, at places experiencing low wall shear stress. These findings are robust and occur for different geometries, degradation intensities, heart rates and pressures. We interpret these findings in the context of recent literature and argue that the degradation-induced hemodynamic changes may lead to a self-amplification of the flow-induced progressive damage of the aneurysmal wall.

Transient numerical simulation of the right coronary artery originating from the left sinus and the effect of its acute take-off angle on hemodynamics

Background

An anomalous origin of the right coronary artery from the left coronary artery sinus is usually characterized by an acute take-off angle. Most affected patients have no clinical symptoms; however, some patients have decreased blood flow into the right coronary artery during exercise, which can lead to symptoms such as myocardial ischemia. Most researchers who have studied an anomalous origin of the right coronary artery from the left coronary artery sinus have done so through clinical cases. In this study, we used numerical simulation to evaluate the hemodynamics of this condition and the effect of an acute take-off angle on hemodynamic parameters. We expect that the results of this study will help in further understanding the clinical symptoms of this anomaly and the hemodynamic impact of an acute take-off angle.

Methods

Three-dimensional models were reconstructed based on the computed tomography images from 16 patients with a normal right coronary artery and 26 patients with an anomalous origin of the right coronary artery from the left coronary artery sinus. A numerical simulation of a two-way fluid-structure interaction was executed with ANSYS Workbench software. The blood was assumed to be an incompressible Newtonian fluid, and the vessel was assumed to be an isotropic, linear elastic material. Hemodynamic parameters and the effect of an acute take-off angle were statistically analyzed.

Results

During the systolic period, the wall pressure in the right coronary artery was significantly reduced in patients with an anomalous origin of the right coronary artery (t =1.32 s, P=0.0001; t =1.34-1.46 s, P<0.0001). The wall shear stress in the abnormal group was higher at the beginning of the systolic period (t =1.24 s, P=0.0473; t =1.26 s, P=0.0193; t =1.28 s, P=0.0441). The acute take-off angle was smaller in patients with clinical symptoms (27.81°±4.406°) than in patients without clinical symptoms (31.86°±2.789°; P=0.017). In the symptomatic group, pressure was negatively correlated with the acute take-off angle (P=0.0185-0.0341, r=-0.459 to -0.4167).

Conclusions

This study shows that an anomalous origin of the right coronary artery from the left coronary artery sinus causes changes in hemodynamic parameters, and that an acute take-off angle in patients with this anomaly is associated with terminal ischemia of the right coronary artery.

Numerical study of hemodynamics in a complete coronary bypass with venous and arterial grafts and different degrees of stenosis

Cardiovascular diseases are among the leading causes of death in the world. The coronary blockage is one of most common types of these diseases that in the majority of cases has been treated by bypass surgery. In the bypass surgery, a graft is implemented to alter the blocked coronary and allow the blood supply process. The hemodynamic characteristics of the bypass strongly depend on the geometry and mechanical properties of the graft. In the present study, the fluid-structure interaction (FSI) analysis is conducted to investigate the bypass performance for a thoracic artery as well as a saphenous vein graft. Blood flow introduces a pressure on the walls of the graft which behaves as a hyperelastic material. A complete coronary bypass with stenosis degrees of 70% and 100% is modeled. To consider the nonlinear stress-strain behavior of the grafts, a five parameter Mooney-Rivlin hyperplastic model is implemented for the structural analysis and blood is assumed to behave as a Newtonian fluid. The simulations are performed for a structured grid to solve the governing equations using finite element method (FEM). The results show that wall shear stress (WSS) for saphenous vein is larger than that of thoracic artery while the total deformation of the thoracic artery is larger compared to the saphenous vein. Also, for the venous grafts or lower stenosis degree, the oscillatory shear index (OSI) is higher at both left and right anastomoses meaning that venous grafts as well as lower degree of stenosis are more critical in terms of restenosis.

Effect of intimal flap motion on flow in acute type B aortic dissection by using fluid-structure interaction

A monolithic, fully coupled fluid-structure interaction (FSI) computational framework was developed to account for dissection flap motion in acute type B aortic dissection (TBAD). Analysis of results included wall deformation, pressure, flow, wall shear stress (WSS), von Mises stress and comparison of hemodynamics between rigid wall and FSI models. Our FSI model mimicked realistic wall deformation that resulted in maximum compression of the distal true lumen (TL) by 21.4%. The substantial movement of intimal flap mostly affected flow conditions in the false lumen (FL). Flap motion facilitated more flow entering the FL at peak systole, with the TL to FL flow split changing from 88:12 in the rigid model to 83:17 in the FSI model. There was more disturbed flow in the FL during systole (5.8% FSI vs 5.2% rigid) and diastole (13.5% FSI vs 9.8% rigid), via a λ₂ -criterion. The flap-induced disturbed flow near the tears in the FSI model caused an increase of local WSS by up to 70.0% during diastole. This resulted in a significant reduction in the size of low time-averaged WSS (TAWSS) regions in the FL (113.11 cm² FSI vs 177.44 cm² rigid). Moreover, the FSI model predicted lower systolic pressure, higher diastolic pressure, and hence lower pulse pressure. Our results provided new insights into the possible impact of flap motion on flow in aortic dissections, which are particularly important when evaluating hemodynamics of acute TBAD. NOVELTY STATEMENT: Our monolithic fully coupled FSI computational framework is able to reproduce experimentally measured range of flap deformation in aortic dissection, thereby providing novel insights into the influence of physiological flap motion on the flow and pressure distributions. The drastic flap movement increases the flow resistance in the true lumen and causes more disturbed flow in the false lumen, as visualized through the λ₂ criterion. The flap-induced luminal pressure is dampened, thereby affecting pressure measures, which may serve as potential prognostic indicators for late complications in acute uncomplicated TBAD patients.

Fluid-Structure Interaction Analysis on Membrane Behavior of a Microfluidic Passive Valve

In this paper, the effect of membrane features on flow characteristics in the microfluidic passive valve (MPV) and the membrane behavior against fluid flow are studied using the fluid-structure interaction (FSI) analysis. Firstly, the microvalve model with different numbers of microholes and pitches of microholes are designed to investigate the flow rate of the MPV. The result shows that the number of microholes on the membrane has a significant impact on the flow rate of the MPV, while the pitch of microholes has little effect on it. The constant flow rate maintained by the microvalve (the number of microholes n = 4) is 5.75 mL/min, and the threshold pressure to achieve the flow rate is 4 kPa. Secondly, the behavior of the membrane against the fluid flow is analyzed. The result shows that as the inlet pressure increases, the flow resistance of the MPV increases rapidly, and the deformation of the membrane gradually becomes stable. Finally, the effect of the membrane material on the flow rate and the deformation of the membrane are studied. The result shows that changes in the material properties of the membrane cause a decrease in the amount of deformation in all stages the all positions of the membrane. This work may provide valuable guidance for the optimization of microfluidic passive valve in microfluidic system.

Fluid-Structure Interaction Based Algorithms for IOP and Corneal Material Behavior

Purpose: This paper presents and clinically validates two algorithms for estimating intraocular pressure (IOP) and corneal material behavior using numerical models that consider the fluid-structure interaction between the cornea and the air-puff used in non-contact tonometry. Methods: A novel multi-physics fluid-structure interaction model of the air-puff test was employed in a parametric numerical study simulating human eyes under air-puff pressure with a wide range of central corneal thickness (CCT = 445-645 μm), curvature (R = 7.4-8.4 mm), material stiffness and IOP (10-25 mmHg). Models were internally loaded with IOP using a fluid cavity, then externally with air-puff loading simulated using a turbulent computational fluid dynamics model. Corneal dynamic response parameters were extracted and used in development of two algorithms for IOP and corneal material behavior; fIOP and fSSI, respectively. The two algorithms were validated against clinical corneal dynamic response parameters for 476 healthy participants. The predictions of IOP and corneal material behavior were tested on how they varied with CCT, R, and age. Results: The present study produced a biomechanically corrected estimation of intraocular pressure (fIOP) and a corneal material stiffness parameter or Stress-Strain Index (fSSI), both of which showed no significant correlation with R (p > 0.05) and CCT (p > 0.05). Further, fIOP had no significant correlation with age (p > 0.05), while fSSI was significantly correlated with age (p = 0.001), which was found earlier to be strongly correlated with material stiffness. Conclusion: The present study introduced two novel algorithms for estimating IOP and biomechanical material behavior of healthy corneas in-vivo. Consideration of the fluid structure interaction between the cornea and the air puff of non-contact tonometry in developing these algorithms led to improvements in performance compared with bIOP and SSI.

Simulation Analysis of Fluid-Structure Interaction of High Velocity Environment Influence on Aircraft Wing Materials under Different Mach Numbers

Aircraft service process is in a state of the composite load of pressure and temperature for a long period of time, which inevitably affects the inherent characteristics of some components in aircraft accordingly. The flow field of aircraft wing materials under different Mach numbers is simulated by Fluent in order to extract pressure and temperature on the wing in this paper. To determine the effect of coupling stress on the wing’s material and structural properties, the fluid-structure interaction (FSI) method is used in ANSYS-Workbench to calculate the stress that is caused by pressure and temperature. Simulation analysis results show that with the increase of Mach number, the pressure and temperature on the wing’s surface both increase exponentially and thermal stress that is caused by temperature will be the main factor in the coupled stress. When compared with three kinds of materials, titanium alloy, aluminum alloy, and Haynes alloy, carbon fiber composite material has better performance in service at high speed, and natural frequency under coupling pre-stressing will get smaller.

Impact of complex blast waves on the human head: a computational study

Head injuries due to complex blasts are not well examined because of limited published articles on the subject. Previous studies have analyzed head injuries due to impact from a single planar blast wave. Complex or concomitant blasts refer to impacts usually caused by more than a single blast source, whereby the blast waves may impact the head simultaneously or consecutively, depending on the locations and distances of the blast sources from the subject, their blast intensities, the sequence of detonations, as well as the effect of blast wave reflections from rigid walls. It is expected that such scenarios will result in more serious head injuries as compared to impact from a single blast wave due to the larger effective duration of the blast. In this paper, the utilization of a head-helmet model for blast impact analyses in Abaqus(TM) (Dassault Systemes, Singapore) is demonstrated. The model is validated against studies published in the literature. Results show that the skull is capable of transmitting the blast impact to cause high intracranial pressures (ICPs). In addition, the pressure wave from a frontal blast may enter through the sides of the helmet and wrap around the head to result in a second impact at the rear. This study recommended better protection at the sides and rear of the helmet through the use of foam pads so as to reduce wave entry into the helmet. The consecutive frontal blasts scenario resulted in higher ICPs compared with impact from a single frontal blast. This implied that blast impingement from an immediate subsequent pressure wave would increase severity of brain injury. For the unhelmeted head case, a peak ICP of 330 kPa is registered at the parietal lobe which exceeds the 235 kPa threshold for serious head injuries. The concurrent front and side blasts scenario yielded lower ICPs and skull stresses than the consecutive frontal blasts case. It is also revealed that the additional side blast would only significantly affect ICPs at the temporal and parietal lobes when compared with results from the single frontal blast case. By analyzing the pressure wave flow surrounding the head and correlating them with the consequential evolution of ICP and skull stress, the paper provides insights into the interaction mechanics between the concomitant blast waves and the biological head model.