MRI-Based Inflow Boundary Conditions for Patient Specific Fluid Structure Interaction Modeling of Abdominal Aortic Aneurysms

Abstract

Rupture of abdominal aortic aneurysm (AAA) is the 10th leading cause of death for men over age of 50 in US. The decision for surgical intervention is currently based on aneurysm diameter or its expansion rate. However, the use of these criteria for all patients is debatable. For example, small aneurysms do rupture or become symptomatic before reaching the critical diameter. Computationally predicted mechanical wall stress is considered a viable alternative criterion for rupture risk assessment. Hence, it is important to evaluate the effect of different modeling approaches on the accuracy of the predicated AAA wall stress. For computational solid stress (CSS) analysis or finite element analysis (FEA), a uniform static or transient intraluminal pressure is generally applied on the wall-lumen surface whereas in fluid-structure interaction (FSI) modeling the wall-lumen surface experiences transient and non-uniform fluid stress. An earlier comparison on idealized AAA models [1] revealed that static and transient CSS underestimate the peak wall stress (PWS) by an average 20–30% for variable wall thickness and 10% for uniform wall thickness when compared to fully coupled FSI. However, FSI-predicted stresses and strains were observed to be sensitive to inflow and outflow boundary conditions, warranting further study on a more accurate approach for FSI modeling. Though significant work has been performed on modeling outflow boundary conditions [2], studies on the sensitivity of computed stress or strain to the type of FSI inflow boundary condition is scarce [2–4]. We hypothesize that a FSI framework with a patient specific velocity boundary condition derived from magnetic resonance imaging (MRI) data applied to patient specific AAA geometry would provide better accuracy of PWS calculations compared to a FEA model. In this work, we present a framework where the AAA geometry is reconstructed from computed tomography (CT) images, on which FSI simulations were performed with inlet velocity components extracted from patient MR images of the abdominal aorta. Fully coupled FSI simulations were performed and results were compared with CSS simulations with uniform transient pressure boundary conditions.