Normative ranges of nasal airflow variables in healthy adults

Algorithm for Nasal Breathing Impairment Evaluation

Assessing patients with complaints of nasal obstruction has traditionally been done by evaluation of the nasal airway looking for fixed or dynamic obstructive locations that could impair nasal airflow. Not infrequently, however, symptoms of nasal obstruction do not match the clinical examination of the nasal airway. Addressing this subset of patients may be a challenge to the surgeon. Evaluation of patients with symptoms of nasal obstruction should include a combination of a patient-reported assessment of nasal breathing and at least one objective method for measuring nasal airflow or nasal airway resistance or dimensions. This will allow distinction between patients with symptoms of nasal obstruction and low airflow or high nasal airway resistance and patients with similar symptoms but whose objective evaluation demonstrates normal nasal airflow or normal airway dimensions or resistance. Patients with low nasal airflow or high nasal airway resistance will require treatment to increase nasal airflow as a necessary step to improve symptoms, whereas patients with normal nasal airflow or nasal airway resistance will require a multidimensional assessment looking for less obvious causes of impaired nasal breathing sensation.

Estimation of Nasal Airway Cross-sectional Area From Endoscopy Using Depth Maps: A Proof-of-Concept Study

OBJECTIVE

Endoscopy is routinely used to diagnose obstructive airway diseases. Currently, endoscopy is only a visualization technique and does not allow quantification of airspace cross-sectional areas (CSAs). This pilot study tested the hypothesis that CSAs can be accurately estimated from depth maps created from virtual endoscopy videos.

STUDY DESIGN

Cross-sectional.

SETTING

Academic tertiary medical center.

METHODS

Virtual endoscopy and depth map videos of the nasal cavity were digitally created based on anatomically accurate three-dimensional (3D) models built from computed tomography scans of 30 subjects. A software tool was developed to outline the airway perimeter and estimate the airspace CSA from the depth maps. Two otolaryngologists used the software tool to estimate the nasopharynx CSA and the nasal valve minimal CSA (mCSA) in the left and right nasal cavities. Model validation statistics were performed.

RESULTS

Nasopharynx CSA had a median percent error of 3.7% to 4.6% when compared to the true values measured in the 3D models. Nasal valve mCSA had a median percent error of 22.7% to 33.6% relative to the true values. Raters successfully used the software tool to identify subjects with nasal valve stenosis (ie, mCSA < 0.20 cm2) with a sensitivity of 83.3%, specificity ≥ 90.7%, and classification accuracy ≥ 90.0%. Interrater and intrarater agreements were high.

CONCLUSION

This study demonstrates that airway CSAs in 3D models can be accurately estimated from depth maps. The development of artificial intelligence algorithms to compute depth maps may soon allow the quantification of airspace CSAs from clinical endoscopies.

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.

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

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

Effects of 3D Airway Geometry on the Airflow of Adults with Cleft Lip and Palate and Obstructive Sleep Apnea: A Functional Imaging Study

Objective  Individuals with cleft lip and palate (CLP) are at a high risk of developing obstructive sleep apnea (OSA). Hypothetically, the severity of OSA might be associated with the morphology of the upper airway (UAW) and the characteristics of the airflow. Thus, the present study aimed to assess and compare, in adults with CLP and skeletal class-III discrepancy, with or without OSA, simulations of airflow resistance and pressure according to the geometrical characteristics of the UAW and cephalometric parameters. Materials and Methods  According to the results of type-I polysomnography tests, the sample ( n  = 21) was allocated in 2 groups: 1) without OSA (N-OSA; n  = 6); and 2) with OSA (OSA; n  = 15). Cephalometric measurements were performed on the cone-beam computed tomography (CBCT) scans of the groups. After three-dimensional (3D) reconstructions, the volume (V) and minimal cross-sectional area (mCSA) of the UAW were generated. Computational fluid dynamics (CFD) simulations were used to assess key airflow characteristics. The results were presented at a significance level of 5%. Results  The UAW pressure values and airway resistance did not differ between the groups, but there was a tendency for more negative pressures (26%) and greater resistance (19%) in the OSA group. Volume and mCSA showed a moderate negative correlation with resistance and pressure. The more inferior the hyoid bone, the more negative the pressures generated on the pharyngeal walls. Conclusion  The position of the hyoid bone and the geometry of the UAW (V and mCSA) exerted effects on the airway-airflow resistance and pressure. However, key airflow characteristics did not differ among subjects with CLP, were they affected or not by OSA.

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.

CT analysis of the anterior nasal airway based on the direction of nasal airflow in patients with nasal obstruction and trauma controls

PURPOSE

The anterior nose is the nasal segment with the highest resistance to airflow. In a hospital-based case-control study, we compared cross-sectional areas of the nasal cavities anterior to the piriform aperture determined by computed tomography (CT-CSA) in patients with nasal obstruction (cases) and unselected patients with trauma unrelated to the head and face (controls).

METHODS

CT-CSA could be reproducibly identified at angles of 0^o, 30°, 60°, and 90° to the nasal floor approximately perpendicular to the arcuate direction of nasal airflow using bony landmarks. CT-CSA were manually segmented and compared in cases and controls. In cases, we compared CT-CSA at 30° (CT-CSA_30-narrow) with the minimum cross-sectional area determined by acoustic rhinometry (AR-MCA_1-narrow), each on the narrower side.

RESULTS

CT-CSA ranged from 7 to 250 mm² with an average of 100 mm² per nasal side. Side differences of the nasal airways indicating asymmetry of the nasal airways were greater in 40 cases than in 44 controls (p < 0.003). Moreover, bilateral CT-CSA were significantly smaller in cases than in controls (p < 0.001). CT-CSA_30-narrow did not significantly correlate with AR-MCA_1-narrow (r = 0.33; p = 0.07) and on average was 58% smaller than AR-MCA_1-narrow.

CONCLUSIONS

Cross-sectional areas of the anterior nose perpendicular to the direction of nasal airflow, which is considered relevant in terms of flow physics, can be reliably measured using CT. Anterior nasal cavities in patients with nasal obstruction were more asymmetric and, as a whole, narrower than in controls, the latter of which is not corrected by routine septoplasty.

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

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

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.

Acoustic Aerosol Delivery: Assessing of Various Nasal Delivery Techniques and Medical Devices on Intrasinus Drug Deposition

This study aims to evaluate the impact of the nasal delivery technique and nebulizing technologies (using different frequencies of oscillating airflow) for acoustic aerosol targeting of maxillary sinuses. Sodium fluoride (chemical used as a marker), tobramycin (drug used as a marker) and ^99mTc-DTPA (radiolabel aerosol) were used to assess the intrasinus aerosol deposition on a nasal cast. Two commercial medical devices (PARI SINUS nebulizer and NL11SN ATOMISOR nebulizer) and various nasal delivery techniques (one or two nostrils connected to the aerosol inlet, the patient with the soft palate closed or open during the acoustic administration of the drug, the presence or not of flow resistance in the nostril opposite to the one allowing the aerosol to be administered) were evaluated. The closed soft palate condition showed a significant increase in drug deposition even though no significant difference in the rest of the nasal fossae was noticed. Our results clearly demonstrated a higher intrasinus aerosol deposition (by a factor 2-3; respectively 0.03 ± 0.007% vs. 0.003 ± 0.0002% in the right maxillary sinus and 0.027 ± 0.006% vs. 0.013 ± 0.004% in the left maxillary sinus) using the acoustic airflow generated by the PARI SINUS compared to the NL11SN ATOMISOR. The results clearly demonstrated that the optimal conditions for aerosol deposition in the maxillary sinuses were obtained with a closed soft palate. Thus, the choice of the nebulizing technology (and mainly the frequency of the pulsating aerosol generated) and also the recommendation of the best nasal delivery technique are key factors to improve intrasinus aerosol deposition.

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

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

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

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

Spatial aerosol deposition correlated to anatomic feature development in 6-year-old upper airway computational models

The upper airways of children undergo developmental changes around age 6, yielding differences between adult and pediatric anatomies. These differences include the cricoid ring area shape, the location of narrowest constriction, and the angle of the epiglottis, all of which are expected to alter local fluid dynamic profiles and subsequent upper airway deposition and downstream aerosol delivery of inhaled therapeutics. In this work, we quantify "pediatric"-like and "adult"-like geometric and fluid dynamic features of two computed tomography (CT)-scan derived models of 6-year-old upper airways in healthy subjects and compare to an idealized model. The two CT-scan models had a mixture of "adult"- and "pediatric"-like anatomic features, with Subject B exhibiting more "pediatric"-like features than Subject A, while the idealized model exhibited entirely "adult"-like features. By computational fluid-particle dynamics, these differences in anatomical features yielded distinct local fluid profiles with altered aerosol deposition between models. Notably, the idealized model better predicted deposition characteristics of Subject A, the more "adult"-like model, including the relationship between the impaction parameter, dp2Q and the fraction of deposition across a range of flow rates and particle diameters, as well as deposition of an approximate pharmaceutical particle size distribution model. Our results with even this limited dataset suggest that there are key personalized metrics that are influenced by anatomical development, which should be considered when developing pediatric inhalable therapeutics. Quantifying anatomical development and correlating to aerosol deposition has the potential for high-throughput developmental characterization and informing desired aerosol characteristics for pediatric applications.

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.

CFD SIMULATIONS AND ANALYSES OF ASYMPTOMATIC AND SYMPTOMATIC NASAL AIRWAY OBSTRUCTIONS

The numerical simulations of the flow in nasal airways were performed for two different clinical cases. The results comprised the distributions of scalars at five different sections and included contours of pressure, velocity magnitude, turbulent kinetic energy and vorticity magnitude. Simulations showed the air branching occurring at the inferior meatus is unaffected by the variations in the volumetric flow rate or the changes in the flow regime through the olfactory cleft. However, the contractions and the following rapid change in the cross-section of the nasopharynx preclude the upward penetration of the vacuum field set by the lungs during the inhalation process. As a result, considerably low velocities and significant cross-sectional nonuniformities are observed, which lead to the appearances of the secondary flow structures and strong unsteadiness. Increased interactions between the airflow and the walls of the nasal cavity resulted in an increase in the vorticity on the right middle meatus and upper inferior meatus. The vorticity was also very high in the nostrils, where the flow was not fully developed.

A systematic analysis of surgical interventions for the airway in the mature unilateral cleft lip nasal deformity: a single case study

PURPOSE

Individuals with unilateral cleft lip nasal deformity (uCLND) often require rhinoplasty in adolescence to correct nasal obstruction. The intent of this study is to identify sites of greatest nasal obstruction and evaluate the effects of isolated and combinations of simulated surgical procedures on these sites using computational fluid dynamics (CFD).

METHODS

Computed tomography imaging of an adolescent subject with uCLND was converted to an anatomically accurate three-dimensional nasal airway model. Initial analysis was performed to identify anatomic sites of obstruction based on CFD computed resistance values. Virtual surgery procedures corresponding to common uCLND surgical interventions were simulated. Resulting airspace models were then analyzed after conducting airflow and heat transfer simulations.

RESULTS

The preoperative model had 21 obstructed sites with a nasal resistance of 0.075 Pa s/mL. Following simulated surgical procedures with functional interventions alone and in combinations, the three virtual surgery models with most improved nasal airflow were inferior turbinate reduction (ITR) with posterior septoplasty (resistance = 0.054 Pa s/ml, reduction in 14 of 21 obstructed sites), ITR with anterior septoplasty (resistance = 0.058 Pa s/ml, reduction in 8 of 21 obstructed sites), and ITR with both anterior and posterior septoplasty (resistance = 0.052 Pa s/ml, reduction in 17 of 21 obstructed sites).

CONCLUSION

This study introduces a new technique for analysis of the impact of different simulated surgical interventions on uCLND-induced nasal obstruction. In this subject, simulated septoplasty with ITR on the non-cleft side provided maximal relief of nasal obstruction. The proposed technique can be further studied for possible utility in analyzing potential surgical interventions for optimal relief of nasal obstruction in patients with uCLND.

Nasally inhaled therapeutics and vaccination for COVID-19: Developments and challenges

The nose is the initial site of viral infection, replication, and transmission in the human body. Nasally inhaled vaccines may act as a promising alternative for COVID-19 management in addition to intramuscular vaccination. In this review, the latest developments of nasal sprays either as repurposed or antiviral formulations were presented. Nasal vaccines based on traditional medicines, such as grapefruit seed extract, algae-isolated carrageenan, and Yogurt-fermenting Lactobacillus, are promising and under active investigations. Inherent challenges that hinder effective intranasal delivery were discussed in detail, which included nasal device issues and human nose physiological complexities. We examined factors related to nasal spray administration, including the nasal angiotensin I converting enzyme 2 (ACE2) locations as the delivery target, nasal devices, medication translocation after application, delivery methods, safety issues, and other nasal delivery options. The effects of human factors on nasal spray efficacy, such as nasal physiology, disease-induced physiological modifications, intersubject variability, and mucociliary clearance, were also examined. Finally, the potential impact of nasal vaccines on COVID-19 management in the developing world was discussed. It is concluded that effective delivery of nasal sprays to ACE2-rich regions is urgently needed, especially in the context that new variants may become unresponsive to current vaccines and more refractory to existing therapies.

Toward automatic atlas-based surgical planning for septoplasty

PURPOSE

Surgery for nasal airway obstruction (NAO) has a high failure rate, with up to 50% of patients reporting persistent symptoms postoperatively. Virtual surgery planning has the potential to improve surgical outcomes, but current manual methods are too labor-intensive to be adopted on a large scale. This manuscript introduces an automatic atlas-based approach for performing virtual septoplasties.

METHODS

A cohort of 47 healthy subjects and 26 NAO patients was investigated. An atlas of healthy nasal geometry was constructed. The automatic virtual septoplasty method consists of a multi-stage registration approach to fit the atlas to a target NAO patient, automatically segment the patient's septum and airway, and deform the patient image to have a non-deviated septum.

RESULTS

Our automatic virtual septoplasty method straightened the septum successfully in 18 out of 26 NAO patients (69% of cases). In these cases, the ratio of the higher to the lower airspace cross-sectional areas in the left and right nasal cavities improved from 1.47 ± 0.45 to 1.16 ± 0.33 in the region surrounding the septal deviation, showing that the nasal airway became more symmetric after virtual septoplasty.

CONCLUSION

This automated virtual septoplasty technique has the potential to greatly reduce the effort required to perform computational fluid dynamics (CFD) analysis of nasal airflow for NAO surgical planning. Future studies are needed to investigate if virtual surgery planning using this method is predictive of subjective symptoms in NAO patients after septoplasty.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

CLINICAL QUESTION/LEVEL OF EVIDENCE

Risk, V.

Which Nasal Airway Dimensions Correlate with Nasal Airflow and with Nasal Breathing Sensation?

Background: Rhinoplasty modifies the nasal pyramid, thereby also modifying the nasal airway. Objectives: To correlate the sensation of nasal breathing, as measured by patient-reported outcome measures, and nasal airflow, as assessed by peak nasal inspiratory flow (PNIF), with nasal airway dimensions, as measured on computed tomography (CT) images. Methods: Fifty Caucasian patients were studied through visual analogue scale (VAS), nasal obstruction symptom evaluation (NOSE) and PNIF. Measurements of the nasal airway were made on CT images: minimal distance between septum and inferior and middle turbinates, nasal valve angle, and nasal valve area. Results: There was a significant association between PNIF and nasal valve area, between VAS and the narrower nasal valve angle and between NOSE and minimal distance between septum and middle turbinate of the narrower side. Conclusions: This study suggests that the dimensions of the nasal valve and of the middle nasal airway have a substantial impact on nasal breathing capacity. It also highlights the importance of unilateral nasal airway obstruction to nasal breathing.

Nasal airflow of patient with septal deviation and allergy rhinitis

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

Dynamic characteristics of heat capacity of the human nasal cavity during a respiratory cycle

The dynamic characteristics of air-conditioning in the human nasal cavity during a respiratory cycle were investigated using unsteady numerical simulations to assess whether inhaled air is sufficiently conditioned by the nasal cavity. Variations in the epithelial surface temperature, surface heat, and water vapor fluxes were found to vary significantly during inspiration while providing substantial air conditioning to the inhaled air, but variations and magnitudes were significantly reduced during the expiration period. Air temperature (31.3-35.3 °C) and relative humidity (85.1-100 %) in the nasopharynx exhibited significant variations during inspiration. Flow rate-weighted average values of the air temperature and relative humidity during inspiration were estimated to be 32.0 °C and 89.1 %, respectively. Inhaled air did not attain alveolar conditions before reaching the nasopharynx, and was therefore thought to be insufficiently conditioned by the nasal cavity alone. A steady flow of approximately 250 mL/s appears to be useful for evaluating the accumulated thermal state of air in the nasopharynx during inspiration.

The Recumbent Position Affects Nasal Resistance: A Systematic Review and Meta-Analysis

OBJECTIVE

Nasal diseases are among the main motives for the early discontinuation of continuous positive airway pressure therapy and for long-term therapeutic compliance with mandibular advancement device. Although our clinical experience leads us to the belief that recumbency impacts nasal airflow in some patient populations, there is no consensus regarding the magnitude of this effect and the specific group of patients who are the most affected by this condition. In this study, we conducted a meta-analysis to assess the effect of the recumbent position on nasal resistance and nasal airflow.

REVIEW METHODS

PubMed (Medline), Cochrane Library, EMBASE, Scopus, and SciELO databases were checked for relevant studies by two members of the YO-IFOS study group. The two authors extracted the data. The main outcome was expressed as the difference between nasal resistance and nasal airflow before and after recumbency.

RESULTS

Nine studies with a total population of 291 individuals were included in the meta-analysis for nasal resistance after recumbency. We found a statistically significant difference in nasal airway resistance of -0.18 Pa sec/cm³ as compared to before and after recumbency through rhinomanometry (RMM) analysis. A subgroup analysis revealed a variation of -0.20 Pa sec/cm³ for patients with snoring or sleep apnea and - 0.10 Pa sec/cm³ for healthy individuals. Regarding nasal airflow measured with RMM, three studies (n = 32) in asymptomatic controls revealed a statistically significant difference of 47.33 ml/sec.

CONCLUSIONS

Recumbency increases nasal resistance and diminishes nasal airflow. This finding is of utmost importance in snorers and sleep apnea patients. Laryngoscope, 2021.

Existence of a Neutral-Impact Maxillo-Mandibular Displacement on Upper Airways Morphology

Current scientific evidence on how orthognathic surgery affects the airways morphology remains contradictory. The aim of this study is to investigate the existence and extension of a neutral-impact interval of bony segments displacement on the upper airways morphology. Its upper boundary would behave as a skeletal displacement threshold differentiating minor and major jaw repositioning, with impact on the planning of the individual case. Pre- and post-operative cone beam computed tomographies (CBCTs) of 45 patients who underwent maxillo-mandibular advancement or maxillary advancement/mandibular setback were analysed by means of a semi-automated three-dimensional (3D) method; 3D models of skull and airways were produced, the latter divided into the three pharyngeal subregions. The correlation between skeletal displacement, stacked surface area and volume was investigated. The displacement threshold was identified by setting three ∆Area percentage variations. No significant difference in area and volume emerged from the comparison of the two surgical procedures with bone repositioning below the threshold (approximated to +5 mm). A threshold ranging from +4.8 to +7 mm was identified, varying in relation to the three ∆Area percentages considered. The ∆Area increased linearly above the threshold, while showing no consistency in the interval ranging from -5 mm to +5 mm.

Regional deposition of the allergens and micro-aerosols in the healthy human nasal airways

The nasal cavity is the inlet to the human respiratory system and is responsible for the olfactory sensation, filtering pollutant particulate matter, and humidifying the air. Many research studies have been performed to numerically predict allergens, contaminants, and/or drug particle deposition in the human nasal cavity; however, the majority of these investigations studied only one or a small number of nasal passages. In the present study, a series of Computed Tomography (CT) scan images of the nasal cavities from ten healthy subjects were collected and used to reconstruct accurate 3D models. All models were divided into twelve anatomical regions in order to study the transport and deposition features of different regions of the nasal cavity with specific functions. The flow field and micro-particle transport equations were solved, and the total and regional particle deposition fractions were evaluated for the rest and low activity breathing conditions. The results show that there are large variations among different subjects. The standard deviation of the total deposition fraction in the nasal cavities was the highest for 5 × 10 4 <impaction parameter (IP)< 1.125 × 10 5 with a maximum of 20%. The achieved results highlighted the nasal cavity sections that are more involved in the particle deposition. Particles with IP = 30,000 deposit more in the middle turbinate and nasopharynx areas, while for particles with IP = 300,000, deposition is mainly in the anterior parts (kiesselbach and vestibule regions). For small IP values, the amounts of deposition fractions in different regions of the nasal cavity are more uniform.

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.

Computational Fluid Dynamic Analysis of the Pharyngeal Airway after Bimaxillary Orthognathic Surgery in Patients with Mandibular Prognathism

This study aimed to analyze pharyngeal airflow characteristics and their relationship with the skeletal movement of the maxilla and mandible after bimaxillary orthognathic surgery in patients with skeletal class III (mandibular prognathism) malocclusion. Cone-beam computed tomography (CBCT) was conducted before surgery (T0), immediately after surgery (T1), and at least six months after surgery (T2). Digital imaging and communications in medicine files were transferred to InVivo (Anatomage) software to measure the skeletal changes after surgery. The changes in the maxillary and mandibular position, tongue position, and hyoid bone position were analyzed. Patient-specific models were reconstructed using 3D-Doctor software. The models after converting to the stereolithography (STL) file for Ansys integrated computer engineering and manufacturing code for computational fluid dynamics (ICEM CFD), commercial software were used for calculating the geometry, pressure drop and adjusted pressure coefficient value. The total volume of the upper airway including nasal cavity was reduced by 23% immediately after surgery and recovered to 92.2% of the initial volume six months after surgery. The airflow computation analysis showed a decrease in the pressure drop values immediately after surgery and six months after surgery. The adjusted pressure coefficients were slightly different but the change was statistically insignificant. The airflow characteristics computed using the computational fluid dynamics were correlated to the surgical changes. The surgical changes can affect the aerodynamics of the pharyngeal airway. In clinical practice, this knowledge is useful for developing a suitable orthognathic surgery treatment plan.

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

PURPOSE OF REVIEW

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

RECENT FINDINGS

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

SUMMARY

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

Extracting signature responses from respiratory flows: Low-dimensional analyses on Direct Numerical Simulation-predicted wakes of a flapping uvula

Uvula-induced snoring and associated obstructive sleep apnea is a complex phenomenon characterized by vibrating structures and highly transient vortex dynamics. This study aimed to extract signature features of uvula wake flows of different pathological origins and develop a linear reduced-order surrogate model for flow control. Six airway models were developed with two uvula kinematics and three pharynx constriction levels. A direct numerical simulation (DNS) flow solver based on the immersed boundary method was utilized to resolve the wake flows induced by the flapping uvula. Key spatial and temporal responses of the flow to uvula kinematics and pharynx constriction were investigated using continuous wavelet transform (CWT), proper orthogonal decomposition (POD), and dynamic mode decomposition (DMD). Results showed highly complex patterns in flow topologies. CWT analysis revealed multiscale correlations in both time and space between the flapping uvular and wake flows. POD analysis successfully separated the flows among the six models by projecting the datasets in the vector space spanned by the first three eigenmodes. Perceivable differences were also captured in the time evolution of the DMD modes among the six models. A linear reduced-order surrogate model was constructed from the predominant eigenmodes obtained from the DMD analysis and predicted vortex patterns from this surrogate model agreed well with the corresponding DNS simulations. The computational and analytical platform presented in this study could bring a variety of applications in breathing-related disorders and beyond. The computational efficiency of surrogate modeling makes it well suited for flow control, forecasting, and uncertainty analyses.

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.

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.

Narrowed Posterior Nasal Airway Limits Efficacy of Anterior Septoplasty

Background: Predicting symptomatic relief after septoplasty has been difficult. Minimal cross-sectional area (mCSA) measured by acoustic rhinometry and airflow resistance (R) measured by rhinomanometry have been used to select surgical candidates with mixed success. An important assumption is that mCSA and resistance are tightly coupled, but studies have reported weak or no correlation. Recently, we proposed the Bernoulli Obstruction Theory as an explanation, where tight coupling between mCSA and R is only predicted below a critical mCSA (A_crit). Methods: The nasal airway and septum of 10 healthy subjects were reconstructed from computed tomography scans. Simulated anterior septal deviations of increasing severity were created. Computational fluid dynamics simulations were performed to quantify mCSA, resistance, and flow in the healthy septum model and four simulated septal deviation models for each subject (total of 50 models). Results: A tighter coupling between mCSA and resistance was found below A_crit, estimated to be 0.20 cm² (a very severe deviation). Above A_crit, enlarging the mCSA had a smaller effect in patients with narrower cross-sectional area in the postvalve region (CSA_PV). Conclusions: Two patterns of flow increase are expected with septoplasty. Below A_crit, enlarging mCSA predictably increases flow. Above A_crit, the effect size of increasing mCSA depends on CSA_PV. Unrecognized small CSA_PV may explain persistent sensation of nasal obstruction after septoplasty. Our data suggest that inferior turbinate reduction ipsilateral to a septal deviation may amplify airflow benefits after septoplasty in patients with a narrow CSA_PV.

Side asymmetry in nasal resistance correlate with nasal obstruction severity in patients with septal deformities: Computational fluid dynamics study

OBJECTIVES

The objective of this study was to investigate the relationship between side asymmetry in nasal resistance (NR) and severity of the nasal airway obstruction (NAO) in patients with different types of nasal septal deformity (NSD).

DESIGN

Computational fluid dynamics (CFD) study.

SETTING

The study was conducted in a tertiary medical centre.

PARTICIPANTS

The study included 232 patients, who were referred to the CT examination of the paranasal sinuses. Exclusion criteria were sinonasal and respiratory diseases that may interfere with the nasal obstruction. The presence and the type of NSD were recorded according to the Mladina's classification.

MAIN OUTCOME MEASURES

The presence and severity of NAO in each patient were assessed by NOSE questionnaire. Eight computational models of the nasal cavity were created from CT scans. Models represented seven Mladina's NSD types and a straight septum of a symptomless patient. CFD calculated airflow partitioning and NR for each nasal passage. Side differences in NR were calculated by the equation ∆NR = NR_left  - NR_right . The relationship between NOSE scores, airflow partitioning and side differences in NR was explored using Spearman's correlation analysis.

RESULTS

Mladina's types of NSD showed differences in airflow partitioning and the degree of side asymmetry in NR. A significant positive correlation was detected between side differences in NR and NOSE scores (R = .762, P = .028). A significant negative correlation was found between the per cent of unilateral airflow and NR (R = -.524, P = .037).

CONCLUSIONS

Our results demonstrated that side asymmetry in NR could explain differences in NAO severity related to the NSD type.

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.

Virtual septoplasty: a method to predict surgical outcomes for patients with nasal airway obstruction

PURPOSE

A deviated nasal septum is the most common etiology for nasal airway obstruction (NAO), and septoplasty is the most common surgical procedure performed by ear-nose-throat surgeons in adults. However, quantitative criteria are rarely adopted to select patients for surgery, which may explain why up to 50% of patients report persistent or recurrent symptoms of nasal obstruction postoperatively. This study reports a systematic virtual surgery method to identify patients who may benefit from septoplasty.

METHODS

One patient with symptoms of NAO due to a septal deviation was selected to illustrate the virtual surgery concept. Virtual septoplasty was implemented in three steps: (1) determining if septal geometry is abnormal preoperatively, (2) virtually correcting the deviation while preserving the anatomical shape of the septum, and (3) estimating the post-surgical improvement in airflow using computational fluid dynamics. Anatomical and functional changes predicted by the virtual surgery method were compared to a standard septoplasty performed independently from the computational analysis.

RESULTS

A benchmark healthy nasal septum geometry was obtained by averaging the septum dimensions of 47 healthy individuals. A comparison of the nasal septum geometry in the NAO patient with the benchmark geometry identified the precise locations where septal deviation and thickness exceeded the healthy range. Good agreement was found between the virtual surgery predictions and the actual surgical outcomes for both airspace minimal cross-sectional area (0.05 cm² pre-surgery, 0.54 cm² virtual surgery, 0.50 cm² actual surgery) and nasal resistance (0.91 Pa.s/ml pre-surgery, 0.08 Pa.s/ml virtual surgery, 0.08 Pa.s/ml actual surgery).

CONCLUSIONS

Previous virtual surgery methods for NAO were based on manual edits and subjective criteria. The virtual septoplasty method proposed in this study is objective and has the potential to be fully automated. Future implementation of this method in virtual surgery planning software has the potential to improve septoplasty outcomes.