Obstructions in the lower airways lead to altered airflow patterns in the central airway

Structural and functional stenosis of the upper airway in Crouzon syndrome patients: A computational fluid dynamics analysis

OBJECTIVES

This study aimed to simulate the aerodynamics and to identify the spatial correlation between anatomical and functional stenoses in Crouzon syndrome patients.

METHODS

Six patients of Crouzon syndrome were included. Computational fluid dynamics (CFD) was utilized to simulate airflow dynamics, and characteristics, including the velocity, pressure intensity, wall shear stress, airflow resistance and streamline, were extracted for quantitative analysis both in overall and regionally. Structural stenosis was defined at the minimum cross-sectional area, while functional stenosis was identified at the point of maximum airflow velocity. The spatial distances between the Frankfurt plane and structural/functional stenosis were calculated and compared.

RESULTS

Structural stenosis occurred in the palatopharynx, while the highest inspiratory resistance and peak airflow velocity during expiration identified the glossopharynx as the functional stenosis site. A steep increase in negative pressure and a significant increase in wall shear stress could be observed surrounding the functional stenosis. The intensity and diffusion range of wall shear stress are positively correlated with age. Notably, the functional stenosis was consistently 5 mm below the structural stenosis (P < 0.05).

CONCLUSIONS

CFD effectively visualized both overall and regional aerodynamics of Crouzon syndrome, providing a novel method for functional airway evaluation. The spatial distributions of structural and functional stenoses did not strictly correspond; the structural stenosis was located on the palatopharynx, while the functional stenosis was on the glossopharynx. The wall shear stress worsens pathologically with age, aggravating functional stenosis to structural stenosis. Therefore, functional stenosis should also be addressed in airway management to ensure therapeutic effectiveness.

Modeling Realistic Geometries in Human Intrathoracic Airways

Geometrical models of the airways offer a comprehensive perspective on the complex interplay between lung structure and function. Originating from mathematical frameworks, these models have evolved to include detailed lung imagery, a crucial enhancement that aids in the early detection of morphological changes in the airways, which are often the first indicators of diseases. The accurate representation of airway geometry is crucial in research areas such as biomechanical modeling, acoustics, and particle deposition prediction. This review chronicles the evolution of these models, from their inception in the 1960s based on ideal mathematical constructs, to the introduction of advanced imaging techniques like computerized tomography (CT) and, to a lesser degree, magnetic resonance imaging (MRI). The advent of these techniques, coupled with the surge in data processing capabilities, has revolutionized the anatomical modeling of the bronchial tree. The limitations and challenges in both mathematical and image-based modeling are discussed, along with their applications. The foundation of image-based modeling is discussed, and recent segmentation strategies from CT and MRI scans and their clinical implications are also examined. By providing a chronological review of these models, this work offers insights into the evolution and potential future of airway geometry modeling, setting the stage for advancements in diagnosing and treating lung diseases. This review offers a novel perspective by highlighting how advancements in imaging techniques and data processing capabilities have significantly enhanced the accuracy and applicability of airway geometry models in both clinical and research settings. These advancements provide unique opportunities for developing patient-specific models.

Aerodynamic Simulation of Small Airway Resistance: A New Imaging Biomarker for Chronic Obstructive Pulmonary Disease

Purpose

To develop a novel method for calculating small airway resistance using computational fluid dynamics (CFD) based on CT data and evaluate its value to identify COPD.

Patients and Methods

24 subjects who underwent chest CT scans and pulmonary function tests between August 2020 and December 2020 were enrolled retrospectively. Subjects were divided into three groups: normal (10), high-risk (6), and COPD (8). The airway from the trachea down to the sixth generation of bronchioles was reconstructed by a 3D slicer. The small airway resistance (RSA) and RSA as a percentage of total airway resistance (RSA%) were calculated by CFD combined with airway resistance and FEV1 measured by pulmonary function test. A correlation analysis was conducted between RSA and pulmonary function parameters, including FEV1/FVC, FEV1% predicted, MEF50% predicted, MEF75% predicted and MMEF75/25% predicted.

Results

The RSA and RSA% were significantly different among the three groups (p<0.05) and related to FEV1/FVC (r = -0.70, p < 0.001; r = -0.67, p < 0.001), FEV1% predicted (r = -0.60, p = 0.002; r = -0.57, p = 0.004), MEF50% predicted (r = -0.64, p = 0.001; r = -0.64, p = 0.001), MEF75% predicted (r = -0.71, p < 0.001; r = -0.60, p = 0.002) and MMEF 75/25% predicted (r = -0.64, p = 0.001; r = -0.64, p = 0.001).

Conclusion

Airway CFD is a valuable method for estimating the small airway resistance, where the derived RSA will aid in the early diagnosis of COPD.

Targeted delivery of inhalable drug particles in the tracheobronchial tree model of a pediatric patient with bronchopneumonia: A numerical study

Pneumonia is a common cause of hospitalization and death in children worldwide. Inhalation therapy is one of the methods treating pneumonia However, there are limited studies that distinguish between the physiology of children and adults, especially with respect to targeted drug delivery. A tracheobronchial (TB) tree model of an 11-year-old child with bronchopneumonia is selected as a testbed for in silico trials of targeted drug delivery. The airflow and particle transport are solved by the computational fluid dynamics method at an airflow rate of 15 LPM. The results indicate that the distribution of deposited particles shows aggregation on the particle release map. Point-source aerosol release (PSAR) method can significantly reduce the deposition efficiency (DE) of particles in the TB tree model. Specifically, the PSAR method can reduce the DE of large particles (i.e., 7.5 µm and 10 µm) by 7.57% and 9.61%, respectively. This enables rapid design of patient-specific treatment for different population age groups and different airway diseases.

Evaluation of computational fluid dynamics models for predicting pediatric upper airway airflow characteristics

Computational fluid dynamics (CFD) has the potential for use as a clinical tool to predict the aerodynamics and respiratory function in the upper airway (UA) of children; however, careful selection of validated computational models is necessary. This study constructed a 3D model of the pediatric UA based on cone beam computed tomography (CBCT) imaging. The pediatric UA was 3D printed for pressure and velocity experiments, which were used as reference standards to validate the CFD simulation models. Static wall pressure and velocity distribution inside of the UA under inhale airflow rates from 0 to 266.67 mL/s were studied by CFD simulations based on the large eddy simulation (LES) model and four Reynolds-averaged Navier-Stokes (RANS) models. Our results showed that the LES performed best for pressure prediction; however, it was much more time-consuming than the four RANS models. Among the RANS models, the Low Reynolds number (LRN) SST k-ω model had the best overall performance at a series of airflow rates. Central flow velocity determined by particle image velocimetry was 3.617 m/s, while velocities predicted by the LES, LRN SST k-ω, and k-ω models were 3.681, 3.532, and 3.439 m/s, respectively. All models predicted jet flow in the oropharynx. These results suggest that the above CFD models have acceptable accuracy for predicting pediatric UA aerodynamics and that the LRN SST k-ω model has the most potential for clinical application in pediatric respiratory studies.

Evaluation of Airflow Sensitivity to the Truncation Level of a Realistic Human Airway Model in an Accurate Numerical Simulation

BACKGROUND

The truncation level of human airways is an influential factor in the analysis of respiratory flow in numerical simulations. Due to computational limitations and limited resolution of diagnostic medical imaging equipment, a truncated geometry of airways is always investigated.

OBJECTIVE

This study aimed to employ image-based geometries with zero generation and 5^th-generation truncation levels and assess bronchial airways truncation's effect on tracheal airflow characteristics.

MATERIAL AND METHODS

In this numerical study, computational fluid dynamics was employed to solve the respiratory flow in a realistic human airway model using the large eddy simulation technique coupling with the wall-adapting local eddy-viscosity (WALE) sub-grid scale model. The accuracy of numerical simulations was ensured by examining the large eddy simulation index of quality and Kolmogorov's K^-5/3 law.

RESULTS

The turbulent kinetic energy along the trachea has increased abnormally in the geometry with the zero-generation truncation level, and more severe fluctuations occurred in the velocity field of this geometry, which increased the tendency of each point to rotate. Compared to the extended model, the airflow's more chaotic behavior prevented larger-scale vortices from forming in the geometry with the zero-generation truncation level. Larger-scale vortices in the extended model caused the primary flow passing next to the vortices to accelerate more intensely, increasing the wall shear stress peaks in this geometry.

CONCLUSION

Eliminating the bronchial airways caused changes in tracheal airflow characteristics.

Structural and functional alterations of the tracheobronchial tree after left upper pulmonary lobectomy for lung cancer

Background

Pulmonary lobectomy has been a well-established curative treatment method for localized lung cancer. After left upper pulmonary lobectomy, the upward displacement of remaining lower lobe causes the distortion or kink of bronchus, which is associated with intractable cough and breathless. However, the quantitative study on structural and functional alterations of the tracheobronchial tree after lobectomy has not been reported. We sought to investigate these alterations using CT imaging analysis and computational fluid dynamics (CFD) method.

Methods

Both preoperative and postoperative CT images of 18 patients who underwent left upper pulmonary lobectomy are collected. After the tracheobronchial tree models are extracted, the angles between trachea and bronchi, the surface area and volume of the tree, and the cross-sectional area of left lower lobar bronchus are investigated. CFD method is further used to describe the airflow characteristics by the wall pressure, airflow velocity, lobar flow rate, etc.

Results

It is found that the angle between the trachea and the right main bronchus increases after operation, but the angle with the left main bronchus decreases. No significant alteration is observed for the surface area or volume of the tree between pre-operation and post-operation. After left upper pulmonary lobectomy, the cross-sectional area of left lower lobar bronchus is reduced for most of the patients (15/18) by 15–75%, especially for 4 patients by more than 50%. The wall pressure, airflow velocity and pressure drop significantly increase after the operation. The flow rate to the right lung increases significantly by 2–30% (but there is no significant difference between each lobe), and the flow rate to the left lung drops accordingly. Many vortices are found in various places with severe distortions.

Conclusions

The favorable and unfavorable adaptive alterations of tracheobronchial tree will occur after left upper pulmonary lobectomy, and these alterations can be clarified through CT imaging and CFD analysis. The severe distortions at left lower lobar bronchus might exacerbate postoperative shortness of breath.