ABSTRACT
The surgical creation of vascular accesses for renal failure patients provides an abnormally high flow rate conduit in the patient's upper arm vasculature that facilitates the hemodialysis treatment. These vascular accesses, however, are very often associated with complications that lead to access failure and thrombotic incidents, mainly due to excessive neointimal hyperplasia (NH) and subsequently stenosis. Development of a framework to monitor and predict the evolution of the venous system post access creation can greatly contribute to maintaining access patency. Computational fluid dynamics (CFD) has been exploited to inspect the non-homeostatic wall shear stress (WSS) distribution that is speculated to trigger NH in the patient cohort under investigation. Thereafter, CFD in liaison with a gradient-free shape optimization method has been employed to analyze the deformation modes of the venous system enduring non-physiological hemodynamics. It is observed that the optimally evolved shapes and their corresponding hemodynamics strive to restore the homeostatic state of the venous system to a normal, pre-surgery condition. It is concluded that a CFD-shape optimization coupling that seeks to regulate the WSS back to a well-defined physiological WSS target range can accurately predict the mode of patient-specific access failure.
ABSTRACT
Given the current emphasis on accurate computational fluid dynamics (CFD) modeling of cardiovascular flows, which incorporates realistic blood vessel geometries and cardiac waveforms, it is necessary to revisit the conventional wisdom regarding the influences of non-Newtonian effects. In this study, patient-specific reconstructed 3D geometries, whole blood viscosity data, and venous pulses postdialysis access surgery are used as the basis for the hemodynamic simulations of renal failure patients with native fistula access. Rheological analysis of the viscometry data initially suggested that the correct choice of constitutive relations to capture the non-Newtonian behavior of blood is important because the end-stage renal disease (ESRD) patient cohort under observation experience drastic variations in hematocrit (Hct) levels and whole blood viscosity throughout the hemodialysis treatment. For this purpose, various constitutive relations have been tested and implemented in CFD practice, namely Quemada and Casson. Because of the specific interest in neointimal hyperplasia and the onset of stenosis in this study, particular attention is placed on differences in nonhomeostatic wall shear stress (WSS) as that drives the venous adaptation process that leads to venous geometric evolution over time in ESRD patients. Surprisingly, the CFD results exhibit no major differences in the flow field and general flow characteristics of a non-Newtonian simulation and a corresponding identical Newtonian counterpart. It is found that the vein's geometric features and the dialysis-induced flow rate have far greater influence on the WSS distribution within the numerical domain.