Prediction of aerosol deposition in 90∘ bends using LES and an efficient Lagrangian tracking method

Abstract Aiming at the better prediction of pharmaceutical aerosol deposition in extrathoracic airways, a simpler test case, namely a 90 ∘ bend flow (tubular cross-section) laden with monodisperse particles, is adopted here and studied numerically. The continuous phase is calculated using a large-eddy simulation technique along with a finite-volume method for block-structured curvilinear grids. The particulate phase is simulated using a Lagrangian approach where hundred thousands of individual monodisperse particles with varying particle diameters are released and tracked throughout the computational domain. To allow such a large number of particles, a highly efficient tracking algorithm is applied, where particle paths are predicted in an orthogonal computational domain, avoiding time-consuming search algorithm, normally required when particles are tracked in the actual physical domain of a curvilinear body-fitted block-structured grid. Both simulation algorithms, for the continuous and particulate phases, are completely parallelized using domain decomposition. Additionally, the in-house code applied supports vector processing allowing efficient usage of nearly all kinds of high-performance architectures. Two different Reynolds numbers Re D are considered where Re D is based on the bend diameter and mean flow velocity. The first case is within the laminar regime at Re D = 1000 and serves for the purpose of verification and validation. The second, more challenging case comprises the turbulent regime at Re D = 10 , 000 , which is the intrinsic objective of the present study. Depending on the Stokes number of the particles, 0.001 ⩽ St ⩽ 1.5 , and the releasing locations at the entrance of the bend, the particles will either deposit on the wall or penetrate and exit the computational domain. Simulation results of aerosol deposition efficiency, over the entire range of particle diameters considered here, show an excellent agreement when compared to experimental values obtained by Pui, Romay-Novas, and Liu [(1987). Experimental study of particle deposition in bends of circular cross-section. Aerosol Science and Technology, 7, 301].

[1]  Michael Breuer,et al.  Comparison of c‐space and p‐space particle tracing schemes on high‐performance computers: accuracy and performance , 2002 .

[2]  J. McLaughlin,et al.  On the role of the lift force in turbulence simulations of particle deposition , 1997 .

[3]  Andrew R. McFarland,et al.  Aerosol deposition in bends with turbulent flow , 1997 .

[4]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[5]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[6]  K. Lilly The representation of small-scale turbulence in numerical simulation experiments , 1966 .

[7]  Benjamin Y. H. Liu,et al.  Experimental Study of Particle Deposition in Bends of Circular Cross Section , 1987 .

[8]  G. Rudolf,et al.  Intercomparison of Experimental Regional Aerosol Deposition Data , 1989 .

[9]  Ernst Heinrich Hirschel,et al.  Flow Simulation with High-Performance Computers II , 1996 .

[10]  J. McLaughlin,et al.  A New Correlation for the Aerosol Deposition Rate in Vertical Ducts , 1995 .

[11]  Bean T. Chen,et al.  Particle Deposition in a Cast of Human Oral Airways , 1999 .

[12]  C. Rhie,et al.  A numerical study of the turbulent flow past an isolated airfoil with trailing edge separation , 1982 .

[13]  D. Lilly,et al.  A proposed modification of the Germano subgrid‐scale closure method , 1992 .

[14]  Zhong-hua Yang,et al.  Multiple laminar flows through curved pipes , 1986 .

[15]  D. Pui,et al.  Numerical study of particle deposition in bends of a circular cross-section-laminar flow regime , 1990 .

[16]  Warren H. Finlay,et al.  On the suitability of k–ε turbulence modeling for aerosol deposition in the mouth and throat: a comparison with experiment , 2000 .

[17]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[18]  W. Uijttewaal,et al.  Particle dispersion and deposition in direct numerical and large eddy simulations of vertical pipe flows , 1996 .

[19]  W. Finlay,et al.  In vitro monodisperse aerosol deposition in a mouth and throat with six different inhalation devices. , 2001, Journal of aerosol medicine : the official journal of the International Society for Aerosols in Medicine.

[20]  M. Breuer LARGE EDDY SIMULATION OF THE SUBCRITICAL FLOW PAST A CIRCULAR CYLINDER: NUMERICAL AND MODELING ASPECTS , 1998 .

[21]  H. L. Stone ITERATIVE SOLUTION OF IMPLICIT APPROXIMATIONS OF MULTIDIMENSIONAL PARTIAL DIFFERENTIAL EQUATIONS , 1968 .

[22]  M. Breuer A CHALLENGING TEST CASE FOR LARGE EDDY SIMULATION: HIGH REYNOLDS NUMBER CIRCULAR CYLINDER FLOW , 2000, Proceeding of First Symposium on Turbulence and Shear Flow Phenomena.

[23]  Krishnaswamy Nandakumar,et al.  Bifurcation in steady laminar flow through curved tubes , 1982, Journal of Fluid Mechanics.

[24]  W. Finlay,et al.  Improved numerical simulation of aerosol deposition in an idealized mouth-throat , 2004 .

[25]  U. Piomelli,et al.  Effect of the subgrid scales on particle motion , 1999 .

[26]  S. Jayanti,et al.  A numerical study of bifurcation in laminar flow in curved ducts , 1992 .

[27]  Y. Cheng,et al.  Motion of particles in bends of circular pipes , 1981 .

[28]  Nakayama Wataru,et al.  Study on forced convective heat transfer in curved pipes: (1st report, laminar region) , 1965 .