Patient-specific hemodynamic markers co-localise with regions of neointimal hyperplasia in an arteriovenous graft

PurposeMost computational fluid dynamics (CFD) studies on arteriovenous grafts (AVGs) adopt idealised geometries and simplified boundary conditions (BCs), potentially resulting in misleading conclusions when attempting to predict neointimal hyperplasia (NIH) development. Moreover, they often analyse a limited range of hemodynamic indices, lack validation, and fail to link the graft-altered hemodynamics with follow-up data. This study develops a novel patient-specific CFD workflow for AVGs using pathophysiological BCs. It validates the CFD results with patient medical data and assesses the co-localisation between CFD results and NIH regions at follow-up. MethodsContrast-enhanced computed tomography angiography images were used to segment the patients AVG geometry. A uniform Doppler ultrasound (DUS)-derived velocity profile was imposed at the inlet, and three-element Windkessel models were applied at the arterial outlets of the domain. Transient, rigid-wall simulations were performed using the k-{omega} SST turbulence model. The CFD-derived flow waveform was compared with the patients DUS image to ensure validation. Turbulent kinetic energy (TKE), helicity and near-wall hemodynamic descriptors were calculated and linked with regions presenting NIH from a 4-month follow-up fistulogram. ResultsIn the analysed patient, areas presenting high TKE and balanced helical flow structures at baseline exhibit NIH growth at follow-up. Transverse wall shear stress index is a stronger predictor of NIH than other commonly analysed near-wall hemodynamic indices, since luminal areas subjected to high values greatly co-localise with observed areas of remodelling. ConclusionThis patient-specific computational workflow for AVGs could be applied to a larger cohort to unravel the link between altered hemodynamics and NIH progression in vascular access.