1
|
Zorrilla R, Soudah E. An efficient procedure for the blood flow computer simulation of patient-specific aortic dissections. Comput Biol Med 2024; 179:108832. [PMID: 39002313 DOI: 10.1016/j.compbiomed.2024.108832] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Revised: 06/06/2024] [Accepted: 06/29/2024] [Indexed: 07/15/2024]
Abstract
In this work we present a novel methodology for the numerical simulation of patient-specific aortic dissections. Our proposal, which targets the seamless virtual prototyping of customized scenarios, combines an innovative two-step segmentation procedure with a CutFEM technique capable of dealing with thin-walled bodies such as the intimal flap. First, we generate the fluid mesh from the outer aortic wall disregarding the intimal flap, similarly to what would be done in a healthy aorta. Second, we create a surface mesh from the approximate midline of the intimal flap. This approach allows us to decouple the segmentation of the fluid volume from that of the intimal flap, thereby bypassing the need to create a volumetric mesh around a thin-walled body, an operation widely known to be complex and error-prone. Once the two meshes are obtained, the original configuration of the dissection into true and false lumen is recovered by embedding the surface mesh into the volumetric one and calculating a level set function that implicitly represents the intimal flap in terms of the volumetric mesh entities. We then leverage the capabilities of unfitted mesh methods, specifically relying on a CutFEM technique tailored for thin-walled bodies, to impose the wall boundary conditions over the embedded intimal flap. We tested the method by simulating the flow in four patient-specific aortic dissections, all involving intricate geometrical patterns. In all cases, the preprocess is greatly simplified with no impact on the computational times. Additionally, the obtained results are consistent with clinical evidence and previous research.
Collapse
Affiliation(s)
- Rubén Zorrilla
- Departament d'Enginyeria Civil i Ambiental, Universitat Politècnica de Catalunya (UPC), Barcelona, 08034, Spain; International Center for Numerical Methods in Engineering (CIMNE), Barcelona, 08034, Spain.
| | - Eduardo Soudah
- Departament de Resistència de Materials i Estructures a l'Enginyeria, Universitat Politècnica de Catalunya (UPC), Barcelona, 08034, Spain; International Center for Numerical Methods in Engineering (CIMNE), Barcelona, 08034, Spain; Departamento de Ciencia de los Materiales e Ingeniería Metalúrgica, Expresión Gráfica en la Ingeniería, Ingeniería Cartográfica, Geodésica y Fotogrametría, Ingeniería Mecánica e Ingeniería de los Procesos de Fabricación, Universidad de Valladolid (UVA), Valladolid, 47011, Spain.
| |
Collapse
|
2
|
Li M, Ma T, Cai Y, Li J, Meng Z, Dong Z, Wang S. Numerical simulation of the distal stent graft-induced new entry after TEVAR. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2024; 40:e3819. [PMID: 38551141 DOI: 10.1002/cnm.3819] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Revised: 02/16/2024] [Accepted: 03/19/2024] [Indexed: 05/15/2024]
Abstract
The study aimed to investigate the mechanical factors for distal stent graft-induced new entry (dSINE) in aortic dissection patients and discussed these factors in conjunction with aortic morphology. Two patients (one dSINE and one non-dSINE), with the same age, gender, and type of implanted stent, were selected, then aortic morphological parameters were calculated. In addition, the stent material parameters used by the patients were also fitted. Simulations were performed based on the patient's aortic model and the stent graft used. The true lumen segment at the distal stent graft was designated as the "dSINE risk zone," and mechanical parameters (maximum principal strain, maximum principal stress) were computed. When approaching the area with higher mechanical parameters in the dSINE risk zone, dSINE patient exhibited higher values and growth rates in mechanical parameters compared to non-dSINE patient. Furthermore, dSINE patient also presented larger aortic taper ratio, stent oversizing ratio, and expansion mismatch ratio of the distal true lumen (EMRDTR). The larger mechanical parameters and growth rates in dSINE patient corresponded to a greater aortic taper ratio, stent oversizing ratio, and EMRDTR. The failure of dSINE prediction by the stent tortuosity index indicated that mechanical parameters were the fundamental reasons for dSINE development.
Collapse
Affiliation(s)
- Meixuan Li
- Department of Aeronautics and Astronautics, Institute of Biomechanics, Fudan University, Shanghai, China
| | - Tao Ma
- Department of Vascular Surgery, Zhongshan Hospital, Fudan University, and National Clinical Research Center for Interventional Medicine, Shanghai, China
| | - Yunhan Cai
- Department of Aeronautics and Astronautics, Institute of Biomechanics, Fudan University, Shanghai, China
| | - Jianming Li
- Department of Aeronautics and Astronautics, Institute of Biomechanics, Fudan University, Shanghai, China
| | - Zhuangyuan Meng
- Department of Aeronautics and Astronautics, Institute of Biomechanics, Fudan University, Shanghai, China
| | - Zhihui Dong
- Department of Vascular Surgery, Zhongshan Hospital, Fudan University, and National Clinical Research Center for Interventional Medicine, Shanghai, China
| | - Shengzhang Wang
- Department of Aeronautics and Astronautics, Institute of Biomechanics, Fudan University, Shanghai, China
- Institute of Biomedical Engineering Technology, Academy of Engineering and Technology, Fudan University, Shanghai, China
| |
Collapse
|
3
|
Zhu Y, Mirsadraee S, Rosendahl U, Pepper J, Xu XY. Fluid-Structure Interaction Simulations of Repaired Type A Aortic Dissection: a Comprehensive Comparison With Rigid Wall Models. Front Physiol 2022; 13:913457. [PMID: 35774287 PMCID: PMC9237394 DOI: 10.3389/fphys.2022.913457] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 05/17/2022] [Indexed: 11/26/2022] Open
Abstract
This study aimed to evaluate the effect of aortic wall compliance on intraluminal hemodynamics within surgically repaired type A aortic dissection (TAAD). Fully coupled two-way fluid-structure interaction (FSI) simulations were performed on two patient-specific post-surgery TAAD models reconstructed from computed tomography angiography images. Our FSI model incorporated prestress and different material properties for the aorta and graft. Computational results, including velocity, wall shear stress (WSS) and pressure difference between the true and false lumen, were compared between the FSI and rigid wall simulations. It was found that the FSI model predicted lower blood velocities and WSS along the dissected aorta. In particular, the area exposed to low time-averaged WSS ( ≤ 0.2 P a ) was increased from 21 cm2 (rigid) to 38 cm2 (FSI) in patient 1 and from 35 cm2 (rigid) to 144 cm2 (FSI) in patient 2. FSI models also produced more disturbed flow where much larger regions presented with higher turbulence intensity as compared to the rigid wall models. The effect of wall compliance on pressure difference between the true and false lumen was insignificant, with the maximum difference between FSI and rigid models being less than 0.25 mmHg for the two patient-specific models. Comparisons of simulation results for models with different Young's moduli revealed that a more compliant wall resulted in further reduction in velocity and WSS magnitudes because of increased displacements. This study demonstrated the importance of FSI simulation for accurate prediction of low WSS regions in surgically repaired TAAD, but a rigid wall computational fluid dynamics simulation would be sufficient for prediction of luminal pressure difference.
Collapse
Affiliation(s)
- Yu Zhu
- Department of Chemical Engineering, Imperial College London, London, United Kingdom
| | - Saeed Mirsadraee
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
- Department of Radiology, Royal Brompton and Harefield Hospitals, London, United Kingdom
| | - Ulrich Rosendahl
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
- Department of Cardiac Surgery, Royal Brompton and Harefield Hospitals, London, United Kingdom
| | - John Pepper
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
- Department of Cardiac Surgery, Royal Brompton and Harefield Hospitals, London, United Kingdom
| | - Xiao Yun Xu
- Department of Chemical Engineering, Imperial College London, London, United Kingdom
| |
Collapse
|
4
|
Stokes C, Bonfanti M, Li Z, Xiong J, Chen D, Balabani S, Díaz-Zuccarini V. A novel MRI-based data fusion methodology for efficient, personalised, compliant simulations of aortic haemodynamics. J Biomech 2021; 129:110793. [PMID: 34715606 PMCID: PMC8907869 DOI: 10.1016/j.jbiomech.2021.110793] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 08/24/2021] [Accepted: 09/30/2021] [Indexed: 01/24/2023]
Abstract
We present a novel, cost-efficient methodology to simulate aortic haemodynamics in a patient-specific, compliant aorta using an MRI data fusion process. Based on a previously-developed Moving Boundary Method, this technique circumvents the high computational cost and numerous structural modelling assumptions required by traditional Fluid-Structure Interaction techniques. Without the need for Computed Tomography (CT) data, the MRI images required to construct the simulation can be obtained during a single imaging session. Black Blood MR Angiography and 2D Cine-MRI data were used to reconstruct the luminal geometry and calibrate wall movement specifically to each region of the aorta. 4D-Flow MRI and non-invasive pressure measurements informed patient-specific inlet and outlet boundary conditions. Luminal area closely matched 2D Cine-MRI measurements with a mean error of less than 4.6% across the cardiac cycle, while physiological pressure and flow distributions were simulated to within 3.3% of patient-specific targets. Moderate agreement with 4D-Flow MRI velocity data was observed. Despite lower peak velocity, an equivalent rigid-wall simulation predicted a mean Time-Averaged Wall Shear Stress (TAWSS) 13% higher than the compliant simulation. The agreement observed between compliant simulation results and MRI data is testament to the accuracy and efficiency of this MRI-based simulation technique.
Collapse
Affiliation(s)
- Catriona Stokes
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| | - Mirko Bonfanti
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| | - Zeyan Li
- School of Life Science, Beijing Institute of Technology, Beijing, China.
| | - Jiang Xiong
- Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China.
| | - Duanduan Chen
- School of Life Science, Beijing Institute of Technology, Beijing, China.
| | - Stavroula Balabani
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| | - Vanessa Díaz-Zuccarini
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| |
Collapse
|
5
|
Li D, Zheng T, Liu Z, Li Y, Yuan D, Fan Y. Influence of Distal Re-entry Tears on False Lumen Thrombosis After Thoracic Endovascular Aortic Repair in Type B Aortic Dissection Patients: A Computational Fluid Dynamics Simulation. Cardiovasc Eng Technol 2021; 12:426-437. [PMID: 33768445 DOI: 10.1007/s13239-021-00532-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Accepted: 03/05/2021] [Indexed: 02/05/2023]
Abstract
PURPOSE Distal re-entry tears play a significant role in false lumen (FL) thrombosis, which will strongly affect the postoperative long-term survival of patients with type B aortic dissection (TBAD) after thoracic endovascular aortic repair (TEVAR). This study aimed to investigate the influence of a peculiar morphological parameter of the residual re-entry tears in TBAD patients after TEVAR on long-term FL thrombosis using the computational fluid dynamics. METHODS Ideal population-based three-dimensional models of post-operative TBAD were established. Numerical simulation was performed to investigate the hemodynamic differences caused by different tear features, including the tear count, the maximum distance between tears, and the tear area. RESULTS Although the low relative residence time (RRT) area did not change significantly when the tear distance was fixed, the area of oscillatory shear index (OSI) > 0.45 and endothelial cell activation potential (ECAP) > 1.5 decreased significantly with the tear count and area increased and a dramatic increase in blood flow into the FL was also observed. When tear count and total area were fixed, for each 10-mm increase in the maximum distance between tears, the area of low RRT in the FL increased significantly, while the average pressure difference increased by 10.85%. CONCLUSION The different morphology of the re-entry tears had different effects on the thrombosis-related hemodynamic parameters in FL following TEVAR. and the number of re-entry tears was most crucial to the potential thrombosis in the post-TEVAR FL of TBAD patients.
Collapse
Affiliation(s)
- Da Li
- Department of Applied Mechanics, Sichuan University, Chengdu, 610065, China
| | - Tinghui Zheng
- Department of Applied Mechanics, Sichuan University, Chengdu, 610065, China
| | - Zhan Liu
- Department of Applied Mechanics, Sichuan University, Chengdu, 610065, China
| | - Yan Li
- Department of Applied Mechanics, Sichuan University, Chengdu, 610065, China
| | - Ding Yuan
- Department Vascular Surgery of West China Hospital, Sichuan University, No. 37 Guo Xue Xiang, Chengdu, 610041, China.
| | - Yubo Fan
- Beijing Advanced Innovation Center of Biomedical Engineering, Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, Beijing, 100191, China.
| |
Collapse
|
6
|
Spatial Configuration of Abdominal Aortic Aneurysm Analysis as a Useful Tool for the Estimation of Stent-Graft Migration. Diagnostics (Basel) 2020; 10:diagnostics10100737. [PMID: 32977588 PMCID: PMC7598279 DOI: 10.3390/diagnostics10100737] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2020] [Revised: 09/11/2020] [Accepted: 09/21/2020] [Indexed: 11/17/2022] Open
Abstract
The aim of this study was to prepare a self-made mathematical algorithm for the estimation of risk of stent-graft migration with the use of data on abdominal aortic aneurysm (AAA) size and geometry of blood flow through aneurysm sac before or after stent-graft implantation. AngioCT data from 20 patients aged 50–60 years, before and after stent-graft placement in the AAA was analyzed. In order to estimate the risk of stent-graft migration for each patient we prepared an opposite spatial configuration of virtually reconstructed stent-graft with long body or short body. Thus, three groups of 3D geometries were analyzed: 20 geometries representing 3D models of aneurysm, 20 geometries representing 3D models of long body stent-grafts, and 20 geometries representing 3D models of short body stent-graft. The proposed self-made algorithm demonstrated its efficiency and usefulness in estimating wall shear stress (WSS) values. Comparison of the long or short type of stent-graft with AAA geometries allowed to analyze the implants’ spatial configuration. Our study indicated that short stent-graft, after placement in the AAA sac, generated lower drug forces compare to the long stent-graft. Each time shape factor was higher for short stent-graft compare to long stent-graft.
Collapse
|
7
|
The effect of the entry and re-entry size in the aortic dissection: a two-way fluid-structure interaction simulation. Biomech Model Mechanobiol 2020; 19:2643-2656. [PMID: 32621161 DOI: 10.1007/s10237-020-01361-0] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Accepted: 06/20/2020] [Indexed: 12/16/2022]
Abstract
Aortic dissection (AD) is one of the most catastrophic cardiovascular diseases. AD occurs when a layer inside the aorta is disrupted and gives rise to the formation of a true lumen and a false lumen. These lumens can be connected through tears in the intimal flap which are known as entries. Despite being known for about two centuries, the effects of many factors on the morbidity and mortality of this disease are still unknown. As the blood interaction with the aorta is crucial in the severity and the progression of the aortic dissection, a biomechanical approach is chosen to investigate the influence of different morphologies on the severity of this disease. Using the finite element method (FEM) and the fluid-structure interaction (FSI) approach, we have evaluated the blood flow characteristics along the diseased aorta, in conjunction with the deformation of the aortic wall. In this study, an idealized geometry of a dissected descending aorta (type B) with two entries has been studied. The values for the diameter of the entry tear were chosen to be 5 mm and 10 mm. Therefore, a total of four conditions were investigated. According to our results, the retrograde flow through the proximal tear is dependent on the size of the distal re-entry and vice versa. Our results revealed that when both entry and re-entry tears are 10 mm in diameter, the flow passes through the true and false lumens with smaller resistance, resulting in a smaller flutter of the intimal flap, and therefore more stable intimal flap. Major oscillation frequencies of 2.5 Hz and 7.4 Hz were observed for the oscillation of the intimal flap, and amplitudes of the waves with higher frequencies were negligible.
Collapse
|
8
|
Shape and Enhancement Analysis as a Useful Tool for the Presentation of Blood Hemodynamic Properties in the Area of Aortic Dissection. J Clin Med 2020; 9:jcm9051330. [PMID: 32370301 PMCID: PMC7290319 DOI: 10.3390/jcm9051330] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2020] [Revised: 04/19/2020] [Accepted: 04/28/2020] [Indexed: 11/17/2022] Open
Abstract
The aim of this study was to create a mathematical approach for blood hemodynamic description with the use of brightness analysis. Medical data was collected from three male patients aged from 45 to 65 years with acute type IIIb aortic dissection that started proximal to the left subclavian artery and involved the renal arteries. For the recognition of wall dissection areas Digital Imaging and Communications in Medicine (DICOM) data were applied. The distance from descending aorta to the diaphragm was analyzed. Each time Feret (DF) and Hydraulic (DHy) diameter were calculated. Moreover, an average brightness (BAV) was analyzed. Finally, to describe blood hemodynamic in the area of aortic wall dissection, mathematical function combining difference in brightness value and diameter for each computed tomography (CT) scan was calculated. The results indicated that DF described common duct more accurately compare to DHy. While, DHy described more accurately true and false ducts. Each time when connection of true and false duct appeared, true duct had lower brightness compare to common duct and false duct. Moreover, false duct characterized with higher brightness compare to common duct. In summary, the proposed algorithm mimics changes in brightness value for patients with acute type IIIb aortic dissection.
Collapse
|
9
|
Computational modeling of the fluid flow in type B aortic dissection using a modified finite element embedded formulation. Biomech Model Mechanobiol 2020; 19:1565-1583. [PMID: 31974816 DOI: 10.1007/s10237-020-01291-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2019] [Accepted: 01/14/2020] [Indexed: 10/25/2022]
Abstract
This work explores the use of an embedded computational fluid dynamics method to study the type B aortic dissection. The use of the proposed technique makes it possible to easily test different intimal flap configurations without any need of remeshing. To validate the presented methodology, we take as reference test case an in vitro experiment present in the literature. This experiment, which considers several intimal flap tear configurations (number, size and location), mimics the blood flow in a real type B aortic dissection. We prove the correctness and suitability of the presented approach by comparing the pressure values and waveform. The obtained results exhibit a remarkable similarity with the experimental reference data. Complementary, we present a feasible surgical application of the presented computer method. The aim is to help the clinicians in the decision making before the type B aortic dissection surgical fenestration. The capabilities of the proposed technique are exploited to efficiently create artificial reentry tear configurations. We highlight that only the radius and center of the reentry tear need to be specified by the clinicians, without any need to modify neither the model geometry nor the mesh. The obtained computational surgical fenestration results are in line with the medical observations in similar clinical studies.
Collapse
|