Open Access

ARTICLE

Patient-Specific Carotid Plaque Progression Simulation Using 3D Meshless Generalized Finite Difference Models with Fluid-Structure Interactions Based on Serial In Vivo MRI Data

Chun Yang^1,2, Dalin Tang², Satya Atluri³

1 Corresponding author, chunyang@bnu.edu.cn, School of Mathematical Sciences, Beijing Normal University, Key Laboratory of Mathematics and Complex Systems, Ministry of Education, Beijing,100875, China
2 Worcester Polytechnic Institute, Worcester, MA 01609
3 Center of Aerospace Research & Education, University of California, Irvine, CA 92612

Computer Modeling in Engineering & Sciences 2011, 72(1), 53-78. https://doi.org/10.3970/cmes.2011.072.053

Abstract

Previously, we introduced a computational procedure based on three-dimensional meshless generalized finite difference (MGFD) method and serial magnetic resonance imaging (MRI) data to quantify patient-specific carotid atherosclerotic plaque growth functions and simulate plaque progression. Structure-only models were used in our previous report. In this paper, fluid-stricture interaction (FSI) was added to improve on prediction accuracy. One participating patient was scanned three times (T1, T2, and T3, at intervals of about 18 months) to obtain plaque progression data. Blood flow was assumed to laminar, Newtonian, viscous and incompressible. The Navier-Stokes equations with arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations. Plaque material was assumed to be uniform, homogeneous, isotropic, linear, and nearly incompressible. The linear elastic model was used. The 3D FSI plaque model was discretized and solved using a meshless generalized finite difference (GFD) method. Growth functions with a) morphology alone; b) morphology and plaque wall stress (PWS); morphology and flow shear stress (FSS), and d) morphology, PWS and FSS were introduced to predict future plaque growth based on previous time point data. Starting from the T2 plaque geometry, plaque progression was simulated by solving the FSI model and adjusting plaque geometry using plaque growth functions iteratively until T3 is reached. Numerically simulated plaque progression agreed very well with the target T3 plaque geometry with errors ranging from 8.62%, 7.22%, 5.77% and 4.39%, with the growth function including morphology, plaque wall stress and flow shear stress terms giving the best predictions. Adding flow shear stress term to the growth function improved the prediction error from 7.22% to 4.39%, a 40% improvement. We believe this is the first time 3D plaque progression FSI simulation based on multi-year patient-tracking data was reported. Serial MRI-based progression simulation adds time dimension to plaque vulnerability assessment and will improve prediction accuracy for potential plaque rupture risk.

---

Keywords

---