Scale resolving simulations of the effect of glottis motion and the laryngeal jet on flow dynamics during respiration

BACKGROUND AND OBJECTIVE

The movement of the respiratory walls has a significant impact on airflow through the respiratory tract. The majority of computational fluid dynamics (CFD) studies assume a static geometry which may not provide a realistic flow field. Furthermore, many studies use Reynolds Averaged Navier-Stokes (RANS) turbulence models that do not resolve turbulence structure. Combining the application of advanced scale-resolving turbulence models with moving respiratory walls using CFD will provide detailed insights into respiratory flow structures.

METHODS

This study simulated a complete breathing cycle involving inhalation and exhalation in a nasal cavity to trachea geometry that incorporated moving glottis walls. A second breathing cycle was simulated with static glottis walls for comparison. A recently developed hybrid RANS-LES turbulence model, the Stress-Blended Eddy Simulation (SBES), was incorporated to resolve turbulent flow structures in fine detail for both transient simulations. Transient results were compared with steady-state RANS simulations for the same respiratory geometry.

RESULTS

Glottis motion caused substantial effects on flow structure through the complete breathing cycle. Significant flow structure and velocity variations were observed due to glottal motion, primarily in the larynx and trachea. Resolved turbulence structures using SBES showed an intense mixing section in the glottis region during inhalation and in the nasopharynx during expiration, which was not present in the RANS simulations.

CONCLUSION

Transient simulations of a realistic breathing cycle uncovered flow structures absent in simulations with a constant flow rate. Furthermore, the incorporation of glottis motion impacted airflow characteristics that suggest rigid respiratory walls do not accurately describe respiratory flow. Future research in respiratory airflow should be conducted using transient scale-resolving models in conjunction with moving respiratory walls to capture flow structures in detail.

Modelling the effects of post-FESS middle turbinate synechiae on sinonasal physiology: A computational fluid dynamics study

OBJECTIVE(S)

Chronic rhinosinusitis (CRS) is common and often requires surgical intervention. Surgical failure may lead to persistent symptoms and recalcitrant disease, often secondary to synechiae between the middle turbinate (MT) and lateral nasal wall. Synechiae prevention techniques have been extensively investigated, however evidence for the effect of synechiae on sinonasal physiology is lacking. We aimed to model the effects of MT synechiae on a post-functional endoscopic sinus surgery (FESS) sinonasal cavity using computational fluid dynamics (CFD).

METHODS

DICOM data from a CT-sinus of a healthy 25-year-old female was segmented to create a three-dimensional model. Virtual surgery was performed to simulate a "full-house" FESS procedure. Multiple models were created, each with a single unilateral virtual MT synechia of varying extent. CFD analysis was performed on each model and compared with a post-FESS control model without synechiae. Airflow velocity, humidity and mucosal surface and air temperature values were calculated.

RESULTS

All synechiae models demonstrated aberrant downstream sinonasal airflow. There was reduced ventilation of the ipsilateral frontal, ethmoid and sphenoid sinuses, with a concentrated central "jet" in the middle meatus region. Effects were proportionate to the size of synechiae. The impact on bulk inspired airflow was negligible.

CONCLUSION

Post-FESS synechiae between the MT and lateral nasal wall significantly disrupt local downstream sinus ventilation and nasal airflow. These findings may explain the persistent symptoms seen in post-FESS CRS patients with MT synechiae, reinforcing the importance of prevention and adhesiolysis. Larger cohort studies with multiple models of actual post-FESS patients with synechiae are required to validate these findings.

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

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

In-silico decongested trial effects on the impaired breathing function of a bulldog suffering from severe brachycephalic obstructive airway syndrome

BACKGROUND AND OBJECTIVE

Brachycephalic obstructive airway syndrome (BOAS) susceptible dogs (e.g., French bulldog), suffer health complications related to deficient breathing primarily due to anatomical airway geometry. Surgical interventions are known to provide acceptable functional and cosmetic results; however, the long-term post-surgery outcome is not well known. In silico analysis provides an objective measure to quantify the respiratory function in postoperative dogs which is critical for successful long-term outcomes. A virtual surgery to open the airway can explore the ability for improved breathing in an obstructed airway of a patient dog, thus supporting surgeons in pre-surgery planning using computational fluid dynamics.

METHODS

In this study five surgical interventions were generated with a gradual increment of decongested levels in a bulldog based on computed tomography images. The effects of the decongested airways on the breathing function of a patient bulldog, i.e., airflow characteristics, pressure drop, wall shear stress, and air-conditioning capacity, were quantified by benchmarking against a clinically healthy bulldog using computational fluid dynamics (CFD) method.

RESULTS

Our findings demonstrated a promising decrease in excessive airstream velocity, pressure drop, and wall shear stress in virtual surgical scenarios, while constantly preserving adequate air-conditioning efficiency. A linear fit curve was proposed to correlate the reduction in the pressure drop and decongested level.

CONCLUSIONS

The in silico analysis is a viable tool providing visual and quantitative insight into new unexplored surgical techniques.

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

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

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.

Flow Patterns and Particle Residence Times in the Oral Cavity during Inhaled Drug Delivery

Pulmonary drug delivery aims to deliver particles deep into the lungs, bypassing the mouth−throat airway geometry. However, micron particles under high flow rates are susceptible to inertial impaction on anatomical sites that serve as a defense system to filter and prevent foreign particles from entering the lungs. The aim of this study was to understand particle aerodynamics and its possible deposition in the mouth−throat airway that inhibits pulmonary drug delivery. In this study, we present an analysis of the aerodynamics of inhaled particles inside a patient-specific mouth−throat model generated from MRI scans. Computational Fluid Dynamics with a Discrete Phase Model for tracking particles was used to characterize the airflow patterns for a constant inhalation flow rate of 30 L/min. Monodisperse particles with diameters of 7 μm to 26 μm were introduced to the domain within a 3 cm-diameter sphere in front of the oral cavity. The main outcomes of this study showed that the time taken for particle deposition to occur was 0.5 s; a narrow stream of particles (medially and superiorly) were transported by the flow field; larger particles > 20 μm deposited onto the oropharnyx, while smaller particles < 12 μm were more disperse throughout the oral cavity and navigated the curved geometry and laryngeal jet to escape through the tracheal outlet. It was concluded that at a flow rate of 30 L/min the particle diameters depositing on the larynx and trachea in this specific patient model are likely to be in the range of 7 μm to 16 μm. Particles larger than 16 μm primarily deposited on the oropharynx.

Nasal Irrigation Delivery in Three Post-FESS Models From a Squeeze-bottle Using CFD

Purpose

Nasal saline irrigation is highly recommended in patients following functional endoscopic sinus surgery (FESS) to aid the postoperative recovery. Post-FESS patients have significantly altered anatomy leading to markedly different flow dynamics from those found in pre-op or non-diseased airways, resulting in unknown flow dynamics.

Methods

This work investigated how the liquid stream disperses through altered nasal cavities following surgery using Computational Fluid Dynamics (CFD). A realistic squeeze profile was determined from physical experiments with a 27-year-old male using a squeeze bottle with load sensors. The administration technique involved a head tilt of 45-degrees forward to represent a head position over a sink. After the irrigation event that lasted 4.5 s, the simulation continued for an additional 1.5 s, with the head orientation returning to an upright position.

Results

The results demonstrated that a large maxillary sinus ostium on the right side allows saline penetration into this sinus. The increased volume of saline entering the maxillary sinus limits the saline volume available to the rest of the sinonasal cavity and reduces the surface coverage of the other paranasal sinuses. The average wall shear stress was higher on the right side than on the other side for two patients. The results also revealed that head position alters the sinuses’ saline residual, especially the frontal sinuses.

Conclusion

While greater access to sinuses is achieved through FESS surgery, patients without a nasal septum limits posterior sinus penetration due to the liquid crossing over to the contralateral cavity and exiting the nasal cavity early.

Correlation of Nasal Mucosal Temperature and Nasal Patency—A Computational Fluid Dynamics Study

OBJECTIVES

Recent evidence suggests that detection of nasal mucosal temperature, rather than direct airflow detection, is the primary determinant of subjective nasal patency. This study examines the role of nasal mucosal temperature in the perception of nasal patency using in vivo and computational fluid dynamics (CFD) measurements.

METHODS

Healthy adult participants completed Nasal Obstruction Symptom Evaluation (NOSE) and Visual Analogue Scale (VAS) questionnaires. A temperature probe measured nasal mucosal temperature at the vestibule, inferior turbinate, middle turbinate, and nasopharynx bilaterally. Participants underwent a CT scan, used to create a 3D nasal anatomy model to perform CFD analysis of nasal mucosal and inspired air temperature and heat flux along with mucosal surface area where heat flux >50 W/m² (SAHF50).

RESULTS

Eleven participants with a median age of 27 (IQR 24; 48) were recruited. Probe-measured temperature values correlated strongly with CFD-derived values (r = 0.87, p < 0.05). Correlations were seen anteriorly in the vestibule and inferior turbinate regions between nasal mucosal temperature and unilateral VAS (r = 0.42-0.46; p < 0.05), between SAHF50 and unilateral VAS (r = -0.31 to -0.36; p < 0.05) and between nasal mucosal temperature and SAHF50 (r = -0.37 to -0.41; p < 0.05). Subjects with high patency (VAS ≤10) had increased heat flux anteriorly compared with lower patency subjects (VAS >10; p < 0.05).

CONCLUSION

Lower nasal mucosal temperature and higher heat flux within the anterior nasal cavity correlates with a perception of improved unilateral nasal patency in healthy individuals.

LEVEL OF EVIDENCE

4 Laryngoscope, 133:1328-1335, 2023.

Wet surface wall model for latent heat exchange during evaporation

Air conditioning is a dual heat and mass transfer process, and the human nasal cavity achieves this through the mucosal wall surface, which is supplied with an energy source through the sub-epithelial network of capillaries. Computational studies of air conditioning in the nasal cavity have included temperature and humidity, but most studies solved these flow parameters separately, and in some cases, a constant mucosal surface temperature was used. Recent developments demonstrated that both heat and mass transfer need to be modeled. This work expands on existing modeling efforts in accounting for the nasal cavity's dual heat and mass transfer process by introducing a new subwall model, given in the Supplementary Materials. The model was applied to a pipe geometry, and a human nasal cavity was recreated from CT-scans, and six inhalation conditions were studied. The results showed that when the energy transfer from the latent heat of evaporation is included, there is a cooling effect on the mucosal surface temperature.

Combining experimental and computational techniques to understand and improve dry powder inhalers

INTRODUCTION

: Dry Powder Inhalers (DPIs) continue to be developed to deliver an expanding range of drugs to treat an ever-increasing range of medical conditions; with each drug and device combination needing a specifically designed inhaler. Fast regulatory approval is essential to be first to market, ensuring commercial profitability.

AREAS COVERED

: In vitro deposition, particle image velocimetry, and computational modelling using the physiological geometry and representative anatomy can be combined to give complementary information to determine the suitability of a proposed inhaler design and to optimise its formulation performance. In combination they allow the entire range of questions to be addressed cost-effectively and rapidly.

EXPERT OPINION

: Experimental techniques and computational methods are improving rapidly, but each needs a skilled user to maximize results obtained from these techniques. Multidisciplinary teams are therefore key to making optimal use of these methods and such qualified teams can provide enormous benefits to pharmaceutical companies to improve device efficacy and thus time to market. There is already a move to integrate the benefits of Industry 4.0 into inhaler design and usage, a trend that will accelerate.

Large eddy simulation of cough jet dynamics, droplet transport, and inhalability over a ten minute exposure

High fidelity simulations of expiratory events such as coughing provide the opportunity to predict the fate of the droplets from the turbulent jet cloud produced from a cough. It is well established that droplets carrying infectious pathogens with diameters of 1 - 5   μ m remain suspended in the air for several hours and transported by the air currents over considerable distances (e.g., in meters). This study used a highly resolved mesh to capture the multiphase turbulent buoyant cloud with suspended droplets produced by a cough. The cough droplets' dispersion was subjected to thermal gradients and evaporation and allowed to disperse between two humans standing 2 m apart. A nasal cavity anatomy was included inside the second human to determine the inhaled droplets. Three diameter ranges characterized the droplet cloud, < 5   μ m , which made up 93% of all droplets by number; 5 to 100 μm comprised 3%, and > 100   μ m comprising 4%. The results demonstrated the temporal evolution of the cough event, where a jet is first formed, followed by a thermally driven puff cloud with the latter primarily composed of droplets under 5 μm diameter, moving with a vortex string structure. After the initial cough, the data were interpolated onto a more coarse mesh to allow the simulation to cover ten minutes, equivalent to 150 breathing cycles. We observe that the critical diameter size susceptible to inhalation was 0.5   μ m , although most inhaled droplets after 10 min by the second human were approximately 0.8   μ m . These observations offer insight into the risk of airborne transmission and numerical metrics for modeling and risk assessment.

Inhalation and deposition of spherical and pollen particles after middle turbinate resection in a human nasal cavity

Middle turbinate resection significantly alters the anatomy and redistributes the inhaled air. The superior half of the main nasal cavity is opened up, increasing accessibility to the region. This is expected to increase inhalation dosimetry to the region during exposure to airborne particles. This study investigated the influence of middle turbinate resection on the deposition of inhaled pollutants that cover spherical and non-spherical particles (e.g. pollen). A computational model of the nasal cavity from CT scans, and its corresponding post-operative model with virtual surgery performed was created. Two constant flow rates of 5 L/min, and 15 L/min were simulated under a laminar flow field. Inhaled particles including pollen (non-spherical), and a spherical particle with reference density of 1000 kg/m³ were introduced in the surrounding atmosphere. The effect of surgery was most prominent in the less patent cavity side, since the change in anatomy was proportionally greater relative to the original airway space. The left cavity produced an increase in particle deposition at a flow rate of 15 L/min. The main particle deposition mechanisms were inertial impaction, and to a lesser degree gravitational sedimentation. The results are expected to provide insight into inhalation efficiency of different aerosol types, and the likelihood of deposition in different nasal cavity surfaces.

N95 respirator mask breathing leads to excessive carbon dioxide inhalation and reduced heat transfer in a human nasal cavity

Face masks and respirators are used to filter inhaled air, which may contain airborne droplets and high particulate matter (PM) concentrations. The respirators act as a barrier to the inhaled and exhaled air, which may change the nasal airflow characteristics and air-conditioning function of the nose. This study aims to investigate the nasal airflow dynamics during respiration with and without an N95 respirator driven by airflow through the nasal cavity to assess the effect of the respirator on breathing conditions during respiration. To achieve the objective of this study, transient computational fluid dynamics simulations have been utilized. The nasal geometry was reconstructed from high-resolution Computed Tomography scans of a healthy 25-year-old female subject. The species transport method was used to analyze the airflow, temperature, carbon dioxide (CO₂), moisture content (H₂O), and temperature distribution within the nasal cavity with and without an N95 respirator during eight consecutive respiration cycles with a tidal volume of 500 ml. The results demonstrated that a respirator caused excessive CO₂ inhalation by approximately 7 × greater per breath compared with normal breathing. Furthermore, heat and mass transfer in the nasal cavity was reduced, which influences the perception of nasal patency. It is suggested that wearers of high-efficiency masks that have minimal porosity and low air exchange for CO₂ regulation should consider the amount of time they wear the mask.

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.

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

Computational fluid dynamics comparison of impaired breathing function in French bulldogs with nostril stenosis and an examination of the efficacy of rhinoplasty

BACKGROUND

Brachycephalic obstructive airway syndrome (BOAS) in dogs indicates a particular set of upper airway abnormalities found in brachycephalic dogs (e.g., French bulldogs). Stenotic nares is one of the primary BOAS-related abnormalities restricting the functional breathing of affected dogs. For severe stenosis, rhinoplasty is required to increase the accessibility of the external nostril to air; however, the specific improvement from surgery in terms of respiratory physiology and uptake of inhaled air has not been fully elucidated METHOD: This study employed Computational Fluid Dynamics (CFD) simulations to evaluate the effects of different stenotic intensities on airflow patterns in a total of eight French bulldog upper airways. A bulldog with severe stenosis after surgery was included to examine the efficacy of the surgical intervention.

RESULTS

The results showed homogeneous airflow distributions in healthy and mild stenosis cases and significantly accelerated airstreams at the constricted positions in moderate and severe stenosis bulldogs. The airflow resistance was over 20-fold greater in severe stenosis cases than the healthy cases. After surgery, a decrease in airflow velocity was observed in the surgical region, and the percentage of reduced airflow resistance was approximately 4%.

CONCLUSIONS

This study suggests impaired breathing function in brachycephalic dogs with moderate and severe stenosis. The results also serve as a reference for veterinarians in surgical planning and monitoring bulldogs' recuperation after surgery.

Ultra-sensitive photoelectrochemical aptamer biosensor for detecting E. coli O157:H7 based on nonmetallic plasmonic two-dimensional hydrated defective tungsten oxide nanosheets coupling with nitrogen-doped graphene quantum dots (dWO3•H2O@N-GQDs)

Light absorption and interfacial engineering of photoactive materials play vital roles in photoexcited electron generation and electron transport, and ultimately boost the performance of photoelectrochemical (PEC) biosensing. In this work, a novel high-performance photoelectrochemical (PEC) biosensing platform was fabricated based on nonmetallic plasmonic tungsten oxide hydrate nanosheets (WO₃•H₂O) coupling with nitrogen doped graphene quantum dots (N-GQDs) by a facile one-step hydrothermal approach. The localized surface plasmon resonance (LSPR) properties were achieved by oxygen vacancy engineered WO₃·H₂O (dWO₃•H₂O), which could greatly extend the light absorption from visible light to near-infrared light. Moreover, by coupling with N-GQDs, the as-fabricated heterojunction (dWO₃•H₂O@N-GQD) provided a much enhanced photoelectric response due to the efficient charge transfer. By conjugation with E.coli O157:H7 aptamer, a novel PEC aptasensor based on dWO₃•H₂O@N-GQD heterojunction was fabricated with a high sensitivity for detection of E.coli O157:H7. The limit of detection (LOD) of this PEC aptasensor is 0.05 CFU/mL with a linear detection range from 0.1 to 10⁴ CFU/mL. Moreover, high reproducibility and good accuracy could also be achieved for analysis in milk samples. This work could provide a promising platform for the development of PEC bioanalysis and offer an insight into the non-metallic plasmonic materials based heterojunctions for high-performances PEC biosensing.

Numerical analysis of nanoparticle transport and deposition in a cynomolgus monkey nasal passage

Environmental exposure to toxic agents is commonly encountered by occupational and residential populations. However, in vivo exposure data in human subjects is limited by measurement and ethical restrictions. Monkey represents a suitable surrogate for human exposure studies, but the particle transport and deposition features in monkey airways are still not well understood. As a response to this research challenge, this paper presents a virtual exposure study that numerically investigated the nanoparticle transport process through a realistic cynomolgus monkey nasal airway. Particles with size of 1 nm to 1 μm were considered and the transport process was modelled by the Lagrangian discrete phase model. Overall and local deposition as well as particle dispersion along the airway were examined by using a variety of non-dimensional parameters including combined diffusion parameter, deposition enhancement factor and particle flux enhancement factor. Consistent deposition patterns were observed in present and literature nasal models. Most particles tended to pass the nasal airway through certain spatial regions, including the middle section of the nasal valve, the lower half of the middle coronal plane, and the central regions of the choana. While naturally inhaled nanoparticles can hardly be delivered to the olfactory region as it is located apart from the mainstream with high particle flux. Research findings provide insight into nanoparticle inhalation exposure characteristics in the monkey airway and can contribute in formulating data extrapolation schemes between monkey and human airways.

Prediction of nasal spray drug absorption influenced by mucociliary clearance

Evaluation of nasal spray drug absorption has been challenging because deposited particles are consistently transported away by mucociliary clearance during diffusing through the mucus layer. This study developed a novel approach combining Computational Fluid Dynamics (CFD) techniques with a 1-D mucus diffusion model to better predict nasal spray drug absorption. This integrated CFD-diffusion approach comprised a preliminary simulation of nasal airflow, spray particle injection, followed by analysis of mucociliary clearance and drug solute diffusion through the mucus layer. The spray particle deposition distribution was validated experimentally and numerically, and the mucus velocity field was validated by comparing with previous studies. Total and regional drug absorption for solute radius in the range of 1 − 110nm were investigated. The total drug absorption contributed by the spray particle deposition was calculated. The absorption contribution from particles that deposited on the anterior region was found to increase significantly as the solute radius became larger (diffusion became slower). This was because the particles were consistently moved out of the anterior region, and the delayed absorption ensured more solute to be absorbed by the posterior regions covered with respiratory epithelium. Future improvements in the spray drug absorption model were discussed. The results of this study are aimed at working towards a CFD-based integrated model for evaluating nasal spray bioequivalence.

Aerodynamic impact of total inferior turbinectomy versus inferior turbinoplasty - a computational fluid dynamics study

BACKGROUND

The aim of this study was to investigate using computational fluid dynamics (CFD) the effects on nasal aerodynamics of two different techniques for reducing the inferior turbinate. This may assist in surgical planning to select the optimal procedure.

METHODS

Virtual surgery using two techniques of turbinate reduction was performed in eight nasal airway obstruction patients. Three bilateral nasal airway models for each patient were compared: 1) Pre-operative 2) Bilateral inferior turbinoplasty 3) Bilateral total inferior turbinate resection (ITR). Two representative healthy models were included. CFD modeling of airflow was performed under steady-state, laminar, inspiratory conditions.

RESULTS

Nasal airway resistance was slightly more reduced following ITR compared to turbinoplasty due to loss of the pressure gradient at the head of the IT. Turbinoplasty resulted in ventilation, pressure and wall shear stress profiles closer to those of healthy models. A more prominent jet-like course of the main flow stream was observed inferiorly in the ITR group.

CONCLUSIONS

Nasal air conditioning was significantly altered following IT surgery. Overall differences between the groups were small and are unlikely to bear influence on nasal function in normal environments. Further studies using a larger number of patients and healthy subjects are required, attempting to establish a clinical correlation with long-term outcomes such as the perception of nasal patency, mucosal crusting and drying, and air conditioning in different environments. Since a large proportion of IT mucosa remains following turbinoplasty, future dependence on topical therapy should also be considered.

Computational modelling of nasal respiratory flow

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

Quantification of airflow in the sinuses following functional endoscopic sinus surgery

BACKGROUND

Despite functional endoscopic sinus surgery (FESS) being the standard of care in medically recalcitrant chronic rhinosinusitis (CRS), its effect on sinus ventilation has not been fully characterized. Airflow simulations can help improve our understanding of how surgical strategies affect post-surgical sinus ventilation.

METHODS

Eight postoperative sinonasal cavity models were reconstructed from a wide spectrum of CRS patients who had undergone FESS. Computational fluid dynamics modeling of steady-state, laminar, inspiratory airflow was performed. Ventilation was quantified and observed for all the sinuses in each model.

RESULTS

Sinus aeration was enhanced following FESS, particularly in the maxillary and ethmoid sinuses. The degree of improvement was related to the extent of surgery performed. This finding was accentuated at a higher inhalational flow rate of 15L/min. The relationship between ostium size and corresponding sinus inflow was stronger for the maxillary and sphenoid sinuses. Maxillary inflow reached 50% in a mega-antrostomy patient, while negligible flow occurred in the frontal sinuses for except one whom had undergone a modified Lothrop procedure.

CONCLUSIONS

This study has quantified sinus airflow in the largest set of post-FESS patients to date, to show that with increasing extensive surgery, the sinus and nasal cavity become more interconnected and functionally interdependent. Accordingly, sinus ventilation is improved. This may have important consequences for pre- and post-surgical assessment and planning, and on predicting how drug delivery treatments and devices can be designed to target the postoperative sinuses.

Primary break-up and atomization characteristics of a nasal spray

The primary objective of this research was to extract the essential information needed for setting atomization break up models, specifically, the Linear Instability Sheet Atomization (LISA) breakup model, and alternative hollow cone models. A secondary objective was to gain visualization and insight into the atomization break up mechanism caused by the effects of viscosity and surface tension on primary break-up, sheet disintegration, ligament and droplet formation. High speed imaging was used to capture the near-nozzle characteristics for water and drug formulations. This demonstrated more rapid atomization for lower viscosities. Image processing was used to analyze the near-nozzle spray characteristics during the primary break-up of the liquid sheet into ligament formation. Edges of the liquid sheet, spray break-up length, break-up radius, cone angle and dispersion angle were obtained. Spray characteristics pertinent for primary breakup modelling were determined from high speed imaging of multiple spray actuations. The results have established input data for computational modelling involving parametrical analysis of nasal drug delivery.

Magnetic resonance imaging evaluation of the distribution of spray and irrigation devices within the sinonasal cavities

BACKGROUND

Optimizing intranasal distribution and retention of nasal sprays is essential in the management of patients with chronic rhinosinusitis (CRS), including those that have had functional endoscopic sinus surgery (FESS). Despite multiple existing distribution studies, there remains a need for a technique that allows regionalization of particle deposition within a patient's unique 3-dimensional (3D) geometry without exposing the patient to radiation.

METHODS

Seven participants delivered normal saline containing a gadolinium-based contrast agent (GBCA) by either saline irrigation or nasal sprays on 1 side of the nasal cavity. The saline irrigation group included 2 participants (both healthy) while the nasal spray group included 5 participants (2 healthy, 2 post-FESS patients, 1 CRS patient without any sinus surgery). The distribution of new signal enhancement was assessed on each participant using magnetic resonance imaging (MRI). Serial scans were performed over an interval of 4 minutes in the nasal spray group to assess changes in intranasal distribution over time.

RESULTS

Signal enhancement was widespread within the nasal cavities and maxillary sinuses of participants (both healthy) that underwent sinus irrigation. For the nasal spray participants, the hotspots for signal enhancement were similar regardless of disease status or previous history of surgery. These included the internal nasal valve, anterior septum, inferior surface of the inferior turbinate, nasal floor, and nasopharynx. No signal enhancement was detected with nasal sprays in either unoperated or operated paranasal sinuses.

CONCLUSION

A technique has been developed using MRI evaluation of radioopaque contrast to characterize the temporospatial distribution of topical drug delivery within the sinonasal cavities.

Comment on the European position paper on diagnostic tools in rhinology â€" computational fluid dynamics

We read with interest the excellent paper authored by Rimmer et al. (1) entitled The European Position Paper on Diagnostic Tools in Rhinology. The authors are to be commended for their comprehensive, up-to-date and thorough summary of available tests in the assessment of rhinologic function. The paper describes traditional tools used in nasal airway obstruction (NAO); including subjective patient reported outcome measure questionnaires (PROMs) such as the NOSE scale; tests that are subjective to the clinician, such as clinical examination, nasendoscopy and imaging; and objective tests such peak nasal inspiratory flow (PNIF), rhinomanometry (RM) and acoustic rhinometry (AR). Unfortunately, each of these readily available, traditional tools fail to meet several accepted criteria of an ideal diagnostic test (Table 1), as outlined (2,3). In our opinion, these limitations restrict the capacity of Rhinology to develop as a discipline founded on sound evidence-based science.

Characterization of nasal irrigation flow from a squeeze bottle using computational fluid dynamics

BACKGROUND

Nasal saline irrigation has become standard of care in various sinonasal conditions, including allergic and nonallergic rhinitis, chronic rhinosinusitis, and in the postoperative patient. Evidence regarding the mechanisms and dynamics of liquid flow through the sinonasal cavity remains limited due to inadequate experimental models (cadaveric, 3-dimensional [3D] printed, imaging of labeled dyes and radioisotopes). We aimed to develop a computational fluid dynamics (CFD) model of nasal irrigation to demonstrate sinonasal surface coverage, residence times across the mucosal surfaces, and shearing force of irrigation.

METHODS

A nasal cavity geometry derived from high-resolution paranasal sinus computed tomography (CT) scans of a healthy, unoperated, 25-year-old patient was created. CFD analysis was performed to assess the distribution of nasal irrigation from a tapered nozzle bottle at a forward head-tilt position of 45 degrees with a 2-second burst at 35 mL/second.

RESULTS

The model demonstrates nasal irrigation from ipsilateral to contralateral with precise measures of velocity, pressure, wall shear stress, and mapping of surface coverage and residence times at specific locations and times. The nasal cavity experiences almost complete coverage of irrigation, while overflow from the nasal cavity facilitates moderate coverage of the ipsilateral maxillary (40%) and anterior ethmoid sinuses (30%). Negligible coverage of the sphenoid and frontal sinuses was noted.

CONCLUSION

Detailed physical mechanisms of liquid irrigation injected from a commonly used squeeze bottle were shown. Ipsilateral maxillary and ethmoid sinus penetration are primarily due to overflow rather than direct jet entry, confirming the recommendation of larger volumes of irrigation to "flood" the sinus ostia.

In silico approaches to respiratory nasal flows: A review

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

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

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

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

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

Computational investigation of dust mite allergens in a realistic human nasal cavity

Aim: Inhaled allergens from house dust mite (HDM) are a major source of allergic disease such as allergic rhinitis and asthma. It has been a challenge to properly evaluate health risks caused by HDM related allergens including mite bodies, eggs and fecal pellets. This paper presents a numerical study on particle deposition of dust mite allergens in a human nasal cavity. Materials and methods: A realistic nasal cavity model was reconstructed from CT scans and a Computational Fluid Dynamics analysis of steady airflow was simulated. The discrete phase model was used to trace particle trajectories of three dust mite related particles. Results: The flow and particle model were validated by comparing with nasal resistance measurement and previous literature respectively. Aerodynamic characteristics and deposition of dust mite allergens in the nasal cavity were analyzed under different breathing conditions including rest and exercising conditions. Conclusions: The numerical results revealed the roles of different nasal cavity regions in filtering various types of dust mite allergens with consideration of breathing conditions.

Ultrafine particle deposition in a realistic human airway at multiple inhalation scenarios

The scarcity of regional deposition data in distal respiratory airways represents an important challenge for current toxicology and pharmacology research. To bridge this gap, a realistic airway model extending from nasal and oral openings to distal bronchial airways with varying pathway length was built in this study. Transport and deposition characteristics of naturally inhaled ultrafine particles (UFPs) ranging from 1 to 100 nm were numerically investigated, and effects of different inhalation scenarios were considered. To enable intercase particle deposition comparison, an adjusted parameter, unified deposition enhancement factor (UDEF), was proposed for quantifying the localised deposition concentration. Results show that compartment particle deposition peaked around the ultrafine end of the considered size range, and it dropped rapidly with the increase of particle size. Different inhalation modes caused notable deposition changes in the extrathoracic region, while its effects in the TB airway are much less. For UFPs larger than 10 nm, predicted deposition efficiencies in all compartments are all at lowest levels among considered particle size range, implying UFPs ranging from 10 to 100 nm can travel through the whole respiratory airway model and escape to the alveolar region. Furthermore, high enhancement factors were observed at the vicinity of most bifurcation apexes, and more even UDEF distribution was observed from 1-nm particle cases. While for 100-nm cases, the deposited particles tend to concentrate at few "hot spots" (areas of high deposition concentration in relation to surrounding surfaces) with greater UDEF in the tracheobronchial airway.

Correlation of regional deposition dosage for inhaled nanoparticles in human and rat olfactory

Background

Nose-to-brain transport of airborne ultrafine particles (UFPs) via the olfactory pathway has been verified as a possible route for particle translocation into the brain. The exact relationship between increased airborne toxicant exposure and neurological deterioration in the human central nervous system, is still unclear. However, the nasal olfactory is undoubtedly a critical junction where the time course and toxicant dose dependency might be inferred.

Method

Computational fluid-particle dynamics modeling of inhaled nanoparticles (1 to 100 nm) under low to moderate breathing conditions (5 to 14 L/min – human; and 0.14 to 0.40 L/min – rat) were performed in physiologically realistic human and rat nasal airways. The simulation emphasized olfactory deposition, and variations in airflow and particle flux caused by the inter-species airway geometry differences. Empirical equations were developed to predict regional deposition rates of inhaled nanoparticles on human and rat olfactory mucosa in sedentary breathing. Considering, breathing and geometric differences, quantified correlations between human and the rat olfactory deposition dose against a variety of metrics were proposed.

Results

Regional deposition of nanoparticles in human and the rat olfactory was extremely low, with the highest deposition (< 3.5 and 8.1%) occurring for high diffusivity particles of 1.5 nm and 5 nm, respectively. Due to significant filtering of extremely small particles (< 2 nm) by abrupt sharp turns at front of the rat nose, only small fractions of the inhaled nanoparticles (in this range) reached rat olfactory than that in human (1.25 to 45%); however, for larger sizes (> 3 nm), significantly higher percentage of the inhaled nanoparticles reached rat nasal olfactory than that in human (2 to 32 folds). Taking into account the physical and geometric features between human and rat, the total deposition rate (#/min) and deposition rate per unit surface area (#/min/mm²) were comparable for particles> 3 nm. However, when body mass was considered, the normalized deposition rate (#/min/kg) in the rat olfactory region exceeded that in the human. Nanoparticles < 1.5 nm were filtered out by rat anterior nasal cavity, and therefore deposition in human olfactory region exceeded that in the rat model.

Conclusion

Regional deposition dose of inhaled nanoparticles in a human and rat olfactory region was governed by particle size and the breathing rate. Interspecies correlation was determined by combining the effect of deposition dosage, physical\geometric features, and genetic differences. Developed empirical equations provided a tool to quantify inhaled nanoparticle dose in human and rat nasal olfactory regions, which lay the ground work for comprehensive interspecies correlation between the two species. Furthermore, this study contributes to the fields in toxicology, i.e., neurotoxicity evaluation and risk assessment of UFPs, in long-term and low-dose inhalation exposure scenarios.

Antennal scales improve signal detection efficiency in moths

The elaborate bipectinate antennae of male moths are thought to increase their sensitivity to female sex pheromones, and so should be favoured by selection. Yet simple filamentous antennae are the most common structure among moths. The stereotypic arrangements of scales on the surface of antennae may resolve this paradox. We use computational fluid dynamics techniques to model how scales on the filamentous antennae of moths affect the passage of different particles in the airflow across the flagellum in both small and large moths. We found that the scales provide an effective solution to improve the efficacy of filamentous antennae, by increasing the concentration of nanoparticles, which resemble pheromones, around the antennae. The smaller moths have a greater increase in antennal efficiency than larger moths. The scales also divert microparticles, which resemble dust, away from the antennal surface, thereby reducing contamination. The positive correlations between antennal scale angles and sensilla number across Heliozelidae moths are consistent with the predictions of our model.

Development of a computational fluid dynamics model for mucociliary clearance in the nasal cavity

Intranasal drug delivery has attracted significant attention because of the opportunity to deliver systemic drugs directly to the blood stream. However, the mucociliary clearance poses a challenge in gaining high efficacy of intranasal drug delivery because cilia continuously carry the mucus blanket towards the laryngeal region. To better understand mucus flow behaviour on the human nasal cavity wall, we present computational model development, and evaluation of mucus motion on a realistic nasal cavity model reconstructed from CT-scans. The model development involved two approaches based on the actual nasal cavity geometry namely: (i) unwrapped-surface model in 2D domain and (ii) 3D-shell model. Conservation equations of fluid motion were applied to the domains, where a mucus production source term was used to initiate the mucus motion. The analysis included mucus flow patterns, virtual saccharin tests and quantitative velocity magnitude analysis, which demonstrated that the 3D-shell model results provided better agreement with experimental data. The unwrapped-surface model also suffered from mesh-deformations during the unwrapping stage and this led to higher mucus velocity compared to experimental data. Therefore, the 3D-shell model was recommended for future mucus flow simulations. As a first step towards mucus motion modelling this study provides important information that accurately simulates a mucus velocity field on a human nasal cavity wall, for assessment of toxicology and efficacy of intranasal drug delivery.

Detailed computational analysis of flow dynamics in an extended respiratory airway model

BACKGROUND

Understanding respiratory physiology can aid clinicians in diagnosing the cause of respiratory symptoms or shed light on drug delivery inhaler device optimisation. However, the sheer complexity of the human lung prohibits a full-scale study.

METHODS

In this study, a realistic respiratory airway model including large-to-small conducting airways was built. This airway model consists of subject-specific upper and lower airways, extending from nasal and oral openings to terminal bronchioles (up to the 15th generation). Based on the subject-specific airway model, topological information was extracted and a digital reference model that exhibits strong asymmetry and multi-fractal properties was provided. Inhalation flow rates 18 L/min and 50 L/min were adopted to understand inspiratory conditions subjecting to resting and light exercise inhalation modes. Regional airflow in terms of axial velocity and secondary flow vortices along the lung airway model was extracted.

FINDINGS

Obvious secondary flow currents were seen in the larynx-trachea segment and left main bronchus, while for the terminal conducting airway in the right lower lobe, the airflow tends to be much smoother with no secondary flow currents.

INTERPRETATION

This paper provides insights on respiratory physiology, especially in the lower lung airways, and will be potentially useful for diagnosis of lower airway diseases.

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.

Air conditioning analysis among human nasal passages with anterior anatomical variations

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

Detailed deposition analysis of inertial and diffusive particles in a rat nasal passage

Rats have been widely used as surrogates for evaluating the health effects of inhaled airborne particulate matter. To provide a thorough understanding of particle transport and deposition mechanisms in the rat nasal airway, this article presents a computational fluid dynamics (CFD) study of particle exposure in a realistic rat nasal passage under a resting flow condition. Particles covering a diameter range from 1 nm to 4 µm were passively released in front of the rat's breathing zone, and the Lagrangian particle tracking approach was used to calculate individual particle trajectories. Detailed particle deposition analysis shows the deposition of inertial particles >2 µm is high in the rat nasal vestibule and more than 70% of all inhaled inertial particles were trapped in this region. While for diffusive nanoparticles, the vestibule filtration effect is reduced, only less than 60% of inhaled nanoparticles were blocked by the anterior nasal structures. The particle exposure in the olfactory region only shows notable deposition for diffusive nanoparticles, which peaks at 9.4% for 5 nm particles. Despite the olfactory deposition remains at a low level, the ratio between the olfactory and the main passage is kept around 30-40% for 10-800 nm particles, which indicates a particle-size-independent distribution pattern in the main nasal passage and olfactory. This study provides a deep understanding of particles deposition features in a rat nasal passage, and the research findings can aid toxicologist in inter-species exposure-response extrapolation study.

Effect of morphology on nanoparticle transport and deposition in human upper tracheobronchial airways

Aerosol transport and deposition in human lungs has attracted considerable attention in the past few years, as it has significant value to the study of toxicity consequence as well as therapeutic potential in occupational health and medical applications. In reproducing human tracheobronchial airways, two approaches were frequently taken: (1) anatomical realistic reconstruction through image scans (e.g. CT and MRI) or cadaver casts and (2) mathematical description using simplified models. Strengths and limitations are primarily focused on accuracy, resolution, repeatability, and computational\physical expenses. While both approaches were reported in literature, detailed comparison of aerosol transport and deposition in the two representations were scarcely performed, largely due to the challenge to acquire comprehensive data from the irregular structured airway replicas (approach 1). To fill the gap, the current study performed a numerical comparison of nanoparticle transport and deposition in human upper tracheobronchial airways by using an anatomical realistic reconstruction (through CT scans) and a mathematically simplified airway model. As the first step, the current study was focused on the variation in breathing airflow pattern and the effect towards fate of the inhaled nanoparticles in human upper tracheobronchial airways. The study provided important information to understand geometric sensitivity of nanoparticle modeling in the human tracheobronchial tree and is of significant value to predict the whole lung uptake of inhaled nanoparticles in the human respiratory system.

Numerical Comparison of Nasal Aerosol Administration Systems for Efficient Nose-to-Brain Drug Delivery

PURPOSE

Nose-to-brain drug administration along the olfactory and trigeminal nerve pathways offers an alternative route for the treatment of central nervous system (CNS) disorders. The characterization of particle deposition remains difficult to achieve in experiments. Alternative numerical approach is applied to identify suitable aerosol particle size with maximized inhaled doses.

METHODS

This study numerically compared the drug delivery efficiency in a realistic human nasal cavity between two aerosol drug administration systems targeting the olfactory region: the aerosol mask system and the breath-powered bi-directional system. Steady inhalation and exhalation flow rates were applied to both delivery systems. The discrete phase particle tracking method was employed to capture the aerosol drug transport and deposition behaviours in the nasal cavity. Both overall and regional deposition characteristics were analysed in detail.

RESULTS

The results demonstrated the breath-powered drug delivery approach can produce superior olfactory deposition with peaking olfactory deposition fractions for diffusive 1 nm particles and inertial 10 μm. While for particles in the range of 10 nm to 2 μm, no significant olfactory deposition can be found, indicating the therapeutic agents should avoid this size range when targeting the olfactory deposition.

CONCLUSIONS

The breath-powered bi-directional aerosol delivery approach shows better drug delivery performance globally and locally, and improved drug administration doses can be achieved in targeted olfactory region.