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

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.

Effects of Hypoglossal Nerve Stimulation on Upper Airway Structure and Function Using Moving Wall Computational Fluid Dynamics Simulations: A Pilot Study

Hypoglossal nerve stimulation (HGNS) is an innovative alternative treatment option for obstructive sleep apnoea (OSA) in patients unable to tolerate continuous positive airway pressure. However, the success rate of HGNS is variable, but the reasons underlying variation in treatment efficacy are not well understood. In this pilot study of three male subjects, we propose an innovative, non-invasive method to quantify the structural and functional changes to the upper airway that occur with HGNS. We used four-dimensional computed tomography (4DCT) and computational fluid dynamics (CFD) simulations of respiratory airflow to quantify how HGNS changes: (1) airway cross-sectional area (CSA), (2) work done by muscles and air pressure in dilating and collapsing the airway and (3) airway resistance. Subjects underwent 4DCT under natural stage non-REM 2 (N2) sleep with and without HGNS. Each patient had concurrent electroencephalograms and airflow measurements. CFD simulations were performed based on anatomy and airway motion from 4DCT images and airflow data. HGNS was associated with an increase in neuromuscular work done in dilating the airway (up to 490%); airway CSA increased by up to 300%. Most motion with HGNS occurred in the oropharynx; changes in the nasopharynx and hypopharynx varied between subjects. Minute ventilation increased in all subjects (15%-36%). Airway resistance decreased across the three subjects (73%-97%). Quantifying the parameters measured in this study may help explain variable responses to HGNS as a treatment for OSA. These procedures may, in future, help predict non-responders to HGNS, isolate reasons for poor responses, or inform device titration.

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.

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.

The interaction between neuromuscular forces, aerodynamic forces, and anatomical motion in the upper airway predicts the severity of pediatric OSA

Upper airway neuromuscular response to air pressure during inhalation is an important factor in assessing pediatric subjects with obstructive sleep apnea (OSA). The neuromuscular response's strength, timing, and duration all contribute to the potential for airway collapses and the severity of OSA. This study quantifies these factors at the soft palate, tongue, and epiglottis to assess the relationship between neuromuscular control and OSA severity in 20 pediatric subjects with and without trisomy 21, under dexmedetomidine-induced sedation. The interaction between neuromuscular force and airflow pressure force was assessed based on power transferred between the airway wall and airflow calculated from airway wall motion (from cine magnetic resonance images) and air pressure acting on the airway wall (from computational fluid dynamics simulations). Airway wall motion could be asynchronous with pressure forces due to neuromuscular activation, or synchronous with pressure forces, indicating a passive response to airflow. The obstructive apnea-hypopnea index (oAHI) quantified OSA severity. During inhalation, the normalized work done through asynchronous dilation of the airway at the soft palate, tongue, and epiglottis correlated significantly with oAHI (Spearman's ρ = 0.54, 0.50, 0.64; P = 0.03, 0.03, 0.003). Synchronous collapse at the epiglottis correlated significantly with oAHI (ρ = 0.52; P = 0.02). Temporal order of synchronous and asynchronous epiglottis motion during inhalation predicted the severity of OSA (moderate vs. severe) with 100% sensitivity and 70% specificity. Subjects with severe OSA and/or trisomy 21 have insufficient neuromuscular activation during inhalation, leading to collapse and increased neuromuscular activation. Airflow-driven airway wall motion during late inhalation likely is the main determinant of OSA severity.NEW & NOTEWORTHY This is the first study that combines cine MRI and computational fluid dynamics with in vivo synchronous respiratory flow measurement to quantify the interaction between airway neuromuscular forces, aerodynamic forces, and airway anatomy noninvasively in pediatric patients with obstructive sleep apnea (OSA). The results indicate power transfer predicts OSA severity.

Computational fluid-particle dynamics modeling of ultrafine to coarse particles deposition in the human respiratory system, down to the terminal bronchiole

BACKGROUND AND OBJECTIVES

Suspended respirable airborne particles are associated with human health risks and especially particles within the range of ultrafine (< 0.1 μm) or fine (< 2.5 μm) have a high possibility of penetrating the lung region, which is concerned to be closely related to the bronchial or alveoli tissue dosimetry. Nature complex structure of the respiratory system requires much effort to explore and comprehend the flow and the inhaled particle dynamics for precise health risk assessment. Therefore, this study applied the computational fluid-particle dynamics (CFPD) method to elucidate the deposition characteristics of ultrafine-to-coarse particles in the human respiratory tract from nostrils to the 16^th generation of terminal bronchi.

METHODS

The realistic bronchi up to the 8^th generation are precisely and perfectly generated from computed tomography (CT) images, and an artificial model compensates for the 9^th-16^th bronchioles. Herein, the steady airflow is simulated at constant breathing flow rates of 7.5, 15, and 30 L/min, reproducing human resting-intense activity. Then, trajectories of the particle size ranging from 0.002 - 10 μm are tracked using a discrete phase model.

RESULTS

Here, we report reliable results of airflow patterns and particle deposition efficiency in the human respiratory system validated against experimental data. The individual-related focal point of ultrafine and fine particles deposition rates was actualized at the 8^th generation; whilst the hot-spot of the deposited coarse particles was found in the 6^th generation. Lobar deposition characterizes the dominance of coarse particles deposited in the right lower lobe, whereas the left upper-lower and right lower lobes simultaneously occupy high deposition rates for ultrafine particles. Finally, the results indicate a higher deposition in the right lung compared to its counterpart.

CONCLUSIONS

From the results, the developed realistic human respiratory system down to the terminal bronchiole in this study, in coupling with the CFPD method, delivers the accurate prediction of a wide range of particles in terms of particle dosimetry and visualization of site-specific in the consecutive respiratory system. In addition, the series of CFPD analyses and their results are to offer in-depth information on particle behavior in human bronchioles, which may benefit health risk assessment or drug delivery studies.

Nichtinvasive funktionelle Lungenbildgebung mit hyperpolarisiertem Xenon

Hintergrund

Die Magnetresonanztomographie (MRT) ist ein nichtinvasives Verfahren mit hervorragendem Weichteilkontrast. Aufgrund der geringen Protonendichte und vielen Luft-Gewebe-Übergängen ist die Anwendung in der Lunge jedoch eingeschränkt, so dass hier häufig röntgenbasierte Methoden eingesetzt werden (mit den bekannten Nachteilen ionisierender Strahlung).

Fragestellung

In dieser Übersichtsarbeit wird die Lungen-MRT mit hyperpolarisiertem Xenon-129 (Xe-MRT) dargestellt. Die Xe-MRT erlaubt einzigartige wertvolle Einblicke in die Mikrostruktur und Funktion der Lunge, einschließlich des Gasaustauschs mit roten Blutkörperchen – Parameter, die mit klinischen Standardmethoden nicht zugänglich sind.

Material und Methoden

Durch die magnetische Markierung, die Hyperpolarisierung, wird das Signal von Xenon-129 um bis zu 100.000-fach verstärkt. Hierbei werden die Elektronen von Rubidium mittels Laserlicht zunächst auf 100 % polarisiert und dann durch Stöße auf Xenon übertragen. Danach wird das hyperpolarisierte Gas in einem Beutel zum Patienten gebracht und eingeatmet, kurz bevor die MRT-Aufnahmen beginnen.

Ergebnisse

Durch spezielle Programmierungen (Sequenzen) in der MRT kann die Ventilation, Mikrostruktur oder der Gasaustausch der Lunge in 3‑D dargestellt werden. Dies ermöglicht z. B. die quantitative Darstellung von Belüftungsdefekten, der Größe der Alveolen, der Gasaufnahme im Gewebe und des Gastransfers ins Blut.

Schlussfolgerung

Die Xe-MRT liefert einzigartige Informationen über den Zustand der Lunge – nichtinvasiv, in vivo und in weniger als einer Minute.

Biomedical and biophysical limits to mathematical modeling of pulmonary system mechanics: a scoping review on aerosol and drug delivery

Undoubtedly, the construction of the biomechanical geometry systems with the help of computer tomography (CT) and magnetic resonance imaging (MRI) has made a significant advancement in studying in vitro numerical models as accurately as possible. However, some simplifying assumptions in the computational studies of the respiratory system have caused errors and deviations from the in vivo actual state. The most important of these hypotheses is how to generate volume from the point cloud exported from CT or MRI images, not paying attention to the wall thickness and its effect in computational fluid dynamic method, statistical logic of aerosol trap in software; and most importantly, the viscoelastic effect of respiratory tract wall in living tissue pointed in the fluid-structure interaction method. So that applying the viscoelastic dynamic mesh effect in the form of the moving deforming mesh can be very effective in achieving more appropriate response quality. Also, changing the volume fraction of the pulmonary extracellular matrix constituents leads to changes in elastic modulus (storage modulus) and the viscous modulus (loss modulus) of lung tissue. Therefore, in the biomedical computational methods where the model wall is considered flexible, the viscoelastic properties of the texture must be considered correctly.