Numerical modelling of the fluid-structure interaction in complex vascular geometries

A complex network of vessels is responsible for the transportation of blood throughout the body and back to the heart. Fluid mechanics and solid mechanics play a fundamental role in this transport phenomenon and are particularly suited for computer simulations. The latter may contribute to a better comprehension of the physiological processes and mechanisms leading to cardiovascular diseases, which are currently the leading cause of death in the western world. In case these computational models include patient-specific geometries and/or the interaction between the blood flow and the arterial wall, they become challenging to develop and to solve, increasing both the operator time and the computational time. This is especially true when the domain of interest involves vascular pathologies such as a local narrowing (stenosis) or a local dilatation (aneurysm) of the arterial wall. To overcome these issues of high operator times and high computational times when addressing the bio(fluid)mechanics of complex geometries, this PhD thesis focuses on the development of computational strategies which improve the generation and the accuracy of image-based, fluid-structure interaction (FSI) models. First, a robust procedure is introduced for the generation of hexahedral grids, which allows for local grid refinements and automation. Secondly, a straightforward algorithm is developed to obtain the prestress which is implicitly present in the arterial wall of a – by the blood pressure – loaded geometry at the moment of medical image acquisition. Both techniques are validated, applied to relevant cases, and finally integrated into a fluid-structure interaction model of an abdominal mouse aorta, based on in vivo measurements

FSI procedures for Civil Engineering Applications

Implementation of a FSI procedure involving the interaction of shell structures and fluids

Lattice-Boltzmann simulations of cerebral blood flow

Computational haemodynamics play a central role in the understanding of blood behaviour in the cerebral vasculature, increasing our knowledge in the onset of vascular diseases and their progression, improving diagnosis and ultimately providing better patient prognosis. Computer simulations hold the potential of accurately characterising motion of blood and its interaction with the vessel wall, providing the capability to assess surgical treatments with no danger to the patient. These aspects considerably contribute to better understand of blood circulation processes as well as to augment pre-treatment planning. Existing software environments for treatment planning consist of several stages, each requiring significant user interaction and processing time, significantly limiting their use in clinical scenarios. The aim of this PhD is to provide clinicians and researchers with a tool to aid in the understanding of human cerebral haemodynamics. This tool employs a high performance fluid solver based on the lattice-Boltzmann method (coined HemeLB), high performance distributed computing and grid computing, and various advanced software applications useful to efficiently set up and run patient-specific simulations. A graphical tool is used to segment the vasculature from patient-specific CT or MR data and configure boundary conditions with ease, creating models of the vasculature in real time. Blood flow visualisation is done in real time using in situ rendering techniques implemented within the parallel fluid solver and aided by steering capabilities; these programming strategies allows the clinician to interactively display the simulation results on a local workstation. A separate software application is used to numerically compare simulation results carried out at different spatial resolutions, providing a strategy to approach numerical validation. This developed software and supporting computational infrastructure was used to study various patient-specific intracranial aneurysms with the collaborating interventionalists at the National Hospital for Neurology and Neuroscience (London), using three-dimensional rotational angiography data to define the patient-specific vasculature. Blood flow motion was depicted in detail by the visualisation capabilities, clearly showing vortex fluid ow features and stress distribution at the inner surface of the aneurysms and their surrounding vasculature. These investigations permitted the clinicians to rapidly assess the risk associated with the growth and rupture of each aneurysm. The ultimate goal of this work is to aid clinical practice with an efficient easy-to-use toolkit for real-time decision support

Comparison of existing aneurysm models and their path forward

The two most important aneurysm types are cerebral aneurysms (CA) and abdominal aortic aneurysms (AAA), accounting together for over 80\% of all fatal aneurysm incidences. To minimise aneurysm related deaths, clinicians require various tools to accurately estimate its rupture risk. For both aneurysm types, the current state-of-the-art tools to evaluate rupture risk are identified and evaluated in terms of clinical applicability. We perform a comprehensive literature review, using the Web of Science database. Identified records (3127) are clustered by modelling approach and aneurysm location in a meta-analysis to quantify scientific relevance and to extract modelling patterns and further assessed according to PRISMA guidelines (179 full text screens). Beside general differences and similarities of CA and AAA, we identify and systematically evaluate four major modelling approaches on aneurysm rupture risk: finite element analysis and computational fluid dynamics as deterministic approaches and machine learning and assessment-tools and dimensionless parameters as stochastic approaches. The latter score highest in the evaluation for their potential as clinical applications for rupture prediction, due to readiness level and user friendliness. Deterministic approaches are less likely to be applied in a clinical environment because of their high model complexity. Because deterministic approaches consider underlying mechanism for aneurysm rupture, they have improved capability to account for unusual patient-specific characteristics, compared to stochastic approaches. We show that an increased interdisciplinary exchange between specialists can boost comprehension of this disease to design tools for a clinical environment. By combining deterministic and stochastic models, advantages of both approaches can improve accessibility for clinicians and prediction quality for rupture risk.Comment: 46 pages, 5 figure

Improving Cardiovascular Stent Design Using Patient-Specific Models and Shape Optimization

Stent geometry influences local hemodynamic alterations (i.e. the forces moving blood through the cardiovascular system) associated with adverse clinical outcomes. Computational fluid dynamics (CFD) is frequently used to quantify stent-induced hemodynamic disturbances, but previous CFD studies have relied on simplified device or vascular representations. Additionally, efforts to minimize stent-induced hemodynamic disturbances using CFD models often only compare a small number of possible stent geometries. This thesis describes methods for modeling commercial stents in patient-specific vessels along with computational techniques for determining optimal stent geometries that address the limitations of previous studies. An efficient and robust method was developed for virtually implanting stent models into patient-specific vascular geometries derived from medical imaging data. Models of commercial stent designs were parameterized to allow easy control over design features. Stent models were then virtually implanted into vessel geometries using a series of Boolean operations. This approach allowed stented vessel models to be automatically regenerated for rapid analysis of the contribution of design features to resulting hemodynamic alterations. The applicability of the method was demonstrated with patient-specific models of a stented coronary artery bifurcation and basilar trunk aneurysm to reveal how it can be used to investigate differences in hemodynamic performance in complex vascular beds for a variety of clinical scenarios. To identify hemodynamically optimal stents designs, a computational framework was constructed to couple CFD with a derivative-free optimization algorithm. The optimization algorithm was fully-automated such that solid model construction, mesh generation, CFD simulation and time-averaged wall shear stress (TAWSS) quantification did not require user intervention. The method was applied to determine the optimal number of circumferentially repeating stent cells (NC) for a slotted-tube stents and various commercial stents. Optimal stent designs were defined as those minimizing the area of low TAWSS. It was determined the optimal value of NC is dependent on the intrastrut angle with respect to the primary flow direction. Additionally, the geometries of current commercial stents were found to generally incorporate a greater NC than is hemodynamically optimal. The application of the virtual stent implantation and optimization methods may lead to stents with superior hemodynamic performance and the potential for improved clinical outcomes. Future in vivo studies are needed to validate the findings of the computational results obtained from the methods developed in this thesis