Optical coherence tomography-based patient-specific coronary artery reconstruction and fluid–structure interaction simulation

Patient-Specific Numerical Simulations of Coronary Artery Hemodynamics and Biomechanics: A Pathway to Clinical Use

PURPOSE

Numerical models that simulate the behaviors of the coronary arteries have been greatly improved by the addition of fluid-structure interaction (FSI) methods. Although computationally demanding, FSI models account for the movement of the arterial wall and more adequately describe the biomechanical conditions at and within the arterial wall. This offers greater physiological relevance over Computational Fluid Dynamics (CFD) models, which assume the walls do not move or deform. Numerical simulations of patient-specific cases have been greatly bolstered by the use of imaging modalities such as Computed Tomography Angiography (CTA), Magnetic Resonance Imaging (MRI), Optical Coherence Tomography (OCT), and Intravascular Ultrasound (IVUS) to reconstruct accurate 2D and 3D representations of artery geometries. The goal of this study was to conduct a comprehensive review on CFD and FSI models on coronary arteries, and evaluate their translational potential.

METHODS

This paper reviewed recent work on patient-specific numerical simulations of coronary arteries that describe the biomechanical conditions associated with atherosclerosis using CFD and FSI models. Imaging modality for geometry collection and clinical applications were also discussed.

RESULTS

Numerical models using CFD and FSI approaches are commonly used to study biomechanics within the vasculature. At high temporal and spatial resolution (compared to most cardiac imaging modalities), these numerical models can generate large amount of biomechanics data.

CONCLUSIONS

Physiologically relevant FSI models can more accurately describe atherosclerosis pathogenesis, and help to translate biomechanical assessment to clinical evaluation.

Fully automated construction of three-dimensional finite element simulations from Optical Coherence Tomography

Despite recent advances in diagnosis and treatment, atherosclerotic coronary artery diseases remain a leading cause of death worldwide. Various imaging modalities and metrics can detect lesions and predict patients at risk; however, identifying unstable lesions is still difficult. Current techniques cannot fully capture the complex morphology-modulated mechanical responses that affect plaque stability, leading to catastrophic failure and mute the benefit of device and drug interventions. Finite Element (FE) simulations utilizing intravascular imaging OCT (Optical Coherence Tomography) are effective in defining physiological stress distributions. However, creating 3D FE simulations of coronary arteries from OCT images is challenging to fully automate given OCT frame sparsity, limited material contrast, and restricted penetration depth. To address such limitations, we developed an algorithmic approach to automatically produce 3D FE-ready digital twins from labeled OCT images. The 3D models are anatomically faithful and recapitulate mechanically relevant tissue lesion components, automatically producing morphologies structurally similar to manually constructed models whilst including more minute details. A mesh convergence study highlighted the ability to reach stress and strain convergence with average errors of just 5.9% and 1.6% respectively in comparison to FE models with approximately twice the number of elements in areas of refinement. Such an automated procedure will enable analysis of large clinical cohorts at a previously unattainable scale and opens the possibility for in-silico methods for patient specific diagnoses and treatment planning for coronary artery disease.

Comparison of multilayer and single-layer coronary plaque models on stress/strain calculations based on optical coherence tomography images

Mechanical stress and strain conditions are closely related to atherosclerotic plaque progression and rupture and have been under intensive investigations in recent years. It is well known that arteries have a three-layer structure: intima, media and adventitia. However, in vivo image-based multilayer plaque models are not available in the current literature due to lack of multilayer image segmentation data. A multilayer segmentation and repairing technique was introduced to segment coronary plaque optical coherence tomography (OCT) image to obtain its three-layer vessel structure. A total of 200 OCT slices from 20 patients (13 male; 7 female) were used to construct multilayer and single-layer 3D thin-slice models to calculate plaque stress and strain and compare model differences. Our results indicated that the average maximum plaque stress values of 20 patients from multilayer and single-layer models were 385.13 ± 110.09 kPa and 270.91 ± 95.86 kPa, respectively. The relative difference was 42.2%, with single-layer stress serving as the base value. The average mean plaque stress values from multilayer and single-layer models were 129.59 ± 32.77 kPa and 93.27 ± 18.20 kPa, respectively, with a relative difference of 38.9%. The maximum and mean plaque strain values obtained from the multilayer models were 11.6% and 19.0% higher than those from the single-layer models. Similarly, the maximum and mean cap strains showed increases of 9.6% and 12.9% over those from the single-layer models. These findings suggest that use of multilayer models could improve plaque stress and strain calculation accuracy and may have large impact on plaque progression and vulnerability investigation and potential clinical applications. Further large-scale studies are needed to validate our findings.

Impact of cyclic bending on coronary hemodynamics

It remains unknown that the degree of bias in computational fluid dynamics results without considering coronary cyclic bending. This study aims to investigate the influence of different rates of coronary cyclic bending on coronary hemodynamics. To model coronary bending, a multi-ring-controlled fluid-structural interaction model was designed. A coronary artery was simulated with various cyclic bending rates (0.5, 0.75 and 1 s, corresponding to heart rates of 120, 80 and 60 bpm) and compared against a stable model. The simulated results show that the hemodynamic parameters of vortex Q-criterion, temporal wall shear stress (WSS), time-averaged WSS (TaWSS) and oscillatory shear index (OSI) were sensitive to the changes in cyclic rate. A higher heart rate resulted in higher magnitude and larger variance in the hemodynamic parameters. Whereas, the values and distributions of flow velocity and relative residence time (RRT) did not show significant differences between different bending periods. This study suggests that a stable coronary model is not sufficient to represent the hemodynamics in a bending coronary artery. Different heart rate conditions were found to have significant impact on the hemodynamic parameters. Thus, cyclic bending should be considered to mimic the realistic hemodynamics in future patient-specific coronary hemodynamics studies.

How getting twisted in scaffold design can promote bone regeneration: A fluid-structure interaction evaluation

Bone tissue engineering (BTE) uses engineering principles to repair large bone defects, which requires effective mass transport ability of scaffolds to support cellular activities during bone regeneration. Since the implanted BTE scaffolds keep deforming under physiological loading which influences the fluid flow and mass transport within the scaffold and surrounding tissue, thus, scaffold design needs to consider the mass transport behavior under the physiological loading. This work proposed a novel twist scaffold, and its mass transport efficiency under physiological loading conditions was evaluated by a fluid-structure interaction analysis. The results showed that compared to the non-twist scaffold, the twist scaffold could form a rotating flow under the physiological loading, which enhanced the mass transport and generated more appropriate wall shear stress (WSS) to promote bone regeneration. This highlighted the better mass transport efficiency of the twist scaffold. Therefore, getting twist may be a promising design strategy for future BTE scaffolds, and the fluid-structure interaction approach may be a more reliable method for bone regeneration studies in either in vivo or in vitro systems.

Automated finite element approach to generate anatomical patient-specific biomechanical models of atherosclerotic arteries from virtual histology-intravascular ultrasound

Despite advancements in early detection and treatment, atherosclerosis remains the leading cause of death across all cardiovascular diseases (CVD). Biomechanical analysis of atherosclerotic lesions has the potential to reveal biomechanically instable or rupture-prone regions. Treatment decisions rarely consider the biomechanics of the stenosed lesion due in-part to difficulties in obtaining this information in a clinical setting. Previous 3D FEA approaches have incompletely incorporated the complex curvature of arterial geometry, material heterogeneity, and use of patient-specific data. To address these limitations and clinical need, herein we present a user-friendly fully automated program to reconstruct and simulate the wall mechanics of patient-specific atherosclerotic coronary arteries. The program enables 3D reconstruction from patient-specific data with heterogenous tissue assignment and complex arterial curvature. Eleven arteries with coronary artery disease (CAD) underwent baseline and 6-month follow-up angiographic and virtual histology-intravascular ultrasound (VH-IVUS) imaging. VH-IVUS images were processed to remove background noise, extract VH plaque material data, and luminal and outer contours. Angiography data was used to orient the artery profiles along the 3D centerlines. The resulting surface mesh is then resampled for uniformity and tetrahedralized to generate the volumetric mesh using TetGen. A mesh convergence study revealed edge lengths between 0.04 mm and 0.2 mm produced constituent volumes that were largely unchanged, hence, to save computational resources, a value of 0.2 mm was used throughout. Materials are assigned and finite element analysis (FEA) is then performed to determine stresses and strains across the artery wall. In a representative artery, the highest average effective stress was in calcium elements with 235 kPa while necrotic elements had the lowest average stress, reaching as low as 0.79 kPa. After applying nodal smoothening, the maximum effective stress across 11 arteries remained below 288 kPa, implying biomechanically stable plaques. Indeed, all atherosclerotic plaques remained unruptured at the 6-month longitudinal follow up diagnosis. These results suggest our automated analysis may facilitate assessment of atherosclerotic plaque stability.

The Need to Shift from Morphological to Structural Assessment for Carotid Plaque Vulnerability

Degree of luminal stenosis is generally considered to be an important indicator for judging the risk of atherosclerosis burden. However, patients with the same or similar degree of stenosis may have significant differences in plaque morphology and biomechanical factors. This study investigated three patients with carotid atherosclerosis within a similar range of stenosis. Using our developed fluid-structure interaction (FSI) modelling method, this study analyzed and compared the morphological and biomechanical parameters of the three patients. Although their degrees of carotid stenosis were similar, the plaque components showed a significant difference. The distribution range of time-averaged wall shear stress (TAWSS) of patient 2 was wider than that of patient 1 and patient 3. Patient 2 also had a much smaller plaque stress compared to the other two patients. There were significant differences in TAWSS and plaque stresses among three patients. This study suggests that plaque vulnerability is not determined by a single morphological factor, but rather by the combined structure. It is necessary to transform the morphological assessment into a structural assessment of the risk of plaque rupture.

In Vivo Intravascular Optical Coherence Tomography (IVOCT) Structural and Blood Flow Imaging Based Mechanical Simulation Analysis of a Blood Vessel

INTRODUCTION

Computer modelling of blood vessels based on biomedical imaging provides important hemodynamic and biomechanical information for vascular disease studies and diagnosis. However due to lacking well-defined physiological blood flow profiles, the accuracy of the simulation results is often not guaranteed. Flow velocity profiles of a specific cross section of a blood vessel were obtained by in vivo Doppler intravascular optical coherence tomography (IVOCT) lately. However due to the influence of the catheter, the velocity profile imaged by IVOCT can't be applied to simulation directly.

METHODS

A simulation-experiment combined method to determine the inlet flow boundary based on in vivo porcine carotid Doppler IVOCT imaging is proposed. A single conduit carotid model was created from the 3D IVOCT structural images using an image processing-computer aided design combined method.

RESULTS

With both high- resolution arterial model and near physiological blood flow profile, stress analysis by fluid-structure interaction and computational fluid dynamics were performed. The influence of the catheter to the wall shear stress, the hemodynamics and the stresses of the carotid wall under the measured inlet flow and various outlet pressure boundary conditions, are analyzed.

CONCLUSION

This study provides a solution to the difficulty of getting inlet flow boundary for numerical simulation of arteries. It paves the way for developing IVOCT based vascular stress analysis and imaging methods for the studies and diagnosis of vascular diseases.

A Study of the Fluid–Structure Interaction of the Plaque Circumferential Distribution in the Left Coronary Artery

Atherosclerotic plaques within the coronary arteries can prevent blood from flowing to downstream tissues, causing coronary heart disease and a myocardial infarction over time. The degree of stenosis is an important reference point during percutaneous coronary intervention (PCI). However, clinically, patients with the same degree of stenosis exhibit different degrees of disease severity. To investigate the connection between this phenomenon and the plaque circumferential distribution, in this paper, four models with different plaque circumferential locations were made based on the CT data. The blood in the coronary arteries was simulated using the fluid–structure interaction method in ANSYS Workbench software. The results showed that the risk of plaque rupture was less affected by the circumferential distribution of plaque, and the distribution of blood in each branch was affected by the circumferential distribution of plaque. Low TAWSS areas were found posterior to the plaque, and the TAWSS < 0.4 Pa area was ranked from highest to lowest in each model species: plaque on the side away from the left circumflex branch, plaque on the side away from the heart; plaque on the side close to the heart; and plaque on the side close to the left circumflex branch. The same trend was also found in the OSI. It was concluded that the circumferential distribution of plaques affects their further development. This finding will be useful for clinical treatment.

Automated Coronary Optical Coherence Tomography Feature Extraction with Application to Three-Dimensional Reconstruction

Coronary optical coherence tomography (OCT) is an intravascular, near-infrared light-based imaging modality capable of reaching axial resolutions of 10–20 µm. This resolution allows for accurate determination of high-risk plaque features, such as thin cap fibroatheroma; however, visualization of morphological features alone still provides unreliable positive predictive capability for plaque progression or future major adverse cardiovascular events (MACE). Biomechanical simulation could assist in this prediction, but this requires extracting morphological features from intravascular imaging to construct accurate three-dimensional (3D) simulations of patients’ arteries. Extracting these features is a laborious process, often carried out manually by trained experts. To address this challenge, numerous techniques have emerged to automate these processes while simultaneously overcoming difficulties associated with OCT imaging, such as its limited penetration depth. This systematic review summarizes advances in automated segmentation techniques from the past five years (2016–2021) with a focus on their application to the 3D reconstruction of vessels and their subsequent simulation. We discuss four categories based on the feature being processed, namely: coronary lumen; artery layers; plaque characteristics and subtypes; and stents. Areas for future innovation are also discussed as well as their potential for future translation.

Accurate Calculation of FFR Based on a Physics-Driven Fluid‐Structure Interaction Model

Background: The conventional FFRct numerical calculation method uses a model with a multi-scale geometry based upon CFD, and rigid walls. Therefore, important interactions between the elastic vessel wall and blood flow are not routinely considered. Changes in the resistance of coronary microcirculation during hyperaemia are likewise not typically incorporated using a fluid–structure interaction (FSI) algorithm. It is likely that both have resulted in FFRct calculation errors.

Objective: In this study we incorporated both the influence of vascular elasticity and coronary microcirculatory structure on FFR, to improve the accuracy of FFRct calculation. Thus, in this study, a physics-driven 3D–0D coupled model including fluid–structure interaction was established to calculate accurate FFRct values.

Methods: Based upon a novel geometric multi-scale modeling technology, a FSI simulation approach was used. A lumped parameter model (0D) was used as the outlet boundary condition for the 3D FSI coronary artery model to incorporate physiological microcirculation, with bidirectional coupling between the two models.

Results: The accuracy, sensitivity, specificity, and both positive and negative predictive values of FFR_DC calculated based upon the coupled 3D–0D model were 86.7, 66.7, 84.6, 66.7, and 91.7%, respectively. Compared to the calculated value using the basic CFD model (MSE = 5.9%, accuracy rate = 80%), the FFR_CFD calculated based on the coupled 3D–0D model has a smaller MSE of 1.9%.

Conclusion: The physics-driven coupled 3D–0D model that incorporates fluid–structure interactions not only consider the influence of the elastic vessel wall on blood flow, but also provides reliable microvascular resistance boundary conditions for the 3D FSI model. This allows for a calculation that is based upon conditions that are closer to the physiological environment, and thus improves the accuracy of FFRct calculation. It is likely that more accurate information will provide an enhanced recommendation regarding percutaneous coronary intervention (PCI) in the clinic.

Automatic Detection of Age-Related Macular Degeneration Based on Deep Learning and Local Outlier Factor Algorithm

Age-related macular degeneration (AMD) is a retinal disorder affecting the elderly, and society’s aging population means that the disease is becoming increasingly prevalent. The vision in patients with early AMD is usually unaffected or nearly normal but central vision may be weakened or even lost if timely treatment is not performed. Therefore, early diagnosis is particularly important to prevent the further exacerbation of AMD. This paper proposed a novel automatic detection method of AMD from optical coherence tomography (OCT) images based on deep learning and a local outlier factor (LOF) algorithm. A ResNet-50 model with L2-constrained softmax loss was retrained to extract features from OCT images and the LOF algorithm was used as the classifier. The proposed method was trained on the UCSD dataset and tested on both the UCSD dataset and Duke dataset, with an accuracy of 99.87% and 97.56%, respectively. Even though the model was only trained on the UCSD dataset, it obtained good detection accuracy when tested on another dataset. Comparison with other methods also indicates the efficiency of the proposed method in detecting AMD.

Image-based biomechanical modeling for coronary atherosclerotic plaque progression and vulnerability prediction

Atherosclerotic plaque progression and rupture play an important role in cardiovascular disease development and the final drastic events such as heart attack and stroke. Medical imaging and image-based computational modeling methods advanced considerably in recent years to quantify plaque morphology and biomechanical conditions and gain a better understanding of plaque evolution and rupture process. This article first briefly reviewed clinical imaging techniques for coronary thin-cap fibroatheroma (TCFA) plaques used in image-based computational modeling. This was followed by a summary of different types of biomechanical models for coronary plaques. Plaque progression and vulnerability prediction studies based on image-based computational modeling were reviewed and compared. Much progress has been made and a reasonable high prediction accuracy has been achieved. However, there are still some inconsistencies in existing literature on the impact of biomechanical and morphological factors on future plaque behavior, and it is very difficult to perform direct comparison analysis as differences like image modality, biomechanical factors selection, predictive models, and progression/vulnerability measures exist among these studies. Encouraging data and model sharing across the research community would partially resolve these differences, and possibly lead to clearer assertive conclusions. In vivo image-based computational modeling could be used as a powerful tool for quantitative assessment of coronary plaque vulnerability for potential clinical applications.

Mathematical modeling of plaque progression and associated microenvironment: How far from predicting the fate of atherosclerosis?

Mathematical modeling contributes to pathophysiological research of atherosclerosis by helping to elucidate mechanisms and by providing quantitative predictions that can be validated. In turn, the complexity of atherosclerosis is well suited to quantitative approaches as it provides challenges and opportunities for new developments of modeling. In this review, we summarize the current 'state of the art' on the mathematical modeling of the effects of biomechanical factors and microenvironmental factors on the plaque progression, and its potential help in prediction of plaque development. We begin with models that describe the biomechanical environment inside and outside the plaque and its influence on its growth and rupture. We then discuss mathematical models that describe the dynamic evolution of plaque microenvironmental factors, such as lipid deposition, inflammation, smooth muscle cells migration and intraplaque hemorrhage, followed by studies on plaque growth and progression using these modelling approaches. Moreover, we present several key questions for future research. Mathematical models can complement experimental and clinical studies, but also challenge current paradigms, redefine our understanding of mechanisms driving plaque vulnerability and propose future potential direction in therapy for cardiovascular disease.

A platform for high-fidelity patient-specific structural modelling of atherosclerotic arteries: from intravascular imaging to three-dimensional stress distributions

The pathophysiology of atherosclerotic lesions, including plaque rupture triggered by mechanical failure of the vessel wall, depends directly on the plaque morphology-modulated mechanical response. The complex interplay between lesion morphology and structural behaviour can be studied with high-fidelity computational modelling. However, construction of three-dimensional (3D) and heterogeneous models is challenging, with most previous work focusing on two-dimensional geometries or on single-material lesion compositions. Addressing these limitations, we here present a semi-automatic computational platform, leveraging clinical optical coherence tomography images to effectively reconstruct a 3D patient-specific multi-material model of atherosclerotic plaques, for which the mechanical response is obtained by structural finite-element simulations. To demonstrate the importance of including multi-material plaque components when recovering the mechanical response, a computational case study was conducted in which systematic variation of the intraplaque lipid and calcium was performed. The study demonstrated that the inclusion of various tissue components greatly affected the lesion mechanical response, illustrating the importance of multi-material formulations. This platform accordingly provides a viable foundation for studying how plaque micro-morphology affects plaque mechanical response, allowing for patient-specific assessments and extension into clinically relevant patient cohorts.

A Review on the Effect of Temporal Geometric Variations of the Coronary Arteries on the Wall Shear Stress and Pressure Drop

Temporal variations of the coronary arteries during a cardiac cycle are defined as the superposition of the changes in the position, curvature, and torsion of the coronary artery axis markers and the variations in the lumen cross-sectional shape due to the distensible wall motion induced by the pulse pressure and contraction of the myocardium in a cardiac cycle. This review discusses whether modeling of the temporal variations of the coronary arteries is needed for the investigation of hemodynamics specifically in time-critical applications such as a clinical environment. The numerical modelings in the literature that model or disregard the temporal variations of the coronary arteries on the hemodynamic parameters are discussed. The results in the literature show that neglecting the effects of temporal geometric variations is expected to result in about 5% deviation of the time-averaged pressure drop and wall shear stress values and also about 20% deviation of the temporal variations of hemodynamic parameters, such as time-dependent wall shear stress and oscillatory shear index. This review study can be considered as a guide for future studies to outline the conditions in which temporal variations of the coronary arteries can be neglected while providing a reliable estimation of hemodynamic parameters.

Impact of stent malapposition on intracoronary flow dynamics: An optical coherence tomography-based patient-specific study

Percutaneous coronary intervention with stent implantation has emerged as a popular approach to treat coronary artery stenosis. Stent malapposition (SM), also referred as incomplete stent apposition, could reduce stent tissue coverage and hence increase the risk of late stent thrombosis. The objective of this study was to investigate the impact of SM on intracoronary flow dynamics by combining optical coherence tomography (OCT) image-based model reconstruction and computational analysis. Firstly, a stenosed coronary artery model was reconstructed from OCT and angiography imaging data of a patient. Two structural analyses were carried out to simulate two types of coronary artery stent implantations: a fully-apposed (FA) case and a SM case. Then, based on the two deformed coronary geometries, two computational fluid dynamics (CFD) analyses were performed to evaluate the differences of hemodynamic metrics between the FA and the SM cases, including wall shear stress (WSS), time-averaged WSS (TWSS), oscillatory shear index (OSI), WSS gradient (WSSG), time-averaged WSSG (TWSSG), and relative residence time (RRT). The results indicated that maximum flow velocity was higher in the SM case than that of the FA case, due to the incomplete expansion of the stent and artery. Moreover, the SM case had a lower percentage of areas of adverse WSS (< 0.5 Pa) and RRT (> 10/Pa) but a higher percentage of areas of adverse OSI (> 0.1) and WSSG (> 5000 Pa/m). Specifically, the differences of OSI, WSSG, and RRT between the two cases were relatively small. It was suggested that SM might not be responsible for negative hemodynamic metrics which would further result in stent thrombosis on the basis of the present specific model.

Impact of left atrial appendage location on risk of thrombus formation in patients with atrial fibrillation

Most strokes in patients with atrial fibrillation (AF) are thought to arise from thrombus formation in the left atrial appendage (LAA). Assessing the hemodynamics in LAA and left atrium (LA) may provide some insights in the evaluation of the risk of thrombus formation. This study aims to find out the impact of different LAA locations with respect of LA on the risk of thrombus formation within LAA in patients with AF. Three different LAA locations at LA were modeled and a fully coupled fluid-structure interaction analysis was performed. A discrete phase method was used for particle residence analysis to evaluate risk of the thrombus formation. The results showed that LAA positions on the LA affected the LAA flow velocity distribution, passive contraction ability, and particle residence. In particular, the left pulmonary veins (PVs) had a greater influence on the LAA hemodynamics than the right PVs. The LAA had the lowest contractibility when it was located between left superior and left inferior PVs, and in this case, a larger number of particles were resided, which indicated a higher risk of thrombus formation. The present work provides a quantitative way to evaluate the risk of thrombus formation within LAA in patients with AF.

Case Report: Evaluating Biomechanical Risk Factors in Carotid Stenosis by Patient-Specific Fluid-Structural Interaction Biomechanical Analysis

BACKGROUND

Carotid atherosclerosis is one of the main underlying inducements of stroke, which is a leading cause of disability. The morphological feature and biomechanical environment have been found to play important roles in atherosclerotic plaque progression. However, the biomechanics in each patient's blood vessel is complicated and unique.

METHOD

To analyse the biomechanical risk of the patient-specific carotid stenosis, this study used the fluid-structure interaction (FSI) computational biomechanical model. This model coupled both structural and hemodynamic analysis. Two patients with carotid stenosis planned for carotid endarterectomy were included in this study. The 3D models of carotid bifurcation were reconstructed using our in-house-developed protocol based on multisequence magnetic resonance imaging (MRI) data. Patient-specific flow and pressure waveforms were used in the computational analysis. Multiple biomechanical risk factors including structural and hemodynamic stresses were employed in post-processing to assess the plaque vulnerability.

RESULTS

Significant difference in morphological and biomechanical conditions between 2 patients was observed. Patient I had a large lipid core and serve stenosis at carotid bulb. The stenosis changed the cross-sectional shape of the lumen. The blood flow pattern changed consequently and led to a complex biomechanical environment. The FSI results suggested a potential plaque progression may lead to a high-risk plaque, if no proper treatment was performed. The patient II had significant tandem stenosis at both common and internal carotid artery (CCA and ICA). From the results of biomechanical factors, both stenoses had a high potential of plaque progression. Especially for the plaque at ICA branch, the current 2 small plaques might further enlarge and merge as a large vulnerable plaque. The risk of plaque rupture would also increase.

CONCLUSIONS

Computational biomechanical analysis is a useful tool to provide the biomechanical risk factors to help clinicians assess and predict the patient-specific plaque vulnerability. The FSI computational model coupling the structural and hemodynamic computational analysis, better replicates the in vivo biomechanical condition, which can provide multiple structural and flow-based risk factors to assess plaque vulnerability.

Application of physics-based flow models in cardiovascular medicine: Current practices and challenges

Personalized physics-based flow models are becoming increasingly important in cardiovascular medicine. They are a powerful complement to traditional methods of clinical decision-making and offer a wealth of physiological information beyond conventional anatomic viewing using medical imaging data. These models have been used to identify key hemodynamic biomarkers, such as pressure gradient and wall shear stress, which are associated with determining the functional severity of cardiovascular diseases. Importantly, simulation-driven diagnostics can help researchers understand the complex interplay between geometric and fluid dynamic parameters, which can ultimately improve patient outcomes and treatment planning. The possibility to compute and predict diagnostic variables and hemodynamics biomarkers can therefore play a pivotal role in reducing adverse treatment outcomes and accelerate development of novel strategies for cardiovascular disease management.

Computational fluid dynamics simulation of time-resolved blood flow in Budd-Chiari syndrome with inferior vena cava stenosis and its implication for postoperative efficacy assessment

BACKGROUND

This study aimed to adopt computational fluid dynamics to simulate the blood flow dynamics in inferior vena cava stenosis based on time-dependent patient-specific models of Budd-Chiari syndrome as well as a normal model. It could offer valuable references for a retrospective insight into the underlying mechanisms of Budd-Chiari syndrome pathogenesis as well as more accurate evaluation of postoperative efficacy.

METHODS

Three-dimensional inferior vena cava models of Budd-Chiari syndrome patient-specific (preoperative and postoperative) and normal morphology model were reconstructed as per magnetic resonance images using Simpleware. Moreover, computational fluid dynamics of time-resolved inferior vena cava blood flow were simulated using actual patient-specific measurements to reflect time-dependent flow rates.

FINDINGS

The assessment of the preoperative model revealed the dramatic variations of hemodynamic parameters of the stenotic inferior vena cava. Moreover, the comparison of the preoperative and postoperative models with the normal model as benchmark showed that postoperative hemodynamic parameters were markedly ameliorated via stenting, with the attenuation of overall velocity and wall shear stress, and the increase of pressure. However, the comparative analysis of the patient-specific simulations revealed that some postoperative hemodynamic profiles still bore some resemblance to the preoperative ones, indicating potential risks of restenosis.

INTERPRETATION

Computational fluid dynamics simulation of time-resolved blood flow could reveal the tight correlation between the hemodynamic characteristics and the pathological mechanisms of inferior vena cava stenosis. Furthermore, such time-resolved hemodynamic profiles could provide a quantitative approach to diagnosis, operative regimen and postoperative evaluation of Budd-Chiari syndrome with inferior vena cava stenosis.

Evaluating the Impact of Calcification on Plaque Vulnerability from the Aspect of Mechanical Interaction Between Blood Flow and Artery Based on MRI

Acute cerebral ischemic events and thrombosis are associated with the rupture/erosion of carotid atherosclerotic plaques. The aim of the present study was to determine the impact of calcification deposition on the wall shear stress (WSS) and stresses within the plaques using 3D fluid-structure interaction (FSI) models. Six patients with calcified carotid atherosclerosis underwent multisequence magnetic resonance imaging (MRI) and were divided into three groups according to the calcification volume. To evaluate the role of the calcification deposition on the stresses, the calcification content was replaced by lipids and arterial tissue, respectively. By comparing the results from the simulation with calcification, and when changing it to lipids there was a significant increment in the stresses at the fibrous cap (p = 0.004). Instead, by changing it to arterial tissue, there was no significant difference (p = 0.07). The calcification shapes that presented the highest stresses were thin concave arc-shaped (AS1) and thin convex arc-shaped (AS3), with mean stress values of 107 ± 54.2 and 99.6 ± 23.4 kPa, respectively. It was also observed that, the calcification shape has more influence on the level of stress than its distance to the lumen. Higher WSS values were associated with the presence of calcification. Calcification shape plays an important role in producing high stresses in the plaque. This work further clarifies the impact of calcification on plaque vulnerability.

Stress-Relaxation and Cyclic Behavior of Human Carotid Plaque Tissue

Atherosclerotic plaque rupture is a catastrophic event that contributes to mortality and long-term disability. A better understanding of the plaque mechanical behavior is essential for the identification of vulnerable plaques pre-rupture. Plaque is subjected to a natural dynamic mechanical environment under hemodynamic loading. Therefore, it is important to understand the mechanical response of plaque tissue under cyclic loading conditions. Moreover, experimental data of such mechanical properties are fundamental for more clinically relevant biomechanical modeling and numerical simulations for risk stratification. This study aims to experimentally and numerically characterize the stress-relaxation and cyclic mechanical behavior of carotid plaque tissue. Instron microtester equipped with a custom-developed setup was used for the experiments. Carotid plaque samples excised at endarterectomy were subjected to uniaxial tensile, stress-relaxation, and cyclic loading protocols. Thirty percent of the underlying load level obtained from the uniaxial tensile test results was used to determine the change in mechanical properties of the tissue over time under a controlled testing environment (Control tests). The stress-relaxation test data was used to calibrate the hyperelastic (neo-Hookean, Ogden, Yeoh) and linear viscoelastic (Prony series) material parameters. The normalized relaxation force increased initially and slowly stabilized toward the end of relaxation phase, highlighting the viscoelastic behavior. During the cyclic tests, there was a decrease in the peak force as a function of the cycle number indicating mechanical distension due to repeated loading that varied with different frequencies. The material also accumulated residual deformation, which increased with the cycle number. This trend showed softening behavior of the samples. The results of this preliminary study provide an enhanced understanding of in vivo stress-relaxation and cyclic behavior of the human atherosclerotic plaque tissue.

Carotid Bifurcation With Tandem Stenosis-A Patient-Specific Case Study Combined in vivo Imaging, in vitro Histology and in silico Simulation

A patient-specific carotid bifurcation with tandem stenosis found at both internal carotid artery (ICA) and common carotid artery (CCA) was studied. The in vivo pre-carotid endarterectomy (pre-CEA) multi-spectral magnetic resonance imaging (MRI) were performed and in vitro post-CEA carotid plaque tissue sample was collected. MR imaging data and tissue sample staining histology were used to recognize the plaque components. Further, the computational fluid dynamics (CFD) were performed on four MR-based reconstructed 3D carotid bifurcation models (the patient-specific geometry with tandem stenosis and three presumptive geometries by removing the stenosis part). The flow and shear stress behavior affected by the tandem stenosis was analyzed. From the results of MR segmentation and histology analysis, plaque lipid pool and calcification were found at both ICA and CCA. From the result of CFD simulation, the flow shear stress behavior suggested the tandem stenosis as a more “dangerous” situation than a single-stenosis artery. Besides, the CFD results deduced that the stenosis at the CCA location formed initially and led to the subsequent formation of stenosis at ICA. This study suggests that when planning CEA, CFD simulation on the presumptive models could help clinicians to estimate the blood flow behavior after surgery. Particular attention should be paid to the case of tandem stenosis, as the local hemodynamic environment is more complex and treatment of one stenosis may lead to a variation in the hemodynamic loading on the second plaque, which may result in either a higher risk of plaque rupture or restenosis.