Inhalability of micron particles through the nose and mouth

Numerical study of transient indoor airflow and virus-laden droplet dispersion: Impact of interactive human movement

Human movement affects indoor airflow and the airborne transmission of respiratory infectious diseases, which has attracted scholars. However, the interactive airflow between moving and stationary people has yet to be studied in detail. This study used the numerical method validated by experimental data to explore the transient indoor airflow and virus-laden droplet dispersion in scenes with interactive human movement. Human-shaped numerical models and the dynamic mesh method were adopted to realize human movement in scenes with different lateral distances (0.2-1.2 m) between a moving person and stationary (standing/sitting) persons. The interactive human movement obviously impacts other persons' respiratory airflow, and the lateral fusion ranged about 0.6 m. The interactive human movement strengthens the indoor airflow convection, and some exhaled virus-laden droplets were carried into wake flow and enhanced long-range airborne transmission. The impact of interactive human movement on sitting patients' exhalation airflow seems more evident than on standing patients. The impact might last over 2 min after movement stopped, and people in the affected area might be at a higher exposure. The results can provide a reference for epidemic control in indoor environments.

Detailed comparison of anatomy and airflow dynamics in human and cynomolgus monkey nasal cavity

Nonhuman primates are occasionally used as laboratory models for sophisticated medical research as they bear the closest resemblance to humans in morphometry and physiological functions. A range of nonhuman primate species have been employed in the inhalation toxicity, nasal drug delivery and respiratory viral infection studies, and they provided valuable insight to disease pathogenesis while other laboratory animals such as rodents cannot recapitulate due to the lesser degree of similarity in metabolism, anatomy and cellular response to that of humans. It is anticipated that nonhuman primate models of respiratory diseases will continue to be instrumental for translating biomedical research for improvement of human health, and the confidence in laboratory data extrapolation between species will play a pivotal role. From the morphometry and flow dynamics point of view, this study performed a detailed comparative analysis between human and a cynomolgus monkey nasal airway, with intention to provide high-fidelity qualitative and quantitative linkage between the two species for more effective laboratory data extrapolation. The study revealed that cynomolgus monkey could be a good human surrogate in nasal inhalation studies; however, care should be given for interspecies data extrapolation as subtle differences in anatomy and airflow dynamics were present between the two species.

Particle inhalability of a standing mannequin with large airways in a ventilated room

This study presents a series of numerical simulations for airflow field and particle dispersion and deposition around a mannequin inside a ventilated room. A 3-D airway system of a volunteer subject with a large respiratory system was reconstructed from the nostril inlet to the end of the tracheobronchial tree 4th generation and was integrated into a standing mannequin at the center of a room. The room ventilation system supplied air through a diffuser and expelled air via a damper in three modes. The airflow field was first evaluated by solving the governing equations and the k-ω SST transitional turbulence model using the Ansys-Fluent software. Then spherical particles with various diameters were released into the room, and their trajectories were evaluated using the Lagrangian approach. Aspiration fraction and particle deposition for inhalation flow rates of 15 and 30 L/min were analyzed using a modified discrete random walk (DRW) stochastic model using a user-defined function (UDF) coupled to the Ansys-Fluent discrete phase model. For the first ventilation mode, a recirculation flow region formed behind the mannequin that led the airflow streamlines to the breathing zone. A recirculation flow formed in front of the face for the second ventilation mode that led the airflow streamlines out of the mannequin breathing zone. For the third mode, however, there was no strong recirculation flow zone around the mannequin. Simulation results showed that the aspiration fraction in the first ventilation mode was higher than the other modes. In addition, the regional deposition rates and deposition patterns of particles inside the respiratory system were presented for each region. Accordingly, most large particles were trapped in the nasal passage; however, some large particles penetrated deeper into the airway due to the large airway size. For the higher breathing rate, the percentage of large escaped particles from the lobe branches dropped by a factor of 7 compared to the lower breathing rate.

New Approach Methodology for Assessing Inhalation Risks of a Contact Respiratory Cytotoxicant: Computational Fluid Dynamics-Based Aerosol Dosimetry Modeling for Cross-Species and In Vitro Comparisons

Regulatory agencies are considering alternative approaches to assessing inhalation toxicity that utilizes in vitro studies with human cells and in silico modeling in lieu of additional animal studies. In support of this goal, computational fluid-particle dynamics models were developed to estimate site-specific deposition of inhaled aerosols containing the fungicide, chlorothalonil, in the rat and human for comparisons to prior rat inhalation studies and new human in vitro studies. Under bioassay conditions, the deposition was predicted to be greatest at the front of the rat nose followed by the anterior transitional epithelium and larynx corresponding to regions most sensitive to local contact irritation and cytotoxicity. For humans, simulations of aerosol deposition covering potential occupational or residential exposures (1-50 µm diameter) were conducted using nasal and oral breathing. Aerosols in the 1-5 µm range readily penetrated the deep region of the human lung following both oral and nasal breathing. Under actual use conditions (aerosol formulations >10 µm), the majority of deposited doses were in the upper conducting airways. Beyond the nose or mouth, the greatest deposition in the pharynx, larynx, trachea, and bronchi was predicted for aerosols in the 10-20 µm size range. Only small amounts of aerosols >20 µm penetrated past the pharyngeal region. Using the ICRP clearance model, local retained tissue dose metrics including maximal concentrations and areas under the curve were calculated for each airway region following repeated occupational exposures. These results are directly comparable with benchmark doses from in vitro toxicity studies in human cells leading to estimated human equivalent concentrations that reduce the reliance on animals for risk assessments.

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.

Identification of refractory zirconia from catalytic converters in dust: An emerging pollutant in urban environments

Using catalytic converters is one of the most effective methods to control vehicle emissions. A washcoat of cerium oxide-zirconia (CeO₂-ZrO₂) has been used to enhance the performance of the catalytic converter device. To date, the prevalence of this material in the environment has not been assessed. In this study, we present evidence of the existence of inhalable zirconia in urban dust. Samples of the washcoat, exhaust pipe, topsoil, and road dust were analyzed by X-ray fluorescence, X-ray diffraction (XRD), scanning electron microscopy (SEM), Raman spectroscopy, photoluminescence (PL) spectroscopy, and thermally stimulated luminescence (TSL). The results showed a CeO₂-ZrO₂ phase separation after sintering. This causes the emission of ZrO₂, CeO₂, and CeZrO_x particles smaller than 1 μm, which can likely reach the alveolar macrophages in the lungs. The Ce-Zr content in road dust exceeds geogenic levels, and a significant correlation of 0.87 (p < 0.05) reflects a common anthropic source. Chronic exposure to such refractory particles may result in the development of non-occupational respiratory diseases. The inhalable crystalline compounds emitted by vehicles are a significant environmental health hazard, revealing the need for further investigation and assessment of zirconia levels generated by automobiles in urban areas worldwide.

Regional deposition of the allergens and micro-aerosols in the healthy human nasal airways

The nasal cavity is the inlet to the human respiratory system and is responsible for the olfactory sensation, filtering pollutant particulate matter, and humidifying the air. Many research studies have been performed to numerically predict allergens, contaminants, and/or drug particle deposition in the human nasal cavity; however, the majority of these investigations studied only one or a small number of nasal passages. In the present study, a series of Computed Tomography (CT) scan images of the nasal cavities from ten healthy subjects were collected and used to reconstruct accurate 3D models. All models were divided into twelve anatomical regions in order to study the transport and deposition features of different regions of the nasal cavity with specific functions. The flow field and micro-particle transport equations were solved, and the total and regional particle deposition fractions were evaluated for the rest and low activity breathing conditions. The results show that there are large variations among different subjects. The standard deviation of the total deposition fraction in the nasal cavities was the highest for 5 × 10 4 <impaction parameter (IP)< 1.125 × 10 5 with a maximum of 20%. The achieved results highlighted the nasal cavity sections that are more involved in the particle deposition. Particles with IP = 30,000 deposit more in the middle turbinate and nasopharynx areas, while for particles with IP = 300,000, deposition is mainly in the anterior parts (kiesselbach and vestibule regions). For small IP values, the amounts of deposition fractions in different regions of the nasal cavity are more uniform.

Host-to-host airborne transmission as a multiphase flow problem for science-based social distance guidelines

The COVID-19 pandemic has strikingly demonstrated how important it is to develop fundamental knowledge related to the generation, transport and inhalation of pathogen-laden droplets and their subsequent possible fate as airborne particles, or aerosols, in the context of human to human transmission. It is also increasingly clear that airborne transmission is an important contributor to rapid spreading of the disease. In this paper, we discuss the processes of droplet generation by exhalation, their potential transformation into airborne particles by evaporation, transport over long distances by the exhaled puff and by ambient air turbulence, and their final inhalation by the receiving host as interconnected multiphase flow processes. A simple model for the time evolution of droplet/aerosol concentration is presented based on a theoretical analysis of the relevant physical processes. The modeling framework along with detailed experiments and simulations can be used to study a wide variety of scenarios involving breathing, talking, coughing and sneezing and in a number of environmental conditions, as humid or dry atmosphere, confined or open environment. Although a number of questions remain open on the physics of evaporation and coupling with persistence of the virus, it is clear that with a more reliable understanding of the underlying flow physics of virus transmission one can set the foundation for an improved methodology in designing case-specific social distancing and infection control guidelines.

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.

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.

Geometry and airflow dynamics analysis in the nasal cavity during inhalation

BACKGROUND

A major issue among computational respiratory studies is the wide variety of nasal morphologies being studied, caused by both inter-population and inter-subject variations.

METHOD

Six nasal cavity geometries exhibiting diverse geometry variations were subjected to steady inhalation flow rate of 15L/min. to determine if any consistent flow behaviour could be found.

FINDINGS

Despite vastly different geometries we were able to identify consistent flow patterns including relatively high velocity in the nasal valve region, followed by flow continuing predominantly in the inferior half of the airway. We also found conformity among models where the inhaled air reached a near-conditioned state by the middle of the nasal cavity. Air from the front of the face reached the olfactory regions while air from the lateral sides of the face moved through the inferior half of the nasal cavity.

INTERPRETATION

The ability to predict gross flow features provides a baseline flow field to compare against. This contributes towards establishing well defined flow predictions and be used as a comparison for future larger studies.

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.

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.

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.

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.

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

BACKGROUND

There are methodological ambiguities in the literature on mesh refinement analysis for computational fluid dynamics (CFD) modeling of physiologically realistic airflow dynamics and particle transport in the human sinonasal cavity. To investigate grid independence in discretization of the (sino)nasal geometry, researchers have considered CFD variables such as pressure drop, velocity profile, wall shear, airflow, and particle deposition fractions. Standardization in nasal geometry is also lacking: unilateral or bilateral nasal cavities with and without paranasal sinuses have been used. These methodological variants have led to inconsistencies in establishing grid-independent mesh densities. The aim of this study is to provide important insight in the role of mesh refinement analysis on airflow and particle deposition in sinonasal airway modeling.

METHODS

A three-dimensional reconstruction of the complete sinonasal cavity was created from computed tomography images of a subject who had functional endoscopic sinus surgery. To investigate airflow grid independence, nine different tetrahedral mesh densities were generated. For particle transport mesh refinement analysis, hybrid tetrahedral-prism elements with near-wall prisms ranging from 1 to 6 layers were implemented. Steady-state, laminar inspiratory airflow simulations under physiologic pressure-driven conditions and nebulized particle transport simulations were performed with particle sizes ranging from 1-20 μm.

RESULTS

Mesh independence for sinonasal airflow was achieved with approximately 4 million unstructured tetrahedral elements. The hybrid mesh containing 4 million tetrahedral cells with three prism layers demonstrated asymptotic behavior for sinonasal particle deposition. Inclusion of boundary prism layers reduced deposition fractions relative to tetrahedral-only meshes.

CONCLUSIONS

To ensure numerically accurate simulation results, mesh refinement analyses should be performed for both airflow and particle transport simulations. Tetrahedral-only meshes overpredict particle deposition and are less accurate than hybrid tetrahedral-prism meshes.

Effects of the facial interface on inhalation and deposition of micrometer particles in calm air in a child airway model

CONTEXT

How the facial interface affects particle inhalability and depositions within the airway is not well understood. Previous studies of inhalation dosimetry are limited to either inhalability or deposition, rather than the two studied in a systematic way.

OBJECTIVE

To systematically evaluate the effects of the facial interface on aerosol inhalability, nasal deposition and thoracic dose in a 5-year-old child airway model using a coupled imaging-computational fluid dynamics approach.

METHODS

A face-nose-throat model was developed from magnetic resonance imaging scans of a 5-year-old boy. Respiration airflows and particle transport were simulated with the low Reynolds number k-ω turbulence model and the Lagrangian tracking approach. Particles ranging from 1 to 70 µm were considered in a calm air.

RESULTS

Retaining the facial interface in the computational model induced substantial variations in flow dynamics, aerosol inhalability and thoracic doses. The nasal and thoracic deposition fractions were much lower with the facial interface due to the low inhalability into downward-facing nostrils and facial deposition losses. For a given inhalation rate of 10 L/min, including the facial interface reduced the thoracic dose by 5% for 2.5-µm particles and by 50% for 10 µm particles in the child model. Considering localized conditions, facial interface substantially increased depositions at the turbinate region and dorsal pharynx.

CONCLUSION

This study highlighted the need to include facial interface in future numerical and in vitro studies. Findings of this study have practical implications in the design of aerosol samplers and interpretation of deposition data from studies without facial interfaces.

An empirical model of human aspiration in low-velocity air using CFD investigations

Computational fluid dynamics (CFD) modeling was performed to investigate the aspiration efficiency of the human head in low velocities to examine whether the current inhaled particulate mass (IPM) sampling criterion matches the aspiration efficiency of an inhaling human in airflows common to worker exposures. Data from both mouth and nose inhalation, averaged to assess omnidirectional aspiration efficiencies, were compiled and used to generate a unifying model to relate particle size to aspiration efficiency of the human head. Multiple linear regression was used to generate an empirical model to estimate human aspiration efficiency and included particle size as well as breathing and freestream velocities as dependent variables. A new set of simulated mouth and nose breathing aspiration efficiencies was generated and used to test the fit of empirical models. Further, empirical relationships between test conditions and CFD estimates of aspiration were compared to experimental data from mannequin studies, including both calm-air and ultra-low velocity experiments. While a linear relationship between particle size and aspiration is reported in calm air studies, the CFD simulations identified a more reasonable fit using the square of particle aerodynamic diameter, which better addressed the shape of the efficiency curve's decline toward zero for large particles. The ultimate goal of this work was to develop an empirical model that incorporates real-world variations in critical factors associated with particle aspiration to inform low-velocity modifications to the inhalable particle sampling criterion.

Influence of secondary aspiration on human aspiration efficiency

Computational fluid dynamics (CFD) was used to evaluate the contribution of secondary aspiration to human aspiration efficiency estimates using a humanoid model with realistic facial features. This study applied coefficient of restitution (CoR) values for working-aged human facial skin to the facial regions on the humanoid CFD model. Aspiration efficiencies for particles ranging from 7 to 116 μm were estimated for bounce (allowing for secondary aspiration) and no-bounce (CoR=0) simulations. Fluid simulations used the standard k-epsilon turbulence model over a range of test conditions: three freestream velocities, two breathing modes (mouth and nose breathing, using constant inhalation), three breathing velocities, and five orientations relative to the oncoming wind. Laminar particle trajectory simulations were used to examine inhaled particle transport and estimate aspiration efficiencies. Aspiration efficiency for the realistic CoR simulations, for both mouth- and nose-breathing, decreased with increasing particle size, with aspiration around 50% for 116 μm particles. For the CoR=0 simulations, aspiration decreased more rapidly with increasing particle size and approached zero for 116 μm compared to realistic CoR models (differences ranged from 0% to 80% over the particle sizes and velocity conditions). Differences in aspiration efficiency were larger with increasing particle size (>52 μm) and increased with decreasing freestream velocity and decreasing breathing rate. Secondary aspiration was more important when the humanoid faced the wind, but these contributions to overall aspiration estimates decreased as the humanoid rotated through 90°. There were minimal differences in aspiration between uniform CoR values of 0.5, 0.8, 1.0 and realistic regionally-applied CoR values, indicating differences between mannequin surfaces and between mannequin and human skin will have negligible effect on aspiration for facing-the-wind orientation.

Numerical investigation of regional particle deposition in the upper airway of a standing male mannequin in calm air surroundings

A 3-D realistic computational model of the airway system integrated into a standing male mannequin was developed. The computational domain includes the regions around the mannequin and the inside of the airway passages. The simulation was performed for low activity breathing rates with calm air around the mannequin. The flowfield of the inhaled air was first obtained from solving the Navier-Stokes and continuity equations. Then the particles were released in the domain around the mannequin and their trajectories were evaluated by using the Lagrangian approach for solving the particle equation of motion. The regional aerosols deposition was evaluated for different parts of the human airway system and the results were compared with those obtained from the separate modeling of the airway system without the interaction of the airflow with the mannequin external face. The results showed when the upper airway is integrated into the mannequin, the regional deposition of inhaled particles mainly changes in the airway system.

Computational fluid dynamics investigation of human aspiration in low velocity air: orientation effects on nose-breathing simulations

An understanding of how particles are inhaled into the human nose is important for developing samplers that measure biologically relevant estimates of exposure in the workplace. While previous computational mouth-breathing investigations of particle aspiration have been conducted in slow moving air, nose breathing still required exploration. Computational fluid dynamics was used to estimate nasal aspiration efficiency for an inhaling humanoid form in low velocity wind speeds (0.1-0.4 m s(-1)). Breathing was simplified as continuous inhalation through the nose. Fluid flow and particle trajectories were simulated over seven discrete orientations relative to the oncoming wind (0, 15, 30, 60, 90, 135, 180°). Sensitivities of the model simplification and methods were assessed, particularly the placement of the recessed nostril surface and the size of the nose. Simulations identified higher aspiration (13% on average) when compared to published experimental wind tunnel data. Significant differences in aspiration were identified between nose geometry, with the smaller nose aspirating an average of 8.6% more than the larger nose. Differences in fluid flow solution methods accounted for 2% average differences, on the order of methodological uncertainty. Similar trends to mouth-breathing simulations were observed including increasing aspiration efficiency with decreasing freestream velocity and decreasing aspiration with increasing rotation away from the oncoming wind. These models indicate nasal aspiration in slow moving air occurs only for particles <100 µm.

Computational fluid dynamics investigation of human aspiration in low-velocity air: orientation effects on mouth-breathing simulations

Computational fluid dynamics was used to investigate particle aspiration efficiency in low-moving air typical of occupational settings (0.1-0.4 m s(-1)). Fluid flow surrounding an inhaling humanoid form and particle trajectories traveling into the mouth were simulated for seven discrete orientations relative to the oncoming wind (0°, 15°, 30°, 60°, 90°, 135° and 180°). Three continuous inhalation velocities (1.81, 4.33, and 12.11 m s(-1)), representing the mean inhalation velocity associated with sinusoidal at-rest, moderate, and heavy breathing (7.5, 20.8, and 50.3 l min(-1), respectively) were simulated. These simulations identified a decrease in aspiration efficiency below the inhalable particulate mass (IPM) criterion of 0.5 for large particles, with no aspiration of particles 100 µm and larger for at-rest breathing and no aspiration of particles 116 µm for moderate breathing, over all freestream velocities and orientations relative to the wind. For particles smaller than 100 µm, orientation-averaged aspiration efficiency exceeded the IPM criterion, with increased aspiration efficiency as freestream velocity decreased. Variability in aspiration efficiencies between velocities was low for small (<22 µm) particles, but increased with increasing particle size over the range of conditions studied. Orientation-averaged simulation estimates of aspiration efficiency agree with the linear form of the proposed linear low-velocity inhalable convention through 100 µm, based on laboratory studies using human mannequins.

Uncertainty in aspiration efficiency estimates from torso simplifications in computational fluid dynamics simulations

Computational fluid dynamics (CFD) has been used to report particle inhalability in low velocity freestreams, where realistic faces but simplified, truncated, and cylindrical human torsos were used. When compared to wind tunnel velocity studies, the truncated models were found to underestimate the air's upward velocity near the humans, raising questions about aspiration estimation. This work compares aspiration efficiencies for particles ranging from 7 to 116 µm using three torso geometries: (i) a simplified truncated cylinder, (ii) a non-truncated cylinder, and (iii) an anthropometrically realistic humanoid body. The primary aim of this work is to (i) quantify the errors introduced by using a simplified geometry and (ii) determine the required level of detail to adequately represent a human form in CFD studies of aspiration efficiency. Fluid simulations used the standard k-epsilon turbulence models, with freestream velocities at 0.1, 0.2, and 0.4 m s(-1) and breathing velocities at 1.81 and 12.11 m s(-1) to represent at-rest and heavy breathing rates, respectively. Laminar particle trajectory simulations were used to determine the upstream area, also known as the critical area, where particles would be inhaled. These areas were used to compute aspiration efficiencies for facing the wind. Significant differences were found in both vertical velocity estimates and the location of the critical area between the three models. However, differences in aspiration efficiencies between the three forms were <8.8% over all particle sizes, indicating that there is little difference in aspiration efficiency between torso models.