Assessing Changes in Airflow and Energy Loss in a Progressive Tracheal Compression Before and After Surgical Correction

Influence of catheter thickness on respiratory physiology during less invasive surfactant administration in extremely preterm infants

Introduction

Delivering surfactant via thin catheters (minimal-invasive surfactant therapy (MIST); less invasive surfactant administration (LISA)) has become a common procedure. However, the effect of tracheal obstruction caused by catheters of different sizes on tracheal resistance in extremely low gestational age newborns (ELGANs) is unknown.

Methods

To investigate the effect of catheters size 3.5, 5 and 6 French on airway resistance in ELGANs of 23-28 weeks gestational age during LISA, we performed calculations based on Hagen-Poiseuille's law and compared these with a clinically and physically more accurate method: computational fluid dynamics (CFD) simulations of respiratory airflow, performed in 3D virtual airway models derived from MRI.

Results

The presence of the above catheters decreased the cross-sectional area of the infants' tracheal entrance (the cricoid ring) by 13-53%. Hagen-Poiseuille's law predicted an increase in resistance by 1.5-4.5 times and 1.3-2.6 times in ELGANs born at 23 and 28 weeks, respectively. However, CFD simulations demonstrated an even higher increase in resistance of 3.4-85.1 and 1.1-3.5 times, respectively. The higher calculated resistances were due to the extremely narrow remaining lumen at the glottis and cricoid with the catheter inserted, resulting in a stronger glottal jet and turbulent airflow, which was not predicted by Hagen-Poiseuille.

Conclusion

Catheter thickness can greatly increase tracheal resistance during LISA-procedures in ELGANs. Based on these models, it is recommended to use the thinnest catheter possible during LISA in ELGANs to avoid unnecessary increases in airway resistance in infants already experiencing dyspnea due to respiratory distress syndrome.

The effect of including dynamic imaging derived airway wall motion in CFD simulations of respiratory airflow in patients with OSA

Obstructive sleep apnea (OSA) is an airway disease caused by periodic collapse of the airway during sleep. Imaging-based subject-specific computational fluid dynamics (CFD) simulations allow non-invasive assessment of clinically relevant metrics such as total pressure loss (TPL) in patients with OSA. However, most of such studies use static airway geometries, which neglect physiological airway motion. This study aims to quantify how much the airway moves during the respiratory cycle, and to determine how much this motion affects CFD pressure loss predictions. Motion of the airway wall was quantified using cine MRI data captured over a single respiratory cycle in three subjects with OSA. Synchronously-measured respiratory airflow was used as the flow boundary condition for all simulations. Simulations were performed for full respiratory cycles with 5 different wall boundary conditions: (1) a moving airway wall, and static airway walls at (2) peak inhalation, (3) end inhalation, (4) peak exhalation, and (5) end exhalation. Geometric analysis exposed significant local airway cross-sectional area (CSA) variability, with local CSA varying as much as 300%. The comparative CFD simulations revealed the discrepancies between dynamic and static wall simulations are subject-specific, with TPL differing by up to 400% between static and dynamic simulations. There is no consistent pattern to which static wall CFD simulations overestimate or underestimate the airway TPL. This variability underscores the complexity of accurately modeling airway physiology and the importance of considering dynamic anatomical factors to predict realistic respiratory airflow dynamics in patients with OSA.

Effects of decongestion on nasal cavity air conditioning efficiency: a CFD cohort study

Decongestion reduces blood flow in the nasal turbinates, enlarging the airway lumen. Although the enlarged airspace reduces the trans-nasal inspiratory pressure drop, symptoms of nasal obstruction may relate to nasal cavity air-conditioning. Thus, it is necessary to quantify the efficiency of nasal cavity conditioning of the inhaled air. This study quantifies both overall and regional nasal air-conditioning in a cohort of 10 healthy subjects using computational fluid dynamics simulations before and after nasal decongestion. The 3D virtual geometry model was segmented from magnetic resonance images (MRI). Each subject was under two MRI acquisitions before and after the decongestion condition. The effects of decongestion on nasal cavity air conditioning efficiency were modelled at two inspiratory flowrates: 15 and 30 L min-1 to represent restful and light exercise conditions. Results show inhaled air was both heated and humidified up to 90% of alveolar conditions at the posterior septum. The air-conditioning efficiency of the nasal cavity remained nearly constant between nostril and posterior septum but dropped significantly after posterior septum. In summary, nasal cavity decongestion not only reduces inhaled air added heat by 23% and added moisture content by 19%, but also reduces the air-conditioning efficiency by 35% on average.

Tracheomalacia Reduces Aerosolized Drug Delivery to the Lung

Rationale: Neonates with respiratory issues are frequently treated with aerosolized medications to manage lung disease or facilitate airway clearance. Dynamic tracheal collapse (tracheomalacia [TM]) is a common comorbidity in these patients, but it is unknown whether the presence of TM alters the delivery of aerosolized drugs. Objectives: To quantify the effect of neonatal TM on the delivery of aerosolized drugs. Methods: Fourteen infant subjects with respiratory abnormalities were recruited; seven with TM and seven without TM. Respiratory-gated 3D ultrashort echo time magnetic resonance imaging (MRI) was acquired covering the central airway and lungs. For each subject, a computational fluid dynamics simulation modeled the airflow and particle transport in the central airway based on patient-specific airway anatomy, motion, and airflow rates derived from MRI. Results: Less aerosolized drug reached the distal airways in subjects with TM than in subjects without TM: of the total drug delivered, less particle mass passed through the main bronchi in subjects with TM compared with subjects without TM (33% vs. 47%, p = 0.013). In subjects with TM, more inhaled particles were deposited on the surface of the airway (48% vs. 25%, p = 0.003). This effect becomes greater with larger particle sizes and is significant for particles with a diameter >2 μm (2-5 μm, p ≤ 0.025 and 5-15 μm, p = 0.004). Conclusions: Neonatal patients with TM receive less aerosolized drug delivered to the lungs than subjects without TM. Currently, infants with lung disease and TM may not be receiving adequate and/or expected medication. Particles >2 μm in diameter are likely to deposit on the surface of the airway due to anatomical constrictions such as reduced tracheal and glottal cross-sectional area in neonates with TM. This problem could be alleviated by delivering smaller aerosolized particles.

Analysis of Upper Airway Flow Dynamics in Robin Sequence Infants Using 4-D Computed Tomography and Computational Fluid Dynamics

Robin Sequence (RS) is a potentially fatal craniofacial condition characterized by undersized jaw, posteriorly displaced tongue, and resultant upper airway obstruction (UAO). Accurate assessment of UAO severity is crucial for management and diagnosis of RS, yet current evaluation modalities have significant limitations and no quantitative measures of airway resistance exist. In this study, we combine 4-dimensional computed tomography and computational fluid dynamics (CFD) to assess, for the first time, UAO severity using fluid dynamic metrics in RS patients. Dramatic intrapopulation differences are found, with the ratio between most and least severe patients in breathing resistance, energy loss, and peak velocity equal to 40:1, 20:1, and 6:1, respectively. Analysis of local airflow dynamics characterized patients as presenting with primary obstructions either at the location of the tongue base, or at the larynx, with tongue base obstructions resulting in a more energetic stenotic jet and greater breathing resistance. Finally, CFD-derived flow metrics are found to correlate with the level of clinical respiratory support. Our results highlight the large intrapopulation variability, both in quantitative metrics of UAO severity (resistance, energy loss, velocity) and in the location and intensity of stenotic jets for RS patients. These results suggest that computed airflow metrics may significantly improve our understanding of UAO and its management in RS.

Computational fluid dynamics assessment of congenital tracheal stenosis

PURPOSE

The severity of congenital tracheal stenosis (CTS) is commonly evaluated based on the degree of stenosis. However, it does not always reflect the clinical respiratory status. We applied computational fluid dynamics (CFD) to the assessment of CTS. The aim of this study was to evaluate its validity.

METHODS

CFD models were constructed on 15 patients (12 preoperative models and 15 postoperative models) with CTS before and after surgery, using the computed tomographic data. Energy flux, needed to drive airflow, measured by CFD and the minimum cross-sectional area of the trachea (MCAT) were quantified and evaluated retrospectively.

RESULTS

The energy flux correlated positively with the clinical respiratory status before and after surgery (rs = 0.611, p = 0.035 and rs = 0.591, p = 0.020, respectively). Although MCAT correlated negatively with the clinical respiratory status before surgery (rs = -0.578, p = 0.044), there was not significant correlation between the two after surgery (p = 0.572).

CONCLUSIONS

The energy flux measured by CFD assessment reflects the respiratory status in CTS before and after surgery. CFD can be an additional objective and quantitative evaluation tool for CTS.

Fluid dynamic assessment of positive end-expiratory pressure in a tracheostomy tube connector during respiration

High-flow oxygen therapy using a tracheostomy tube is a promising clinical approach to reduce the work of breathing in tracheostomized patients. Positive end-expiratory pressure (PEEP) is usually applied during oxygen inflow to improve oxygenation by preventing end-expiratory lung collapse. However, much is still unknown about the geometrical effects of PEEP, especially regarding tracheostomy tube connectors (or adapters). Quantifying the degree of end-expiratory pressure (EEP) that takes patient-specific spirometry into account would be useful in this regard, but no such framework has been established yet. Thus, a platform to assess PEEP under respiration was developed, wherein three-dimensional simulation of airflow in a tracheostomy tube connector is coupled with a lumped lung model. The numerical model successfully reflected the magnitude of EEP measured experimentally using a lung phantom. Numerical simulations were further performed to quantify the effects of geometrical parameters on PEEP, such as inlet angles and rate of stenosis in the connector. Although sharp inlet angles increased the magnitude of EEP, they cannot be expected to achieve clinically reasonable PEEP. On the other hand, geometrical constriction in the connector can potentially result in PEEP as obtained with conventional nasal cannulae.

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.

CFD simulations of inhalation through a subject-specific human larynx - Impact of the unilateral vocal fold immobility

BACKGROUND

The larynx of the human respiratory tract plays a vital role in breathing and voice production. Both can be influenced by functional and/or morphological changes of the larynx, e.g., immobility of one or both vocal folds (VF). The immobile VF can become stationary in different positions such as the median, paramedian, intermediate or lateral position. The impact of unilateral vocal fold immobility (UVFI) on inhalation is the focus of this study.

METHODS

Transient numerical simulations of the inhalation process in patient-specific airways are performed. Five configurations are considered: paramedian and intermediate VF positions on the left and right, and healthy. Large eddy simulations are used to describe the complex laryngeal turbulent flow. Airway resistance, power loss, and spectral entropy are calculated to quantify the work of inspiration and evaluate flow regimes.

RESULTS

The laryngeal jet intensity and flow disturbance increase with the severity of immobility. In comparison to the healthy configuration, UVFI with right/left intermediate and right/left paramedian VF position increases the airway resistance over the oropharynx to the trachea by 69%/58% and 310%/285%, respectively. When the entire respiratory system is considered, an increase of up to 48% is estimated. Spectral entropy increases of up to 2.5 times indicate higher turbulence levels due to UVFI.

CONCLUSIONS

Surgery of immobile VF aims to improve glottis closure. However, this can have a negative impact on breathing efficiency. To that end, this study provides initial insights into the conflicting objectives of open versus closed VFs.

Neonates With Tracheomalacia Generate Auto-Positive End-Expiratory Pressure via Glottis Closure

BACKGROUND

In pediatrics, tracheomalacia is an airway condition that causes tracheal lumen collapse during breathing and may lead to the patient requiring respiratory support. Adult patients can narrow their glottis to self-generate positive end-expiratory pressure (PEEP) to raise the pressure in the trachea and prevent collapse. However, auto-PEEP has not been studied in newborns with tracheomalacia. The objective of this study was to measure the glottis cross-sectional area throughout the breathing cycle and to quantify total pressure difference through the glottis in patients with and without tracheomalacia.

RESEARCH QUESTION

Do neonates with tracheomalacia narrow their glottises? How does the glottis narrowing affect the total pressure along the airway?

STUDY DESIGN AND METHODS

Ultrashort echo time MRI was performed in 21 neonatal ICU patients (11 with tracheomalacia, 10 without tracheomalacia). MRI scans were reconstructed at four different phases of breathing. All patients were breathing room air or using noninvasive respiratory support at the time of MRI. Computational fluid dynamics simulations were performed on patient-specific virtual airway models with airway anatomic features and motion derived via MRI to quantify the total pressure difference through the glottis and trachea.

RESULTS

The mean glottis cross-sectional area at peak expiration in the patients with tracheomalacia was less than half that in patients without tracheomalacia (4.0 ± 1.1 mm2 vs 10.3 ± 4.4 mm2; P = .002). The mean total pressure difference through the glottis at peak expiration was more than 10 times higher in patients with tracheomalacia compared with patients without tracheomalacia (2.88 ± 2.29 cm H2O vs 0.26 ± 0.16 cm H2O; P = .005).

INTERPRETATION

Neonates with tracheomalacia narrow their glottises, which raises pressure in the trachea during expiration, thereby acting as auto-PEEP.

Human upper-airway respiratory airflow: In vivo comparison of computational fluid dynamics simulations and hyperpolarized 129Xe phase contrast MRI velocimetry

Computational fluid dynamics (CFD) simulations of respiratory airflow have the potential to change the clinical assessment of regional airway function in health and disease, in pulmonary medicine and otolaryngology. For example, in diseases where multiple sites of airway obstruction occur, such as obstructive sleep apnea (OSA), CFD simulations can identify which sites of obstruction contribute most to airway resistance and may therefore be candidate sites for airway surgery. The main barrier to clinical uptake of respiratory CFD to date has been the difficulty in validating CFD results against a clinical gold standard. Invasive instrumentation of the upper airway to measure respiratory airflow velocity or pressure can disrupt the airflow and alter the subject's natural breathing patterns. Therefore, in this study, we instead propose phase contrast (PC) velocimetry magnetic resonance imaging (MRI) of inhaled hyperpolarized 129Xe gas as a non-invasive reference to which airflow velocities calculated via CFD can be compared. To that end, we performed subject-specific CFD simulations in airway models derived from 1H MRI, and using respiratory flowrate measurements acquired synchronously with MRI. Airflow velocity vectors calculated by CFD simulations were then qualitatively and quantitatively compared to velocity maps derived from PC velocimetry MRI of inhaled hyperpolarized 129Xe gas. The results show both techniques produce similar spatial distributions of high velocity regions in the anterior-posterior and foot-head directions, indicating good qualitative agreement. Statistically significant correlations and low Bland-Altman bias between the local velocity values produced by the two techniques indicates quantitative agreement. This preliminary in vivo comparison of respiratory airway CFD and PC MRI of hyperpolarized 129Xe gas demonstrates the feasibility of PC MRI as a technique to validate respiratory CFD and forms the basis for further comprehensive validation studies. This study is therefore a first step in the pathway towards clinical adoption of respiratory CFD.

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

Malignant central airway obstruction can cause severe breathing difficulty in a patient that requires surgical intervention or stent implantation to alleviate it. A predictive model to identify the onset of this event as the central airway is progressively compressed by tumor growth will be helpful for clinicians to plan for medical intervention. We present such a model to simulate tumor compression of the trachea and the resulting change in airflow dynamics to estimate the level of stenosis that will cause severe breathing difficulties. A patient-specific model of trachea was generated from acquired Computed Tomography (CT) scans for the simulations. The compression of this trachea due to tumor growth is modeled using nonlinear contact simulations of ellipsoidal tumors with the trachea. Computational fluid dynamics (CFD) is employed to simulate the turbulent airflow during inhalation in the stenosed trachea. From the CFD simulated flow fields, the power loss due to airflow through the domain is calculated. The results show that when the obstruction in the trachea reaches 50%, compared to the undeformed model, the power loss can rise to more than 66%. A measure of breathing difficulty can be derived by correlating it with the power loss. Thus, medical intervention can be predicted based on the degree of stenosis if the induced power loss exceeds a threshold that causes severe breathing discomfort.

Correlating Local Volumetric Tissue Strains with Global Lung Mechanics Measurements

The mechanics of breathing is a fascinating and vital process. The lung has complexities and subtle heterogeneities in structure across length scales that influence mechanics and function. This study establishes an experimental pipeline for capturing alveolar deformations during a respiratory cycle using synchrotron radiation micro-computed tomography (SR-micro-CT). Rodent lungs were mechanically ventilated and imaged at various time points during the respiratory cycle. Pressure-Volume (P-V) characteristics were recorded to capture any changes in overall lung mechanical behaviour during the experiment. A sequence of tomograms was collected from the lungs within the intact thoracic cavity. Digital volume correlation (DVC) was used to compute the three-dimensional strain field at the alveolar level from the time sequence of reconstructed tomograms. Regional differences in ventilation were highlighted during the respiratory cycle, relating the local strains within the lung tissue to the global ventilation measurements. Strains locally reached approximately 150% compared to the averaged regional deformations of approximately 80-100%. Redistribution of air within the lungs was observed during cycling. Regions which were relatively poorly ventilated (low deformations compared to its neighbouring region) were deforming more uniformly at later stages of the experiment (consistent with its neighbouring region). Such heterogenous phenomena are common in everyday breathing. In pathological lungs, some of these non-uniformities in deformation behaviour can become exaggerated, leading to poor function or further damage. The technique presented can help characterize the multiscale biomechanical nature of a given pathology to improve patient management strategies, considering both the local and global lung mechanics.

Subglottic Stenosis Position Affects Work of Breathing

OBJECTIVES

Subglottic stenosis (SGS) is the most common type of laryngeal stenosis in neonates. SGS severity is currently graded based on percent area of obstruction (%AO) via the Myer-Cotton grading scale. However, patients with similar %AO can have widely different clinical courses. Computational fluid dynamics (CFD) based on patient-specific imaging can quantify the relationship between airway geometry and flow dynamics. We investigated the effect of %AO and axial position of SGS on work of breathing (WOB) in neonates using magnetic resonance imaging.

METHODS

High-resolution ultrashort echo-time MRI of the chest and airway was obtained in three neonatal patients with no suspected airway abnormalities; images were segmented to construct three-dimensional (3D) models of the neonatal airways. These models were then modified with virtual SGSs of varying %AO and axial positioning. CFD simulations of peak inspiratory flow were used to calculate patient-specific WOB in nonstenotic and artificially stenosed airway models.

RESULTS

CFD simulations demonstrated a relationship between stenosis geometry and WOB increase. WOB rapidly increased with %AO greater than about 70%. Changes in axial position could also increase WOB by approximately the same amount as a 10% increase in %AO. Increased WOB was particularly pronounced when the SGS lumen was misaligned with the glottic jet.

CONCLUSION

The results indicate a strong, predictable relationship between WOB and axial position of the stenotic lumen relative to the glottis, which has not been previously reported. These findings may lead to precision diagnosis and treatment prediction tools in individual patients.

LEVEL OF EVIDENCE

4 Laryngoscope, 131:E1220-E1226, 2021.

Increased Work of Breathing due to Tracheomalacia in Neonates

Rationale: Dynamic collapse of the tracheal lumen (tracheomalacia) occurs frequently in premature neonates, particularly in those with common comorbidities such as bronchopulmonary dysplasia. The tracheal collapse increases the effort necessary to breathe (work of breathing [WOB]). However, quantifying the increased WOB related to tracheomalacia has previously not been possible. Therefore, it is also not currently possible to separate the impact of tracheomalacia on patient symptoms from parenchymal abnormalities.Objectives: To measure the increase in WOB due to airway motion in individual subjects with and without tracheomalacia and with different types of respiratory support.Methods: Fourteen neonatal intensive care unit subjects not using invasive mechanical ventilation were recruited. In eight, tracheomalacia was diagnosed via clinical bronchoscopy, and six did not have tracheomalacia. Self-gated three-dimensional ultrashort-echo-time magnetic resonance imaging (MRI) was performed on each subject with clinically indicated respiratory support to obtain cine images of tracheal anatomy and motion during the respiratory cycle. The component of WOB due to resistance within the trachea was then calculated via computational fluid dynamics (CFD) simulations of airflow on the basis of the subject's anatomy, motion, and respiratory airflow rates. A second CFD simulation was performed for each subject with the airway held static at its largest (i.e., most open) position to determine the increase in WOB due to airway motion and collapse.Results: The tracheal-resistive component of WOB was increased because of airway motion by an average of 337% ± 295% in subjects with tracheomalacia and 24% ± 14% in subjects without tracheomalacia (P < 0.02). In the tracheomalacia group, subjects who were treated with continuous positive airway pressure (CPAP) using a RAM cannula expended less energy for breathing compared with the subjects who were breathing room air or on a high-flow nasal cannula.Conclusions: Neonatal subjects with tracheomalacia have increased energy expenditure compared with neonates with normal airways, and CPAP may be able to attenuate the increase in respiratory work. Subjects with tracheomalacia expend more energy on the tracheal-resistive component of WOB alone than nontracheomalacia patients expend on the resistive WOB for the entire respiratory system, according to previously reported values. CFD may be able to provide an objective measure of treatment response for children with tracheomalacia.