Assessing the relationship between movement and airflow in the upper airway using computational fluid dynamics with motion determined from magnetic resonance imaging

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.

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.

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.

Hypopharyngeal geometry impact on air-induced loads on the supraglottis

Exercise-induced laryngeal obstruction (EILO) describes paradoxical laryngeal closure during inspiration at high-intensity exercise. It is hypothesised that during intense activity, the air-induced loads on supraglottic walls overcome their internal stiffness, leading to the obstruction. Recent investigations have revealed that the air-induced loads on the supraglottic walls vary nonlinearly with increasing flow rate. It is, however, unclear whether certain geometric configurations of the hypopharynx and larynx may contribute to the predisposition to EILO. This study investigates the influence of hypopharyngeal and laryngeal geometry on upper respiratory tract airflow and air-induced forces. A computational fluid dynamics model is developed to study airflow through larynx. Four real, adult upper respiratory tracts with variable configurations are considered. Two steady, uniform inspiratory flow rates of 60 L/min and 180 L/min are considered. The analysis shows that geometries with a space lateral to the epiglottis (EpiS) and piriform fossae (PF) directs the hypopharyngeal and supraglottic pressure field to remain positive and increase with the flow rate. In geometries with EpiS and PF, pressure differential occurs around the aryepiglottic fold producing a net inward force over the region. The three-fold increase in flow rate induces near ten-fold increases in force over the region which may facilitate the closure. It is concluded that hypopharyngeal anatomy, particularly the piriform fossae, play a significant role in the obstruction of the supraglottic airway and should be considered in research and clinical assessment of EILO.

Computational fluid dynamics comparison of the upper airway velocity, pressure, and resistance in cats using an endotracheal tube or a supraglottic airway device

Intoduction

In veterinary medicine, airway management of cats under general anesthesia is performed with an endotracheal tube (ETT) or supraglottic airway device (SGAD). This study aims to describe the use of computational fluid dynamics (CFD) to assess the velocities, pressures, and resistances of cats with ETT or SGAD.

Methods

A geometrical reconstruction model of the device, trachea, and lobar bronchi was carried out from computed tomography (CT) scans that include the head, neck, and thorax. Twenty CT scans of cats under general anesthesia using ETT (n = 10) and SGAD (n = 10) were modeled and analyzed. An inspiratory flow of 2.4 L/min was imposed in each model and velocity (m/s), general and regional pressures (cmH2O) were computed. General resistance (cmH2O/L/min) was calculated using differential pressure differences between the device inlet and lobar bronchi. Additionally, regional resistances were calculated at the device's connection with the breathing circuit (region A), at the glottis area for the SGAD, and the area of the ETT exit (bevel) (region B) and the device itself (region C).

Results

Recirculatory flow and high velocities were found at the ETT's bevel and at the glottis level in the SGAD group. The pressure gradient (Δp) was more enhanced in the ETT cases compared with the SGAD cases, where the pressure change was drastic. In region A, the Δp was higher in the ETT group, while in regions B and C, it was higher in the SGAD group. The general resistance was not statistically significant between groups (p = 0.48). Higher resistances were found at the region A (p = <0.001) in the ETT group. In contrast, the resistance was higher in the SGAD cases at the region B (p = 0.001).

Discussion

Overall, the provided CT-based CFD analysis demonstrated regional changes in airway pressure and resistance between ETT and SGAD during anesthetic flow conditions. Correct selection of the airway device size is recommended to avoid upper airway obstruction or changes in flow parameters.

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.

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.

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.

A review of upper airway physiology relevant to the delivery and deposition of inhalation aerosols

Developing effective oral inhaled drug delivery treatment strategies for respiratory diseases necessitates a thorough knowledge of the respiratory system physiology, such as the differences in the airway channel's structure and geometry in health and diseases, their surface properties, and mechanisms that maintain their patency. While respiratory diseases, such as chronic obstructive pulmonary disease (COPD) and asthma and their implications on the lower airways have been the core focus of most of the current research, the role of the upper airway in these diseases is less known, especially in the context of inhaled drug delivery. This is despite the fact that the upper airway is the passageway for inhaled drugs to be delivered to the lower airways, and their replicas are indispensable in current standards, such as the cascade impactor experiments for testing inhaled drug delivery technology. This review provides an overview of upper airway collapsibility and their mechanical properties, the effects of age and gender on upper airway geometry, and surface properties. The review also discusses how COPD and asthma affect the upper airway and the typical inhalation flow characteristics exhibited by the patients with these diseases.

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.

Evaluation of the Polyhedral Mesh Style for Predicting Aerosol Deposition in Representative Models of the Conducting Airways

A critical factor affecting the accuracy of Computational Fluid Dynamic (CFD) simulations and the time required to conduct them is construction of the computational mesh. This study aimed to evaluate the relatively new polyhedral mesh style for simulating aerosol deposition in the upper conducting airways compared with established meshing techniques and experimental data. Hexahedral and polyhedral mesh solutions were compared in two benchmark geometries: 1) a 90°-bend with flow characteristics similar to the extrathoracic airways of an adolescent child, and 2) a double bifurcation representing bifurcations B3-B5 in an adult. Both 4-block and 5-block hexahedral meshes were used in the 90°-bend to capture the potential of fully-structured hexahedral meshes. In the 90°-bend, polyhedral elements matched polydisperse in vitro deposition data with 20% relative error (RE; averaged across the particle sizes considered), which is an improvement on the accuracy of the 4-block hexahedral mesh (35% RE) and is similar to the accuracy of the 5-block hexahedral mesh (19% RE). In the double bifurcation, deposition fraction relative differences evaluated between polyhedral and hexahedral meshes ranged from 0.3% to 28.6% for the different particle sizes assessed, which is an order of magnitude improvement compared with previous studies that considered hexahedral vs. hybrid tetrahedral-prism meshes for the same flow field. Solution convergence time with polyhedral elements was found to be 50% to 140% higher than with hexahedral meshes of comparable size. While application dependent, the increase in simulation time observed with polyhedral meshes will likely be outweighed by the ease and convenience of polyhedral mesh construction. It was concluded that the polyhedral mesh style, with sufficient resolution especially near the walls, is an excellent alternative to the highly regarded hexahedral mesh style for predicting upper airway aerosol transport and deposition and provides a powerful new tool in the assessment of respiratory aerosol dosimetry.

A computational model of upper airway respiratory function with muscular coupling

A two dimensional finite element model of upper airway respiratory function was developed emphasizing the effects of dilator muscular activation on the human retro-lingual airway. The model utilized an upright mid-sagittal computed tomography of the human head and neck to reconstruct relevant structures of the tongue, mandible, and the hyoid-related soft tissues, along with the retro-lingual airway. The reconstructed geometry was divided into fluid and solid domains and discretized into finite element (FE) meshes used for the computational model. Three cases were investigated: standing position; supine position; and supine position coupled with dilator muscle activation. Computations were performed for the inspiration stage of the breathing cycle, utilizing a fluid-structure interaction (FSI) method to couple structural deformation with airflow dynamics. The spatio-temporal deformation of the structures surrounding the airway wall were predicted to be in general agreement with known changes from upright to supine posture on luminal opening, as well as the distribution of airflow. The model effectively captured the effects of muscular stimulation on the upper airway anatomical changes, the flow characteristics relevant to airway reduction in the supine position and airway enlargement with muscle activation. The smallest airway opening in the retro-lingual section is predicted to occur at the epiglottic region in all the three cases considered, an unexpected vulnerable location of airway obstruction. The model also predicted that hyoid displacement would be associated with recovery from airway collapse. This information may be useful for building more complex models relevant to mechanisms and clinical interventions for obstructive sleep apnea.

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.

Effect of tube length on the buckling pressure of collapsible tubes

BACKGROUND

The higher incidence of obstructive sleep apnea (OSA) in men than in women has been attributed to the upper airway being longer in men. The Starling resistor is the paradigm biomechanical model of upper airway collapse in OSA where a collapsible tube (representing the pharynx) is located between two rigid tubes (representing the nasal cavity and trachea). While the Starling resistor has been extensively studied due to its relevance to many physiological phenomena, the effect of tube length on tube collapsibility has not been quantified yet.

METHODS

Finite element analysis of a 3-dimensional collapsible tube subjected to a transmural pressure was performed in ANSYS Workbench. The numerical methods were validated with in vitro experiments in a silicone tube whose modulus of elasticity (361 ± 28 kPa) and dimensions (length = 100 mm, diameter = 22.2 mm, and wall thickness = 1.59 mm) were selected so that tube compliance was similar to pharyngeal compliance in humans during sleep. The buckling pressure (transmural pressure at which the tube collapses) was quantified in tubes of three different diameters (10 mm, 16 mm, and 22.2 mm) and ten length-to-diameter ratios (L/D = 4 to 13), while keeping the wall-thickness-to-radius ratio constant at 0.143.

RESULTS

The absolute value of the buckling pressure decreased from 4.7 to 3.3 cmH₂O (461-324 Pa) when L/D increased from 4 to 13. The buckling pressure was nearly independent from tube length for L/D >10.

CONCLUSIONS

Our finding that longer tubes are more collapsible than shorter tubes is consistent with the higher incidence of obstructive sleep apnea in males than females.

Advancements in Methods and Camera-Based Sensors for the Quantification of Respiration

Assessment of respiratory function allows early detection of potential disorders in the respiratory system and provides useful information for medical management. There is a wide range of applications for breathing assessment, from measurement systems in a clinical environment to applications involving athletes. Many studies on pulmonary function testing systems and breath monitoring have been conducted over the past few decades, and their results have the potential to broadly impact clinical practice. However, most of these works require physical contact with the patient to produce accurate and reliable measures of the respiratory function. There is still a significant shortcoming of non-contact measuring systems in their ability to fit into the clinical environment. The purpose of this paper is to provide a review of the current advances and systems in respiratory function assessment, particularly camera-based systems. A classification of the applicable research works is presented according to their techniques and recorded/quantified respiration parameters. In addition, the current solutions are discussed with regards to their direct applicability in different settings, such as clinical or home settings, highlighting their specific strengths and limitations in the different environments.

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.

Computational modelling of nasal respiratory flow

CFD has emerged as a promising diagnostic tool for clinical trials, with tremendous potential. However, for real clinical applications to be useful, overall statistical findings from large population samples (e.g., multiple cases and models) are needed. Fully resolved solutions are not a priority, but rather rapid solutions with fast turn-around times are desired. This leads to the issue of what are the minimum modelling criteria for achieving adequate accuracy in respiratory flows for large-scale clinical applications, with a view to rapid turnaround times. This study simulated a highly-resolved solution using the large eddy simulation (LES) method as a reference case for comparison with lower resolution models that included larger time steps and no turbulence modelling. Differences in solutions were quantified by pressure loss, flow resistance, unsteadiness, turbulence intensity, and hysteresis effects from multiple cycles. The results demonstrated that sufficient accuracy could be achieved with lower resolution models if the mean flow was considered. Furthermore, to achieve an established transient result unaffected by the initial start-up quiescent effects, the results need to be taken from at least the second respiration cycle. It was also found that the exhalation phase exhibited strong turbulence. The results are expected to provide guidance for future modelling efforts for clinical and engineering applications requiring large numbers of cases using simplified modelling approaches.

Computational fluid dynamics modelling of human upper airway: A review

BACKGROUND AND OBJECTIVE

Human upper airway (HUA) has been widely investigated by many researchers covering various aspects, such as the effects of geometrical parameters on the pressure, velocity and airflow characteristics. Clinically significant obstruction can develop anywhere throughout the upper airway, leading to asphyxia and death; this is where recognition and treatment are essential and lifesaving. The availability of advanced computer, either hardware or software, and rapid development in numerical method have encouraged researchers to simulate the airflow characteristics and properties of HUA by using various patient conditions at different ranges of geometry and operating conditions. Computational fluid dynamics (CFD) has emerged as an efficient alternative tool to understand the airflow of HUA and in preparing patients to undergo surgery. The main objective of this article is to review the literature that deals with the CFD approach and modeling in analyzing HUA.

METHODS

This review article discusses the experimental and computational methods in the study of HUA. The discussion includes computational fluid dynamics approach and steps involved in the modeling used to investigate the flow characteristics of HUA. From inception to May 2020, databases of PubMed, Embase, Scopus, the Cochrane Library, BioMed Central, and Web of Science have been utilized to conduct a thorough investigation of the literature. There had been no language restrictions in publication and study design of the database searches. A total of 117 articles relevant to the topic under investigation were thoroughly and critically reviewed to give a clear information about the subject. The article summarizes the review in the form of method of studying the HUA, CFD approach in HUA, and the application of CFD for predicting HUA obstacle, including the type of CFD commercial software are used in this research area.

RESULTS

This review found that the human upper airway was well studied through the application of computational fluid dynamics, which had considerably enhanced the understanding of flow in HUA. In addition, it assisted in making strategic and reasonable decision regarding the adoption of treatment methods in clinical settings. The literature suggests that most studies were related to HUA simulation that considerably focused on the aspects of fluid dynamics. However, there is a literature gap in obtaining information on the effects of fluid-structure interaction (FSI). The application of FSI in HUA is still limited in the literature; as such, this could be a potential area for future researchers. Furthermore, majority of researchers present the findings of their work through the mechanism of airflow, such as that of velocity, pressure, and shear stress. This includes the use of Navier-Stokes equation via CFD to help visualize the actual mechanism of the airflow. The above-mentioned technique expresses the turbulent kinetic energy (TKE) in its result to demonstrate the real mechanism of the airflow. Apart from that, key result such as wall shear stress (WSS) can be revealed via turbulent kinetic energy (TKE) and turbulent energy dissipation (TED), where it can be suggestive of wall injury and collapsibility tissue to the HUA.

The perspective of fluid flow behavior of respiratory droplets and aerosols through the facemasks in context of SARS-CoV-2

In the unfortunate event of the current ongoing pandemic COVID-19, where vaccination development is still in the trial phase, several preventive control measures such as social distancing, hand-hygiene, and personal protective equipment have been recommended by health professionals and organizations. Among them, the safe wearing of facemasks has played a vital role in reducing the likelihood and severity of infectious respiratory disease transmission. The reported research in facemasks has covered many of their material types, fabrication techniques, mechanism characterization, and application aspects. However, in more recent times, the focus has shifted toward the theoretical investigations of fluid flow mechanisms involved in the virus-laden particles' prevention by using facemasks. This exciting research domain aims to address the complex fluid transport that led to designing a facemask with a better performance. This Review discusses the recent updates on fluid flow dynamics through the facemasks. Key design aspects such as thermal comfort and flow resistance are discussed. Furthermore, the recent progress in the investigations on the efficacy of facemasks for the prevention of COVID-19 spread and the impact of wearing facemasks is presented.

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.

Imaging the breastfeeding swallow: Pilot study utilizing real-time MRI

OBJECTIVE

Knowledge of the breastfeeding swallow is limited by practical challenges. Radiation exposure to both mother and infant and the radiolucent properties of breastmilk make videofluoroscopy an unsuitable imaging modality. Furthermore, ultrasound is not ideal for capturing the complex 3-dimensional functional anatomy of swallowing. In this study we explore the feasibility of using real-time MRI to capture the breastfeeding swallow.

METHODS

Prospective observational study: Review of imaging from 12 normal infants (<5 months of age) and their mothers while breastfeeding using real-time MRI.

RESULTS

Static images were successfully captured in 11 infants and dynamic images in nine infants. This imaging modality confirms the dorsal surface of the infant's tongue elevates the maternal nipple to the hard palate, closing the space around the nipple with no air visible in the oral cavity during sucking and swallowing. We obtained dynamic imaging of mandibular movement with sucking, palatal elevation and pharyngeal constriction with swallowing, diaphragm movement with breathing and milk entering the stomach. Breastmilk was easily visualized, being high intensity on T2 sequences. Technical challenges were encountered secondary to infant movement and difficulties acquiring and maintaining midsagittal orientation. The similarity in tissue densities of the lips, tongue, nipple and hard palate limited definition between these structures.

CONCLUSION

Real-time MRI imaging was successful in capturing dynamic images of the breastfeeding swallow. However, technical and practical challenges make real-time MRI unlikely at present to be suitable for swallow assessment in clinical practice. Advances in technology and expertise in dynamic image capture may improve the feasibility of using MRI to understand and assess the breastfeeding swallow in the near future.

LEVEL OF EVIDENCE

4.

Targeted Drug Delivery of Magnetic Nano-Particle in the Specific Lung Region

Aerosolized drug inhalation plays an important role in the treatment of respiratory diseases. All of the published in silico, in vivo, and in vitro studies have improved the knowledge of aerosol delivery in the human respiratory system. However, aerosolized magnetic nano-particle (MNP) transport and deposition (TD) for the specific position of the human lung are still unavailable in the literature. Therefore, this study is aimed to provide an understanding of the magnetic nano-particle TD in the targeted region by imposing an external magnetic field for the development of future therapeutics. Uniform aerosolized nano-particle TD in the specific position of the lung airways will be modelled by adopting turbulence k–ω low Reynolds number simulation. The Euler–Lagrange (E–L) approach and the magneto hydrodynamics (MHD) model are incorporated in the ANSYS fluent (18.0) solver to investigate the targeted nano-particle TD. The human physical activity conditions of sleeping, resting, light activity and fast breathing are considered in this study. The aerosolized drug particles are navigated to the targeted position under the influence of external magnetic force (EMF), which is applied in two different positions of the two-generation lung airways. A numerical particle tracing model is also developed to predict the magnetic drug targeting behavior in the lung. The numerical results reveal that nano-particle deposition efficiency (DE) in two different magnetic field position is different for various physical activities, which could be helpful for targeted drug delivery to a specific region of the lung after extensive clinical trials. This process will also be cost-effective and will minimize unwanted side effects due to systemic drug distribution in the lung.

A Review of Respiratory Anatomical Development, Air Flow Characterization and Particle Deposition

The understanding of complex inhalation and transport processes of pollutant particles through the human respiratory system is important for investigations into dosimetry and respiratory health effects in various settings, such as environmental or occupational health. The studies over the last few decades for micro- and nanoparticle transport and deposition have advanced the understanding of drug-aerosol impacts in the mouth-throat and the upper airways. However, most of the Lagrangian and Eulerian studies have utilized the non-realistic symmetric anatomical model for airflow and particle deposition predictions. Recent improvements to visualization techniques using high-resolution computed tomography (CT) data and the resultant development of three dimensional (3-D) anatomical models support the realistic representation of lung geometry. Yet, the selection of different modelling approaches to analyze the transitional flow behavior and the use of different inlet and outlet conditions provide a dissimilar prediction of particle deposition in the human lung. Moreover, incorporation of relevant physical and appropriate boundary conditions are important factors to consider for the more accurate prediction of transitional flow and particle transport in human lung. This review critically appraises currently available literature on airflow and particle transport mechanism in the lungs, as well as numerical simulations with the aim to explore processes involved. Numerical studies found that both the Euler–Lagrange (E-L) and Euler–Euler methods do not influence nanoparticle (particle diameter ≤50 nm) deposition patterns at a flow rate ≤25 L/min. Furthermore, numerical studies demonstrated that turbulence dispersion does not significantly affect nanoparticle deposition patterns. This critical review aims to develop the field and increase the state-of-the-art in human lung modelling.

Nasal sprayed particle deposition in a human nasal cavity under different inhalation conditions

Deposition of polydisperse particles representing nasal spray application in a human nasal cavity was performed under transient breathing profiles of sniffing, constant flow, and breath hold. The LES turbulence model was used to describe the fluid phase. Particles were introduced into the flow field with initial spray conditions, including spray cone angle, insertion angle, and initial velocity. Since nasal spray atomizer design determines the particle conditions, fifteen particle size distributions were used, each defined by a log-normal distribution with a different volume mean diameter (Dv50). Particle deposition in the anterior region was approximately 80% when Dv50 > 50μm, and this decreased to 45% as Dv50 decreased to 10μ m for constant and sniff breathing conditions. The decrease in anterior deposition was countered with increased deposition in the middle and posterior regions. The significance of increased deposition in the middle region for drug delivery shows there is potential for nasal delivered drugs to reach the highly vascularised mucosal walls in the main nasal passages. For multiple targeted deposition sites, an optimisation equation was introduced where deposition results of any two targeted sites could be combined and a weighting between 0 to 1 was applied to each targeted site, representing the relative importance of each deposition site.

Quantitative Assessment of Regional Dynamic Airway Collapse in Neonates via Retrospectively Respiratory-Gated 1 H Ultrashort Echo Time MRI

BACKGROUND

Neonatal dynamic tracheal collapse (tracheomalacia, TM) is a common and serious comorbidity in infants, particularly those with chronic lung disease of prematurity (bronchopulmonary dysplasia, BPD) or congenital airway or lung-related conditions such as congenital diaphragmatic hernia (CDH), but the underlying pathology, impact on clinical outcomes, and response to therapy are not well understood. There is a pressing clinical need for an accurate, objective, and safe assessment of neonatal TM.

PURPOSE

To use retrospectively respiratory-gated ultrashort echo-time (UTE) MRI to noninvasively analyze moving tracheal anatomy for regional, quantitative evaluation of dynamic airway collapse in quiet-breathing, nonsedated neonates.

STUDY TYPE

Prospective.

POPULATION/SUBJECTS

Twenty-seven neonatal subjects with varying respiratory morbidities (control, BPD, CDH, abnormal polysomnogram).

FIELD STRENGTH/SEQUENCE

High-resolution 3D radial UTE MRI (0.7 mm isotropic) on 1.5T scanner sited in the neonatal intensive care unit.

ASSESSMENT

Images were retrospectively respiratory-gated using the motion-modulated time-course of the k-space center. Tracheal surfaces were generated from segmentations of end-expiration/inspiration images and analyzed geometrically along the tracheal length to calculate percent-change in luminal cross-sectional area (A % ) and ratio of minor-to-major diameters at end-expiration (r D,exp ). Geometric results were compared to clinically available bronchoscopic findings (n = 14).

STATISTICAL TESTS

Two-sample t-test.

RESULTS

Maximum A % significantly identified subjects with/without a bronchoscopic TM diagnosis (with: 46.9 ± 10.0%; without: 27.0 ± 5.8%; P < 0.001), as did minimum r D,exp (with: 0.346 ± 0.146; without: 0.671 ± 0.218; P = 0.008). Subjects with severe BPD exhibited a far larger range of minimum r D,exp than subjects with mild/moderate BPD or controls (0.631 ± 0.222, 0.782 ± 0.075, and 0.776 ± 0.030, respectively), while minimum r D,exp was reduced in CDH subjects (0.331 ± 0.171) compared with controls (P < 0.001).

DATA CONCLUSION

Respiratory-gated UTE MRI can quantitatively and safely evaluate neonatal dynamic tracheal collapse, as validated with the clinical standard of bronchoscopy, without requiring invasive procedures, anesthesia, or ionizing radiation.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2019;49:659-667.

Effect of upper airway on tracheobronchial fluid dynamics

The upper airways play a significant role in the tracheal flow dynamics. Despite many previous studies, however, the effect of the upper airways on the ventilation distribution in distal airways has remained a challenge. The aim of this study is to experimentally and computationally investigate the dynamic behaviour in the intratracheal flow induced by the upper respiratory tract and to assess its influence on the subsequent tributaries. Patient-specific images from 2 different modalities (magnetic resonance imaging of the upper airways and computed tomography of the lower airways) were segmented and combined. An experimental phantom of patient-specific airways (including the oral cavity, larynx, trachea, down to generations 6-8) was generated using 3D printing. The flow velocities in this phantom model were measured by the flow-sensitised phase contrast magnetic resonance imaging technique and compared with the computational fluid dynamics simulations. Both experimental and computational results show a good agreement in the time-averaged velocity fields as well as fluctuating velocity. The flows in the proximal trachea were complex and unsteady under both lower- and higher-flow rate conditions. Computational fluid dynamics simulations were also performed with an airways model without the upper airways. Although the flow near the carina remained unstable only when the inflow rate was high, the influence of the upper airways caused notable changes in distal flow distributions when the 2 airways models were compared with and without the upper airways. The results suggest that the influence of the upper airways should be included in the respiratory flow assessment as the upper airways extensively affect the flows in distal airways and consequent ventilation distribution in the lungs.

Magnetic resonance imaging of obstructive sleep apnea in children

Sleep-disordered breathing has a spectrum of severity that spans from snoring and partial airway collapse with increased upper airway resistance, to complete upper airway obstruction with obstructive sleep apnea during sleeping. While snoring occurs in up to 20% of children, obstructive sleep apnea affects approximately 1-5% of children. The obstruction that occurs in obstructive sleep apnea is the result of the airway collapsing during sleep, which causes arousal and impairs restful sleep. Adenotonsillectomy is the first-line treatment of obstructive sleep apnea and is usually effective in otherwise healthy nonsyndromic children. However, there are subgroups in which this surgery is less effective. These subgroups include children with obesity, severe obstructive sleep apnea preoperatively, Down syndrome, craniofacial anomalies and polycystic ovarian disease. Continuous positive airway pressure (CPAP) is the first-line therapy for persistent obstructive sleep apnea despite previous adenotonsillectomy, but it is often poorly tolerated by children. When CPAP is not tolerated or preferred by the family, surgical options beyond adenotonsillectomy are discussed with the parent and child. Dynamic MRI of the airway provides a means to identify and localize the site or sites of obstruction for these children. In this review the authors address clinical indications for imaging, ideal team members to involve in an effective multidisciplinary program, basic anesthesia requirements, MRI protocol techniques and interpretation of the findings on MRI that help guide surgery.

Molecular Binding Contributes to Concentration Dependent Acrolein Deposition in Rat Upper Airways: CFD and Molecular Dynamics Analyses

Existing in vivo experiments show significantly decreased acrolein uptake in rats with increasing inhaled acrolein concentrations. Considering that high-polarity chemicals are prone to bond with each other, it is hypothesized that molecular binding between acrolein and water will contribute to the experimentally observed deposition decrease by decreasing the effective diffusivity. The objective of this study is to quantify the probability of molecular binding for acrolein, as well as its effects on acrolein deposition, using multiscale simulations. An image-based rat airway geometry was used to predict the transport and deposition of acrolein using the chemical species model. The low Reynolds number turbulence model was used to simulate the airflows. Molecular dynamic (MD) simulations were used to study the molecular binding of acrolein in different media and at different acrolein concentrations. MD results show that significant molecular binding can happen between acrolein and water molecules in human and rat airways. With 72 acrolein embedded in 800 water molecules, about 48% of acrolein compounds contain one hydrogen bond and 10% contain two hydrogen bonds, which agreed favorably with previous MD results. The percentage of hydrogen-bonded acrolein compounds is higher at higher acrolein concentrations or in a medium with higher polarity. Computational dosimetry results show that the size increase caused by the molecular binding reduces the effective diffusivity of acrolein and lowers the chemical deposition onto the airway surfaces. This result is consistent with the experimentally observed deposition decrease at higher concentrations. However, this size increase can only explain part of the concentration-dependent variation of the acrolein uptake and acts as a concurrent mechanism with the uptake-limiting tissue ration rate. Intermolecular interactions and associated variation in diffusivity should be considered in future dosimetry modeling of high-polarity chemicals such as acrolein.