A 4DCT imaging-based breathing lung model with relative hysteresis

• Published in: Journal of computational physics Vol. 326; pp. 76 - 90
• Main Authors: Miyawaki, Shinjiro, Choi, Sanghun, Hoffman, Eric A., Lin, Ching-Long
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 01.12.2016
• Subjects: 4DCT; ALGORITHMS; BIOMEDICAL RADIOGRAPHY; BOUNDARY CONDITIONS; CAT SCANNING; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; DEFORMATION; DISTRIBUTION; FLUID MECHANICS; FLUIDS; FOUR-DIMENSIONAL CALCULATIONS; HYSTERESIS; Image registration; INTERPOLATION; LUNGS; MATHEMATICAL METHODS AND COMPUTING; Multiscale; PRESSURE DROP; Pulmonary airflow; Regional ventilation; Relative hysteresis; RESPIRATION; SURFACES; Image registration; 4DCT; Pulmonary airflow; Relative hysteresis; Multiscale; Regional ventilation; pulmonary airflow; regional ventilation; relative hysteresis; multiscale; image registration
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2016.08.039

To reproduce realistic airway motion and airflow, the authors developed a deforming lung computational fluid dynamics (CFD) model based on four-dimensional (4D, space and time) dynamic computed tomography (CT) images. A total of 13 time points within controlled tidal volume respiration were used to account for realistic and irregular lung motion in human volunteers. Because of the irregular motion of 4DCT-based airways, we identified an optimal interpolation method for airway surface deformation during respiration, and implemented a computational solid mechanics-based moving mesh algorithm to produce smooth deforming airway mesh. In addition, we developed physiologically realistic airflow boundary conditions for both models based on multiple images and a single image. Furthermore, we examined simplified models based on one or two dynamic or static images. By comparing these simplified models with the model based on 13 dynamic images, we investigated the effects of relative hysteresis of lung structure with respect to lung volume, lung deformation, and imaging methods, i.e., dynamic vs. static scans, on CFD-predicted pressure drop. The effect of imaging method on pressure drop was 24 percentage points due to the differences in airflow distribution and airway geometry. •We developed a breathing human lung CFD model based on 4D-dynamic CT images.•The 4DCT-based breathing lung model is able to capture lung relative hysteresis.•A new boundary condition for lung model based on one static CT image was proposed.•The difference between lung models based on 4D and static CT images was quantified.