Influence of Mesh Density on Airflow and Particle Deposition in Sinonasal Airway Modeling

In silico investigation of inhalation condition impacts on hygroscopic growth and deposition of salbutamol sulphate in human airways

The objective of this study is to explore the transport, size growth, and deposition of Salbutamol Sulphate (SS) using Computational Fluid Dynamics (CFD). A CT-based realistic model of human airways from the oral cavity to the 5th generation of the lung was utilized as the computational domain. Four Test Cases (TC) with varying temperature and relative humidity (RH) under two inspiratory waveforms were considered to completely evaluate the impact of inhalation conditions on particle growth. Salbutamol Sulphate (SS) is a β2-adrenergic agonist and has been extensively used for asthma treatment. A monodispersed distribution of SS particles with an initial diameter of 167 nm was considered at the mouth inlet based on pharmaceutical data. Results indicated that inhalation of saturated/supersaturated air (RH>100%) leads to significant hygroscopic growth of SS particles with a factor of 10. In addition, the deposition efficiency of SS particles under the Quick and Deep (QD) inhalation profile was enhanced as the flow temperature and humidity increased. However, the implementation of Slow and Deep (SD) inspiratory waveform revealed that the same particle size growth is achieved in the respiratory system with lower deposition efficiency in the mouth-throat (less than 3%) and tracheobronchial airway (less than 2.18%). For the escaped particles form the right lung, in the SD waveform under TC 3, the maximum particle size distribution was for 600 nm particles with 25% probability. In the left lung, 30% of the particles were increased up to 950 nm in size. For the QD waveform in TC 3 and TC4, the most frequent particles were 800 nm with 36% probability. This holds practical significance in the context of deep lung delivery for asthmatic patients with enhanced deposition efficiency and large particle size. The findings of the present study can contribute to the development of targeted drug delivery strategies for the treatment of pulmonary diseases using hygroscopic dry powder formulations.

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

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

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

BACKGROUND

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

METHODS

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

FINDINGS

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

INTERPRETATION

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

Computational Rhinology: Unraveling Discrepancies between In Silico and In Vivo Nasal Airflow Assessments for Enhanced Clinical Decision Support

Computational rhinology is a specialized branch of biomechanics leveraging engineering techniques for mathematical modelling and simulation to complement the medical field of rhinology. Computational rhinology has already contributed significantly to advancing our understanding of the nasal function, including airflow patterns, mucosal cooling, particle deposition, and drug delivery, and is foreseen as a crucial element in, e.g., the development of virtual surgery as a clinical, patient-specific decision support tool. The current paper delves into the field of computational rhinology from a nasal airflow perspective, highlighting the use of computational fluid dynamics to enhance diagnostics and treatment of breathing disorders. This paper consists of three distinct parts-an introduction to and review of the field of computational rhinology, a review of the published literature on in vitro and in silico studies of nasal airflow, and the presentation and analysis of previously unpublished high-fidelity CFD simulation data of in silico rhinomanometry. While the two first parts of this paper summarize the current status and challenges in the application of computational tools in rhinology, the last part addresses the gross disagreement commonly observed when comparing in silico and in vivo rhinomanometry results. It is concluded that this discrepancy cannot readily be explained by CFD model deficiencies caused by poor choice of turbulence model, insufficient spatial or temporal resolution, or neglecting transient effects. Hence, alternative explanations such as nasal cavity compliance or drag effects due to nasal hair should be investigated.

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Olfactory drug delivery with intranasal sprays after nasal midvault reconstruction

Conductive olfaction and nose to brain drug delivery are important processes that remain limited by inadequate odorant or drug delivery to the olfactory airspace. Primary challenges include anatomic barriers and poor targeting to the olfactory region. This study uses computational fluid dynamics to investigate the effects of nasal midvault surgery on olfactory drug delivery with intranasal sprays. Soft tissue elevation, spreader flaps, and spreader grafts were performed on two fresh cadaveric specimens, using computed tomography for airway reconstruction. Nasal airflow and drug particle transport simulations were performed under these conditions: inhalation rate (15, 30 L/min), spray velocity (1, 5, 10 m/s), spray location (top, bottom, center, medial, lateral), head position (upright, supine, forward, backward), and particle size (1-100 µm). Simulation results were used to calculate drug particle deposition to the olfactory airspaces and bulbs. Total olfactory deposition was < 5% but attained a maximum of 36.33% when sorted by particle size. There was no association between nasal midvault surgery and olfactory deposition. No single parameter or technique demonstrated superior olfactory deposition, but smaller particle size, slower spray velocity, and higher inhalation rate tended to optimize olfactory deposition, providing important implications for future intranasal spray and drug design to target the olfactory airspace.

Computational investigation of flow characteristics and particle deposition patterns in a realistic human airway model under different breathing conditions

Airborne particle pollution causes a range of respiratory and cardiovascular disorders by entering the human respiratory system through the breathing process. The administration of pharmaceutical particles by inhalation is another effective way to treat pulmonary illnesses. Studying particle deposition in the respiratory system during human breathing is crucial to maintaining human health. This necessity served as the impetus for this work, which aims to investigate how the airflow and particles' deposition are influenced by constant inhalation and circulatory breathing, particle diameter, and changes in airflow rate. The focus of this paper is to compare the particle deposition results of circulatory respiration with constant respiration. Based on computed tomography (CT) scan pictures, a precise human airway model from the mouth cavity to the fifth-generation bronchi was created. Flow fields and particle deposition inside the respiratory tract were examined at varied breathing rates (30, 60, and 90L/min of constant and circulatory breathing) and varying haled particle sizes (5 and 10 μm). The results showed that the oropharyngeal area is often where the majority of particles are deposited. The particle distribution fraction is more significant in the bronchial area than the oropharyngeal region due to lower inhalation velocities and smaller particle sizes. For particles with a diameter of 5µm, constant respiration and circulatory respiration have virtually identical particle distribution fractions in each region. For particles with a diameter of 10µm, the particle distribution fraction for circulatory respiration is slightly higher than for constant respiration in the bronchial region as the flow rate increases. For both constant and circulatory respiration, particles are deposited more in the right lung and less in the left. These results contribute to further research on respiratory diseases caused by inhaled particles and guide inhalation therapy for better treatment outcomes.

Intranasal Spray Characteristics for Best Drug Delivery in Patients With Chronic Rhinosinusitis

OBJECTIVES

To determine parameter combinations for effective drug delivery of intranasal spray steroids to the ostiomeatal complex (OMC) and maxillary sinus (MS) in patients with chronic rhinosinusitis (CRS).

METHODS

Each patient's sinonasal cavity was reconstructed from computed tomography scans. Intranasal airflow and drug particle transport were simulated using computational fluid dynamic modeling. Airflow simulations were performed at 15 Pascal inhalation pressure. Intranasal spray particles of 1-100 μm were simulated at release speeds of 1, 5, and 10 m/s from 6 release locations (Bottom, Center, Top, Lateral, Lateral-Bottom, and Lateral-Top) at a nozzle insertion depth of 15 mm. Drug delivery simulations were performed in the head tilted forward position.

RESULTS

Maximal OMC deposition was 0.78%-12.44%, while maximal MS deposition was 0.02%-1.03% across all simulations. In general, particles between 6 and 10 μm had the best OMC (at 1 m/s particle velocity) and MS (at 10 m/s particle velocity) deposition. Particles ranging from 21 to 30 μm also had superior OMC deposition. The lateral and lateral-top spray release locations produced maximum OMC deposition, but no one release location demonstrated an increase in MS deposition.

CONCLUSION

This preliminary study suggests that it is challenging to determine a common set of intranasal spray parameter combinations for effective drug delivery to the OMC and MSs. Although drug particle size and spray particle velocity seem to impact particle deposition patterns, spray release location appears to vary with anatomical differences between subjects, particularly when the MS is the target location for particle deposition. Laryngoscope, 133:1036-1043, 2023.

A computational analysis on the impact of multilevel laryngotracheal stenosis on airflow and drug particle dynamics in the upper airway

Laryngotracheal stenosis (LTS) is a type of airway narrowing that is frequently caused by intubation-related trauma. LTS can occur at one or multiple locations in the larynx and/or trachea. This study characterizes airflow dynamics and drug delivery in patients with multilevel stenosis. Two subjects with multilevel stenosis (S1 = glottis + trachea, S2 = glottis + subglottis) and one normal subject were retrospectively selected. Computed tomography scans were used to create subject-specific upper airway models. Computational fluid dynamics modeling was used to simulate airflow at inhalation pressures of 10, 25, and 40 Pa, and orally inhaled drug transport with particle velocities of 1, 5, and 10 m/s, and particle size range of 100 nm-40 µm. Subjects had increased airflow velocity and resistance at stenosis with decreased cross-sectional area (CSA): S1 had the smallest CSA at trachea (0.23 cm2) and resistance = 0.3 Pa·s/mL; S2 had the smallest CSA at glottis (0.44 cm2), and resistance = 0.16 Pa·s/mL. S1 maximal stenotic deposition was 4.15% at trachea; S2 maximal deposition was 2.28% at glottis. Particles of 11-20 µm had the greatest deposition, 13.25% (S1-trachea) and 7.81% (S2-subglottis). Results showed differences in airway resistance and drug delivery between subjects with LTS. Less than 4.2% of orally inhaled particles deposited at stenosis. Particle sizes with most stenotic deposition were 11-20 µm and may not represent typical particle sizes emitted by current-use inhalers.

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.

Computational fluid dynamics model of laryngotracheal stenosis and correlation to pulmonary function measures

3D models of airway lumens were created from CT scans of 19 patients with laryngotracheal stenosis. Computational fluid dynamics (CFD) simulations were completed for each, and results were compared to measured peak inspiratory flow rate, grade of lumen constriction, and measures of airway geometry. Results demonstrate flow resistance and shear stress correlate with degree of lumen constriction and absolute cross-sectional area as well as flow rate. Flow recirculation depends on airway constriction but does not vary with flow rate. Resistance and wall shear stress did not correlate well with functional measures. Flow recirculation did differ between subjects with higher functional measures and subjects with lower functional measures. This analysis provides mathematical models to predict airway resistance, wall shear stress, and flow reversal according lumen constriction and inspiratory flow rate. It suggests aerodynamic factors such as flow recirculation play a role in differences in functional performance between patients with similar airway measures.

Comparison of Inhaled Drug Delivery in Patients With One- and Two-level Laryngotracheal Stenosis

OBJECTIVES/HYPOTHESIS

Laryngotracheal stenosis (LTS) is a functionally devastating condition with high respiratory morbidity and mortality. This preliminary study investigates airflow dynamics and stenotic drug delivery in patients with one- and two-level LTS.

STUDY DESIGN

A Computational Modeling Restropective Cohort Study.

METHODS

Computed tomography scans from seven LTS patients, five with one-level (three subglottic, two tracheal), and two with two-level (glottis + trachea, glottis + subglottis) were used to reconstruct patient-specific three-dimensional upper airway models. Airflow and orally inhaled drug particle transport were simulated using computational fluid dynamics modeling. Drug particle transport was simulated for 1-20 μm particles released into the mouth at velocities of 0 m/s, 1 m/s, 3 m/s, and 10 m/s for metered dose inhaler (MDI) and 0 m/s for dry powder inhaler (DPI) simulations. Airflow resistance and stenotic drug deposition in the patients' airway models were compared.

RESULTS

Overall, there was increased airflow resistance at stenotic sites in subjects with two-level versus one-level stenosis (0.136 Pa s/ml vs. 0.069 Pa s/ml averages). Subjects with two-level stenosis had greater particle deposition at sites of stenosis compared to subjects with one-level stenosis (average deposition 2.31% vs. 0.96%). One-level stenosis subjects, as well as one two-level stenosis subject, had the greatest deposition using MDI with a spacer (0 m/s): 2.59% and 4.34%, respectively. The second two-level stenosis subject had the greatest deposition using DPI (3.45%). Maximum deposition across all stenotic subtypes except one-level tracheal stenosis was achieved with particle sizes of 6-10 μm.

CONCLUSIONS

Our results suggest that patients with two-level LTS may experience a more constricted laryngotracheal airflow profile compared to patients with one-level LTS, which may enhance overall stenotic drug deposition.

LEVEL OF EVIDENCE

NA Laryngoscope, 133:366-374, 2023.

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.

Computational investigation of nasal surface coverage from squeeze bottle and Neti Pot saline irrigation flow

BACKGROUND AND OBJECTIVE

Nasal saline irrigation is a common therapy for inflammatory nasal and paranasal disease or for managing post nasal and sinus surgery recovery. Two common irrigation devices include the netipot and squeeze bottles, where anecdotally, these devices alleviate congestion, facial pain, and pressure. However, a quantitative evaluation of these devices' performance and the fluid dynamics responsible for the irrigation distribution through the nose is lacking. This study tracked the liquid surface coverage and wall shear stresses during nasal saline irrigation produced from a Neti Pot and squeeze bottle.

METHODS

This study used transient computational fluid dynamics (CFD) simulations to investigate the saline irrigation flow field in a subject-specific sinonasal model. The computational nasal cavity model was constructed from a high-resolution computed tomography scan (CT). The irrigation procedure applied a head position tilted at 90° forward using an 80 ml squeeze bottle and 120 ml Neti Pot.

RESULTS

The results from a single sinonasal model demonstrated that the Neti Pot irrigation was more effective in delivering saline solution to the nasal cavity on the contralateral side of irrigation due to typically larger volumes but at the expense of reduced flow and shearing rates, as the flow entered under gravitational forces. The squeeze bottle irrigation provided greater surface coverage on the side of irrigation.

CONCLUSIONS

The results from the single patient model, demonstrated the Neti Pot increased surface coverage in the paranasal sinuses. Reducing the jet diameter may aid the direct targeting of a specific region at the side of irrigation by preventing the impingement of the jet to the nasal passage surface and redirection of the flow. Evaluating this performance across a wider cohort of patients can strengthen the findings.

The Impact of Adhesions on Nasal Airflow: A Quantitative Analysis Using Computational Fluid Dynamics

Background

Nasal adhesions (NAs) are a known complication of nasal airway surgery. Even minor NAs can lead to significant postoperative nasal airway obstruction (NAO). Division of such NAs often provides much greater relief than anticipated.

Objective

We examine the impact of NAs at various anatomical sites on nasal airflow and mucosal cooling using computational fluid dynamics (CFD) and multiple test subjects.

Methods

CT scans of healthy adult subjects were used to construct three-dimensional nasal airway computational models. A single virtual 2.5 mm diameter NA was placed at one of five sites commonly seen following NAO surgery within each nasal cavity bilaterally, resulting in 10 NA models and 1 NA-free control for each subject. CFD analysis was performed on each NA model and compared with the subject's NA-free control model.

Results

4 subjects were recruited to create 44 computational models. The NAs caused the airflow streamlines to separate, leading to a statistically significant increase in mucosal temperature immediately downstream to the NAs (wake region). Changes in the mucosal temperature in the wake region of the NAs were most prominent in anteriorly located NAs with a mean increase of 1.62 °C for the anterior inferior turbinate NAs ( P < .001) and 0.63 °C for the internal valve NAs ( P < .001).

Conclusion

NAs result in marked disruption to airflow patterns and reduced mucosal cooling on critical surfaces, particularly in the wake region. Reduced wake region mucosal cooling may be a contributing factor to the exaggerated perception of nasal obstruction experienced by patients with NAs.

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

OBJECTIVE

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

MATERIALS AND METHODS

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

RESULTS AND DISCUSSION AND CONCLUSION

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

Analysis of improved oral drug delivery with different helical stream inhalation modes

A challenging aspect of pulmonary drug delivery devices, e.g., metered dose inhalers (MDIs), is to deliver therapeutic drugs to prescribed target locations at the required dosage level. In this study, validated computer simulations of micron-drug inhalation with angled or radially positioned helical fluid-particle streams are simulated and analyzed. For a suitable swirl number significant improvements in drug delivery, especially to deeper lung regions, have been achieved. Specifically, considering realistic polydisperse particle distributions at the mouth inlet for a subject-specific upper lung airway geometry, a 10-degree angled helical stream increased the local efficacy by up to 26% in comparison to a conventional helical stream, causing an overall dosage of about 60% to the deep lung. Considering lobe-specific drug targeting scenarios, while using an off-center, i.e., radially well positioned, helical-flow mouthpiece, the local particle-deposition efficacy increased from 9% to 24% in the left lobe and from 25% to 38% in the right lobe in comparison to conventional drug-aerosol stream released from the central position. The efficacy of helical streams for pulmonary drug delivery applications has been established.

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

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

Evaluation of the Polyhedral Mesh Style for Predicting Aerosol Deposition in Representative Models of the Conducting Airways

A critical factor affecting the accuracy of Computational Fluid Dynamic (CFD) simulations and the time required to conduct them is construction of the computational mesh. This study aimed to evaluate the relatively new polyhedral mesh style for simulating aerosol deposition in the upper conducting airways compared with established meshing techniques and experimental data. Hexahedral and polyhedral mesh solutions were compared in two benchmark geometries: 1) a 90°-bend with flow characteristics similar to the extrathoracic airways of an adolescent child, and 2) a double bifurcation representing bifurcations B3-B5 in an adult. Both 4-block and 5-block hexahedral meshes were used in the 90°-bend to capture the potential of fully-structured hexahedral meshes. In the 90°-bend, polyhedral elements matched polydisperse in vitro deposition data with 20% relative error (RE; averaged across the particle sizes considered), which is an improvement on the accuracy of the 4-block hexahedral mesh (35% RE) and is similar to the accuracy of the 5-block hexahedral mesh (19% RE). In the double bifurcation, deposition fraction relative differences evaluated between polyhedral and hexahedral meshes ranged from 0.3% to 28.6% for the different particle sizes assessed, which is an order of magnitude improvement compared with previous studies that considered hexahedral vs. hybrid tetrahedral-prism meshes for the same flow field. Solution convergence time with polyhedral elements was found to be 50% to 140% higher than with hexahedral meshes of comparable size. While application dependent, the increase in simulation time observed with polyhedral meshes will likely be outweighed by the ease and convenience of polyhedral mesh construction. It was concluded that the polyhedral mesh style, with sufficient resolution especially near the walls, is an excellent alternative to the highly regarded hexahedral mesh style for predicting upper airway aerosol transport and deposition and provides a powerful new tool in the assessment of respiratory aerosol dosimetry.

Effect of swirling flow and particle-release pattern on drug delivery to human tracheobronchial airways

The present study aims to investigate the effect of swirling flow on particle deposition in a realistic human airway. A computational fluid dynamic (CFD) model was utilized for the simulation of oral inhalation and particle transport patterns, considering the k-ω turbulence model. Lagrangian particle tracking was used to track the particles' trajectories. A normal breathing condition (30 L/min) was applied, and two-micron particles were injected into the mouth, considering swirling flow to the oral inhalation airflow. Different cases were considered for releasing the particles, which evaluated the impacts of various parameters on the deposition efficiency (DE), including the swirl intensity, injection location and pattern of the particle. The work's novelty is applying several injection locations and diameters simultaneously. The results show that the swirling flow enhances the particle deposition efficiency (20-40%) versus no-swirl flow, especially in the mouth. However, releasing particles inside the mouth, or injecting them randomly with a smaller injection diameter (d_inj) reduced DE in swirling flow condition, about 50 to 80%. Injecting particles inside the mouth can decrease DE by about 20%, and releasing particles with smaller d_inj leads to 50% less DE in swirling flow. In conclusion, it is indicated that the airflow condition is an important parameter for a reliable drug delivery, and it is more beneficial to keep the inflow uniform and avoid swirling flow.

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

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

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

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

Effects of head tilt on squeeze-bottle nasal irrigation - A computational fluid dynamics study

Nasal irrigation is a widely recognized treatment for several sinonasal diseases. However, there is a lack of clear evidence-based guidelines for optimal irrigation delivery to improve lavage and topical drug delivery. This study uses computational fluid dynamics (CFD) to assess the effects of different head tilt positions on sinonasal coverage, residence time and shear stresses in squeeze-bottle nasal irrigation. A sinonasal cavity computational model was constructed from a high-resolution CT scan of a healthy, 25-year-old Asian female. The Volume of Fluid method was used to track the interface between the two immiscible fluids (air and water). The direction of gravity was varied to simulate different head tilt-positions (0° Straight, 45° Forward, 45° Left, 45° Right and 45° Backward) during nasal irrigation with 150 mL liquid via a squeeze bottle through the left nostril for 2 s with a 0.1 s acceleration/deceleration time. The results showed that the 45° backward head tilt position was the most effective in delivering irrigation to the ethmoid, frontal and sphenoid sinuses. Altering head tilt had minimal impact on irrigation delivery to the maxillary sinuses. Maximum wall shear stresses seen in localized areas of the sinus mucosa varied significantly with different head tilt angles. However, the difference in mean wall shear stress on the sinus surfaces was marginal with changing head tilt position. The findings suggest that an optimized head tilt position can be identified to improve liquid irrigation to targeted sinuses, as per treatment requirements (lavage and topical drug delivery).

Orally Inhaled Drug Particle Transport in Computerized Models of Laryngotracheal Stenosis

OBJECTIVE

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

STUDY DESIGN

Descriptive research.

SETTING

Academic medical center.

METHODS

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

RESULTS

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

CONCLUSION

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

Computational characterization of inhaled droplet transport to the nasopharynx

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

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

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

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.

Rhinomanometry Versus Computational Fluid Dynamics: Correlated, but Different Techniques

BACKGROUND

Past studies reported a low correlation between rhinomanometry and computational fluid dynamics (CFD), but the source of the discrepancy was unclear. Low correlation or lack of correlation has also been reported between subjective and objective measures of nasal patency.

OBJECTIVE

This study investigates (1) the correlation and agreement between nasal resistance derived from CFD (R_CFD) and rhinomanometry (R_RMN), and (2) the correlation between objective and subjective measures of nasal patency.

METHODS

Twenty-five patients with nasal obstruction underwent anterior rhinomanometry before and after mucosal decongestion with oxymetazoline. Subjective nasal patency was assessed with a 0-10 visual analog scale (VAS). CFD simulations were performed based on computed tomography scans obtained after mucosal decongestion. To validate the CFD methods, nasal resistance was measured in vitro (R_EXPERIMENT) by performing pressure-flow experiments in anatomically accurate plastic nasal replicas from 6 individuals.

RESULTS

Mucosal decongestion was associated with a reduction in bilateral nasal resistance (0.34 ± 0.23 Pa.s/ml to 0.19 ± 0.24 Pa.s/ml, p = 0.003) and improved sensation of nasal airflow (bilateral VAS decreased from 5.2 ± 1.9 to 2.6 ± 1.9, p < 0.001). A statistically significant correlation was found between VAS in the most obstructed cavity and unilateral airflow before and after mucosal decongestion (r = -0.42, p = 0.003). Excellent correlation was found between R_CFD and R_EXPERIMENT (r = 0.96, p < 0.001) with good agreement between the numerical and in vitro values (R_CFD/R_EXPERIMENT = 0.93 ± 0.08). A weak correlation was found between R_CFD and R_RMN (r = 0.41, p = 0.003) with CFD underpredicting nasal resistance derived from rhinomanometry (R_CFD/R_RMN = 0.65 ± 0.63). A stronger correlation was found when unilateral airflow at a pressure drop of 75 Pa was used to compare CFD with rhinomanometry (r = 0.76, p < 0.001).

CONCLUSION

CFD and rhinomanometry are moderately correlated, but CFD underpredicts nasal resistance measured in vivo due in part to the assumption of rigid nasal walls. Our results confirm previous reports that subjective nasal patency correlates better with unilateral than with bilateral measurements and in the context of an intervention.

Nasal Airflow Changes With Bioabsorbable Implant, Butterfly, and Spreader Grafts

OBJECTIVES/HYPOTHESIS

Internal nasal valve compromise is a major cause of nasal obstruction, with a growing number of ways to treat this condition. In this study, we compared the effects of butterfly graft, spreader graft, and the bioabsorbable nasal implant on nasal airflow resistance.

STUDY DESIGN

Cadaver study.

METHODS

Computational fluid dynamics (CFD) simulations were completed from nine preoperative and postoperative cadaveric subjects. Each cadaveric head underwent placement of a bioabsorbable nasal implant (BNI) (Spirox Latera; Stryker ENT, Plymouth, MN), butterfly graft, or spreader graft. Pre- and postoperative computed tomography (CT) scans were used to generate three-dimensional models of the nasal airway used in steady-state CFD simulations of airflow and heat transfer during inspiration.

RESULTS

Butterfly graft placement resulted in a mean improvement in nasal airway resistance of 24.9% (±7.3), whereas BNI placement resulted in a 6.7% (±1.2) improvement, and spreader graft placement also resulted in a consistent improvement of 2.6% (±13.5). Pressure within the main nasal cavity was consistently lower following butterfly graft placement versus a spreader graft or BNI. Butterfly and spreader graft placement also resulted in modest improvements in airflow allocation, whereas BNI demonstrated more variation (-1% to 12%). Heat flux was not significantly different; however, a small improvement in total heat flux was seen with all three interventions.

CONCLUSIONS

The results of this study demonstrate reduction in nasal airway resistance in all three surgical interventions, with the butterfly graft demonstrating superiority to the other two techniques. However, these data only reflect a static environment and not dynamic changes in airflow seen during respiration.

LEVEL OF EVIDENCE

NA Laryngoscope, 130:E817-E823, 2020.

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.

Comparison of Airflow Between Spreader Grafts and Butterfly Grafts Using Computational Flow Dynamics in a Cadaveric Model

Importance

Nasal valve compromise is a major cause of nasal obstruction, and multiple methods have been developed to address it.

Objective

To compare nasal airflow resistance, airflow partitioning, and mucosal cooling (heat flux) before and after 2 surgical interventions, butterfly and spreader graft placement, used to treat nasal valve compromise.

Design, Setting, and Participants

In this cadaveric tissue study, 4 fresh cadaveric heads underwent both spreader graft and butterfly graft surgical procedures in alternating sequence in March 2016. Preoperative and postoperative computed tomographic scans were used to generate 3-dimensional (3-D) models of the nasal airway. These models were then used in steady state computational fluid dynamics simulations of airflow and heat transfer during inspiration.

Intervention

Butterfly and spreader graft techniques.

Main Outcomes and Measures

Nasal airflow resistance, airflow partitioning, and heat flux.

Results

Donors 1, 2, and 3 were white males; donor 4, a white female. Computational fluid dynamics simulations during inspiration in 3-D models generated from preoperative and postoperative computed tomographic scans of the 4 cadaveric heads indicated reductions from preoperative values in nasal airflow resistance associated with both butterfly grafts (range, 20%-51%) and spreader grafts (range, 2%-29%). Butterfly grafts were associated with a greater reduction in nasal airflow resistance in models of all 4 cadaveric heads. Changes from preoperative values for heat flux, a biophysical variable that correlates with the subjective sensation of nasal patency, were more variable, ranging from -11% to 4% following butterfly grafts and -9% to 10% following spreader grafts. The preoperative airflow allocation in the left and right nostrils improved consistently with the butterfly graft. With the spreader graft, there were improvements for donors 1 and 4, but the allocations were worse for donors 2 and 3.

Conclusions and Relevance

The results of this study suggest that the more recently developed butterfly graft technique may be associated with a similar level of improved nasal airflow as that observed with the use of a spreader graft in nasal valve compromise. Both interventions were associated with comparable changes in heat flux. Because this study addressed only static internal nasal valve stenosis, even greater differences in air flow and heat flux between the 2 techniques may be anticipated in a dynamic model. Further investigation in patients is warranted.

Level of Evidence

NA.

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

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

Computational Analysis of the Mature Unilateral Cleft Lip Nasal Deformity on Nasal Patency

Background

Nasal airway obstruction (NAO) due to nasal anatomic deformities is known to be more common among cleft patients than the general population, yet information is lacking regarding severity and variability of cleft-associated nasal obstruction relative to other conditions causing NAO. This preliminary study compares differences in NAO experienced by unilateral cleft lip nasal deformity (uCLND) subjects with noncleft subjects experiencing NAO.

Methods

Computational modeling techniques based on patient-specific computed tomography images were used to quantify the nasal airway anatomy and airflow dynamics in 21 subjects: 5 healthy normal subjects; 8 noncleft NAO subjects; and 8 uCLND subjects. Outcomes reported include Nasal Obstruction Symptom Evaluation (NOSE) scores, cross-sectional area, and nasal resistance.

Results

uCLND subjects had significantly larger cross-sectional area differences between the left and right nasal cavities at multiple cross sections compared with normal and NAO subjects. Median and interquartile range (IQR) NOSE scores between NAO and uCLND were 75 (IQR = 22.5) and 67.5 (IQR = 30), respectively. Airflow partition difference between both cavities were: median = 9.4%, IQR = 10.9% (normal); median = 31.9%, IQR = 25.0% (NAO); and median = 29.9%, IQR = 44.1% (uCLND). Median nasal resistance difference between left and right nasal cavities were 0.01 pa.s/ml (IQR = 0.03 pa.s/ml) for normal, 0.09 pa.s/ml (IQR = 0.16 pa.s/ml) for NAO and 0.08 pa.s/ml (IQR = 0.25 pa.s/ml) for uCLND subjects.

Conclusions

uCLND subjects demonstrated significant asymmetry between both sides of the nasal cavity. Furthermore, there exists substantial disproportionality in flow partition difference and resistance difference between cleft and noncleft sides among uCLND subjects, suggesting that both sides may be dysfunctional.

Effect of scaffold architecture on cell seeding efficiency: A discrete phase model CFD analysis

Within perfusion cell culture systems, scaffold architecture is able to control important biological parameters such as permeability and fluid flow-induced shear stress. As well, one of the main factors affecting the final fate of this process as well as optimal cell differentiation and proliferation in these systems is initial adhesion of cells to scaffolds. In this study, the effect of scaffold architecture on the adhesion of the cells was computationally investigated. For this purpose, four scaffold models including double-diamond, gyroid, FR-D, and Schwarz-primitive were designed using triply periodic minimal surface (TPMS) geometry with a constant porosity of 80%. As well, the inlet velocity of zero to simulate static cell culture and three different inlet velocities for modeling the dynamic cell culture conditions were also selected. The results showed that cell culture efficiency of scaffolds could be changed up to seven times from architecture to architecture under the same conditions. The efficiency of cell culture in scaffolds with tortuous architecture was also reported higher than those with relatively straight microchannels. In terms of culture methods, unlike dynamic cell culture model in which almost a homogeneous cell distribution was observed in static cell culture simulation, more cells adhered, but they had agglomerated in the scaffold entrance regions and had failed to reach all regions. The results of this study shed more light on the selection and design of scaffold architecture for optimal cell culture in tissue engineering.

Impact of endoscopic craniofacial resection on simulated nasal airflow and heat transport

BACKGROUND

Endoscopic craniofacial resections (CFR) are performed for extensive anterior skull base lesions. This surgery involves removal of multiple intranasal structures, potentially leading to empty nose syndrome (ENS). However, many patients remain asymptomatic postoperatively. Our objective was to analyze the impact of CFR on nasal physiology and airflow using computational fluid dynamics (CFD). This is the first CFD analysis of post-CFR patients.

METHODS

Three-dimensional sinonasal models were constructed from 3 postoperative images using Mimics^TM . Hybrid computational meshes were created. Steady inspiratory airflow and heat transport were simulated at patient-specific flow rates using shear stress transport k-omega turbulent flow modeling in Fluent^TM . Simulated average heat flux (HF) and surface area where HF exceeded 50 W/m² (SAHF50) were compared with laminar simulations in 9 radiographically normal adults.

RESULTS

Three adults underwent CFR without developing ENS. Average HF (W/m² ) were 132.70, 134.84, and 142.60 in the CFR group, ranging from 156.24 to 234.95 in the nonoperative cohort. SAHF50 (m² ) values were 0.0087, 0.0120, and 0.0110 in the CFR group, ranging from 0.0082 to 0.0114 in the radiographically normal cohort. SAHF50 was distributed throughout the CFR cavities, with increased HF at the roof and walls compared with the nonoperative cohort.

CONCLUSION

Average HF was low in the CFR group compared with the nonoperative group. However, absence of ENS in most CFR patients may be due to large stimulated mucosal surface area, commensurate with the nonoperative cohort. Diffuse distribution of stimulated area may result from turbulent mixing after CFR. To better understand heat transport post-CFR, a larger cohort is necessary.

Comparative study of simulated nebulized and spray particle deposition in chronic rhinosinusitis patients

BACKGROUND

Topical intranasal drugs are widely prescribed for chronic rhinosinusitis (CRS), although delivery can vary with device type and droplet size. The study objective was to compare nebulized and sprayed droplet deposition in the paranasal sinuses and ostiomeatal complex (OMC) across multiple droplet sizes in CRS patients using computational fluid dynamics (CFD).

METHODS

Three-dimensional models of sinonasal cavities were constructed from computed tomography (CT) scans of 3 subjects with CRS refractory to medical therapy using imaging software. Assuming steady-state inspiratory airflow at resting rate, CFD was used to simulate 1-µm to 120-µm sprayed droplet deposition in the left and right sinuses and OMC with spray nozzle positioning as in current nasal spray use instructions. Zero-velocity nebulization simulations were performed for 1-µm to 30-µm droplet sizes, maximal sinus and OMC deposition fractions (MSDF) were obtained, and sizes that achieved at least 50% of MSDF were identified. Nebulized MSDF was compared to sprayed droplet deposition. We also validated CFD framework through in vitro experiments.

RESULTS

Among nebulized droplet sizes, 11-µm to 14-µm droplets achieved at least 50% of MSDF in all 6 sinonasal cavities. Four of 6 sinonasal cavities had greater sinus and OMC deposition with nebulized droplets than with sprayed droplets at optimal sizes.

CONCLUSION

Nebulized droplets may target the sinuses and OMC more effectively than sprayed particles at sizes achieving best deposition. Further studies are needed to confirm our preliminary findings. Several commercial nasal nebulizers have average particle sizes outside the optimal nebulized droplet size range found here, suggesting potential for product enhancement.

A hierarchical stepwise approach to evaluate nasal patency after virtual surgery for nasal airway obstruction

BACKGROUND

Despite advances in medicine and expenditures associated in treatment of nasal airway obstruction, 25-50% of patients undergoing nasal surgeries complain of persistent obstructive symptoms. Our objective is to develop a "stepwise virtual surgery" method that optimizes surgical outcomes for treatment of nasal airway obstruction.

METHODS

Pre-surgery radiographic images of two subjects with nasal airway obstruction were imported into Mimics imaging software package for three-dimension reconstruction of the airway. A hierarchical stepwise approach was used to create seven virtual surgery nasal models comprising individual (inferior turbinectomy or septoplasty) procedures and combined inferior turbinectomy and septoplasty procedures via digital modifications of each subject's pre-surgery nasal model. To evaluate the effects of these procedures on nasal patency, computational fluid dynamics modeling was used to perform steady-state laminar inspiratory airflow and heat transfer simulations in every model, at resting breathing. Airflow-related variables were calculated for virtual surgery models and compared with dataset containing results of healthy subjects with no symptoms of nasal obstruction.

FINDINGS

For Subject 1, nasal models with virtual septoplasty only and virtual septoplasty plus inferior turbinectomy on less obstructed side were within the healthy reference thresholds on both sides of the nasal cavity and across all three computed variables. For Subject 2, virtual septoplasty plus inferior turbinectomy on less obstructed side model produced the best result.

INTERPRETATION

The hierarchical stepwise approach implemented in this preliminary report demonstrates computational fluid dynamics modeling ability to evaluate the efficiency of different surgical procedures for nasal obstruction in restoring nasal patency to normative level.

High-Efficiency Nose-to-Lung Aerosol Delivery in an Infant: Development of a Validated Computational Fluid Dynamics Method

Background: Computational fluid dynamics (CFD) provides a powerful tool for developing new high-efficiency aerosol delivery strategies, such as nose-to-lung (N2L) aerosol administration to infants and children using correctly sized aerosols. The objective of this study was to establish numerically efficient CFD solution methods and guidelines for simulating N2L aerosol administration to an infant based on comparisons with concurrent in vitro experiments. Materials and Methods: N2L administration of a micrometer-sized aerosol (mass median aerodynamic diameter [MMAD] = 1.4 μm) was evaluated using concurrent CFD simulations and in vitro experiments. Aerosol transport and deposition was assessed in a new nasal airway geometry of a 6-month-old infant with a streamlined nasal cannula interface, which was constructed as a CFD mesh and three-dimensionally printed to form an identical physical prototype. CFD meshes explored were a conventional tetrahedral approach with near-wall (NW) prism elements and a new polyhedral mesh style with an equally refined NW layer. The presence of turbulence in the model was evaluated using a highly efficient low-Reynolds number (LRN) k-ω turbulence model, with previously established NW corrections that accounted for anisotropic wall-normal turbulence as well as improved NW velocity interpolations and hydrodynamic particle damping. Results: Use of the new polyhedral mesh was found to improve numerical efficiency by providing more rapid convergence and requiring fewer control volumes. Turbulent flow was found in the nasal geometry, generated by the inlet jets from the nasal cannula interface. However, due to the small particle size, turbulent dispersion was shown to have little effect on deposition. Good agreement was established between the CFD predictions using the numerically efficient LRN k-ω model with appropriate NW corrections and in vitro deposition data. Aerosol transmission efficiencies through the delivery tube, nasal cannula, and infant nasal model, based on experimental and CFD predictions, were 93.0% and 91.5%, respectively. Conclusions: A numerically efficient CFD approach was established to develop transnasal aerosol administration to infants and children. Small particle aerosols with aerodynamic diameters of ∼1.5 μm were confirmed to have low inertial depositional loss, and have low deposition from turbulent dispersion, making them ideal for high-efficiency lung delivery through an infant nasal cannula interface.

Examining mesh independence for flow dynamics in the human nasal cavity

Increased computational resources provide new opportunities to explore sophisticated respiratory modelling. A survey of recent publications showed a steady increase in the number of mesh elements used in computational models over time. Complex geometries such as the nasal cavity exhibit sharp gradients and irregular curvatures, leading to abnormal flow development across their surfaces. As such, a robust method for examining the near-wall mesh resolution is required. The non-dimensional wall unit y⁺ (often used in turbulent flows) was used as a parameter to evaluate the near-wall mesh in laminar flows. Mesh independence analysis from line profiles showed that the line location had a significant influence on the result. Furthermore, using a single line profile as a measure for mesh convergence was unsuitable for representing the entire flow field. To improve this, a two-dimensional (2D) cross-sectional plane subtraction method where scalar values (such as the velocity magnitude) on a cross-sectional plane were interpolated onto a regularly spaced grid was proposed. The new interpolated grid values from any two meshed models could then be compared for changes caused by the different meshed models. The application of this method to three-dimensional (3D) volume subtraction was also demonstrated. The results showed that if the near-wall mesh was sufficiently refined, then narrow passages were less reliant on the overall mesh size. However, in wider passages, velocity magnitudes were sensitive to mesh size, requiring a more refined mesh.

Upper airway reconstruction using long-range optical coherence tomography: Effects of airway curvature on airflow resistance

OBJECTIVES

Adenotonsillectomy (AT) is commonly used to treat upper airway obstruction in children, but selection of patients who will benefit most from AT is challenging. The need for diagnostic evaluation tools without sedation, radiation, or high costs has motivated the development of long-range optical coherence tomography (LR-OCT), providing real-time cross-sectional airway imaging during endoscopy. Since the endoscope channel location is not tracked in conventional LR-OCT, airway curvature must be estimated and may affect predicted airway resistance. The study objective was to assess effects of three realistic airway curvatures on predicted airway resistance using computational fluid dynamics (CFD) in LR-OCT reconstructions of the upper airways of pediatric patients, before and after AT.

METHODS

Eight subjects (five males, three females, aged 4-9 years) were imaged using LR-OCT before and after AT during sedated endoscopy. Three-dimensional (3D) airway reconstructions included three airway curvatures. Steady-state, inspiratory airflow simulations were conducted under laminar conditions, along with turbulent simulations for one subject using the k-ω turbulence model. Airway resistance (pressure drop/flow) was compared using two-tailed Wilcoxon signed rank tests.

RESULTS

Regardless of the airway curvatures, CFD findings corroborate a surgical end-goal with computed post-operative airway resistance significantly less than pre-operative (P < 0.01). The individual resistances did not vary significantly for different airway curvatures (P > 0.25). Resistances computed using turbulent simulations differed from laminar results by less than ∼5%.

CONCLUSIONS

The results suggest that reconstruction of the upper airways from LR-OCT imaging data may not need to account for airway curvature to be predictive of surgical effects on airway resistance. Lasers Surg. Med. 51:150-160, 2019. © 2018 Wiley Periodicals, Inc.

Numerical simulation of thermal water delivery in the human nasal cavity

This work describes an extensive numerical investigation of thermal water delivery for the treatment of inflammatory disorders in the human nasal cavity. The numerical simulation of the multiphase air-droplets flow is based upon the Large Eddy Simulation (LES) technique, with droplets of thermal water described via a Lagrangian approach. Droplet deposition is studied for different sizes of water droplets, corresponding to two different thermal treatments, i.e. aerosol and inhalation. Numerical simulations are conducted on a patient-specific anatomy, employing two different grid sizes, under steady inspiration at two breathing intensities. The results are compared with published in vivo and in vitro data. The effectiveness of the various thermal treatments is then assessed qualitatively and quantitatively, by a detailed analysis of the deposition patterns of the droplets. Discretization effects on the deposition dynamics are addressed. The level of detail of the present work, together with the accuracy afforded by the LES approach, leads to an improved understanding of how the mixture of air-water droplets is distributed within the nose and the paranasal sinuses.

On computational fluid dynamics models for sinonasal drug transport: Relevance of nozzle subtraction and nasal vestibular dilation

Generating anatomically realistic 3-dimensional (3D) models of the human sinonasal cavity for numerical investigations of sprayed drug transport presents a host of methodological ambiguities. For example, subject-specific radiographic images used for 3D reconstructions typically exclude spray bottles. Subtracting a bottle contour from the 3D airspace and dilating the anterior nasal vestibule for nozzle placement augment the complexity of model building. So we explored the question: how essential are these steps to adequately simulate nasal airflow and identify the optimal delivery conditions for intranasal sprays? In particular, we focused on particle deposition patterns in the maxillary sinus, a critical target site for chronic rhinosinusitis. The models were reconstructed from postsurgery computed tomography scans for a 39-year-old Caucasian male, with chronic rhinosinusitis history. Inspiratory airflow patterns during resting breathing are reliably tracked through computational fluid dynamics-based steady-state laminar-viscous modeling, and such regimes portray relative lack of sensitivity to inlet perturbations. Consequently, we hypothesized that the posterior airflow transport and the particle deposition trends should not be radically affected by the nozzle subtraction and vestibular dilation. The study involved 1 base model and 2 derived models; the latter 2 with nozzle contours (2 different orientations) subtracted from the dilated anterior segment of the left vestibule. We analyzed spray transport in the left maxillary sinus for multiple release conditions. Similar release points, localized on an approximately 2 mm × 4.5 mm contour, facilitated improved maxillary deposition in all 3 test cases. This suggests functional redundancy of nozzle insertion in a 3D numerical model for identifying the optimal spray release locations.

Investigating the effects of laryngotracheal stenosis on upper airway aerodynamics

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

LEVEL OF EVIDENCE

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

Flow perfusion rate modulates cell deposition onto scaffold substrate during cell seeding

The combination of perfusion bioreactors with porous scaffolds is beneficial for the transport of cells during cell seeding. Nonetheless, the fact that cells penetrate into the scaffold pores does not necessarily imply the interception of cells with scaffold substrate and cell attachment. An in vitro perfusion system was built to relate the selected flow rate with seeding efficiency. However, the in vitro model does not elucidate how the flow rate affects the transport and deposition of cells onto the scaffold. Thus, a computational model was developed mimicking in vitro conditions to identify the mechanisms that bring cells to the scaffold from suspension flow. Static and dynamic cell seeding configurations were investigated. In static seeding, cells sediment due to gravity until they encounter the first obstacle. In dynamic seeding, 12, 120 and 600 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu \hbox {l/min}$$\end{document}μl/min flow rates were explored under the presence or the absence of gravity. Gravity and secondary flow were found to be key factors for cell deposition. In vitro and in silico seeding efficiencies are in the same order of magnitude and follow the same trend with the effect of fluid flow; static seeding results in higher efficiency than dynamic perfusion although irregular spatial distribution of cells was found. In dynamic seeding, 120 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu \hbox {l/min}$$\end{document}μl/min provided the best seeding results. Nevertheless, the perfusion approach reports low efficiencies for the scaffold used in this study which leads to cell waste and low density of cells inside the scaffold. This study suggests gravity and secondary flow as the driving mechanisms for cell-scaffold deposition. In addition, the present in silico model can help to optimize hydrodynamic-based seeding strategies prior to experiments and enhance cell seeding efficiency.

Ideal Particle Sizes for Inhaled Steroids Targeting Vocal Granulomas: Preliminary Study Using Computational Fluid Dynamics

Objectives Vocal fold granulomas are benign lesions of the larynx commonly caused by gastroesophageal reflux, intubation, and phonotrauma. Current medical therapy includes inhaled corticosteroids to target inflammation that leads to granuloma formation. Particle sizes of commonly prescribed inhalers range over 1 to 4 µm. The study objective was to use computational fluid dynamics to investigate deposition patterns over a range of particle sizes of inhaled corticosteroids targeting the larynx and vocal fold granulomas. Study Design Retrospective, case-specific computational study. Setting Tertiary academic center. Subjects/Methods A 3-dimensional anatomically realistic computational model of a normal adult airway from mouth to trachea was constructed from 3 computed tomography scans. Virtual granulomas of varying sizes and positions along the vocal fold were incorporated into the base model. Assuming steady-state, inspiratory, turbulent airflow at 30 L/min, computational fluid dynamics was used to simulate respiratory transport and deposition of inhaled corticosteroid particles ranging over 1 to 20 µm. Results Laryngeal deposition in the base model peaked for particle sizes 8 to 10 µm (2.8%-3.5%). Ideal sizes ranged over 6 to 10, 7 to 13, and 7 to 14 µm for small, medium, and large granuloma sizes, respectively. Glottic deposition was maximal at 10.8% for 9-µm-sized particles for the large posterior granuloma, 3 times the normal model (3.5%). Conclusion As the virtual granuloma size increased and the location became more posterior, glottic deposition and ideal particle size generally increased. This preliminary study suggests that inhalers with larger particle sizes, such as fluticasone propionate dry-powder inhaler, may improve laryngeal drug deposition. Most commercially available inhalers have smaller particles than suggested here.

Estimates of nasal airflow at the nasal cycle mid-point improve the correlation between objective and subjective measures of nasal patency

INTRODUCTION

The nasal cycle represents a significant challenge when comparing pre- and post-surgery objective measures of nasal airflow.

METHODS

Computational fluid dynamics (CFD) simulations of nasal airflow were conducted in 12 nasal airway obstruction patients showing significant nasal cycling between pre- and post-surgery computed tomography scans. To correct for the nasal cycle, mid-cycle models were created virtually. Subjective scores of nasal patency were obtained via the Nasal Obstruction Symptom Evaluation (NOSE) and unilateral visual analog scale (VAS).

RESULTS

The correlation between objective and subjective measures of nasal patency increased after correcting for the nasal cycle. In contrast to biophysical variables in individual patients, cohort averages were not significantly affected by the nasal cycle correction.

CONCLUSIONS

The ability to correct for the confounding effect of the nasal cycle is a key element that future virtual surgery planning software for nasal airway obstruction will need to account for when using anatomic models based on single instantaneous imaging.

Understanding the mechanisms underlying pulsating aerosol delivery to the maxillary sinus: In vitro tests and computational simulations

BACKGROUND

Pulsating aerosol delivery has been demonstrated in depositing medications into paranasal sinuses. However, its mechanisms are not fully understood. Influences of the nasal anatomy and sound frequency on intrasinus delivery are not yet clear.

OBJECTIVES

This study aimed to gain a better understanding of the mechanisms for enhanced intrasinus delivery with pulsating sound. Specifically, effects of the pulsation frequency, ostium size, and sinus shape on the intrasinus dosage and resonance frequency would be examined.

METHODS AND MATERIALS

Both experiments and computational modeling were conducted to understand the pulsating aerosol delivery in both idealized (two-bottle) and realistic nose-sinus models. A computational model of intrasinus pulsation delivery was developed using COMSOL and was cross-validated with both experimental and theoretical results.

RESULTS

In contrast to previous studies, seemingly erratic relations between the intrasinus dosage and ostium diameter were observed in experiments, which suggested a more complicated particle transport mechanism. Improved agreement was achieved when grouping the ostium size and sinus volume into the resonance frequency, and therefore, validated the hypothesis that intrasinus deposition strongly depends on the resonance frequency. Extensive computational simulations revealed that the deposition was highest at the resonance frequency and decreased gradually at off-resonance frequencies. The resonance frequency depended on the ostium and sinus morphology, but was independent of the nasal cavity.

CONCLUSION

Results of this study verified the hypothesis of resonance being the mechanism for enhanced particle deposition in the maxillary sinus. A better knowledge of the relationship between sinus dosages, pulsating frequency, and nasal morphometry is essential for improving the design of intrasinus delivery devices.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Characterizing human nasal airflow physiologic variables by nasal index

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