A novel method to generate dynamic boundary conditions for airway CFD by mapping upper airway movement with non-rigid registration of dynamic and static MRI

Effects of the deformation and size of the upper airway on the deposition of aerosols

Aerosol drug delivery in the human airway is significantly affected by the morphology and size of the airway. This work developed a CFD-DEM model to simulate and analyze air flow and powder dynamics in combined inhaler-airway systems with different degrees of airway deformation (non-deformed, 50%, and 75% deformed) and sizes (adult, 0.80, and 0.62 scaled). The airways were generated based on a regular airway constructed from the MRI images through finite element method (for deformed airways) or scaling-down (for smaller airways). The airways were connected to Turbuhaler® through a connector. The results showed that under the same flow rate, the variation in the airway geometry and size had a minimum impact on the flow field and powder deposition in the device and the connector. However, deformation caused more particle deposition in the deformed region. Notably, the airway with 50% deformation had the most particles passing through the airway with the largest particle sizes due to its lower air velocity in the deformed area. Reducing airway size resulted in more powder deposition on the airway, particularly at the pharynx and mouth regions. This was because, with the same flow rate, the flow velocity in the smaller airway was higher, causing more particle-wall collisions in the mouth and pharynx regions. More importantly, the deposition efficiency in the 0.62-scaled airway was significantly higher than the other two airways, highlighting the importance of the different administration of aerosol drugs for young children.

Comparison of weighting algorithms to mitigate respiratory motion in free-breathing neonatal pulmonary radial UTE-MRI

Background. Thoracoabdominal MRI is limited by respiratory motion, especially in populations who cannot perform breath-holds. One approach for reducing motion blurring in radially-acquired MRI is respiratory gating. Straightforward 'hard-gating' uses only data from a specified respiratory window and suffers from reduced SNR. Proposed 'soft-gating' reconstructions may improve scan efficiency but reduce motion correction by incorporating data with nonzero weight acquired outside the specified window. However, previous studies report conflicting benefits, and importantly the choice of soft-gated weighting algorithm and effect on image quality has not previously been explored. The purpose of this study is to map how variable soft-gated weighting functions and parameters affect signal and motion blurring in respiratory-gated reconstructions of radial lung MRI, using neonates as a model population.Methods. Ten neonatal inpatients with respiratory abnormalities were imaged using a 1.5 T neonatal-sized scanner and 3D radial ultrashort echo-time (UTE) sequence. Images were reconstructed using ungated, hard-gated, and several soft-gating weighting algorithms (exponential, sigmoid, inverse, and linear weighting decay outside the period of interest), with %Nprojrepresenting the relative amount of data included. The apparent SNR (aSNR) and motion blurring (measured by the maximum derivative of image intensity at the diaphragm, MDD) were compared between reconstructions.Results. Soft-gating functions produced higher aSNR and lower MDD than hard-gated images using equivalent %Nproj, as expected. aSNR was not identical between different gating schemes for given %Nproj. While aSNR was approximately linear with %Nprojfor each algorithm, MDD performance diverged between functions as %Nprojdecreased. Algorithm performance was relatively consistent between subjects, except in images with high noise.Conclusion. The algorithm selection for soft-gating has a notable effect on image quality of respiratory-gated MRI; the timing of included data across the respiratory phase, and not simply the amount of data, plays an important role in aSNR. The specific soft-gating function and parameters should be considered for a given imaging application's requirements of signal and sharpness.

Impact of palatopharyngeal sizes changing on pharyngeal airflow fluctuation and airway vibration in a pediatric airway

Snoring is common in children and is associated with many adverse consequences. One must study the relationships between pharyngeal morphology and snoring physics to understand snoring progression. Although some model studies have provided fluid-structure interaction dynamic descriptions for the correlation between airway size and snoring physics, the descriptions still need to be further investigated in patient-specific airway models. Fluid-structure interaction studies using patient-specific airway structures complement the above model studies. Based on reported cephalometric measurement methods, this study quantified and preset the size of the palatopharynx airway in a patient-specific airway and investigated how the palatopharynx size affects the pharyngeal airflow fluctuation, soft palate vibration, and glossopharynx vibration with the help of a verified FSI method. The results showed that the stenosis anterior airway of the soft palate increased airway resistance and airway resistance fluctuations, which can lead to increased sleep effort and frequent snoring. Widening of the anterior airway can reduce airflow resistance and avoid obstructing the anterior airway by the soft palate vibration. The pharyngeal airflow resistance, mouth inflow proportion, and soft palate apex displacement have components at the same frequencies in all airway models, and the glossopharynx vibration and instantaneous inflow rate have components at the same frequencies, too. The mechanism of this same frequency fluctuation phenomenon can be explained by the fluid-structure interaction dynamics of an ideal coupled model consisting of a flexible plate model and a collapsible tube model. The results of this study demonstrate the potential of FSI in studying snoring physics and clarify to some degree the mechanism of airway morphology affecting airway vibration physics.

Scale resolving simulations of the effect of glottis motion and the laryngeal jet on flow dynamics during respiration

BACKGROUND AND OBJECTIVE

The movement of the respiratory walls has a significant impact on airflow through the respiratory tract. The majority of computational fluid dynamics (CFD) studies assume a static geometry which may not provide a realistic flow field. Furthermore, many studies use Reynolds Averaged Navier-Stokes (RANS) turbulence models that do not resolve turbulence structure. Combining the application of advanced scale-resolving turbulence models with moving respiratory walls using CFD will provide detailed insights into respiratory flow structures.

METHODS

This study simulated a complete breathing cycle involving inhalation and exhalation in a nasal cavity to trachea geometry that incorporated moving glottis walls. A second breathing cycle was simulated with static glottis walls for comparison. A recently developed hybrid RANS-LES turbulence model, the Stress-Blended Eddy Simulation (SBES), was incorporated to resolve turbulent flow structures in fine detail for both transient simulations. Transient results were compared with steady-state RANS simulations for the same respiratory geometry.

RESULTS

Glottis motion caused substantial effects on flow structure through the complete breathing cycle. Significant flow structure and velocity variations were observed due to glottal motion, primarily in the larynx and trachea. Resolved turbulence structures using SBES showed an intense mixing section in the glottis region during inhalation and in the nasopharynx during expiration, which was not present in the RANS simulations.

CONCLUSION

Transient simulations of a realistic breathing cycle uncovered flow structures absent in simulations with a constant flow rate. Furthermore, the incorporation of glottis motion impacted airflow characteristics that suggest rigid respiratory walls do not accurately describe respiratory flow. Future research in respiratory airflow should be conducted using transient scale-resolving models in conjunction with moving respiratory walls to capture flow structures in detail.

Numerical simulation of flow behavior in basilar bifurcation computed tomography angiography

In this study, a moving boundary deformation model based on four-dimensional computed tomography angiography (4D-CTA) with high temporal resolution is constructed, and blood flow dynamics of cerebral aneurysms are investigated by numerical simulation. A realistic moving boundary deformation model of a cerebral aneurysm was constructed based on 4D-CTA in each phase. Four hemodynamic factors (wall shear stress [WSS], wall shear stress divergence [WSSD], oscillatory shear index [OSI], and residual residence time [RRT]) were obtained from numerical simulations, and these factors were evaluated in basilar artery aneurysms. Comparison of the rigid body condition and the moving boundary condition investigating the relationship between wall displacement and hemodynamic factors clarified that the spatial-averaged WSS and maximum WSSD considering only the aneurysmal dome has a large difference between conditions during the peak systole, and there were also significant differences in OSI and RRT.

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.

A flexible 16-channel custom coil array for accelerated imaging of upper and infraglottic airway at 3 T

PURPOSE

To develop a custom coil and evaluate its utility for accelerated upper and infraglottic airway MRI at 3 T.

METHODS

A 16-channel flexible and anatomy-conforming coil was developed to provide localized sensitivity over upper and infraglottic airway regions of interest. Parallel-imaging capabilities were compared against existing head and head-neck coils. SENSE geometry factor losses were quantified for retrospectively accelerating 3D MRI. Blinded image-quality ratings from two experts were performed. Spiral GRAPPA reconstructions were evaluated for a speaking task at a time resolution of 40 ms. Contrast-to-noise ratios between air and tissue at key landmarks along the vocal tract were compared. SENSE imaging with the custom coil in the lateral recumbent posture was evaluated. Multislice imaging was performed to image swallowing at 17 ms/frame via constrained reconstruction.

RESULTS

The custom coil showed improved SENSE imaging up to 3-fold acceleration when accelerated along either the anterior-posterior or the superior-inferior direction and a net 4-fold acceleration when accelerated along both directions. Spiral GRAPPA reconstructions with the custom coil showed higher contrast-to-noise ratio when compared with existing coils. In the lateral posture, robust SENSE imaging was achieved at up to 2-fold and 3-fold acceleration levels in the superior-inferior and anterior-posterior directions, respectively. Key events of swallowing in the multislice dynamic images were identified by an otolaryngologist.

CONCLUSION

The coil provided improved parallel imaging of upper and infraglottic airway in both supine and lateral recumbent postures. It enabled efficient accelerated dynamic imaging of speaking and swallowing.

Impact of sleep posture and breathing pattern on soft palate flutter and pharynx vibration in a pediatric airway using fluid-structure interaction

Snoring is a common condition in the general population, and the management of snoring requires a better understanding of its mechanism through a fluid-structure interaction (FSI) perspective. Despite the recent popularity of numerical FSI techniques, outstanding challenges are accurately predicting airway deformation and its vibration during snoring due to complex airway morphology. In addition, there still needs to be more understanding of snoring inhibition when lying on the side, and the possible effect of airflow rates, as well as nose or mouth-nose breathing, on snoring remains to be investigated. In this study, an FSI method verified against in vitro models was introduced to predict upper airway deformation and vibration. The technique was applied to predict airway aerodynamics, soft palate flutter, and airway vibration in four sleep postures (supine, left/right lying, and sitting positions) and four breathing patterns (mouth-nose, nose, mouth, and unilateral nose breathing). It was found that, at given elastic properties of soft tissues, the evaluated flutter frequency of 19.8 Hz in inspiration was in good agreement with the reported frequency of snoring sound in literature. Reduction in flutter and vibrations due to the mouth-nose airflow proportion changes were also noticed when having side-lying and sitting positions. Breathing through the mouth results in larger airway deformation than breathing through the nose or mouth-nose. These results collectively demonstrate the potential of FSI for studying the physics of airway vibration and clarify to some degree the reason for snoring inhibition during sleep postures and breathing patterns.

Computational assessment of upper airway muscular activity in obstructive sleep apnea - In vitro validation

Neuromuscular control of the upper airway contributes to obstructive sleep apnea (OSA). An accurate, non-invasive method to assess neuromuscular function is needed to improve surgical treatment outcomes. Currently, surgical approaches for OSA are based on airway anatomy and are often not curative. When the airway surface moves, the power transferred between air in the airway lumen and the structures of the upper airway may be a measure of airway neuromuscular activity. The aim of this study was to validate power transfer as a measure of externally applied forces, representing neuromuscular activity, through cine computed tomography (CT) imaging and computational fluid dynamics (CFD) analysis in a 3D-printed airway model. A hollow elastic airway model was manufactured. An insufflation/exsufflation device generated airflow within the model lumen. The model was contained in an airtight chamber that could be positively or negatively pressurized to represent muscular forces. These forces were systematically applied to dilate and collapse the model. Cine CT imaging captured airway wall movement during respiratory cycles with and without externally applied forces. Power transfer was calculated from the product of wall movement and internal aerodynamic pressure forces using CFD simulations. Cross-correlation peaks between power transfer and changes in externally applied pressure during exhalation and inhalation were -0.79 and 0.95, respectively. Power transfer calculated via cine CT imaging and CFD was an accurate surrogate measure of externally applied forces representing airway muscular activity. In the future, power transfer may be used in clinical practice to phenotype patients with OSA and select personalized therapies.

Predicting tracheal work of breathing in neonates based on radiological and pulmonary measurements

Tracheomalacia is an airway condition in which the trachea excessively collapses during breathing. Neonates diagnosed with tracheomalacia require more energy to breathe, and the effect of tracheomalacia can be quantified by assessing flow-resistive work of breathing (WOB) in the trachea using computational fluid dynamics (CFD) modeling of the airway. However, CFD simulations are computationally expensive; the ability to instead predict WOB based on more straightforward measures would provide a clinically useful estimate of tracheal disease severity. The objective of this study is to quantify the WOB in the trachea using CFD and identify simple airway and/or clinical parameters that directly relate to WOB. This study included 30 neonatal intensive care unit subjects (15 with tracheomalacia and 15 without tracheomalacia). All subjects were imaged using ultrashort echo time (UTE) MRI. CFD simulations were performed using patient-specific data obtained from MRI (airway anatomy, dynamic motion, and airflow rates) to calculate the WOB in the trachea. Several airway and clinical measurements were obtained and compared with the tracheal resistive WOB. The maximum percent change in the tracheal cross-sectional area (ρ = 0.560, P = 0.001), average glottis cross-sectional area (ρ = -0.488, P = 0.006), minute ventilation (ρ = 0.613, P < 0.001), and lung tidal volume (ρ = 0.599, P < 0.001) had significant correlations with WOB. A multivariable regression model with three independent variables (minute ventilation, average glottis cross-sectional area, and minimum of the eccentricity index of the trachea) can be used to estimate WOB more accurately (R2 = 0.726). This statistical model may allow clinicians to estimate tracheal resistive WOB based on airway images and clinical data.NEW & NOTEWORTHY The work of breathing due to resistance in the trachea is an important metric for quantifying the effect of tracheal abnormalities such as tracheomalacia, but currently requires complex dynamic imaging and computational fluid dynamics simulation to calculate it. This study produces a method to predict the tracheal work of breathing based on readily available imaging and clinical metrics.

Nasal Irrigation Delivery in Three Post-FESS Models From a Squeeze-bottle Using CFD

Purpose

Nasal saline irrigation is highly recommended in patients following functional endoscopic sinus surgery (FESS) to aid the postoperative recovery. Post-FESS patients have significantly altered anatomy leading to markedly different flow dynamics from those found in pre-op or non-diseased airways, resulting in unknown flow dynamics.

Methods

This work investigated how the liquid stream disperses through altered nasal cavities following surgery using Computational Fluid Dynamics (CFD). A realistic squeeze profile was determined from physical experiments with a 27-year-old male using a squeeze bottle with load sensors. The administration technique involved a head tilt of 45-degrees forward to represent a head position over a sink. After the irrigation event that lasted 4.5 s, the simulation continued for an additional 1.5 s, with the head orientation returning to an upright position.

Results

The results demonstrated that a large maxillary sinus ostium on the right side allows saline penetration into this sinus. The increased volume of saline entering the maxillary sinus limits the saline volume available to the rest of the sinonasal cavity and reduces the surface coverage of the other paranasal sinuses. The average wall shear stress was higher on the right side than on the other side for two patients. The results also revealed that head position alters the sinuses’ saline residual, especially the frontal sinuses.

Conclusion

While greater access to sinuses is achieved through FESS surgery, patients without a nasal septum limits posterior sinus penetration due to the liquid crossing over to the contralateral cavity and exiting the nasal cavity early.

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.

Combining experimental and computational techniques to understand and improve dry powder inhalers

INTRODUCTION

: Dry Powder Inhalers (DPIs) continue to be developed to deliver an expanding range of drugs to treat an ever-increasing range of medical conditions; with each drug and device combination needing a specifically designed inhaler. Fast regulatory approval is essential to be first to market, ensuring commercial profitability.

AREAS COVERED

: In vitro deposition, particle image velocimetry, and computational modelling using the physiological geometry and representative anatomy can be combined to give complementary information to determine the suitability of a proposed inhaler design and to optimise its formulation performance. In combination they allow the entire range of questions to be addressed cost-effectively and rapidly.

EXPERT OPINION

: Experimental techniques and computational methods are improving rapidly, but each needs a skilled user to maximize results obtained from these techniques. Multidisciplinary teams are therefore key to making optimal use of these methods and such qualified teams can provide enormous benefits to pharmaceutical companies to improve device efficacy and thus time to market. There is already a move to integrate the benefits of Industry 4.0 into inhaler design and usage, a trend that will accelerate.

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.

Oral Positive Expiratory Pressure Device for Excessive Dynamic Airway Collapse Caused by Emphysema

Excessive dynamic airway collapse (EDAC) contributes to breathlessness and reduced quality of life in individuals with emphysema. We tested a novel, portable, oral positive expiratory pressure (o-PEP) device in a patient with emphysema and EDAC. MRI revealed expiratory tracheal narrowing to 80 mm² that increased to 170 mm² with the o-PEP device. After 2-weeks use of the o-PEP device for 33% to 66% of activities, breathlessness, quality of life, and exertional dyspnea improved compared with minimal clinically important differences (MCID): University of California-San Diego Shortness of Breath questionnaire score declined 69 to 42 (MCID, ≥5), St. George's Respiratory Questionnaire score decreased 71 to 27 (MCID, ≥4), and before and after the 6-minute walk test Borg score difference improved from Δ3 to Δ2 (MCID, ≥1). During the 6-minute walk test on room air without the use of the o-PEP device, oxyhemoglobin saturation declined 91% to 83%; whereas, with the o-PEP device, the nadir was 90%. Use of the o-PEP device reduced expiratory central airway collapse and improved dyspnea, quality of life, and exertional desaturation in a patient with EDAC and emphysema.

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.

The effect of airway motion and breathing phase during imaging on CFD simulations of respiratory airflow

RATIONALE

Computational fluid dynamics (CFD) simulations of respiratory airflow can quantify clinically useful information that cannot be obtained directly, such as the work of breathing (WOB), resistance to airflow, and pressure loss. However, patient-specific CFD simulations are often based on medical imaging that does not capture airway motion and thus may not represent true physiology, directly affecting those measurements.

OBJECTIVES

To quantify the variation of respiratory airflow metrics obtained from static models of airway anatomy at several respiratory phases, temporally averaged airway anatomies, and dynamic models that incorporate physiological motion.

METHODS

Neonatal airway images were acquired during free-breathing using 3D high-resolution MRI and reconstructed at several respiratory phases in two healthy subjects and two with airway disease (tracheomalacia). For each subject, five static (end expiration, peak inspiration, end inspiration, peak expiration, averaged) and one dynamic CFD simulations were performed. WOB, airway resistance, and pressure loss across the trachea were obtained for each static simulation and compared with the dynamic simulation results.

RESULTS

Large differences were found in the airflow variables between the static simulations at various respiratory phases and the dynamic simulation. Depending on the static airway model used, WOB, resistance, and pressure loss varied up to 237%, 200%, and 94% compared to the dynamic simulation respectively.

CONCLUSIONS

Changes in tracheal size and shape throughout the breathing cycle directly affect respiratory airflow dynamics and breathing effort. Simulations incorporating realistic airway wall dynamics most closely represent airway physiology; if limited to static simulations, the airway geometry must be obtained during the respiratory phase of interest for a given pathology.

Pediatric obstructive sleep apnea revisited: Perioperative considerations for the pediatric Anesthesiologist

Pediatric obstructive sleep apnea presents in up to 7% of children and represents a constellation from nasal turbulence to cessation in gas exchange. There are numerous end organ sequelae including neurocognitive morbidity associated with persistent OSA. Adenotonsillectomy (AT), the first line therapy for pediatric OSA, has not been demonstrated to reduce all end organ morbidity, specifically neurological and behavioral morbidity. Furthermore, certain at-risk populations are at higher risk from neurocognitive morbidity. Precise knowledge and perioperative planning is required to ensure optimal evidence-based practices in children with OSA. This comprehensive review covers the seminal perioperative implications of OSA, including preoperative polysomnography, pharmacotherapeutics, and postoperative risk stratification.

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.

Ultrashort Echo-Time MRI for the Assessment of Tracheomalacia in Neonates

BACKGROUND

Bronchoscopy is the gold standard for evaluating tracheomalacia; however, reliance on an invasive procedure limits understanding of normal airway dynamics. Self-gated ultrashort echo-time MRI (UTE MRI) can assess tracheal dynamics but has not been rigorously evaluated.

METHODS

This study was a validation of UTE MRI diagnosis of tracheomalacia in neonates using bronchoscopy as the gold standard. Bronchoscopies were reviewed for the severity and location of tracheomalacia based on standardized criteria. The percent change in cross-sectional area (CSA) of the trachea between end-inspiration and end-expiration was determined by UTE MRI, and receiver-operating curves were used to determine the optimal cutoff values to predict tracheomalacia and determine positive and negative predictive values.

RESULTS

Airway segments with tracheomalacia based on bronchoscopy had a more than threefold change in CSA measured from UTE MRI (54.4 ± 56.1% vs 14.8 ± 19.5%; P < .0001). UTE MRI correlated moderately with bronchoscopy for tracheomalacia severity (ρ = 0.39; P = .0001). Receiver-operating curves, however, showed very good ability of UTE MRI to identify tracheomalacia (area under the curve, 0.78). A "loose" definition (> 20% change in CSA) of tracheomalacia had good sensitivity (80%) but low specificity (64%) for identifying tracheomalacia based on UTE MRI, whereas a "strict" definition (> 40% change in CSA) was poorly sensitive (48%) but highly specific (93%).

CONCLUSIONS

Self-gated UTE MRI can noninvasively assess tracheomalacia in neonates without sedation, ionizing radiation, or increased risk. This technique overcomes major limitations of other diagnostic modalities and may be suitable for longitudinal population studies of tracheal dynamics.