Modeling the pharyngeal anatomical effects on breathing resistance and aerodynamically generated sound

Optimized gravity-driven intranasal drop administration delivers significant doses to the ostiomeatal complex and maxillary sinus

Chronic and allergic rhinosinusitis impacts approximately 12% of the global population. Challenges in rhinosinusitis treatment include paranasal sinus inaccessibility and variability in delivery efficiency among individuals. This study addresses these challenges of drug delivery by developing a high-efficiency, low-variability protocol for nasal drop delivery to the ostiomeatal complex (OMC) and maxillary sinus. Patient-specific nasal casts were dissected to reveal the configurations of conchae and meatus, providing insights into anatomical features amendable for sinus delivery. Fluorescent dye-enhanced videos visualized the dynamic liquid translocation in transparent nasal casts, allowing real-time assessment and quick adjustment to delivery parameters. Dosimetry to the OMC and maxillary sinus were quantified as drop count and mass using a precision scale. Key delivery factors, including the device type, formulation, and head-chin orientation, were systematically investigated in a cohort of ten nasal casts. Results show that both the squeeze bottle and soft-mist nasal pump yielded notably low doses to the OMC with high variability, and no dose from these two devices was detected within the maxillary sinuses. In contrast, the proposed approach, which included a curved nozzle surpassing the nasal valve and leveraged gravity-driven liquid translocation along the lateral nasal wall, delivered significant doses to the OMC and maxillary sinus. Iterative experimentations identified the optimal head tilt to be 40° and chin tilt to be° from the lateral recumbent position. Statistical analyses established the drop count required for effective OMC/sinus delivery. The proposed delivery protocol holds the potential to enhance chronic rhinosinusitis treatment outcomes with low variability. The dual role of nasal anatomy in posing challenges and offering opportunities highlights the need for future investigations using diverse formulations in a larger cohort of nasal models.

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.

Study on the flow mechanism and frequency characteristics of rales in lower respiratory tract

The frequency characteristics of lung sounds have great significance for noninvasive diagnosis of respiratory diseases. The rales in the lower respiratory tract region that can provide rich information about symptoms of respiratory diseases are not clear. In this paper, a three-dimensional idealized bifurcated lower respiratory tract geometric model, which contains 3rd to 13th generation (G3-G13) bronchi is constructed, where Re ∼ 10 1 - 10 3 , and then the large eddy simulation and volume of fluid are used to study the fluid flow characteristics. Ffowcs Williams and Hawkings model are subsequently used to study the frequency characteristics of rale of different generations of bronchi. The results showed that bronchial blockage and sputum movement will enhance the turbulence intensity and vortex shedding intensity of flow. The dominant frequency and highest value of sound pressure level (SPL) of rhonchi/moist crackles decrease with the increase of bronchial generation. The change rates of dominant frequency of rhonchi / moist crackles in adjacent generations were 5.0 ± 0.1 ~ 9.1 ± 0.2% and 3.1 ± 0.1 ~ 11.9 ± 0.3%, respectively, which is concentrated in 290 ~ 420 Hz and 200 ~ 300 Hz, respectively. The change rates of SPL of rhonchi/moist crackles were 8.8 ± 0.1 ~ 15.7 ± 0.1% and 7.1 ± 0.1 ~ 19.5 ± 0.2%, respectively, which is concentrated in 28 ~ 50 dB and 16 ~ 32 dB, respectively. In the same generation of bronchus (e.g., G8, G9) with the same degree of initial blockage, the dominant frequency and SPL of moist crackles can be 3.7 ± 0.2% and 4.5 ± 0.3% slightly higher than that of rhonchi, respectively. This research is conducive to the establishment of a rapid and accurate noninvasive diagnosis system for respiratory diseases.

Simulating the effect of individual upper airway anatomical features on drug deposition

This study aims to systematically isolate different anatomical features of the human pharynx with the goal to investigate their independent influence on airflow dynamics and particle deposition characteristics in a geometrically realistic human airway. Specifically, the effects of the uvula, epiglottis and soft palate on drug particle deposition are studied systematically, by carefully removing each of these anatomical features from reconstructed models based on MRI data and comparing them to a benchmark realistic airway model. Computational Fluid Dynamics using established turbulence models is employed to simulate the transport of mono-dispersed particles (3 µm) in the airway at two flow-rates. The simulations suggest three findings: 1) widening the space between the oral cavity and oropharynx and where the soft palate is situated leads to the most dramatic reduction in drug deposition in the upper airway; 2) exclusion of the uvula and epiglottis: a) affects flow dynamics in the airway; b) alters regional deposition behaviour; c) does not significantly affect the total number of particles deposited in the pharynx; and 3) the space adjacent to the soft palate is a key determinant for aerosol deposition in the extrathoracic region and is related to mechanisms of flow acceleration, diversion and recirculation.

Evaluation of Impulse Oscillometry in Respiratory Airway Casts with Varying Obstruction Phenotypes, Locations, and Complexities

The use of impulse oscillometry (IOS) for lung function testing does not need patient cooperation and has gained increasing popularity among both young and senior populations, as well as in patients with breathing difficulties. However, studies of the IOS sensitivity to regional lung obstructions are limited and have shown mixed results. The objective of this study was to evaluate the performance of an IOS system in 3D-printed lung models with structural abnormalities at different locations and with different severities. Lung trees of two complexity levels were tested, with one extending to the sixth generation (G6) and the other to G12. The IOS responses to varying glottal apertures, carina ridge tumors, and segmental bronchial constrictions were quantified in the G6 lung geometry. Both the G6 and G12 lung casts were prepared using high-resolution 3D printers. Overall, IOS detected the progressive airway obstructions considered in this study. The resonant frequency dropped with increasing obstructions for all three disease phenotypes in the G6 lung models. R20Hz increased with the increase in airway obstructions. Specifically, R20Hz in the airway model with varying glottal apertures agreed reasonably well with complementary measurements using TSI VelociCalc. In contrast to the high-resistance (R) sensitivity to the frequency in G6 lung models, R was nearly independent of frequency in G12 lung models. IOS R20Hz demonstrated adequate sensitivity to the structural remodeling in the central airways. However, the changes of R5Hz and X5Hz vs. airway obstructions were inconclusive in this study, possibly due to the rigid lung casts and the difference of a container–syringe system from human lungs.

Inhalation dosimetry of nasally inhaled respiratory aerosols in the human respiratory tract with locally remodeled conducting lungs

Objective: Respiratory diseases are often accompanied by alterations to airway morphology. However, inhalation dosimetry data in remodeled airways are scarce due to the challenges in reconstructing diseased respiratory morphologies. This study aims to study the airway remodeling effects on the inhalation dosimetry of nasally inhaled nanoparticles in a nose-lung geometry that extends to G9 (ninth generation).Materials and methods: Statistical shape modeling was used to develop four diseased lung models with varying levels of bronchiolar dilation/constriction in the left-lower (LL) lobe (i.e. M1-M4). Respiratory airflow and particle deposition were simulated using a low Reynolds number k-ω turbulence model and a Lagrangian tracking approach.Results: Significant discrepancies were observed in the flow partitions between the left and right lungs, as well as between the lower and upper lobes of the left lung, which changed by 10-fold between the most dilated and constricted models.Much lower doses were predicted on the surface of the constricted LL bronchioles G4-G9, as well as into the peripheral airways beyond G9 of the LL lung. However, the LL lobar remodeling had little effect on the dosimetry in the nasopharynx, as well as on the total dosimetry in the nose-lung geometry (up to G9).Conclusion: It is suggested that airway remodeling may pose a higher viral infection risk to the host by redistributing the inhaled viruses to healthy lung lobes. Airway remodeling effects should also be considered in the treatment planning of inhalation therapies, not only because of the dosimetry variation from altered lung morphology but also its evolution as the disease progresses.

Septal destruction enhances chaotic mixing and increases cellular doses of nanoparticles in emphysematous acinus

One hallmark of emphysema is the breakdown of inter-alveolar septal walls in pulmonary acini. How the acinar dosimetry of environmental aerosols varies at different stages of emphysema remains unclear; this is specifically pertinent to users of tobacco products, which is the leading cause of emphysema. The objective of this study is to systematically assess the impacts of septal destruction on the behavior and fate of nanoparticles (1–800 nm) in a pyramid-shaped sub-acinar model consisting of 496 alveoli. Four diseased geometry variants were created by gradually removing the septal walls from the base model. Particle motions within the acinar region were tracked for particles raging 1–800 nm at four emphysema stages using a well-tested Lagrangian tracking model. Both spatial profile and temporal variation of particle deposition were predicted in healthy and diseased sub-acinar geometries on both a total and regional basis. Results show large differences in airflow and particle dynamics among different emphysema stages. Large differences in particle dynamics are also observed among different particle sizes, with one order of magnitude’s variation in the speeds of particles of 1, 10, and 200 nm. The destruction of septal walls also changed the deposition mechanisms, shifting from connective diffusion to chaotic mixing with emphysema progression. The sub-acinar dosimetry became less sensitive to particle size variation with more septal destructions. The lowest retention rate was found at 200–500 nm in the healthy sub-acinar geometry, but at 800 nm in all emphysematous models considered. The acinus-averaged dose for nanoparticles (1–800 nm) increases with aggravating septal destructions, indicating an even higher risk to the acinus at later emphysema stages.

SARS COV-2 virus-laden droplets coughed from deep lungs: Numerical quantification in a single-path whole respiratory tract geometry

When an infected person coughs, many virus-laden droplets will be exhaled out of the mouth. Droplets from deep lungs are especially infectious because the alveoli are the major sites of coronavirus replication. However, their exhalation fraction, size distribution, and exiting speeds are unclear. This study investigated the behavior and fate of respiratory droplets (0.1-4 μm) during coughs in a single-path respiratory tract model extending from terminal alveoli to mouth opening. An experimentally measured cough waveform was used to control the alveolar wall motions and the flow boundary conditions at lung branches from G2 to G18. The mouth opening was modeled after the image of a coughing subject captured using a high-speed camera. A well-tested k-ω turbulence model and Lagrangian particle tracking algorithm were applied to simulate cough flow evolutions and droplet dynamics under four cough depths, i.e., tidal volume ratio (TVR) = 0.13, 0.20. 0.32, and 0.42. The results show that 2-μm droplets have the highest exhalation fraction, regardless of cough depths. A nonlinear relationship exists between the droplet exhalation fraction and cough depth due to a complex deposition mechanism confounded by multiscale airway passages, multiregime flows, and drastic transient flow effects. The highest exhalation fraction is 1.6% at the normal cough depth (TVR = 0.32), with a mean exiting speed of 20 m/s. The finding that most exhaled droplets from deep lungs are 2 μm highlights the need for more effective facemasks in blocking 2-μm droplets and smaller both in infectious source control and self-protection from airborne virus-laden droplets.

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

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

Effects of mask-wearing on the inhalability and deposition of airborne SARS-CoV-2 aerosols in human upper airway

Even though face masks are well accepted as tools useful in reducing COVID-19 transmissions, their effectiveness in reducing viral loads in the respiratory tract is unclear. Wearing a mask will significantly alter the airflow and particle dynamics near the face, which can change the inhalability of ambient particles. The objective of this study is to investigate the effects of wearing a surgical mask on inspiratory airflow and dosimetry of airborne, virus-laden aerosols on the face and in the respiratory tract. A computational model was developed that comprised a pleated surgical mask, a face model, and an image-based upper airway geometry. The viral load in the nose was particularly examined with and without a mask. Results show that when breathing without a mask, air enters the mouth and nose through specific paths. When wearing a mask, however, air enters the mouth and nose through the entire surface of the mask at lower speeds, which favors the inhalation of ambient aerosols into the nose. With a 65% filtration efficiency (FE) typical for a three-layer surgical mask, wearing a mask reduces dosimetry for all micrometer particles except those of size 1 µm-3 µm, for which equivalent dosimetry with and without a mask in the upper airway was predicted. Wearing a mask reduces particle penetration into the lungs, regardless of the FE of the mask. The results also show that mask-wearing protects the upper airway (particularly the nose and larynx) best from particles larger than 10 µm while protecting the lungs best from particles smaller than 10 µm.

Alveolar size effects on nanoparticle deposition in rhythmically expanding-contracting terminal alveolar models

Significant differences in alveolar size exist in humans of different ages, gender, health, and among different species. The effects of alveolar sizes, as well as the accompanying breathing frequencies, on regional and local dosimetry of inhaled nanoparticles have not been sufficiently studied. Despite a well-accepted qualitative understanding of the advection-diffusion-sedimentation mechanism in the acinar region, a quantitative picture of the interactions among these factors remains inchoate. The objective of this study is to quantify the effects of alveolar size on the regional and local deposition of inhaled nanoparticles in alveolar models of varying complexities and to understand the dynamic interactions among different deposition mechanisms. Three different models were considered that retained 1, 4, and 45 alveoli, respectively. For each model, the baseline geometry was scaled by ¼, ½, 2, 4, and 8 times by volume. Temporal evolution and spatial distribution of particle deposition were tracked using a discrete-phase Lagrangian model. Lower retentions of inhaled nanoparticles were observed in the larger alveoli under the same respiration frequency, while similar retentions were found among different geometrical scales if breathing frequencies allometrically matched the alveolar size. Dimensional analysis reveals a manifold deposition mechanism with tantamount contributions from advection, diffusion, and gravitational sedimentation, each of which can become dominant depending on the location in the alveoli. Results of this study indicate that empirical correlations obtained from one sub-population cannot be directly applied to others, nor can they be simply scaled as a function of the alveolar size or respiration frequency due to the regime-transiting deposition mechanism that is both localized and dynamic.

Ventilation Modulation and Nanoparticle Deposition in Respiratory and Olfactory Regions of Rabbit Nose

The rabbit nose's ability to filter out inhaled agents is directly related to its defense to infectious diseases. The knowledge of the rabbit nose anatomy is essential to appreciate its functions in ventilation regulation, aerosol filtration and olfaction. The objective of this study is to numerically simulate the inhalation and deposition of nanoparticles in a New Zealand white (NZW) rabbit nose model with an emphasis on the structure-function relation under normal and sniffing conditions. To simulate the sniffing scenario, the original nose model was modified to generate new models with enlarged nostrils or vestibules based on video images of a rabbit sniffing. Ventilations into the maxilloturbinate and olfactory region were quantified with varying nostril openings, and deposition rates of inhaled aerosols ranging from 0.5 nm to 1000 nm were characterized on the total, sub-regional and local basis. Results showed that particles which deposited in the olfactory region came from a specific area in the nostril. The spiral vestibule played an essential role in regulating flow resistance and flow partition into different parts of the nose. Increased olfactory doses were persistently predicted in models with expanded nostrils or vestibule. Particles in the range of 5-50 nm are more sensitive to the geometry variation than other nanoparticles. It was also observed that exhaled aerosols occupy only the central region of the nostril, which minimized the mixing with the aerosols close to the nostril wall, and potentially allowed the undisruptive sampling of odorants. The results of this study shed new light on the ventilation regulation and inhalation dosimetry in the rabbit nose, which can be further implemented to studies of infectious diseases and immunology in rabbits.

Nanoparticle Deposition in Rhythmically Moving Acinar Models with Interalveolar Septal Apertures

Pulmonary delivery of nanomedicines has been extensively studied in recent years because of their enhanced biocompatibility, sustained-release properties, and surface modification capability. The lung as a target also offers many advantages over other routers, such as large surface area, noninvasive, quick therapeutic onset, and avoiding first-pass metabolism. However, nanoparticles smaller than 0.26 µm typically escape phagocytosis and remain in the alveoli for a long time, leading to particle accumulation and invoking tissue responses. It is imperative to understand the behavior and fates of inhaled nanoparticles in the alveoli to reliably assess therapeutic outcomes of nanomedicines or health risk of environmental toxins. The objective of this study is to numerically investigate nanoparticle deposition in a duct-alveolar model with varying sizes of inter-alveolar septal apertures (pores). A discrete phase Lagrangian model was implemented to track nanoparticle trajectories under the influence of rhythmic wall expansion and contraction. Both temporal and spatial dosimetry in the alveoli were computed. Wall motions are essential for nanoparticles to penetrate the acinar region and deposit in the alveoli. The level of aerosol irreversibility (i.e., mixing of inhaled nanoparticles with residual air in the alveolar airspace) is determined by the particle diffusivity, which in turn, dictates the fraction of particles being exhaled out. When deposition in the upper airways was not considered, high alveolar deposition rates (74–95%) were predicted for all nanoparticles considered (1–1000 nm), which were released into the alveoli at the beginning of the inhalation. The pore size notably affects the deposition pattern of inhaled nanoparticles but exerts a low impact upon the total deposition fractions. This finding indicates that consistent pulmonary doses of nanomedicine are possible in emphysema patients if breathing maneuver with the same tidal volume can be performed.

Effects of Increased Pharyngeal Tissue Mass Due to Fluid Accumulation in the Neck on the Acoustic Features of Snoring Sounds in Men

STUDY OBJECTIVES

Snoring sounds are generated by the vibration of pharyngeal tissue due to the upper airway narrowing. While recorded by a microphone placed over the neck, snoring can pass through the pharyngeal tissue surrounding the upper airway. Thus, changes in the pharyngeal tissue content may change the acoustic properties of the snoring sounds. Rostral fluid shift and the consequent increases in neck fluid volume (NFV) and neck circumference (NC) can increase pharyngeal tissue mass. Therefore, the goal of this study was to investigate the relationship between increases in pharyngeal tissue mass, as assessed by increased NFV and NC, and snoring sounds features.

METHODS

We obtained data from a previous study where 20 males who were not obese participated in a daytime polysomnography and their NC and NFV were measured before and after sleep. During sleep, snoring sounds were recorded with a microphone placed over the neck. Spectral centroid of the snoring sounds was estimated. Then, the first five snoring segments were selected from the first and last 30 minutes of stage N2 sleep.

RESULTS

We found a significant decrease in the snoring spectral centroid from the beginning to end of sleep. We also found that spectral centroid from the end of sleep in frequency ranges below 200 Hz was inversely correlated with the increases in NFV and NC from before to after sleep.

CONCLUSIONS

These results suggest that snoring spectral centroid can be used as a noninvasive and convenient method to assess variations in the pharyngeal tissue mass.

Airflow and Particle Deposition in Acinar Models with Interalveolar Septal Walls and Different Alveolar Numbers

Unique features exist in acinar units such as multiple alveoli, interalveolar septal walls, and pores of Kohn. However, the effects of such features on airflow and particle deposition remain not well quantified due to their structural complexity. This study aims to numerically investigate particle dynamics in acinar models with interalveolar septal walls and pores of Kohn. A simplified 4-alveoli model with well-defined geometries and a physiologically realistic 45-alveoli model was developed. A well-validated Lagrangian tracking model was used to simulate particle trajectories in the acinar models with rhythmically expanding and contracting wall motions. Both spatial and temporal dosimetries in the acinar models were analyzed. Results show that collateral ventilation exists among alveoli due to pressure imbalance. The size of interalveolar septal aperture significantly alters the spatial deposition pattern, while it has an insignificant effect on the total deposition rate. Surprisingly, the deposition rate in the 45-alveoli model is lower than that in the 4-alveoli model, indicating a stronger particle dispersion in more complex models. The gravity orientation angle has a decreasing effect on acinar deposition rates with an increasing number of alveoli retained in the model; such an effect is nearly negligible in the 45-alveoli model. Breath-holding increased particle deposition in the acinar region, which was most significant in the alveoli proximal to the duct. Increasing inhalation depth only slightly increases the fraction of deposited particles over particles entering the alveolar model but has a large influence on dispensing particles to the peripheral alveoli. Results of this study indicate that an empirical correlation for acinar deposition can be developed based on alveolar models with reduced complexity; however, what level of geometry complexity would be sufficient is yet to be determined.

Deposition of bolus and continuously inhaled aerosols in rhythmically moving terminal alveoli

The particle dynamics in an oscillating alveolus under tidal breathing can be dramatically different from those in a static alveolus. Despite its close relevance to pulmonary drug delivery and health risk from airborne exposure, quantifications of alveolar deposition are scarce due to its inaccessibility to in vivo measurement instruments, tiny size to replicate in vitro, and dynamic wall motions to model. The objective of this study is to introduce a numerical method to quantify alveolar deposition with continuous particle release in a rhythmically oscillating alveolus by integrating the deposition curves for bolus aerosols and use this method to develop correlations applicable in assessing alveolar drug delivery efficiency or dosimetry of inhaled toxicants. An idealized blind-end terminal alveolus model was developed with rhythmically moving alveolar boundary conditions in phase with tidal breathing. The dynamic wall expansion mode and magnitude were based on experimentally measured chest wall motions and tidal volumes. A well-validated Lagrangian tracking model was used to simulate the transport and deposition of inhaled micrometer particles. Large differences were observed between dynamic and static alveoli in particle motion, deposition onset, and final alveolar deposition fraction. Alveolar deposition of bolus aerosols is highly sensitive to breath-holding duration, particle release time, and alveolar dimension. For 1 µm particles, there exists a cut-off release time (zero bolus deposition), which decreases with alveolar size (i.e., 1.0 s in a 0.2-mm-diameter alveolus and 0.56 s in a 0.8-mm-diameter alveolus). The cumulative alveolar deposition was predicted to be 39% for a 0.2-mm-diameter alveolus, 22% for a 0.4-mm-diameter alveolus, and 10% for a 0.8-mm-diameter alveolus. A cumulative alveolar deposition correlation was developed for inhalation delivery with a prescribed period of drug release and the second correlation for the time variation of alveolar deposition of ambient aerosols, both of which captured the relative dependence of the particle release time and alveolar dimension.

Numerical study of dynamic glottis and tidal breathing on respiratory sounds in a human upper airway model

BACKGROUND

Human snores are caused by vibrating anatomical structures in the upper airway. The glottis is a highly variable structure and a critical organ regulating inhaled flows. However, the effects of the glottis motion on airflow and breathing sound are not well understood, while static glottises have been implemented in most previous in silico studies. The objective of this study is to develop a computational acoustic model of human airways with a dynamic glottis and quantify the effects of glottis motion and tidal breathing on airflow and sound generation.

METHODS

Large eddy simulation and FW-H models were adopted to compute airflows and respiratory sounds in an image-based mouth-lung model. User-defined functions were developed that governed the glottis kinematics. Varying breathing scenarios (static vs. dynamic glottis; constant vs. sinusoidal inhalations) were simulated to understand the effects of glottis motion and inhalation pattern on sound generation. Pressure distributions were measured in airway casts with different glottal openings for model validation purpose.

RESULTS

Significant flow fluctuations were predicted in the upper airways at peak inhalation rates or during glottal constriction. The inhalation speed through the glottis was the predominating factor in the sound generation while the transient effects were less important. For all frequencies considered (20-2500 Hz), the static glottis substantially underestimated the intensity of the generated sounds, which was most pronounced in the range of 100-500 Hz. Adopting an equivalent steady flow rather than a tidal breathing further underestimated the sound intensity. An increase of 25 dB in average was observed for the life condition (sine-dynamic) compared to the idealized condition (constant-rigid) for the broadband frequencies, with the largest increase of approximately 40 dB at the frequency around 250 Hz.

CONCLUSION

Results show that a severely narrowing glottis during inhalation, as well as flow fluctuations in the downstream trachea, can generate audible sound levels.

Optimization of magnetophoretic-guided drug delivery to the olfactory region in a human nose model

Magnetophoretic-guided delivery has been shown to be able to improve the olfactory doses. However, due to the complex nasal structure and quick decay of magnetic intensity, precise control of particle motion in the human nose remains a challenge. In this study, an optimization model was developed for magnetophoretic olfactory delivery systems. The performance of the model was evaluated using a baseline device design in an MRI-based human nose geometry. Three key components of the delivery system were examined, which included the particle release position, the front magnet to minimize nasal valve depositions, and the top magnet to attract particles into the olfactory region. Results show that the magnetophoretic olfactory delivery device can be significantly improved by optimizing the product and operational parameters. The olfactory delivery efficiency was increased by 1.5-fold compared to the baseline design. The top magnet height and strength were shown to be the most influential factor in olfactory delivery, followed by the drug release position and the front magnet strength. The optimization framework developed in this study can be easily adapted for the optimization of intranasal drug delivery to other regions such as paranasal sinuses.

Improving intranasal delivery of neurological nanomedicine to the olfactory region using magnetophoretic guidance of microsphere carriers

BACKGROUND

Although direct nose-to-brain drug delivery has multiple advantages, its application is limited by the extremely low delivery efficiency (<1%) to the olfactory region where drugs can enter the brain. It is crucial to developing new methods that can deliver drug particles more effectively to the olfactory region.

MATERIALS AND METHODS

We introduced a delivery method that used magnetophoresis to improve olfactory delivery efficiency. The performance of the proposed method was assessed numerically in an image-based human nose model. Influences of the magnet layout, magnet strength, drug-release position, and particle diameter on the olfactory dosage were examined.

RESULTS AND DISCUSSION

Results showed that particle diameter was a critical factor in controlling the motion of nasally inhaled ferromagnetic drug particles. The optimal particle size was found to be approximately 15 μm for effective magnetophoretic guidance while avoiding loss of particles to the walls in the anterior nose. Olfactory delivery efficiency was shown to be sensitive to the position and strength of magnets and the release position of drug particles. The results of this study showed that clinically significant olfactory doses (up to 45%) were feasible using the optimal combination of magnet layout, selective drug release, and microsphere-carrier diameter. A 64-fold-higher delivery of dosage was predicted in the magnetized nose compared to the control case, which did not have a magnetic field. However, the sensitivity of olfactory dosage to operating conditions and the unstable nature of magnetophoresis make controlled guidance of nasally inhaled aerosols still highly challenging.