Laminar airflow and nanoparticle or vapor deposition in a human nasal cavity model

The transport and deposition of nanoparticles, i.e., dp = 1-2 nm, or equivalent vapors, in the human nasal cavities is of interest to engineers, scientists, air-pollution regulators, and healthcare officials alike. Tiny ultrafine particles, i.e., dp < or = 5 nm, are of special interest because they are most rapidly absorbed and hence have an elevated toxic or therapeutic impact when compared to larger particles. Assuming transient laminar 3-D incompressible flow in a representative human nasal cavity, the cyclic airflow pattern as well as local and overall nanoparticle depositions were computationally simulated and analyzed. The focus was on transient effects during inhalation/exhalation as compared to the steady-state assumption typically invoked. Then, an equation for a matching steady-state inhalation flow rate was developed that generates the same deposition results as cyclic inhalation. Of special interest is the olfactory region where the narrow channel surfaces receive only about one-half of a percent of the inhaled nanoparticles because the airflow bypasses these recesses located in the superior-most portions in the geometrically complex nasal cavities.     

Particle Deposition in Oral-tracheal Airway Models with Very Low Inhalation Profiles

Considerable progress has been made on modeling particle deposition in the oral-tracheal airway under some normal breathing conditions,i.e.,resting,light activity and moderate exercise.None of these standard breathing patterns correspond to very low inhalation profiles.It is known that particle deposition in the oral-tracheal airway is greatly influenced by flow and particle inlet conditions.In this work,very low inhalation flow rates are considered.Particle deposition is numerically investigated in different oral-tracheal airway models,i.e.,circular,elliptic and realistic oral-tracheal airway models.Both micro- and nano-particles that are normally present in cigarette smoke are considered.Results show that inhalation profiles greatly influence the particle deposition.Due to relatively low flow rate,for ultra-fine particles,the oral deposition is enhanced due to longer residence time in oral cavity and stronger Brownian motion.However,for larger particles,less particles deposit in the oral-tracheal airway due to the weaker impaction.The transition happens when particle size changes from 0.01 μm to 0.1 μm.The influence of the limited entrance area is shown and discussed.Under the low inhalation profiles,the highest deposition fraction could be in either circular or realistic models depending on the particle property and the geometric characteristic of oral cavity.The knowledge obtained in this study may be beneficial for the design of bionic inhaler and understanding of health effect from smoke particle on human being.     

Numerical modeling of turbulent and laminar airflow and odorant transport during sniffing in the human and rat nose

Human sniffing behavior usually involves bouts of short, high flow rate inhalation (>300 ml/s through each nostril) with mostly turbulent airflow. This has often been characterized as a factor enabling higher amounts of odorant to deposit onto olfactory mucosa than for laminar airflow and thereby aid in olfactory detection. Using computational fluid dynamics human nasal cavity models, however, we found essentially no difference in predicted olfactory odorant flux (g/cm2 s) for turbulent versus laminar flow for total nasal flow rates between 300 and 1000 ml/s and for odorants of quite different mucosal solubility. This lack of difference was shown to be due to the much higher resistance to lateral odorant mass transport in the mucosal nasal airway wall than in the air phase. The simulation also revealed that the increase in airflow rate during sniffing can increase odorant uptake flux to the nasal/olfactory mucosa but lower the cumulative total uptake in the olfactory region when the inspired air/odorant volume was held fixed, which is consistent with the observation that sniff duration may be more important than sniff strength for optimizing olfactory detection. In contrast, in rats, sniffing involves high-frequency bouts of both inhalation and exhalation with laminar airflow. In rat nose odorant uptake simulations, it was observed that odorant deposition was highly dependent on solubility and correlated with the locations of different types of receptors.     

Analytical model of hygroscopic particle behavior in human airways

A model is developed to calculate the deposition of hygroscopic aerosols in the human tracheobronchial (TB) tree. The TB airflow pattern assumed is consistent with experimental observations and accounts for anatomical features such as the larynx and cartilaginous rings in large airways. Some original deposition efficiency formulae are presented for laminar and turbulent airstreams. Stepwise growth is simulated by changes in particle size and density at each TB generation. The dose distribution of NaCl aerosols is studied as a function of inhaled particle size and flow rate. Two NaCl growth rate curves are used which differ in the mode of aerosol-air mixing in the trachea. The initial rate of aerosol mixing in the human due to the laryngeal jet is shown to be an important factor affecting the deposition of hygroscopic aerosols. Total TB deposition of NaCl exceeds that for nonhygroscopic particles of the same inhaled aerodynamic size. Hygroscopic growth can also influence the regional TB distribution of dose when submicron NaCl particles grow rapidly enough to deposit by impaction and sedimentation.     

The effect of inlet flow profile and nozzle diameter on drug delivery to the maxillary sinus

In this paper, the effect of the turbulence and swirling of the inlet flow and the diameter of the nozzle on the flow characteristics and the particles' transport/deposition patterns in a realistic combination of the nasal cavity (NC) and the maxillary sinus (MS) were examined. A computational fluid dynamics (CFD) model was developed in ANSYS® Fluent using a hybrid Reynolds averaged Navier–Stokes–large-eddy simulation algorithm. For the validation of the CFD model, the pressure distribution in the NC was compared with the experimental data available in the literature. An Eulerian–Lagrangian approach was employed for the prediction of the particle trajectories using a discrete phase model. Different inlet flow conditions were investigated, with turbulence intensities of 0.15 and 0.3, and swirl numbers of 0.6 and 0.9 applied to the inlet flow at a flow rate of 7 L/min. Monodispersed particles with a diameter of 5 µm were released into the nostril for various nozzle diameters. The results demonstrate that the nasal valve plays a key role in nasal resistance, which damps the turbulence and swirl intensities of the inlet flow. Moreover, it was found that the effect of turbulence at the inlet of the NC on drug delivery to the MS is negligible. It was also demonstrated that increasing the flow swirl at the inlet and decreasing the nozzle diameter improves the total particle deposition more than threefold due to the generation of the centrifugal force, which acts on the particles in the nostril and vestibule. The results also suggest that the drug delivery efficiency to the MS can be increased by using a swirling flow with a moderate swirl number of 0.6. It was found that decreasing the nozzle diameter can increase drug delivery to the proximity of the ostium in the middle meatus by more than 45%, which subsequently increases the drug delivery to the MS. The results can help engineers design a nebulizer to improve the efficiency of drug delivery to the maxillary sinuses.

     

Numerical simulation of airflow and microparticle deposition in a synchrotron micro-CT-based pulmonary acinus model

The acinus consists of complex, branched alveolar ducts and numerous surrounding alveoli, and so in this study, we hypothesized that the particle deposition can be much influenced by the complex acinar geometry, and simulated the airflow and particle deposition (density = 1.0 g/cm³, diameter = 1 and 3 μm) numerically in a pulmonary acinar model based on synchrotron micro-CT of the mammalian lung. We assumed that the fluid–structure interaction was neglected and that alveolar flow was induced by the expansion and contraction of the acinar model with the volume changing sinusoidally with time as the moving boundary conditions. The alveolar flow was dominated by radial flows, and a weak recirculating flow was observed at the proximal side of alveoli during the entire respiratory cycle, despite the maximum Reynolds number at the inlet being 0.029. Under zero gravity, the particle deposition rate after single breathing was less than 0.01, although the particles were transported deeply into the acinus after inspiration. Under a gravitational field, the deposition rate and map were influenced strongly by gravity orientation. In the case of a particle diameter of 1 μm, the rate increased dramatically and mostly non-deposited particles remained in the model, indicating that the rate would increase further after repeated breathing. At a particle diameter of 3 μm, the rate was 1.0 and all particles were deposited during single breathing. Our results show that the particle deposition rate in realistic pulmonary acinar model is higher than in an idealized model.     

Theoretical study of nasal deposition

A theory is derived to calculate the regional and total deposition of aerosol particles in the nasal passages during inhalation. The particle size studied range from 0.2 to 10.0 μm diameter. The deposition is calculated in five regions; (I) the region filled with nasal hair, (II) the nasal valve, (III) the expansion region, (IV) the turbinate region and (V) the posterior bend. Equations are derived to determine the deposition caused by direct impaction on the nasal hairs and bends of the passages. The calculations show the deposition due to direct impaction does not account for the amount or location of deposited particles measured in experiments. Secondary flows have been speculated to exist in the expansion region after the nasal valve and an equation is derived to estimate the deposition caused by the secondary flows. The calculated deposition, due to direct impaction and secondary flows, shows general agreement with the experiment as to the predicted amount and location of deposited particles.     

Nasal aerodynamics protects brain and lung from inhaled dust in subterranean diggers,Ellobius talpinus

Inhalation of air-dispersed sub-micrometre and nano-sized particles presents a risk factor for animal and human health. Here, we show that nasal aerodynamics plays a pivotal role in the protection of the subterranean mole vole Ellobius talpinus from an increased exposure to nano-aerosols. Quantitative simulation of particle flow has shown that their deposition on the total surface of the nasal cavity is higher in the mole vole than in a terrestrial rodent Mus musculus (mouse), but lower on the olfactory epithelium. In agreement with simulation results, we found a reduced accumulation of manganese in olfactory bulbs of mole voles in comparison with mice after the inhalation of nano-sized MnCl₂ aerosols. We ruled out the possibility that this reduction is owing to a lower transportation from epithelium to brain in the mole vole as intranasal instillations of MnCl₂ solution and hydrated nanoparticles of manganese oxide MnO · (H₂O)_x revealed similar uptake rates for both species. Together, we conclude that nasal geometry contributes to the protection of brain and lung from accumulation of air-dispersed particles in mole voles.     

Particle deposition and sensory drive

The mutualism between chemical cues emitted into the air and variations in how primates respond to them using olfaction has demonstrated aspects of species‐specific adaptations. Building on this mutualism we can look at particle deposition as another means to understanding how various environments may have elicited biological changes that enable efficient communication. Research on particle movement and deposition within the nasal cavity is largely based on questions about health as it relates to drug delivery systems and overall olfactory function in modern humans. With increased access to 3D models and the use of computational fluid dynamic analysis, researchers have been able to simulate site‐specific deposition, to determine what particles are making it through the nasal cavity to the main olfactory epithelium, which ultimately leads to processing in the olfactory bulb. Here we discuss particle deposition research, sensory drive and their potential applications to evolutionary anthropology.     

Detailed flow patterns in the nasal cavity.

The human nasal cavity filters and conditions inspired air while providing olfactory function. Detailed experimental study of nasal airflow patterns has been limited because of the complex geometry of the nasal cavity. In this work, particle image velocimetry was used to determine two-dimensional instantaneous velocity vector fields in parallel planes throughout a model of the nasal cavity that was subjected to a nonoscillatory flow rate of 125 ml/s. The model, which was fabricated from 26 computed tomography scans by using rapid prototyping techniques, is a scaled replica of a human right nasal cavity. The resulting vector plots show that the flow is laminar and regions of highest velocity are in the nasal valve and in the inferior airway. The relatively low flow in the olfactory region appears to protect the olfactory bulb from particulate pollutants. Low flows were also observed in the nasal meatuses, whose primary function has been the subject of debate. Comparison of sequentially recorded data suggests a steady flow.     

Modeling inspiratory and expiratory steady-state velocity fields in the Sprague-Dawley rat nasal cavity

Distribution patterns of odorant molecules in the rat nasal olfactory region depend in large part on the detailed airflow patterns in the nasal cavity, which in turn depend on the anatomical structure. To investigate these flow patterns, we constructed an anatomically accurate finite element model of the right nasal cavity of the Sprague-Dawley rat based on horizontal (anterior-posterior) nasal cast cross sections. By numerically solving the fluid mechanical momentum and continuity equations using the finite element method, we studied the flow distribution and the complete velocity field for both inspiration and expiration throughout the nasal cavity under physiological flow rates of resting breathing and sniffing. Detailed velocity profiles, volumetric flow distributions, and streamline patterns for quasi-steady airflow are presented. S-shaped streamlines passing through the olfactory region are found to be less prevalent during expiratory than inspiratory flow leading to trapping and an increase in odorant molecule retention in the olfactory region during sniffing. The rat nasal velocity calculations will be used to study the distribution of odorant uptake onto the rat olfactory mucosa and compare it with the known anatomic location of some types of rat olfactory receptors.     

Numerical Simulation of Airway Dimension Effects on Airflow Patterns and Odorant Deposition Patterns in the Rat Nasal Cavity

The sense of smell is largely dependent on the airflow and odorant transport in the nasal cavity, which in turn depends on the anatomical structure of the nose. In order to evaluate the effect of airway dimension on rat nasal airflow patterns and odorant deposition patterns, we constructed two 3-dimensional, anatomically accurate models of the left nasal cavity of a Sprague-Dawley rat: one was based on high-resolution MRI images with relatively narrow airways and the other was based on artificially-widening airways of the MRI images by referencing the section images with relatively wide airways. Airflow and odorant transport, in the two models, were determined using the method of computational fluid dynamics with finite volume method. The results demonstrated that an increase of 34 µm in nasal airway dimension significantly decreased the average velocity in the whole nasal cavity by about 10% and in the olfactory region by about 12% and increased the volumetric flow into the olfactory region by about 3%. Odorant deposition was affected to a larger extent, especially in the olfactory region, where the maximum odorant deposition difference reached one order of magnitude. The results suggest that a more accurate nasal cavity model is necessary in order to more precisely study the olfactory function of the nose when using the rat.     

Flow simulation in the human upper respiratory tract

Computer simulations of airflow patterns within the human upper respiratory tract (URT) are presented. The URT model includes airways of the head (nasal and oral), throat (pharyngeal and laryngeal), and lungs (trachea and main bronchi). The head and throat morphology was based on a cast of a medical school teaching model; tracheobronchial airways were defined mathematically. A body-fitted three-dimensional curvilinear grid system and a multiblock method were employed to graphically represent the surface geometries of the respective airways and to generate the corresponding mesh for computational fluid dynamics simulations. Our results suggest that for a prescribed phase of breath (i.e., inspiration or expiration), convective respiratory airflow patterns are highly dependent on flow rate values. Moreover, velocity profiles were quite different during inhalation and exhalation, both in terms of the sizes, strengths, and locations of localized features such as recirculation zones and air jets. Pressure losses during inhalation were 30-35% higher than for exhalation and were proportional to the square of the flow rate. Because particles are entrained and transported within airstreams, these results may have important applications to the targeted delivery of inhaled drugs.     

Respiratory deposition model of an inhaled aerosol bolus

The bolus delivery method is designed to deliver a dose to the desired location in the lung, and it has the advantage of fewer side effects and a more efficient way of delivery. Based upon the lung deposition model developed for continuously inhaling aerosols of constant concentration, a mathematical model of aerosol bolus deposition is proposed. The calculated results show that the recovery depends on the bolus penetration depth, flow rate, particle size, breath holding time and bolus volume. Three sets of published experimental data with different controlling factors (particle size, flow rate and breath holding time) are adopted to make the quantitative comparisons with the calculated results. The predictions and data for the low intrinsic motion particles (～1 μm) have good agreement, as do the coarse particles in the shallow airways region. For females, the recovery was found to be consistently lower than that for males.     

A Simplified Geometric Model to Predict Nasal Spray Deposition in Children and Adults

A mathematical approach was developed to estimate spray deposition patterns in the nasal cavity based on the geometric relationships between the emitted spray plume and the anatomical dimensions of the nasal valve region of the nasal cavity. Spray plumes were assumed to be spherical cones and the nasal valve region was approximated as an ellipse. The effect of spray plume angle (15–85°) on the fraction of the spray able to pass through the nasal valve (deposition fraction) was tested for a variety of nasal valve (ellipse) shapes and cross-sectional areas based on measured dimensions from pediatric and adult nasal cavities. The effect of the distances between the tip of the nasal spray device and the nasal valve (0.2–1.9 cm) on the deposition fraction was also tested. Simulation results show that (1) decreasing spray plume angles resulted in higher deposition fractions, (2) deposition fraction was inversely proportional to the spray distance and the nasal valve (ellipse) major/minor axis ratio, and (3) for fixed major/minor axis ratios, improved deposition occurred with larger nasal valve cross-sectional areas. For a typical adult nasal valve, plume angles of less than 40° emitted from a distance of 1 cm resulted depositions greater than 90% within the main nasal cavity, whereas for a 12-year-old child, only the most narrow plume angles (<?20°) administered resulted in significant deposition beyond the nasal valve.     

Deposit feeding in Abra tenuis (Bivalvia: Tellinacea)

Deposit feeding in Abra tenuis is described in terms of the size of particles utilized.
Material is collected by the inhalant siphon performing a circular motion sucking in sediment from beneath and on the surface.
The size distribution of silica admitted into the mantle cavity is described and shown to be controlled by physical parameters. The density of a particle does not affect its uptake by the inhalant siphon. The size distribution of the sediment affects the size distribution of particles admitted to the mantle cavity.
No selection of material for ingestion in terms of size occurs after it has been taken into the mantle cavity. Thus the range of material ingested is ultimately controlled by a physical parameter, the inhalant opening.     

Computational Analysis of Micron-particle Deposition in a Human Triple Bifurcation Airway Model

Steady laminar axisymmetric inhalation flow and wall deposition of micron-size particles in representative triple bifurcation airways have been simulated using a commercial finite-volume code with user-enhanced programs. Assuming spherical non-interacting particles (3 μm≤ d p ≤7 μm), various inlet Reynolds numbers (Re=500-2000) and Stokes numbers (St=0.02-0.23) were considered. The resulting particle deposition patterns were analyzed and then summarized in terms of deposition efficiencies, i.e. DE=DE(Re,St) Surprisingly high DE-values occur at relatively low Reynolds numbers (e.g., Re=500 ) in the third bifurcation. The quantitative results are of interest to researchers either conducting health risk assessment studies for inhaled particulate pollutants or analyzing drug aerosol inhalation and deposition at desired lung target sites.     

Determining the basic characteristics of aerosols suitable for studies of deposition in the respiratory tract

Studies of aerosol particle deposition in the respiratory tract requires experimental inhalation of artificial model aerosols. The paper formulates some of the most important requirements for the properties of such aerosols. Several suitable fractions were prepared as part of a research project dealing with the use of microporous polymers for diagnostic purposes. 5 fractions of the polymer designated G-gel 60 with the particle size as stated by the manufacturer, ranging from 3 to 7 micron were evaluated using a 16-channel particle dispersity analyzer HIAC/ROYCO MT 3210 with the sensor 1200 and operated by a microprocessor, the equipment being coupled to an APPLE IIe computer. G-gel 60 particles introduced into the aerosol were characterized by the parameters CMAD, MMAD and sg both numerically and graphically. The measurement procedure was found to be very sensitive with respect to all fractions in evaluating the subtile differences between different lot numbers of the aerosol. G-gel 60 fractions characterized both numerically and graphically were compared with the known aerosols from paraffin oil and atmospheric air. The equipment MT 3210 enables prompt determination of the percentages of aerosol particles distribution by size class. The authors conclude that the procedure, both in its numerical and graphical versions, is particularly suitable for the diagnosis of aerosol particles deposition in the respiratory tract, offering a new application for HIAC/ROYCO in the field of medicine. In evaluating atmospheric aerosol in exhaled air, the number of particles was found to be below that in inhaled air, the difference being dependent on the choice of investigation methods. Percentual distribution of deposited particles following one minute ventilation proved to be at its maximum, as regards atmospheric aerosol, in the 0.30-0.50 micron range. The deposition curve was similar to already published curves, being characterized by an S-shaped pattern with maximum deposition in the greater size classes. An analysis of inhaled, exhaled and deposited aerosol suggested that deposited aerosol is more polydisperse and has particles of greater sizes than inhaled aerosol. Investigation of the effect of apnoe on deposition indicated that deposition increased as a function of apnoeic pause.     

Computational analysis of micron-particle deposition in a human triple bifurcation airway model

Steady laminar axisymmetric inhalation flow and wall deposition of micron-size particles in representative triple bifurcation airways have been simulated using a commercial finite-volume code with user-enhanced programs. Assuming spherical non-interacting particles (3 microm相似文献   

Computational fluid dynamics simulations of particle deposition in large-scale, multigenerational lung models

Computational fluid dynamics (CFD) has emerged as a useful tool for the prediction of airflow and particle transport within the human lung airway. Several published studies have demonstrated the use of Eulerian finite-volume CFD simulations coupled with Lagrangian particle tracking methods to determine local and regional particle deposition rates in small subsections of the bronchopulmonary tree. However, the simulation of particle transport and deposition in large-scale models encompassing more than a few generations is less common, due in part to the sheer size and complexity of the human lung airway. Highly resolved, fully coupled flowfield solution and particle tracking in the entire lung, for example, is currently an intractable problem and will remain so for the foreseeable future. This paper adopts a previously reported methodology for simulating large-scale regions of the lung airway (Walters, D. K., and Luke, W. H., 2010, "A Method for Three-Dimensional Navier-Stokes Simulations of Large-Scale Regions of the Human Lung Airway," ASME J. Fluids Eng., 132(5), p. 051101), which was shown to produce results similar to fully resolved geometries using approximate, reduced geometry models. The methodology is extended here to particle transport and deposition simulations. Lagrangian particle tracking simulations are performed in combination with Eulerian simulations of the airflow in an idealized representation of the human lung airway tree. Results using the reduced models are compared with those using the fully resolved models for an eight-generation region of the conducting zone. The agreement between fully resolved and reduced geometry simulations indicates that the new method can provide an accurate alternative for large-scale CFD simulations while potentially reducing the computational cost of these simulations by several orders of magnitude.