A computational model for false lumen thrombosis in type B aortic dissection following thoracic endovascular repair

Thoracic endovascular repair (TEVAR) has recently been established as the preferred treatment option for complicated type B dissection. This procedure involves covering the primary entry tear to stimulate aortic remodelling and promote false lumen thrombosis thereby restoring true lumen flow. However, complications associated with incomplete false lumen thrombosis, such as aortic dilatation and stent graft induced new entry tears, can arise after TEVAR. This study presents the application and validation of a recently developed mathematical model for patient-specific prediction of thrombus formation and growth under physiologically realistic flow conditions. The model predicts thrombosis through the evaluation of shear rates, fluid residence time and platelet distribution, based on convection-diffusion-reaction transport equations. The model was applied to 3 type B aortic dissection patients: two TEVAR cases showing complete and incomplete false lumen thrombosis respectively, and one medically treated dissection with no signs of thrombosis. Predicted thrombus growth over time was validated against follow-up CT scans, showing good agreement with in vivo data in all cases with a maximum difference between predicted and measured false lumen reduction below 8%. Our results demonstrate that TEVAR-induced thrombus formation in type B aortic dissection can be predicted based on patient-specific anatomy and physiologically realistic boundary conditions. Our model can be used to identify anatomical or stent graft related factors that are associated with incomplete false lumen thrombosis following TEVAR, which may help clinicians develop personalised treatment plans for dissection patients in the future.     

构建基于患者个体化的Stanford B型主动脉夹层计算流体力学数值模拟分析模型

目的:探索简便、高效、精确的构建基于真实人体解剖形态结构的Stanford B型主动脉夹层计算流体力学数值模拟分析模型的方法.方法:利用Siemens Sensation Cardiac 64层螺旋CT薄层扫描技术,基于1mm层厚获取6例Stanford B型主动脉夹层连续断层DICOM格式图像,导入Materialise MIMICS v12.11软件,界定目标区域后生成三维动脉模型,经网格优化处理去除低质量及相交面网格,保存结果输出,导入TGrid 5.0软件,对主动脉面网格模型进行几何修复,使面网格扭曲率<0.75,采用自由分网方式生成Stanford B型主动脉夹层计算流体力学分析体网格模型,并对所构建模型进行血流属性、流场边界等界定,初步验证模型的有效性.结果:通过初步计算求解,确定所构建的6例Stanford B型主动脉夹层计算流体力学分析模型分别包含1857030,1820501,1844181,1849651,1858246及1814914个四面体单元.结论:利用64层螺旋CT薄层扫描技术获取DICOM格式连续断层CT图像可快速、准确地构建Stanford B型主动脉夹层计算流体力学数值模拟分析模型,为进一步的计算流体力学分析奠定了良好的基础.     

Numerical simulation of two-phase non-Newtonian blood flow with fluid-structure interaction in aortic dissection

The behavior of blood cells and vessel compliance significantly influence hemodynamic parameters, which are closely related to the development of aortic dissection. Here the two-phase non-Newtonian model and the fluid-structure interaction (FSI) method are coupled to simulate blood flow in a patient-specific dissected aorta. Moreover, three-element Windkessel model is applied to reproduce physiological pressure waves. Important hemodynamic indicators, such as the spatial distribution of red blood cells (RBCs) and vessel wall displacement, which greatly influence the hemodynamic characteristics are analyzed. Results show that the proximal false lumen near the entry tear appears to be a vortex zone with a relatively lower volume fraction of RBCs, a low time-averaged wall shear stress (TAWSS) and a high oscillatory shear index (OSI), providing a suitable physical environment for the formation of atherosclerosis. The highest TAWSS is located in the narrow area of the distal true lumen which might cause further dilation. TAWSS distributions in the FSI model and the rigid wall model show similar trend, while there is a significant difference for the OSI distributions. We suggest that an integrated model is essential to simulate blood flow in a more realistic physiological environment with the ultimate aim of guiding clinical treatment.     

Evaluation of 4D flow MRI-based non-invasive pressure assessment in aortic coarctations

Severity of aortic coarctation (CoA) is currently assessed by estimating trans-coarctation pressure drops through cardiac catheterization or echocardiography. In principle, more detailed information could be obtained non-invasively based on space- and time-resolved magnetic resonance imaging (4D flow) data. Yet the limitations of this imaging technique require testing the accuracy of 4D flow-derived hemodynamic quantities against other methodologies.With the objective of assessing the feasibility and accuracy of this non-invasive method to support the clinical diagnosis of CoA, we developed an algorithm (4DF-FEPPE) to obtain relative pressure distributions from 4D flow data by solving the Poisson pressure equation. 4DF-FEPPE was tested against results from a patient-specific fluid-structure interaction (FSI) simulation, whose patient-specific boundary conditions were prescribed based on 4D flow data. Since numerical simulations provide noise-free pressure fields on fine spatial and temporal scales, our analysis allowed to assess the uncertainties related to 4D flow noise and limited resolution.4DF-FEPPE and FSI results were compared on a series of cross-sections along the aorta. Bland-Altman analysis revealed very good agreement between the two methodologies in terms of instantaneous data at peak systole, end-diastole and time-averaged values: biases (means of differences) were +0.4 mmHg, −1.1 mmHg and +0.6 mmHg, respectively. Limits of agreement (2 SD) were ±0.978 mmHg, ±1.06 mmHg and ±1.97 mmHg, respectively. Peak-to-peak and maximum trans-coarctation pressure drops obtained with 4DF-FEPPE differed from FSI results by 0.75 mmHg and −1.34 mmHg respectively. The present study considers important validation aspects of non-invasive pressure difference estimation based on 4D flow MRI, showing the potential of this technology to be more broadly applied to the clinical practice.     

Validation of numerical flow simulations against in vitro phantom measurements in different type B aortic dissection scenarios

An aortic dissection (AD) is a serious condition defined by the splitting of the arterial wall, thus generating a secondary lumen [the false lumen (FL)]. Its management, treatment and follow-up are clinical challenges due to the progressive aortic dilatation and potentially severe complications during follow-up. It is well known that the direction and rate of dilatation of the artery wall depend on haemodynamic parameters such as the local velocity profiles, intra-luminal pressures and resultant wall stresses. These factors act on the FL and true lumen, triggering remodelling and clinical worsening. In this study, we aimed to validate a computational fluid dynamic (CFD) tool for the haemodynamic characterisation of chronic (type B) ADs. We validated the numerical results, for several dissection geometries, with experimental data obtained from a previous in vitro study performed on idealised dissected physical models. We found a good correlation between CFD simulations and experimental measurements as long as the tear size was large enough so that the effect of the wall compliance was negligible.     

Comparison between a generalized Newtonian model and a network-type multiscale model for hemodynamic behavior in the aortic arch: Validation with 4D MRI data for a case study

Blood is a complex fluid in which the presence of the various constituents leads to significant changes in its rheological properties. Thus, an appropriate non-Newtonian model is advisable; and we choose a Modified version of the rheological model of Phan-Thien and Tanner (MPTT). The different parameters of this model, derived from the rheology of polymers, allow characterization of the non-Newtonian nature of blood, taking into account the behavior of red blood cells in plasma. Using the MPTT model that we implemented in the open access software OpenFOAM, numerical simulations have been performed on blood flow in the thoracic aorta for a healthy patient. We started from a patient-specific model which was constructed from medical images. Exiting flow boundary conditions have been developped, based on a 3-element Windkessel model to approximate physiological conditions. The parameters of the Windkessel model were calibrated with in vivo measurements of flow rate and pressure. The influence of the selected viscosity of red blood cells on the flow and wall shear stress (WSS) was investigated. Results obtained from this model were compared to those of the Newtonian model, and to those of a generalized Newtonian model, as well as to in vivo dynamic data from 4D MRI during a cardiac cycle. Upon evaluating the results, the MPTT model shows better agreement with the MRI data during the systolic and diastolic phases than the Newtonian or generalized Newtonian model, which confirms our interest in using a complex viscoelastic model.