Inflow boundary profile prescription for numerical simulation of nasal airflow

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.

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.

Reducing variability in nasal surgery outcomes through computational fluid dynamics and advanced 3D virtual surgery techniques

Objectives

This study aims to delineate the specific impact of using computational fluid dynamics (CFD) and 3D virtual surgery techniques in otolaryngology surgery, focusing on their roles in enhancing the precision of nasal surgery and optimizing future patient outcomes. The central objective was to assess whether these advanced technologies could reduce variability in surgical approaches and decision-making among specialists, thereby improving the consistency and efficacy of patient care in cases of nasal obstruction.

Methods and results

Our methodology involved a detailed analysis of pre- and post-operative scenarios using CFD feedback. Six otolaryngologists participated, employing virtual surgery techniques on two patients with diagnosed nasal obstruction. The CFD analysis focused on quantifying key airflow parameters: right nasal flow rate (QR), left nasal flow rate (QL), flow symmetry (Ф), and bilateral nasal resistance (R). These parameters were meticulously compared before and after the application of CFD feedback to evaluate changes in surgical planning and outcomes. Quantitative analysis revealed a notable decrease in the standard deviation of the measured parameters among the specialists post-CFD feedback, indicating reduced variability in surgical approaches. Specifically, for Patient #1 the standard deviation for QR values dropped from 0.694 L/min to 0.602 L/min, and for QL values from 0.676 L/min to 0.584 L/min, and for Patient #2, the standard deviation for QR values decreased from 2.204 L/min to 0.958 L/min, and for QL values from 2.295 L/min to 1.014 L/min. Moreover, the variability range, represented by the differences between the maximum and minimum values for Ф and R, diminished significantly. Post-operative average values for all parameters showed a convergence towards ideal basal levels, suggesting a more uniform and effective surgical strategy across different surgeons.

Conclusions

Both integration of CFD and 3D virtual surgery techniques in otolaryngology can substantially reduce variability in surgical planning and decision-making, ultimately leading to improved patient outcomes. These advanced tools have the potential to standardize the diagnosis and treatment of nasal pathologies, contributing to more effective and consistent care. Future research in this area should focus on larger patient cohorts and further exploration of the potential benefits and applications of CFD and virtual surgery in otolaryngology.

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.

Morphology of pig nasal structure and modulation of airflow and basic thermal conditioning

Mammals have presumably evolved to adapt to a diverse range of ambient environmental conditions through the optimized heat and mass exchange. One of the crucial biological structures for survivability is the nose, which efficiently transports and thermally preconditions the external air before reaching the internal body. Nasal mucosa and cavity help warm and humidify the inhaled air quickly. Despite its crucial role, the morphological features of mammal noses and their effect in modulating the momentum of the inhaled air, heat transfer dynamics, and particulate trapping remain poorly understood. Tortuosity of the nasal cavity in high-olfactory mammalian species, such as pigs and opossum, facilitates the formation of complex airflow patterns inside the nasal cavity, which leads to the screening of particulates from the inhaled air. We explored basic nasal features in anatomically realistic nasal pathways, including tortuosity, radius of curvature, and gap thickness; they show strong power-law correlations with body weight. Complementary inspection of tortuosity with idealized conduits reveals that this quantity is central in particle capture efficiency. Mechanistic insights into such nuances can serve as a tipping point to transforming nature-based designs into practical applications. In-depth characterization of the fluid-particle interactions in nasal cavities is necessary to uncover nose mechanistic functionalities. It is instrumental in developing new devices and filters in a number of engineering processes.

Machine learning and sensitivity analysis for predicting nasal drug delivery for targeted deposition

Targeted nasal drug delivery can provide improved efficacy for drug formulations to be delivered at high efficacy rates. Some parameters that influence drug delivery have a dependency on the patient's technique of administration and the spray device itself. When the different parameters, each having a specific range of values are combined, the combinatory permutations for studying its effects on particle deposition become large. In this study, we combine six input spray parameters (the spray half-cone angle, the mean spray exit velocity, the breakup length from the nozzle exit, the diameter of the nozzle spray device, the particle size, and the sagittal angle of the spray) with a range of values to produce 384 combinations of spray characteristics. This was repeated for three inhalation flow rates of 20, 40, and 60 L/min. To reduce the computational costs of a full transient Large Eddy Simulation flow field, we create a time-averaged frozen field and perform the time integration of particle trajectories through the flow field to determine the particle deposition in four anatomical regions of the nasal cavity (anterior, middle, olfactory and posterior) for each of the 384 spray field. A sensitivity analysis determined the significance of each input variable on the deposition. It was found the particle size distribution significantly affected deposition in the olfactory and posterior regions, while the spray device insertion angle was significant for deposition in the anterior and middle regions. Five machine learning models were evaluated based on 384 cases and it was found that despite the small sample dataset the simulation data was sufficient to provide accurate machine-learning predictions.

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.

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 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.

Three-dimensional form and function of the nasal cavity and nasopharynx in humans and chimpanzees

The facial differences between recent Pan troglodytes and Homo sapiens can be used as a proxy for the reduction of facial prognathism that happened during evolutionary transition between Australopithecines and early Homo. The projecting nasal morphology of Homo has been considered both a passive consequence of anatomical reorganization related to brain and integrated craniofacial evolution as well as an adaptation related to air-conditioning during physiological and behavioral shifts in human evolution. Yet, previous research suggested impaired air-conditioning in Homo challenging respiratory adaptations based on computational fluid dynamics (CFD) and airflow simulations. Here we improved CFD model at the inflow region and also carried out three-dimensional (3D) geometric morphometrics to address the hypothesis of impaired air-conditioning in humans and species differences in airway shape. With the new CFD model we simulated pressure, velocity, and temperature changes in airflow of six adult humans and six chimpanzees and analyzed 164 semi-landmarks of 10 humans and 10 chimpanzees for 3D size and shape comparisons. Our finding shows significantly different internal 3D nasal airways. Also, species means of pressure, velocity, and temperature differed statistically significantly. However, form-related differences in temperature exchanges seem subtle and may question adaptive disadvantages. We rather support a hypothesis of craniofacial changes in the Australopithecus-Homo transition that are related to brain evolution and craniofacial integration with facial and nasal modifications that contribute to maintain respiratory adaptations related to air conditioning.

Linking Chronic Otitis Media and Nasal Obstruction: A CFD Approach

OBJECTIVES

To investigate a possible relationship between altered nasal flow and chronic otitis media (COM) using computational fluid dynamics (CFD).

STUDY DESIGN

Retrospective case series.

METHODS

Retrospective cohort sample of CT scans from patients with COM and controls without COM to compare the results of various nasal airflow parameters determined by CFD between a group of patients with COM (N = 60) and a control group of subjects without any evidence of ear disease (N = 81). The CT were subjected to various procedures to carry out CFD studies, determining the resistance to nasal flow, the proportion of flow through the right and left nasal cavity, and two nondimensional estimators. The results of CFD studies between patients with COM and controls were compared.

RESULTS

Whereas only 12.3% of the controls had CFD alteration (10 out of 81), 43.3% of the patients suffering COM displayed alterations of our nondimensional parameters R - ϕ (26 out of 60).

CONCLUSIONS

According to our results, the incidence of alterations in nasal airflow by studying with CFD is significantly higher in patients with COM than in controls. To our knowledge, this is the first article linking nasal cavity and COM using a CFD approach. Our results support the hypothesis that nasal flow alterations could be implicated in the etiopathogenesis of the COM.

LEVEL OF EVIDENCE

Level 4 Laryngoscope, 2021.

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.

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

RATIONALE

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Computational modelling of nasal respiratory flow

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

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.

Evaluation of airflow characteristics before and after septoplasty in unilateral cleft patients with a deviated nasal septum: a computational fluid dynamics study

The aim of this study was to evaluate changes in airflow characteristics before and after septoplasty in unilateral cleft lip and palate (UCLP) patients using computational fluid dynamics (CFD) models. The study was designed as a prospective cohort study involving pre- and postoperative computed tomography data from 12 UCLP patients with a deviated nasal septum who underwent septoplasty. CFD analysis of nasal airflow was conducted to study changes in velocity, pressure, volume, nasal resistance, and wall shear stress of the nasal domain before and after surgery. The study results demonstrated a statistically significant difference in pressure drop after septoplasty: median 116.10Pa (interquartile range (IQR) 749.02Pa) preoperative compared with 43.39Pa (IQR 349.01Pa) postoperative (P= 0.004). Maximum wall shear stress was found to be approximately three times lower after septoplasty: median 6.15 Pa (IQR 1908.62 Pa) preoperative versus median 2.51 Pa (IQR 540.06 Pa) postoperative (P=0.002). Changes in nasal resistance were also found to be statistically significant: median 460.59 Pa·s/l (IQR 1946.99 Pa·s/l) preoperative versus median 166.61 Pa·s/l (IQR 694.08 Pa·s/l) postoperative (P=0.04). These values demonstrate significant changes in flow dynamics after surgery indicative of a more uniform airflow pattern and stabilization of the nasal mucosa.

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.

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.

From indoor exposure to inhaled particle deposition: A multiphase journey of inhaled particles

Indoor air quality and its effect on respiratory health are reliant on understanding the level of inhalation exposure, particle inhalability, and particle deposition in the respiratory airway. In the indoor environment, controlling airflow through different ventilation systems can reduce inhalation exposure. This produces a wide variety of complex flow phenomena, such as recirculation, coanda flow, separation, and reattachment. Airborne particles drifting through the air, that move within the breathing region become inhaled into nasal cavity the nostrils. Studies have developed the aspiration efficiency to assist in predicting the fraction of inhaled particles. Inside the nasal cavity, micron and submicron particle deposition occurs in very different ways (inertial impaction, sedimentation, diffusion) and different locations. In addition, fibrous particles such as asbestos are influenced by tumbling effects and its deposition mechanism can include interception. Indoor fluid-particle dynamics related to inhalation exposure and eventual deposition in the respiratory airway is presented. This study involves multi-disciplinary fields involving building science, fluid dynamics, computer science, and medical imaging disciplines. In the future, an integrated approach can lead to digital/in-silico representations of the human respiratory airway able to predict the inhaled particle exposure and its toxicology effect.

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

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

Review: The role of computational simulation in understanding the postoperative sinonasal environment

Nasal surgery improves symptoms in a majority of patients for whom medical treatment has failed. In rhinosinusitis patients, endoscopic sinus surgery aims to alleviate obstruction and re-establish mucociliary clearance. Surgery alters the structure-function relationship within the nasal passage, which is difficult to assess clinically. Computational modelling has been used to investigate this relationship by simulating air flow and environmental variables inside realistic three-dimensional models of the human nasal airway but many questions remain unanswered and need further investigation. The application of computational models to improve pre-surgical planning and post-surgical treatment may not be currently possible due to the absence of knowledge correlating the model-predicted parameters to physiological variables. Links between these parameters to patient outcomes are yet to be established. This article reviews the recent application of computational modelling to understand the nasal structure-function relationship following surgery in patients with sinusitis and nasal obstruction.

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.

Review: The role of computational simulation in understanding the postoperative sinonasal environment

Nasal surgery improves symptoms in a majority of patients for whom medical treatment has failed. In rhinosinusitis patients, endoscopic sinus surgery aims to alleviate obstruction and re-establish mucociliary clearance. Surgery alters the structure-function relationship within the nasal passage, which is difficult to assess clinically. Computational modelling has been used to investigate this relationship by simulating air flow and environmental variables inside realistic three-dimensional models of the human nasal airway but many questions remain unanswered and need further investigation. The application of computational models to improve pre-surgical planning and post-surgical treatment may not be currently possible due to the absence of knowledge correlating the model-predicted parameters to physiological variables. Links between these parameters to patient outcomes are yet to be established. This article reviews the recent application of computational modelling to understand the nasal structure-function relationship following surgery in patients with sinusitis and nasal obstruction.

Sensitivity of nasal airflow variables computed via computational fluid dynamics to the computed tomography segmentation threshold

Computational fluid dynamics (CFD) allows quantitative assessment of transport phenomena in the human nasal cavity, including heat exchange, moisture transport, odorant uptake in the olfactory cleft, and regional delivery of pharmaceutical aerosols. The first step when applying CFD to investigate nasal airflow is to create a 3-dimensional reconstruction of the nasal anatomy from computed tomography (CT) scans or magnetic resonance images (MRI). However, a method to identify the exact location of the air-tissue boundary from CT scans or MRI is currently lacking. This introduces some uncertainty in the nasal cavity geometry. The radiodensity threshold for segmentation of the nasal airways has received little attention in the CFD literature. The goal of this study is to quantify how uncertainty in the segmentation threshold impacts CFD simulations of transport phenomena in the human nasal cavity. Three patients with nasal airway obstruction were included in the analysis. Pre-surgery CT scans were obtained after mucosal decongestion with oxymetazoline. For each patient, the nasal anatomy was reconstructed using three different thresholds in Hounsfield units (-800HU, -550HU, and -300HU). Our results demonstrate that some CFD variables (pressure drop, flowrate, airflow resistance) and anatomic variables (airspace cross-sectional area and volume) are strongly dependent on the segmentation threshold, while other CFD variables (intranasal flow distribution, surface area) are less sensitive to the segmentation threshold. These findings suggest that identification of an optimal threshold for segmentation of the nasal airway from CT scans will be important for good agreement between in vivo measurements and patient-specific CFD simulations of transport phenomena in the nasal cavity, particularly for processes sensitive to the transnasal pressure drop. We recommend that future CFD studies should always report the segmentation threshold used to reconstruct the nasal anatomy.

Investigation of flow characteristics in the maxillary sinus where polypoid changes develop

Flow behavior in the maxillary sinus where polypoid changes develop was investigated using computational fluid dynamics. A nasal cavity model was constructed, after performing a virtual polypectomy based on computed tomography images of a patient, using a computer-aided design software to artificially remove polypoid changes inside the maxillary sinus. Local flow characteristics in the maxillary sinus were examined for one full respiration period. The results showed that the epithelial surfaces where polypoid changes occur are located in the lower part of the maxillary sinus which contains a protruding zone of the sinus and are characterized by stagnation of air during the entire respiration period. Due to the geometric characteristics, a very slow recirculating motion was found to occur in the bulging area for approximately half of the respiration period as a result of interaction with a larger-scale, counter-rotating vortex filling the middle of the maxillary sinus. With a much smaller velocity inside the maxillary sinus compared to that typically found in the airway passage through the middle meatus, both wall shear and pressure changes were found to be vanishingly small along the epithelial surface of the maxillary sinus where polypoid changes were found.

Nasal surgery handled by CFD tools

Annually, hundreds of thousands of surgical interventions to correct nasal airway obstruction are performed throughout the world. Recent studies have noted that a significant number of patients have persistent symptoms of nasal obstruction postoperatively. In the present work, we introduce a new methodology that raises the success rate of nasal cavity surgery. In this procedure, the surgeon performs virtual surgery on a 3D nasal model of a patient prior to the real surgery. The main goal of the methodology is to guide the surgeon throughout the virtual operation using mathematical estimators based on CFD results. The virtual surgery intervention ends as soon as the estimators fall into a region of a Cartesian coordinate system with a high success probability. This region is defined according to a statistical analysis of estimators corresponding to sets of healthy and diseased cavities. As examples of this application, this study includes 2 surgical operations performed with this innovative methodology on patients with severe nasal obstruction. The patients underwent nasal surgery according to the final nasal geometry revealed by CFD-guided virtual surgery. Currently, both subjects show high degrees of satisfaction.

A CFD approach to understand nasoseptal perforations

INTRODUCTION

Nasoseptal perforations (NSP) are becoming common in the modern world, and can cause a wide variety of symptoms, including a sensation of nasal obstruction, epistaxis, crusting, dryness, headache, nasal pain and a whistling sound. There is an extensive range of surgical treatment techniques, but reported results were rarely statistically significant. The lack of consistent surgical results may be related to the lack of knowledge about the pathophysiology of NSP and how they affect the nasal flow. Computational fluid dynamics (CFD) has proved to be a very useful tool to study nasal function.

METHODS

We have used CFD software (the program MECOMLAND® and the Digbody® tool for virtual surgery) to investigate the behaviour of the parameters R-[Formula: see text] based on CFD results, when four subjects underwent virtual surgery to induce a septal perforation: two subjects with healthy noses and two patients suffering from nasal airway obstruction. For each case a CFD study was performed, before and after creating an anterior (close to nostrils) or a posterior (close to choanae) NSP.

RESULTS

In all cases analyzed, a posterior septal perforation did not result in a significant volumetric flow rate [Formula: see text] through the perforation between nasal passages. However, for anterior defects only in those nasal cavities considered diseased or unhealthy, high values of [Formula: see text] were found.

CONCLUSION

The induced NSP only rendered significant flow alterations in noses with preexisting nasal airway obstruction alterations, whereas in nasal cavities considered as normal the creation of a NSP did not produce significant differences between both sides. We strongly suggest that this finding can explain the variety of symptoms and the number of asymptomatic patients bearing NSP.

Flow and air conditioning simulations of computer turbinectomized nose models

Air conditioning for the human respiratory system is the most important function of the nose. When obstruction occurs in the nasal airway, turbinectomy is used to correct such pathology. However, mucosal atrophy may occur sometime after this surgery when it is overdone. There is not enough information about long-term recovery of nasal air conditioning performance after partial or total surgery. The purpose of this research was to assess if, based on the flow and temperature/humidity characteristics of the air intake to the choana, partial resection of turbinates is better than total resection. A normal nasal cavity geometry was digitized from tomographic scans and a model was printed in 3D. Dynamic (sinusoidal) laboratory tests and computer simulations of airflow were conducted with full agreement between numerical and experimental results. Computational adaptations were subsequently performed to represent six turbinectomy variations and a swollen nasal cavity case. Streamlines along the nasal cavity and temperature and humidity distributions at the choana indicated that the middle turbinate partial resection is the best alternative. These findings may facilitate the diagnosis of nasal obstruction and can be useful both to plan a turbinectomy and to reduce postoperative discomfort. Graphical Abstract ᅟ.

Flow visualization through particle image velocimetry in realistic model of rhesus monkey's upper airway

Studies concerning inhalation toxicology and respiratory drug-delivery systems require biological testing involving experiments performed on animals. Particle image velocimetry (PIV) is an effective in vitro technique that reveals detailed inhalation flow patterns, thereby assisting analyses of inhalation exposure to various substances. A realistic model of a rhesus-monkey upper airway was developed to investigate flow patterns in its oral and nasal cavities through PIV experiments performed under steady-state constant inhalation conditions at various flow rates-4, 10, and 20 L/min. Flow rate of the fluid passing through the inlet into the trachea was measured to obtain characteristic flow mechanisms, and flow phenomena in the model were confirmed via characterized flow fields. It was observed that increase in flow rate leads to constant velocity profiles in upper and lower trachea regions. It is expected that the results of this study would contribute to future validation of studies aimed at developing in silico models, especially those involving computational fluid dynamic (CFD) analysis.

Robust nondimensional estimators to assess the nasal airflow in health and disease

There are significant variations of both human nose shapes and airflow patterns inside nasal cavities, so it is difficult to provide a comprehensive medical identification using a universal template for what otolaryngologists consider normal breathing at rest. In addition, airflow patterns present even more random characteristics in diseased nasal cavities. To give a medical assessment to differentiate the nasal cavities in health and disease, we propose 2 nondimensional estimators obtained from both medical images and computational fluid dynamics. The first mathematical estimator ϕ is a function of geometric features and potential asymmetries between nasal passages, while the second estimator R represents in fluid mechanics terms the total nasal resistance that corresponds to the atmosphere-choana pressure drop. These estimators only require global information such as nasal geometry and magnitudes of flow determined by simulations under laminar conditions. We find that these estimators take low and high values for healthy and diseased nasal cavities, respectively. Our study, based on 24 healthy and 25 diseased Caucasian subjects, reveals that there is an interval of values associated with healthy cavities that clusters in a small region of the plane ϕ-R. Therefore, these estimators can be seen as a first approximation to provide nasal airflow data to the clinician in a noninvasive method, as the computed tomography scan that provides the required images is routinely obtained as a result of the preexisting naso-sinusal condition.

New CFD tools to evaluate nasal airflow

Computational fluid dynamics (CFD) is a mathematical tool to analyse airflow. As currently CFD is not a usual tool for rhinologists, a group of engineers in collaboration with experts in Rhinology have developed a very intuitive CFD software. The program MECOMLAND^® only required snapshots from the patient's cross-sectional (tomographic) images, being the output those results originated by CFD, such as airflow distributions, velocity profiles, pressure, temperature, or wall shear stress. This is useful complementary information to cover diagnosis, prognosis, or follow-up of nasal pathologies based on quantitative magnitudes linked to airflow. In addition, the user-friendly environment NOSELAND^® helps the medical assessment significantly in the post-processing phase with dynamic reports using a 3D endoscopic view. Specialists in Rhinology have been asked for a more intuitive, simple, powerful CFD software to offer more quality and precision in their work to evaluate the nasal airflow. We present MECOMLAND^® and NOSELAND^® which have all the expected characteristics to fulfil this demand and offer a proper assessment with the maximum of quality plus safety for the patient. These programs represent a non-invasive, low-cost (as the CT scan is already performed in every patient) alternative for the functional study of the difficult rhinologic case. To validate the software, we studied two groups of patients from the Ear Nose Throat clinic, a first group with normal noses and a second group presenting septal deviations. Wall shear stresses are lower in the cases of normal noses in comparison with those for septal deviation. Besides, velocity field distributions, pressure drop between nasopharynx and the ambient, and flow rates in each nostril were different among the nasal cavities in the two groups. These software modules open up a promising future to simulate the nasal airflow behaviour in virtual surgery intervention scenarios under different pressure or temperature conditions to understand the effects on nasal airflow.

Numerical study on the air conditioning characteristics of the human nasal cavity

The air-conditioning characteristics of the human nasal cavity were investigated using computational fluid dynamics. The wall layer was modeled as a heat conducting layer consisting of water with constant thickness placed on top of epithelial cells. By assuming constant tissue temperature, prescribed to be 36 ^°C, which is close to the alveolar condition, the proposed wall model yielded a spatially varying surface temperature distribution that is in reasonable agreement with the measurement studies in the literature. The results show that the regions of the main airway between the nasal valve, and the anterior of the middle turbinate were shown to have relatively low temperatures, whereas the superior meatus exhibited relatively high temperature. Water vapor flux evaluated at the surface of the mucus layer was found to be quite large in the region between the posterior of the vestibule and the anterior of the middle turbinate. Comparing the results obtained from the present model to those obtained with a constant surface temperature boundary condition of 32.6 °C or 34 °C revealed that temperature, and absolute humidity of the airflow increased faster through the turbinated airway passage. Even in the presence of sizable differences in the distributions of surface temperature and water vapor concentration, distributions of relative humidity of the air were found to be quite similar regardless of temperature boundary conditions.

The incompressibility assumption in computational simulations of nasal airflow

Most of the computational works on nasal airflow up to date have assumed incompressibility, given the low Mach number of these flows. However, for high temperature gradients, the incompressibility assumption could lead to a loss of accuracy, due to the temperature dependence of air density and viscosity. In this article we aim to shed some light on the influence of this assumption in a model of calm breathing in an Asian nasal cavity, by solving the fluid flow equations in compressible and incompressible formulation for different ambient air temperatures using the OpenFOAM package. At low flow rates and warm climatological conditions, similar results were obtained from both approaches, showing that density variations need not be taken into account to obtain a good prediction of all flow features, at least for usual breathing conditions. This agrees with most of the simulations previously reported, at least as far as the incompressibility assumption is concerned. However, parameters like nasal resistance and wall shear stress distribution differ for air temperatures below [Formula: see text]C approximately. Therefore, density variations should be considered for simulations at such low temperatures.

The impact of endoscopic sinus surgery on paranasal physiology in simulated sinus cavities

BACKGROUND

Surgery improves symptoms for the majority of chronic rhinosinusitis (CRS) patients; however, physiological changes in the sinus cavities remain poorly characterized. Direct measurement of changes in airflow, pressure, temperature, humidity, and intranasal spray distribution following surgery is technically challenging. Accordingly, we have used computational fluid dynamic modeling to quantify how these parameters change postoperatively.

METHODS

Computed tomography images from a normal control, a patient with CRS preoperatively and postoperatively, and a patient following an endoscopic Lothrop procedure (ELP) were used to create 4 three-dimensional models of the sinus cavities. Changes in physiologic parameters and topical drug distribution were modeled (inhaled air at 16°C and 10% humidity) at the maxillary ostium, frontal recess, and sphenoid ostium.

RESULTS

Large differences were seen between models. Following surgery, the maxillary ostia were found on average to be cooler (by 2.4°C), with an increased airflow (0.26 m/second; from 0 m/second), and a 9% reduction in absolute humidity. Sphenoid ostial parameters followed a similar trend. Significant changes in frontal recess physiology were seen following ELP in which the recess was 4.2°C cooler, had increased airflow (0.76 m/second) and a 17% reduction in absolute humidity. Topical drug distribution increased with surgery, particularly after ELP.

CONCLUSION

Surgery changes the geometry and physiology of the paranasal sinuses. These changes are likely to have an impact on wound healing, mucociliary function, and microbial ecology in postoperative cavities. Application of this model to further understand the effects of surgery may help to optimize surgical techniques and improve topical drug delivery.

Assessment of nasal resistance using computational fluid dynamics

Anterior rhinomanometry is the current gold standard for the objective assessment of nasal breathing by determining the nasal resistance. However, computational fluid dynamics would allow spatially and temporally well- resolved investigation of additional flow parameters. In this study, measured values of nasal resistance are compared with measured values. An unclear discrepancy between the two methods was found, suggesting further investigation.

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.

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.

A computational analysis of nasal vestibule morphologic variabilities on nasal function

Although advances in computational modeling have led to increased understanding of nasal airflow, not much is known about the effects of normal sinonasal anatomic variabilities on nasal function. In this study, three distinct variations in the human nasal vestibule airspace that have not been previously described were identified. Computational fluid dynamics modeling of nasal airflow profile in each identified variation of nasal vestibule phenotype was conducted to assess the role of these phenotypes on nasal physiology. Three-dimensional reconstructions of the nasal geometry in sixteen subjects with normal radiographic sinonasal images were created and each respective unilateral nasal cavity was classified as Notched, Standard, or Elongated phenotype based nasal vestibule morphology. Steady state, laminar and incompressible flow simulations were performed in the nasal geometries under physiological, pressure-driven conditions with constant inspiratory pressure. Results showed that at localized regions of the unilateral nasal cavity, average resistance was significantly different among nasal vestibule phenotypes. However, global comparison from nostril to choana showed that average resistance was not significantly different across phenotypes; suggesting that with normal anatomic variations, the nose has a natural compensatory mechanism that modulates localized airflow in order to achieve a desired amount of global airflow.

Large-scale CFD simulations of the transitional and turbulent regime for the large human airways during rapid inhalation

The dynamics of unsteady flow in the human large airways during a rapid inhalation were investigated using highly detailed large-scale computational fluid dynamics on a subject-specific geometry. The simulations were performed to resolve all the spatial and temporal scales of the flow, thanks to the use of massive computational resources. A highly parallel finite element code was used, running on two supercomputers, solving the transient incompressible Navier-Stokes equations on unstructured meshes. Given that the finest mesh contained 350 million elements, the study sets a precedent for large-scale simulations of the respiratory system, proposing an analysis strategy for mean flow, fluctuations and wall shear stresses on a rapid and short inhalation (a so-called sniff). The geometry used encompasses the exterior face and the airways from the nasal cavity, through the trachea and up to the third lung bifurcation; it was derived from a contrast-enhanced computed tomography (CT) scan of a 48-year-old male. The transient inflow produces complex flows over a wide range of Reynolds numbers (Re). Thanks to the high fidelity simulations, many features involving the flow transition were observed, with the level of turbulence clearly higher in the throat than in the nose. Spectral analysis revealed turbulent characteristics persisting downstream of the glottis, and were captured even with a medium mesh resolution. However a fine mesh resolution was found necessary in the nasal cavity to observe transitional features. This work indicates the potential of large-scale simulations to further understanding of airway physiological mechanics, which is essential to guide clinical diagnosis; better understanding of the flow also has implications for the design of interventions such as aerosol drug delivery.

Dynamics of airflow in a short inhalation

During a rapid inhalation, such as a sniff, the flow in the airways accelerates and decays quickly. The consequences for flow development and convective transport of an inhaled gas were investigated in a subject geometry extending from the nose to the bronchi. The progress of flow transition and the advance of an inhaled non-absorbed gas were determined using highly resolved simulations of a sniff 0.5 s long, 1 l s⁻¹ peak flow, 364 ml inhaled volume. In the nose, the distribution of airflow evolved through three phases: (i) an initial transient of about 50 ms, roughly the filling time for a nasal volume, (ii) quasi-equilibrium over the majority of the inhalation, and (iii) a terminating phase. Flow transition commenced in the supraglottic region within 20 ms, resulting in large-amplitude fluctuations persisting throughout the inhalation; in the nose, fluctuations that arose nearer peak flow were of much reduced intensity and diminished in the flow decay phase. Measures of gas concentration showed non-uniform build-up and wash-out of the inhaled gas in the nose. At the carina, the form of the temporal concentration profile reflected both shear dispersion and airway filling defects owing to recirculation regions.

Variability in Nose-to-Lung Aerosol Delivery

Nasal delivery of lung targeted pharmaceutical aerosols is ideal for drugs that need to be administered during high flow nasal cannula (HFNC) gas delivery, but based on previous studies losses and variability through both the delivery system and nasal cavity are expected to be high. The objective of this study was to assess the variability in aerosol delivery through the nose to the lungs with a nasal cannula interface for conventional and excipient enhanced growth (EEG) delivery techniques. A database of nasal cavity computed tomography (CT) scans was collected and analyzed, from which four models were selected to represent a wide range of adult anatomies, quantified based on the nasal surface area-to-volume ratio (SA/V). Computational fluid dynamics (CFD) methods were validated with existing in vitro data and used to predict aerosol delivery through a streamlined nasal cannula and the four nasal models at a steady state flow rate of 30 L/min. Aerosols considered were solid particles for EEG delivery (initial 0.9 μm and 1.5 μm aerodynamic diameters) and conventional droplets (5 μm) for a control case. Use of the EEG approach was found to reduce depositional losses in the nasal cavity by an order of magnitude and substantially reduce variability. Specifically, for aerosol deposition efficiency in the four geometries, the 95% confidence intervals (CI) for 0.9 and 5 μm aerosols were 2.3-3.1 and 15.5-66.3%, respectively. Simulations showed that the use of EEG as opposed to conventional methods improved delivered dose of aerosols through the nasopharynx, expressed as penetration fraction (PF), by approximately a factor of four. Variability of PF, expressed by the coefficient of variation (CV), was reduced by a factor of four with EEG delivery compared with the control case. Penetration fraction correlated well with SA/V for larger aerosols, but smaller aerosols showed some dependence on nasopharyngeal exit hydraulic diameter. In conclusion, results indicated that the EEG technique not only improved lung aerosol delivery, but largely eliminated variability in both nasal depositional loss and lung PF in a newly developed set of nasal airway models.

Numerical investigation of regional particle deposition in the upper airway of a standing male mannequin in calm air surroundings

A 3-D realistic computational model of the airway system integrated into a standing male mannequin was developed. The computational domain includes the regions around the mannequin and the inside of the airway passages. The simulation was performed for low activity breathing rates with calm air around the mannequin. The flowfield of the inhaled air was first obtained from solving the Navier-Stokes and continuity equations. Then the particles were released in the domain around the mannequin and their trajectories were evaluated by using the Lagrangian approach for solving the particle equation of motion. The regional aerosols deposition was evaluated for different parts of the human airway system and the results were compared with those obtained from the separate modeling of the airway system without the interaction of the airflow with the mannequin external face. The results showed when the upper airway is integrated into the mannequin, the regional deposition of inhaled particles mainly changes in the airway system.

Accuracy evaluation of a numerical simulation model of nasal airflow

CONCLUSION

Our numerical simulation model provides an accurate reflection of nasal airflow, and the results were validated by clinical measurements.

OBJECTIVES

To evaluate the accuracy of a numerical simulation model of nasal airflow.

METHODS

Ten volunteers with normal nasal cavities underwent CT, acoustic rhinometry, and rhinomanometry. CT data were uploaded into Mimics, ICEM-CFD, Fluent, and CFD-Post software for three-dimensional modeling, finite element grid division, transient calculations, and analysis, respectively. Velocity and pressure data of airflow were obtained during the normal respiratory cycle. The accuracy of the simulation was evaluated by two methods: acoustic rhinometry measurements were used to evaluate the accuracy of the anatomic model, and rhinomanometry measurements were used to evaluate the accuracy of the nasal resistance values obtained by numerical simulation.

RESULTS

There were no significant differences between the values describing the model and the acoustic rhinometry measurements, the nasal resistance values obtained by numerical simulation. The airflow through the nasal cavity was mainly laminar. The maximum velocities were measured at the nasal valve, the amplitudes of all velocity curves at locations beyond the nasal valve were reduced. The amplitudes of the pressure curves increased from the front to the back of the airway.

Effects of the ambient temperature on the airflow across a Caucasian nasal cavity

We analyse the effects of the air ambient temperature on the airflow across a Caucasian nasal cavity under different ambient temperatures using CFD simulations. A three-dimensional nasal model was constructed from high-resolution computed tomography images for a nasal cavity from a Caucasian male adult. An exhaustive parametric study was performed to analyse the laminar-compressible flow driven by two different pressure drops between the nostrils and the nasopharynx, which induced calm breathing flow rates ࣈ 5.7 L/min and ࣈ 11.3 L/min. The inlet air temperature covered the range - 10(o) C ⩽ To ⩽50(o) C. We observed that, keeping constant the wall temperature of the nasal cavity at 37(o) C, the ambient temperature affects mainly the airflow velocity into the valve region. Surprisingly, we found an excellent linear relationship between the ambient temperature and the air average temperature reached at different cross sections, independently of the pressure drop applied. Finally, we have also observed that the spatial evolution of the mean temperature data along the nasal cavity can be collapsed for all ambient temperatures analysed with the introduction of suitable dimensionless variables, and this evolution can be modelled with the help of hyperbolic functions, which are based on the heat exchanger theory.