1
|
Decroocq M, Frindel C, Rougé P, Ohta M, Lavoué G. Modeling and hexahedral meshing of cerebral arterial networks from centerlines. Med Image Anal 2023; 89:102912. [PMID: 37549612 DOI: 10.1016/j.media.2023.102912] [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: 12/31/2021] [Revised: 06/12/2023] [Accepted: 07/25/2023] [Indexed: 08/09/2023]
Abstract
Computational fluid dynamics (CFD) simulation provides valuable information on blood flow from the vascular geometry. However, it requires extracting precise models of arteries from low-resolution medical images, which remains challenging. Centerline-based representation is widely used to model large vascular networks with small vessels, as it encodes both the geometric and topological information and facilitates manual editing. In this work, we propose an automatic method to generate a structured hexahedral mesh suitable for CFD directly from centerlines. We addressed both the modeling and meshing tasks. We proposed a vessel model based on penalized splines to overcome the limitations inherent to the centerline representation, such as noise and sparsity. The bifurcations are reconstructed using a parametric model based on the anatomy that we extended to planar n-furcations. Finally, we developed a method to produce a volume mesh with structured, hexahedral, and flow-oriented cells from the proposed vascular network model. The proposed method offers better robustness to the common defects of centerlines and increases the mesh quality compared to state-of-the-art methods. As it relies on centerlines alone, it can be applied to edit the vascular model effortlessly to study the impact of vascular geometry and topology on hemodynamics. We demonstrate the efficiency of our method by entirely meshing a dataset of 60 cerebral vascular networks. 92% of the vessels and 83% of the bifurcations were meshed without defects needing manual intervention, despite the challenging aspect of the input data. The source code is released publicly.
Collapse
Affiliation(s)
- Méghane Decroocq
- CREATIS, Université Lyon1, CNRS UMR5220, INSERM U1206, INSA-Lyon, 69621 Villeurbanne, France; LIRIS, CNRS UMR 5205, F-69621, France; ELyTMaX IRL3757, CNRS, INSA Lyon, Centrale Lyon, Université Claude Bernard Lyon 1, Tohoku University, 980-8577, Sendai, Japan; Institute of Fluid Science, Tohoku University, 2-1-1, Katahira, Aoba-ku, Sendai, Miyagi 980-8577, Japan; Graduate School of Biomedical Engineering, Tohoku University, 6-6 Aramaki-aza-aoba, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Carole Frindel
- CREATIS, Université Lyon1, CNRS UMR5220, INSERM U1206, INSA-Lyon, 69621 Villeurbanne, France; Institute of Fluid Science, Tohoku University, 2-1-1, Katahira, Aoba-ku, Sendai, Miyagi 980-8577, Japan.
| | - Pierre Rougé
- ELyTMaX IRL3757, CNRS, INSA Lyon, Centrale Lyon, Université Claude Bernard Lyon 1, Tohoku University, 980-8577, Sendai, Japan; Université de Reims Champagne Ardenne, CReSTIC, 51100 Reims, France
| | - Makoto Ohta
- ELyTMaX IRL3757, CNRS, INSA Lyon, Centrale Lyon, Université Claude Bernard Lyon 1, Tohoku University, 980-8577, Sendai, Japan; Institute of Fluid Science, Tohoku University, 2-1-1, Katahira, Aoba-ku, Sendai, Miyagi 980-8577, Japan
| | - Guillaume Lavoué
- LIRIS, CNRS UMR 5205, F-69621, France; Ecole Centrale de Lyon, France
| |
Collapse
|
2
|
Ho H, Means S, Safaei S, Hunter PJ. In silico modeling for the hepatic circulation and transport: From the liver organ to lobules. WIREs Mech Dis 2023; 15:e1586. [PMID: 36131627 DOI: 10.1002/wsbm.1586] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2021] [Revised: 08/08/2022] [Accepted: 08/15/2022] [Indexed: 11/12/2022]
Abstract
The function of the liver depends critically on its blood supply. Numerous in silico models have been developed to study various aspects of the hepatic circulation, including not only the macro-hemodynamics at the organ level, but also the microcirculation at the lobular level. In addition, computational models of blood flow and bile flow have been used to study the transport, metabolism, and clearance of drugs in pharmacokinetic studies. These in silico models aim to provide insights into the liver organ function under both healthy and diseased states, and to assist quantitative analysis for surgical planning and postsurgery treatment. The purpose of this review is to provide an update on state-of-the-art in silico models of the hepatic circulation and transport processes. We introduce the numerical methods and the physiological background of these models. We also discuss multiscale frameworks that have been proposed for the liver, and their linkage with the large context of systems biology, systems pharmacology, and the Physiome project. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering Cardiovascular Diseases > Computational Models.
Collapse
Affiliation(s)
- Harvey Ho
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Shawn Means
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Soroush Safaei
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Peter John Hunter
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
3
|
T A, Narayan R, Shenoy PA, Nayak UY. Computational modeling for the design and development of nano based drug delivery systems. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
4
|
Albadawi M, Abuouf Y, Elsagheer S, Sekiguchi H, Ookawara S, Ahmed M. Influence of Rigid-Elastic Artery Wall of Carotid and Coronary Stenosis on Hemodynamics. Bioengineering (Basel) 2022; 9:708. [PMID: 36421109 PMCID: PMC9687628 DOI: 10.3390/bioengineering9110708] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Revised: 11/15/2022] [Accepted: 11/16/2022] [Indexed: 12/18/2023] Open
Abstract
Cardiovascular system abnormalities can result in serious health complications. By using the fluid-structure interaction (FSI) procedure, a comprehensive realistic approach can be employed to accurately investigate blood flow coupled with arterial wall response. The hemodynamics was investigated in both the coronary and carotid arteries based on the arterial wall response. The hemodynamics was estimated based on the numerical simulation of a comprehensive three-dimensional non-Newtonian blood flow model in elastic and rigid arteries. For stenotic right coronary artery (RCA), it was found that the maximum value of wall shear stress (WSS) for the FSI case is higher than the rigid wall. On the other hand, for the stenotic carotid artery (CA), it was found that the maximum value of WSS for the FSI case is lower than the rigid wall. Moreover, at the peak systole of the cardiac cycle (0.38 s), the maximum percentage of arterial wall deformation was found to be 1.9%. On the other hand, for the stenotic carotid artery, the maximum percentage of arterial wall deformation was found to be 0.46%. A comparison between FSI results and those obtained by rigid wall arteries is carried out. Findings indicate slight differences in results for large-diameter arteries such as the carotid artery. Accordingly, the rigid wall assumption is plausible in flow modeling for relatively large diameters such as the carotid artery. Additionally, the FSI approach is essential in flow modeling in small diameters.
Collapse
Affiliation(s)
- Muhamed Albadawi
- Department of Energy Resources Engineering, Egypt-Japan University of Science and Technology (E-JUST), P.O. Box 179, New Borg El-Arab City 5221241, Egypt
- Biomedical Flow Dynamics Laboratory, Institute of Fluid Science, Tohoku University, Sendai 980-8577, Japan
- Engineering Mathematics and Physics Department, Faculty of Engineering, Alexandria University, Alexandria 5424041, Egypt
| | - Yasser Abuouf
- Mechanical Engineering Department, Faculty of Engineering, Alexandria University, Alexandria 5424041, Egypt
| | - Samir Elsagheer
- Department of Energy Resources Engineering, Egypt-Japan University of Science and Technology (E-JUST), P.O. Box 179, New Borg El-Arab City 5221241, Egypt
- Faculty of Engineering, Aswan University, Aswan 81528, Egypt
| | - Hidetoshi Sekiguchi
- Department of Energy Resources Engineering, Egypt-Japan University of Science and Technology (E-JUST), P.O. Box 179, New Borg El-Arab City 5221241, Egypt
- Department of Chemical Engineering, Graduate School of Science and Engineering, Tokyo Institute of Technology, Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Shinichi Ookawara
- Department of Energy Resources Engineering, Egypt-Japan University of Science and Technology (E-JUST), P.O. Box 179, New Borg El-Arab City 5221241, Egypt
- Department of Chemical Engineering, Graduate School of Science and Engineering, Tokyo Institute of Technology, Ookayama, Meguro-ku, Tokyo 152-8552, Japan
| | - Mahmoud Ahmed
- Department of Energy Resources Engineering, Egypt-Japan University of Science and Technology (E-JUST), P.O. Box 179, New Borg El-Arab City 5221241, Egypt
- Mechanical Engineering Department, Assiut University, Assiut 71516, Egypt
| |
Collapse
|
5
|
Taebi A, Janibek N, Goldman R, Pillai R, Vu CT, Roncali E. The Impact of Injection Distance to Bifurcations on Yttrium-90 Distribution in Liver Cancer Radioembolization. J Vasc Interv Radiol 2022; 33:668-677.e1. [PMID: 35301128 PMCID: PMC9156550 DOI: 10.1016/j.jvir.2022.03.006] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2021] [Revised: 02/22/2022] [Accepted: 03/06/2022] [Indexed: 01/26/2023] Open
Abstract
PURPOSE To model the effect of the injection location on the distribution of yttrium-90 (90Y) microspheres in the liver during radioembolization using computational simulation and to determine the potential effects of radial movements of the catheter tip. MATERIALS AND METHODS Numerical studies were conducted using images from a representative patient with hepatocellular carcinoma. The right hepatic artery (RHA) was segmented from contrast-enhanced cone-beam computed tomography scans. The blood flow was investigated in the trunk of the RHA using numerical simulations for 6 injection position scenarios at 2 sites located at a distance of approximately 5 and 20 mm upstream of the first bifurcation (RHA diameters of approximately 4.6 mm). The 90Y delivery to downstream vessels was calculated from the simulated hepatic artery hemodynamics. RESULTS Varying the injection location along the RHA and across the vessel cross-section resulted in different simulated microsphere distributions in the downstream vascular bed. When the catheter tip was 5 mm upstream of the bifurcation, 90Y distribution in the downstream branches varied by as much as 53% with a 1.5-mm radial movement of the tip. However, the catheter radial movement had a weaker effect on the microsphere distribution when the injection plane was farther from the first bifurcation (20 mm), with a maximum delivery variation of 9% to a downstream branch. CONCLUSIONS An injection location far from bifurcations is recommended to minimize the effect of radial movements of the catheter tip on the microsphere distribution.
Collapse
Affiliation(s)
- Amirtaha Taebi
- Department of Agricultural and Biological Engineering, Mississippi State University
| | - Nursultan Janibek
- Department of Mechanical and Aerospace Engineering, University of California Davis
| | - Roger Goldman
- Department of Radiology, University of California Davis
| | - Rex Pillai
- Department of Radiology, University of California Davis
| | | | - Emilie Roncali
- Department of Radiology, University of California Davis,Department of Biomedical Engineering, University of California Davis
| |
Collapse
|
6
|
Mahmoudi M, Jennings C, Pereira K, Hall AF, Arzani A. Guiding the prostatic artery embolization procedure with computational uid dynamics. J Biomech Eng 2022; 144:1140781. [DOI: 10.1115/1.4054515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Indexed: 11/08/2022]
Abstract
Abstract
Benign prostatic hyperplasia (BPH) is a common disease associated with lower urinary tract symptoms and the most frequent benign tumor in men. To reduce BPH therapy complications, prostatic artery embolization (PAE) was developed to replace the surgical options. PAE is a minimally invasive technique in which emboli are injected into the prostate arteries (PA), obstructing the blood flow in the hypervascular nodules. In this work, a personalized PAE treatment strategy was proposed using patient-specific computational fluid dynamics (CFD). First, the hemodynamics environment in the iliac arterial tree considering a large network of bifurcations was studied. The results showed complex blood flow patterns in the iliac arterial network. Subsequently, the transport of embolic particulates during PAE for the standard horizontal and a hypothetical vertical patient positioning was simulated using Lagrangian particle tracking. Emboli with different sizes were released at various locations across the iliac arterial tree. The emboli entering the PA were mapped back to their initial location to create emboli release maps (ERMs). The obtained ERMs during the standard patient positioning for smaller emboli at certain release locations showed distinct regions in which if the emboli were released within these regions, all of them would reach the PA without non-target embolization. During the hypothetical vertical patient positioning, the larger emboli formed a larger coherent region in the ERMs. Our patient-specific model can be used to find the best spatial location for emboli injection and perform the embolization procedure with minimal off-target delivery.
Collapse
Affiliation(s)
- Mostafa Mahmoudi
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, United States
| | - Chadrick Jennings
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, United States
| | - Keith Pereira
- Department of Radiology, Saint Louis University, St Louis, MO, United States
| | - Andrew F. Hall
- Department of Biomedical Engineering, Saint Louis University, St Louis, MO, United States
| | - Amirhossein Arzani
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, United States
| |
Collapse
|
7
|
Aramburu J, Antón R, Rodríguez-Fraile M, Sangro B, Bilbao JI. Computational Fluid Dynamics Modeling of Liver Radioembolization: A Review. Cardiovasc Intervent Radiol 2021; 45:12-20. [PMID: 34518913 PMCID: PMC8716346 DOI: 10.1007/s00270-021-02956-5] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/29/2021] [Accepted: 08/25/2021] [Indexed: 12/16/2022]
Abstract
Yttrium-90 radioembolization (RE) is a widely used transcatheter intraarterial therapy for patients with unresectable liver cancer. In the last decade, computer simulations of hepatic artery hemodynamics during RE have been performed with the aim of better understanding and improving the therapy. In this review, we introduce the concept of computational fluid dynamics (CFD) modeling with a clinical perspective and we review the CFD models used to study RE from the fluid mechanics point of view. Finally, we show what CFD simulations have taught us about the hemodynamics during RE, the current capabilities of CFD simulations of RE, and we suggest some future perspectives.
Collapse
Affiliation(s)
- Jorge Aramburu
- Universidad de Navarra, TECNUN Escuela de Ingeniería, 20018, Donostia-San Sebastián, Spain.
| | - Raúl Antón
- Universidad de Navarra, TECNUN Escuela de Ingeniería, 20018, Donostia-San Sebastián, Spain.,IdiSNA, Instituto de Investigación Sanitaria de Navarra, 31008, Pamplona, Spain
| | - Macarena Rodríguez-Fraile
- IdiSNA, Instituto de Investigación Sanitaria de Navarra, 31008, Pamplona, Spain.,Department of Nuclear Medicine, Clínica Universidad de Navarra, 31008, Pamplona, Spain
| | - Bruno Sangro
- IdiSNA, Instituto de Investigación Sanitaria de Navarra, 31008, Pamplona, Spain.,Liver Unit, Clínica Universidad de Navarra and CIBEREHD, 31008, Pamplona, Spain
| | - José Ignacio Bilbao
- IdiSNA, Instituto de Investigación Sanitaria de Navarra, 31008, Pamplona, Spain.,Department of Radiology, Clínica Universidad de Navarra, 31008, Pamplona, Spain
| |
Collapse
|
8
|
Taebi A, Berk S, Roncali E. Realistic boundary conditions in SimVascular through inlet catheter modeling. BMC Res Notes 2021; 14:215. [PMID: 34103097 PMCID: PMC8186195 DOI: 10.1186/s13104-021-05631-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2021] [Accepted: 05/22/2021] [Indexed: 01/25/2023] Open
Abstract
OBJECTIVE This study aims at developing a pipeline that provides the capability to include the catheter effect in the computational fluid dynamics (CFD) simulations of the cardiovascular system and other human vascular flows carried out with the open-source software SimVascular. This tool is particularly useful for CFD simulation of interventional radiology procedures such as tumor embolization where estimation of a therapeutic agent distribution is of interest. RESULTS A pipeline is developed that generates boundary condition files which can be used in SimVascular CFD simulations. The boundary condition files are modified such that they simulate the effect of catheter presence on the flow field downstream of the inlet. Using this pipeline, the catheter flow, velocity profile, radius, wall thickness, and deviation from the vessel center can be defined. Since our method relies on the manipulation of the boundary condition that is imposed on the inlet, it is sensitive to the mesh density. The finer the mesh is (especially around the catheter wall), the more accurate the velocity estimations are. In this study, we also utilized this pipeline to qualitatively investigate the effect of catheter presence on the flow field in a truncated right hepatic arterial tree of a liver cancer patient.
Collapse
Affiliation(s)
- Amirtahà Taebi
- Department of Biomedical Engineering, University of California, Davis, One Shields Ave, Davis, CA, 95616-5270, USA.
| | - Selin Berk
- Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, 420 Westwood Plaza, Los Angeles, CA, 90095-1597, USA
| | - Emilie Roncali
- Department of Biomedical Engineering, University of California, Davis, One Shields Ave, Davis, CA, 95616-5270, USA.,Department of Radiology, University of California, Davis, 4860 Y Street, Suite 3100, Sacramento, CA, 95817, USA
| |
Collapse
|