Computational modeling and validation of human nasal airflow under various breathing conditions

Experimental and numerical investigation of pressure distribution in human upper airway passage before and after maxillary sinus surgery

Sinusitis, a common disease of the maxillary sinus, is initially managed with saline solution and medication, resulting in the resolution of symptoms within a few days in most cases. However, Functional Endoscopic Sinus Surgeries are recommended if pharmacological treatments prove ineffective. This research aims to investigate the effects of maxillary sinus surgery on the airflow field, pressure distribution within the nasal cavity, and overall ventilation. This study utilized a three-dimensional realistic nasal cavity model constructed from CT images of a healthy adult. Virtual surgery including uncinectomy with Middle Meatal Antrostomy, two standard procedures performed during such surgeries, was performed on the model under the supervision of a clinical specialist. Two replicas representing pre- and post-operative cases were created using 3D printing for experimental purposes. Various breathing rates ranging from 3.8 to 42.6 L/min were examined through experimental and numerical simulations. To ensure the accuracy of the numerical simulations, the results were compared to measured pressure data, showing a reasonable agreement between the two. The findings demonstrate that uncinectomy and Middle Meatal Antrostomy significantly enhance the ventilation of the maxillary sinuses. Furthermore, increasing inspiratory rates leads to further improvements in ventilation. The static pressure distribution within the maxillary sinuses remains relatively uniform, except in regions close to the sinus ostium, even after surgical intervention.

Is the mouse nose a miniature version of a rat nose? A computational comparative study

BACKGROUND AND OBJECTIVE

Although the mouse is a widely used animal model in biomedical research, there are few published studies on its nasal aerodynamics, potentially due to its small size. It is not appropriate to assume that mice and rats' nasal structure and airflow characteristics are the same because the ratio of nasal surface area to nasal volume and body weight is much higher in a mouse than in a rat. The aim of this work is to use anatomically accurate image-based computational fluid dynamic modeling to quantitatively reveal the characteristics of mouse nasal airflow and mass transport that haven't been detailed before and find key differences to that of rat nose, which will deepen our understanding of the mouse's physiological functions.

METHODS

We created an anatomically accurate 3D computational nasal model of a B6 mouse using postmortem high-resolution micro-CT scans and simulated the airflow distribution and odor transport patterns under restful breathing conditions. The deposition pattern of airborne particles was also simulated and validated against experimental data. In addition, we calculated the gas chromatograph efficiency of odor transport in the mouse employing the theoretical plate concept and compared it with previous studies involving cat and rat models.

RESULTS

Similar to the published rat model, respiratory and olfactory flow regimes are clearly separated in the mouse nasal cavity. A high-speed dorsal medial (DM) stream was observed, which enhances the delivery speed and efficiency of odor to the ethmoid (olfactory) recess (ER). The DM stream split into axial and secondary paths in the ER. However, the secondary flow in the mouse is less extensive than in the rat. The gas chromatograph efficiency calculations suggest that the rat may possess a moderately higher odorant transport efficiency than that of the mouse due to its more complex ethmoid recess structure and extensive secondary flow. However, the mouse's nasal structure seems to adapt better to varying airflow velocity.

CONCLUSIONS

Due to the inherent structural disparities, the rat and mouse models exhibit moderate differences in airflow and mass transport patterns, potentially impacting their olfaction and other behavioral habits.

Oxymetazoline as a predictor of turbinate reduction surgery outcomes: Objective support from a prospective, single-blinded, computational fluid dynamics study

BACKGROUND

A patient's subjective response to topical nasal decongestant is often used to screen for turbinate reduction surgery suitability. However, this anecdotal strategy has not been objectively and quantitatively evaluated.

METHODS

Prospective, longitudinal, and single-blinded cohort study employing computational fluid dynamic modeling based on computed tomography scans at baseline, 30 min postoxymetazoline, and 2 months postsurgery on 11 patients with chronic turbinate hypertrophy.

RESULTS

Nasal obstruction symptom evaluation (NOSE) and visual analogue scale (VAS) obstruction scores significantly improved from baseline to postoxymetazoline and again to postsurgery (NOSE: 71.82 ± 14.19 to 42.27 ± 25.26 to 22.27 ± 21.04; VAS: 6.09 ± 2.41 to 4.14 ± 2.20 to 2.08 ± 1.56; each interaction p < 0.05), with significant correlation between the latter two states (r∼0.37-0.69, p < 0.05). Oxymetazoline had a broader anatomical impact throughout inferior and middle turbinates than surgery (many p < 0.05); however, the improvement in regional airflow is similar (most p > 0.05) and predominantly surrounding the inferior turbinate. Strong postoxymetazoline to postsurgery correlations were observed in decreased nasal resistance (r = 0.79, p < 0.05), increased regional airflow rates (r = -0.47 to -0.55, p < 0.05) and regional air/mucosa shear force and heat flux (r = 0.43 to 0.58, p < 0.05); however, only increasing peak heat flux significantly correlated to symptom score improvement (NOSE: r = 0.48, p < 0.05).

CONCLUSION

We present the first objective evidence that the "topical decongestant test" can help predict turbinate reduction surgery outcomes. The predictive effect is driven by similar improvementin regional airflow that leading to improved air/mucosa stimulations (peak heat flux) rather than through reduced nasal resistance.

Computational analysis of nasal airflow and its alteration by a nasal dilator

Nasal airflow obstruction correlates with several ailments, such as higher patency, increased friction at the mucosal wall or the so-called Little's area, improper air conditioning, and snoring. Nasal dilators are frequently employed, mainly due to their ease of access and use, combined with their non-permanent and non-surgical nature. Their overall efficacy, however, has not been clearly demonstrated so far, with some studies reporting conflicting outcomes, mainly because being based on subjective evaluations. This study employs Computational Fluid Dynamics simulations to analyze the flow inside a real nose, performs an objective assessment of a nasal dilator's effect in terms of airflow and air conditioning, reporting flow paths, friction levels, heat and water fluxes and detailed temperature and humidity distributions. Coincidentally, the studied nose presents a septal deviation, with one nostril being wider than the other. The tubes of the dilator used in both nostrils are identical, as with any standard commercial dilator. Consequently, the dilator widens one nostril, as intended, but results in an obstruction in the other. This allows simultaneously addressing two situations, the nominal function of the dilator, as well as an off-design case. Results indicate a 24 % increase in nasal patency in the design situation. The effect, however, is limited, as quantified by appropriate measures, such as the flow-generated friction at the nose surfaces and the temperature fluxes. Hence, the effect of such a dilator in nominal conditions is perhaps not as large as might be hoped. In the off-design situation, nasal resistance increases by 62 %, an undesirable effect, illustrating the consequences of using an inappropriate dilator.

Flow Analysis of Central Venous Outflow Tract: A New Approach to Understanding Pulse-Synchronous Tinnitus

OBJECTIVE

Pulse-synchronous tinnitus (PST) has been linked to multiple anatomical variants of the central venous outflow tract (CVOT) including sigmoid sinus (SS) dehiscence and diverticulum. This study investigates flow turbulence, pressure, and wall shear stress along the CVOT and proposes a mechanism that results in SS dehiscence and PST.

STUDY DESIGN

Case series.

SETTING

Tertiary Academic Center.

METHODS

Venous models were reconstructed from computed tomography scans of 3 patients with unilateral PST. Two models for each patient are obtained: a symptomatic and contralateral asymptomatic side. A turbulent model-enabled commercial flow solver was used to simulate the pulsatile blood flow over the cardiac cycle through the models. Fluid flow through the transverse and SS junction was analyzed to observe the velocity, pressure, turbulent kinetic energy (TKE), and shear stress over a simulated cardiac cycle.

RESULTS

Fluid flow on the symptomatic side showed increased vorticity in the presence of an SS diverticulum. Higher TKE with periodicity following the cardiac cycle was observed on the symptomatic side, and a sharp increase was observed if SS diverticulum was present. Shear stress was highest near the narrowest segments of the vessel. Pressure was observed to be lower on the symptomatic side at the transverse-SS junction for all 3 patients.

CONCLUSION

Computational fluid dynamics modeling of blood flow through the CVOT in PST suggests that low pressure may be the cause of dehiscence, and tinnitus may result from periodic increases in TKE.

An adjoint-based approach for the surgical correction of nasal septal deviations

Deviations of the septal wall are widespread anatomic anomalies of the human nose; they vary significantly in shape and location, and often cause the obstruction of the nasal airways. When severe, septal deviations need to be surgically corrected by ear-nose-throat (ENT) specialists. Septoplasty, however, has a low success rate, owing to the lack of suitable standardized clinical tools for assessing type and severity of obstructions, and for surgery planning. Moreover, the restoration of a perfectly straight septal wall is often impossible and possibly unnecessary. This paper introduces a procedure, based on advanced patient-specific Computational Fluid Dynamics (CFD) simulations, to support ENT surgeons in septoplasty planning. The method hinges upon the theory of adjoint-based optimization, and minimizes a cost function that indirectly accounts for viscous losses. A sensitivity map is computed on the mucosal wall to provide the surgeon with a simple quantification of how much tissue removal at each location would contribute to easing the obstruction. The optimization procedure is applied to three representative nasal anatomies, reconstructed from CT scans of patients affected by complex septal deviations. The computed sensitivity consistently identifies all the anomalies correctly. Virtual surgery, i.e. morphing of the anatomies according to the computed sensitivity, confirms that the characteristics of the nasal airflow improve significantly after small anatomy changes derived from adjoint-based optimization.

Computational Fluid Dynamics and Its Potential Applications for the ENT Clinician

This article is an examination of computational fluid dynamics in the field of otolaryngology, specifically rhinology. The historical development and subsequent application of computational fluid dynamics continues to enhance our understanding of various sinonasal conditions and surgical planning in the field today. This article aims to provide a description of computational fluid dynamics, the methods for its application, and the clinical relevance of its results. Consideration of recent research and data in computational fluid dynamics demonstrates its use in nonhistological disease pathology exploration, accompanied by a large potential for surgical guidance applications. Additionally, this article defines in lay terms the variables analyzed in the computational fluid dynamic process, including velocity, wall shear stress, area, resistance, and heat flux.

The air conditioning in the nose of mammals depends on their mass and on their maximal running speed

The nose of the mammals is responsible for filtering, humidifying, and heating the air before entering the lower respiratory tract. This conditioning avoids, notably, dehydration of the bronchial and alveolar mucosa. However, since this conditioning is not perfect, exercising in cold air can induce lung inflammation, both for human and non-human mammals. This work aims to compare the air conditioning in the noses of various mammals during inspiration. We build our study on computational fluid dynamics simulations of the heat exchanges in the lumen of the upper respiratory tract of these mammals. These simulations show that the efficiency of the air conditioning in the nose during inspiration does not relate only to the mass m of the mammal but also to its maximal running speed v. More precisely, the results allow establishing a scaling law relating the efficiency of air conditioning in the nose of mammals to the ratio v / log 10 ( m ) . The simulations also correlate the resistance to the flow in the nose to the efficiency of this air conditioning. The obtained scaling law allows predicting the air temperature at the top of the trachea during inspiration for nasal-breathing mammals, and thus notably for humans of various ages.

Aerodynamic evaluation of surgical design for the stenosis correction of airway

Introduction

Congenital tracheal stenosis (CTS) is a rare but life-threatening disease that can lead to respiratory dysfunction in children. Obstructive sleep apnea syndrome (OSAS) in children is characterized by prolonged partial upper airway obstruction and/or intermittent complete obstruction. Both of the diseases require surgical intervention. Although respective treatments of these two diseases are clear, there is a lack of literature discussing the surgical treatment of patients with CTS complicated by OSAS.

Methods

We conducted a patient-specific study of patient with CTS complicated by OSAS. Computer-aided design was used to simulate surgical correction under different surgical sequences. Computational fluid dynamics was used to compare the outcomes of different sequences.

Results

Aerodynamic parameters, pressure drop, velocity streamlines, wall shear stress (WSS), and the ratio of airflow distribution and energy loss rate were evaluated. An obvious interaction was found between the two diseases in different surgical sequences. The order of correction for CTS or OSAS greatly affected the aerodynamic parameters and turbulence flows downstream of tracheal stenosis and upstream of epiglottis. The CTS and OSAS had mutual influences on each other on the aerodynamic parameters, such as pressure drops and WSS.

Discussion

When evaluating the priority of surgical urgency of CTS and OSAS, surgeons need to pay attention to the state of both CTS and OSAS and the physiological conditions of patients. The aerodynamic performance of the uneven airflow distribution and the potential impact caused by the correction of CTS should be considered in surgical planning and clinical management.

How Does Oxymetazoline Change Nasal Aerodynamics and Symptomatology in Patients with Turbinate Hypertrophy?

OBJECTIVES

Oxymetazoline relieves nasal obstructive symptoms via vasoconstriction, however, the changes in nasal structures and aerodynamics that impact symptoms the most remain unclear.

METHODS

This prospective, longitudinal, and single blinded cohort study applied Computational Fluid Dynamic (CFD) modeling based on CT scans at baseline and post-oxymetazoline on 13 consecutive patients with chronic nasal obstruction secondary to inferior turbinate hypertrophy from a tertiary medical center. To account for placebo effect, a sham saline spray was administered with subject blindfolded prior to oxymetazoline, with 30 min rest in between. Nasal Obstruction Symptom Evaluation (NOSE) and unilateral Visual Analogue Scale (VAS) scores of nasal obstructions were collected at baseline, after sham, and 30 min after oxymetazoline.

RESULTS

Both VAS and NOSE scores significantly improved from baseline to post-oxymetazoline (NOSE: 62.3 ± 12.4 to 31.5 ± 22.5, p < 0.01; VAS: 5.27 ± 2.63 to 3.85 ± 2.59, p < 0.05), but not significantly from baseline to post-sham. The anatomical effects of oxymetazoline were observed broadly throughout the entire length of the inferior and middle turbinates (p < 0.05). Among many variables that changed significantly post-oxymetazoline, only decreased nasal resistance (spearman r = 0.4, p < 0.05), increased regional flow rates (r = -0.3 to -0.5, p < 0.05) and mucosal cooling heat flux (r = -0.42, p < 0.01) in the inferior but not middle turbinate regions, and nasal valve Wall Shear Stress (WSS r = -0.43, p < 0.05) strongly correlated with symptom improvement.

CONCLUSION

Oxymetazoline broadly affects the inferior and middle turbinates, however, symptomatic improvement appears to be driven more by global nasal resistance and regional increases in airflow rate, mucosal cooling, and WSS, especially near the head of the inferior turbinate.

LEVEL OF EVIDENCE

3: Well-designed, prospective, single blinded cohort trial. Laryngoscope, 134:1100-1106, 2024.

Computational Rhinology: Unraveling Discrepancies between In Silico and In Vivo Nasal Airflow Assessments for Enhanced Clinical Decision Support

Computational rhinology is a specialized branch of biomechanics leveraging engineering techniques for mathematical modelling and simulation to complement the medical field of rhinology. Computational rhinology has already contributed significantly to advancing our understanding of the nasal function, including airflow patterns, mucosal cooling, particle deposition, and drug delivery, and is foreseen as a crucial element in, e.g., the development of virtual surgery as a clinical, patient-specific decision support tool. The current paper delves into the field of computational rhinology from a nasal airflow perspective, highlighting the use of computational fluid dynamics to enhance diagnostics and treatment of breathing disorders. This paper consists of three distinct parts-an introduction to and review of the field of computational rhinology, a review of the published literature on in vitro and in silico studies of nasal airflow, and the presentation and analysis of previously unpublished high-fidelity CFD simulation data of in silico rhinomanometry. While the two first parts of this paper summarize the current status and challenges in the application of computational tools in rhinology, the last part addresses the gross disagreement commonly observed when comparing in silico and in vivo rhinomanometry results. It is concluded that this discrepancy cannot readily be explained by CFD model deficiencies caused by poor choice of turbulence model, insufficient spatial or temporal resolution, or neglecting transient effects. Hence, alternative explanations such as nasal cavity compliance or drag effects due to nasal hair should be investigated.

Numerical and experimental analysis of pollen inhalation exposure in nasal airways following various middle turbinectomy

BACKGROUND

Middle turbinectomy (MT) has always been controversial. MT significantly alters the anatomy and redistributes the inhaled air. The current study is designed to quantify the effect of MT with varying resection volumes on airflow and associated pollen inhalation exposure characteristics in the nasal airways.

METHODS

Six realistic models following bilateral comprehensive Functional Endoscopic Sinus Surgery (FESS) deriving from CT images were constructed and their corresponding post-MT models with four types of MT procedures were virtually conducted. Inhalation exposure to pollen particles was simulated by the Computational Fluid-Particle Dynamics (CFPD) approach and validated through in vitro experiments.

RESULTS

Following the excision of the middle turbinate, a significant escalation in airflow was observed within the upper-middle region of the nasal cavities. Pollen deposition was observed to be more prominent in the nasal septum, laryngopharynx, and maxillary sinus, varying with the types of MT procedures. Notably, particles with diameters smaller than 50 μm exhibited two distinct "high peaks" and three "small peaks" within the nasal airways.

CONCLUSION

MT resulted in increased airflow volume within the upper-middle region of the nasal cavities. Following MT, notable shifts in pollen particle deposition hot spots were observed, transitioning from the nasal vestibule, nasal septum, and middle meatus to the nasal septum and laryngopharynx. These findings are anticipated to contribute valuable perspectives on pollen inhalation exposure risk assessments following diverse MT surgical interventions.

Numerical simulation of upper airway heat transfer in children with mandibular retrognathia during inspiratory process

Introduction

The human upper airway regulates temperature, but its heating capacity remains unclear when the ambient temperature is low and the airway structure is abnormal. Therefore, the purpose of this study was to evaluate the heat transfer characteristics of the upper airway in patients with mandibular retrognathia, and to quantitatively evaluate the influence of ambient temperature on the temperature field of the upper airway, which could provide a valuable reference for the prediction, diagnosis and treatment of respiratory tract related diseases.

Methods

Two typical ambient temperatures of -10 °C and 20 °C were selected to numerically simulate the air flow and heat transfer synchronization in the upper airway model of mandibular retrognathia under quiet breathing and heavy breathing.

Results and discussion

The results showed that the inspired air could not be sufficiently heated after flowing through the upper airway and main trachea in the two breathing states under low temperature conditions, and the inferior bronchus was more stimulated under the state of heavy breathing. In addition, the vortex flow structure in the upper airway can enhance the convective heat transfer effect in the corresponding area.

Toward a Mesoscopic Modeling Approach of Magnetohydrodynamic Blood Flow in Pathological Vessels: A Comprehensive Review

The investigation of magnetohydrodynamic (MHD) blood flow within configurations that are pertinent to the human anatomy holds significant importance in the realm of scientific inquiry because of its practical implications within the medical field. This article presents an exhaustive appraisal of the diverse applications of magnetohydrodynamics and their computational modeling in biological contexts. These applications are classified into two categories: simple flow and pulsatile flow. An alternative approach of traditional CFD methods called Lattice Boltzmann Method (LBM), a mesoscopic method based on kinetic theory, is introduced to solve complex problems, such as hemodynamics. The results show that the flow velocity reduces considerably by increasing the magnetic field intensity, and the flow separation area is minimized by the increase of magnetic field strength. The LBM with BGK collision model has shown good results in terms of precision. Finally, this literature review has revealed a number of potential avenues for further research. Suggestions for future works are proposed accordingly.

Effect of ambient temperature and respiration rate on nasal dominance: preliminary findings from a nostril-specific wearable

The nasal dominance (ND) determination is crucial for nasal synchronized ventilator, optimum nasal drug delivery, identifying brain hemispheric dominance, nasal airway obstruction surgery, mindfulness breathing, and for possible markers of a conscious state. Given these wider applications of ND, it is interesting to understand the patterns of ND with varying temperature and respiration rates. In this paper, we propose a method which measures peak-to-peak temperature oscillations (difference between end-expiratory and end-inspiratory temperature) for the left and right nostrils during nasal breathing. These nostril-specific temperature oscillations are further used to calculate the nasal dominance index (NDI), nasal laterality ratio (NLR), inter-nostril correlation, and mean of peak-to-peak temperature oscillation for inspiratory and expiratory phase at (1) different ambient temperatures of 18 °C, 28 °C, and 38 °C and (2) at three different respiration rate of 6 bpm, 12 bpm, and 18 bpm. The peak-to-peak temperature (Tpp) oscillation range (averaged across participants;n= 8) for the left and right nostril were 3.80 ± 0.57 °C and 2.34 ± 0.61 °C, 2.03 ± 0.20 °C and 1.40 ± 0.26 °C, and 0.20 ± 0.02 °C and 0.29 ± 0.03 °C at the ambient temperature of 18 °C, 28 °C, and 38 °C respectively (averaged across participants and respiration rates). The NDI and NLR averaged across participants and three different respiration rates were 35.67 ± 5.53 and 2.03 ± 1.12; 8.36 ± 10.61 and 2.49 ± 3.69; and -25.04 ± 14.50 and 0.82 ± 0.54 at the ambient temperature of 18 °C, 28 °C, and 38 °C respectively. The Shapiro-Wilk test, and non-parametric Friedman test showed a significant effect of ambient temperature conditions on both NDI and NLR. No significant effect of respiration rate condition was observed on both NDI and NLR. The findings of the proposed study indicate the importance of ambient temperature while determining ND during the diagnosis of breathing disorders such as septum deviation, nasal polyps, nosebleeds, rhinitis, and nasal fractions, and in the intensive care unit for nasal synchronized ventilator.

Domestic cat nose functions as a highly efficient coiled parallel gas chromatograph

The peripheral structures of mammalian sensory organs often serve to support their functionality, such as alignment of hair cells to the mechanical properties of the inner ear. Here, we examined the structure-function relationship for mammalian olfaction by creating an anatomically accurate computational nasal model for the domestic cat (Felis catus) based on high resolution microCT and sequential histological sections. Our results showed a distinct separation of respiratory and olfactory flow regimes, featuring a high-speed dorsal medial stream that increases odor delivery speed and efficiency to the ethmoid olfactory region without compromising the filtration and conditioning purpose of the nose. These results corroborated previous findings in other mammalian species, which implicates a common theme to deal with the physical size limitation of the head that confines the nasal airway from increasing in length infinitely as a straight tube. We thus hypothesized that these ethmoid olfactory channels function as parallel coiled chromatograph channels, and further showed that the theoretical plate number, a widely-used indicator of gas chromatograph efficiency, is more than 100 times higher in the cat nose than an "amphibian-like" straight channel fitting the similar skull space, at restful breathing state. The parallel feature also reduces airflow speed within each coil, which is critical to achieve the high plate number, while feeding collectively from the high-speed dorsal medial stream so that total odor sampling speed is not sacrificed. The occurrence of ethmoid turbinates is an important step in the evolution of mammalian species that correlates to their expansive olfactory function and brain development. Our findings reveal novel mechanisms on how such structure may facilitate better olfactory performance, furthering our understanding of the successful adaptation of mammalian species, including F. catus, a popular pet, to diverse environments.

Changes in nasolabial angle may alter nasal valve morphology and airflow: a computational fluid dynamics study

Aim: Nasal valve (NV) dysfunctions are a significant cause of nasal obstruction. Changes in the nasolabial angle (NLA) may also cause changes in NV morphology. The effect of changes in the 3D structure of the nasal valve region (NVR) on nasal airflow has yet to be studied sufficiently. The accuracy of computational fluid dynamics (CFD) simulation results of nasal airflow has been confirmed by in vitro tests. Therefore, this study aimed to evaluate the effect of changes in NV structure and volume on nasal airflow based on the CFD method. Material and Method: We used CT images to create a 3D structural model of the NVR. First, CT images were transferred to MIMICS® software, and the nasal air passage was modeled. A solid reference model of the NVR was then created using SolidWorks software. Five different solid 3D nasal valve models were created with nasolabial angles of 85˚ in Model 1, 90˚ in Model 2, 95˚ in Model 3, 100˚ in Model 4, and 105˚ in Model 5. To simulate breathing during rest and exercise using the CFD method, the unilateral nasal airflow rates were set at 150 ml/s and 500 ml/s, respectively. The CFD method was then used to calculate each model’s airflow properties. Finally, the volumes of the models, pressure at the NV outlet, and airflow velocity were evaluated and calculated to investigate each model’s NV airflow characteristics. Results: Our study found a significant correlation between the nasolabial angle (NLA) and NVR volume (r=-0.998, p=0.000), flow rate and velocity (r=0.984, p=0.000), velocity and maximum pressure (r=0.920, p=0.000), velocity and minimum pressure (r=-0.969, p=0.000), flow rate and maximum pressure (r=0.974, p=0.000), and flow rate and minimum pressure (r=-0.950, p=0.000). There was no correlation between NLA increase and nasal airflow velocity. We determined that the highest pressure and lowest airflow velocity values were in the upper angle region and that the lowest pressure and highest airflow velocity values were at the bottom of the NVR in all models. Conclusion: Using the CFD method, we found a decrease in NVR volume and an increase in airflow velocity with an increase in NLA. In addition, we found that the pressure values in the NVR did not change significantly with the increase in NLA.

Numerical and Experimental Analysis of Drug Inhalation in Realistic Human Upper Airway Model

The demand for a more efficient and targeted method for intranasal drug delivery has led to sophisticated device design, delivery methods, and aerosol properties. Due to the complex nasal geometry and measurement limitations, numerical modeling is an appropriate approach to simulate the airflow, aerosol dispersion, and deposition for the initial assessment of novel methodologies for better drug delivery. In this study, a CT-based, 3D-printed model of a realistic nasal airway was reconstructed, and airflow pressure, velocity, turbulent kinetic energy (TKE), and aerosol deposition patterns were simultaneously investigated. Different inhalation flowrates (5, 10, 15, 30, and 45 L/min) and aerosol sizes (1, 1.5, 2.5, 3, 6, 15, and 30 µm) were simulated using laminar and SST viscous models, with the results compared and verified by experimental data. The results revealed that from the vestibule to the nasopharynx, the pressure drop was negligible for flow rates of 5, 10, and 15 L/min, while for flow rates of 30 and 40 L/min, a considerable pressure drop was observed by approximately 14 and 10%, respectively. However, from the nasopharynx and trachea, this reduction was approximately 70%. The aerosol deposition fraction alongside the nasal cavities and upper airway showed a significant difference in pattern, dependent on particle size. More than 90% of the initiated particles were deposited in the anterior region, while just under 20% of the injected ultrafine particles were deposited in this area. The turbulent and laminar models showed slightly different values for the deposition fraction and efficiency of drug delivery for ultrafine particles (about 5%); however, the deposition pattern for ultrafine particles was very different.

New approach to establish a surgical planning in infantile vallecular cyst synchronous with laryngomalacia based on aerodynamic analysis

BACKGROUND AND OBJECTIVES

A large proportion of infants with vallecular cyst (VC) have coexisting laryngomalacia (LM). Feeding difficulties, regurgitation, occasional cough, and sleep-disordered breathing are the common symptoms in moderate to severe cases. The surgical management of these cases is more challenging and remains controversial. The purpose of this study is to help surgeons select the effective surgical strategies by computer-aided design (CAD) and computational fluid dynamics (CFD) simulations of the upper airway flow characteristics.

METHODS

The three dimensional (3D) geometric model of the upper airway was reconstructed based on two dimensional (2D) medical images of the patient with VC accompanied with LM. Virtual surgeries were carried out preoperatively to simulate three possible post-operative states in silico. The different outcomes of virtual surgical strategies were predicted based on computational evaluations of airway fluid dynamics including pressure, resistance, velocity, and wall shear stress (WSS).

RESULTS

The CFD results of this study suggested the importance of the angle between the rim of epiglottis and arytenoid epiglottic (AE) fold. There was a small impact on the upper airway flow field while the VC was removed and the angle of epiglottis was unchanged. The partial lifting of epiglottis can further improve the flow field. With performing supraglottoplasty (SGP) and the marsupialization of VC, epiglottis was completely recovered, and the flow field was significantly improved. The clinical symptoms of this patient improved greatly after surgeries and no recurrence or growth retardation were noted during 1-year follow-up. The clinical prognosis was consistent with the prediction of the CFD results.

CONCLUSIONS

The state of epiglottis needs to be carefully checked to evaluate the necessity of performing further SGP in the patients with VC accompanied with LM. CFD and CAD could be developed as a new approach to help surgeons predict the post-operative outcomes through quantification of the airflow dynamics, and make the optimal and individualized surgical approaches for patients with airway obstruction.

In Silico Quantification of Intersubject Variability on Aerosol Deposition in the Oral Airway

The extrathoracic oral airway is not only a major mechanical barrier for pharmaceutical aerosols to reach the lung but also a major source of variability in lung deposition. Using computational fluid dynamics, deposition of 1−30 µm particles was predicted in 11 CT-based models of the oral airways of adults. Simulations were performed for mouth breathing during both inspiration and expiration at two steady-state flow rates representative of resting/nebulizer use (18 L/min) and of dry powder inhaler (DPI) use (45 L/min). Consistent with previous in vitro studies, there was a large intersubject variability in oral deposition. For an optimal size distribution of 1−5 µm for pharmaceutical aerosols, our data suggest that >75% of the inhaled aerosol is delivered to the intrathoracic lungs in most subjects when using a nebulizer but only in about half the subjects when using a DPI. There was no significant difference in oral deposition efficiency between inspiration and expiration, unlike subregional deposition, which shows significantly different patterns between the two breathing phases. These results highlight the need for incorporating a morphological variation of the upper airway in predictive models of aerosol deposition for accurate predictions of particle dosimetry in the intrathoracic region of the lung.

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.

On the Design and Testing of an Origami Inspired Nasal Cover: Mitigating Aerosol Risks During Endoscopic Sinus Procedures

Aerosols generated during endoscopic sinus procedures present a concern to the health and safety of healthcare personnel especially with the presence of COVID-19. The purpose of this study is to describe the design and testing of a nasal cover to restrict aerosolized pathogens. The nasal cover was designed to sit overtop the nose with conformal slits for insertion of endoscopic instrumentation. Polydimethylsiloxane (PDMS) was chosen as the nasal mask material and its composition, thickness, and slit geometry were selected using a Taguchi experimental design and survey with clinical collaborators at Penn State Milton S. Hershey Medical Center. The nasal cover was designed using principles of origami engineering to be manufactured flat then folded into its operating state. Form and functionality were evaluated by surgeons, fellows, and residents in the aforementioned survey. Aerosol containment was evaluated by measuring smoke, representative of surgical aerosols, with an optical particle counter. A 25:1 composition PDMS with 3mm thickness and vertical slit geometry was chosen for the nasal cover design. Survey results demonstrated that the origami cover sat well on the nose and did not significantly impact the surgical conditions with single instrumentation. On average, this nasal cover was found to restrict more than 93% of 0.3µm aerosols, and more than 99% of all aerosols larger than 0.5µm in size. Use of a patient worn nasal cover has the potential to drastically reduce the risk to hospital personnel during endonasal surgeries by reducing aerosol generation and potential pathogen spread.

Computational Modeling of Nasal Drug Delivery Using Different Intranasal Corticosteroid Sprays for the Treatment of Eustachian Tube Dysfunction

Eustachian tube dysfunction (ETD) is a common otolaryngologic condition associated with decreased quality of life. The first-line treatment of ETD is intranasal corticosteroid sprays (INCS). Computational fluid dynamics (CFD) was used to study particle deposition on the Eustachian tube (ET) using two commercial INCS (Flonase and Sensimist). Simulations also considered the effects of nostril side, insertion depth, insertion angle, cone spray angle, inhaling rates, wall impingement treatment, and fluid film. Flonase and Sensimist produced different particle size distributions and sizes. Sensimist droplets are smaller, less sensitive to asymmetry in nostrils anatomy and variation in insertion angle, and therefore can reach the posterior nasopharynx more readily. Flonase produces larger particles with greater inertia. Its particles deposition is more sensitive to intrasubject variation in nasal anatomy and insertion angles. The particle deposition on the ET was sensitive to the wall impingement model. The deposition on the ET was insignificant with adherence only <0.15% but increased up to 1-4% when including additional outcomes rebound and splash effects when droplets impact with the wall. The dose redistribution with the fluid film is significant but plays a secondary effect on the ET deposition. Flonase aligned parallel with the hard palate produced 4% deposition efficiency on the ET, but this decreased <0.14% at the higher insertion angle. INCS with larger droplet sizes with a small insertion angle may be more effective at targeting droplet deposition on the ET opening.

Comparison of rhinomanometric and computational fluid dynamic assessment of nasal resistance with respect to measurement accuracy

PURPOSE

Computational fluid dynamics (CFD)-based calculation of intranasal airflow became an important method in rhinologic research. Current evidence shows weak to moderate correlation as well as a systematic underprediction of nasal resistance by numerical simulations. In this study, we investigate whether these differences can be explained by measurement uncertainties caused by rhinomanometric devices and procedures. Furthermore, preliminary findings regarding the impact of tissue movements are reported.

METHODS

A retrospective sample of 17 patients, who reported impaired nasal breathing and for which rhinomanometric (RMM) measurements using two different devices as well as computed tomography scans were available, was investigated in this study. Three patients also exhibited a marked collapse of the nasal valve. Agreement between both rhinomanometric measurements as well as between rhinomanometry and CFD-based calculations was assessed using linear correlation and Bland-Altman analyses. These analyses were performed for the volume flow rates measured at trans-nasal pressure differences of 75 and 150 Pa during inspiration and expiration.

RESULTS

The correlation between volume flow rates measured using both RMM devices was good (R² > 0.72 for all breathing states), and no relevant differences in measured flow rates was observed (21.6 ml/s and 14.8 ml/s for 75 and 150 Pa, respectively). In contrast, correlation between RMM and CFD was poor (R² < 0.5) and CFD systematically overpredicted RMM-based flow rate measurements (231.8 ml/s and 328.3 ml/s). No differences between patients with and without nasal valve collapse nor between inspiration and expiration were observed.

CONCLUSION

Biases introduced during RMM measurements, by either the chosen device, the operator or other aspects as for example the nasal cycle, are not strong enough to explain the gross differences commonly reported between RMM- and CFD-based measurement of nasal resistance. Additionally, tissue movement during breathing is most likely also no sufficient explanation for these differences.

Numerical Simulation of Aspergillus Niger Spore Deposition in Nasal Cavities of a Population in Northwest China

Background: As common pathogens in the human respiratory tract, fungal-spore-related health risks have been challenging to evaluate properly. This paper presents numerical simulations of particle deposition of Aspergillus niger spores in human nasal cavities. Methods: 30 healthy adults (including 60 nasal chambers) who lived in northwest China were recruited to conduct a nasal cavity numerical simulation using computational fluid dynamics–discrete phase model (CFD-DPM). The deposition rate in each anatomic area and its influencing variables, such as body position and respiratory flow rate, were analyzed. Results: (1) Under a resting condition, only about 5.57% ± 1.51% Aspergillus niger spores were deposited in the nasal cavity, while most of them escaped from the nasopharynx, and 0.31% ± 0.20% spores entered the maxillary sinus; (2) under an exercising condition, spores deposited in the nasal cavity were about 2.09 times as many as that in the resting state; (3) in a lying position, the A. niger spores deposited evenly on the lateral wall of the nasal cavity and the sinus when compared with a standing position. However, the deposition rate in each anatomic area did not change significantly.

Upper airway flow characteristics of childhood obstructive sleep apnea-hypopnea syndrome

Revealing the structural morphology and inner flow field of the upper airway is important for understanding obstructive sleep apnea-hypopnea syndrome (OSAHS) incidence phenomena and pathological diagnosis in children. However, prior work on this topic has been focused on adults and the findings cannot be directly extrapolated to children because of different inducing factors. Therefore, this paper employs a simulation method to investigate upper airway flow characteristics of childhood OSAHS. It is found that the Reynold number changes highly throughout the whole upper airway, and the laminar assumption is no longer suitable for low Reynold number flow, which is much unlike classic fluid mechanics. Turbulent models of Standard k-ω and Spalart-Allmaras were developed prior to suggestion. The simulation is validated by experiments with an error of approximately 20%. Additionally, carried out in this analysis is the influence of adenoidal hypertrophy with different narrow levels. The cross-sectional area, flow velocity, pressure drop and volume rate will change greatly when the narrow level is above 64% of the upper airway, which can be a quantitative explanation for medical intervention if adenoid hypertrophy blocks 2/3 of the upper airway in the common clinical judgment of otorhinolaryngology. It is expected that this paper can be a meaningful instruction on OSAHS surgery plan making as well as recovery evaluation postoperatively.

Human Multi-Compartment Airways-on-Chip Platform for Emulating Respiratory Airborne Transmission: From Nose to Pulmonary Acini

The past decade has witnessed tremendous endeavors to deliver novel preclinical in vitro lung models for pulmonary research endpoints, including foremost with the advent of organ- and lung-on-chips. With growing interest in aerosol transmission and infection of respiratory viruses within a host, most notably the SARS-CoV-2 virus amidst the global COVID-19 pandemic, the importance of crosstalk between the different lung regions (i.e., extra-thoracic, conductive and respiratory), with distinct cellular makeups and physiology, are acknowledged to play an important role in the progression of the disease from the initial onset of infection. In the present Methods article, we designed and fabricated to the best of our knowledge the first multi-compartment human airway-on-chip platform to serve as a preclinical in vitro benchmark underlining regional lung crosstalk for viral infection pathways. Combining microfabrication and 3D printing techniques, our platform mimics key elements of the respiratory system spanning (i) nasal passages that serve as the alleged origin of infections, (ii) the mid-bronchial airway region and (iii) the deep acinar region, distinct with alveolated airways. Crosstalk between the three components was exemplified in various assays. First, viral-load (including SARS-CoV-2) injected into the apical partition of the nasal compartment was detected in distal bronchial and acinar components upon applying physiological airflow across the connected compartment models. Secondly, nebulized viral-like dsRNA, poly I:C aerosols were administered to the nasal apical compartment, transmitted to downstream compartments via respiratory airflows and leading to an elevation in inflammatory cytokine levels secreted by distinct epithelial cells in each respective compartment. Overall, our assays establish an in vitro methodology that supports the hypothesis for viral-laden airflow mediated transmission through the respiratory system cellular landscape. With a keen eye for broader end user applications, we share detailed methodologies for fabricating, assembling, calibrating, and using our multi-compartment platform, including open-source fabrication files. Our platform serves as an early proof-of-concept that can be readily designed and adapted to specific preclinical pulmonary research endpoints.

Evidence of Nasal Cooling and Sensory Impairments Driving Patient Symptoms With Septal Deviation

OBJECTIVES/HYPOTHESIS

About 260,000 septoplasties are performed annually in the US to address nasal septal deviation (NSD). Yet, we do not consistently understand what aspects of NSD result in symptoms.

STUDY DESIGN

Blinded cohort study.

METHODS

Two fellowship-trained surgeons blindly reviewed computerized tomography (CTs) of 10 confirmed NSD patients mixed with 36 healthy controls. All patients were correctly identified, however, 24/36 controls were falsely identified by both surgeons as patients (33.3% specificity), which were grouped as asymptomatic NSD (aNSD), while the remaining controls as non-NSD (healthy). Acoustic rhinometry, rhinomanometry, individual CT-based computational fluid dynamics and nasal sensory testing were applied to address the puzzling questions of why these aNSD had no symptoms and, more fundamentally, what caused symptoms in sNSD patients.

RESULTS

aNSD reported no nasal symptoms - Nasal Obstruction Symptom Evaluation score (sNSD: 60.50 ± 13.00; aNSD: 5.20 ± 5.41; non-NSD: 6.66 ± 7.17, P < .05); 22-item Sino-Nasal Outcome Test score (sNSD: 32.60 ± 14.13; aNSD: 10.04 ± 10.10; non-NSD: 9.08 ± 12.42, P < .001). No significant differences in measured nasal resistance, minimum cross-sectional area (MCA), degree of septal deviation, and nasal airflow distributions were found between sNSD and aNSD groups. Only three variables differentiate sNSD versus aNSD: anterior averaged heat flux on deviated side, inferior turbinate peak heat flux on non-deviated side, and nasal cool sensitivity measured by menthol lateralization threshold, with no significant differences among these variables found between the two healthy groups (aNSD vs. non-NSD). These variables by themselves or combined can differentiate sNSD from controls with higher specificity than the physicians (ROC area under the curve = 0.84 with 70% sensitivity and 91.6% specificity).

CONCLUSIONS

This study sheds light on the potential mechanisms of NSD symptomatology: distorted nasal cooling due to NSD exacerbated by poorer nasal mucosal sensitivity. It further supports our previous hypothesis that nasal obstruction complaints do not result directly from obstruction, rather from the capacity of our nose to subjectively sense airflow cooling.

LEVEL OF EVIDENCE

3 Laryngoscope, 132:509-517, 2022.

Numerical comparison of inspiratory airflow patterns in human nasal cavities with distinct age differences

As a primary determinant of nasal physiological functions, the nasal morphology and its effects on the airflow dynamics have been extensively studied in literature. However, gross flow features reported in literature are mostly obtained from subjects at similar ages, while studies focusing on nasal subjects with distinct age differences are significantly less. To advance current understandings of nasal airflow dynamics in the context of age diversity, this study employed three anatomically accurate nasal cavity models with distinct age features (5-, 24- and 77-year-old models) and numerically compared the physiological nasal airflow fields within these nasal cavity models. To demonstrate the validity of the present numerical models, in vivo rhinomanometry measurement was conducted on the 24-year-old female nasal model, and key anatomical features and pressure-flow curves of all three models were compared with models with similar age features in literature work. Apart from results comparison based on conventional velocity flow fields and wall shear stress distributions, a method for quantifying flow partitions in confined airway spaces was developed to reveal the proportions of fractional flow that enters the olfactory region. Our results revealed dramatic intersubject discrepancies between considered nasal cavity models, especially for the fractional flow that enters the olfactory region. Specifically, the 5-year-old girl nasal model received the highest proportion of fractional flow, which accounts for 13.3% ~ 15% of overall inhalation flow rates under different activity levels. For the 24-year-old female model, on the contrary, the olfactory fractional flow was dramatically reduced (with a local to overall percentage around 4.3%-7.7%). Finally, for the elderly subject-77-year-old male model, minimum level of olfactory flux was observed with a local to overall percentage ranging between 3.1% and 4.9% for considered wide range of inhalation flow rates. Therefore, the local flow intersubject variation can reach nearly fourfold. The vast local flow difference is mainly due to the inherent anatomical features (e.g., immature nasal turbinate structure in the child model, the partial narrowing superior nasal valve in the elder model). The results may further lead to discrepant health effects associated with inhalation exposure to airborne particles.

In Silico Study to Enhance Delivery Efficiency of Charged Nanoscale Nasal Spray Aerosols to the Olfactory Region Using External Magnetic Fields

Various factors and challenges are involved in efficiently delivering drugs using nasal sprays to the olfactory region to treat central nervous system diseases. In this study, computational fluid dynamics was used to simulate nasal drug delivery to (1) examine effects on drug deposition when various external magnetic fields are applied to charged particles, (2) comprehensively study effects of multiple parameters (i.e., particle aerodynamic diameter; injection velocity magnitude, angle, and position; magnetic force strength and direction), and (3) determine how to achieve the optimal delivery efficiency to the olfactory epithelium. The Reynolds-averaged Navier–Stokes equations governed airflow, with a realistic inhalation waveform implemented at the nostrils. Particle trajectories were modeled using the one-way coupled Euler–Lagrange model. A current-carrying wire generated a magnetic field to apply force on charged particles and direct them to the olfactory region. Once drug particles reached the olfactory region, their diffusion through mucus to the epithelium was calculated analytically. Particle aerodynamic diameter, injection position, and magnetic field strength were found to be interconnected in their effects on delivery efficiency. Specific combinations of these parameters achieved over 65-fold higher drug delivery efficiency compared with uniform injections with no magnetic fields. The insight gained suggests how to integrate these factors to achieve the optimal efficiency.

A machine-learning-based method for automatizing lattice-Boltzmann simulations of respiratory flows

Many simulation workflows require to prepare the data for the simulation manually. This is time consuming and leads to a massive bottleneck when a large number of numerical simulations is requested. This bottleneck can be overcome by an automated data processing pipeline. Such a novel pipeline is developed for a medical use case from rhinology, where computer tomography recordings are used as input and flow simulation data define the results. Convolutional neural networks are applied to segment the upper airways and to detect and prepare the in- and outflow regions for accurate boundary condition prescription in the simulation. The automated process is tested on three cases which have not been used to train the networks. The accuracy of the pipeline is evaluated by comparing the network-generated output surfaces to those obtained from a semi-automated procedure performed by a medical professional. Except for minor deviations at interfaces between ethmoidal sinuses, the network-generated surface is sufficiently accurate. To further analyze the accuracy of the automated pipeline, flow simulations are conducted with a thermal lattice-Boltzmann method for both cases on a high-performace computing system. The comparison of the results of the respiratory flow simulations yield averaged errors of less than 1% for the pressure loss between the in- and outlets, and for the outlet temperature. Thus, the pipeline is shown to work accurately and the geometrical deviations at the ethmoidal sinuses to be negligible.

Use of computational fluid dynamics (CFD) to model observed nasal nitric oxide levels in human subjects

BACKGROUND

Upper airway nitric oxide (NO) is physiologically important in airway regulation and defense, and nasal NO (nNO) levels typically exceed those in exhaled breath (fractional exhaled NO [FeNO]). Elevated concentrations of NO sampled from the nose, in turn, reflect even higher concentrations in the paranasal sinuses, suggesting a "reservoir" role for the latter. However, the dynamics of NO flux within the sinonasal compartment are poorly understood.

METHODS

Data from 10 human subjects who had previously undergone both real-time nNO sampling and computed tomography (CT) scanning of the sinuses were analyzed using computational fluid dynamics (CFD) methods. Modeled and observed nNO values during the initial 2-s transient ("spike") during nasal exhalation were then compared.

RESULTS

Examining the initial 2-s transient spike for each subject (as well as the pooled group), there was a statistically significant correlation between modeled and observed nNO levels, with r values ranging from 0.43 to 0.89 (p values ranging from <0.05 to <0.0001). Model performance varied between subjects, with weaker correlations evident in those with high background (FeNO) levels. In addition, the CFD simulation suggests that ethmoid sinuses (>60%) and diffusion process (>54%) contributed most to total nasal NO emissions.

CONCLUSION

Analysis of this dataset confirms that CFD is a valuable modeling tool for nNO dynamics, and highlights the importance of the ethmoid sinuses, as well as the role of diffusion as an initiating step in sinonasal NO flux. Future model iterations may apply more generally if baseline FeNO is taken into account.

Detailed comparison of anatomy and airflow dynamics in human and cynomolgus monkey nasal cavity

Nonhuman primates are occasionally used as laboratory models for sophisticated medical research as they bear the closest resemblance to humans in morphometry and physiological functions. A range of nonhuman primate species have been employed in the inhalation toxicity, nasal drug delivery and respiratory viral infection studies, and they provided valuable insight to disease pathogenesis while other laboratory animals such as rodents cannot recapitulate due to the lesser degree of similarity in metabolism, anatomy and cellular response to that of humans. It is anticipated that nonhuman primate models of respiratory diseases will continue to be instrumental for translating biomedical research for improvement of human health, and the confidence in laboratory data extrapolation between species will play a pivotal role. From the morphometry and flow dynamics point of view, this study performed a detailed comparative analysis between human and a cynomolgus monkey nasal airway, with intention to provide high-fidelity qualitative and quantitative linkage between the two species for more effective laboratory data extrapolation. The study revealed that cynomolgus monkey could be a good human surrogate in nasal inhalation studies; however, care should be given for interspecies data extrapolation as subtle differences in anatomy and airflow dynamics were present between the two species.

Septal deviation in the nose of the longest faced crocodylian: A description of nasal anatomy and airflow in the Indian gharial (Gavialis gangeticus) with comments on acoustics

The remarkably thin rostrum in the Indian gharial (Gavialis gangeticus) imparts challenges to nasal physiology. Competition for space in the slim jaws necessitates a thin nasal septum, leaving this taxon susceptible to nasal passage abnormalities such as septal deviation. Here we describe the nasal anatomy of gharials based on multiple individuals including one that showcases an extreme instance of nasal septum deviation. We found that gharials have both confluent nostrils and choanae, which may be important for their unique nasal acoustics. The deviated nasal septum in the female showed distinct waviness that affected the nasal passages by alternately compressing them. We performed a computational fluid dynamic analysis on the nasal passages to visualize the effects of septal deviation on airflow. Our analysis found the deviated septum increased nasal resistance and wall shear stress during respiration, resulting in unequal distribution of the air field between both sides of the nasal passage. Our findings indicate that gharials-and potentially other longirostrine crocodylians-may be particularly susceptible to septal deviations. Lastly, we observed pterygoid bullae to be present in both sexes, though their morphology differed. Airflow in the male pterygoid bullae produced a Bernoulli effect which may be responsible for the unique "pop" sounds recorded in this species.

Computational Fluid Dynamic Modeling Reveals Nonlinear Airway Stress during Trachea Development

Normative trachea dimensions and aerodynamic information during development was collected to establish clinical benchmarks and showed that airway development seems to outpace respiratory demands. Infants and toddlers' trachea exhibit higher aerodynamic stress that significantly decreases by teenage years. This implies large airway pathology in younger children may have a more substantial clinical impact.

Inhalation and deposition of spherical and pollen particles after middle turbinate resection in a human nasal cavity

Middle turbinate resection significantly alters the anatomy and redistributes the inhaled air. The superior half of the main nasal cavity is opened up, increasing accessibility to the region. This is expected to increase inhalation dosimetry to the region during exposure to airborne particles. This study investigated the influence of middle turbinate resection on the deposition of inhaled pollutants that cover spherical and non-spherical particles (e.g. pollen). A computational model of the nasal cavity from CT scans, and its corresponding post-operative model with virtual surgery performed was created. Two constant flow rates of 5 L/min, and 15 L/min were simulated under a laminar flow field. Inhaled particles including pollen (non-spherical), and a spherical particle with reference density of 1000 kg/m³ were introduced in the surrounding atmosphere. The effect of surgery was most prominent in the less patent cavity side, since the change in anatomy was proportionally greater relative to the original airway space. The left cavity produced an increase in particle deposition at a flow rate of 15 L/min. The main particle deposition mechanisms were inertial impaction, and to a lesser degree gravitational sedimentation. The results are expected to provide insight into inhalation efficiency of different aerosol types, and the likelihood of deposition in different nasal cavity surfaces.

Computational modeling of nasal nitric oxide flux from the paranasal sinuses: Validation against human experiment

BACKGROUND

Nitric oxide (NO) is important in respiratory physiology and airway defense. Although the paranasal sinuses are the major source of nasal NO, transport dynamics between the sinuses and nasal cavities are poorly understood.

METHODS

Exhaled nasal NO tracings were measured in two non-asthmatic subjects (one with allergic rhinitis, one without) using NO analyzer connected via face mask. We subsequently performed computational fluid dynamics NO emission simulations based on individual CT scans and compared to the experimental data.

RESULTS

Simulated exhaled NO tracings match well with experimental data (r > 0.84, p < 0.01) for both subjects, with measured peaks reaching 319.6 ppb in one subject (allergic-rhinitis), and 196.9 ppb in the other. The CFD simulation accurately captured the peak differences, even though the initial sinus NO concentration for both cases was set to the same 9000 ppb based on literature value. Further, the CFD simulation suggests that ethmoid sinuses contributed the most (>67%, other sinuses combined <33%) to total nasal NO emission in both cases and that diffusion contributes more than convective transport. By turning off diffusion (setting NO diffusivity to ~0), the NO emission peaks for both cases were reduced by >70%.

CONCLUSION

Historically, nasal NO emissions were thought to be contributed mostly by the maxillary sinuses (the largest sinuses) and active air movement (convection). Here, we showed that the ethmoid sinuses and diffusive transport dominate the process. These findings may have a substantial impact on our view of nasal NO emission mechanisms and sinus physiopathology in general.

The effect of decongestion on nasal airway patency and airflow

Nasal decongestant reduces blood flow to the nasal turbinates, reducing tissue volume and increasing nasal airway patency. This study maps the changes in nasal anatomy and measures how these changes affect nasal resistance, flow partitioning between superior and inferior cavity, flow patterns and wall shear stress. High-resolution MRI was applied to capture nasal anatomy in 10 healthy subjects before and after application of a topical decongestant. Computational fluid dynamics simulated nasal airflow at steady inspiratory flow rates of 15 L.min[Formula: see text] and 30 L.min[Formula: see text]. The results show decongestion mainly increases the cross-sectional area in the turbinate region and SAVR is reduced (median approximately 40[Formula: see text] reduction) in middle and lower parts of the cavity. Decongestion reduces nasal resistance by 50[Formula: see text] on average, while in the posterior cavity, nasal resistance decreases by a median factor of approximately 3 after decongestion. We also find decongestant regularises nasal airflow and alters the partitioning of flow, significantly decreasing flow through the superior portions of the nasal cavity. By comparing nasal anatomies and airflow in their normal state with that when pharmacologically decongested, this study provides data for a broad range of anatomy and airflow conditions, which may help characterize the extent of nasal variability.

Leveraging statistical shape modeling in computational respiratory dynamics: Nanomedicine delivery in remodeled airways

BACKGROUND AND OBJECTIVE

Accurate knowledge of the delivered doses to the diseased site in the respiratory tract is crucial to elicit desired therapeutic outcomes. However, such information is still difficult to obtain due to inaccessibility for measurement or visualization, complex network structure, and challenges in reconstructing lung geometries with disease-invoked airway remodeling. This study presents a novel method to simulate the airway remodeling in a mouth-lung geometry extending to G9.

METHODS

Statistical shape modeling was used to extract morphological features from a lung geometry database and four new models (i.e., M1-M4) were generated with parameter-controlled dilated/constricted bronchioles in the left-lower (LL) lung. The variations in airflow and particle deposition due to the airway remodeling were simulated using a well-tested k-ω turbulence model and a Lagrangian tracking approach.

RESULTS

Significant variations in flow partitions between the lower and upper lobes of the left lung, as well as between the left and right lungs. The flow partition into the LL lobe varied by 10-fold between the most dilated and constricted models in this study. Significantly lower doses were also predicted on the surface of the constricted LL bronchioles G4-G9, as well as into the peripheral airways beyond G9. However, the total dosimetry in the mouth-lung geometry (up to G9) exhibited low sensitivity to the LL lobar remodeling. Results in this study suggest that the optimal nanomedicine should be 2-10 nm in diameter if targeted at the constricted bronchioles G4-G9 as in topical inhalation therapy but should be larger than 20 nm if targeted at the alveolar region as in systemic therapy.

CONCLUSIONS

This study highlights the large dose variability from local airway remodeling and the need to consider these variations in the treatment planning for pneumonia and other obstructive respiratory diseases.

Numerical simulation of unsteady airflow in a nasal cavity for various sizes of maxillary sinus opening in a virtual endoscopic surgery

Functional endoscopic sinus surgery (FESS) is performed to treat sinusitis when treatment with medication fails. In the present study, three different virtual maxillary sinus endoscopic surgeries were performed on a realistic 3-D computational model of the nasal cavity of an adult male under the supervision of a specialist. They included only uncinectomy, uncinectomy + 8mm Middle Meatal Antrostomy (MMA) and uncinectomy + 18 mm MMA. Simulations were performed for two human activity respiratory rates, including rest and moderate activities, and effects of different surgeries and respiratory rates on maxillary sinus were investigated. It was found that after endoscopic sinus surgery, the volume of air entering the maxillary sinus increased significantly, and as the size of the MMA increased, or the breathing condition changed from rest to moderate activity, this volume of air increased. For the rest condition, on average for both nasal passages, for uncinectomy +8 mm MMA, around 15 % of the inhaled flow and 7 % of the exhaled flow enter the maxillary sinuses. For uncinectomy +18 mm MMA, these values are 24 % and 14 %, respectively. As human activity increases, a lower portion of inhaled and exhaled air enters the maxillary sinuses. For the moderate activity condition, on average for both nasal passages, for uncinectomy +8 mm MMA, around 11 % of the inhaled flow and 6 % of the exhaled flow rate enters the maxillary sinus. For uncinectomy +18 mm MMA, these values are 16 % and 8%, respectively. Comparing the steady and unsteady simulation results showed that the quasi-steady flow assumption could predict the flow in the maxillary sinus and the volume of air entering the sinuses, almost at any moment of respiration, with acceptable accuracy.

The Virtual Nose: Assessment of Static Nasal Airway Obstruction Using Computational Simulations and 3D-Printed Models

Background: The use of virtual noses to predict the outcome of surgery is of increasing interests, particularly, as detailed and objective pre- and postoperative assessments of nasal airway obstruction (NAO) are difficult to perform. The objective of this article is to validate predictions using virtual noses against their experimentally measured counterpart in rigid 3D-printed models. Methods: Virtual nose models, with and without NAO, were reconstructed from patients' cone beam computed tomography scans, and used to evaluate airflow characteristics through computational fluid dynamics simulations. Prototypes of the reconstructed models were 3D printed and instrumented experimentally for pressure measurements. Results: Correlation between the numerical predictions and experimental measurements was shown. Analysis of the flow field indicated that the NAO in the nasal valve increases significantly the wall pressure, shear stress, and incremental nasal resistance behind the obstruction. Conclusions: Airflow predictions in static virtual noses correlate well with detailed experimental measurements on 3D-printed replicas of patient airways.

Inhalation dosimetry of nasally inhaled respiratory aerosols in the human respiratory tract with locally remodeled conducting lungs

Objective: Respiratory diseases are often accompanied by alterations to airway morphology. However, inhalation dosimetry data in remodeled airways are scarce due to the challenges in reconstructing diseased respiratory morphologies. This study aims to study the airway remodeling effects on the inhalation dosimetry of nasally inhaled nanoparticles in a nose-lung geometry that extends to G9 (ninth generation).Materials and methods: Statistical shape modeling was used to develop four diseased lung models with varying levels of bronchiolar dilation/constriction in the left-lower (LL) lobe (i.e. M1-M4). Respiratory airflow and particle deposition were simulated using a low Reynolds number k-ω turbulence model and a Lagrangian tracking approach.Results: Significant discrepancies were observed in the flow partitions between the left and right lungs, as well as between the lower and upper lobes of the left lung, which changed by 10-fold between the most dilated and constricted models.Much lower doses were predicted on the surface of the constricted LL bronchioles G4-G9, as well as into the peripheral airways beyond G9 of the LL lung. However, the LL lobar remodeling had little effect on the dosimetry in the nasopharynx, as well as on the total dosimetry in the nose-lung geometry (up to G9).Conclusion: It is suggested that airway remodeling may pose a higher viral infection risk to the host by redistributing the inhaled viruses to healthy lung lobes. Airway remodeling effects should also be considered in the treatment planning of inhalation therapies, not only because of the dosimetry variation from altered lung morphology but also its evolution as the disease progresses.

Pre-surgery planning tool for estimation of resection volume to improve nasal breathing based on lattice Boltzmann fluid flow simulations

Purpose

State-of-the-art medical examination techniques (e.g., rhinomanometry and endoscopy) do not always lead to satisfactory postoperative outcome. A fully automatized optimization tool based on patient computer tomography (CT) data to calculate local pressure gradient regions to reshape pathological nasal cavity geometry is proposed.

Methods

Five anonymous pre- and postoperative CT datasets with nasal septum deviations were used to simulate the airflow through the nasal cavity with lattice Boltzmann (LB) simulations. Pressure gradient regions were detected by a streamline analysis. After shape optimization, the volumetric difference between the two shapes of the nasal cavity yields the estimated resection volume.

Results

At LB rhinomanometry boundary conditions (bilateral flow rate of 600 ml/s), the preliminary study shows a critical pressure gradient of −1.1 Pa/mm as optimization criterion. The maximum coronal airflow ΔA  := cross-section ratio \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{\mathrm{virtual surgery }}{\mathrm{post}-\mathrm{surgery}}$$\end{document}virtualsurgerypost-surgery found close to the nostrils is 1.15. For the patients a pressure drop ratio ΔΠ  := (pre-surgery − virtual surgery)/(pre-surgery − post-surgery) between nostril and nasopharynx of 1.25, 1.72, −1.85, 0.79 and 1.02 is calculated.

Conclusions

LB fluid mechanics optimization of the nasal cavity can yield results similar to surgery for air-flow cross section and pressure drop between nostril and nasopharynx. The optimization is numerically stable in all five cases of the presented study. A limitation of this study is that anatomical constraints (e.g. mucosa) have not been considered.

Agreement between rhinomanometry and computed tomography-based computational fluid dynamics

Purpose

Active anterior rhinomanometry (AAR) and computed tomography (CT) are standardized methods for the evaluation of nasal obstruction. Recent attempts to correlate AAR with CT-based computational fluid dynamics (CFD) have been controversial. We aimed to investigate this correlation and agreement based on an in-house developed procedure.

Methods

In a pilot study, we retrospectively examined five subjects scheduled for septoplasty, along with preoperative digital volume tomography and AAR. The simulation was performed with Sailfish CFD, a lattice Boltzmann code. We examined the correlation and agreement of pressure derived from AAR (RhinoPress) and simulation (SimPress) and these of resistance during inspiration by 150 Pa pressure drop derived from AAR (RhinoRes150) and simulation (SimRes150). For investigation of correlation between pressures and between resistances, a univariate analysis of variance and a Pearson’s correlation were performed, respectively. For investigation of agreement, the Bland–Altman method was used.

Results

The correlation coefficient between RhinoPress and SimPress was r = 0.93 (p < 0.001). RhinoPress was similar to SimPress in the less obstructed nasal side and two times greater than SimPress in the more obstructed nasal side. A moderate correlation was found between RhinoRes150 and SimRes150 (r = 0.65; p = 0.041).

Conclusion

The simulation of rhinomanometry pressure by CT-based CFD seems more feasible with the lattice Boltzmann code in the less obstructed nasal side. In the more obstructed nasal side, error rates of up to 100% were encountered. Our results imply that the pressure and resistance derived from CT-based CFD and AAR were similar, yet not same.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11548-021-02332-1.

[Analysis of upper airway flow field between obstructive sleep apnea and normal children based on computational fluid dynamics]

Objective:In this study, the characteristics of the upper airway flow field were analyzed by using computational fluid dynamics（CFD）. The study analyze the differences in the upper airway flow field between normal children and children with obstructive sleep apnea（OSA）, and the pathological characteristics of children with OSA were elaborated from the perspective of airway fluid dynamics. Methods:The upper airway models of a normal child and a child with OSA were constructed. The differences in the same inspiration pressure, such as airflow velocity, airflow pattern, ventilation volume, and pressure, were analyzed. To verify CFD results, rhinomanometry was carried out and an experimental bench based 3D technology was also built. Results:The CFD results are consistent with the in vitro 3D model experiments and clinical measurement results. The adenoid area of nasopharynx is only 11.274 mm²of the child with OSA, about 1/6 of that of normal children. At the area of nasopharyngeal in OSA children, the flow velocity increased but the pressure dropped sharply, which was 69.197% of the total pressure drop, and the resistance value was 6.59 times of that of normal children. Streamline of nasopharyngeal is more disorder. Normal children's inspiratory flow was 116.139 mL/s, while OSA children's inspiratory flow was 47.055 mL/s, with a difference rate as high as 59.48%. Conclusion:The airflow of OSA children in nasopharynx is significantly different from that of normal children. The airflow characteristics of upper airway were discussed in detail with the use of CFD, which can help clinicians intuitively understand the abnormal flow behavior of children with OSA.

SARS COV-2 virus-laden droplets coughed from deep lungs: Numerical quantification in a single-path whole respiratory tract geometry

When an infected person coughs, many virus-laden droplets will be exhaled out of the mouth. Droplets from deep lungs are especially infectious because the alveoli are the major sites of coronavirus replication. However, their exhalation fraction, size distribution, and exiting speeds are unclear. This study investigated the behavior and fate of respiratory droplets (0.1-4 μm) during coughs in a single-path respiratory tract model extending from terminal alveoli to mouth opening. An experimentally measured cough waveform was used to control the alveolar wall motions and the flow boundary conditions at lung branches from G2 to G18. The mouth opening was modeled after the image of a coughing subject captured using a high-speed camera. A well-tested k-ω turbulence model and Lagrangian particle tracking algorithm were applied to simulate cough flow evolutions and droplet dynamics under four cough depths, i.e., tidal volume ratio (TVR) = 0.13, 0.20. 0.32, and 0.42. The results show that 2-μm droplets have the highest exhalation fraction, regardless of cough depths. A nonlinear relationship exists between the droplet exhalation fraction and cough depth due to a complex deposition mechanism confounded by multiscale airway passages, multiregime flows, and drastic transient flow effects. The highest exhalation fraction is 1.6% at the normal cough depth (TVR = 0.32), with a mean exiting speed of 20 m/s. The finding that most exhaled droplets from deep lungs are 2 μm highlights the need for more effective facemasks in blocking 2-μm droplets and smaller both in infectious source control and self-protection from airborne virus-laden droplets.