Modeling the bifurcating flow in a CT-scanned human lung airway

The inspiratory flow characteristics in a CT-scanned human lung model were numerically investigated using low Reynolds number (LRN) kappa-omega turbulent model. The five-generation airway is extracted from the trachea to segmental bronchi of a 60-year-old Chinese male patient. Computations were carried out in the Reynolds number range of 900-2100, corresponding to mouth-air breathing rates of 190-440ml/s. Flow patterns on the Re=2100 and flow rate distribution were presented. In this model, the flow pattern is very complex. To count the effect of laryngeal jet on trachea inlet, the trachea was extended and modified to simulate the larynx, consequently the inlet velocity profile is biased towards the rear wall. In the inferior lobar bronchi, there are two stems in which the axial velocity is stronger but secondary velocity is weaker. Secondary flow in the lateral bronchi is stronger than the medial ones. With increasing Re, the air flow increases in the middle, inferior lobes and left main bronchus, i.e., flow biases to left and downward.     

Effects of oral airway geometry characteristics on the diffusional deposition of inhaled nanoparticles

The deposition of ultrafine aerosols in the respiratory tract presents a significant health risk due to the increased cellular-level response that these particles may invoke. However, the effects of geometric simplifications on local and regional nanoparticle depositions remain unknown for the oral airway and throughout the respiratory tract. The objective of this study is to assess the effects of geometric simplifications on diffusional transport and deposition characteristics of inhaled ultrafine aerosols in models of the extrathoracic oral airway. A realistic model of the oral airway with the nasopharynx (NP) included has been constructed based on computed tomography scans of a healthy adult in conjunction with measurements reported in the literature. Three other geometries with descending degrees of physical realism were then constructed with successive geometric simplifications of the realistic model. A validated low Reynolds number k-omega turbulence model was employed to simulate laminar, transitional, and fully turbulent flow regimes for the transport of 1-200 nm particles. Results of this study indicate that the geometric simplifications considered did not significantly affect the total deposition efficiency or maximum local deposition enhancement of nanoparticles. However, particle transport dynamics and the underlying flow characteristics such as separation, turbulence intensity, and secondary motions did show an observable sensitivity to the geometric complexity. The orientation of the upper trachea was shown to be a major factor determining local deposition downstream of the glottis and should be retained in future models of the respiratory tract. In contrast, retaining the NP produced negligible variations in airway dynamics and could be excluded for predominantly oral breathing conditions. Results of this study corroborate the use of existing diffusion correlations based on a circular oral airway model. In comparison to previous studies, an improved correlation for the deposition of nanoparticles was developed based on a wider range of particle sizes and flow rates, which captures the dependence of the Sherwood number on both Reynolds and Schmidt numbers.     

Validating CFD predictions of respiratory aerosol deposition: effects of upstream transition and turbulence

A number of computational fluid dynamics (CFD) studies have explored local deposition patterns of inhaled aerosols in the respiratory tract. These studies have highlighted the effects of multiple physiologic, geometric, and particle characteristics on deposition. However, very few studies have reported local or sub-branch quantitative comparisons to in vitro particle deposition data. The objective of this study is to numerically investigate the effects of transition and turbulence on highly localized particle deposition in a respiratory double bifurcation model in order to quantitatively validate CFD results. To perform the validations, local comparisons have been made to a specific in vitro case study of 10 microm particles depositing in a model of respiratory generations G3-G5. To achieve this objective, two geometric cases have been considered. The first case includes only the double bifurcation model. The second case includes a portion of the experimental particle delivery geometry, where transitional flow is expected. To evaluate the effectiveness of two-equation turbulence models in this system, the flow field solutions have been computed using laminar, standard k-omega, and low Reynolds number (LRN) k-omega approximations. Results indicate that even though the Reynolds number remained below the critical limit required for full turbulence, transition and turbulence have a significant impact on the flow field and local particle deposition patterns. For the experimental case considered, turbulence impacted the local deposition of 10 microm particles primarily by influencing the initial velocity and particle profiles. As such, both the laminar and LRN k-omega flow models provided good local quantitative matches to the in vitro deposition data, provided that the correct initial particle profile was specified. Implications of this study include the need for local quantitative validations of particle deposition results, the importance of correct inlet conditions, and the need to consider upstream effects in experimental and computational studies of the respiratory tract. Applications of these results to realistic respiratory geometries will require consideration on upstream flow conditions in the lung, transient flow, and intermittent turbulent structures.     

Computational simulations of airflow in an in vitro model of the pediatric upper airways

In order to understand mechanisms of gas and aerosol transport in the human respiratory system airflow in the upper airways of a pediatric subject (male aged 5) was calculated using Computational Fluid Dynamic techniques. An in vitro reconstruction of the subject's anatomy was produced from MRI images. Flow fields were solved for steady inhalation at 6.4 and 8 LPM. For validation of the numerical solution, airflow in an adult cadaver based trachea was solved using identical numerical methods. Comparisons were made between experimental results and computational data of the adult model to determine solution validity. It was found that numerical simulations can provide an accurate representation of axial velocities and turbulence intensity. Data on flow resistance, axial velocities, secondary velocity vectors, and turbulent kinetic energy are presented for the pediatric case. Turbulent kinetic energy and axial velocities were heavily dependant on flow rate, whereas turbulence intensity varied less over the flow rates studied. The laryngeal jet from an adult model was compared to the laryngeal jet in the pediatric model based on Tracheal Reynolds number. The pediatric case indicated that children show axial velocities in the laryngeal jet comparable to adults, who have much higher tracheal Reynolds numbers than children due to larger characteristic dimensions. The intensity of turbulence follows a similar trend, with higher turbulent kinetic energy levels in the pediatric model than would be expected from measurements in adults at similar tracheal Reynolds numbers. There was reasonable agreement between the location of flow structures between adults and children, suggesting that an unknown length scale correlation factor could exist that would produce acceptable predictions of pediatric velocimetry based off of adult data sets. A combined scale for turbulent intensity as well may not exist due to the complex nature of turbulence production and dissipation.     

The effect of inlet flow profile and nozzle diameter on drug delivery to the maxillary sinus

In this paper, the effect of the turbulence and swirling of the inlet flow and the diameter of the nozzle on the flow characteristics and the particles' transport/deposition patterns in a realistic combination of the nasal cavity (NC) and the maxillary sinus (MS) were examined. A computational fluid dynamics (CFD) model was developed in ANSYS® Fluent using a hybrid Reynolds averaged Navier–Stokes–large-eddy simulation algorithm. For the validation of the CFD model, the pressure distribution in the NC was compared with the experimental data available in the literature. An Eulerian–Lagrangian approach was employed for the prediction of the particle trajectories using a discrete phase model. Different inlet flow conditions were investigated, with turbulence intensities of 0.15 and 0.3, and swirl numbers of 0.6 and 0.9 applied to the inlet flow at a flow rate of 7 L/min. Monodispersed particles with a diameter of 5 µm were released into the nostril for various nozzle diameters. The results demonstrate that the nasal valve plays a key role in nasal resistance, which damps the turbulence and swirl intensities of the inlet flow. Moreover, it was found that the effect of turbulence at the inlet of the NC on drug delivery to the MS is negligible. It was also demonstrated that increasing the flow swirl at the inlet and decreasing the nozzle diameter improves the total particle deposition more than threefold due to the generation of the centrifugal force, which acts on the particles in the nostril and vestibule. The results also suggest that the drug delivery efficiency to the MS can be increased by using a swirling flow with a moderate swirl number of 0.6. It was found that decreasing the nozzle diameter can increase drug delivery to the proximity of the ostium in the middle meatus by more than 45%, which subsequently increases the drug delivery to the MS. The results can help engineers design a nebulizer to improve the efficiency of drug delivery to the maxillary sinuses.

     

Analytical model of hygroscopic particle behavior in human airways

A model is developed to calculate the deposition of hygroscopic aerosols in the human tracheobronchial (TB) tree. The TB airflow pattern assumed is consistent with experimental observations and accounts for anatomical features such as the larynx and cartilaginous rings in large airways. Some original deposition efficiency formulae are presented for laminar and turbulent airstreams. Stepwise growth is simulated by changes in particle size and density at each TB generation. The dose distribution of NaCl aerosols is studied as a function of inhaled particle size and flow rate. Two NaCl growth rate curves are used which differ in the mode of aerosol-air mixing in the trachea. The initial rate of aerosol mixing in the human due to the laryngeal jet is shown to be an important factor affecting the deposition of hygroscopic aerosols. Total TB deposition of NaCl exceeds that for nonhygroscopic particles of the same inhaled aerodynamic size. Hygroscopic growth can also influence the regional TB distribution of dose when submicron NaCl particles grow rapidly enough to deposit by impaction and sedimentation.     

Estimation of turbulent shear stresses in pulsatile flow immediately downstream of two artificial aortic valves in vitro

Measuring turbulent shear stresses is of major importance in artificial heart valve evaluation. Bi- and unidirectional fluid velocity measurements enable calculation of Reynolds shear stress ( ) and Reynolds normal stress ( ). τ is important due to the relation to hemolysis and thrombus formation, but σ is the only obtainable parameter in vivo. Therefore, determination of a correlation factor between τ and σ is pertinent.In a pulsatile flow model, laser Doppler (LDA) and hot-film (HFA) anemometry were used for simultaneous bi- and unidirectional fluid velocity measurements downstream of a Hall Kaster and a Hancock Porcine aortic valve. Velocities were registered in two flow field locations and at four cardiac outputs. The velocity signals were subjected to analog signal processing prior to digital turbulence analysis, as a basis for calculation of τ and σ.A correlation factor of 0.5 with a correlation coefficient of 0.97 was found between the maximum Reynolds shear stress and Reynolds normal stress, implying . In vitro estimation of turbulent shear stresses downstream of artificial aortic valves, based on the axial velocity component alone, seems possible.     

Effects of the laryngeal jet on nano- and microparticle transport and deposition in an approximate model of the upper tracheobronchial airways.

The extent to which laryngeal-induced flow features penetrate into the upper tracheobronchial (TB) airways and their related impact on particle transport and deposition are not well understood. The objective of this study was to evaluate the effects of including the laryngeal jet on the behavior and fate of inhaled aerosols in an approximate model of the upper TB region. The upper TB model was based on a simplified numerical reproduction of a replica cast geometry used in previous in vitro deposition experiments that extended to the sixth respiratory generation along some paths. Simulations with and without an approximate larynx were performed. Particle sizes ranging from 2.5 nm to 12 mum were considered using a well-tested Lagrangian tracking model. The model larynx was observed to significantly affect flow dynamics, including a laryngeal jet skewed toward the right wall of the trachea and a significant reverse flow in the left region of the trachea. Inclusion of the laryngeal model increased the tracheal deposition of nano- and micrometer particles by factors ranging from 2 to 10 and significantly reduced deposition in the first three bronchi of the model. Considering localized conditions, inclusion of the laryngeal approximation decreased deposition at the main carina and produced a maximum in local surface deposition density in the lobar-to-segmental bifurcations (G2-G3) for both 40-nm and 4-microm aerosols. These findings corroborate previous experiments and highlight the need to include a laryngeal representation in future computational and in vitro models of the TB region.     

Effects of curved inlet tubes on air flow and particle deposition in bifurcating lung models

In vivo bifurcating airways are complex and the airway segments leading to the bifurcations are not always straight, but curved to various degrees. How do such curved inlet tubes influence the motion as well as local deposition and hence the biological responses of inhaled particulate matter in lung airways? In this paper steady laminar dilute suspension flows of micron-particles are simulated in realistic double bifurcations with curved inlet tubes, i.e., 0 degrees < or =theta< or =90 degrees, using a commercial finite-volume code with user-enhanced programs. The resulting air-flow patterns as well as particle transport and wall depositions were analyzed for different flow inlet conditions, i.e., uniform and parabolic velocity profiles, and geometric configurations. The curved inlet segments have quite pronounced effects on air-flow, particle motion and wall deposition in the downstream bifurcating airways. In contrast to straight double bifurcations, those with bent parent tubes also exhibit irregular variations in particle deposition efficiencies as a function of Stokes number and Reynolds number. There are fewer particles deposited at mildly curved inlet segments, but the particle deposition efficiencies at the downstream sequential bifurcations vary much when compared to those with straight inlets. Under certain flow conditions in sharply curved lung airways, relatively high, localized particle depositions may take place. The findings provide necessary information for toxicologic or therapeutic impact assessments and for global lung dosimetry models of inhaled particulate matter.     

Study on biomass circulation and gasification performance in a clapboard-type internal circulating fluidized bed gasifier

We investigated the solid particle flow characteristics and biomass gasification in a clapboard-type internal circulating fluidized bed reactor. The effect of fluidization velocity on particle circulation rate and pressure distribution in the bed showed that fluidization velocities in the high and low velocity zones were the main operational parameters controlling particle circulation. The maximum internal circulation rates in the low velocity zone came almost within the range of velocities in the high velocity zone, when u_H/u_mf = 2.2–2.4 for rice husk and u_H/u_mf = 3.5–4.5 for quartz sand. In the gasification experiment, the air equvalence ratio (ER) was the main controlling parameter. Rice husk gasification gas had a maximum heating value of around 5000 kJ/m³ when ER = 0.22–0.26, and sawdust gasification gas reached around 6000–6500 kJ/m³ when ER = 0.175–0.24. The gasification efficiency of rice husk reached a maximum of 77% at ER = 0.28, while the gasification efficiency of sawdust reached a maximum of 81% at ER = 0.25.     

Large eddy simulation of the unsteady flow-field in an idealized human mouth-throat configuration

The present study concerns the simulation and analysis of the flow field in the upper human respiratory system in order to gain an improved understanding of the complex flow field with respect to the process affecting drug delivery for medical treatment of the human air system. For this purpose, large eddy simulation (LES) is chosen because of its powerful performance in the transitional range of laminar and turbulent flow fields. The average gas velocity in a constricted tube is compared with experimental data (Ahmed and Giddens, 1983) and numerical data from Reynolds-averaged Navier-Stokes (RANS) equations coupled with low Reynolds number (LRN) κ-ω model (Zhang and Kleinstreuer, 2003) and LRN shear-stress transport κ-ω model (Jayaraju et al., 2007), for model validation. The present study emphasizes on the instantaneous flow field, where the simulations capture different scales of secondary vortices in different flow zones including recirculation zones, the laryngeal jet zone, the mixing zone, and the wall shear layer. It is observed that the laryngeal jet tail breaks up, and the unsteady motion of laryngeal jet is coupled with the unsteady distribution of secondary vortices in the jet boundary. The present results show that it is essential to study the unsteady flow field since it strongly affects the particle flow in the human upper respiratory system associated with drug delivery for medical treatment.     

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相似文献   

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 μm≤ d p ≤7 μm), 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.     

Aerosol transport and deposition in sequentially bifurcating airways

Deposition patterns and efficiencies of a dilute suspension of inhaled particles in three-dimensional double bifurcating airway models for both in-plane and 90 deg out-of-plane configurations have been numerically simulated assuming steady, laminar, constant-property air flow with symmetry about the first bifurcation. Particle diameters of 3, 5, and 7 microns were used in the simulation, while the inlet Stokes and Reynolds numbers varied from 0.037 to 0.23 and 500 to 2000, respectively. Comparisons between these results and experimental data based on the same geometric configuration showed good agreement. The overall trend of the particle deposition efficiency, i.e., an exponential increase with Stokes number, was somewhat similar for all bifurcations. However, the deposition efficiency of the first bifurcation was always larger than that of the second bifurcation, while in general the particle efficiency of the out-of-plane configuration was larger than that of the in-plane configuration. The local deposition patterns consistently showed that the majority of the deposition occurred in the carinal region. The distribution pattern in the first bifurcation for both configurations were symmetric about the carina, which was a direct result of the uniaxial flow at the inlet. The deposition patterns about the second carina showed increased asymmetry due to highly nonuniform flow generated by the first bifurcation and were extremely sensitive to bifurcation orientation. Based on the deposition variations between bifurcation levels and orientations, the use of single bifurcation models was determined to be inadequate to resolve the complex fluid-particle interactions that occur in multigenerational airways.     

The drift of zooplankton and microzoobenthos in the river Strandaelva,western Norway

The drift of zooplankton (rotifers, cladocerans, cyclopoid copepods) and microscopical zoobenthos (mainly bdelloid rotifers and small chironomid larvae) was investigated by filtering samples of river water. The number of drifting benthic rotifers varied between 1 000 and 6 000 ind. m^–3 in the lake inlet, and between 30 and 500 ind. m^–3 in the lake outlet, without any seasonal trend. The number of drifting insect larvae was approx. equal in the lake inlet and outlet, with a maximum in summer (250–300 ind. m^–3) and minimum in winter (ca. 10 ind. m^–3). Increasing water flow resulted in an increasing number of drifting zoobenthos. Downstream from the lake, the number of drifting benthic rotifers was increasing from approx. 300 ind. m^–3 in the outlet to 6 500 ind. m^–3 3.4 km downstream, while the number of insect larvae was ca. 100 ind. m^–3 in the outlet and leveled off at approx. 300 ind. m^–3 after 200 m. The number of drifting zooplankton in the lake outlet varied between 20 and 2 000 ind. m^–3 for crustaceans, and between 300 and 20 000 ind. m^–3 for rotifers, both with a maximum in late summer/autumn and a minimum in winter. The number of drifting zooplankton decreased by some 45% in the first 200 m from the lake outlet, but some zooplankton was still found in the drift 3.4 km downstream. The largest species was removed first from the drift. The diurnal variation in the number of drifting zooplankton in lake outlets appear to be related to the vertical migration in the lake, i.e. the largest number drifting when most animals are in the upper water layers.Contribution from the Voss Project, University of OsloContribution from the Voss Project, University of Oslo     

Measurement of steady-flow instability and turbulence levels in Dacron vascular grafts.

Fluid dynamic properties of Dacron vascular grafts were studied under controlled steady-flow conditions over a Reynolds number range of 800 to 4500. Knitted and woven grafts having nominal diameters of 6 mm and 10 mm were studied. Thermal anemometry was used to measure centerline velocity at the downstream end of the graft; pressure drop across the graft was also measured. Transition from laminar flow to turbulent flow was observed, and turbulence intensity and turbulent stresses (Reynolds normal stresses) were measured in the turbulent regime. Knitted grafts were found to have greater pressure drop than the woven grafts, and one sample was found to have a critical Reynolds number (Rc) of less than one-half the value of Rc for a smooth-walled tube.     

Viscous energy dissipation for steady flow in models of arterial bifurcations

In this paper the viscous energy dissipation in a series of glass model symmetric bifurcations—typical of human vascular branching—was studied. The bifurcations studied have included angles of 75°, 100° and 125° and total output/input area ratios of 0.73, 1.07 and 1.33. The flowrate range studied corresponded to parent tube Reynolds numbers in the range 100–1000.Pressure and flow measurements were made using a highly sensitive variable reluctance pressure transducer and electromagnetic flowmeter. The measurements were made in such a way as to indicate the net effect of the bifurcation.It was found that a dimensionless form of the viscous dissipation had a constant value for all the geometries investigated up to a Reynolds number of 800. Above this, for an angle of 125°, the measured energy dissipation increased. An analytical model based on entry flow principles showed good agreement with measured values except at an area ratio of 0.73. The reasons for this are discussed with the physiological implications of the results.     

Sources,fates, and impacts of nitrogen inputs to terrestrial ecosystems: review and synthesis

The relative importance of nitrogen inputs from atmospheric deposition and biological fixation is reviewed in a number of diverse, non-agricultural terrestrial ecosystems. Bulk precipitation inputs of N (l–l2 kg N ha^–1 yr^–1) are the same order of magnitude as, or frequently larger than, the usual range of inputs from nonsymbiotic fixation (< 1=" –=" 5=" kg=" n=">^–1 yr^–1), especially in areas influenced by industrial activity. Bulk precipitation measurements may underestimate total atmospheric deposition by 30–40% because they generally do not include all forms of wet and dry deposition. Symbiotic fixation generally ranges from 10–160 kg N ha^–1 yr^–1) in ecosystems where N-fixing species are present during early successional stages, and may exceed the range under unusual conditions.Rates of both symbiotic and nonsymbiotic fixation appear to be greater during early successional stages of forest development, where they have major impacts on nitrogen dynamics and ecosystem productivity. Fates and impacts of these nitrogen inputs are important considerations that are inadequately understood. These input processes are highly variable in space and time, and few sites have adequate comparative information on both nitrogen deposition and fixation.

Modelling of parameters for optimization of maturity in composting trimming residues

In this paper, changes in physico-chemical parameters during trimmings residue composting (cation exchange capacity, germination index, self-heated, NH₄/NO₃ ratio and C_FA/C_HA ratio) in relation to environmental composting parameters (time, aeration, moisture and particle size) of the composting process were studied. A central composite experimental design was used to obtain the polynomial model for each dependent variable. Results of the modelling showed that among the studied range, moisture was the highest influenced parameter in maturity evaluation, with respect to aeration and particle size. An exception was found for CEC evolution. In this parameter, the highest influence was found for particle size. Moreover, a product with acceptable chemical properties entails operating at medium moisture content (55%) and medium-to-high particle size (3–5 cm). Moderate to low aeration (0.2 m³ air kg⁻¹ d⁻¹) would be the best compromise to composting this residue, due to the scarce statistical influence of this independent variable.     

A method for low volume and low Se concentration samples and application to paired cerebrospinal fluid and serum samples

The well-known beneficial health effects of Se have demanded the development of rapid and accurate methods for its analysis. A flow injection (FI) method with inductively coupled plasma mass spectrometry (ICP-MS) as a selenium-selective detector was optimized. Flow injection was carried out using a Knauer 1100 smartline inert series liquid chromatograph coupled with a Perkin Elmer DRC II ICP-mass spectrometer. For sample injection a Perkin Elmer electronic valve equipped with a 25 μL sample loop was employed. Before measurement, standards or samples were administered with 1 μg/L rhodium as internal standard for correction of changes in detector response according to changes in sample electrolyte concentration. The method characterization parameters are: LOD (3σ criterion): 26 ng/L, LOQ (10σ criterion): 86 ng/L, linearity: 0.05–>10 μg/L, r²=0.9999, serial or day-to-day precision at 2 μg/L: 4.48% or 5.6%. Accuracy was determined by (a) recovery experiments (CSF spiked with 2 μg/L Se); (b) comparison of FI-ICP-MS measurement with graphite furnace atomic absorption (GFAAS) measurements of 1:10 diluted serum samples; (c) Se determination in urine and serum control materials. Recovery (a) was 101.4%, measurement comparison with GFAAS (b) showed 98.8% (5 serum samples, 1:10 diluted in the range of 0.5–1.3 μg/L, compared to GFAAS determination, which was set to 100%), and accuracy was 96.8% or 105.6% for the serum or urine control material. Analysis time per sample was short and typically below 2 min for the complete measurement, including sample introduction, sample-line purge and quadruplicate Se determination.This method was used to determine Se in cerebrospinal fluid (CSF) and plasma (here parallel to GFAAS) in 35 paired serum and CSF samples. Se determination gave values in the range of 42–130 μg/L for serum and 1.63–6.66 μg/L for CSF. The median for Se in 35 individual CSF samples was 3.28 μg/L, the mean (±SD) was 3.67 (1.35) μg/L, whilst for individual serum samples the median was 81 μg/L and the mean (±SD) was 85 (26) μg/L. When relating the paired Se concentrations of CSF samples to respective serum samples it turned out that Se-CSF (behind blood brain barrier (BBB)) is independent on Se-serum concentration (before BBB).