Effects of asymmetric branch flow rates on aerosol deposition in bifurcating airways

Visualizing particle/flow structure interactions in the small bronchial tubes

Particle deposition in the small bronchial tubes (generations six through twelve) is strongly influenced by the vortex-dominated secondary flows that are induced by axial curvature of the tubes. In this paper, we employ particle destination maps in conjunction with two-dimensional, finite-time Lyapunov exponent maps to illustrate how the trajectories of finite-mass particles are influenced by the presence of vortices. We consider two three-generation bronchial tube models: a planar, asymmetric geometry and a non-planar, asymmetric geometry. Our visualizations demonstrate that these techniques, coupled with judiciously seeded particle trajectories, are effective tools for studying particle/flow structure interactions.

Three-Dimensional Simulations of Reactive Gas Uptake in Single Airway Bifurcations

The pattern of lung injury induced by the inhalation of ozone (O(3)) depends on the dose delivered to different tissues in the airways. This study examined the distribution of O(3) uptake in a single, symmetrically branched airway bifurcation. Reaction in the epithelial lining fluid was assumed to be so rapid that O(3) concentration was negligible along the entire surface of the bifurcation wall. Three-dimensional numerical solutions of the continuity, Navier-Stokes and convection-diffusion equations were obtained for steady inspiratory and expiratory flows at Reynolds numbers ranging from 100 to 500. The total rate of O(3) uptake was found to increase with increasing flow rate during both inspiration and expiration. Hot spots of O(3) flux appeared at the carina of the bifurcation for virtually all inspiratory and expiratory Reynolds numbers considered in the simulations. At the lowest expiratory Reynolds number, however, the location of the maximum flux was shifted to the outer wall of the daughter branch. For expiratory flow, additional hot spots of flux were found on the parent branch wall just downstream of the branching region. In all cases, O(3) uptake in the single bifurcation was larger than that in a straight tube of equal inlet radius and wall surface area. This study provides insight into the effect of flow conditions on O(3) uptake and dose distribution in individual bifurcations.

Effects of mesh style and grid convergence on particle deposition in bifurcating airway models with comparisons to experimental data

A number of research studies have employed a wide variety of mesh styles and levels of grid convergence to assess velocity fields and particle deposition patterns in models of branching biological systems. Generating structured meshes based on hexahedral elements requires significant time and effort; however, these meshes are often associated with high quality solutions. Unstructured meshes that employ tetrahedral elements can be constructed much faster but may increase levels of numerical diffusion, especially in tubular flow systems with a primary flow direction. The objective of this study is to better establish the effects of mesh generation techniques and grid convergence on velocity fields and particle deposition patterns in bifurcating respiratory models. In order to achieve this objective, four widely used mesh styles including structured hexahedral, unstructured tetrahedral, flow adaptive tetrahedral, and hybrid grids have been considered for two respiratory airway configurations. Initial particle conditions tested are based on the inlet velocity profile or the local inlet mass flow rate. Accuracy of the simulations has been assessed by comparisons to experimental in vitro data available in the literature for the steady-state velocity field in a single bifurcation model as well as the local particle deposition fraction in a double bifurcation model. Quantitative grid convergence was assessed based on a grid convergence index (GCI), which accounts for the degree of grid refinement. The hexahedral mesh was observed to have GCI values that were an order of magnitude below the unstructured tetrahedral mesh values for all resolutions considered. Moreover, the hexahedral mesh style provided GCI values of approximately 1% and reduced run times by a factor of 3. Based on comparisons to empirical data, it was shown that inlet particle seedings should be consistent with the local inlet mass flow rate. Furthermore, the mesh style was found to have an observable effect on cumulative particle depositions with the hexahedral solution most closely matching empirical results. Future studies are needed to assess other mesh generation options including various forms of the hybrid configuration and unstructured hexahedral meshes.

The effect of inlet velocity profile on the bifurcation COPD airway flow

The effect of inlet velocity profile on the flow features in obstructed airways is investigated in this study. In reality, the inlet velocity distributions on such models, which are extracted from medial branches of natural human lung, should be neither uniform, nor symmetric parabolic, but skewed-parabolic due to having been skewed by the upper carina ridges. Four different three-dimensional three-generation models based on the 23 generations model of Weibel have been considered, respectively. The fully three-dimensional incompressible laminar Navier-Stokes equations and continuity equation have been solved using CFD solver on unstructured tetrahedral meshes. To reduce the complexity of the simulations, only one Reynolds number of 900 was used in this calculation. Four types of inlet boundary conditions, namely uniform, parabolic, positive-skewed parabolic (skewed to the positive x-direction), and negative-skewed parabolic, were imposed on the obstructed airway models, which were considered to be obstructed at either the second generation or the third generation airways, respectively. The results show that the inlet velocity profile has significant influence on the flow patterns, mass distributions, and pressure drops in either the symmetric model, or the three obstructed models. The three generation airways may not be enough to study the bifurcation flow in chronic obstructive pulmonary disease (COPD) airways, and a four-generation or more airway model is necessary to get better predictive results.

Modeling the bifurcating flow in an asymmetric human lung airway

In a former paper, the inspiratory flow characteristics in a three-generation symmetric bifurcation airway have been numerically investigated using a control volume method to solve the fully three-dimensional laminar Navier-Stokes equations. The present paper extends the work to deal with asymmetric airway extracted from the 5th-11th branches of the model of Weibel (Morphometry of the Human Lung. New York Academic Press, Verlag, 1963) in order to more appropriately model human air passage. Computations are carried out in the Reynolds number range 200-1600, corresponding to mouth-air breathing rates ranging from 0.27 to 2.16l/s, representative for an averaged height man breathing from quiet to vigorous state. Particular attentions are paid to establishing relations between the Reynolds number and the overall flow characteristics, including flow patterns and pressure drop. The study shows that the ratios of airflow rate through the medial branches to that of their mother branches are the same, and this is also true for the ratios of airflow rate through the lateral branches. This partially explains why regular human breathing is not affected by airways of different sizes.

Targeted drug aerosol deposition analysis for a four-generation lung airway model with hemispherical tumors

One important research area of broad interest is the development of highly efficient drug delivery systems for desired site deposition and uptake. For example, controlled drug aerosol release and targeting to specific regions of the lung is a novel way to combat lung diseases, diabetes, virus infections, cancers, etc. Determination of feasible air-particle streams is a prerequisite for the development of such delivery devices, say, smart inhalers. The concept of "controlled particle release and targeting" is introduced and results are discussed for a representative model of bronchial lung airways afflicted with hemispherical tumors of different sizes and locations. It is shown that under normal particle inlet conditions a particle mass fraction of only up to 11% may deposit on the surface of a specific tumor with critical radius r/R approximately 1.25, while a controlled particle release achieves deposition fractions of 35 to 92% for a realistic combination of inlet Stokes and Reynolds numbers, depending mainly on tumor size. Furthermore, with the controlled release and targeting approach nearby healthy tissue is hardly impacted by the typically aggressive drug aerosols. Assuming laminar, quasi-steady, three-dimensional air flow and spherical non-interacting micron-particles in sequentially bifurcating rigid airways, the results were obtained using a validated commercial finite-volume code with user-enhanced programs on a high-end engineering workstation. The new concept is generic and hence should be applicable to other regions of the respiratory system as well.

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</=d(p)</=7 microm), 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.