A computational analysis of nasal vestibule morphologic variabilities on nasal function

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.

Parameter characteristics in intranasal drug delivery: A key to targeting medications to the olfactory airspace

BACKGROUND

The nose is a viable pathway for topical drug delivery to the olfactory cleft for treatment of obstructive smell loss and nose-to-brain drug delivery. This study investigates how variations in nasal vestibule morphology influence intranasal spray drug transport to the olfactory cleft and olfactory roof/bulb regions.

METHODS

The unilateral nasal vestibule morphology in three healthy subjects with healthy normal nasal anatomy was classified as Elongated (Subject DN001), Notched (Subject DN002), and Standard (Subject DN003). Computational fluid and particle dynamics modelling were used to simulate nasal airflow and drug particle transport to the olfactory cleft and olfactory roof/bulb regions in each subject-specific nasal cavity. To evaluate highest drug depositions in these regions, the particle transport simulations involved extensive parameter combination analyses: 6 inspiratory flow rates mimicking resting to sniffing (10-50 L/min); 5 spray release locations (Top, Bottom, Central, Lateral, and Medial); 5 head positions (Upright, Tilted Forward, Tilted Back, Supine, and Mygind); 3 particle velocities (1, 5, and 10 m/s); 350,000 μm-particles (1-100 μm) and 346,500 nanoparticles (10-990 nm).

FINDINGS

Particle size groups with highest depositions in olfactory cleft: DN001 left = 28.4% at 11-20 μm, right = 75.3% at 6-10 μm; DN002 left = 16.8% at 1-5 μm, right = 45.3% at 30-40 nm; DN003 left = 29.1% at 21-30 μm, right = 15.9% at 6-10 μm. Highest depositions in olfactory roof/bulb: DN001 left = 6.5% at 11-20 μm, right = 26.4% at 11-20 μm; DN002 left = 3.6% at 1-5 μm, right = 2.6% at 1-5 μm; DN003 left = 2.8% at 21-30 μm, right = 1.7% at 31-40 μm.

INTERPRETATION

DN001 (Elongated nasal vestibule) had the most deposition in the olfactory regions. Micron-particles size groups generally had better deposition in the olfactory regions.

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.

Insights into exercise-induced rhinitis based on nasal aerodynamics induced by airway morphology

BACKGROUND

Exercise-induced rhinitis (EIR) is a poorly understood phenomenon that may be related to increased inspiratory airflow. Characterization of the development of EIR is important to understand contributing factors.

OBJECTIVE

To characterize how different nasal morphologies respond to airflow-related variables during rapid/deep inspiratory conditions.

METHODS

Subject-specific nasal airways were reconstructed from radiographic images. Unilateral airways were classified as Standard, Notched, or Elongated accord to their distinct nasal vestibule morphology. Computational fluid dynamics simulations were performed at various airflow rates.

RESULTS

For all simulated flow rates, average resistance at the nasal vestibule, airflow velocity and wall sheer stress were highest in Notched. Average mucosal heat flux was highest in Standard. Notched phenotypes showed lower mean percent increases from 10 L/min to 50 L/min in all computed variables.

CONCLUSION

Resistance values and airflow velocities depicted a more constricted nasal vestibule in the Notched phenotypes, while perception of nasal mucosal cooling (heat flux) favored the Standard phenotypes. Different nasal phenotypes may predispose to EIR.

Characterizing Olfactory Dysfunction in Patients with Unilateral Cleft Lip Nasal Deformities

Background: Unilateral cleft lip nasal deformity (uCLND) is associated with olfactory dysfunction, but the underlying etiology remains poorly understood. Objective: To investigate the etiology of uCLND-associated olfactory dysfunction using clinical, computational, and histologic assessments. Methods: Inclusion criteria: uCLND patients >16 years undergoing septorhinoplasty. Exclusion criteria: prior septoplasty or rhinoplasty, pregnancy, sinusitis. Measured outcomes: patient-reported scores, rhinomanometry, smell identification and threshold tests, computational fluid dynamics (CFD) airflow simulations, and histologic analysis of olfactory epithelium. Results: Five uCLND subjects were included: 18-23 years, three male and two female, four left-sided cleft and one right-sided cleft. All subjects reported moderate to severe nasal obstruction. Smell identification and threshold tests showed varying degrees of hyposmia. Nasal resistance was higher on the cleft side versus noncleft side measured by rhinomanometry (median 3.85 Pa-s/mL, interquartile range [IQR] = 21.96, versus 0.90 Pa-s/mL, IQR = 5.17) and CFD (median 1.04 Pa-s/mL, IQR = 0.94 vs. 0.11 Pa-s/mL, IQR = 0.12). Unilateral olfaction varied widely and was dependent on unilateral percentage olfactory airflow. Biopsies revealed intact olfactory neuroepithelium. Conclusions: uCLND-associated olfactory dysfunction appears to be primarily conductive in etiology and highly susceptible to variations in nasal anatomy. Clinical Trial Registration number: NCT04150783.

Exploring the influence of nasal vestibule structure on nasal obstruction using CFD and Machine Learning method

Motivated by clinical findings about the nasal vestibule, this study analyzes the aerodynamic characteristics of the nasal vestibule and attempt to determine anatomical features which have a large influence on airflow through a combination of Computational Fluid Dynamics (CFD) and machine learning method. Firstly, the aerodynamic characteristics of the nasal vestibule are detailedly analyzed using the CFD method. Based on CFD simulation results, we divide the nasal vestibule into two types with distinctly different airflow patterns, which is consistent with clinical findings. Secondly, we explore the relationship between anatomical features and aerodynamic characteristics by developing a novel machine learning model which could predict airflow patterns based on several anatomical features. Feature mining is performed to determine the anatomical feature which has the greatest impact on respiratory function. The method is developed and validated on 41 unilateral nasal vestibules from 26 patients with nasal obstruction. The correctness of the CFD analysis and the developed model is verified by comparing them with clinical findings.

Gender Differences in Nasal Anatomy and Function Among Caucasians

Background: Anatomical variations influence nasal physiology, yet sex differences in physiology remains unclear. Objective: To investigate sex differences among Caucasians using computational fluid dynamics. Methods: Adult subjects were selected with normal nasal cone beam computed tomography (CBCT) images and Nasal Obstruction Symptom Evaluation scores ≤30. The CBCT images were used to create subject-specific airway models. Nasal surface area (SA) and volume were computed, and airflow and heat transfer were simulated. Results: The CBCT scans were taken from 23 females and 12 males. The SA and volume (males: mean = 25.0 cm3; females: mean = 19.5 cm3; p < 0.001; Cohen's d = 1.51) were significantly larger for males, but SA-to-volume ratio did not differ significantly. Although unilateral nasal resistance did not vary greatly, females had higher bilateral resistance (males: mean = 0.04 Pa.s/mL; females: mean = 0.05 Pa.s/mL; p = 0.044; Cohen's d = 0.37). Females had higher heat flux (males: mean = 158.5 W/m2; females: mean = 191.8 W/m2; p = 0.012; Cohen's d = 0.79), but males had larger SA where mucosal heat flux exceeds 50 W/m2. Conclusions: These findings suggest differences in normal nasal anatomy and physiology between Caucasian males and females, which may be useful when assessing sex-specific functional outcomes after nasal surgery.

Detailed comparative analysis of environmental microparticle deposition characteristics between human and monkey nasal cavities using a surface mapping technique

Inhaled particulate matter is associated with nasal diseases such as allergic rhinitis, rhinosinusitis and neural disorders. Its health risks on humans are usually evaluated by measurements on monkeys as they share close phylogenetic relationship. However, the reliability of cross-species toxicological extrapolation is in doubt due to physiological and anatomical variations, which greatly undermine the reliability of these expensive human surrogate models. This study numerically investigated in-depth microparticle transport and deposition characteristics on human and monkey (Macaca fuscata) nasal cavities that were reconstructed from CT-images. Deposition characteristics of 1-30μm particles were investigated under resting and active breathing conditions. Similar trends were observed for total deposition efficiencies and a single correlation using Stokes Number was fitted for both species and both breathing conditions, which is convenient for monkey-human extrapolation. Regional deposition patterns were carefully compared using the surface mapping technique. Deposition patterns of low, medium and high inertial particles, classified based on their total deposition efficiencies, were further analyzed in the 3D view and the mapped 2D view, which allows locating particle depositions on specific nasal regions. According to the particle intensity contours and regional deposition profiles, the major differences were observed at the vestibule and the floor of the nasal cavity, where higher deposition intensities of medium and high inertial particles were shown in the monkey case than the human case. Comparisons of airflow streamlines indicated that the cross-species variations of microparticle deposition patterns are mainly contributed by two factors. First, the more oblique directions of monkey nostrils result in a sharper airflow turn in the vestibule region. Second, the monkey's relatively narrower nasal valves lead to higher impaction of medium and high inertial particles on the nasal cavity floor. The methods and findings in this study would contribute to an improved cross-species toxicological extrapolation between human and monkey nasal cavities.

Analyses on the influence of normal nasal morphological variations on odorant transport to the olfactory cleft

OBJECTIVE

Olfaction requires a combination of sensorineural components and conductive components, but conductive mechanisms have not typically received much attention. This study investigates the role of normal nasal vestibule morphological variations in ten healthy subjects on odorant flux in the olfactory cleft.

MATERIALS AND METHODS

Computed tomography images were used to create subject-specific nasal models. Each subject's unilateral nasal cavity was classified according to its nasal vestibule shape as Standard or Notched. Inspiratory airflow simulations were performed at 15 L/min, simulating resting inspiration using computational fluid dynamics modeling. Odorant transport simulations for three odorants (limonene, 2,4-dinitrotoluene, and acetaldehyde) were then performed at concentrations of 200 ppm for limonene and acetaldehyde, and 0.2 ppm for dinitrotoluene. Olfactory cleft odorant flux was computed for each simulation.

RESULTS AND DISCUSSION AND CONCLUSION

Simulated results showed airflow in the olfactory cleft was greater in the Standard phenotype compared to the Notched phenotype. For Standard, median airflow was greatest in the anterior region (0.5006 L/min) and lowest in the posterior region (0.1009 L/min). Median airflow in Notched was greatest in the medial region (0.3267 L/min) and lowest in the posterior region (0.0756 L/min). Median olfactory odorant flux for acetaldehyde and limonene was greater in Standard (Acetaldehyde: Standard = 140.45 pg/cm2-s; Notched = 122.20 pg/cm2-s. Limonene: Standard = 0.67 pg/cm2-s; Notched = 0.65 pg/cm2-s). Median dinitrotoluene flux was greater in Notched (Standard = 2.86 × 10-4pg/cm2-s; Notched = 4.29 × 10-4 pg/cm2-s). The impact of nasal vestibule morphological variations on odorant flux at the olfactory cleft may have implications on individual differences in olfaction, which should be investigated further.

The role of normal nasal morphological variations from race and gender differences on respiratory physiology

This study identifies anatomical and airflow-induced relationships based on nasal morphological variations due to inter- and intra-racial differences and gender. Subject-specific nasal airway reconstruction was created from computed tomography images in 16 subjects: 4 subjects from each ethnic group (Black, East Asian, Caucasian, and Latino) comprising of 2 males and 2 females. Volume, surface area and nasal index were calculated, as well as airflow rate and nasal resistance after computational fluid dynamics simulations in the nasal airway. Results showed that nasal airspace surface area (p = 0.0499) and volume (p = 0.0281) were significantly greater in males than in females. Nasal volume was greatest in East Asians (Median = 20.38cm³, Interquartile Range [IQR] = 4.58 cm³), Latinos had the greatest surface area (Median = 219.70cm², IQR = 29.56cm²). On average, East Asian and Black females had larger nasal index than their male counterparts. Caucasians had the highest median nasal resistance (0.050 Pa.s/mL, IQR = 0.025 Pa.s/mL). Results indicate that there exist anatomical variabilities based on race and gender. However, these variabilities may not significantly influence nasal function.

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.

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.

Analysis of nasal air conditioning in subjects with unilateral cleft lip nasal deformity

This study evaluated the impact of unilateral cleft lip nasal deformity (uCLND) on the ability of the nasal passages to warm and humidify inspired environmental air using computational fluid dynamics (CFD) modeling. Nasal air conditioning was simulated at resting inspiration in ten individuals with uCLND and seven individuals with normal anatomy. The overall heat and water transfer through nasal mucosa was significantly greater (p = 0.02 for both heat and moisture fluxes) on the non-cleft side than on the cleft side. Unilateral median and interquartile range (IQR) for heat flux (W/m²) was 190.3 (IQR 59.9) on the non-cleft side, 160.9 (IQR 105.0) on the cleft side, and 170.7 (IQR 87.8) for normal subjects. For moisture flux (mg/(s·m²), they were 357.4 (IQR 112.9), 298.7 (IQR 200.3) and 320.8 (IQR 173.0), respectively. Significant differences of SAHF50 between cleft side of uCLND and normal existed except for anterior region. Nevertheless, air conditioning ability in subjects with uCLND was generally comparable to that of normal subjects.

Orally Inhaled Drug Particle Transport in Computerized Models of Laryngotracheal Stenosis

OBJECTIVE

Adjuvant management for laryngotracheal stenosis (LTS) may involve inhaled corticosteroids, but metered dose inhalers are designed for pulmonary drug delivery. Comprehensive analyses of drug particle deposition efficiency for orally inhaled corticosteroids in the stenosis of LTS subjects are lacking.

STUDY DESIGN

Descriptive research.

SETTING

Academic medical center.

METHODS

Anatomically realistic 3-dimensional reconstructions of the upper airway were created from computed tomography images of 4 LTS subjects-2 subglottic stenosis and 2 tracheal stenosis subjects. Computational fluid dynamics modeling was used to simulate airflow and drug particle transport in each airway. Three inhalation pressures were simulated, 10 Pa, 25 Pa, and 40 Pa. Drug particle transport was simulated for 100 to 950 nanoparticles and 1 to 50 micron-particles. Particles were released into the airway to mimic varying inhaler conditions with and without a spacer chamber.

RESULTS

Based on smallest to largest cross-sectional area ratio, the laryngotracheal stenotic segment shrunk by 57% and 47%, respectively, for subglottic stenosis models and by 53% for both tracheal stenosis models. Airflow resistance at the stenotic segment was lower in subglottic stenosis models than in tracheal stenosis models: 0.001 to 0.011 Pa.s/mL vs 0.024 to 0.082 Pa.s/mL. Drug depositions for micron-particles and nanoparticles at stenosis were 0.06% to 2.48% and 0.10% to 2.60% for subglottic stenosis and tracheal stenosis models, respectively. Particle sizes with highest stenotic deposition were 6 to 20 µm for subglottic stenosis models and 1 to 10 µm for tracheal stenosis models.

CONCLUSION

This study suggests that at most, 2.60% of inhaled drug particles deposit at the stenosis. Particle size ranges with highest stenotic deposition may not represent typical sizes emitted by inhalers.

Role of nasal vestibule morphological variations on olfactory airflow dynamics

BACKGROUND

The conductive mechanisms of olfaction are typically given little priority in the evaluation of olfactory function. The objective of this study is to investigate the role of nasal vestibule morphological variations on airflow volume at the olfactory recess in healthy subjects.

METHODS

Anatomically realistic three-dimensional nasal airway models were constructed from computed tomography scans in five subjects. Each individual's unilateral nasal cavity (10 total) was classified according to the shape of their nasal vestibule: Standard, Notched, or Elongated. Nasal airflow simulations were performed using computational fluid dynamics modeling at two inspiratory flow rates (15 L/min and 30 L/min) to reflect resting and moderate breathing rates. Olfactory airflow volume and cross-sectional flow resistance were computed.

FINDINGS

Average olfactory airflow volumes (and percent airflow in olfactory) were: 0.25 L/min to 0.64 L/min (3.0%-7.7%; 15 L/min simulations) and 0.53 L/min to 1.30 L/min (3.2%-7.8%; 30 L/min simulations) for Standard; 0.13 L/min - 0.47 L/min (2.0%-6.8%; 15 L/min simulations) and 0.06 L/min - 0.82 L/min (1.7%-6.1%; 30 L/min simulations) for Notched; and 0.07 L/min - 0.39 L/min (1.2%-5.4%; 15 L/min simulations) and 0.30 L/min - 0.99 L/min (2.1%-6.7%; 30 L/min simulations) for Elongated. On average, relative difference in olfactory resistance between left and right sides was 141.5% for patients with different unilateral phenotypes and 82.2% for patients with identical unilateral phenotype.

INTERPRETATION

Olfactory cleft airflow volume was highest in the Standard nasal vestibule phenotype, followed by Notched phenotype for 15 L/min simulations and Elongated phenotype for 30 L/min simulations. Further, intra-patient variation in olfactory cleft airflow resistance differs greatly for patients with different unilateral phenotypes compared to patients with identical unilateral phenotype.

Details of the physiology of the aerodynamic and heat and moisture transfer in the normal nasal cavity

Anatomically accurate 3D models of 10 healthy nasal cavities are developed from computerized tomography (CT) scan images. Considering anatomical and physiological importance of different parts of the nasal cavity, the surface of each nasal passage is divided to eleven anatomical surfaces. Also the coronal cross sections in the nasal passage are divided to six sub-sections that share the total nasal passage airflow. The details of the flow field, heat transfer and water-vapor transport are numerically investigated for resting and low activity conditions. The mean and standard deviation of the different anatomical and air conditioning parameters such as: surface area, wall shear stress, heat and moisture transfer on different parts of the nasal passage surfaces and volume flow rates through different sections are presented. Results show that the percentages of airflow for inferior, middle and superior meatuses are 11.3 ± 6.4, 36.5 ± 9.5, 1.9 ± 0.81 % respectively and 4.1 ± 2.1 % of air passes through olfactory area. The inhaled air passing from the remaining surface (main passage) is 46.2 ± 10 %. Heat and moisture fluxes are highest in the anterior part of the nasal cavity, turbinates and lower part of the septum respectively. The percentage of the heat transfer from turbinates is 25.7 ± 3.9 % of total nasal heat transfer.

Olfactory Targeting of Microparticles Through Inhalation and Bi-directional Airflow: Effect of Particle Size and Nasal Anatomy

Background: Targeting drugs to the olfactory region in the nasal cavity can bypass the restrictive blood-brain barrier and enhance their direct delivery to the brain. However, complex nasal geometry and its demographical variations can pose challenges for targeted drug deposition in the olfactory region. Deposition of particles in the nasal cavity is influenced by particle size, airflow rate, and nasal geometry. Therefore, this study investigated the effect of these parameters on regional microparticle deposition with the view to provide insights into the nose-to-brain delivery of drugs. Methods: In this study, three anatomically accurate human nasal cavities were reconstructed in silico and deposition of microparticles under nebulization and bi-directional airflow conditions was simulated. Microparticle deposition data were analyzed to gain insight into the effect of particle size and nasal geometry. Results: Maximum olfactory deposition was observed with particles in the size range of 8 to 12 μm under nebulization and 14 to 18 μm under bi-directional airflow condition. Geometric differences between subjects were shown to significantly impact overall and regional particle deposition and introduced inter-subject variability. Significant intra-subject variability in microparticle deposition was also observed in the bi-directional delivery cases. Conclusions: The data from this study suggest that tailoring particle size, combined with a delivery protocol, may provide a unique and pragmatic way to target drugs to the olfactory region. Differences in nasal anatomy among humans can cause variability in particle deposition and need to be considered in any future applications.

In silico approaches to respiratory nasal flows: A review

The engineering discipline of in silico fluid dynamics delivers quantitative information on airflow behaviour in the nasal regions with unprecedented detail, often beyond the reach of traditional experiments. The ability to provide visualisation and analysis of flow properties such as velocity and pressure fields, as well as wall shear stress, dynamically during the respiratory cycle may give significant insight to clinicians. Yet, there remains ongoing challenges to advance the state-of-the-art further, including for example the lack of comprehensive CFD modelling on varied cohorts of patients. The present article embodies a review of previous and current in silico approaches to simulating nasal airflows. The review discusses specific modelling techniques required to accommodate physiologically- and clinically-relevant findings. It also provides a critical summary of the reported results in the literature followed by an outlook on the challenges and topics anticipated to drive research into the future.

Normative ranges of nasal airflow variables in healthy adults

PURPOSE

Virtual surgery planning based on computational fluid dynamics (CFD) simulations of nasal airflow has the potential to improve surgical outcomes for patients with nasal airway obstruction (NAO). Virtual surgery planning requires normative ranges of airflow variables, but few studies to date have quantified inter-individual variability of nasal airflow among healthy subjects. This study reports CFD simulations of nasal airflow in 47 healthy adults.

METHODS

Anatomically accurate three-dimensional nasal models were reconstructed from cone beam computed tomography scans and used for steady-state inspiratory airflow simulations with a bilateral flowrate of 250 ml/s. Normal subjective sensation of nasal patency was confirmed using the nasal obstruction symptom evaluation and visual analog scale. Healthy ranges for several CFD variables known to correlate with subjective nasal patency were computed, including unilateral airflow, nasal resistance, airspace minimal cross-sectional area (mCSA), heat flux (HF), and surface area stimulated by mucosal cooling (defined as the area where HF > 50 W/m²). The normative ranges were targeted to contain 95% of the healthy population and computed using a nonparametric method based on order statistics.

RESULTS

A wide range of inter-individual variability in nasal airflow was observed among healthy subjects. Unilateral airflow varied from 60 to 191 ml/s, airflow partitioning ranged from 23.8 to 76.2%, and unilateral mCSA varied from 0.24 to 1.21 cm². These ranges are in good agreement with rhinomanometry and acoustic rhinometry data from the literature. A key innovation of this study are the normative ranges of flow variables associated with mucosal cooling, which recent research suggests is the primary physiological mechanism of nasal airflow sensation. Unilateral HF ranged from 94 to 281 W/m², while the surface area stimulated by cooling ranged from 27.4 to 64.3 cm².

CONCLUSIONS

These normative ranges may serve as targets in future virtual surgery planning for patients with NAO.

Geometry and airflow dynamics analysis in the nasal cavity during inhalation

BACKGROUND

A major issue among computational respiratory studies is the wide variety of nasal morphologies being studied, caused by both inter-population and inter-subject variations.

METHOD

Six nasal cavity geometries exhibiting diverse geometry variations were subjected to steady inhalation flow rate of 15L/min. to determine if any consistent flow behaviour could be found.

FINDINGS

Despite vastly different geometries we were able to identify consistent flow patterns including relatively high velocity in the nasal valve region, followed by flow continuing predominantly in the inferior half of the airway. We also found conformity among models where the inhaled air reached a near-conditioned state by the middle of the nasal cavity. Air from the front of the face reached the olfactory regions while air from the lateral sides of the face moved through the inferior half of the nasal cavity.

INTERPRETATION

The ability to predict gross flow features provides a baseline flow field to compare against. This contributes towards establishing well defined flow predictions and be used as a comparison for future larger studies.

Nasal Structural and Aerodynamic Features That May Benefit Normal Olfactory Sensitivity

Nasal airflow that effectively transports ambient odors to the olfactory receptors is important for human olfaction. Yet, the impact of nasal anatomical variations on airflow pattern and olfactory function is not fully understood. In this study, 22 healthy volunteers were recruited and underwent computed tomographic scans for computational simulations of nasal airflow patterns. Unilateral odor detection thresholds (ODT) to l-carvone, phenylethyl alcohol (PEA) and d-limonene were also obtained for all participants. Significant normative variations in both nasal anatomy and aerodynamics were found. The most prominent was the formation of an anterior dorsal airflow vortex in some but not all subjects, with the vortex size being significantly correlated with ODT of l-carvone (r = 0.31, P < 0.05). The formation of the vortex is likely the result of anterior nasal morphology, with the vortex size varying significantly with the nasal index (ratio of the width and height of external nose, r = -0.59, P < 0.001) and nasal vestibule "notch" index (r = 0.76, P < 0.001). The "notch" is a narrowing of the upper nasal vestibule cartilage region. The degree of the notch also significantly correlates with ODT for PEA (r = 0.32, P < 0.05) and l-carvone (r = 0.33, P < 0.05). ODT of d-limonene, a low mucosal soluble odor, does not correlate with any of the anatomical or aerodynamic variables. The current study revealed that nasal anatomy and aerodynamics might have a significant impact on normal olfactory sensitivity, with greater airflow vortex and a narrower vestibule region likely intensifying the airflow vortex toward the olfactory region and resulting in greater olfactory sensitivity to high mucosal soluble odors.

The navigational nose: a new hypothesis for the function of the human external pyramid

One of the outstanding questions in evolution is why Homo erectus became the first primate species to evolve the external pyramid, i.e. an external nose. The accepted hypothesis for this trait has been its role in respiration, to warm and humidify air as it is inspired. However, new studies testing the key assumptions of the conditioning hypothesis, such as the importance of turbulence to enhance heat and moisture exchange, have called this hypothesis into question. The human nose has two functions, however, respiration and olfaction. It is thus also possible that the external nose evolved in response to selection for olfaction. The genus Homo had many adaptations for long-distance locomotion, which allowed Homo erectus to greatly expand its species range, from Africa to Asia. Long-distance navigation in birds and other species is often accomplished by orientation to environmental odors. Such olfactory navigation, in turn, is enhanced by stereo olfaction, made possible by the separation of the olfactory sensors. By these principles, the human external nose could have evolved to separate olfactory inputs to enhance stereo olfaction. This could also explain why nose shape later became so variable: as humans became more sedentary in the Neolithic, a decreasing need for long-distance movements could have been replaced by selection for other olfactory functions, such as detecting disease, that would have been critical to survival in newly dense human settlements.

High quality statistical shape modelling of the human nasal cavity and applications

The human nose is a complex organ that shows large morphological variations and has many important functions. However, the relation between shape and function is not yet fully understood. In this work, we present a high quality statistical shape model of the human nose based on clinical CT data of 46 patients. A technique based on cylindrical parametrization was used to create a correspondence between the nasal shapes of the population. Applying principal component analysis on these corresponded nasal cavities resulted in an average nasal geometry and geometrical variations, known as principal components, present in the population with a high precision. The analysis led to 46 principal components, which account for 95% of the total geometrical variation captured. These variations are first discussed qualitatively, and the effect on the average nasal shape of the first five principal components is visualized. Hereafter, by using this statistical shape model, two application examples that lead to quantitative data are shown: nasal shape in function of age and gender, and a morphometric analysis of different anatomical regions. Shape models, as the one presented here, can help to get a better understanding of nasal shape and variation, and their relationship with demographic data.

Air conditioning analysis among human nasal passages with anterior anatomical variations

A major functional role of the nasal cavity is air conditioning of the inspired environmental air to near alveolar conditions. It is well known that the anatomical disparities among nasal passages can change airflow patterns to a great extent. However, its effect on nasal air conditioning performance remains largely unexplored. This research investigated the nasal air conditioning performance among nasal models with distinct vestibule phenotypes, including subjects with and without vestibule notches. For the mass transfer, we used a two-film theory model to determine the species transport. Airflow patterns, heat and mass transfer between the inhaled airflow and the nasal mucosa were analysed and compared. Results showed that the nasal air conditioning performance is closely related to nasal passage structures. The anatomical variations, especially the geometry changes in the anterior vestibule region, can increase both heat and mass transfer rate between nasal mucous and respiratory air at the vicinity of the notched regions, while for other regions such as the anterior superior nasal cavity, the heat transfer is greatly reduced to even zero heat flux due to lack of active airflow passing.

Investigating the effects of laryngotracheal stenosis on upper airway aerodynamics

OBJECTIVE

Very little is known about the impact of laryngotracheal stenosis (LTS) on inspiratory airflow and resistance, especially in air hunger states. This study investigates the effect of LTS on airway resistance and volumetric flow across three different inspiratory pressures.

METHODS

Head-and-neck computed tomography scans of 11 subjects from 2010 to 2016 were collected. Three-dimensional reconstructions of the upper airway from the nostrils to carina, including the oral cavity, were created for one subject with a normal airway and for 10 patients with LTS. Airflow simulations were conducted using computational fluid dynamics modeling at three different inspiratory pressures (10, 25, 40 pascals [Pa]) for all subjects under two scenarios: 1) inspiration through nostrils only (MC), and 2) through both nostrils and mouth (MO).

RESULTS

Volumetric flows in the normal subject at the three inspiratory pressures were considerably higher (MC: 11.8-26.1 L/min; MO: 17.2-36.9 L/min) compared to those in LTS (MC: 2.86-6.75 L/min; MO: 4.11-9.00 L/min). Airway resistances in the normal subject were 0.051 to 0.092 pascal seconds per milliliter (Pa.s)/mL (MC) and 0.035-0.065 Pa.s/mL (MO), which were approximately tenfold lower than those of subjects with LTS: 0.39 to 0.89 Pa.s/mL (MC) and 0.45 to 0.84 Pa.s/mL (MO). Furthermore, subjects with glottic stenosis had the greatest resistance, whereas subjects with subglottic stenosis had the greatest variability in resistance. Subjects with tracheal stenosis had the lowest resistance.

CONCLUSION

This pilot study demonstrates that LTS increases resistance and decreases airflow. Mouth breathing significantly improved airflow and resistance but cannot completely compensate for the effects of stenosis. Furthermore, location of stenosis appears to modulate the effect of the stenosis on resistance differentially.

LEVEL OF EVIDENCE

NA. Laryngoscope, 128:E141-E149, 2018.

Nasal dilation effects on olfactory deposition in unilateral and bi-directional deliveries: In vitro tests and numerical modeling

The human nose can expand either actively or passively to increase airflow. Nasal dilation may alter drug delivery efficiencies in the nasal airway or olfactory region. However, the dosage enhancement from nasal dilations has not been quantified. The mechanisms underlying the dilation-induced deposition variation are also not clear. This study aims to quantify the nasal dilation effects on drug delivery in the nasal airway and olfactory region using in vitro tests and numerical analysis. Two variants of an existing normal nasal airway model were developed with different levels of airway dilation. Airway dimensions were quantified in terms of hydraulic diameter, cross-sectional area, and surface area to volume ratio. Sectional nose casts were prepared using a 3-D printer for visualizing deposition patterns and quantifying delivered dosages. A well-validated computational fluid-particle dynamics (CFPD) model was utilized to understand the underlying mechanisms in the unilateral and bi-directional deliveries. In vitro tests show that nasal dilation lowered the total dosage in the nose but increased the dosage to the olfactory region in both the unilateral and bi-directional deliveries. Compared to the normal nose with unilateral delivery, nasal dilation enhanced the olfactory deposition by a factor of 2.2, while nasal dilatation with the bi-directional delivery increased by a factor of 4. Complementary numerical analyses revealed the growth of a recirculation zone in the middle meatus of dilated noses, which induced lower pressure and increased ventilation to the upper nose. In bi-directional deliveries, a significantly higher fraction of airflow was ventilated to the upper airway in the outflow side of the nose and contributed to the elevated olfactory dosage. Nasal dilation in combination with the bi-directional delivery is recommended over the conventional unilateral method for olfactory targeting.

Characterizing human nasal airflow physiologic variables by nasal index

Although variations in nasal index (NI) have been reported to represent adaptation to climatic conditions, assessments of NI with airflow variables have not been rigorously investigated. This study uses computational fluid dynamics modeling to investigate the relationship between NI and airflow variables in 16 subjects with normal nasal anatomy. Airflow simulations were conducted under constant inspiratory pressure. Nasal resistance (NR) against NI showed weak association from nostrils to anterior inferior turbinate (R(2)=0.26) and nostril to choanae (R(2)=0.12). NI accounted for 38% and 41% of the respective variation in wall shear stress (WSS) and heat flux (HF) at the nasal vestibule, and 52% and 49% of variability in WSS and HF across the entire nose. HF and WSS had strong correlation with NI<80, and weakly correlated with NI>80; these differences in HF and WSS for NI<80 and NI>80 were not statistically significant. Results suggest strong relationship between NI and both WSS and HF but not NR, particularly in subjects with NI<80.

The relationship between nasal resistance to airflow and the airspace minimal cross-sectional area

The relationship between nasal resistance (R) and airspace minimal cross-sectional area (mCSA) remains unclear. After the introduction of acoustic rhinometry, many otolaryngologists believed that mCSA measurements would correlate with subjective perception of nasal airway obstruction (NAO), and thus could provide an objective measure of nasal patency to guide therapy. However, multiple studies reported a low correlation between mCSA and subjective nasal patency, and between mCSA and R. This apparent lack of correlation between nasal form and function has been a long-standing enigma in the field of rhinology. Here we propose that nasal resistance is described by the Bernoulli Obstruction Theory. This theory predicts two flow regimes. For mCSA>Acrit, the constriction is not too severe and there is not a tight coupling between R and mCSA. In contrast, when mCSAAcrit (estimated to be 0.37cm(2)), this theory suggests that airway constrictions are rarely an exclusive contributor to nasal resistance, which may explain the weak correlation between mCSA and subjective nasal patency.

Airflow in the Human Nasal Passage and Sinuses of Chronic Rhinosinusitis Subjects

Endoscopic surgery is performed on patients with chronic inflammatory disease of the paranasal sinuses to improve sinus ventilation. Little is known about how sinus surgery affects sinonasal airflow. In this study nasal passage geometry was reconstructed from computed tomographic imaging from healthy normal, pre-operative, and post-operative subjects. Transient air flow through the nasal passage during calm breathing was simulated. Subject-specific differences in ventilation of the nasal passage were observed. Velocity magnitude at ostium was different between left and right airway. In FESS, airflow in post-surgical subjects, airflow at the maxillary sinus ostium was upto ten times higher during inspiration. In a Lothrop procedure, airflow at the frontal sinus ostium can be upto four times higher during inspiration. In both post-operative subjects, airflow at ostium was not quasi-steady. The subject-specific effect (of surgery) on sinonasal interaction evaluated through airflow simulations may have important consequences for pre- and post-surgical assessment and surgical planning, and design for improvement of the delivery efficiency of nasal therapeutics.