Computational Fluid Dynamics Simulations of Particle Deposition in Large-Scale, Multigenerational Lung Models

2010 ◽  
Vol 133 (1) ◽  
Author(s):  
D. Keith Walters ◽  
William H. Luke

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.

Author(s):  
D. Keith Walters ◽  
William H. Luke

Computational fluid dynamics (CFD) has evolved as a useful tool for the prediction of airflow and particle transport within the human lung airway. A large number of published studies have demonstrated the use of CFD coupled with Lagrangian particle tracking methods to determine local and regional deposition rates in small subsections of the bronchopulmonary tree. However, simulation of particle transport and deposition in large-scale models encompassing more than a few generations is less common, due primarily to the sheer size and complexity of the human lung airway geometry. 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 [1], which was shown to produce results similar to fully resolved geometries using approximate, reduced geometry models. The methodology is here extended to particle transport and deposition simulations. Lagrangian particle-tracking simulations are performed in combination with Eulerian simulations of the air flow in an idealized representation of the human lung airway tree. Results using the reduced models are compared to fully resolved models for an eight-generation region of the conducting zone. Agreement between fully resolved and reduced geometry simulations indicates that the new method can provide an accurate alternative for large-scale CFD simulations while reducing the computational cost of these simulations by an order of magnitude or more.


Author(s):  
Kevin T. Shanley ◽  
Goodarz Ahmadi ◽  
Philip K. Hopke ◽  
Yung-Sung Cheng

As the interface between the human respiratory system and the environment, the nose plays many vital roles. Not the least of which is filter. Resulting from numerous natural and anthropogenic processes, particulate matter becomes airborne. Should particulate matter reach the lower portions of the respiratory tract, a host of maladies may occur. In an attempt to further understand the physics behind particulate matter transitioning from the environment into humans a computational model has been developed to predict the efficiency with which human noses can remove particles before they reach the lungs. To this end computational fluid dynamics and Lagrangian particle tracking simulations have been run to gather information on the deposition behavior of both fibrous and spherical particles. MRI data was collected from the left and right passages of a 181.6 cm, 120.2 kg, human male. The two passages were constructed into separate computational volumes consisting of approximately 950,000 unstructured tetrahedral cells each. A steady laminar flow model was used to simulate the inhalation portion of a human breathing cycle. Volumetric flow rates were varied to represent the full range of human nasal breathing. General agreement was shared quantitatively and qualitatively with previously published in vitro studies on other nasal models. Lagrangian particle tracking was performed for varying sizes of fibrous and spherical particles. Deposition efficiency was shown to increase with fiber aspect ratio, particle size, and flow rate. Anatomy was also identified as effecting deposition.


Author(s):  
Johannes Bosbach ◽  
Daniel Schanz ◽  
Phillip Godbersen ◽  
Andreas Schröder

We present spatially and temporally resolved velocity and acceleration measurements of turbulent RayleighBénard convection spanning the whole volume (~ 1 m³) of a cylindrical sample with aspect ratio one. With the "Shake-The-Box" (STB) Lagrangian particle tracking (LPT) algorithm, we were able to instantaneously track up to 560,000 particles, corresponding to mean inter-particle distances down to 6 - 8 Kolmogorov lengths. We used the data assimilation scheme ‘FlowFit’, which involves continuity and Navier-Stokesconstraints, to map the scattered velocity and acceleration data on cubic grids, herewith recovering the smallest flow scales. Lagrangian and Eulerian visualizations reveal the dynamics of the large-scale circulation and its interplay with small scale structures, such as thermal plumes and turbulent background fluctuations. As a result, the complex time-dependent behavior of the LSC comprising azimuthal rotations, torsional oscillation and sloshing can be extracted from the data. Further, we found more seldom dynamic events, such as spontaneous reorientations of the LSC in the data from long-term measurements.


Author(s):  
Andreas Schröder ◽  
Daniel Schanz ◽  
Johannes Bosbach ◽  
Matteo Novara ◽  
Reinhard Geisler ◽  
...  

Exhalation of small aerosol droplets and their transport, dispersion and (local) accumulation in closed rooms have been identified as the main pathway for indirect or airborne respiratory virus transmission from person to person, e.g. for SARS-CoV 2 or measles (Morawska and Cao 2020). Understanding airborne transport mechanisms of viruses via small bio-aerosol particles inside closed populated rooms is an important key factor for optimizing various mitigation strategies (Morawska et al. 2020), which can play an important role for damping the infection dynamics of any future and the ongoing present pandemic scenario, which unfortunately, is still threatening due to the spreading of several SARS-CoV2 variants of concern, e.g. delta (Kupferschmidt and Wadman 2021). Therefore, a large-scale 3D Lagrangian Particle Tracking experiment using up to 3 million long lived and nearly neutrally buoyant helium-filled soap bubbles (HFSB) with a mean diameter of ~ 370 µm as passive tracers in a 12 m³ generic test room has been performed, which allows to fully resolve the Lagrangian transport properties and flow field inside the whole room around a cyclically breathing thermal manikin (Lange et al. 2012) with and without mouth-nose-masks and shields applied. Six high-resolution CMOS streaming cameras, a large array of powerful pulsed LEDs have been used and the Shake-The-Box (STB) (Schanz et al. 2016) Lagrangian particle tracking algorithm has been applied in this experimental study of internal flows in order to gain insight into the complex transient and turbulent aerosol particle transport and dispersion processes around seated breathing persons.


2021 ◽  
Vol 67 (5) ◽  
pp. 1497-1513
Author(s):  
Jana Wedel ◽  
Paul Steinmann ◽  
Mitja Štrakl ◽  
Matjaž Hriberšek ◽  
Jure Ravnik

AbstractTo respond to the ongoing pandemic of SARS-CoV-2, this contribution deals with recently highlighted COVID-19 transmission through respiratory droplets in form of aerosols. Unlike other recent studies that focused on airborne transmission routes, this work addresses aerosol transport and deposition in a human respiratory tract. The contribution therefore conducts a computational study of aerosol deposition in digital replicas of human airways, which include the oral cavity, larynx and tracheobronchial airways down to the 12th generation of branching. Breathing through the oral cavity allows the air with aerosols to directly impact the larynx and tracheobronchial airways and can be viewed as one of the worst cases in terms of inhalation rate and aerosol load. The implemented computational model is based on Lagrangian particle tracking in Reynolds-Averaged Navier–Stokes resolved turbulent flow. Within this framework, the effects of different flow rates, particle diameters and lung sizes are investigated to enable new insights into local particle deposition behavior and therefore virus loads among selected age groups. We identify a signicant increase of aerosol deposition in the upper airways and thus a strong reduction of virus load in the lower airways for younger individuals. Based on our findings, we propose a possible relation between the younger age related fluid mechanical protection of the lower lung regions due to the airway size and a reduced risk of developing a severe respiratory illness originating from COVID-19 airborne transmission.


2016 ◽  
Vol 46 (10) ◽  
pp. 2995-3010 ◽  
Author(s):  
Edward W. Doddridge ◽  
David P. Marshall ◽  
Andrew McC. Hogg

AbstractThe presence of large-scale Ekman pumping associated with the climatological wind stress curl is the textbook explanation for low biological activity in the subtropical gyres. Using an idealized, eddy-resolving model, it is shown that Eulerian-mean Ekman pumping may be opposed by an eddy-driven circulation, analogous to the way in which the atmospheric Ferrel cell and the Southern Ocean Deacon cell are opposed by eddy-driven circulations. Lagrangian particle tracking, potential vorticity fluxes, and depth–density streamfunctions are used to show that, in the model, the rectified effect of eddies acts to largely cancel the Eulerian-mean Ekman downwelling. To distinguish this effect from eddy compensation, it is proposed that the suppression of Eulerian-mean downwelling by eddies be called “eddy cancellation.”


1996 ◽  
Vol 33 (9) ◽  
pp. 163-170 ◽  
Author(s):  
Virginia R. Stovin ◽  
Adrian J. Saul

Research was undertaken in order to identify possible methodologies for the prediction of sedimentation in storage chambers based on computational fluid dynamics (CFD). The Fluent CFD software was used to establish a numerical model of the flow field, on which further analysis was undertaken. Sedimentation was estimated from the simulated flow fields by two different methods. The first approach used the simulation to predict the bed shear stress distribution, with deposition being assumed for areas where the bed shear stress fell below a critical value (τcd). The value of τcd had previously been determined in the laboratory. Efficiency was then calculated as a function of the proportion of the chamber bed for which deposition had been predicted. The second method used the particle tracking facility in Fluent and efficiency was calculated from the proportion of particles that remained within the chamber. The results from the two techniques for efficiency are compared to data collected in a laboratory chamber. Three further simulations were then undertaken in order to investigate the influence of length to breadth ratio on chamber performance. The methodology presented here could be applied to complex geometries and full scale installations.


Sign in / Sign up

Export Citation Format

Share Document