Wall Shear Stress Gradient Research Articles

4D-Flow MRI and Vector Ultrasound in the In-Vitro Evaluation of Surgical Aortic Heart Valves - a Pilot Study.

The aim of this study was the initial investigation of 4D-Flow MRI and Vector Ultrasound as novel imaging techniques in the in-vitro analysis of hemodynamics in anatomical models. Specifically, by looking at the hemodynamic performance of state-of-the-art surgical heart valves in a 3D-printed aortic arch. The mock circulatory loop simulated physiological, pulsatile flow. Two mechanical and three biological aortic valves prostheses were compared in a 3D-printed aortic arch. 4D magnetic resonance imaging and vector flow Doppler ultrasound served as imaging methods. Hemodynamic parameters such as wall shear stress, flow velocities and pressure gradients were analyzed. The flow analysis revealed characteristic flow-patterns in the 3D-printed aortic arch. The blood-flow in the arch presented complex patterns, including the formation of helixes and vortices. Higher proximal peak velocities and lower flow volumes were found for biological valves. The mock circulatory loop in combination with modern radiological imaging provides a sufficient basis for the hemodynamic comparison of aortic valves.

Can β-blockers prevent intracranial aneurysm rupture? - insights from Computational Fluid Dynamics analysis.

Hypertension is a risk factor for intracranial aneurysm rupture. We analyzed whether the intake of drugs from specific classes of anti-hypertensive medications affects hemodynamic parameters of intracranial aneurysm dome. We recorded medical history including medications and the in-hospital blood pressure values. We then obtained 3D reconstruction of each patients' aneurysm dome and the feeding artery. Using OpenFOAM software we performed Computational Fluid Dynamics analysis of blood flow through the modeled structures. Blood was modeled as Newtonian fluid, using the incompressible transient solver. As the inlet boundary condition we used the patient-specific Internal Carotid Artery blood velocity waves obtained with Doppler ultrasound. We calculated haemodynamic parameters of the aneurysm dome. All presented analyses are cross-sectional.We included 72 patients with a total of 91 unruptured intracranial aneurysms. The history of β-blocker intake significantly influenced hemodynamic parameters of aneurysm dome. The patients on β-blockers had significantly smaller aneurysm domes (5.09 ± 2.11 mm vs. 7.41 ± 5.89 mm; p = 0.03) and did not have aneurysms larger than 10 mm (0% vs 17.0%; p = 0.01). In the Computational Fluid Dynamics analysis, walls of aneurysms in patients who took β-blockers were characterized by lower Wall Shear Stress Gradient (1.67 ± 1.85 Pa vs. 4.3 ± 6.06 Pa; p = 0.03), Oscillatory Shear Index (0.03 ± 0.02 vs. 0.07 ± 0.10; p = 0.04) and Surface Vortex Fraction (16.2% ± 5.2% vs. 20.0% ± 6.8%; p<0.01). After controlling for covariates, we demonstrated difference of Surface Vortex Fraction (F[1, 48] = 4.36; p = 0.04) and Oscillatory Shear Index (F[1, 48] = 6.51; p = 0.01) between patients taking and not taking β-blockers, respectively. Intake of β-blockers might contribute to more favorable hemodynamics inside aneurysmal sac. Other antihypertensive medication classes were not associated with differences in intracranial aneurysm parameters.

Numerical investigation on flow-induced wall shear stress variation of metastatic cancer cells in lymphatics with elastic valves

Hematogenous metastasis occurs when cancer cells detach from the extracellular matrix in the primary tumor into the bloodstream or lymphatic system. Elucidating the response of metastatic tumor cells in suspension to the flow conditions in lymphatics with valves from a mechanical/fluidic perspective is necessary. A physiologically relevant computational model of a lymphatic vessel with valves was constructed using fully coupled fluid–cell–vessel interactions to investigate the effects of lymphatic vessel contractility, valve properties, and cell size and stiffness on the variations in magnitude and gradient of the flow-induced wall shear stress (WSS) experienced by suspended tumor cells. Results indicated that the maximum WSSmax increased with the increments in cell diameter, vessel contraction amplitude, and valve stiffness. The decrease in vessel contraction period and valve aspect ratio also increased the maximum WSSmax. The influence of the properties of the valve on the WSS was more significant among the factors mentioned above. The maximum WSSmax acting on the cancer cell when the cell reversed the direction of its motion in the valve region increased by 0.5–1.4 times that before the cell entered the valve region. The maximum change in WSS was in the range of 0.004–0.028 Pa/µm depending on the factors studied. They slightly exceeded the values associated with breast cancer cell apoptosis. The results of this study provide biofluid mechanics-based support for mechanobiological research on the metastasis of metastatic cancer cells in suspension within the lymphatics.

Comparative analysis of mechanical wall shear stress and hemodynamics to study the influence of asymmetry in abdominal aortic aneurysm and descending thoracic aortic aneurysm

An aneurysm's rupture is commonly associated with its maximum diameter, yet biomechanical studies emphasize the significant influence of mechanical wall shear stress (WSS) in this process. This study focuses on two models of aortic aneurysms: abdominal aortic aneurysm and descending thoracic aortic aneurysm. Five cases, comprising two for model 1 and three for model 2, are examined to explore both axisymmetric and asymmetric shapes, as patient geometry may manifest as either fusiform (axisymmetric) or saccular (asymmetric), while maintaining a consistent aneurysm diameter and adjusting the bulge shape factor to induce asymmetry. Hemodynamic factors, including WSS and wall shear stress gradient, are computed to evaluate thrombus formation and rupture risk within the aneurysms. Our results indicate the presence of recirculation zones in both the medial and transverse planes, generating vortices within the aneurysm. These vortices are more prominent in asymmetric cases compared to axisymmetric cases, leading to increased blood residence time within the aneurysm and a higher likelihood of thrombus formation. Thrombus formation can further impede blood flow, heightening the risk of embolism or ischemic events. Rupture occurs when the WSS surpasses tissue strength; thus, if the tissue strength of all aneurysms is same, our findings suggest that rupture risk varies according to asymmetry. In the transverse direction, our results demonstrate that in model 1, case 1 exhibits uniform WSS on both sides, while in case 2, WSS is higher at the posterior sides of the aneurysm sac. Conversely, in model 2, WSS is higher at the anterior side of the aneurysm. In the medial direction of the aneurysm, WSS is highest for case 5, followed by case 3, case 4, case 2, and case 1, respectively, indicating elevated WSS when the anterior bulge dominates over the posterior bulge for each model. Overall, a higher rupture risk is observed in model 2 compared to model 1 due to increased mechanical stresses.

Identifying hemodynamic factors associated with the rupture of anterior communicating artery aneurysms based on global modeling of blood flow in the cerebral artery network

Anterior communicating artery (ACoA) aneurysms are more prone to rupture compared to aneurysms present in other cerebral arteries. We hypothesize that systemic blood flow in the cerebral artery network plays an important role in shaping intra-aneurysmal hemodynamic environment thereby affecting the rupture risk of ACoA aneurysms. The majority of existing numerical studies in this field employed local modeling methods where the physical boundaries of a model are confined to the aneurysm region, which, though having the benefit of reducing computational cost, may compromise the physiological fidelity of numerical results due to insufficient account of systemic cerebral arterial hemodynamics. In the present study, we firstly carried out numerical experiments to address the difference between the outcomes of local and global modeling methods, demonstrating that local modeling confined to the aneurysm region results in inaccurate predictions of hemodynamic parameters compared with global modeling of the ACoA aneurysm as part of the cerebral artery network. Motivated by this finding, we built global hemodynamic models for 40 ACoA aneurysms (including 20 ruptured and 20 unruptured ones) based on medical image data. Statistical analysis of the computed hemodynamic data revealed that maximum wall shear stress (WSS), minimum WSS divergence, and maximum WSS gradient differed significantly between the ruptured and unruptured ACoA aneurysms. Optimal threshold values of high/low WSS metrics were determined through a series of statistical tests. In the meantime, some morphological parameters of aneurysms, such as large nonsphericity index, aspect ratio, and bottleneck factor, were found to be associated closely with aneurysm rupture. Furthermore, multivariate logistic regression analyses were performed to derive models combining hemodynamic and morphological parameters for discriminating the rupture status of aneurysms. The capability of the models in rupture status discrimination was high, with the area under the receiver operating characteristic curve reaching up to 0.9. The findings of the study suggest that global modeling of the cerebral artery network is essential for reliable quantification of hemodynamics in ACoA aneurysms, disturbed WSS and irregular aneurysm morphology are associated closely with aneurysm rupture, and multivariate models integrating hemodynamic and morphological parameters have high potential for assessing the rupture risk of ACoA aneurysms.

Predicting red blood cell traffcking and capillary hemodynamics in angiogenic and tumor microcirculation in silico

Objective: Angiogenic and tumor microvasculatures are known to have abnormal topology due to the presence of frequent vessel junctions, irregular and deflated blood vessels, multi-furcations, and tessellated vessel organization. Although recent advances in imaging techniques in vivo have enabled mapping such vasculatures at high spatial resolution, simultaneous measurements of hemodynamic parameters, such as the wall shear stress (WSS) with full 3D details, remain a challenge. Theoretical network flow models, often used for hemodynamic predictions in such experimentally acquired images, cannot provide the full 3D hemodynamic details either, as these models treat each blood vessel as 1D segment and do not explicitly model red blood cells (RBCs). To overcome this limitation, we have developed a high-fidelity, 3D Computational Fluid Dynamics modeling to predict the flow of a large number of deformable RBCs through physiologically realistic tumor/angiogenic microvascular networks in silico. Methods: We use in vivo images to create such vascular networks in silico and then predict RBC traffcking and capillary hemodynamics. Deformation of every flowing RBC is considered with high accuracy, and 3D geometry of each vessel is accurately modeled. Flow is driven by specifying physiological pressure boundary conditions. Model predictions have been validated against in vivo data. This in-house predictive tool is versatile, can be applied to any microvascular network image obtained in vivo in any organ, and can predict trajectories of diverse cell types including leukocytes, platelets and circulating tumor cells, drug and molecular transport in capillary blood, and cell-vessel adhesion. Results: We provide quantitative differences between healthy microvascular networks and tumor/angiogenic networks in terms of RBC distribution, perfusion, and wall shear stress. Our model shows increased heterogeneity in RBC and flow distribution in both tumor and angiogenic vasculatures than the healthy one. Also, we predict reduced flow and hematocrit in several vessels in both tumor and angiogenic vasculatures. Interestingly, several vessels in the angiogenic vasculature are predicted to have higher flow than the healthy one, while most vessels in the tumor vasculature show flow reduction. This in silico prediction is consistent with a recent in vivo study which showed higher flow in peri-tumor region and reduced flow in tumor. We further predict a significant heterogeneity in WSS and WSS gradient, blood velocity profiles, and near-wall RBC-depleted region. Conclusion: In conclusion, we have developed a versatile, in silico model that allows high-fidelity prediction of capillary hemodynamics in tumor microcirculation and provide information on hemodynamic variables that are not readily measurable in vivo but have physiological significance in tumor progression and treatment. NIH (R01EY033003) and NSF (CBET1804591). This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Experimental study on hemodynamics of an end-to-side anastomosis

A three-dimensional and three-component velocity measurement on the flow field in a 45° end-to-side anastomosis model is conducted to investigate the hemodynamics, which is an important factor to the intimal hyperplasia formation and graft failure after surgery. Thanks to the advanced volumetric measurement technology of tomographic particle image velocimetry, the recirculation zone, low-speed region, and the spiral flow structures can be visualized. As a result, the flow field of three cases with the local maximum velocity of 0.15, 0.8, and 1.4 m/s are visible and the inlet velocity profile tends to be skewed as the flow rate increases. The mean vorticity contours indicate that the positive vortex center rotates 6.47°, 50.23°, and 90.4° and the negative vortex center rotates 20.44°, 15.73°, and 68.47°, respectively, in three cases. The instantaneous vortex structures identified by the λci criterion demonstrate two large-scale vortex structures in the distal section. The two vortices have the tendency to intertwine while one of them decays earlier. The wall shear stress (WSS) distributions on the entire model with the local maximum of 0.8, 5.8, and 13.8 Pa in three cases have been quantitatively achieved. The abnormal WSS and WSS gradient can help localize risk areas and understand the intimal hyperplasia formation. A detailed illustration of hemodynamics inside the 45° end-to-side anastomosis model has been provided, which demonstrates more comprehensive large-scale flow structures and abnormal WSS regions. Combined with the information of flow structures and WSS distribution, the understanding of the hemodynamics in the anastomosis can be strengthened.

A parametric study of the effect of 3D plaque shape on local hemodynamics and implications for plaque instability

The vast majority of heart attacks occur when vulnerable plaques rupture, releasing their lipid content into the blood stream leading to thrombus formation and blockage of a coronary artery. Detection of these unstable plaques before they rupture remains a challenge. Hemodynamic features including wall shear stress (WSS) and wall shear stress gradient (WSSG) near the vulnerable plaque and local inflammation are known to affect plaque instability. In this work, a computational workflow has been developed to enable a comprehensive parametric study detailing the effects of 3D plaque shape on local hemodynamics and their implications for plaque instability. Parameterized geometric 3D plaque models are created within a patient-specific coronary artery tree using a NURBS (non-uniform rational B-splines)-based vascular modeling pipeline. Realistic blood flow features are simulated by using a Navier–Stokes solver within an isogeometric finite-element analysis framework. Near wall hemodynamic quantities such as WSS and WSSG are quantified, and vascular distribution of an inflammatory marker (VCAM-1) is estimated. Results show that proximally skewed eccentric plaques have the most vulnerable combination of high WSS and high positive spatial WSSG, and the presence of multiple lesions increases risk of rupture. The computational tool developed in this work, in conjunction with clinical data, -could help identify surrogate markers of plaque instability, potentially leading to a noninvasive clinical procedure for the detection of vulnerable plaques before rupture.

Hemodynamic factors of spontaneous vertebral artery dissecting aneurysms assessed with numerical and deep learning algorithms: Role of blood pressure and asymmetry

Background and ObjectivesThe pathophysiology of spontaneous vertebral artery dissecting aneurysms (SVADA) is poorly understood. Our goal is to investigate the hemodynamic factors contributing to their formation using computational fluid dynamics (CFD) and deep learning algorithms. MethodsWe have developed software that can use patient imagery as input to recreate the vertebrobasilar arterial system, both with and without SVADA, which we used in a series of three patients. To obtain the kinematic blood flow data before and after the aneurysm forms, we utilized numerical methods to solve the complex Navier-Stokes partial differential equations. This was accomplished through the application of a finite volume solver (OpenFoam/Helyx OS). Additionally, we trained a neural ordinary differential equation (NODE) to learn and replicate the dynamical streamlines obtained from the computational fluid dynamics (CFD) simulations. ResultsIn all three cases, we observed that the equilibrium of blood pressure distributions across the VAs, at a specific vertical level, accurately predicted the future SVADA location. In the two cases where there was a dominant VA, the dissection occurred on the dominant artery where blood pressure was lower compared to the contralateral side. The SVADA sac was characterized by reduced wall shear stress (WSS) and decreased velocity magnitude related to increased turbulence. The presence of a high WSS gradient at the boundary of the SVADA may explain its extension. Streamlines generated by CFD were learned with a neural ordinary differential equation (NODE) capable of capturing the system’s dynamics to output meaningful predictions of the flow vector field upon aneurysm formation. ConclusionIn our series, asymmetry in the vertebrobasilar blood pressure distributions at and proximal to the site of the future SVADA accurately predicted its location in all patients. Deep learning algorithms can be trained to model blood flow patterns within biological systems, offering an alternative to the computationally intensive CFD. This technology has the potential to find practical applications in clinical settings.

Standardized viscosity as a source of error in computational fluid dynamic simulations of cerebral aneurysms.

Computational fluid dynamics (CFD) simulations are a powerful tool for studying cerebral aneurysms, capable of evaluating hemodynamics in a way that is infeasible with imaging alone. However, the difficulty of incorporating patient-specific information and inherent obstacles of in vivo validation have limited the clinical usefulness of CFD of cerebral aneurysms. In this work we investigate the effect of using standardized blood viscosity values in CFD simulations of cerebral aneurysms when compared to simulations of the same aneurysms using patient-specific viscosity values derived from hematocrit measurements. The objective of this work is to determine the level of error, on average, that is, caused by using standardized values of viscosity in CFD simulations of cerebral aneurysms. By quantifying this error, we demonstrate the need for incorporating patient-specific viscosity in future CFD investigations of cerebral aneurysms. CFD simulations of forty-one cerebral aneurysms were conducted using patient-specific boundary conditions. For each aneurysm two simulations were conducted, one utilizing patient-specific blood viscosity derived from hematocrit measurements and another using a standardized value for blood viscosity. Hemodynamic parameters such as wall shear stress (WSS), wall shear stress gradient (WSSG), and the oscillatory shear index (OSI) were calculated for each of the simulations for each aneurysm. Paired t-tests for differences in the time-averaged maps of these hemodynamic parameters between standardized and patient-specific viscosity simulations were conducted for each aneurysm. Bland-Altman analysis was used to examine the cohort-wide changes in the hemodynamic parameters. Subjects were broken into two groups, those with higher than standard viscosity and those with lower than standard viscosity. An unpaired t-test was used to compare the percent change in WSS, WSSG, and OSI between patient-specific and standardized viscosity simulations for the two cohorts. The percent changes in hemodynamic parameters were correlated against the direction and magnitude of percent change in viscosity, aneurysm size, and aneurysm location. For all t-tests, a Bonferroni-corrected significance level of 0.0167 was used. 63.2%, 41.5%, and 48.7% of aneurysms showed statistically significant differences between patient-specific and standardized viscosity simulations for WSS, WSSG, and OSI respectively. No statistically significant difference was found in the percent changes in WSS, WSSG, and OSI between the group with higher than standard viscosity and those with lower than standard viscosity, indicating an increase in viscosity can cause either an increase or decrease in each of the hemodynamic parameters. On a study-wide level no significant bias was found in either direction for WSS, WSSG, or OSI between the simulation groups due to the bidirectional effect of changing viscosity. No correlation was found between percent change of viscosity and percent change of WSS, WSSG, or OSI, meaning an after-the-fact correction for patient-specific viscosity is not feasible. Standardizing viscosity values in CFD of cerebral aneurysms has a large and unpredictable impact on the calculated WSS, WSSG, and OSI when compared to CFD simulations of the same aneurysms using a patient-specific viscosity. We recommend implementing hematocrit-based patient-specific blood viscosity values for all CFD simulations of cerebral aneurysms.

A cell-and-plasma numerical model reveals hemodynamic stress and flow adaptation in zebrafish microvessels after morphological alteration.

The development of a functional cardiovascular system ensures a sustainable oxygen, nutrient and hormone delivery system for successful embryonic development and homeostasis in adulthood. While early vessels are formed by biochemical signaling and genetic programming, the onset of blood flow provides mechanical cues that participate in vascular remodeling of the embryonic vascular system. The zebrafish is a prolific animal model for studying the quantitative relationship between blood flow and vascular morphogenesis due to a combination of favorable factors including blood flow visualization in optically transparent larvae. In this study, we have developed a cell-and-plasma blood transport model using computational fluid dynamics (CFD) to understand how red blood cell (RBC) partitioning affect lumen wall shear stress (WSS) and blood pressure in zebrafish trunk blood vascular networks with altered rheology and morphology. By performing live imaging of embryos with reduced hematocrit, we discovered that cardiac output and caudal artery flow rates were maintained. These adaptation trends were recapitulated in our CFD models, which showed reduction in network WSS via viscosity reduction in the caudal artery/vein and via pressure gradient weakening in the intersegmental vessels (ISVs). Embryos with experimentally reduced lumen diameter showed reduced cardiac output and caudal artery flow rate. Factoring in this trend into our CFD models, simulations highlighted that lumen diameter reduction increased vessel WSS but this increase was mitigated by flow reduction due to the adaptive network pressure gradient weakening. Additionally, hypothetical network CFD models with different vessel lumen diameter distribution characteristics indicated the significance of axial variation in lumen diameter and cross-sectional shape for establishing physiological WSS gradients along ISVs. In summary, our work demonstrates how both experiment-driven and hypothetical CFD modeling can be employed for the study of blood flow physiology during vascular remodeling.

Multi-effect analysis of nanofluid flow in stenosed arteries with variable pressure gradient: analytical study

This study advances the understanding of nanofluid behaviour within stenosed arteries, highlighting the importance of considering multifaceted effects in the modelling process. It investigates the combined impact of pressure gradient variation, heat transfer, chemical reactions, and magnetic field effects on nano-blood flow in stenosed arteries. Unlike previous studies that made the assumption that the pulsatile pressure gradient remains constant during channel narrowing, this novel investigation introduces a variable pressure gradient. This, in turn, significantly impacts several associated parameters. The mathematical model describing nano-blood flow in a horizontally stenosed artery is solved using perturbation techniques. Analytical solutions for key variables, including velocity, temperature, concentration, wall shear stress, flow rate, and pressure gradient, are visually presented for various physical parameter values.

Numerical investigation of unruptured middle cerebral artery bifurcation aneurysms: influence of aspect ratio

An aneurysm is a disease condition, which is due to the pathological weakening of an arterial wall. These aneurysms are often found in various branch points and bifurcations of an artery in the cerebral circulation. Most aneurysms come to medical attention, either due to brain hemorrhages caused by rupture or found unruptured. To consider surgically invasive treatment modalities, clinicians need scientific methods such as, hemodynamic analysis to assess rupture risk. The arterial wall loses its structural integrity when wall shear stress (WSS) and other hemodynamic parameters exceed a certain threshold. In the present study, numerical simulations are carried out for unruptured middle cerebral artery (MCA) aneurysms. Three distinct representative sizes are chosen from a larger patient pool of 26 MCA aneurysms. Logically, these aneurysms represent three growth stages of any patient with similar anatomical structure. Simulations are performed to compare the three growth phases (with different aspect ratios) of an aneurysm and correlate their hemodynamic parameters. Simulations with patient specific boundary conditions reveal that, aneurysms with a higher aspect ratio (AR) correspond to an attendant decrease in both time-averaged wall shear stress (TAWSS) and spatial wall shear stress gradients (WSSG). Smaller MCAs were observed to have higher positive wall shear stress divergence (WSSD), exemplifying the tensile nature of arterial wall stretching. Present study identifies positive wall shear stress divergence (PWSSD) to be a potential biomarker for evaluating the growth of an aneurysm.

Hemodynamic analysis of coil filled patient-specific middle cerebral artery aneurysm using porous medium approach

Cerebral aneurysms are bulges of an artery, which could be life-threatening when ruptured. Depending on their size, shape, and location, they need to be managed either through clipping or an endovascular coiling intervention. When coiled, reduced hemodynamic activity enables the coil to get thrombosed and achieve flow stasis. However, some coils delivered into the aneurysm tend to prolapse into the parent vessel and cause stroke due to obstruction and embolization. The recurrence of an aneurysm after endovascular coiling is of concern in the treatment of wide necked aneurysms. The initial packing density or improper coiling of the aneurysm and its relation to recurrence remains uncertain. This study investigates the influence of reduction in coil fill volume and packing density on the aneurysm recurrence using hemodynamic parameters by analyzing its flow features. Finite element method based commercial computational fluid dynamics solver is employed for performing patient-specific simulations for the coil filled aneurysm. The present approach uses porous medium based formulation. The numerical simulations show that any reduction below the optimal coil fill volume and packing density inside the aneurysm increases the velocity magnitude, wall shear stress, time-averaged wall shear stress, and spatial gradient of wall shear stress and reduces the relative residence time. The hemodynamic parameters and flow features suggest that a reduction in the coil packing density inside the aneurysm increases the chances of aneurysm recurrence. Hence, an assessment on how to achieve optimal coil fill volume and packing density is critical in reducing the risk of aneurysm recurrence.

Association of Arterial Tortuosity with Hemodynamic Parameters—A Computational Fluid Dynamics Study

Tortuosity of intracranial arteries has been proven to be associated with the risk of intracranial aneurysm development. We decided to analyze which aspects of tortuosity are correlated with hemodynamics parameters promoting intracranial aneurysm development. We constructed 73 idealized models of tortuous artery (length: 25 mm, diameter: 2.5 mm) with single bifurcation. For each model, on the course of segment before bifurcation, we placed 1-3 angles with measures 15, 30, 45, 60, or 75 degrees and arc lengths 2, 5, 7, 10, or 15 mm. We performed computational fluid dynamics analysis. Blood was modeled as Newtonian fluid. We have set velocity wave of 2 cardiac cycles. After performing simulation we calculated following hemodynamic parameters at the bifurcation: time average wall shear stress (TAWSS), time average wall shear stress gradient (TAWSSG), oscillatory shear index (OSI), and relative residence time (RRT). We found a significant positive correlation with number of angles and TAWSS (R= 0.329; P < 0.01), TAWSSG (R= 0.317; P < 0.01), and negative with RRT (R=-0.335; P < 0.0.01). Similar results were obtained in terms of arcs lengths. On the other hand, mean angle measure was negatively correlated to TAWSS (R=-0.333; P < 0.01), TAWSSG (R=-0.473 P < 0.01), OSI (R=-0.463; P < 0.01), and positively to RRT (R= 0.332; P < 0.01). On the basis of the obtained results, we developed new tortuosity descriptor, which considered angle measures normalized to its arc length and distance from bifurcation. For such descriptor we found strong negative correlation with TAWSS (R=-0.701; P < 0.01), TAWSSG (R= 0.778; P < 0.01), OSI (R=-0.776; P < 0.01), and positive with RRT (R= 0.747; P < 0.01). Hemodynamic parameters promoting aneurysm development are correlated with larger number of smaller angles located on larger arcs.

Pathophysiology of intracranial aneurysms in monozygotic twins: A rare case study from hemodynamic perspectives

Hemodynamic mechanisms of the formation and growth of intracranial aneurysms (IA) in monozygotic twins (MTs) are still under–reported. To partially fill such knowledge gap, this study employed an experimentally validated numerical model to compare hemodynamics in 3 anatomical and 5 ablation study neurovascular models from a rare pair of MTs in terms of 7 critical hemodynamic parameters. Numerical results showed significant differences in hemodynamics between the MTs, although they share the same genes, indicating that genetic mutation and environmental factors might affect neurovascular morphologies and cause hemodynamic changes. After virtual removals of IAs in the ablation study, the locations where the aneurysmal sac/bleb generated in bifurcated anterior cerebral arteries (ACAs) register a locally high instantaneous wall shear stress (IWSS) of 52.9 and 70.1 Pa at the systolic peak in twin A and twin B, respectively. Same scenario can be observed in the distribution of instantaneous wall shear stress gradient (IWSSG), with 571.1 Pa/mm for twin A and 301.3 Pa/mm for twin B due to aggressive blood impingements, leading to IA generation. The fenestrated complex approaching ACA bifurcations in twin A may assist IA growth and rupture, via. Causing abnormal IWSS of 116.3 Pa, IWSSG of 832.5 Pa/mm, and oscillatory shear index (OSI) of 0.49. The bleb in twin B has high risks of progression and possible rupture as the IA suffers relatively low IWSS and high OSI. Additionally, IA generation can change blood flow rates in each connected artery, then affecting blood supplies to associated tissues and organs.

Wall enhancement predictive of abnormal hemodynamics and ischemia in vertebrobasilar non-saccular aneurysms: a pilot study.

To analyze how wall enhancement affects hemodynamics and cerebral ischemic risk factors in vertebrobasilar non-saccular intracranial aneurysms (VBNIAs). Ten consecutive non-saccular aneurysms were collected, including three transitional vertebrobasilar dolichoectasia (TVBD). A wall enhancement model was quantitatively constructed to analyze how wall enhancement interacts with hemodynamics and cerebral ischemic factors. Enhanced area revealed low wall shear stress (WSS) and wall shear stress gradient (WSSG), with high oscillatory shear index (OSI), relative residence time (RRT), and gradient oscillatory number (GON) while the vortex and slow flow region in fusiform aneurysms are similar to TVBD fusiform aneurysms. With low OSI, high RRT and similar GON in the dilated segment, the enhanced area still manifests low WSS and WSSG in the slow flow area with no vortex. In fusiform aneurysms, wall enhancement was negatively correlated with WSS (except for case 71, all p values < 0.05, r = -0.52 ~ -0.95), while wall enhancement was positively correlated with OSI (except for case 5, all p values < 0.05, r = 0.50 ~ 0.83). For the 10 fusiform aneurysms, wall enhancement is significantly positively correlated with OSI (p = 0.0002, r = 0.75) and slightly negatively correlated with WSS (p = 0.196, r = -0.30) throughout the dataset. Aneurysm length, width, low wall shear stress area (LSA), high OSI, low flow volume (LFV), RRT, and high aneurysm-to-pituitary stalk contrast ratio (CRstalk) area plus proportion may be predictive of cerebral ischemia. A wall enhancement quantitative model was established for vertebrobasilar non-saccular aneurysms. Low WSS was negatively correlated with wall enhancement, while high OSI was positively correlated with wall enhancement. Fusiform aneurysm hemodynamics in TVBD are similar to simple fusiform aneurysms. Cerebral ischemia risk appears to be correlated with large size, high OSI, LSA, and RRT, LFV, and wall enhancement.

A spatiotemporal analysis of the left coronary artery biomechanics using fluid-structure interaction models.

Biomechanics plays a critical role in coronary artery disease development. FSI simulation is commonly used to understand the hemodynamics and mechanical environment associated with atherosclerosis pathology. To provide a comprehensive characterization of patient-specific coronary biomechanics, an analysis of FSI simulation in the spatial and temporal domains was performed. In the current study, a three-dimensional FSI model of the LAD coronary artery was built based on a patient-specific geometry using COMSOL Multiphysics. The effect of myocardial bridging was simulated. Wall shear stress and its derivatives including time-averaged wall shear stress, wall shear stress gradient, and OSI were calculated across the cardiac cycle in multiple locations. Arterial wall strain (radial, circumferential, and longitudinal) and von Mises stress were calculated. To assess perfusion, vFFR was calculated. The results demonstrated the FSI model could identify regional and transient differences in biomechanical parameters within the coronary artery. The addition of myocardial bridging caused a notable change in von Mises stress and an increase in arterial strain during systole. The analysis performed in this manner takes greater advantage of the information provided in the space and time domains and can potentially assist clinical evaluation.

Comprehensive morphomechanical analysis of brain aneurysms.

Brain aneurysms comprise different compartments that undergo unique biological processes. A detailed multimodal analysis incorporating 3D aneurysm wall enhancement (AWE), computational fluid dynamics (CFD), and finite element analysis (FEA)data can provide insights into the aneurysm wall biology. Unruptured aneurysms were prospectively imaged with 7T high-resolution MRI (HR-MRI). 3D AWE color maps of the entire aneurysm wall were generated and co-registered with contour plots of morphomechanical parameters derived from CFD and FEA. A multimodal analysis of the entire aneurysm was performed using 3D circumferential AWE (3D-CAWE), wall tension (WT), time-averaged wall shear stress (TAWSS), wall shear stress gradient (WSSG), and oscillatory shear index (OSI). A detailed compartmental analysis of each aneurysm's dome, bleb, and neck was also performed. Twenty-six aneurysms were analyzed. 3D-CAWE + aneurysms had higher WT (p = 0.03) and higher TAWSS (p = 0.045) than 3D-CAWE- aneurysms. WT, TAWSS, and WSSG were lower in areas of focal AWE in the aneurysm dome compared to the neck (p = 0.009, p = 0.049, and p = 0.040, respectively), whereas OSI was higher in areas of focal AWE compared to the neck (p = 0.020). When compared to areas of no AWE of the aneurysm sac (AWE = 0.92 vs. 0.49, p = 0.001), blebs exhibited lower WT (1.6 vs. 2.45, p = 0.010), lower TAWSS (2.6 vs. 6.34), lower OSI (0.0007 vs. 0.0010), and lower WSSG (2900 vs. 5306). Fusiform aneurysms had a higher 3D-CAWE and WT than saccular aneurysms (p = 0.046 and p = 0.003, respectively). Areas of focal high AWE in the sac and blebs are associated with low wall tension, low wall shear stress, and low flow conditions (TAWSS and WSSG). Conversely, the neck had average AWE, high wall tension, high wall shear stress, and high flow conditions. The aneurysm dome and the aneurysm neck have different morphomechanical environments, with increased mechanical load at the neck.

Analysis of Hemodynamic Effect on Different Stent Strut Configuration in Femoral Popliteal Artery during the Physical and Physiological Conditions

Peripheral arterial disease (PAD) is a narrowing of the peripheral arteries that might result in blockage if not immediately treated. Normally, an invasive technique called stenting is used at the stenosed arterial region to restore normal blood flow. Thus, it promotes the formation of thrombosis on the stented artery due to the presenting flow recirculation. However, the rate of thrombosis growth was reported to be different for both genders. This is due to an increase in body surface area, body mass index, and weight of the body. Thus, this study aims to investigate the effect of the physiological and physical conditions of men and women with different hemodynamic parameters on the strut configuration in FPA. Five different stent strut configurations were modelled and inserted into the FPA. The computational fluid dynamic (CFD) method was implemented to solve the continuity and N-S equations. The hemodynamic performance of the stent was analyzed based on hemodynamic parameters consisting of time- averaged wall shear stress (TAWSS), time-averaged wall shear stress gradient (TAWSSG), oscillatory shear index (OSI), and relative residence time (RRT). According to the observations, the distal region of the stented FPA had more dominant flow re-circulation than the proximal region. The high void area contributed to less growth of the thrombosis.