Modeling the effect of tumor compression on airflow dynamics in trachea using contact simulation and CFD analysis

Stenotic geometry effects on airflow dynamics and respiration for central airway obstruction

BACKGROUND AND OBJECTIVE

The quantitative relationship between tracheal anatomy and ventilation function can be analyzed by using engineering-derived methods, including mathematical modeling and numerical simulations. In order to provide quantitative functional evaluation for patients with tracheobronchial stenosis, we here propose an aerodynamics-based assessment method by applying computational fluid dynamics analysis on synthetic and patient-specific airway models.

METHODS

By using 3D reconstruction of tracheobronchial tree and computational fluid dynamics simulations, the aerodynamic environment from the stenotic central airway down to the 4th-6th bifurcation of the tracheobronchial tree is examined in both synthetic and patient-derived models. The effects of stenotic anatomy (the degree of stenosis, stenotic length and location) on the aerodynamic parameters, including pressure drop, area-average velocity, volume flow rate, wall shear stress and airflow resistance, are investigated on three-dimensional models of tracheobronchial tree.

RESULTS

The results from 36 synthetic models demonstrate that 70% constriction marks the onset of a precipitous decrease in airflow relative to a normal airway. The analyses of simulation results of 8 patient-specific models indicate that the Myer-Cotton stenosis grading system can be interpreted in terms of aerodynamics-derived description, such as flow resistance. The tracheal stenosis significantly influences the resistance of peripheral bronchi, especially for patients with severe stenosis.

CONCLUSIONS

The present study forms a systematic framework for future development of more robust, bioengineering-informed evaluation methods for quantitative assessment of respiratory function of patients with central airway obstruction.

The impact of endotracheal intubation on oxygen delivery, trachea pressure and wall deformation

This paper concerns improving endotracheal tube (ETT) insertion through advanced computational science modelling. The study aims to better understand endotracheal intubation (ETI) and reduce medical errors in intensive and critical care units since ETT insertion is unique for each patient, depending on age, gender, size, physiology, and underlying health conditions. We have employed computational fluid dynamics and biomechanics modelling to investigate the effect of ETT for three ventilation modes on (a) local oxygen delivery to the lungs, (b) air pressure and wall shear stress at the tracheal walls, and (c) oscillatory elastic deformation of the tracheal tissues and muscle. For the first time, we reveal how the ventilation mode and ETT insertion in the trachea may induce major complications, especially in long periods of ETT. We show that rotating the ETT or displacing it by 2 mm only can induce a significant rise in the tracheal pressure up to 177 cmH2O. This study, for the first time, shows the vital role of computers in biology and medicine to provide enhanced decision-making-support to clinicians and medical doctors dealing with ETI.

Effect of cartilaginous rings in tracheal flow with stenosis

BACKGROUND

In respiratory fluid dynamics research, it is typically assumed that the wall of the trachea is smooth. However, the trachea is structurally supported by a series of cartilaginous rings that create undulations on the wall surface, which introduce perturbations into the flow. Even though many studies use realistic Computer Tomography (CT) scan data to capture the complex geometry of the respiratory system, its limited spatial resolution does not resolve small features, including those introduced by the cartilaginous rings.

RESULTS

Here we present an experimental comparison of two simplified trachea models with Grade II stenosis (70% blockage), one with smooth walls and second with cartilaginous rings. The use a unique refractive index-matching method provides unprecedented optical access and allowed us to perform non-intrusive velocity field measurements close to the wall (e.g., Particle Image Velocimetry (PIV)). Measurements were performed in a flow regime comparable to a resting breathing state (Reynolds number Re_D = 3350). The cartilaginous rings induce velocity fluctuations in the downstream flow, enhancing the near-wall transport of momentum flux and thus reducing flow separation in the downstream flow. The maximum upstream velocity in the recirculation region is reduced by 38%, resulting in a much weaker recirculation zone- a direct consequence of the cartilaginous rings.

CONCLUSIONS

These results highlight the importance of the cartilaginous rings in respiratory flow studies and the mechanism to reduce flow separation in trachea stenosis.

LES study on the impact of airway deformation on the airflow structures in the idealized mouth–throat model

To investigate the impacts of upper airway deformation on the airflow structures, the airflow fields in the trachea are simulated using three geometrical models considering three different levels of airway deformations. Structured grids are used to create the high-quality grids. Large eddy simulation with the Smagorinsky sub-grid model is adopted to solve the three-dimensional in-compressible Navier–Stokes equations using the solver pisoFoam in the open-source CFD software OpenFOAM. The numerical results demonstrate that the airway deformation influences the main airflow structures depending on the deformation level. Particularly, it slightly impacts on the laryngeal jet such as the profile and the strength of laryngeal jet. The strength of the laryngeal jet increases slightly for the heavy deformation. In contrast, it impacts on the recirculation zone, secondary vortices, and turbulent kinetic energy more obviously. The increasing airway deformation will produce stronger secondary flow, smaller recirculation zone, and weaker turbulent kinetic energy. The turbulence intensity distribution varies as well. The obviously impacted flow region is mainly within the region of one to six tracheal diameters downstream the glottis.