Abstract
PURPOSE: Most 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 verification, 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 verifies the CFD results with patient medical data and assesses the co-localisation between CFD results and NIH regions at follow-up. METHODS: Contrast-enhanced computed tomography angiography images were used to segment the patient's 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-ω SST turbulence model. The CFD-derived flow waveform was compared with the patient's DUS image to ensure verification. 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. RESULTS: In 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. CONCLUSION: This 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.