FSI modeling on the effect of artery-aneurysm thickness and coil embolization in patient cases

Nonlinear biomechanical behaviour of extracranial carotid artery aneurysms in the framework of Windkessel effect via FSI technique

Extracranial carotid artery aneurysms (ECCA) lead to rupture and neurologic symptoms from embolisation, with potentially fatal outcomes. Investigating the biomechanical behaviour of EECA with blood flow dynamics is crucial for identifying regions more susceptible to rupture. A coupled three-dimensional (3D) Windkessel-framework and hyperelastic fluid-structure interaction (FSI) analysis of ECCAs with patient-specific geometries, was developed in this paper with a particular focus on hemodynamic parameters and the arterial wall's biomechanical response. The blood flow has been modelled as non-Newtonian, pulsatile, and turbulent. The biomechanical characteristics of the aneurysm and artery are characterised employing a 5-parameter Mooney-Rivlin hyperelasticity model. The Windkessel effect is also considered to efficiently simulate pressure profile of the outlets and to capture the dynamic changes over the cardiac cycle. The study found the aneurysm carotid artery exhibited the high levels of pressure, wall shear stress (WSS), oscillatory shear index (OSI), and relative residence time (RRT) compared to the healthy one. The deformation of the arterial wall and the corresponding von Mises (VM) stress were found significantly increased in aneurysm cases, in comparison to that of no aneurysm cases, which strongly correlated with the hemodynamic characteristics of the blood flow and the geometric features of the aneurysms. This escalation would intensify the risk of aneurysm wall rupture. These findings have critical implications for enhancing treatment strategies for patients with extracranial aneurysms.

Patient-specific arterial wall generation for intracranial aneurysms with a variable and a near realistic vessel wall thickness for FSI studies

BACKGROUND AND OBJECTIVE

Imaging methodologies such as, computed tomography (CT) aid in three-dimensional (3D) reconstruction of patient-specific aneurysms. The radiological data is useful in understanding their location, shape, size, and disease progression. However, there are serious impediments in discerning the blood vessel wall thickness due to limitations in the current imaging modalities. This further restricts the ability to perform high-fidelity fluid structure interaction (FSI) studies for an accurate assessment of rupture risk. FSI studies would require the arterial wall mesh to be generated to determine realistic maximum allowable wall stresses by performing coupled calculations for the hemodynamic forces with the arterial walls.

METHODS

In the present study, a novel methodology is developed to geometrically model variable vessel wall thickness for the lumen isosurface extracted from CT scan slices of patient-specific aneurysms based on clinical and histopathological inputs. FSI simulations are carried out with the reconstructed models to assess the importance of near realistic wall thickness model on rupture risk predictions.

RESULTS

During surgery, clinicians often observe translucent vessel walls, indicating the presence of thin regions. The need to generate variable vessel wall thickness model, that embodies the wall thickness gradation, is closer to such clinical observations. Hence, corresponding FSI simulations performed can improve clinical outcomes. Considerable differences in the magnitude of instantaneous wall shear stresses and von Mises stresses in the walls of the aneurysm was observed between a uniform wall thickness and a variable wall thickness model.

CONCLUSION

In the present study, a variable vessel wall thickness generation algorithm is implemented. It was shown that, a realistic wall thickness modeling is necessary for an accurate prediction of the shear stresses on the wall as well as von Mises stresses in the wall. FSI simulations are performed to demonstrate the utility of variable wall thickness modeling.

A new self-adjustable glaucoma valve

Introduction: Glaucoma, the leading cause of irreversible blindness globally, affects more than 70 million people across the world. When initial treatments prove ineffective, especially for cases with high intraocular pressure (IOP), the preferred approach involves employing glaucoma drainage devices (GDDs). Methods: This study introduces a novel self-adjustable glaucoma drainage device (SAGDD) designed to maintain IOP within the desired biological range (10 mmHg < IOP <18 mmHg) by dynamically modulating its fluidic resistance. Inspired by the starling resistor, we designed a circular valve with a thin, flexible membrane placed over the valve's inlet and outlet. To achieve the ideal design for the SAGDD and optimize its parameters, we utilized fluid-solid interaction (FSI) numerical models and conducted parametric studies, wherein simulations demonstrated the validity of the concept. Subsequently, to confirm and validate the numerical results, we fabricated a SAGDD at a 3:1 scale and subjected it to in vitro testing. Results: Our findings demonstrate that, on a 3:1 scale, a circular SAGDD with a diameter of 8.1 mm and a stainless-steel membrane with a thickness of 10 µm effectively maintained IOP within the target range when the membrane exposed to external pressures of 7.5 or 10 mmHg. Discussion: In summary, our study establishes a strong foundation for further exploration of the potential efficacy of SAGDD as a promising treatment for glaucoma. The cost-effectiveness and simplicity of its design, devoid of costly instrumentation, hold considerable promise in addressing the challenges associated with glaucoma.

Influence of circle of Willis modeling on hemodynamic parameters in anterior communicating artery aneurysms and recommendations for model selection

Computational fluid dynamics (CFD) has been utilized to calculate hemodynamic parameters in anterior communicating artery aneurysm (AComA), which is located at a junction between left and right A1 and A2 segments. However, complete or half circle of Willis (CoW) models are used indiscriminately. This study aims to suggest recommendations for determining suitable CoW model. Five patient-specific CoW models with AComA were used, and each model was divided into complete, left-half, and right-half models. After validating the CFD using a flow experiment, the hemodynamic parameters and flow patterns in five AComAs were compared. In four out of five cases, inflow from one A1 side had a dominant influence on the AComA, while both left and right A1 sides affected the AComA in the remaining case. Also, the average difference in time-averaged wall shear stress between the complete and half models for four cases was 4.6%, but it was 62% in the other case. The differences in the vascular resistances of left and right A1 and A2 segments greatly influenced the flow patterns in the AComA. These results may help to enhance clinicians' understanding of blood flow in the brain, leading to improvements in diagnosis and treatment of cerebral aneurysms.

Modeling and evaluation of biomechanics and hemodynamic based on patient-specific small intracranial aneurysm using fluid-structure interaction

BACKGROUND AND OBJECTIVE

Rupture of small intracranial aneurysm (IA) often leads to the development of highly fatal clinical syndromes such as subarachnoid hemorrhage. Due to the patient specificity of small IA, there are many difficulties in evaluating the rupture risk of small IA such as multiple influencing factors, high clinical experience requirements and poor reusability.

METHODS

In this study, clinical methods such as transcranial doppler (TCD) and magnetic resonance imaging (MRI) are used to obtain patient-specific parameters, and the fluid-structure interaction method (FSI) is used to model and evaluate the biomechanics and hemodynamics of patient-specific small IA.

RESULTS

The results show that a spiral vortex stably exists in the patient-specific small IA. Due to the small size of the patient-specific small IA, the blood flow velocity still maintains a high value with maximum reaching 3 m/s. The inertial impact of blood flow and vortex convection have certain influence on hemodynamic and biomechanics parameters. They cause three high value areas of WSSM on the patient-specific small IA with maximum of 180 Pa, 130 Pa and 110 Pa, respectively. They also cause two types of WSS concentration points, positive normal stress peak value areas and negative normal stress peak value areas to appear.

CONCLUSION

This paper found that the factors affecting hemodynamic parameters and biomechanical parameters are different. Unlike hemodynamic parameters, biomechanical parameters are also affected by blood pressure in addition to blood flow velocity. This study reveals the relationship between the flow field distribution and changes of patient-specific small IA, biomechanics and hemodynamics.

A novel gas embolotherapy using microbubbles electrocoalescence for cancer treatment

BACKGROUND AND OBJECTIVE

Embolotherapy has been increasingly used to disrupt tumor growth. Despite its success in the occlusion of microvessels, it has drawbacks such as limited access to the target location, limited control of the blocker size, and inattention to the tumor characteristics, especially high interstitial fluid pressure. The present work introduces a novel numerical method of gas embolotherapy for cancer treatment through tumor vessel occlusion.

METHODS

The gas microbubbles are generated from Levovist bolus injection into the tumor microvessel. The microbubble movement in the blood flow is innovatively controlled by an electric field applied to the tumor-feeding vessel. The interaction between the Levovist microbubbles and the electric field is resolved by developing a fully coupled model using the phase-field model, Carreau model for non-Newtonian blood, Navier-Stokes equations and Maxwell stress tensor. Additionally, the critical effect of high interstitial fluid pressure as a characteristic of solid tumors is included.

RESULTS

The findings of this study indicate that the rates of microbubble deformation and displacement increase with the applied potential intensity to the microvessel wall. Accordingly, the required time for a microbubble to join the upper microvessel wall reduces from 1.97ms to 22 μs with an increase of the electric potential from 3.5V to 12.5V. Additionally, an electric potential of 12.5V causes the microbubbles coalescence and formation of a gas column against the bloodstream.

CONCLUSIONS

Clinically, our novel embolization procedure can be considered a non-invasive targeted therapy, and under a controlled electric field, the blocker size can be precisely controlled. Also, the proposed method has the potential to be used as a gradual treatment in advanced cancers as tumors develop resistance and relapse.

A novel parameter for the prediction of rupture risk of cerebral aneurysms based on morphology

Neurosurgeons often encounter dilemmas in the clinical management of cerebral aneurysms owing to an uncertainty of their rupture status and rupture risk. This study evaluates the influence of natural frequency of an aneurysm, as a novel morphological parameter to understand and analyze rupture status and risk prediction. In this work, we employ the natural frequency of 20 idealized and 50 patient specific aneurysms. The natural frequency of patient specific aneurysms is then compared against their rupture status. A strong correlation was observed between various morphological indicators and natural frequency for ideal and patient specific geometries. A statistical analysis with both Mann Whitney U test and T-test for rupture status against natural frequency has given a p-value less than 0.01 indicating a strong correlation between them. The correlation of morphological parameters with natural frequency from Pearson correlation coefficient and T-test suggests a holistic reflection of their effects on the natural frequency of an aneurysm. Thus, natural frequency could be a good indicator to discern the rupture potential of an aneurysm. The correlation between rupture status and natural frequency makes it a novel parameter that can differentiate between ruptured and unruptured patient specific aneurysms.

Numerical Study of a Thrombus Migration Risk in Aneurysm After Coil Embolization in Patient Cases: FSI Modelling

PURPOSE

There are still many challenges for modelling a thrombus migration process in aneurysms. The main novelty of the present research lies in the modelling of aneurysm clot migration process in a realistic cerebral aneurysm, and the analysis of forces suffered by clots inside an aneurysm, through transient FSI simulations.

METHODS

The blood flow has been modelled using a Womersley velocity profile, and following the Carreau viscosity model. Hyperelastic Ogden model has been used for clot and isotropic linear elastic model for the artery walls. The FSI coupled model was implemented in ANSYS® software. The hemodynamic forces suffered by the clot have been quantified using eight different clot sizes and positions inside a real aneurysm.

RESULTS

The obtained results have shown that it is almost impossible for clots adjacent to aneurysm walls, to leave the aneurysm. Nevertheless, in clots positioned in the centre of the aneurysm, there is a real risk of clot migration. The risk of migration of a typical post-coiling intervention clot in an aneurysm, in contact with the wall and occupying a significant percentage of its volume is very low in the case studied, even in the presence of abnormally intense events, associated with sneezes or impacts.

CONCLUSIONS

The proposed methodology allows evaluating the clot migration risk, vital for evaluating the progress after endovascular interventions, it is a step forward in the personalized medicine, patient follow-up, and helping the medical team deciding the optimal treatment.

A personalized 0D-1D model of cardiovascular system for the hemodynamic simulation of enhanced external counterpulsation

BACKGROUND AND OBJECTIVE

Enhanced external counterpulsation (EECP) is a non-invasive treatment modality capable of treating a variety of ischemic diseases. Currently, no effective methods of predicting the patient-specific hemodynamic effects of EECP are available. In this study, a personalized 0D-1D model of the cardiovascular system was developed for hemodynamic simulation to simulate the changes in blood flow in the EECP state and develop the best treatment protocol for each individual.

METHODS

A 0D-1D closed-loop model of the cardiovascular system was developed for hemodynamic simulation, consisting of a 1D wave propagation model for arteries, a 0D model for veins and capillaries, and a one-fiber model for the heart. Additionally, a simulation model coupling EECP with a 1D model was established. Physiological data, including the blood flow in different arteries, were clinically collected from 22 volunteers at rest and in the EECP state. Sensitivity analysis and a simulated annealing algorithm were used to build personalized 0D-1D models using the clinical data in the rest state as optimization objectives. Then, the clinical data on EECP were used to verify the applicability and accuracy of the personalized models.

RESULTS

The simulation results and clinical data were found to be in agreement for all 22 subjects, with waveform similarity coefficients (r) exceeding 90% for most arteries at rest and 80% for most arteries during EECP.

CONCLUSIONS

The 0D-1D closed-loop model and the optimized method can facilitate personalized modeling of the cardiovascular system using the data in the rest state and effectively predict the hemodynamic changes in the EECP state, which is significant for the numerical simulation of personalized hemodynamics. The model can also potentially be used to make decisions regarding patient-specific treatment.

Hemodynamics and bio-mechanics of morphologically distinct saccular intracranial aneurysms at bifurcations: Idealised vs Patient-specific geometries

Understanding the factors that influence the rupture of aneurysms is of primary concern to the clinicians, who are grappled with patient management. It is important to know how the relation between morphological features of the cerebral aneurysm, and the mechanical stresses on the containing arterial walls are influenced by the hemodynamic forces. Present study investigates three different shapes, which have been identified correspondingly in patient-specific scenarios as well. The primary objective is to categorize the bifurcation aneurysms into standard shapes such as, spherical, beehive and pear-shaped, based on patient-specific clinical studies and further compare and contrast the model aneurysms with the patient specific configurations, for their hemodynamic factors as well as the attendant stresses on the wall. MethodsComputational fluid dynamic simulations are performed accounting for the fluid-structure interaction (FSI) effects between the flowing fluid and the containing vessel wall. Blood is assumed to be Newtonian, while the arterial walls are assumed to be linearly elastic. A commercial solver is used for performing detailed calculations. Hemodynamic and bio-mechanical rupture predictions are carried out for the three different shapes. Observations derived from the idealised simulations are compared and contrasted against their patient-specific counterparts. ResultsFrom detailed numerical simulations, it was observed that pear-shaped aneurysms exhibit large re-circulation bubble and flow stagnation zone, with higher residence time for the particles, which may lead to atherosclerotic lesions. Beehive shape allows for maximum flow into the aneurysmal sac with concentrated jet impinging on the dome, leading to high values of maximum WSS (MWSS) resulting in great propensity to form a secondary bleb. However, flow field inside a spherical aneurysm is found to be stable with fewer vortices, and nearly uniform distribution of wall stresses are observed though-out the sac, which perhaps signifies hemodynamically and bio-mechanically stable condition. ConclusionCategorizing patient-specific intracranial aneurysms into standard shapes viz, spherical, beehive and pear could generalize the process of prediction of hemodynamic and bio-mechanical rupture indicators. Comparative assessment of the flow field and stresses reported from the simulations on idealised models, with corresponding patient-specific simulations reveal that, these studies could aid in understanding the generalised shape dependence of hemodynamic and bio-mechanical behaviour of aneurysms.