Subject-Specific Fully-Coupled and One-Way Fluid-Structure Interaction Models for Modeling of Carotid Atherosclerotic Plaques in Humans

Image-based insilico investigation of hemodynamics and biomechanics in healthy and diabetic human retinas

Retinal hemodynamics and biomechanics play a significant role in understanding the pathophysiology of several ocular diseases. However, these parameters are significantly affected due to changed blood vessel morphology ascribed to pathological conditions, particularly diabetes. In this study, an image-based computational fluid dynamics (CFD) model is applied to examine the effects of changed vascular morphology due to diabetes on blood flow velocity, vorticity, wall shear stress (WSS), and oxygen distribution and compare it with healthy. The 3D patient-specific vascular architecture of diabetic and healthy retina is extracted from Optical Coherence Tomography Angiography (OCTA) images and fundus to extract the capillary level information. Further, Fluid-structure interaction (FSI) simulations have been performed to compare the induced tissue stresses in diabetic and healthy conditions. Results illustrate that most arterioles possess higher velocity, vorticity, WSS, and lesser oxygen concentration than arteries for healthy and diabetic cases. However, an opposite trend is observed for venules and veins. Comparisons show that, on average, the blood flow velocity in the healthy case decreases by 42 % in arteries and 21 % in veins, respectively, compared to diabetic. In addition, the WSS and von Mises stress (VMS) in healthy case decrease by 49 % and 72 % in arteries and by 6 % and 28 % in veins, respectively, when compared with diabetic, making diabetic blood vessels more susceptible to wall rupture and tissue damage. The in-silico results may help predict the possible abnormalities region early, helping the ophthalmologists use these estimates as prognostic tools and tailor patient-specific treatment plans.

Comparative study between 1-way and 2-way coupled fluid-structure interaction in numerical simulation of aortic arch aneurysms

Hemodynamic forces are related to pathological variations of the cardiovascular system, and numerical simulations for fluid-structure interaction have been systematically used to analyze the behavior of blood flow and the arterial wall in aortic aneurysms. This paper proposes a comparative analysis of 1-way and 2-way coupled fluid-structure interaction for aortic arch aneurysm. The coupling models of fluid-structure interaction were conducted using 3D geometry of the thoracic aorta from computed tomography. Hyperelastic anisotropic properties were estimated for the Holzapfel arterial wall model. The rheological behavior of the blood was modeled by the Carreau-Yasuda model. The results showed that the 1-way approach tends to underestimate von Mises stress, displacement, and strain over the entire cardiac cycle, compared to the 2-way approach. In contrast, the behavior of the variables of flow field, velocity, wall shear stress, and Reynolds number when coupled by the 1-way model was overestimated at the systolic moment and tends to be equal at the diastolic moment. The quantitative differences found, especially during the systole, suggest the use of 2-way coupling in numerical simulations of aortic arch aneurysms due to the hyperelastic nature of the arterial wall, which leads to a strong iteration between the fluid and the arterial wall.

Fluid structure interaction analysis to stratify the behavior of different atheromatous carotid plaques

BACKGROUND

In asymptomatic carotid stenosis (ACS), different plaque types, i.e. lipidic (LP), fibrous (FP), and calcific (CP), could have different hemodynamic and structural behaviors.

METHODS

Different carotid plaques, reconstructed from medical imaging of ACS >70%, were analyzed by computing fluid structure interaction (FSI), modeling the spatial distribution of wall shear stresses (WSS), plaque displacements (D), von Mises stresses (VMS), and absorbed elastic energy (AEE) together with their maximum-in-space values at the systole (WSS<inf>syst</inf>, D<inf>syst</inf>, VMS<inf>syst</inf> and AEE<inf>syst</inf>).

RESULTS

WSS resulted significantly higher in CP, whereas D and VMS showed the highest values for LP. Regarding AEE<inf>syst</inf> stored by the plaques, LP absorbed in average 2320 J/m³, FP 408 J/m³ (470%) and CP 99 J/m³ (2240%), (P<0.01, P<0.01, and P<0.01, respectively).

CONCLUSIONS

Depending upon their nature, plaques store different deformations and inner distributions of forces, thus potentially influencing their vulnerability.

Medial Gap: A Structural Factor at the Arterial Bifurcation Aggravating Hemodynamic Insult

Previous studies have reported that intracranial aneurysms frequently occur adjacent to the medial gap. However, the role of the medial gap in aneurysm formation is controversial. We designed this study to explore the potential role of the medial gap in aneurysm formation. Widened artery bifurcations with or without medial gaps were microsurgically created and pathologically stained in the carotid arteries of 30 rats. Numerical artery bifurcation models were constructed, and bidirectional fluid-solid interaction analyses were performed. Animal experiments showed that the apexes of widened bifurcations with a medial gap were prone to being insulted by blood flow compared to those without a medial gap. The bidirectional fluid-solid interaction analyses indicated that artery bifurcations with the medial gap exhibited higher wall shear stress (WSS) and von Mises stress (VMS) at the apex of the bifurcation. The disparity of stress between the gap and no-gap model was larger for widened bifurcations, peaking at 180° with a maximum of 1.9 folds. The maximum VMS and relatively high WSS were located at the junction between the medial gap and the adjacent arterial wall. Our results suggest that the medial gap at the widened arterial bifurcation may promote aneurysm formation.

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.

Patient-specific fluid-structure interaction model of bile flow: comparison between 1-way and 2-way algorithms

Gallbladder disease is one of the most spread pathologies in the world. Despite the number of operations dealing with biliary surgery increases, the number of postoperative complications is also high. The aim of this study is to show the influence of the biliary system pathology on bile flow character and to numerically assess the effect of surgical operation (cholecystectomy) on the fluid dynamics in the extrahepatic biliary tree, and also to reveal the difference between 1-way and 2-way FSI algorithms on the results. Moreover, the bile viscosity and biliary tree geometry influence on the choledynamics were evaluated. Bile velocity, pressure, wall shear stress (WSS), displacements and von Mises stress distributions in the extrahepatic biliary tree are presented, and comparison is made between a healthy and a lithogenic bile. The patient-specific biliary tree model is created using magnetic resonance imaging (MRI) and imported in a commercial finite element analysis software. It is found that in the case of lithogenic bile, velocities have lower magnitude while pressures are higher. Furthermore, stress analysis of the bile ducts shows that the WSS distribution is found mostly prevailing in the common hepatic duct and common bile duct areas. It is shown that when it is necessary to evaluate the bile flow dynamics in urgent medical situations, 1-way analysis is acceptable. Nevertheless, 2-way FSI provides more accurate data, if necessary to evaluate the stress-strain state of bile ducts. The proposed model can be applied to medical practice to reduce the number of post-operative complications.

A Computational Fluid-Structure Interaction Study for Carotids With Different Atherosclerotic Plaques

Atherosclerosis is a systemic disease that leads to accumulation of deposits, known as atherosclerotic plaques, within the walls of the carotids. In particular, three types of plaque can be distinguished: soft, fibrous, and calcific. Most of the computational studies who investigated the interplay between the plaque and the blood flow on patient-specific geometries used nonstandard medical images to directly delineate and segment the plaque and its components. However, these techniques are not so widely available in the clinical practice. In this context, the aim of our work was twofold: (i) to propose a new geometric tool that allowed to reconstruct a plausible plaque in the carotids from standard images and (ii) to perform three-dimensional (3D) fluid-structure interaction (FSI) simulations where we compared some fluid-dynamic and structural quantities among 15 patients characterized by different typologies of plaque. Our results highlighted that both the morphology and the mechanical properties of different plaque components play a crucial role in determining the vulnerability of the plaque.

A surrogate model for plaque modeling in carotids based on Robin conditions calibrated by cine MRI data

We propose a surrogate model for the fluid-structure interaction (FSI) problem for the study of blood dynamics in carotid arteries in presence of plaque. This is based on the integration of a numerical model with subject-specific data and clinical imaging. We propose to model the plaque as part of the tissues surrounding the vessel wall through the application of an elastic support boundary condition. In order to characterize the plaque and other surrounding tissues, such as the close-by jugular vein, the elastic parameters of the boundary condition were spatially differentiated and their values were estimated by minimizing the discrepancies between computed vessel displacements and reference values obtained from CINE Magnetic Resonance Imaging data. We applied the model to three subjects with a degree of stenosis greater than 70%. We found that accounting for both plaque and jugular vein in the estimation of the elastic parameters increases the accuracy. In particular, in all patients, mismatches between computed and in vivo measured wall displacements were one to two orders of magnitude lower than the spatial resolution of the original MRI data. These results confirmed the validity of the proposed surrogate plaque model. We also compared fluid-dynamics results with those obtained in a fixed wall setting and in a full FSI model, used as gold standard, highlighting the better accordance of our results in comparison to the rigid ones.

Influence of coupled hemodynamics-arterial wall interaction on compliance in a realistic pulmonary artery with variable intravascular wall properties

Pulmonary hypertension is characterized by elevation of pulmonary artery (PA) pressure (p) and structural remodeling of the PA wall, leading to reduction in arterial compliance (c). As a step towards improving diagnosis of pulmonary disease, we use the PA branch geometry (main pulmonary artery (MPA) branching into left (LPA) and right (RPA) pulmonary arteries) obtained from MRI in conjunction with an inverse algorithm to obtain the pre-stress level in the artery walls. Next, a coupled blood-wall interaction (BWI) calculation provides hemodynamic information as well as compliance of the PA walls. We show that the computed load-free geometry from the inverse algorithm exhibits a 27.8% lower inner diameter (d) and 18.5% lower outer d compared to the in vivo geometry from MRI. Further, the mean p computed from the BWI computation in the main PA (p_MPA-n) is within 4% of the mean p_MPA-e (n-numerical; e-experimental). Also, the mean Q computed in the left PA (Q_LPA-n) is within 10% of the mean Q_LPA-e. Finally, the compliance c_MPA-n is computed to be 27% lower than c_MPA-e, while the c_LPA-n is computed to be 20.4% lower than c_LPA-e. Importantly, the PA shows significant intra-vascular variation in compliance, with the MPA showing higher overall compliance compared to the LPA (3.5-4 times).