Search

The trend towards realistic numerical models of (pathologic) patient-specific vascular structures brings along larger computational domains and more complex geometries, increasing both the computation time and the operator time. Hexahedral grids effectively lower the computational run time and the required computational infrastructure, but at high cost in terms of operator time and minimal cell quality, especially when the computational analyses are targeting complex geometries such as aneurysm necks, severe stenoses and bifurcations. Moreover, such grids generally do not allow local refinements. As an attempt to overcome these limitations, a novel approach to hexahedral meshing is proposed in this paper, which combines the automated generation of multi-block structures with a grid-based method. The robustness of the novel approach is tested on common complex geometries, such as tree-like structures (including trifurcations), stenoses, and aneurysms. Additionally, the performance of the generated grid is assessed using two numerical examples. In the first example, a grid sensitivity analysis is performed for blood flow simulated in an abdominal mouse aorta and compared to tetrahedral grids with a prismatic boundary layer. In the second example, the fluid-structure interaction in a model of an aorta with aortic coarctation is simulated and the effect of local grid refinement is analyzed.

Carotid artery stenting (CAS) has emerged as a minimally invasive alternative to endarterectomy but its use in clinical treatment is limited due to the post-stenting complications. Haemodynamic actors, related to blood flow in the stented vessel, have been suggested to play a role in the endothelium response to stenting, including adverse reactions such as in-stent restenosis and late thrombosis. Accessing the flow-related shear forces acting on the endothelium in vivo requires space and time resolutions which are currently not achievable with non-invasive clinical imaging techniques but can be obtained from image-based computational analysis. In this study, we present a framework for accurate determination of the wall shear stress (WSS) in a mildly stenosed carotid artery after the implantation of a stent, resembling the commercially available Acculink (Abbott Laboratories, Abbott Park, Illinois, USA). Starting from angiographic CT images of the vessel lumen and a micro-CT scan of the stent, a finite element analysis is carried out in order to deploy the stent in the vessel, reproducing CAS in silico. Then, based on the post-stenting anatomy, the vessel is perfused using a set of boundary conditions: total pressure is applied at the inlet, and impedances that are assumed to be insensitive to the presence of the stent are imposed at the outlets. Evaluation of the CAS outcome from a geometrical and haemodynamic perspective shows the presence of atheroprone regions (low time-average WSS, high relative residence time) colocalised with stent malapposition and stent strut interconnections. Stent struts remain unapposed in the ostium of the external carotid artery disturbing the flow and generating abnormal shear forces, which could trigger thromboembolic events.

Recent numerical studies of abdominal aortic aneurysm (AAA) suggest that intraluminal thrombus (ILT) may reduce the stress loading on the aneurysmal wall. Detailed fluid structure interaction (FSI) in the presence and absence of ILT may help predict AAA rupture risk better. Two patients, with varied AAA geometries and ILT structures, were studied and compared in detail. The patient specific 3D geometries were reconstructed from CT scans, and uncoupled FSI approach was applied. Complex flow trajectories within the AAA lumen indicated a viable mechanism for the formation and growth of the ILT. The resulting magnitude and location of the peak wall stresses was dependent on the shape of the AAA, and the ILT appeared to reduce wall stresses for both patients. Accordingly, the inclusion of ILT in stress analysis of AAA is of importance and would likely increase the accuracy of predicting AAA risk of rupture.