A combined experimental and numerical study on upper airway dosimetry of inhaled nanoparticles from an electrical discharge machine shop

BACKGROUNDS

Exposure to nanoparticles in the workplace is a health concern to occupational workers with increased risk of developing respiratory, cardiovascular, and neurological disorders. Based on animal inhalation study and human lung tumor risk extrapolation, current authoritative recommendations on exposure limits are either on total mass or number concentrations. Effects of particle size distribution and the implication to regional airway dosages are not elaborated.

METHODS

Real time production of particle concentration and size distribution in the range from 5.52 to 98.2 nm were recorded in a wire-cut electrical discharge machine shop (WEDM) during a typical working day. Under the realistic exposure condition, human inhalation simulations were performed in a physiologically realistic nasal and upper airway replica. The combined experimental and numerical study is the first to establish a realistic exposure condition, and under which, detailed dose metric studies can be performed. In addition to mass concentration guided exposure limit, inhalation risks to nano-pollutant were reexamined accounting for the actual particle size distribution and deposition statistics. Detailed dosimetries of the inhaled nano-pollutants in human nasal and upper airways with respect to particle number, mass and surface area were discussed, and empirical equations were developed.

RESULTS

An astonishing enhancement of human airway dosages were detected by current combined experimental and numerical study in the WEDM machine shop. Up to 33 folds in mass, 27 folds in surface area and 8 folds in number dosages were detected during working hours in comparison to the background dosimetry measured at midnight. The real time particle concentration measurement showed substantial emission of nano-pollutants by WEDM machining activity, and the combined experimental and numerical study provided extraordinary details on human inhalation dosimetry. It was found out that human inhalation dosimetry was extremely sensitive to real time particle concentration and size distribution. Averaged particle concentration over 24-h period will inevitably misrepresent the sensible information critical for realistic inhalation risk assessment.

CONCLUSIONS

Particle size distribution carries very important information in determining human airway dosimetry. A pure number or mass concentration recommendation on the exposure limit at workplace is insufficient. A particle size distribution, together with the deposition equations, is critical to recognize the actual exposure risks. In addition, human airway dosimetry in number, mass and surface area varies significantly. A complete inhalation risk assessment requires the knowledge of toxicity mechanisms in response to each individual metric. Further improvements in these areas are needed.

From the Cover: Comparative Numerical Modeling of Inhaled Nanoparticle Deposition in Human and Rat Nasal Cavities

To gain a better understanding of nanoparticle exposure in human nasal cavities, laboratory animals (e.g. rat) are used for in vivo studies. However, due to anatomical differences between human and rodent nasal cavities, direct particle deposition comparisons between species are difficult. This paper presents a comparative nanoparticle (1 nm, 10 nm, and 100 nm) deposition study using anatomically realistic models of a human and rat nasal cavity. The particle deposition fraction was highest consistently in the main nasal passage, for all nanoparticles tested, in the human model; whereas this was only the case for 10 nm, and 100 nm particles for the rodent model, where greater deposition was found in the anterior nose for 1 nm particles. A deposition intensity (DI) term was introduced to represent the accumulated deposition fraction on cross-sectional slices. A common and preferential deposition site in the human model was found for all nanoparticles occurring at a distance of 3.5 cm inside the nasal passage. For the rodent model maximum DI occurred in the vestibule region at a distance of 0.3 cm, indicating that the rodent vestibule produces exceptionally high particle filtration capability. We also introduced a deposition flux which was a ratio of the regional deposition fraction relative to the region's surface area fraction. This value allowed direct comparison of deposition flux between species, and a regional extrapolation scaling factor was found (e.g. 1/10 scale for vestibule region for rat to human comparison). This study bridges the in vitro exposure experiments and in vivo nanomaterials toxicity studies, and can contribute towards improving inter-species exposure extrapolation studies in the future.

Transport and Deposition of Welding Fume Agglomerates in a Realistic Human Nasal Airway

Welding fume is a complex mixture containing ultra-fine particles in the nanometer range. Rather than being in the form of a singular sphere, due to the high particle concentration, welding fume particles agglomerate into long straight chains, branches, or other forms of compact shapes. Understanding the transport and deposition of these nano-agglomerates in human respiratory systems is of great interest as welding fumes are a known health hazard. The neurotoxin manganese (Mn) is a common element in welding fumes. Particulate Mn, either as soluble salts or oxides, that has deposited on the olfactory mucosa in human nasal airway is transported along the olfactory nerve to the olfactory bulb within the brain. If this Mn is further transported to the basal ganglia of the brain, it could accumulate at the part of the brain that is the focal point of its neurotoxicity. Accounting for various dynamic shape factors due to particle agglomeration, the current computational study is focused on the exposure route, the deposition pattern, and the deposition efficiency of the inhaled welding fume particles in a realistic human nasal cavity. Particular attention is given to the deposition pattern and deposition efficiency of inhaled welding fume agglomerates in the nasal olfactory region. For particles in the nanoscale, molecular diffusion is the dominant transport mechanism. Therefore, Brownian diffusion, hydrodynamic drag, Saffman lift force, and gravitational force are included in the model study. The deposition efficiencies for single spherical particles, two kinds of agglomerates of primary particles, two-dimensional planar and straight chains, are investigated for a range of primary particle sizes and a range of number of primary particles per agglomerate. A small fraction of the inhaled welding fume agglomerates is deposited on the olfactory mucosa, approximately in the range 0.1-1%, and depends on particle size and morphology. The strong size dependence of the deposition in olfactory mucosa on particle size implies that the occupation deposition of welding fume manganese can be expected to vary with welding method.

Comparative numerical modeling of inhaled micron-sized particle deposition in human and rat nasal cavities

Micron-sized particle deposition in anatomically realistic models of a rat and human nasal cavity was numerically investigated. A steady laminar inhalation flow rate was applied and particles were released from the outside air. Particles showing equivalent total particle deposition fractions were classified into low, medium and high inertial particle. Typical particle sizes are 2.5, 9 and 20 μm for the human model and 1, 2 and 3 μm for the rat model, respectively. Using a surface-mapping technique the 3D nasal cavity surface was "unwrapped" into a 2D domain and the particle deposition locations were plotted for complete visual coverage of the domain surface. The total surface area comparison showed that the surface area of the human nasal model was about ten times the size of the rat model. In contrast, the regional surface area percentage analysis revealed the olfactory region of the rat model was significantly larger than all other regions making up ∼55.6% of the total surface area, while that of the human nasal model only occupying 10.5%. Flow pattern comparisons showed rapid airflow acceleration was found at the nasopharynx region and the nostril region for the human and rat model, respectively. For the human model, the main passage is the major deposition region for micro-particles. While for the rat model, it is the vestibule. Through comparing the regional deposition flux between human and rat models, this study can contribute towards better extrapolation approach of inhalation exposure data between inter-subject species.

Measurements of droplet size distribution and analysis of nasal spray atomization from different actuation pressure

BACKGROUND

To evaluate the deposition efficiency of spray droplets in a nasal cavity produced from a spray device, it is important to determine droplet size distribution, velocity, and its dispersion during atomization. Due to the limiting geometric dimensions of the nasal cavity airway, the spray plume cannot develop to its full size inside the nasal vestibule to penetrate the nasal valve region for effective drug deposition.

METHODS

Particle/droplet image analysis was used to determine local mean droplet sizes at eight regions within the spray plume under different actuation pressures that represent typical hand operation from pediatric to adult patients.

RESULTS

The results showed that higher actuation pressure produces smaller droplets in the atomization. Stronger actuation pressure typical of adult users produces a longer period of the fully atomized spray stage, despite a shorter overall spray duration. This produces finer droplets when compared with the data obtained by weaker actuation pressure, typical of pediatric users.

CONCLUSION

The experimental technique presented is able to capture a more complete representation of the droplet size distribution and the atomization process during an actuation. The measured droplet size distribution produced can be related to the empirically defined deposition efficiency curve of the nasal cavity, allowing a prediction of the likely deposition.

Fluid-structure interaction analysis of the left coronary artery with variable angulation

The aim of this study is to elucidate the correlation between coronary artery branch angulation, local mechanical and haemodynamic forces at the vicinity of bifurcation. Using a coupled fluid-structure interaction (FSI) modelling approach, five idealized left coronary artery models with various angles ranging from 70° to 110° were developed to investigate the influence of branch angulations. In addition, one CT image-based model was reconstructed to further demonstrate the medical application potential of the proposed FSI coupling method. The results show that the angulation strongly alters its mechanical stress distribution, and the instantaneous wall shear stress distributions are substantially moderated by the arterial wall compliance. As high tensile stress is hypothesized to cause stenosis, the left circumflex side bifurcation shoulder is indicated to induce atherosclerotic changes with a high tendency for wide-angled models.

Nonlinear eddy viscosity modeling and experimental study of jet spreading rates

Indoor airflow pattern is strongly influenced by turbulent shear and turbulent normal stresses that are responsible for entrainment effects and turbulence-driven secondary motion. Therefore, an accurate prediction of room airflows requires reliable modeling of these turbulent quantities. The most widely used turbulence models include RANS-based models that provide quick solutions but are known to fail in turbulent free shear and wall-affected flows. In order to cope with this deficiency, this study presents a nonlinear k-ε turbulence model and evaluates it along with linear k-ε models for an indoor isothermal linear diffuser jet flow measured in two model rooms using PIV. The results show that the flow contains a free jet near the inlet region and a wall-affected region downstream where the jet is pushed toward the ceiling by entrainment through the well-known Coanda effect. The CFD results show that an accurate prediction of the entrainment process is very important and that the nonlinear eddy viscosity model is able to predict the turbulence-driven secondary motions. Furthermore, turbulence models that are calibrated for high Reynolds free shear layer flows were not able to reproduce the measured velocity distributions, and it is suggested that the model constants of turbulence models should be adjusted before they are used for room airflow simulations.

Surface mapping for visualization of wall stresses during inhalation in a human nasal cavity

Airflow analysis can assist in better understanding the physiology however the human nasal cavity is an extremely complicated geometry that is difficult to visualize in 3D space, let alone in 2D space. In this paper, an anatomically accurate 3D surface of the nasal passages derived from CT data was unwrapped and transformed into a 2D space, into a UV-domain (where u and v are the coordinates) to allow a complete view of the entire wrapped surface. This visualization technique allows surface flow parameters to be analyzed with greater precision. A UV-unwrapping tool is developed and a strategy is presented to allow deeper analysis to be performed. This includes (i) the ability to present instant comparisons of geometry and flow variables between any number of different nasal cavity models through normalization of the 2D unwrapped surface; (ii) visualization of an entire surface in one view and; (iii) a planar surface that allows direct 1D and 2D analytical solutions of diffusion of inhaled vapors and particles through the nasal walls. This work lays a foundation for future investigations that correlates adverse and therapeutic health responses to local inhalation of gases and particles.

External characteristics of unsteady spray atomization from a nasal spray device

The nasal route presents an enormous opportunity to exploit the highly vascularized respiratory airway for systemic drug delivery to provide more rapid onset of therapy and reduced drug degradation compared with conventional oral routes. The dynamics of atomization at low injection pressure is less known as typical spray atomization studies have focused on industrial applications such as fuel injection that are performed at much higher pressure. An experimental test station was designed in house and an alternative method to characterize the external spray is presented. This involved the use of high-speed camera to capture the temporal development of the spray as it is atomized through actuation of the spray device. An image-processing technique based on edge detection was developed to automate processing through the large number of images captured. The results showed that there are three main phases of spray development (prestable, stable, and poststable) that can be correlated by examining the spray width. A comparison with a human nasal cavity is made to put into perspective the dimensions and geometry that the spray atomization produces. This study aimed to extend the current existing set of data to contribute toward a better understanding in nasal spray drug delivery.

Local deposition fractions of ultrafine particles in a human nasal-sinus cavity CFD model

Ultrafine particle deposition studies in the human nasal cavity regions often omit the paranasal sinus regions. Because of the highly diffusive nature of nanoparticles, it is conjectured that deposition by diffusion may occur in the paranasal sinuses, which may affect the residual deposition fraction that leaves the nasal cavity. Two identical CFD models of a human nasal cavity, one with sinuses and one without, were reconstructed from CT-scans to determine the uptake of ultrafine particles. In general, there was little flow passing through the paranasal sinuses. However, flow patterns revealed that some streamlines reached the upper nasal cavity near the olfactory regions. These flow paths promote particle deposition in the sphenoid and ethmoid sinuses. It was found that there were some differences in the deposition fractions and patterns for 5 and 10 nm particles between the nasal-sinus and the nasal cavity models. This difference is amplified when the flow rate is decreased and at a flow rate of 4 L/min the maximum difference was 17%. It is suggested that evaluations of nanoparticle deposition should consider some deposition occurring in the paranasal sinuses especially if flow rates are of concern.

Numerical modelling of nanoparticle deposition in the nasal cavity and the tracheobronchial airway

Recent advances in nanotechnology have seen the manufacture of engineered nanoparticles for many commercial and medical applications such as targeted drug delivery and gene therapy. Transport of nanoparticles is mainly attributed to the Brownian force which increases as the nanoparticle decreases to 1 nm. This paper first verifies a Lagrangian Brownian model found in the commercial computational fluid dynamics software Fluent before applying the model to the nasal cavity and the tracheobronchial (TB) airway tree with a focus on drug delivery. The average radial dispersion of the nanoparticles was 9x greater for the user-defined function model over the Fluent in-built model. Deposition in the nasal cavity was high for very small nanoparticles. The particle diameter range in which the deposition drops from 80 to 18% is between 1 and 10 nm. From 10 to 150 nm, however, there is only a small change in the deposition curve from 18 to 15%. A similar deposition curve profile was found for the TB airway.

Unsteady Particle Deposition in a Human Nasal Cavity during Inhalation

The present study investigates the deposition efficiency during the unsteady inhalation cycle by using Computational Fluid Dynamics (CFD). The unsteady inhalation profile was applied at the outlet of nasopharynx, which had a maximum flow rate of 40.3L/min which corresponds to an equivalent steady inhalation tidal volume flow rate of 24.6L/min. Aerodynamic particle sizes of 5μm and 20μm were studied in order to reflect contrasting Stokes numbered particle behaviour. Two particle deposition efficiencies in the nasal cavity versus time are presented. In general, the deposition of 5μm particles was much less than 20μm particles. The first 0.2 second of the inhalation cycle was found to be significant to the particle transport, since the majority of particles were deposited during this period (i.e. its residence time). Comparisons were also made with its equivalent steady inhalation flow rate which found that the unsteady inhalation produced lower deposition efficiency for both particle sizes.

Inhalability of micron particles through the nose and mouth

Aspiration efficiencies from nose and mouth inhalations are investigated at low and high inhalation rates by using the commercial Computational Fluid Dynamics (CFD) software CFX 11. A realistic human head with detailed facial features was constructed. Facial features were matched to represent the 50th percentile of a human male, aged between 20 and 65 years old, based on anthropometric data. The constant freestream velocity was 0.2 ms(-1), normal to the face, and inhalation rates through the mouth and nose were 15 liters per minute (LPM) for light breathing and 40 LPM for heavy breathing. It was found that the flow field in the near breathing region exhibited vertical direction caused by the presence of the torso where the airstream diverges as it flows around and over the body. The critical area concept was used as a tool to determine the aspiration efficiency of particles. Comparisons between critical areas for the nose and mouth inhalations show similar geometric properties such as the area's shape, and its vertical distance location on the x-z plane located at y = 80 cm upstream. The critical area sizes were found to be slightly larger for the mouth inhalation mainly due to the larger mouth area and also the aligned orientation of the mouth to the upstream flow, whereas the nose is perpendicular to the upstream flow. This study was undertaken to establish the flow field in the near breathing region that will help to characterize the flow and particle field for initial boundary conditions leading to a more holistic modeling approach of respiration through the internal nasal cavity and mouth.

Comparison of micron- and nanoparticle deposition patterns in a realistic human nasal cavity

Knowledge regarding particle deposition processes in the nasal cavity is important in aerosol therapy and inhalation toxicology applications. This paper presents a comparative study of the deposition of micron and submicron particles under different steady laminar flow rates using a Lagrangian approach. A computational model of a nasal cavity geometry was developed from CT scans and the simulation of the fluid and particle flow within the airway was performed using the commercial software GAMBIT and FLUENT. The air flow patterns in the nasal cavities and the detailed local deposition patterns of micron and submicron particles were presented and discussed. It was found that the majority of micron particles are deposited near the nasal valve region and some micron particles are deposited on the septum wall in the turbinate region. The deposition patterns of micron particles in the left cavity are different compared with that in the right one especially in the turbinate regions. In contrast, the deposition for nanoparticles shows a moderately even distribution of particles throughout the airway. Furthermore the particles releasing position obviously influences the local deposition patterns. The influence of the particle releasing position is mainly shown near the nasal valve region for micron particle deposition, while for submicron particles deposition, both the nasal valve and turbinate region are influenced. The results of the paper are valuable in aerosol therapy and inhalation toxicology.

Computational Modelling of Gas-Particle Flows with Different Particle Morphology in the Human Nasal Cavity

This paper summarises current studies related to numerical gas-particle flows in the human nasal cavity. Of interest are the numerical modelling requirements to consider the effects of particle morphology for a variety of particle shapes and sizes such as very small particles sizes (nanoparticles), elongated shapes (asbestos fibres), rough shapes (pollen), and porous light density particles (drug particles) are considered. It was shown that important physical phenomena needed to be addressed for different particle characteristics. This included the Brownian diffusion for submicron particles. Computational results for the nasal capture efficiency for nano-particles and various breathing rates in the laminar regime were found to correlate well with the ratio of particle diffusivity to the breathing rate. For micron particles, particle inertia is the most significant property and the need to use sufficient drag laws is important. Drag correlations for fibrous and rough surfaced particles were investigated to enable particle tracking. Based on the simulated results, semi-empirical correlations for particle deposition were fitted in terms of Peclet number and inertial parameter for nanoparticles and micron particles respectively.

Optimising nasal spray parameters for efficient drug delivery using computational fluid dynamics

Experimental images from particle/droplet image analyser (PDIA) and particle image velocimetry (PIV) imaging techniques of particle formation from a nasal spray device were taken to determine critical parameters for the study and design of effective nasal drug delivery devices. The critical parameters found were particle size, diameter of spray cone at a break-up length and a spray cone angle. A range of values for each of the parameters were ascertained through imaging analysis which were then transposed into initial particle boundary conditions for particle flow simulation within the nasal cavity by using Computational Fluid Dynamics software. An Eulerian-Lagrangian scheme was utilised to track mono-dispersed particles (10 and 20 microm) at a breathing rate of 10 L/min. The results from this qualitative study aim to assist the pharmaceutical industry to improve and help guide the design of nasal spray devices.

Numerical simulations for detailed airflow dynamics in a human nasal cavity

Nasal physiology is dependent on the physical structure of the nose. Individual aspects of the nasal cavity such as the geometry and flow rate collectively affect nasal function such as the filtration of foreign particles by bringing inspired air into contact with mucous-coated walls, humidifying and warming the air before it enters the lungs and the sense of smell. To better understand the physiology of the nose, this study makes use of CFD methods and post-processing techniques to present flow patterns between the left and right nasal cavities and compared the results with experimental and numerical data that are available in literature. The CFD simulation adopted a laminar steady flow for flow rates of 7.5 L/min and 15 L/min. General agreement of gross flow features were found that included high velocities in the constrictive nasal valve area region, high flow close to the septum walls, and vortex formations posterior to the nasal valve and olfactory regions. The differences in the left and right cavities were explored and the effects it had on the flow field were discussed especially in the nasal valve and middle turbinate regions. Geometrical differences were also compared with available models.

Deposition of inhaled wood dust in the nasal cavity

Detailed deposition patterns of inhaled wood dust in an anatomically accurate nasal cavity were investigated using computational fluid dynamics (CFD) techniques. Three wood dusts, pine dust, heavy oak dust, and light oak dust, with a particle size distribution generated by machining (Chung et al., 2000), were simulated at an inhalation flow rate of 10 L/min. It was found that the major particle deposition sites were the nasal valve region and anterior section of the middle turbinate. Wood dust depositing in these regions is physiologically removed much more slowly than in other regions. This leads to the surrounding layer of soft tissues being damaged by the deposited particles during continuous exposure to wood dust. Additionally, it was found that pine dust had a higher deposition efficiency in the nasal cavity than the two oak dusts, due to the fact that it comprises a higher proportion of larger sized particles. Therefore, this indicates that dusts with a large amount of fine particles, such as those generated by sanding, may penetrate the nasal cavity and travel further into the lung.