Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D flow-MRI

MacDonald Black, Scott and Maclean, Craig and Hall Barrientos, Pauline and Ritos, Konstantinos and McQueen, Alistair and Kazakidi, Asimina (2023) Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D flow-MRI. Frontiers in Bioengineering and Biotechnology, 11. 1178483. ISSN 2296-4185 (https://doi.org/10.3389/fbioe.2023.1178483)

[thumbnail of Black-etal-FBB-2023-Calibration-of-patient-specific-boundary-conditions-for-coupled-CFD-models-of-the-aorta]
Preview
Text. Filename: Black_etal_FBB_2023_Calibration_of_patient_specific_boundary_conditions_for_coupled_CFD_models_of_the_aorta.pdf
Final Published Version
License: Creative Commons Attribution 4.0 logo

Download (3MB)| Preview

Abstract

Introduction: Patient-specific computational fluid dynamics (CFD) models permit analysis of complex intra-aortic hemodynamics in patients with aortic dissection (AD), where vessel morphology and disease severity are highly individualized. The simulated blood flow regime within these models is sensitive to the prescribed boundary conditions (BCs), so accurate BC selection is fundamental to achieve clinically relevant results. Methods: This study presents a novel reduced-order computational framework for the iterative flow-based calibration of 3-Element Windkessel Model (3EWM) parameters to generate patient-specific BCs. These parameters were calibrated using time-resolved flow information derived from retrospective four-dimensional flow magnetic resonance imaging (4D Flow-MRI). For a healthy and dissected case, blood flow was then investigated numerically in a fully coupled zero dimensional-three dimensional (0D-3D) numerical framework, where the vessel geometries were reconstructed from medical images. Calibration of the 3EWM parameters was automated and required ~3.5 min per branch. Results: With prescription of the calibrated BCs, the computed near-wall hemodynamics (time-averaged wall shear stress, oscillatory shear index) and perfusion distribution were consistent with clinical measurements and previous literature, yielding physiologically relevant results. BC calibration was particularly important in the AD case, where the complex flow regime was captured only after BC calibration. Discussion: This calibration methodology can therefore be applied in clinical cases where branch flow rates are known, for example, via 4D Flow-MRI or ultrasound, to generate patient-specific BCs for CFD models. It is then possible to elucidate, on a case-by-case basis, the highly individualized hemodynamics which occur due to geometric variations in aortic pathology high spatiotemporal resolution through CFD.