Heterogeneous Porous Media Simulation

Intracranial aneurysms are vascular disorders in which weakness in the wall of a cerebral artery or vein causes a localized dilation of the blood vessel. Flow diversion is an endovascular technique where a flow diverter stent is placed in the parent blood vessel to divert blood flow away from the aneurysm itself. Simulation by computational fluid dynamics is an attractive method to study flow diverters, particularly to model the small gaps between stent struts as a porous media. In many cases obstructions are not equal across the free medium and the porous one must be heterogeneous. Finite Volume Method solves numerical problems of computational fluid dynamics, splitting the region of interest in cells of small volumes. Porous media are usually modeled as a set of simulation cells described in a dictionary with constant porosity parameters (Homogeneous medium). An heterogeneous medium can be described as multiple homogeneous media, one by one. However, creating multiple homogeneous porous media is a tedious job if each simulation cell requires different parameters. Also, porous medium sets creates overheads on memory and processor load. The open source tool OpenFOAM is a open source C++ toolbox for field operations and partial differential equations solving using Finite Volume Method, including computational fluid dynamics. The tool is well prepared to describe heterogeneous fields. In this work, porous media coefficients are described as tensor fields. A steady state flow solver considering this fields is developed. The fidelity of the solver is then studied qualitatively and quantitatively.Fil: Dazeo, Nicolás Ignacio. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; ArgentinaFil: Dottori, Javier Alejandro. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; ArgentinaFil: Boroni, Gustavo Adolfo. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; ArgentinaFil: Larrabide, Ignacio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; Argentina. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; Argentin

Intracranial Stents Being Modeled as a Porous Medium: Flow Simulation in Stented Cerebral Aneurysms

Intracranial aneurysms may be treated by flow diverters, alternatively to stents and coils combination. Numerical simulation allows the assessment of the complex nature of aneurismal flow. Endovascular devices present a rather dense and fine strut network, increasing the complexity of the meshing. We propose an alternative strategy, which is based on the modeling of the device as a porous medium. Two patient-specific aneurysm data sets were reconstructed using conventional clinical setups. The aneurysms selection was done so that intra-aneurismal flow was shear driven in one and inertia driven in the other. Stents and their porous medium analog were positioned at the aneurysm neck. Physiological flow and standard boundary conditions were applied. The comparison between both approaches was done by analyzing the velocity, vorticity, and shear rate magnitudes inside the aneurysm as well as the wall shear stress (WSS) at the aneurysm surface. Simulations without device were also computed. The average flow reduction reaches 76 and 41% for the shear and inertia driven flow models, respectively. When comparing the two approaches, results show a remarkable similarity in the flow patterns and magnitude. WSS, iso-velocity surfaces and velocity on a trans-sectional plane are in fairly good agreement. The root mean squared error on the investigated parameters reaches 20% for aneurysm velocity, 30.6% for aneurysm shear rate, and 47.4% for aneurysm vorticity. It reaches 20.6% for WSS computed on the aneurysm surface. The advantages of this approach reside in its facility to implement and in the gain in computational time. Results predicted by the porous medium approach compare well with the real stent geometry model and allow predicting the main effects of the device on intra-aneurismal flow, facilitating thus the analysi

Hemodynamic investigation of the flow diverter treatment of intracranial aneurysm

Flow diverter stents (FDS) are increasingly used for the treatment of complex intracranial aneurysms such as fusiform, giant, or wide-neck aneurysms. The primary goal of these devices is to reconstruct the diseased vascular segment by diverting blood flow from the aneurysm. The resulting intra-aneurysmal flow reduction promotes progressive aneurysm thrombosis and healing of the disease. In the present study, a numerical investigation was performed for modeling blood flow inside a patient-specific intracranial aneurysm virtually treated with FDS. The aim of the study is to investigate the effects of FDS placement prior to the actual endovascular treatment and to compare the effectiveness of devices differing in porosity. Numerical simulations were performed under pulsatile flow conditions, taking into account the non-Newtonian behavior of blood. Two possible post-operative conditions with virtual stent deployment were simulated. Hemodynamic parameters were calculated and compared between the pre-operative (no stent placement) and post-operative (virtual stent placement) aneurysm models. FDS placement significantly reduced intra-aneurysmal flow velocity and increased the Relative Residence Time (RRT) on the aneurysm, thus promoting thrombus formation within the dilatation and aneurysm occlusion. The results highlighted an increase in the effectiveness of FDS as its porosity increased. The proposed analysis provides pre-operative knowledge on the impact of FDS on intracranial hemodynamics, allowing the selection of the most effective treatment for the specific patient

Stenting as porous media in anatomically accurate geometries: A comparison of models and spatial heterogeneity

Modelling intracranial aneurysm blood flow after flow diverter treatment has proven to be of great scientific and clinical interest. One of the reasons for not having CFD as an everyday clinical tool yet is the time required to set-up such simulations plus the required computational time. The speed-up of these simulations can have a considerable impact during treatment planning and device selection. Modelling flow diverters as a porous medium (PM) can considerably improve the computational time. Many models have been presented in literature, but quantitative comparisons between models are scarce.In this study, the untreated case, the explicit definition of the flow diverter wires as no-slip boundary condition and five different porous medium models were chosen for comparison, and evaluated on intracranial aneurysm of 14 patients with different shapes, sizes, and locations. CFD simulations were made using finite volume method on steady flow conditions. Velocities, kinetic energy, wall shear stress, and computational time were assessed for each model. Then, all models are compared against the no-slip boundary condition using non parametric Kolmogorov–Smirnov test.The model with least performance showed a mean K-S statistic of 0.31 and deviance of 0.2, while the model with best values always gave K-S statistics below 0.2. Kinetic energy between PM models varied between an over estimation of 218.3% and an under estimation of 73.06%. Also, speedups were between 4.75x and 5.3x (stdev: 0.38x and 0.15x) when using PM models.Flow diverters can be simulated with PM with a good agreement to standard CFD simulations were FD wires are represented with no-slip boundary condition in less than a quarter of the time. Best results were obtained on PM models based on geometrical properties, in particular, when using a heterogeneous medium based on equations for flat rhomboidal wire frames.Fil: Dazeo, Nicolás Ignacio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; Argentina. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; ArgentinaFil: Dottori, Javier Alejandro. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; Argentina. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; ArgentinaFil: Boroni, Gustavo Adolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; Argentina. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; ArgentinaFil: Narata, Ana Paula. Universite de Tours; FranciaFil: Larrabide, Ignacio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Tandil; Argentina. Universidad Nacional del Centro de la Provincia de Buenos Aires. Facultad de Ciencias Exactas. Grupo de Plasmas Densos Magnetizados. Provincia de Buenos Aires. Gobernación. Comision de Investigaciones Científicas. Grupo de Plasmas Densos Magnetizados; Argentin

Fluid-structure interaction modeling of artery aneurysms with steady-state configurations

This paper addresses numerical simulations of fluid-structure interaction (FSI) problems involving artery aneurysms, focusing on steady-state configurations. Both the fluid flow and the hyperelastic material are incompressible. A monolithic formulation for the FSI problem is considered, where the deformation of the fluid domain is taken into account according to an Arbitrary Lagrangian Eulerian (ALE) scheme. The numerical algorithm is a Newton-Krylov method combined with geometric multigrid preconditioner and smoothing based on domain decomposition. The system is modeled using a specific equation shuffling that aims at improving the row pivoting. Due to the complexity of the operators, the exact Jacobian matrix is evaluated using automatic differentiation tools. We describe benchmark settings which shall help to test and compare different numerical methods and code implementations for the FSI problem in hemodynamics. The configurations consist of realistic artery aneurysms. A case of endovascular stent implantation on a cerebral aneurysm is also presented. Hybrid meshes are employed in such configurations. We show numerical results for the described aneurysm geometries for steady-state boundary conditions. Parallel implementation is also addressed

In-silico clinical trials for assessment of intracranial flow diverters

In-silico trials refer to pre-clinical trials performed, entirely or in part, using individualised computer models that simulate some aspect of drug effect, medical device, or clinical intervention. Such virtual trials reduce and optimise animal and clinical trials, and enable exploring a wider range of anatomies and physiologies. In the context of endovascular treatment of intracranial aneurysms, in-silico trials can be used to evaluate the effectiveness of endovascular devices over virtual populations of patients with different aneurysm morphologies and physiologies. However, this requires (i) a virtual endovascular treatment model to evaluate device performance based on a reliable performance indicator, (ii) models that represent intra- and inter-subject variations of a virtual population, and (iii) creation of cost-effective and fully-automatic workflows to enable a large number of simulations at a reasonable computational cost and time. Flow-diverting stents have been proven safe and effective in the treatment of large wide-necked intracranial aneurysms. The presented thesis aims to provide the ingredient models of a workflow for in-silico trials of flow-diverting stents and to enhance the general knowledge of how the ingredient models can be streamlined and accelerated to allow large-scale trials. This work contributed to the following aspects: 1) To understand the key ingredient models of a virtual treatment workflow for evaluation of the flow-diverter performance. 2) To understand the effect of input uncertainty and variability on the workflow outputs, 3) To develop generative statistical models that describe variability in internal carotid artery flow waveforms, and investigate the effect of uncertainties on quantification of aneurysmal wall shear stress, 4) As part of a metric to evaluate success of flow diversion, to develop and validate a thrombosis model to assess FD-induced clot stability, and 5) To understand how a fully-automatic aneurysm flow modelling workflow can be built and how computationally inexpensive models can reduce the computational costs

Intracranial Stents Being Modeled as a Porous Medium: Flow Simulation in Stented Cerebral Aneurysms

Intracranial aneurysms may be treated by flow diverters, alternatively to stents and coils combination. Numerical simulation allows the assessment of the complex nature of aneurismal flow. Endovascular devices present a rather dense and fine strut network, increasing the complexity of the meshing. We propose an alternative strategy, which is based on the modeling of the device as a porous medium. Two patient-specific aneurysm data sets were reconstructed using conventional clinical setups. The aneurysms selection was done so that intra-aneurismal flow was shear driven in one and inertia driven in the other. Stents and their porous medium analog were positioned at the aneurysm neck. Physiological flow and standard boundary conditions were applied. The comparison between both approaches was done by analyzing the velocity, vorticity, and shear rate magnitudes inside the aneurysm as well as the wall shear stress (WSS) at the aneurysm surface. Simulations without device were also computed. The average flow reduction reaches 76 and 41% for the shear and inertia driven flow models, respectively. When comparing the two approaches, results show a remarkable similarity in the flow patterns and magnitude. WSS, iso-velocity surfaces and velocity on a trans-sectional plane are in fairly good agreement. The root mean squared error on the investigated parameters reaches 20% for aneurysm velocity, 30.6% for aneurysm shear rate, and 47.4% for aneurysm vorticity. It reaches 20.6% for WSS computed on the aneurysm surface. The advantages of this approach reside in its facility to implement and in the gain in computational time. Results predicted by the porous medium approach compare well with the real stent geometry model and allow predicting the main effects of the device on intra-aneurismal flow, facilitating thus the analysis

Virtual endovascular treatment of intracranial aneurysms: models and uncertainty

Virtual endovascular treatment models (VETMs) have been developed with the view to aid interventional neuroradiologists and neurosurgeons to pre-operatively analyze the comparative efficacy and safety of endovascular treatments for intracranial aneurysms. Based on the current state of VETMs in aneurysm rupture risk stratification and in patient-specific prediction of treatment outcomes, we argue there is a need to go beyond personalized biomechanical flow modeling assuming deterministic parameters and error-free measurements. The mechanobiological effects associated with blood clot formation are important factors in therapeutic decision making and models of post-treatment intra-aneurysmal biology and biochemistry should be linked to the purely hemodynamic models to improve the predictive power of current VETMs. The influence of model and parameter uncertainties associated to each component of a VETM is, where feasible, quantified via a random-effects meta-analysis of the literature. This allows estimating the pooled effect size of these uncertainties on aneurysmal wall shear stress. From such meta-analyses, two main sources of uncertainty emerge where research efforts have so far been limited: (1) vascular wall distensibility, and (2) intra/intersubject systemic flow variations. In the future, we suggest that current deterministic computational simulations need to be extended with strategies for uncertainty mitigation, uncertainty exploration, and sensitivity reduction techniques. WIREs Syst Biol Med 2017, 9:e1385. doi: 10.1002/wsbm.138