Mechanics of airflow in the human nasal airways

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.

Assessing nasal airway resistance and symmetry: An approach to global perspective through computational fluid dynamics

This study aimed to explore the variability in nasal airflow patterns among different sexes and populations using computational fluid dynamics (CFD). We focused on evaluating the universality and applicability of dimensionless parameters R (bilateral nasal resistance) and ϕ (nasal flow asymmetry), initially established in a Caucasian Spanish cohort, across a broader spectrum of human populations to assess normal breathing function in healthy airways. In this retrospective study, CT scans from Cambodia (20 males, 20 females), Russia (20 males, 18 females), and Spain (19 males, 19 females) were analyzed. A standardized CFD workflow was implemented to calculate R-ϕ parameters from these scans. Statistical analyses were conducted to assess and compare these parameters across different sexes and populations, emphasizing their distribution and variances. Our results indicated no significant sex-based differences in the R parameter across the populations. However, moderate sexual dimorphism in the ϕ parameter was observed in the Cambodian group. Notably, no geographical differences were found in either R or ϕ parameters, suggesting consistent nasal airflow characteristics across the diverse human groups studied. The study also emphasized the importance of using dimensionless variables to effectively analyze the relationships between form and function in nasal airflow. The observed consistency of R-ϕ parameters across various populations highlights their potential as reliable indicators in both medical practice and further CFD research, particularly in diverse human populations. Our findings suggest the potential applicability of dimensionless CFD parameters in analyzing nasal airflow, highlighting their utility across diverse demographic and geographic contexts. This research advances our understanding of nasal airflow dynamics and underscores the need for additional studies to validate these parameters in broader population cohorts. The approach of employing dimensionless parameters paves the way for future research that eliminates confounding size effects, enabling more accurate comparisons across different populations and sexes. The implications of this study are significant for the advancement of personalized medicine and the development of diagnostic tools that accommodate individual variations in nasal airflow.

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.

Beyond skeletal studies: A computational analysis of nasal airway function in climate adaptation

OBJECTIVES

Ecogeographic variation in human nasal anatomy has historically been analyzed on skeletal morphology and interpreted in the context of climatic adaptations to respiratory air-conditioning. Only a few studies have analyzed nasal soft tissue morphology, actively involved in air-conditioning physiology.

MATERIALS AND METHODS

We used in vivo computer tomographic scans of (N = 146) adult individuals from Cambodia, Chile, Russia, and Spain. We conducted (N = 438) airflow simulations during inspiration using computational fluid dynamics to analyze the air-conditioning capacities of the nasal soft tissue in the inflow, functional, and outflow tract, under three different environmental conditions: cold-dry; hot-dry; and hot-humid. We performed statistical comparisons between populations and sexes.

RESULTS

Subjects from hot-humid regions showed significantly lower air-conditioning capacities than subjects from colder regions in all the three conditions, specifically within the isthmus region in the inflow tract, and the anterior part of the internal functional tract. Posterior to the functional tract, no differences were detected. No differences between sexes were found in any of the tracts and under any of the conditions.

DISCUSSION

Our statistical analyses support models of climatic adaptations of anterior nasal soft tissue morphology that fit with, and complement, previous research on dry skulls. However, our results challenge a morpho-functional model that attributes air-conditioning capacities exclusively to the functional tract located within the nasal cavity. Instead, our findings support studies that have suggested that both, the external nose and the intra-facial soft tissue airways contribute to efficiently warming and humidifying air during inspiration. This supports functional interpretations in modern midfacial variation and evolution.

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.

A Novel Bionebulizer Approach to Study the Effects of Natural Mineral Water on a 3D In Vitro Nasal Model from Allergic Rhinitis Patients

Sulfurous thermal waters (STWs) are used as a complementary treatment for allergic rhinitis. However, there is scant data on the effects of STW on nasal epithelial cells, and in vitro models are warranted. The main aim of this study was to evaluate the dose and time effects of exposure to 3D nasal inserts (MucilAirTM-HF allergic rhinitis model) with STW or isotonic sodium chloride solution (ISCS) aerosols. Transepithelial electrical resistance (TEER) and histology were assessed before and after nebulizations. Chemokine/cytokine levels in the basal supernatants were assessed by enzyme-linked immunosorbent assay. The results showed that more than four daily nebulizations of four or more minutes compromised the normal epithelial integrity. In contrast, 1 or 2 min of STW or ISCS nebulizations had no toxic effect up to 3 days. No statistically significant changes in release of inflammatory chemokines MCP-1/CCL2 > IL-8/CXCL8 > MIP-1α/CCL3, no meaningful release of "alarmins" (IL-1α, IL-33), nor of anti-inflammatory IL-10 cytokine were observed. We have characterized safe time and dose conditions for aerosol nebulizations using a novel in vitro 3D nasal epithelium model of allergic rhinitis patients. This may be a suitable in vitro setup to mimic in vivo treatments of chronic rhinitis with STW upon triggering an inflammatory stimulus in the future.

Impact of nasal septal perforation on the airflow and air-conditioning characteristics of the nasal cavity

We investigated (1) how nasal septal perforations (NSPs) modify nasal airflow and air-conditioning characteristics and (2) how the modifications of nasal airflow are influenced by the size and location of the NSP. Computed tomography scans of 14 subjects with NSPs were used to generate nasal cavity models. Virtual repair of NSPs was conducted to examine the sole effect of NSPs on airflow. The computational fluid dynamics technique was used to assess geometric and airflow parameters around the NSPs and in the nasopharynx. The net crossover airflow rate, the increased wall shear stress (WSS) and the surface water-vapor flux on the posterior surface of the NSPs were not correlated with the size of the perforation. After the virtual closure of the NSPs, the levels in relative humidity (RH), air temperature (AT) and nasal resistance did not improve significantly both in the choanae and nasopharynx. A geometric parameter associated with turbinate volume, the surface area-to-volume ratio (SAVR), was shown to be an important factor in the determination of the RH and AT, even in the presence of NSPs. The levels of RH and AT in the choanae and nasopharynx were more influenced by SAVR than the size and location of the NSPs.

Modelling the effects of post-FESS middle turbinate synechiae on sinonasal physiology: A computational fluid dynamics study

OBJECTIVE(S)

Chronic rhinosinusitis (CRS) is common and often requires surgical intervention. Surgical failure may lead to persistent symptoms and recalcitrant disease, often secondary to synechiae between the middle turbinate (MT) and lateral nasal wall. Synechiae prevention techniques have been extensively investigated, however evidence for the effect of synechiae on sinonasal physiology is lacking. We aimed to model the effects of MT synechiae on a post-functional endoscopic sinus surgery (FESS) sinonasal cavity using computational fluid dynamics (CFD).

METHODS

DICOM data from a CT-sinus of a healthy 25-year-old female was segmented to create a three-dimensional model. Virtual surgery was performed to simulate a "full-house" FESS procedure. Multiple models were created, each with a single unilateral virtual MT synechia of varying extent. CFD analysis was performed on each model and compared with a post-FESS control model without synechiae. Airflow velocity, humidity and mucosal surface and air temperature values were calculated.

RESULTS

All synechiae models demonstrated aberrant downstream sinonasal airflow. There was reduced ventilation of the ipsilateral frontal, ethmoid and sphenoid sinuses, with a concentrated central "jet" in the middle meatus region. Effects were proportionate to the size of synechiae. The impact on bulk inspired airflow was negligible.

CONCLUSION

Post-FESS synechiae between the MT and lateral nasal wall significantly disrupt local downstream sinus ventilation and nasal airflow. These findings may explain the persistent symptoms seen in post-FESS CRS patients with MT synechiae, reinforcing the importance of prevention and adhesiolysis. Larger cohort studies with multiple models of actual post-FESS patients with synechiae are required to validate these findings.

Computational, In Vitro, and In Vivo Models for Nose-to-Brain Drug Delivery Studies

Direct nose-to-brain drug delivery offers the opportunity to treat central nervous system disorders more effectively due to the possibility of drug molecules reaching the brain without passing through the blood-brain barrier. Such a delivery route allows the desired anatomic site to be reached while ensuring drug effectiveness, minimizing side effects, and limiting drug losses and degradation. However, the absorption of intranasally administered entities is a complex process that considerably depends on the interplay between the characteristics of the drug delivery systems and the nasal mucosa. Various preclinical models (in silico, in vitro, ex vivo, and in vivo) are used to study the transport of drugs after intranasal administration. The present review article attempts to summarize the different computational and experimental models used so far to investigate the direct delivery of therapeutic agents or colloidal carriers from the nasal cavity to the brain tissue. Moreover, it provides a critical evaluation of the data available from different studies and identifies the advantages and disadvantages of each model.

Impact of the Location of Nasal Septal Deviation on the Nasal Airflow and Air Conditioning Characteristics

The location of nasal septal deviation (NSD) directly impacts nasal physiology. The objective is to examine, using computational fluid dynamics (CFD), the difference in the airflow and air conditioning characteristics according to the location of NSD. Twenty patients with septal deviation were divided into two: 10 caudal septal deviation (CSD) and 10 posterior septal deviation (PSD). Physiological variables were compared and numerical models for nasal cavity were created with CT scans. Cases with CSD had distinctive features including restricted airflow partition, larger nasal resistance, and decreased surface heat flux in the more obstructed side (MOS), and lower humidity and air temperature in the lesser obstructed side (LOS). Physiological differences were observed according to the location of septal deviation, CSD cases exhibit significantly more asymmetric airflow characteristics and air conditioning capacity between LOS and MOS.

Effect of different degrees of adenoid hypertrophy on pediatric upper airway aerodynamics: a computational fluid dynamics study

To improve the diagnostic accuracy of adenoid hypertrophy (AH) in children and prevent further complications in time, it is important to study and quantify the effects of different degrees of AH on pediatric upper airway (UA) aerodynamics. In this study, based on computed tomography (CT) scans of a child with AH, UA models with different degrees of obstruction (adenoidal-nasopharyngeal (AN) ratio of 0.9, 0.8, 0.7, and 0.6) and no obstruction (AN ratio of 0.5) were constructed through virtual surgery to quantitatively analyze the aerodynamic characteristics of UA with different degrees of obstruction in terms of the peak velocity, pressure drop (△P), and maximum wall shear stress (WSS). We found that two obvious whirlpools are formed in the anterior upper part of the pediatric nasal cavity and in the oropharynx, which is caused by the sudden increase in the nasal cross-section area, resulting in local flow separation and counterflow. In addition, when the AN ratio was ≥ 0.7, the airflow velocity peaked at the protruding area in the nasopharynx, with an increase 1.1-2.7 times greater than that in the nasal valve area; the △P in the nasopharynx was significantly increased, with an increase 1.1-6.8 times greater than that in the nasal cavity; and the maximum WSS of the posterior wall of the nasopharynx was 1.1-4.4 times larger than that of the nasal cavity. The results showed that the size of the adenoid plays an important role in the patency of the pediatric UA.

Effect of transient breathing cycle on the deposition of micro and nanoparticles on respiratory walls

BACKGROUND AND OBJECTIVE

From various perspectives (e.g. inhalation exposure and drug delivery), it is important to provide insights into the behavior of inhaled particles in the human respiratory system. Although most of the experimental and numerical studies have relied on an assumption of steady inhalation, the transient breathing profile is a key factor in particle deposition in the respiratory tract. In this study, particle transportation and deposition were predicted in a realistic human airway model during a breathing cycle and the effects of steady-state and transient flows on the deposition fraction were observed using computational fluid dynamics.

METHODS

Two transient breathing cycles with different respiratory durations were considered to evaluate the effects of respiration duration on particle transport and deposition characteristics. Two types of steady breathing conditions with corresponding steady-state respiratory volumes were reproduced. The Lagrangian discrete phase model approach was used to investigate particle transportation and deposition under transient breathing conditions. Additionally, the Eulerian approach was used to analyze the transport of nanoparticles in the gas phase. A total of >50,000 monodispersed particles with aerodynamic diameters ranging between 2 nm and 10 μm were selected for comprehensive deposition predictions for particle sizes ranging from the nano- to microscale.

RESULTS

The predicted results were compared with the experimental data. The particle deposition fraction in the nasal cavity and tracheal regions showed differences between the steady and transient simulations. In addition, particle analysis under steady inhalation conditions cannot accurately predict particle transportation and deposition in the lower airway. Furthermore, the breathing cycle had a significant effect on the deposition fraction of the particles and the behavior of the inhaled particles.

CONCLUSIONS

Transient simulation mimicking the breathing cycle was observed to be an important factor in accurately predicting the transportation and deposition of particles in the respiratory tract.

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.

Computational analysis of human upper airway aerodynamics

There is a considerable interest in understanding transient human upper airway aerodynamics, especially in view of assessing the effects of various ventilation therapies. Experimental analyses in a patient-specific manner pose challenges as the upper airway consists of a narrow confined region with complex anatomy. Pressure measurements are feasible, but, for example, PIV experiments require special measures to accommodate for the light refraction by the model. Computational fluid dynamics can bridge the gap between limited experimental data and detailed flow features. This work aims to validate the use of combined lattice Boltzmann method and a large eddy scale model for simulating respiration, and to identify clinical features of the flow and show the clinical potential of the method. Airflow was computationally analyzed during a realistic, transient, breathing profile in an upper airway geometry ranging from nose to trachea, and the resulting pressure calculations were compared against in vitro experiments. Simulations were conducted on meshes containing about 1 billion cells to ensure accuracy and to capture intrinsic flow features. Airway pressures obtained from simulations and in vitro experiments are in good agreement both during inhalation and exhalation. High velocity pharyngeal and laryngeal jets and recirculation in the region of the olfactory cleft are observed. Graphical Abstract The Lattice-Boltzmann Method combined with Large Eddy Simulations was used to compute the aerodynamics in a human upper airway geometry. The left side of this graphical abstract shows the velocity and vorticity (middle figure in bottom row, and right figure of the right bottom figure) profiles at peak exhalation. The simulations were validated against experiments on a 3D-print of the geometry (shown in the top figures on the right hand side). The pressure drop (right bottom corner) shows a good agreement between experiments and simulations.

INFLUENCE OF BAD HABITS ON THE DEVELOPMENT OF ACQUIRED DEFORMATIONS IN THE MAXILLOFACIAL AREA

OBJECTIVE

The aim: Identifying the relationship between the presence of oral habit and acquired maxillomandibular anomalies, influence of oral habits on the skeleton and muscular system formation in children.

PATIENTS AND METHODS

Materials and methods: We conducted clinical, radiological methods of examination of 60 patients aged 9-12 with acquired maxillomandibular anomalies, 15 persons aged 9-12 years without maxillomandibular anomalies and acquired deformities (norm group) and 15 persons aged 9-12 years with hereditary syndromes, which are combined with bone deformities in the maxillofacial area (comparison group).

RESULTS

Results: Clinical examination showed that oral habits were manifested in 98.3% of patients. The results of clinical and radiological examination, analysis of cephalometric parameters and data on the thickness of the masticatory muscles on symmetrical areas of the face confirm the relationship between chronic oral habits and formation of acquired maxillomandibular anomalies; confirm the presence of acquired rather than congenital deformity of the facial skeleton, which is associated with changes in the thickness of the masticatory muscles on the part of the deformation та compensatory muscle hypertrophy on the opposite side.

CONCLUSION

Conclusions: The oral habit should be considered as one of the triggers in the development of acquired deformities of the maxillofacial area.

Structural and functional changes of nasal cavity and maxillary sinus in patients with skeletal class III malocclusion 1 year after bimaxillary surgery

PURPOSE

This study aimed to analyse changes in the nasal cavity and maxillary sinus structure and function in patients with skeletal class III malocclusion 1 year after bimaxillary surgery.

MATERIALS AND METHODS

In this study, cone-beam computed tomography (CBCT) images of 20 patients (10 men and 10 women; mean age 24.3 ± 3.4 years) with skeletal class III malocclusion who underwent Le Fort I osteotomy and bilateral sagittal split osteotomy were obtained before and 1 year after the surgery. CBCT data were stored opened with element 3D (E3D) to establish a nasal airway model (the paranasal sinus includes only the maxillary sinus). Ansys (ANSYS) software is used for simulation and analysis.

RESULTS

The maxillary sinus and nasal cavity volumes decreased significantly 1 year after the surgery. After surgery, the volume of nasal cavity decreased by 13.5%, and the average volume of maxillary sinus decreased by 7.8%. There was no significant difference in the degree of deviation of the septum and nasal cavity resistance, and air distribution in the maxillary sinus did not change. The nasal cavity wall shear stress change was similar to that before surgery.

CONCLUSIONS

The maxillary sinus volume and nasal cavity volume of patients with skeletal class III malocclusion changed significantly after bimaxillary surgery, but there was no significant change in nasal ventilation function 1 year after surgery.

Improvements in airflow characteristics and effect on the NOSE score after septoturbinoplasty: A computational fluid dynamics analysis

Septoturbinoplasty is a surgical procedure that can improve nasal congestion symptoms in patients with nasal septal deviation and inferior turbinate hypertrophy. However, it is unclear which physical domains of nasal airflow after septoturbinoplasty are related to symptomatic improvement. This work employs computational fluid dynamics modeling to identify the physical variables and domains associated with symptomatic improvement. Sixteen numerical models were generated using eight patients' pre- and postoperative computed tomography scans. Changes in unilateral nasal resistance, surface heat flux, relative humidity, and air temperature and their correlations with improvement in the Nasal Obstruction Symptom Evaluation (NOSE) score were analyzed. The NOSE score significantly improved after septoturbinoplasty, from 14.4 ± 3.6 to 4.0 ± 4.2 (p < 0.001). The surgery not only increased the airflow partition on the more obstructed side (MOS) from 31.6 ± 9.6 to 41.9 ± 4.7% (p = 0.043), but also reduced the unilateral nasal resistance in the MOS from 0.200 ± 0.095 to 0.066 ± 0.055 Pa/(mL·s) (p = 0.004). Improvement in the NOSE score correlated significantly with the reduction in unilateral nasal resistance in the preoperative MOS (r = 0.81). Also, improvement in the NOSE score correlated better with the increase in surface heat flux in the preoperative MOS region from the nasal valve to the choanae (r = 0.87) than in the vestibule area (r = 0.63). Therefore, unilateral nasal resistance and mucous cooling in the preoperative MOS can explain the perceived improvement in symptoms after septoturbinoplasty. Moreover, the physical domain between the nasal valve and the choanae might be more relevant to patient-reported patency than the vestibule area.

Septoplasty Training During the COVID-19 Era: Development and Validation of a Novel Low-Cost Simulation Model

OBJECTIVE

In a context of increasingly limited surgical exposition, enhanced by the coronavirus disease 2019 (COVID-19) pandemic context, the objective of this article is to explain the development of a novel low-cost and simple replication animal-based septoplasty training model for otolaryngology residents, to assess its face and construct validity, and to validate a specific rating scale for each task.

STUDY DESIGN

Experimental study.

SETTING

Surgical simulation laboratory.

METHODS

Septoplasty experts divided the procedure into key tasks. A simulator model to perform tasks was developed using pig ears to imitate human nasal septum cartilage, and a Specific Rating Scale was constructed. Trainees and faculty performed all tasks in the model. The participants were videotaped, and operative time, hand movements, and path length were recorded using a motion sensor device. Two blinded experts evaluated the videos with Global and Specific Rating Scales. All participants answered a satisfaction survey.

RESULTS

Fifteen subjects were recruited (7 trainees and 8 faculty). Significantly higher Global Rating Scale score, shorter operative time and path length, and fewer hand movements were observed in the faculty group. The satisfaction survey showed high applicability to a real scenario (mean score of 4.6 out of 5). Specific Rating Scale showed construct and concurrent validity and high reliability.

CONCLUSION

This simulation model and its specific rating scale can be accurately used as a validated surgical assessment tool for endonasal septoplasty skills. Its low cost and simple replicability make it a potentially useful tool in any otolaryngology surgical training program.

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.

Wet surface wall model for latent heat exchange during evaporation

Air conditioning is a dual heat and mass transfer process, and the human nasal cavity achieves this through the mucosal wall surface, which is supplied with an energy source through the sub-epithelial network of capillaries. Computational studies of air conditioning in the nasal cavity have included temperature and humidity, but most studies solved these flow parameters separately, and in some cases, a constant mucosal surface temperature was used. Recent developments demonstrated that both heat and mass transfer need to be modeled. This work expands on existing modeling efforts in accounting for the nasal cavity's dual heat and mass transfer process by introducing a new subwall model, given in the Supplementary Materials. The model was applied to a pipe geometry, and a human nasal cavity was recreated from CT-scans, and six inhalation conditions were studied. The results showed that when the energy transfer from the latent heat of evaporation is included, there is a cooling effect on the mucosal surface temperature.

Nasal airflow comparison in neonates, infant and adult nasal cavities using computational fluid dynamics

BACKGROUND AND OBJECTIVE

Neonates are preferential nasal breathers up to 3 months of age. The nasal anatomy in neonates and infants is at developing stages whereas the adult nasal cavities are fully grown which implies that the study of airflow dynamics in the neonates and infants are significant. In the present study, the nasal airways of the neonate, infant and adult are anatomically compared and their airflow patterns are investigated.

METHODS

Computational Fluid Dynamics (CFD) approach is used to simulate the airflow in a neonate, an infant and an adult in sedentary breathing conditions. The healthy CT scans are segmented using MIMICS 21.0 (Materialise, Ann arbor, MI). The patient-specific 3D airway models are analyzed for low Reynolds number flow using ANSYS FLUENT 2020 R2. The applicability of the Grid Convergence Index (GCI) for polyhedral mesh adopted in this work is also verified.

RESULTS

This study shows that the inferior meatus of neonates accounted for only 15% of the total airflow. This was in contrast to the infants and adults who experienced 49 and 31% of airflow at the inferior meatus region. Superior meatus experienced 25% of total flow which is more than normal for the neonate. The highest velocity of 1.8, 2.6 and 3.7 m/s was observed at the nasal valve region for neonates, infants and adults, respectively. The anterior portion of the nasal cavity experienced maximum wall shear stress with average values of 0.48, 0.25 and 0.58 Pa for the neonates, infants and adults.

CONCLUSIONS

The neonates have an underdeveloped nasal cavity which significantly affects their airway distribution. The absence of inferior meatus in the neonates has limited the flow through the inferior regions and resulted in uneven flow distribution.

From SARS-CoV-2 infection to COVID-19 morbidity: an in silico projection of virion flow rates to the lower airway via nasopharyngeal fluid boluses

Background: While the nasopharynx is initially the dominant upper airway infection site for SARS-CoV-2, the physiologic mechanism launching the infection at the lower airway is still not well-understood. Based on the rapidity of infection progression to the lungs, it has been hypothesized that the nasopharynx may be acting as the primary seeding zone for subsequent contamination of the lower airway via aspiration of virus-laden boluses of nasopharyngeal fluids. Methodology: To examine the plausibility of the aspiration-driven mechanism, we have computationally tracked the inhalation process in three anatomic airway reconstructions and have quantified the nasopharyngeal liquid volume transmitted to the lower airspace during each aspiration. Results: Extending the numerical trends on aspiration volume to earlier records on aspiration frequencies indicates a total aspirated nasopharyngeal liquid volume of 0.3 – 0.76 ml/day. Subsequently, for mean sputum viral load, our modeling projects that the number of virions reaching the lower airway will range over 2.1×106 – 5.3×106 /day; for peak viral load, the corresponding number hovers between 7.1×108 – 1.8×109. Conclusions: The virion transmission findings fill in a key piece of the mechanistic puzzle on the systemic progression of SARS-CoV-2, and subjectively point to health conditions like dysphagia, with proclivity to increased aspiration, as some of the potential underlying risk factors for aggressive lung infections.

In-silico investigation of airflow and micro-particle deposition in human nasal airway pre- and post-virtual transnasal sphenoidotomy surgery

Sphenoid sinus, located posterior to the nasal cavity, is difficult to reach for a surgery. Several operation procedures are available for sphenoidotomy, including endoscopic surgeries. Although the endoscopic sinus surgery is minimally invasive with low post-operative side effects, further optimization is required. Transnasal sphenoidotomy is a low invasive alternative to transethmoidal sphenoidotomy, but it still needs to be studied to understand its effects on the airflow pattern and the particle deposition. In this work, we simulated airflow and the micro-particle deposition in the nasal airway of a middle-aged man to investigate the change in particle deposition in the sphenoid sinus after virtual transnasal sphenoidotomy surgery. The results demonstrated that after transnasal sphenoidotomy, particle deposition in the targeted sphenoid sinus was an order of magnitude lower than that observed after virtual transethmoidal sphenoidotomy surgery. In addition, the diameter of the particles for the peak deposition fraction in the targeted sinus was shifted to smaller diameters after the transnasal sphenoidotomy surgery compared with that in the post-transethmoidal condition. These results suggest that the endoscopic transnasal sphenoidotomy can be a better procedure for sphenoid surgeries as it decreases the chance of bacterial contaminations and consequently lowers the surgical side effects and recovery time.

Computational Fluid Dynamics Modeling of Nasal Obstruction and Associations with Patient-Reported Outcomes

BACKGROUND

Nasal obstruction is a common problem, with significant impact on quality of life. Accurate diagnosis may be challenging because of the complex and dynamic nature of the involved anatomy. Computational fluid dynamics modeling has the ability to identify specific anatomical defects, allowing for a targeted surgical approach. The goal of the current study is to better understand nasal obstruction as it pertains to disease-specific quality of life by way of a novel computational fluid dynamics model of nasal airflow.

METHODS

Fifty-three patients with nasal obstruction underwent computational fluid dynamics modeling based on computed tomographic imaging. Nasal resistance was compared to demographic data and baseline subjective nasal patency based on Nasal Obstructive Symptom Evaluation scores.

RESULTS

Mean Nasal Obstructive Symptom Evaluation score among all patients was 72.6. Nasal Obstructive Symptom Evaluation score demonstrated a significant association with nasal resistance in patients with static obstruction (p = 0.03). There was a positive correlation between Nasal Obstructive Symptom Evaluation score and nasal resistance in patients with static bilateral nasal obstruction (R2 = 0.32) and poor correlation in patients with dynamic bilateral obstruction caused by nasal valve collapse (R2 = 0.02). Patients with moderate and severe bilateral symptoms had significantly higher nasal resistance compared to those with unilateral symptoms (p = 0.048).

CONCLUSIONS

Nasal obstruction is a multifactorial condition in most patients. This study shows correlation between simulated nasal resistance and Nasal Obstructive Symptom Evaluation score in a select group of patients. There is currently no standardized diagnostic algorithm or gold standard objective measure of nasal airflow; however, computational fluid dynamics may better inform treatment planning and surgical techniques on an individual basis.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Risk, V.

Heterogeneous distribution of mature olfactory sensory neurons in human olfactory epithelium

BACKGROUND

The olfactory cleft (OC) comprising the olfactory epithelium (OE) is the most important anatomical location for olfactory function. Endoscopic sinus surgery (ESS) is used to treat diseases related to the OC and improve olfactory dysfunction. However, iatrogenic OE injury occasionally occurs. Comprehensive knowledge of the olfactory region is required to avoid damage to the OE during endoscopic procedures.

METHODS

Immunohistochemistry was performed on olfactory mucosa obtained from the unaffected side of olfactory neuroblastoma surgical specimens. The OE was defined as the epithelium containing mature olfactory sensory neurons (OSNs). The distribution and cell kinetics of the OE were examined.

RESULTS

The OE was selectively localized to the anterior two-thirds of the superior turbinate (ST) and in the nasal septum (NS) just opposite to the ST; the OE was not detected within the mucosa of the superior meatus. The density of mature OSNs was high at the ethmoid tegmen but gradually decreased with distance from the ethmoid tegmen. The extent of cell death and proliferation was relatively even across the OE. Analysis of airflow profiles revealed that resection of inferior ST does not decrease airflow to the OC.

CONCLUSION

The results indicate that the distribution and degree of differentiation of mature OSNs are heterogenous throughout the OE. Epithelial resection of the anterior or superior ST has the potential to damage olfactory function. Resection of the inferior or posterior ST or widening of the superior meatus is a safer alternative that does not damage mature OSNs or alter airflow to the OC.

A numerical study of the effects of ambient temperature and humidity on the particle growth and deposition in the human airway

A numerical study was conducted on the effects of ambient temperature and humidity on the transportation of sodium chloride particles (100 nm-1 μm) in a human airway model ranging from the nasal cavity to bronchi. A mucus-tissue structure was adopted to model the mass and heat transfer on the airway surface boundary. The temperature and humidity distributions of the respiratory flow were calculated and then the interaction between the particle and water vapor was further analyzed. It was predicted that the particle size grew to the ratio of 5-6 under subsaturation conditions because of hygroscopicity, which shifted the deposition efficiency in opposite directions on dependence of the initial particle size. However, the particles could be drastically raised to 40 times of the initial 100 nm diameter if the supersaturation-induced condensation was established, that was prone to occur under the cold-dry condition, and consequently promoted the deposition significantly. Such behavior might effectively contribute to the revitalized coronavirus disease 2019 (COVID-19) pandemic in addition to the more active virus itself in winter.

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.

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.

The impact of nasal adhesions on airflow and mucosal cooling - A computational fluid dynamics analysis

Nasal adhesions are a known postoperative complication following surgical procedures for nasal airway obstruction (NAO); and are a common cause of surgical failure, with patients often reporting significant NAO, despite relatively minor adhesion size. Division of such nasal adhesions often provides much greater relief than anticipated, based on the minimal reduction in cross-sectional area associated with the adhesion. The available literature regarding nasal adhesions provides little evidence examining their quantitative and qualitative effects on nasal airflow using objective measures. This study examined the impact of nasal adhesions at various anatomical sites on nasal airflow and mucosal cooling using computational fluid dynamics (CFD). A high-resolution CT scan of the paranasal sinuses of a 25-year-old, healthy female patient was segmented to create a three-dimensional nasal airway model. Virtual nasal adhesions of 2.5 mm diameter were added to various locations within the nasal cavity, representing common sites seen following NAO surgery. A series of models with single adhesions were created. CFD analysis was performed on each model and compared with a baseline no-adhesion model, comparing airflow and heat and mass transfer. The nasal adhesions resulted in no significant change in bulk airflow patterns through the nasal cavity. However, significant changes were observed in local airflow and mucosal cooling around and immediately downstream to the nasal adhesions. These were most evident with anterior nasal adhesions at the internal valve and anterior inferior turbinate. Postoperative nasal adhesions create local airflow disruption, resulting in reduced local mucosal cooling on critical surfaces, explaining the exaggerated perception of nasal obstruction. In particular, anteriorly located adhesions created greater disruption to local airflow and mucosal cooling, explaining their associated greater subjective sensation of obstruction.

Nasal airflow of patient with septal deviation and allergy rhinitis

A numerical simulation of a patient’s nasal airflow was developed via computational fluid dynamics. Accordingly, computerized tomography scans of a patient with septal deviation and allergic rhinitis were obtained. The three-dimensional (3D) nasal model was designed using InVesalius 3.0, which was then imported to (computer aided 3D interactive application) CATIA V5 for modification, and finally to analysis system (ANSYS) flow oriented logistics upgrade for enterprise networks (FLUENT) to obtain the numerical solution. The velocity contours of the cross-sectional area were analyzed on four main surfaces: the vestibule, nasal valve, middle turbinate, and nasopharynx. The pressure and velocity characteristics were assessed at both laminar and turbulent mass flow rates for both the standardized and the patient’s model nasal cavity. The developed model of the patient is approximately half the size of the standardized model; hence, its velocity was approximately two times more than that of the standardized model.

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.

Computational characterization of inhaled droplet transport to the nasopharynx

How human respiratory physiology and the transport phenomena associated with inhaled airflow in the upper airway proceed to impact transmission of SARS-CoV-2, leading to the initial infection, stays an open question. An answer can help determine the susceptibility of an individual on exposure to a COVID-2019 carrier and can also provide a preliminary projection of the still-unknown infectious dose for the disease. Computational fluid mechanics enabled tracking of respiratory transport in medical imaging-based anatomic domains shows that the regional deposition of virus-laden inhaled droplets at the initial nasopharyngeal infection site peaks for the droplet size range of approximately 2.5-19 [Formula: see text]. Through integrating the numerical findings on inhaled transmission with sputum assessment data from hospitalized COVID-19 patients and earlier measurements of ejecta size distribution generated during regular speech, this study further reveals that the number of virions that may go on to establish the SARS-CoV-2 infection in a subject could merely be in the order of hundreds.

An Overview of Computational Fluid Dynamics Preoperative Analysis of the Nasal Airway

Evaluation of the nasal airway is crucial for every patient with symptoms of nasal obstruction as well as for every patient with other nasal symptoms. This assessment of the nasal airway comprises clinical examination together with imaging studies, with the correlation between findings of this evaluation and symptoms reported by the patient being based on the experience of the surgeon. Measuring nasal airway resistance or nasal airflow can provide additional data regarding the nasal airway, but the benefit of these objective measurements is limited due to their lack of correlation with patient-reported evaluation of nasal breathing. Computational fluid dynamics (CFD) has emerged as a valuable tool to assess the nasal airway, as it provides objective measurements that correlate with patient-reported evaluation of nasal breathing. CFD is able to evaluate nasal airflow and measure variables such as heat transfer or nasal wall shear stress, which seem to reflect the activity of the nasal trigeminal sensitive endings that provide sensation of nasal breathing. Furthermore, CFD has the unique capacity of making airway analysis of virtual surgery, predicting airflow changes after trial virtual modifications of the nasal airway. Thereby, CFD can assist the surgeon in deciding surgery and selecting the surgical techniques that better address the features of each specific nose. CFD has thus become a trend in nasal airflow assessment, providing reliable results that have been validated for analyzing airflow in the human nasal cavity. All these features make CFD analysis a mainstay in the armamentarium of the nasal surgeon. CFD analysis may become the gold standard for preoperative assessment of the nasal airway.

A numerical simulation of air flow in the human respiratory system for various environmental conditions

The functions of the nasal cavity are very important for maintaining the internal environment of the lungs since the inner walls of the nasal cavity control the temperature and saturation of the inhaled air with water vapor until the nasopharynx is reached. In this paper, three-dimensional computational studies of airflow transport in the models of the nasal cavity were carried out for the usual inspiratory velocity in various environmental conditions. Three-dimensional numerical results are compared with experimental data and calculations of other authors. Numerical results show that during normal breathing, the human nose copes with heat and relative moisture metabolism in order to balance the intra-alveolar conditions. It is also shown in this paper that a normal nose can maintain balance even in extreme conditions, for example, in cold and hot weather. The nasal cavity accelerates heat transfer by narrowing the air passages and swirls from the nasal concha walls of the inner cavity.

Nasal air conditioning following total inferior turbinectomy compared to inferior turbinoplasty - A computational fluid dynamics study

BACKGROUND

The aim of this study was to use computational fluid dynamics (CFD) to investigate the effects on nasal heat exchange and humidification of two different surgical techniques for reducing the inferior turbinate under different environmental conditions.

METHODS

Virtual surgery using two techniques of turbinate reduction was performed in eight nasal airway obstruction patients. Bilateral nasal airway models for each patient were compared: 1) Pre-operative 2) Post inferior turbinoplasty 3) Post total inferior turbinate resection (ITR). Two representative healthy models were included. Three different environmental conditions were investigated 1) ambient air 2) cold, dry air 3) hot, humid air. CFD modelling of airflow and conditioning was performed under steady-state, laminar, inspiratory conditions.

FINDINGS

Nasal conditioning is significantly altered following inferior turbinate reduction surgery, particularly with ITR under cold, dry inspired air (CDA). The degree of impairment is minor under the simulated range of environmental conditions (temperature = 12-40 °C; relative humidity = 13-80%). Streams of significantly colder air are found in the nasopharynx and more prevalent under CDA in ITR. These are related to high velocity flow streams, which remain cool in their centre throughout the widened inferior nasal cavity.

INTERPRETATION

Reduced air-mucosal heat exchange and moisture carrying capacity occurs under cooler temperatures in patients following inferior turbinate surgery. The clinical impact in extremely cold and dry conditions in groups with poor baseline respiratory function, respiratory illness, or endurance athletes is of special interest.

Effect of inhalation on oropharynx collapse via flow visualisation

Computational fluid dynamics (CFD) modelling has made significant contributions to the analysis and treatment of obstructive sleep apnoea (OSA). While several investigations have considered the flow field within the airway and its effect on airway collapse, the effect of breathing on the pharynx region is still poorly understood. We address this gap via a combined experimental and numerical study of the flow field within the pharynx and its impacts upon airway collapse. Two 3D experimental models of the upper airway were constructed based upon computerised tomography scans of a specific patient diagnosed with severe OSA; (i) a transparent, rigid model for flow visualisation, and (ii) a semi-flexible model for understanding the effect of flow on pharynx collapse. Validated simulation results for this geometry indicate that during inhalation, negative pressure (with respect to atmospheric pressure) caused by vortices drives significant narrowing of the pharynx. This narrowing is strongly dependent upon whether inhalation occurs through the nostrils. Thus, the methodology presented here can be used to improve OSA treatment by improving the design methodology for personalised, mandibular advancement splints (MAS) that minimise OSA during sleep.

Can computational fluid dynamic models help us in the treatment of chronic rhinosinusitis

PURPOSE OF REVIEW

The aim of this study was to review the recent literature (January 2017-July 2020) on computational fluid dynamics (CFD) studies relating to chronic rhinosinusitis (CRS), including airflow within the pre and postoperative sinonasal cavity, virtual surgery, topical drug and saline delivery (sprays, nebulizers and rinses) and olfaction.

RECENT FINDINGS

Novel CFD-specific parameters (heat flux and wall shear stress) are highly correlated with patient perception of nasal patency. Increased ostial size markedly improves sinus ventilation and drug delivery. New virtual surgery tools allow surgeons to optimize interventions. Sinus deposition of nasal sprays is more effective with smaller, low-inertia particles, outside of the range produced by many commercially available products. Saline irrigation effectiveness is improved using greater volume, with liquid entering sinuses via 'flooding' of ostia rather than direct jet entry.

SUMMARY

CFD has provided new insights into sinonasal airflow, air-conditioning function, the nasal cycle, novel measures of nasal patency and the impact of polyps and sinus surgery on olfaction. The deposition efficiency of topical medications on sinus mucosa can be markedly improved through parametric CFD experiments by optimising nasal spray particle size and velocity, nozzle angle and insertion location, while saline irrigation effectiveness can be optimized by modelling squeeze bottle volume and head position. More sophisticated CFD models (inhalation and exhalation, spray particle and saline irrigation) will increasingly provide translational benefits in the clinical management of CRS.

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.

Voxel-based simulation of flow and temperature in the human nasal cavity

The nasal airway is an extremely complex structure, therefore grid generation for numerical prediction of airflow in the nasal cavity is time-consuming. This paper describes the development of a voxel-based model with a Cartesian structured grid, which is characterized by robust and automatic grid generation, and the simulation of the airflow and air-conditioning in an individual human nasal airway. Computed tomography images of a healthy adult nose were used to reconstruct a virtual three-dimensional model of the nasal airway. Simulations of quiet restful inspiratory flow were then performed using a Neumann boundary condition for the energy equation to adequately resolve the flow and heat transfer. General agreements of airflow patterns, which were a high-speed jet posterior to the nasal valve and recirculating flow that occupied the anterior part of the upper cavity, and temperature distributions of the airflow and septum wall were confirmed by comparing in-vivo measurements with numerical simulation results.

Numerical study on the heat-recovery capacity of the human nasal cavity during expiration

The characteristics of the thermal field in the human nasal cavity during the expiration period were investigated using computational fluid dynamics. Heat and water-vapor recovery features were quantitatively investigated under realistic distributions of the epithelial surface and air temperature. A constant expiratory flow rate of 250 mL/s was assumed. The epithelial surface temperature was approximately 34.3-34.4 °C in the nasopharynx and 33.5-33.6 °C in the vestibule region, and these values are in good agreement with the measurement data in the literature. We observed that heat-recovery from the exhaled air mostly occurred in the posterior turbinate region, and the amount of heat recovered is estimated to be approximately 1/3 of the heat supply during inspiration. Because of this heat transfer from the exhaled air to the epithelial surface, the temperature of the epithelial surface increased in this region, and the exhaled air temperature dropped through the turbinate airway. Water-vapor recovery primarily occurs in the posterior segments of the turbinates; however, the amount of water-vapor transfer was approximately 1/5 of that in inspiration. Accordingly, the relative humidity of the exhaled air remained constant throughout the airway.

Effect of nasal airway nonlinearities on oscillometric resistance measurements in infants

Oscillometric measurements of respiratory system resistance (R_rs) in infants are usually made via the nasal pathways, which not only significantly contribute to overall R_rs but also introduce marked flow (V')-dependent changes. We employed intrabreath oscillometry in casts of the upper airways constructed from head CT images of 46 infants. We examined oscillometric nasal resistance (R_n) in upper airway casts with no respiratory flow (R₀) and the effect of varying V' on R_n by simulating tidal breathing. A characteristic nonlinear relationship was found between R_n and V', exhibiting segmental linearity and a prominent breakpoint (V'_bp) after log-log transformation. V'_bp was linearly related to the preceding value of end-expiratory volume acceleration (V″_eE; on average r² = 0.96, P < 0.001). R_n depended on V', and R at end-expiration (R_eE) showed a strong dependence on V″_eE in every cast (r² = 0.994, P < 001) with considerable interindividual variability. The intercept of the linear regression of R_eE versus V″_eE was found to be a close estimate of R₀. These findings were utilized in reanalyzed R_rs data acquired in vivo in a small group of infants (n = 15). Using a graphical method to estimate R₀ from R_eE, we found a relative contribution of V'-dependent nonlinearity to total resistance of up to 33%. In conclusion, we propose a method for correcting the acceleration-dependent nonlinearity error in R_eE. This correction can be adapted to estimate R₀ from a single intrabreath oscillometric measurement, which would reduce the masking effects of the upper airways on the changes in the intrathoracic resistance.NEW & NOTEWORTHY Oscillometric measurements of respiratory system resistance (R_rs) in infants are usually made via the nasal pathways, which not only significantly contribute to overall R_rs but also introduce marked flow acceleration-dependent distortions. Here, we propose a method for correcting flow acceleration-dependent nonlinearity error based on in vitro measurements in 3D-printed upper airway casts of infants as well as in vivo measurements. This correction can be adapted to estimate R_rs from a single intrabreath oscillometric measurement.

Three-dimensional analysis of sexual dimorphism in the soft tissue morphology of the upper airways in a human population

OBJECTIVES

Several studies have analyzed the sexual dimorphism of the skeletal cranial airways. This study aimed to quantify the three-dimensional (3D) morphology of the soft tissues of the upper airways in a human population. We addressed hypotheses about morphological features related to respiratory and energetic aspects of nasal sexual dimorphism.

METHODS

We reconstructed 3D models of 41 male and female soft tissue nasal airways from computed tomography data. We measured 280 landmarks and semilandmarks for 3D-geometric morphometric analyses to test for differences in size and 3D morphology of different functional compartments of the soft tissue airways.

RESULTS

We found statistical evidence for sexual dimorphism: Males were larger than females. 3D features indicated taller and wider inflow tracts, taller outflow tracts and slightly taller internal airways in males. These characteristics are compatible with greater airflow in males.

DISCUSSION

The differences in 3D nasal airway morphology are compatible with the respiratory-energetics hypothesis according to which males differ from females because of greater energetic demands. Accordingly, structures related to inflow and outflow of air show stronger signals than structures relevant for air-conditioning.

Assessment of PIV performance in validating CFD models from nasal cavity CBCT scans

OBJECTIVE

The aim of our study was to investigate how well Particle Image Velocimetry (PIV) measurements could serve Computational Fluid Dynamics (CFD) model validation for nasal airflow.

MATERIAL AND METHODS

For the PIV measurements, a silicone model of the nose based on cone beam computed tomography (CBCT) scans of a patient was made. Corresponding CFD calculations were conducted with laminar and two turbulent models (k-ω and k-ω SST).

RESULTS

CFD and PIV results corresponded well in our study. Especially, the correspondence of CFD calculations between the laminar and turbulent models was found to be even stronger. When comparing CFD with PIV, we found that the results were most convergent in the wider parts of the nasal cavities.

CONCLUSION

PIV measurements in realistically modelled nasal cavities succeed acceptably and CFD calculations produce corresponding results with PIV measurements. Greater model scaling is, however, necessary for better validations with PIV and comparisons of competing CFD models.

Use of anatomically-accurate 3-dimensional nasal airway models of adult human subjects in a novel methodology to identify and evaluate the internal nasal valve

The optimal method for radiographic evaluation of the internal nasal valve (INV) has not been established. The objective of this study was to develop a method to assess the cross-sectional area and the angle of the INV using anatomically-accurate 3D digital nasal airway models. Axial CT images of the paranasal sinuses of twenty adult subjects with healthy nasal airways (50% female and 50% age ≥ 50) were used to create the models. Patients with significant radiographic evidence of sinonasal disease were excluded. A primary cutting plane that passed through the edge of the nasal bone, upper lateral cartilage, and the head of the inferior turbinate was defined in coronal view. This primary coronal cutting plane was then rotated in 5° increments anteriorly while ensuring the anatomic criteria for the INV were still met. The cutting plane resulting in the minimum INV area was identified as the optimal cutting plane and the total cross-sectional area of INV in this plane,198.79 ± 54.57 mm², was significantly less than the areas obtained using the existing methods for radiographic evaluation of the INV. The angle between the optimal cutting plane and nasal dorsum was 75.00 ± 10.26°, and the corresponding INV angle was 10.77 ± 6.02°.

Characterization of the Airflow within an Average Geometry of the Healthy Human Nasal Cavity

This study's objective was the generation of a standardized geometry of the healthy nasal cavity. An average geometry of the healthy nasal cavity was generated using a statistical shape model based on 25 symptom-free subjects. Airflow within the average geometry and these geometries was calculated using fluid simulations. Integral measures of the nasal resistance, wall shear stresses (WSS) and velocities were calculated as well as cross-sectional areas (CSA). Furthermore, individual WSS and static pressure distributions were mapped onto the average geometry. The average geometry featured an overall more regular shape that resulted in less resistance, reduced WSS and velocities compared to the median of the 25 geometries. Spatial distributions of WSS and pressure of the average geometry agreed well compared to the average distributions of all individual geometries. The minimal CSA of the average geometry was larger than the median of all individual geometries (83.4 vs. 74.7 mm²). The airflow observed within the average geometry of the healthy nasal cavity did not equal the average airflow of the individual geometries. While differences observed for integral measures were notable, the calculated values for the average geometry lay within the distributions of the individual parameters. Spatially resolved parameters differed less prominently.