Modelling of the arteriovenous malformation embolization optimal scenario

Rheological Properties of Non-Adhesive Embolizing Compounds-The Key to Fine-Tuning Embolization Process-Modeling in Endovascular Surgery

The study of polymers' rheological properties is of paramount importance both for the problems of their industrial production as well as for their practical application. Two polymers used for embolization of arteriovenous malformations (AVMs) are studied in this work: Onyx-18^® and Squid-12^®. Viscosity curve tests and computational fluid dynamics (CFD) were used to uncover viscosity law as a function of shear rate as well as behavior of the polymers in catheter or pathological tissue models. The property of thermal activation of viscosity was demonstrated, namely, the law of dependence of viscosity on temperature in the range from 20 °C to 37 °C was established. A zone of viscosity nonmonotonicity was identified, and a physical interpretation of the dependence of the embolic polymers' viscosity on the shear rate was given on the basis of Cisco's model. The obtained empirical constants will be useful for researchers based on the CFD of AVMs. A description of the process of temperature activation of the embolic polymers' viscosity is important for understanding the mechanics of the embolization process by practicing surgeons as well as for producing new prospective embolic agents.

A Computational Framework for Pre-Interventional Planning of Peripheral Arteriovenous Malformations

PURPOSE

Peripheral arteriovenous malformations (pAVMs) are congenital lesions characterised by abnormal high-flow, low-resistance vascular connections-the so-called nidus-between arteries and veins. The mainstay treatment typically involves the embolisation of the nidus, however the complexity of pAVMs often leads to uncertain outcomes. This study aims at developing a simple, yet effective computational framework to aid the clinical decision making around the treatment of pAVMs using routinely acquired clinical data.

METHODS

A computational model was developed to simulate the pre-, intra-, and post-intervention haemodynamics of a patient-specific pAVM. A porous medium of varying permeability was employed to simulate the sclerosant effect on the nidus haemodynamics. Results were compared against clinical data (digital subtraction angiography, DSA, images) and experimental flow-visualization results in a 3D-printed phantom of the same pAVM.

RESULTS

The computational model allowed the simulation of the pAVM haemodynamics and the sclerotherapy-induced changes at different interventional stages. The predicted inlet flow rates closely matched the DSA-derived data, although the post-intervention one was overestimated, probably due to vascular system adaptations not accounted for numerically. The nidus embolization was successfully captured by varying the nidus permeability and increasing its hydraulic resistance from 0.330 to 3970 mmHg s ml-1. The nidus flow rate decreased from 71% of the inlet flow rate pre-intervention to 1%: the flow completely bypassed the nidus post-intervention confirming the success of the procedure.

CONCLUSION

The study demonstrates that the haemodynamic effects of the embolisation procedure can be simulated from routinely acquired clinical data via a porous medium with varying permeability as evidenced by the good qualitative agreement between numerical predictions and both in vivo and in vitro data. It provides a fundamental building block towards a computational treatment-planning framework for AVM embolisation.

Optimal control problem arising in mathematical modeling of cerebral vascular pathology embolization

Arteriovenous malformation (AVM) of the brain is a congenital vascular abnormality, in which the arterial and venous blood pools are intertwined and directly connected. This dangerous disease causes a high risk of intracranial hemorrhage and disrupts brain functioning. The preferred method of AVM treating is embolization, which is the endovascular filling of abnormal AVM vessels with a special embolic agent. Despite the fact that this method is widely used in neurosurgery, in some cases its use is accompanied by perioperative AVM vessels rupture. In this regard, the aim of this work is to study the optimal scenarios for multi-stage AVM embolization from the effectiveness and safety of the procedure point of view. Mathematically, the joint movement of blood and embolic agent in the AVM body is described on the basis of a one-dimensional two-phase filtration model, which takes into account the redistribution of blood to surrounding healthy vessels. For the numerical solution of the resulting integro-differential system of equations, a monotonic modification of the CABARET scheme is used. To find optimal embolization scenarios, the optimal control problem with phase constraints arising from medicine is formulated. A modified particle swarm optimization method is used to solve this problem numerically. This technique is used to obtain optimal embolization scenarios on the basis of real patients clinical data collected during neurosurgical operations.