Computational fluid dynamics evaluation of posterior septectomy as a viable treatment option for large septal perforations

Computational Fluid Dynamics and Its Potential Applications for the ENT Clinician

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

On the comparison between pre- and post-surgery nasal anatomies via computational fluid dynamics

Nasal breathing difficulties (NBD) are widespread and difficult to diagnose; the failure rate of their surgical corrections is high. Computational fluid dynamics (CFD) enables diagnosis of NBD and surgery planning, by comparing a pre-operative (pre-op) situation with the outcome of virtual surgery (post-op). An equivalent comparison is involved when considering distinct anatomies in the search for the functionally normal nose. Currently, this comparison is carried out in more than one way, under the implicit assumption that results are unchanged, which reflects our limited understanding of the driver of the respiratory function. The study describes how to set up a meaningful comparison. A pre-op anatomy, derived via segmentation from a CT scan, is compared with a post-op anatomy obtained via virtual surgery. State-of-the-art numerical simulations for a steady inspiration carry out the comparison under three types of global constraints, derived from the field of turbulent flow control: a constant pressure drop (CPG) between external ambient and throat, a constant flow rate (CFR) through the airways and a constant power input (CPI) from the lungs can be enforced. A significant difference in the quantities of interest is observed depending on the type of comparison. Global quantities (flow rate, pressure drop and nasal resistance) as well as local ones are affected. The type of flow forcing affects the outcome of the comparison between pre-op and post-op anatomies. Among the three available options, we argue that CPG is the least adequate. Arguments favouring either CFR or CPI are presented.

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.

Anatomic Optical Coherence Tomography (aOCT) for Evaluation of the Internal Nasal Valve

OBJECTIVES/HYPOTHESIS

To establish the utility of anatomic optical coherence tomography (aOCT) in evaluating internal nasal valve (INV).

STUDY DESIGN

Anatomic specimen imaging study.

METHODS

Fresh-harvested human specimen heads were evaluated using both computed tomography (CT) imaging as well as using aOCT. Scans were performed at three time points: 1) After septoplasty for cartilage harvest, 2) after placement of butterfly graft (BFG), and 3) after placement of bilateral spreader grafts (SG). Imaging data were then converted into 3D models of the nasal airway. CT- and aOCT-generated models were compared by both static volumetric analysis and computational fluid dynamics (CFD) to predict nasal resistance and pressure.

RESULTS

Scans using aOCT showed comparable results to CT in terms of volumetric parameters both before and after intervention. Analysis of aOCT data by CFD demonstrated decrease in pressure after SG or BFG intervention. No statistically significant difference was observed when comparing CT- and aOCT-generated calculations of pressure or resistance.

CONCLUSION

The INV can be imaged in a static fashion using aOCT technology. Advantages over traditional CT imaging include lack of exposure to radiation and rapid scan time. In addition, in-office use is possible as aOCT technology develops. Further investigation will be necessary to define the role of aOCT in the dynamic evaluation of this vital component of the nasal airway.

LEVEL OF EVIDENCE

3 Laryngoscope, 132:2148-2156, 2022.

Is nasal airflow disrupted after endoscopic skull base surgery? A short review

Even the most delicate endonasal surgery for skull base lesion causes changes in the nasal cavity, some of them permanent. Morphological changes in the nasal cavity and their consequences (changes in nasal airflow) are often studied by advanced numerical analysis called computational fluid dynamics. This review summarizes current knowledge of endoscopic transsphenoidal skull base surgery effects on nasal airflow. Several studies have shown that endoscopic skull base surgery changes nasal anatomy to the extent that nasal airflow changes significantly postoperatively. Removing any intranasal structure increases the cross-sectional area of the respective nasal meatus, leading to increased nasal airflow in this area while airflow in the narrower periphery decreases. Middle turbinate resection increases airflow in the middle meatus and reduces airflow in the superior and inferior meatus. Small posterior septectomy does not cause a significant change in nasal airflow. Nasal septum deviation is an important factor in airflow changes. Current studies describe nasal changes after rather extensive procedures (e.g., middle turbinectomy, ethmoidectomy) that are unnecessary in routine pituitary adenoma surgery. No studies have compared changes using pre- and postoperative scans of the same patients after actual surgery.

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

OBJECTIVES/HYPOTHESIS

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

STUDY DESIGN

Blinded cohort study.

METHODS

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

RESULTS

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

CONCLUSIONS

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

LEVEL OF EVIDENCE

3 Laryngoscope, 132:509-517, 2022.

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

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

Inferior meatus augmentation procedure (IMAP) normalizes nasal airflow patterns in empty nose syndrome patients via computational fluid dynamics (CFD) modeling

BACKGROUND

Empty nose syndrome (ENS) is a controversial upper airway disorder most commonly associated with tissue loss from the inferior turbinates. The inferior meatus augmentation procedure (IMAP) has been shown to effectively reduce ENS symptoms in a durable manner, but the precise mechanisms that may govern this symptomatic improvement remain unknown.

METHODS

Five patients with ENS who underwent bilateral IMAP via submucosal costal cartilage implant were assessed. Pre-implant and 6 months post-implant computed tomography (CT) imaging for each ENS patient was analyzed in a blinded fashion using computational fluid dynamics (CFD) modeling to investigate intrapatient changes in airflow parameters.

RESULTS

Following surgery, ENS patients have significantly improved symptoms as indexed by Empty Nose Syndrome 6-Item Questionnaire (ENS6Q) scoring (pre-implant: 14.00 ± 4.06 [mean ± standard deviation]; 95% confidence interval [CI], 10.44 to 17.56; post-implant: 4.8 ± 2.77; 95% CI, 2.37 to 7.23; Cohen's d = 2.64; p = 0.02). Using CFD, a significant shift in nasal airflow patterns was observed, where airflow deviates away from the middle meatus upon hitting the implant (pre-implant: 67.13% ± 11.14%; 95% CI, 60.22% to 74.04%; post-implant: 46.18% ± 12.81%; 95% CI, 38.23% to 54.12%; d = 1.74; p < 0.05) toward the inferior meatus (pre-implant: 30.55% ± 11.29%; 95% CI, 23.55% to 37.55%; post-implant: 42.59% ± 9.60%; 95% CI, 36.63 to 48.54%; d = 1.14; p < 0.05). No significant changes were found in nasal resistance (pre-implant: 0.102 ± 0.015; 95% CI, 0.092 to 0.112 Pa*s/mL; post-implant: 0.105 ± 0.041; 95% CI, 0.081 to 0.130 Pa*s/mL). In addition, the improvement of ENS6Q scoring significantly correlated with percent reduction in aberrant airflow through the middle meatus (R² = 0.60, p = 0.04).

CONCLUSION

This study supports our prior working hypothesis that disordered vectors of nasal airflow congregate in the middle meatus contribute to ENS symptoms, not nasal resistance. Moreover, these data illuminate a paradoxical, but consistent, restoration of nasal airflow to the inferior meatus following the replacement of turbinate tissue volume in the inferior meatus via IMAP surgery, potentially due to the Coandă effect.

Regional Peak Mucosal Cooling Predicts Radiofrequency Treatment Outcomes of Nasal Valve Obstruction

OBJECTIVES/HYPOTHESIS

Low energy radiofrequency may offer effective treatment for narrow or obstructed nasal valve, yet its precise mechanism is not fully understood.

STUDY DESIGN

Prospective, nonrandomized, case series.

METHODS

Twenty prospective patients with internal nasal valve obstruction underwent office-based Vivaer treatment (Aerin Medical, Inc) under local anesthesia. Computational fluid dynamics (CFD) models were constructed based on the pre- and 90 days post-procedure computed tomography (CT) scans to identify salient changes in nasal airflow parameters.

RESULTS

Patients' Nasal Obstruction Symptom Evaluation score (NOSE: pre-treatment 78.89 ± 11.57; post-treatment 31.39 ± 18.30, P = 5e-7) and Visual Analog Scale of nasal obstruction (VAS: pre-treatment 6.01 ± 1.83; post-treatment 3.44 ± 2.11, P = 1e-4) improved significantly at 90 days after the minimally invasive approach. Nasal airway volume in the treatment area increased ~7% 90 days post-treatment (pre-treatment 5.97 ± 1.20, post-treatment 6.38 ± 1.50 cm³ , P = .018), yet there were no statistically significant changes in the measured peak nasal inspiratory flowrate (PNIF, pre-treatment: 60.16 ± 34.49; post-treatment: 72.38 ± 43.66 ml/s; P = .13) and CFD computed nasal resistance (pre-treatment: 0.096 ± 0.065; post-treatment: 0.075 ± 0.026 Pa/(ml/s); P = .063). As validation, PNIF correlated significantly with nasal resistance (r = 0.47, P = .004). Among all the variables, only the peak mucosal cooling posterior to the nasal vestibule significantly correlated with the NOSE at baseline (r = -0.531, P = .023) and with post-treatment improvement (r = 0.659, P = .003).

CONCLUSION

Minimal remodeling of the nasal valve (7% in this study) may have a profound effect on perceived nasal obstruction, despite little effect on nasal resistance, or PNIF. The results corroborated our previous findings that subjective relief of nasal obstruction correlates with regional mucosal cooling rather than nasal resistance or peak flow rate, a potential target for future effective, personalized therapeutic approaches.

LEVEL OF EVIDENCE

4 Laryngoscope, 131:E1760-E1769, 2021.

Optimizing septal perforation repair techniques

PURPOSE OF REVIEW

Multiple successful techniques and approaches for nasal septal perforation repair have been described, yet consistency in perforation and outcome metrics is required to identify the optimal approach to repair. The present article will review the recent literature.

RECENT FINDINGS

Computational fluid dynamic studies continue to expand our understanding of the airflow dynamics in nasal septal perforation and after repair. Combining rhinoplasty and nasal septal perforation repair in appropriately selected patients can be safely done with excellent results. There has been a rise in utilization of a temporoparietal fascia with polydiaxonone plate construct for septal perforation repair with excellent outcomes.

SUMMARY

The present review provides the reconstructive surgeon with an update on nasal septal perforation repair and describes a recently popularized technique of temporoparietal fascia-polydiaxonone plate for perforation reconstruction.

Correlation of sinonasal symptoms with the size and position of nasal septal perforations

OBJECTIVE

To assess the correlation of sinonasal symptoms and quality of life with size and position of nasal septal perforation (NSP).

METHOD

This is a prospective observational study in a tertiary referral center involving adult patients presented with NSP. The Sino-Nasal Outcome Test-22 (SNOT-22) and its clinico-psychometric domains were analyzed, including additional NSP-specific symptoms (nasal crusting, epistaxis, and whistling noise during nasal breathing). The size of NSP was measured radiologically by calculating the area in cm² and anteroposterior (AP) diameter. Position of perforation was determined clinically by distance from columella to the anterior edge of the perforation.

RESULTS

Forty patients were included in this study (22 males). The most common etiology of NSP was iatrogenic. The mean SNOT-22 score was 50.8 (standard deviation 23.8), and mean NSP size was 3.0 cm² . No statistical correlation was observed between the total SNOT-22 score with either position or size of NSP. As for NSP-specific symptoms, there was a statistically significant negative correlation with the size of perforation (A-P diameter) (r = -0.34, P = 0.03) and position of the perforation (r = -0.49, P = 0.0016), suggesting that these symptoms improved with posterior and larger perforations.

CONCLUSION

Reported SNOT-22 scores were comparable to patients with recalcitrant chronic rhinosinusitis, although the scores did not correlate with size and position of NSP. Sinonasal symptoms typically observed in NSP improved with more posteriorly placed and larger sized perforations. This study provides an insight into the quality-of-life impact of NSP and affirms the clinical observation that anterior NSP are more symptomatic.

LEVEL OF EVIDENCE

4 Laryngoscope, 2019.

Peak Sinus Pressures During Sneezing in Healthy Controls and Post-Skull Base Surgery Patients

OBJECTIVES/HYPOTHESIS

Patients are frequently advised to sneeze with an open mouth and avoid nose-blowing following an endoscopic endonasal approache (EEA) to the skull base, despite a lack of quantitative evidence. This study applies computational fluid dynamics (CFD) to quantify sinus pressures along the skull base during sneezing.

STUDY DESIGN

Case-control series.

METHODS

Computed tomography or magnetic resonance imaging scans of four post-EEA patients and four healthy controls were collected and analyzed utilizing CFD techniques. A pressure drop of 6,000 Pa was applied to the nasopharynx based on values in the literature to simulate expiratory nasal airflow during sneezing. Peak pressures along the skull base in frontal, ethmoid, and sphenoid sinuses were collected.

RESULTS

Significant increases in skull base peak pressure was observed during sneezing, with significant individual variations from 2,185 to 5,685 Pa. Interestingly, healthy controls had significantly higher pressures compared to post-EEA patients (5179.37 ± 198.42 Pa vs. patients 3,347.82 ± 1,472.20 Pa, P < .05), which could be related to higher anterior nasal resistance in unoperated healthy controls (0.44 ± 0.22 vs. 0.31 ± 0.16 Pa/mL/sec for patients, P = .38). The sinus pressure buildup may be due to airway resistance functioning as a valve preventing air from being released quickly. Supporting this theory, there was a strong correlation (r = 0.82) between peak skull base pressure and the ratio of anterior resistance to total resistance. Within-subject variation in pressures between different skull base regions was much lower (average = ~5%).

CONCLUSIONS

This study provided the first quantitative analysis of air pressure along the skull base during sneezing in post-EEA patients through CFD, suggesting that pressure buildup may depend on individual anatomy.

LEVEL OF EVIDENCE

3b Laryngoscope, 130:2138-2143, 2020.

Functional relevance of computational fluid dynamics in the field of nasal obstruction: A literature review

BACKGROUND

Nasal airway obstruction (nasal obstruction) is a common symptom affecting the quality of life of patients. It can be estimated by patient perception or physical measurements. Computational fluid dynamics (CFD) can be used to analyse nasal ventilation modalities. There is a lack of comparative studies investigating the correlations between CFD variables and patient perception or physical measurements.

OBJECTIVE OF THE REVIEW

Our goal was to define correlations between CFD variables and patient perception and physical measurements. We also aimed to identify the most reliable CFD variable (heat flux, WSS, total pressure, temperature…) characterising nasal breathing perception.

TYPE OF REVIEW

Systematic literature review using PRISMA guidelines.

SEARCH STRATEGY

The selected studies were obtained from the US National Library of Medicine (PubMed) online database, MEDLINE (Ovid), Google Scholar and the Cochrane Library using a combination of MeSH terms (nose, paranasal sinus, fluid dynamics, rhinology) and non-MeSH terms (CFD, nasal airway, nasal airflow, numerical, nasal symptoms). Studies that did not incorporate objective or subjective clinical assessment were excluded.

EVALUATION METHOD

We compared all results obtained by authors regarding CFD variables and assessment of nasal airway obstruction (clinical or physical).

RESULTS

To compare nasal obstruction with CFD variables, most authors use CFD-calculated nasal resistances, airflow, heat flux, wall shear stress, total pressure, velocities and streamlines. We found that heat flux appears to be the CFD variable most closely correlated with patient perception. Total pressure, wall shear stress and velocities are also useful and show good correlations. Correlations between CFD-calculated nasal resistances and patient perception are stronger after correction of the nasal cycle.

CONCLUSIONS

The growing number of CFD studies on the nose has led to a better understanding of nasal obstruction. The clinical interpretation of previously unknown data, such as WSS and heat flux, is opening up new horizons in the understanding of this symptom. Heat fluxes are among the best CFD values correlated with patient perception. More studies need to be performed including temperature and humidity exchanges.

Asymptomatic vs symptomatic septal perforations: a computational fluid dynamics examination

BACKGROUND

A nasal septal perforation (NSP) can lead to frustrating symptoms for some patients while remaining completely asymptomatic for others, without a clear mechanism differentiating them.

METHODS

We applied individual computed tomography (CT)-based computational fluid dynamics (CFD) to examine the nasal aerodynamics differences between 5 asymptomatic and 15 symptomatic NSP patients. Patients' symptoms were confirmed through interviews, 22-item Sino-Nasal Outcome Test score (asymptomatic, 25 ± 18.8; symptomatic, 53.7 ± 18.2), nasal obstruction symptom evaluation score (asymptomatic, 28.0 ± 32.1; symptomatic, 62.2 ± 32.2), and review of medical history.

RESULTS

No statistical differences were found in perforation location, size (asymptomatic, 1.94 ± 1.88 cm² ; symptomatic, 1.36 ± 1.44 cm² ), nasal resistance (asymptomatic, 0.059 ± 0.012 Pa·s/mL; symptomatic, 0.063 ± 0.022 Pa·s/mL), and computed flow rate shunting across the perforation (asymptomatic, 52.9 ± 30.9 mL/s; symptomatic, 27.4 ± 23.6 mL/s; p > 0.05). However, symptomatic patients had significantly higher wall shear stress (WSS) and heat flux, especially along the posterior perforation margin (WSS, 0.54 ± 0.12 vs 1.15 ± 0.49 Pa, p < 0.001; heat flux, 0.21 ± 0.05 vs 0.37 ± 0.14 W/cm² , p < 0.01). A WSS cutoff at 0.72 Pa can separate asymptomatic vs symptomatic NSP with 87% sensitivity and 100% specificity. Flow visualization showed flow peaks toward the posterior margin that may be responsible for the high WSS and heat flux among symptomatic NSPs.

CONCLUSION

This study is the first CFD examination of asymptomatic and symptomatic NSP with regional aerodynamics and stress abnormalities, beyond size or location, being implicated as the mechanism behind the symptomology of NSP. This finding could serve as an objective basis for future personalized treatment decisions and optimization.

Computational fluid dynamic analysis of aggressive turbinate reductions: is it a culprit of empty nose syndrome?

BACKGROUND

Empty nose syndrome (ENS) remains highly controversial, with aggressive inferior turbinate reduction (ITR) or mucociliary dysfunction frequently implicated. However, the appropriate degree of ITR is highly debatable.

METHODS

We applied individual computed tomography (CT)-based computational fluid dynamics (CFD) to 5 patients receiving relatively aggressive ITR but with no ENS symptoms, and compared them to 27 symptomatic ENS patients who all had histories of aggressive ITRs, and 42 healthy controls. Patients' surgical outcomes were confirmed with 22-item Sino-Nasal Outcome Test (SNOT-22) (ITR: 6.40 ± 4.56; ENS: 58.2 ± 15.9; healthy: 13.2 ± 14.9), Nasal Obstruction Symptom Evaluation (NOSE) scores (ITR: 4.00 ± 2.24; ENS: 69.4 ± 17.1; healthy: 11.9 ± 12.9), and Empty Nose Syndrome 6-Item Questionnaire (ENS6Q) (≥11 for ENS).

RESULTS

Both aggressive ITR without ENS symptoms and symptomatic ENS patients had significantly lower nasal resistance (ITR: 0.059 ± 0.020 Pa·s/mL; ENS: 0.052 ± 0.015 Pa·s/mL; healthy: 0.070 ± 0.021 Pa·s/mL) and higher cross-sectional areas surrounding the inferior turbinate (ITR: 0.94 ± 0.21 cm² ; ENS: 1.19 ± 1.05 cm² ; healthy: 0.42 ± 0.22 cm² ) than healthy controls. The lack of significant differences among patient groups indicated similar degrees of surgeries between ITR with and without ENS symptom cohorts. However, symptomatic ENS patients have paradoxical significantly less airflow in the inferior meatus (ITR: 47.7% ± 23.6%; ENS: 25.8% ± 17.6%; healthy: 36.5 ± 15.9%; both p < 0.01), but higher airflow around the middle meatus (ITR: 49.7% ± 22.6%; ENS: 66.5% ± 18.3%; healthy: 49.9% ± 15.1%, p < 0.0001) than aggressive ITR without symptoms and controls. Aggressive ITR patients have increased inferior meatus airflow as expected (p < 0.05). This imbalanced airflow produced less inferior wall-shear-stress distribution among symptomatic ENS patients only (ITR: 42.45% ± 11.4%; ENS: 32.2% ± 12.6%; healthy: 49.7% ± 9.9%). ENS patients (n = 12) also had impaired nasal trigeminal function, as measured by menthol lateralization detection thresholds (ITR: 15.2 ± 1.2; ENS: 10.3 ± 3.9; healthy: 13.8 ± 3.09, both p < 0.0001). Surprisingly, aggressive ITR patients without ENS symptoms have better menthol lateralization detection thresholds (LDTs) than healthy controls.

CONCLUSION

Although turbinate tissue loss is linked with ENS, the degree of ITR that might distinguish postoperative patient satisfaction in their nasal breathing vs development of ENS symptoms is unclear. Our results suggest that a combination of distorted nasal aerodynamics and loss of mucosal sensory function may potentially lead to ENS symptomology.

Nasal Structural and Aerodynamic Features That May Benefit Normal Olfactory Sensitivity

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

Impact of Middle Turbinectomy on Airflow to the Olfactory Cleft: A Computational Fluid Dynamics Study

BACKGROUND

The impact of middle turbinate resection (MTR) on olfaction remains a point of debate in the current literature. Few studies have objectively evaluated olfactory cleft airflow following MTR; thus, the mechanism by which MTR may impact olfaction is poorly understood. It is not known whether the postsurgical changes in airway volume, flow, and resistance increase odorant transport or disrupt the patterns of normal airflow. Computational fluid dynamics can be used to study the nasal airway and predict responses to surgical intervention.

OBJECTIVE

To evaluate the functional impact of MTR on nasal airflow, resistance, and olfaction.

METHODS

Five maxillofacial computed tomography scans of patients without signs of significant sinusitis or nasal polyposis were used. Control models for each patient were compared to their corresponding model after virtual total MTR. For each model, nasal airway volume, nasal resistance, and air flow rate were determined. Odorant transport of 3 different odorants in the nasal cavity was simulated based on the computed steady airflow field.

RESULTS

Total airflow significantly increased following bilateral MTR in all patient models ( P < .05). Consistent with our airflow results, we found a decrease in nasal resistance following MTR. MTR significantly increased area averaged flux to the olfactory cleft when compared to controls for phenylethyl alcohol (high-sorptive odorant). Results for carvone (medium sorptive) were similarly elevated. MTR impact on limonene, a low flux odorant, was equivocal.

CONCLUSION

MTR increases nasal airflow while decreasing the nasal resistance. Overall, olfactory flux increased for high sorptive (phenylethyl alcohol) and medium sorpitve (l-carvone) odorants. However, the significant variation observed in one of our models suggests that the effects of MTR on the nasal airflow and the resultant olfaction can vary between individuals based on individual anatomic differences.

Computational fluid dynamics after endoscopic endonasal skull base surgery-possible empty nose syndrome in the context of middle turbinate resection

BACKGROUND

Empty nose syndrome (ENS) is a rare and debilitating disease with a controversial definition, etiology, and treatment. One puzzling fact is that patients who undergo an endoscopic endonasal approach (EEA) often have resection of multiple anatomic structures, yet seldom develop ENS. In this pilot study, we analyzed and compared the computational fluid dynamics (CFD) and symptoms among post-EEA patients, ENS patients, and healthy subjects.

METHODS

Computed tomography scans of 4 post-EEA patients were collected and analyzed using CFD techniques. Two patients had significant ENS symptoms based on results of the Empty Nose Syndrome 6-item Questionnaire (score >11), whereas the other 2 were asymptomatic. As a reference, their results were compared with previously published CFD results of 27 non-EEA ENS patients and 42 healthy controls.

RESULTS

Post-EEA patients with ENS symptoms had a similar nasal airflow pattern as non-EEA ENS patients. This pattern differed significantly from that of EEA patients without ENS symptoms and healthy controls. Overall, groups with ENS symptoms exhibited airflow dominant in the middle meatus region and a significantly lower percentage of airflow in the inferior turbinate region (EEA with ENS, 17.74 ± 4.00% vs EEA without ENS, 51.25 ± 3.33% [t test, p < 0.02]; non-EEA ENS, 25.8 ± 17.6%; healthy subjects, 36.5 ± 15.9%) as well as lower peak wall shear stress (EEA with ENS, 0.30 ± 0.13 Pa vs EEA without ENS, 0.61 ± 0.03 Pa [p = 0.003]; non-EEA ENS, 0.58 ± 0.24 Pa; healthy subjects, 1.18 ± 0.81 Pa).

CONCLUSION

These results suggest that turbinectomy and/or posterior septectomy may have a varying functional impact and that ENS symptoms go beyond anatomy and correlate with aerodynamic changes. The findings open the door for CFD as a potential objective diagnosis tool for ENS.

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.

Modeling congenital nasal pyriform aperture stenosis using computational fluid dynamics

OBJECTIVES

Congenital nasal pyriform aperture stenosis (CNPAS) is a rare cause of airway obstruction in the neonate. Computational airway modeling has not been done in neonates and young infants to understand the impact of stenosis on functional nasal airflow. In this study, we 1) applied computational fluid dynamics (CFD) model to the airway of a neonate with CNPAS and 2) compare airflow dynamics of a normal and CNPAS airway.

METHODS

Three-dimensional models of the nasal airway of a normal neonate and a neonate with CNPAS were created using computed tomography scans of the facial bones. Measured anatomic parameters included volume, surface area, and cross-sectional area. CFD simulation was then performed. Simulated flow parameters included pressure, average velocity, and resistance.

RESULTS

The neonate with CNPAS had a lesser volume (2.74 cm³ vs. 4.50 cm³) and surface area (18.8 cm² vs. 45.5 cm²) than the normal airway. The CNPAS airway had a lesser bilateral cross-sectional area and average cross-sectional velocity throughout the length of the model. While there is a large pressure drop in the normal airway immediately after the entry point, the pressure drop in the CNPAS airway occurs more posteriorly. The total nasal resistance was approximately eight-fold greater in the CNPAS airway than the normal.

CONCLUSIONS

CFD analysis can be performed on airways of neonates with nasal obstruction, such as in CNPAS. A CFD model may help characterize severity of airway obstruction as it can predict the three-dimensional pattern of airflow. Determining the role of CFD in clinical management of CNPAS requires further investigation.

Quantification of tissue-engineered trachea performance with computational fluid dynamics

OBJECTIVES/HYPOTHESIS

Current techniques for airway characterization include endoscopic or radiographic measurements that produce static, two-dimensional descriptions. As pathology can be multilevel, irregularly shaped, and dynamic, minimal luminal area (MLA) may not provide the most comprehensive description or diagnostic metric. Our aim was to examine the utilization of computational fluid dynamics (CFD) for the purpose of defining airway stenosis using an ovine model of tissue-engineered tracheal graft (TETG) implantation.

STUDY DESIGN

Animal research model.

METHODS

TETGs were implanted into sheep, and MLA was quantified with imaging and endoscopic measurements. Graft stenosis was managed with endoscopic dilation and stenting when indicated. Geometries of the TETG were reconstructed from three-dimensional fluoroscopic images. CFD simulations were used to calculate peak flow velocity (PFV) and peak wall shear stress (PWSS). These metrics were compared to values derived from a quantitative respiratory symptom score.

RESULTS

Elevated PFV and PWSS derived from CFD modeling correlated with increased respiratory symptoms. Immediate pre- and postimplantation CFD metrics were similar, and implanted sheep were asymptomatic. Respiratory symptoms improved with stenting, which maintained graft architecture similar to dilation procedures. With stenting, baseline PFV (0.33 m/s) and PWSS (0.006 Pa) were sustained for the remainder of the study. MLA measurements collected via bronchoscopy were also correlated with respiratory symptoms. PFV and PWSS found via CFD were correlated (R²  = 0.92 and 0.99, respectively) with respiratory symptoms compared to MLA (R²  = 0.61).

CONCLUSIONS

CFD is valid for informed interventions based on multilevel, complex airflow and airway characteristics. Furthermore, CFD may be used to evaluate TETG functionality.

LEVEL OF EVIDENCE

NA. Laryngoscope, E272-E279, 2018.

Investigation of the abnormal nasal aerodynamics and trigeminal functions among empty nose syndrome patients

BACKGROUND

Abnormal nasal aerodynamics or trigeminal functions have been frequently implicated in the symptomology of empty nose syndrome (ENS), yet with limited evidence.

METHODS

Individual computed tomography (CT)-based computational fluid dynamics (CFD) was applied to 27 ENS patients to simulate their nasal aerodynamics and compared with 42 healthy controls. Patients' symptoms were confirmed with Empty Nose Syndrome 6-item Questionnaire (ENS6Q), 22-item Sino-Nasal Outcome Test (SNOT-22), and Nasal Obstruction Symptom Evaluation (NOSE) scores. Nasal trigeminal sensitivity was measured with menthol lateralization detection thresholds (LDTs).

RESULTS

ENS patients had significantly lower (∼25.7%) nasal resistance and higher (∼2.8 times) cross-sectional areas compared to healthy controls (both p < 0.001). Despite inferior turbinate reductions, CFD analysis demonstrated that ENS patients had increased airflow concentrated in the middle meatus region (66.5% ± 18.3%) compared to healthy controls (49.9% ± 15.1%, p < 0.0001). Significantly less airflow (25.8% ± 17.6%) and lower peak wall shear stress (WSS) (0.58 ± 0.24 Pa) were found in the inferior meatus (vs healthy: 36.5% ± 15.9%; 1.18 ± 0.81 Pa, both p < 0.05), with the latter significantly correlated with the symptom scores of ENS6Q (r = -0.398, p = 0.003). Item-wise, complaints of "suffocation" and "nose feels too open" were also found to be significantly correlated with peak WSS around the inferior turbinate (r = -0.295, p = 0.031; and r = -0.388, p = 0.004, respectively). These correlations were all negative, indicating that less air-mucosal stimulations resulted in worse symptom scores. ENS patients (n = 12) also had impaired menthol LDT when compared to healthy controls (p < 0.0001).

CONCLUSION

This is the first CFD examination of nasal aerodynamics in a large cohort of ENS patients. The results indicated that a combination of loss of neural sensitivity and poorer inferior air-mucosal stimulation may potentially lead to ENS symptomology.

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

The human nose serves vital physiological functions, including warming, filtration, humidification, and olfaction. These functions are based on transport phenomena that depend on nasal airflow patterns and turbulence. Accurate prediction of these airflow properties requires careful selection of computational fluid dynamics models and rigorous validation. The validation studies in the past have been limited by poor representations of the complex nasal geometry, lack of detailed airflow comparisons, and restricted ranges of flow rate. The objective of this study is to validate various numerical methods based on an anatomically accurate nasal model against published experimentally measured data under breathing flow rates from 180 to 1100ml/s. The numerical results of velocity profiles and turbulence intensities were obtained using the laminar model, four widely used Reynolds-averaged Navier-Stokes (RANS) turbulence models (i.e., k-ε, standard k-ω, Shear Stress Transport k-ω, and Reynolds Stress Model), large eddy simulation (LES) model, and direct numerical simulation (DNS). It was found that, despite certain irregularity in the flow field, the laminar model achieved good agreement with experimental results under restful breathing condition (180ml/s) and performed better than the RANS models. As the breathing flow rate increased, the RANS models achieved more accurate predictions but still performed worse than LES and DNS. As expected, LES and DNS can provide accurate predictions of the nasal airflow under all flow conditions but have an approximately 100-fold higher computational cost. Among all the RANS models tested, the standard k-ω model agrees most closely with the experimental values in terms of velocity profile and turbulence intensity.