1
|
Jiang F, Hirano T, Liang C, Zhang G, Matsunaga K, Chen X. Multi-scale simulations of pulmonary airflow based on a coupled 3D-1D-0D model. Comput Biol Med 2024; 171:108150. [PMID: 38367450 DOI: 10.1016/j.compbiomed.2024.108150] [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: 08/11/2023] [Revised: 12/25/2023] [Accepted: 02/12/2024] [Indexed: 02/19/2024]
Abstract
Pulmonary airflow simulation is a valuable tool for studying respiratory function and disease. However, the respiratory system is a complex multiscale system that involves various physical and biological processes across different spatial and temporal scales. In this study, we propose a 3D-1D-0D multiscale method for simulating pulmonary airflow, which integrates different levels of detail and complexity of the respiratory system. The method consists of three components: a 3D computational fluid dynamics model for the airflow in the trachea and bronchus, a 1D pipe model for the airflow in the terminal bronchioles, and a 0D biphasic mixture model for the airflow in the respiratory bronchioles and alveoli coupled with the lung deformation. The coupling between the different components is achieved by satisfying the mass and momentum conservation law and the pressure continuity condition at the interfaces. We demonstrate the validity and applicability of our method by comparing the results with data of previous models. We also investigate the reduction in inhaled air volume due to the pulmonary fibrosis using the developed multiscale model. Our method provides a comprehensive and realistic framework for simulating pulmonary airflow and can potentially facilitate the diagnosis and treatment of respiratory diseases.
Collapse
Affiliation(s)
- Fei Jiang
- Department of Mechanical Engineering, Graduate School of Sciences and Technology for Innovation, Yamaguchi University, Tokiwadai, Ube, 7558611, Yamaguchi, Japan; Biomedical Engineering Center (YUBEC), Tokiwadai, Ube, 7558611, Yamaguchi, Japan.
| | - Tsunahiko Hirano
- Department of Respiratory Medicine and Infectious Disease, Graduate School of Medicine, Yamaguchi University, Minamikogushi, Ube, 7558505, Yamaguchi, Japan
| | - Chenyang Liang
- Department of Mechanical Engineering, Graduate School of Sciences and Technology for Innovation, Yamaguchi University, Tokiwadai, Ube, 7558611, Yamaguchi, Japan
| | - Guangzhi Zhang
- Keisoku Engineering System Co., Ltd., Uchikanda, Chiyoda-ku, Tokyo, 1010047, Japan
| | - Kazuto Matsunaga
- Department of Respiratory Medicine and Infectious Disease, Graduate School of Medicine, Yamaguchi University, Minamikogushi, Ube, 7558505, Yamaguchi, Japan
| | - Xian Chen
- Department of Mechanical Engineering, Graduate School of Sciences and Technology for Innovation, Yamaguchi University, Tokiwadai, Ube, 7558611, Yamaguchi, Japan; Biomedical Engineering Center (YUBEC), Tokiwadai, Ube, 7558611, Yamaguchi, Japan
| |
Collapse
|
2
|
Goodarzi Ardakani V, Gambaruto AM, Silva G, Pereira R. A porosity model for medical image segmentation of vessels. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2022; 38:e3580. [PMID: 35142065 PMCID: PMC9285562 DOI: 10.1002/cnm.3580] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Revised: 02/03/2022] [Accepted: 02/05/2022] [Indexed: 06/14/2023]
Abstract
A physics-based medical image segmentation method is developed. Specifically, the image greyscale intensity is used to infer the voxel partial volumes and subsequently formulate a porous medium analogy. The method involves first translating the medical image volumetric data into a three-dimensional computational domain of a porous material. A velocity field is then obtained from numerical simulations of incompressible fluid flow in the porous material, and finally a velocity iso-surface provides the surface description of the target object. The approach is tested on CT images of eight patient-specific cases, where cerebral aneurysms, nasal cavities (NC), and an aortic arch (AA) are the objects of interest. In the aneurysm cases, the results are compared against constant greyscale thresholding and manual segmentation. The manual segmentations of the aneurysms are validated by a clinical practitioner. Only a qualitative comparison is available for the NC, and the AA geometries. The results show that the proposed method is effective and capable of extracting the target object in a noisy domain. A sensitivity study is carried out to verify the method's performance with respect to modelling or user choices. The segmentation by the proposed method is also evaluated by performing computational fluid dynamics simulation, including a near-wall flow analysis, to ensure that the segmented geometry and the resulting computed solution are representative and meaningful.
Collapse
Affiliation(s)
| | | | - Goncalo Silva
- IDMEC, Mechatronics DepartmentUniversity of ÉvoraÉvoraPortugal
| | - Ricardo Pereira
- Department of NeurosurgeryCoimbra University Hospital Center, Faculty of Medicine, University of CoimbraCoimbraPortugal
| |
Collapse
|
3
|
Volumetric lattice Boltzmann method for wall stresses of image-based pulsatile flows. Sci Rep 2022; 12:1697. [PMID: 35105911 PMCID: PMC8807599 DOI: 10.1038/s41598-022-05269-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Accepted: 12/23/2021] [Indexed: 11/18/2022] Open
Abstract
Image-based computational fluid dynamics (CFD) has become a new capability for determining wall stresses of pulsatile flows. However, a computational platform that directly connects image information to pulsatile wall stresses is lacking. Prevailing methods rely on manual crafting of a hodgepodge of multidisciplinary software packages, which is usually laborious and error-prone. We present a new computational platform, to compute wall stresses in image-based pulsatile flows using the volumetric lattice Boltzmann method (VLBM). The novelty includes: (1) a unique image processing to extract flow domain and local wall normality, (2) a seamless connection between image extraction and VLBM, (3) an en-route calculation of strain-rate tensor, and (4) GPU acceleration (not included here). We first generalize the streaming operation in the VLBM and then conduct application studies to demonstrate its reliability and applicability. A benchmark study is for laminar and turbulent pulsatile flows in an image-based pipe (Reynolds number: 10 to 5000). The computed pulsatile velocity and shear stress are in good agreements with Womersley's analytical solutions for laminar pulsatile flows and concurrent laboratory measurements for turbulent pulsatile flows. An application study is to quantify the pulsatile hemodynamics in image-based human vertebral and carotid arteries including velocity vector, pressure, and wall-shear stress. The computed velocity vector fields are in reasonably well agreement with MRA (magnetic resonance angiography) measured ones. This computational platform is good for image-based CFD with medical applications and pore-scale porous media flows in various natural and engineering systems.
Collapse
|