Dynamics of disk-like particles in turbulent vertical channel flow

The dynamical behavior of inertial disk-like particles in turbulent vertical channel flow is investigated by an Eulerian–Lagrangian point-particle approach. Gravity effects on distribution, translation, rotation and orientation statistics of non-spherical particles modeled as oblate spheroids have been studied both in an upward and a downward flow and compared with results obtained in the absence of gravity. Altogether 12 different particle classes have been studied, with inertia and shape parameterized by means of Stokes number St and aspect ratio λ ≤ 1. The St = 5 disk-like particles distribute more evenly across the channel in upward than in downward flow. The gravity effect on the particle concentration diminishes with large inertia and the spheroid shape has only a modest influence. Although the gravity significantly affects the streamwise and wall-normal mean slip velocities with increasing inertia, the particle shape rarely has any impact on the translational motion, except for the mean wall-normal velocity. The fluctuations of the velocity of disk-like particles are mainly ascribed to inertia, whereas the gravity and shape only have marginal effects. The presence of gravity is moreover found to have a negligible effect on the particles’ orientation and rotation, in spite of the striking effect of λ on the orientation and rotation seen in the near-wall region. The tendency of the disks to align their symmetry axis orthogonal to the fluid vorticity in the channel center is stronger for particles with modest inertia. In the near-wall region, however, oblate spheroids preferentially align with the fluid vorticity for St >> 1. The observed behavior is believed to be caused by the influence of the gravity force on the turbophoresis; i.e. that inertial particles move towards low-turbulence regions.     

Particle deposition model for particulate flows at high temperatures in gas turbine components

This study proposes an improved physical model to predict sand deposition at high temperature in gas turbine components. This model differs from its predecessor (Sreedharan and Tafti, 2011) by improving the sticking probability by accounting for the energy losses during particle-wall collision based on our previous work (Singh and Tafti, 2013). This model predicts the probability of sticking based on the critical viscosity approach and collision losses during a particle–wall collision. The current model is novel in the sense that it predicts the sticking probability based on the impact velocity along with the particle temperature. To test the model, deposition from a sand particle laden jet impacting on a flat coupon geometry is computed and the results obtained from the numerical model are compared with experiments (Delimont et al., 2014) conducted at Virginia Tech, on a similar geometry and flow conditions, for jet temperatures of 950 °C, 1000 °C and 1050 °C. Large Eddy Simulations (LES) are used to model the flow field and heat transfer, and sand particles are modeled using a discrete Lagrangian framework. Results quantify the impingement and deposition for 20–40 μm sand particles. The stagnation region of the target coupon is found to experience most of the impingement and deposition. For 950 °C jet temperature, around 5% of the particle impacting the coupon deposit while the deposition for 1000 °C and 1050 °C is 17% and 28%, respectively. In general, the sticking efficiencies calculated from the model show good agreement with the experiments for the temperature range considered.     

Instantaneous flow field above the free end of finite-height cylinders and prisms

The flow above the free ends of surface-mounted finite-height circular cylinders and square prisms was studied experimentally using particle image velocimetry (PIV). Cylinders and prisms with aspect ratios of AR = 9, 7, 5, and 3 were tested at a Reynolds number of Re = 4.2 × 10⁴. The bodies were mounted normal to a ground plane and were partially immersed in a turbulent zero-pressure-gradient boundary layer, where the boundary layer thickness relative to the body width was δ/D = 1.6. PIV measurements were made above the free ends of the bodies in a vertical plane aligned with the flow centreline. The present PIV results provide insight into the effects of aspect ratio and body shape on the instantaneous flow field. The recirculation zone under the separated shear layer is larger for the square prism of AR = 3 compared to the more slender prism of AR = 9. Also, for a square prism with low aspect ratio (AR = 3), the influence of the reverse flow over the free end surface becomes more significant compared to that for a higher aspect ratio (AR = 9). For the circular cylinder, a cross-stream vortex forms within the recirculation zone. As the aspect ratio of the cylinder decreases, the reattachment point of the separated flow on the free end surface moves closer to the trailing edge. For both the square prism and circular cylinder cases, the instantaneous velocity vector field and associated in-plane vorticity field revealed small-scale structures mostly generated by the separated shear layer.     

Numerical investigation of particle deposition inside aero-shielded solar cyclone reactor: A promising solution for reactor clogging

Solar cracking of methane is considered to be an attractive option due to its CO₂ free hydrogen production process. Carbon particle deposition on the reactor window, walls and exit is a major obstacle to achieve continuous operation of methane cracking solar reactors. As a solution to this problem a novel “aero-shielded solar cyclone reactor” was created. In this present study the prediction of particle deposition at various locations for the aero-shielded reactor is numerically investigated by a Lagrangian particle dispersion model. A detailed three dimensional computational fluid dynamic (CFD) analysis for carbon deposition at the reactor window, walls and exit is presented using a Discrete Phase Model (DPM). The flow field is based on a RNG k–ε model and species transport with methane as the main flow and argon/ hydrogen as window and wall screening fluid. Flow behavior and particle deposition have been observed with the variation of main flow rates from 10–20 L/min and with carbon particle mass flow rate of 7 × 10⁻⁶ and 1.75 × 10⁻⁵ kg/s. In this study the window and wall screening flow rates have been considered to be 1 L/min and 10 L/min by employing either argon or hydrogen. Also, to study the effect of particle size simulations have also been carried out (i) with a variation of particle diameter with a size distribution of 0.5–234 μm and (ii) by taking 40 μm mono sized particles which is the mean value for the considered size distribution. Results show that by appropriately selecting the above parameters, the concept of the aero-shielded reactor can be an attractive option to resolve the problem of carbon deposition at the window, walls and exit of the reactor.     

LES-based filter-matrix lattice Boltzmann model for simulating turbulent natural convection in a square cavity

In this paper, a novel thermal filter-matrix lattice Boltzmann model based on large eddy simulation (LES) is proposed for simulating turbulent natural convection. In this study, the Vreman subgrid-scale eddy-viscosity model is introduced into the present framework of LES to accurately predict the flow in near-wall region. Two dimensional numerical simulations of natural convection in a square cavity were performed at high Rayleigh number varying from 10⁷ to 10¹⁰ with a fixed Prandtl number of Pr = 0.71. The influences of the higher-order terms upon the present results at high Rayleigh numbers are examined, taking Ra = 10⁷ and 10⁸ as the example, revealing that the proper minimization of the higher-order terms can improve numerical accuracy of present model for high Rayleigh convective flow. For the turbulent convective flow, the time-averaged quantities in the median lines are presented and compared against those available results from previous studies. The general structure of turbulent boundary layers is well predicted. All numerical results exhibit good agreement with the benchmark solutions available in the previous literatures.     

Turbulent boundary layers over sparsely-spaced rod-roughened walls

Direct numerical simulations (DNSs) of spatially developing turbulent boundary layers (TBLs) over sparsely-spaced two-dimensional (2D) rod-roughened walls were performed. The rod elements were periodically arranged along the streamwise direction with pitches of p_x/k = 8, 16, 32, 64 and 128, where p_x is the streamwise spacing of the rods, and k is the roughness height. The Reynolds number based on the momentum thickness was varied from Re_θ = 300–1400, and the height of the roughness element was k = 1.5θ_in, where θ_in is the momentum thickness at the inlet. The characteristics of the TBLs, such as the friction velocity, mean velocity, and Reynolds stresses over the rod-roughened walls, were examined by varying the spacing of the roughness features (8 ⩽ p_x/k ⩽ 128). The outer-layer similarity between the rough and smooth walls was established for the sparsely-distributed rough walls (p_x/k ⩾ 32) based on the profiles of the Reynolds stresses, whereas those are not for p_x/k = 8 and 16. Inspection of the interaction between outer-layer large-scale motions and near-wall small-scale motions using two-point amplitude modulation (AM) covariance showed that modulation effect of large-scale motions on near-wall small-scale motions was strongly disturbed over the rough wall for p_x/k = 8 and 16. For p_x/k ⩾ 32, the flow that passed through the upstream roughness element transitioned to a smooth wall flow between the consecutive rods. The strong influence of the surface roughness in the outer layer for p_x/k = 8 and 16 was attributed to large-scale erupting motions by the surface roughness, creating both upward shift of the near-wall turbulent energy and active energy production in the outer layer with little influence on the near-wall region.     

Six-day measurement of size-resolved indoor fluorescent bioaerosols of outdoor origin in an office

Indoor airborne bioaerosols of outdoor origin play an important role in determining the exposure of humans to bioaerosols because people spend most of their time indoors. However, there are few studies focusing on indoor bioaerosols originating from outdoors. In this study, indoor versus outdoor size-resolved concentrations and particle asymmetry factors of airborne fluorescent bioaerosols in an office room were measured continuously for 6 days (144 h) using a fluorescent bioaerosol detector. The windows and door of this room were closed to ensure that there was only air infiltration; moreover, any human activities were ceased during sampling to inhibit effects of indoor sources. We focused on fine particles, since few coarse particles enter indoor environments, when windows and doors are closed. Both indoor and outdoor fluorescent bioaerosol size distributions were fit with two-mode lognormal distributions (indoor R² = 0.935, outdoor R² = 0.938). Asymmetry factor distributions were also fit with lognormal distributions (indoor R² = 0.992, outdoor R² = 0.992). Correlations between indoor and outdoor fluorescent bioaerosol concentrations show significant concentration-attenuation and a time lag during the study period. A two-parameter, semi-empirical model was used to predict concentrations of indoor fluorescent bioaerosols of outdoor origin. The measured and predicted concentrations had a linear relationship for the studied size fractions, with an R² for all size fractions of larger than 0.83.     

Finite-length effects on dynamical behavior of rod-like particles in wall-bounded turbulent flow

Combined Particle Image Velocimetry (PIV) and Particle Tracking Velocimetry (PTV) measurements have been performed in dilute suspensions of rod-like particles in wall turbulence. PIV results for the turbulence field in the water table flow apparatus compared favorably with data from Direct Numerical Simulations (DNS) of channel flow turbulence and the universality of near-wall turbulence justified comparisons with DNS of fiber-laden channel flow. In order to examine any shape effects on the dynamical behavior of elongated particles in wall-bounded turbulent flow, fibers with three different lengths but the same diameter were used. In the logarithmic part of the wall-layer, the translational fiber velocity was practically unaffected by the fiber length l. In the buffer layer, however, the fiber dynamics turned out to be severely constrained by the distance z to the wall. The short fibers accumulated preferentially in low-speed areas and adhered to the local fluid speed. The longer fibers (l/z > 1) exhibited a bi-modal probability distribution for the fiber velocity, which reflected an almost equal likelihood for a long fiber to reside in an ejection or in a sweep. It was also observed that in the buffer region, high-speed long fibers were almost randomly oriented whereas for all size cases the slowly moving fibers preferentially oriented in the streamwise direction. These phenomena have not been observed in DNS studies of fiber suspension flows and suggested l/z to be an essential parameter in a new generation of wall-collision models to be used in numerical studies.     

CFD–DEM simulation of the pneumatic conveying of fine particles through a horizontal slit

Fine particles play a significant role in many industrial processes. To study the dynamic behavior of fine particle and their deposition in rock fractures, the pneumatic conveying of fine particles (approximately 100 μm in diameter) through a small-scale horizontal slit (0.41 m × 0.025 m) was studied, which is useful for the sealing technology of underground gas drainage in coal mining production. The CFD–DEM method was adopted to model the gas-particle two-phase flow; the gas phase was treated as a continuum and modeled using computational fluid dynamics (CFD), particle motion and collisions were simulated using the DEM code. Then, the bulk movement of fine particles through a small-scale horizontal slit was explored numerically, and the flow patterns were further investigated by visual inspection. The simulation results indicated that stratified flow or dune flow can be observed at low gas velocities. For intermediate gas velocities, the flow patterns showed pulsation phenomena, and dune flow reappeared in the tail section. Moreover, periodic flow regimes with alternating thick and sparse stream structures were observed at a high gas velocity. The simulation results of the bulk movement of fine particles were in good agreement with the experimental findings, which were obtained by video-imaging experiments. Furthermore, the calculated pressure drop versus gas velocity profile was investigated and compared with relative experimental findings, and the results showed good agreement. Furthermore, the particle velocity vectors and voidage distribution were numerically simulated. Selected stimulation results are presented and provide a reference for the further study of fine particles.     

Turbulent stresses and particle break-up criteria in particle-laden pipe flows

Three-dimensional particle tracking velocimetry (3D-PTV) is applied to particle-laden pipe flows at Reynolds number 10,300, based on the bulk velocity and the pipe diameter. The effects of flow direction (upward or downward) and mean concentration (in the range 0.5 × 10⁻⁵–3.2 × 10⁻⁵) on the production of turbulence are assessed for inertial particles with a Stokes number equal to 2.3, based on the particle relaxation time and viscous scales. The turbulence production and the Kolmogorov constant, both measured for particle laden flows in upflow and downflow, allowed for the derivation of a break-up criterion as a function of the radial coordinate. This criterion predicts the maximum possible particle size before break-up may occur. It is shown that the maximum particle size is bigger at the pipe centerline than in the near-wall zone by more than a factor of 5. Flow direction affects the particle concentration profile, with wall peaking in downflow and core peaking in upflow. This affects both the residence time and the maximum particle size, the latter by 7%.     

Experiments with a Wire-Mesh Sensor for stratified and dispersed oil-brine pipe flow

Two-phase oil–water flow was studied in a 15 m long horizontal steel pipe, with 8.28 cm internal diameter, using mineral oil (having 830 kg/m³ density and 7.5 mPa s viscosity) and brine (1073 kg/m³ density and of 0.8 mPa s viscosity). Measurements of the holdup and of the cross-sectional phase fraction distribution were obtained for stratified flow and for highly dispersed oil–water flows, applying a capacitive Wire-Mesh Sensor specially designed for that purpose. The applicability of this measurement technique, which uses a circuit for capacitive measurements that is adapted to conductive measurements, where one of the fluids is water with high salinity (mimicking sea water), was assessed. Values for the phase fraction values were derived from the raw data obtained by the Wire-Mesh Sensor using several mixture permittivity models. Two gamma-ray densitometers allowed the accurate measurement of the holdups, which was used to validate the data acquired with the capacitive Wire-Mesh Sensor. The measured time-averaged distribution of the phase fraction over the cross-sectional area was used to investigate the details of the observed two-phase flow patterns, including the interface shape and water height. The experiments were conducted in the multiphase-flow test facility of Shell Global International B.V. in the Netherlands.     

Flow fields and packing states in the discharge flow of noncircular particles—A SIPHPM simulation

Although the discharge flow of spherical materials has been extensively explored, the effect of particle shape on discharge is still poorly understood. The present work explores the two-dimensional discharge flow fields of noncircular particles using the soft-sphere-imbedded pseudo-hard particle model method. Rectangular particles having different aspect ratios (R_a = 1, 1.5, 2–5) and regular polygonal particles having different numbers of sides (N_s = 3–8, 10) are discharged through hopper beds having different orifice widths (D_i = 40, 70.77, 99.13, 125.74, 151.13 mm). The discharge rates of differently shaped particles in different beds are consistent with Beverloo’s relation. Moreover, the flow fields are computed and evaluated to study the effects of R_a, N_s, and D_i on particle discharge. The characteristics of particle–particle connections in the discharge process are evaluated according to the temporal evolution and spatial distribution of the contact points. Additionally, the effect of the initial packing on the discharge profile is investigated. The findings help clarify the discharge of noncircular particles.     

Origin of regime transition to turbulent flow in bubble column: Orifice- and column-induced transitions

The possible events during bubble formation on an orifice were investigated using a rectangular bubble column (30 cm × 30 cm × 100 cm). The gas flow rate through a single orifice was adjusted from 0.1 dm³/min to 5.0 dm³/min covering a high flow rate regime. At the high gas flow rate, the bubble formation process was complicated by diverse events, such as wake effect, channeling, and orifice-induced turbulent flow. The detachment period could be used to discern the bubble formation steps because it was strongly affected by the above events. The bubble size distribution around the orifice was also analyzed to gain a clearer understanding of the bubble formation process. Above the rate of 3.0 dm³/min through a single orifice, the detachment period converged to a value of 25 ms irrespective of the orifice diameter. The bubble size distribution also showed little difference in this range of gas flow rate. This could be explained by the development of turbulent flow around the orifice. A 0.15 m in-diameter bubble column was tested to investigate the effect of orifice-induced turbulent flow on the regime transition in which the homogeneous flow regime is converted into the heterogeneous flow regime in the column. Obvious distinction between the orifice- and column-induced transitions was observed.     

PIV, 2D-LIF and 1D-Raman measurements of flow field,composition and temperature in premixed gas turbine flames

Several laser diagnostic measurement techniques have been applied to study the lean premixed natural gas/air flames of an industrial swirl burner. This was made possible by equipping the burner with an optical combustion chamber that was installed in the high-pressure test rig facility at the DLR Institute of Combustion Technology in Stuttgart. The burner was operated with preheated air at various operating conditions with pressures up to p = 6 bar and a maximum thermal power of P = 1 MW.The instantaneous planar flow field inside the combustor was studied with particle image velocimetry (PIV). Planar laser induced fluorescence (PLIF) of OH radicals on a single-shot basis was used to determine the shape and the location of the flame front as well as the spatial distribution of reaction products. 1D laser Raman spectroscopy was successfully applied for the measurement of the temperature and the concentration of major species under realistic gas turbine conditions.Results of the flow field analysis show the shape and the size of the main flow regimes: the inflow region, the inner and the outer recirculation zone. The highly turbulent flow field of the inner shear layer is found to be dominated by small and medium sized vortices. High RMS fluctuations of the flow velocity in the exhaust gas indicate the existence of a rotating exhaust gas swirl. From the PLIF images it is seen that the primary reactions happened in the shear layers between inflow and the recirculation zones and that the appearance of the reaction zones changed with flame parameters. The results of the multiscalar Raman measurements show a strong variation of the local mixture fraction allowing conclusions to be drawn about the premix quality. Furthermore, mixing effects of unburnt fuel and air with fully reacted combustion products are studied giving insights into the processes of the turbulence–chemistry interaction.     

Emission of organic carbon,elemental carbon and water-soluble ions from crop straw burning under flaming and smoldering conditions

Emissions from major agricultural residues were measured using a self-designed combustion system. Emission factors (EFs) of organic carbon (OC), elemental carbon (EC), and water-soluble ions (WSIs) (K⁺, NH₄⁺, Na⁺, Mg²⁺, Ca²⁺, Cl⁻, NO₃⁻, SO₄^2–) in smoke from wheat and rice straw were measured under flaming and smoldering conditions. The OC₁/TC (total carbon) was highest (45.8% flaming, 57.7% smoldering) among carbon fractions. The mean EFs for OC (EF_OC) and EC (EF_EC) were 9.2 ± 3.9 and 2.2 ± 0.7 g/kg for wheat straw and 6.4 ± 1.9 and 1.1 ± 0.3 g/kg for rice straw under flaming conditions, while they were 40.8 ± 5.6 and 5.8 ± 1.0 g/kg and 37.6 ± 6.3 and 5.0 ± 1.4 g/kg under smoldering conditions, respectively. Higher EC ratios were observed in particulate matter (PM) mass under flaming conditions. The OC and EC for the two combustion patterns were significantly correlated (p < 0.01, R = 0.95 for wheat straw; p < 0.01, R = 0.97 for rice straw), and a higher positive correlation between OC₃ and EC was observed under both combustion conditions. WSIs emitted from flaming smoke were dominated by Cl⁻ and K⁺, which contributed 3.4% and 2.4% of the PM mass for rice straw and 2.2% and 1.0% for wheat straw, respectively. The EFs of Cl⁻ and K⁺ were 0.73 ± 0.16 and 0.51 ± 0.14 g/kg for wheat straw and 0.25 ± 0.15 and 0.12 ± 0.05 g/kg for rice straw under flaming conditions, while they were 0.42 ± 0.28 and 0.12 ± 0.06 g/kg and 0.30 ± 0.27 and 0.05 ± 0.03 g/kg under smoldering conditions, respectively. Na⁺, Mg²⁺, and NH₄⁺ were vital components in PM, comprising from 0.8% (smoldering) to 3.1% (flaming) of the mass. Strong correlations of Cl⁻ with K⁺, NH₄⁺, and Na⁺ ions were observed in rice straw and the calculated diagnostic ratios of OC/EC, K⁺/Na⁺ and Cl⁻/Na⁺ could be useful to distinguishing crop straw burning from other sources of atmospheric pollution.     

Turbulent natural convection scaling in a vertical channel

Using direct numerical simulation (DNS) data, this study appraises existing scaling laws in literature for turbulent natural convection of air in a differentially heated vertical channel. The present data is validated using past DNS studies, and covers a range of Rayleigh number, Ra between 5.4 × 10⁵ and 2.0 × 10⁷. We then appraise and compare the various scaling laws proposed by Versteegh and Nieuwstadt, 1999, Hölling and Herwig, 2005, Shiri and George, 2008, George and Capp, 1979 with the profiles of the mean temperature defect, mean streamwise velocity, normal velocity fluctuations, temperature fluctuations and Reynolds shear stress. Based on the arguments of an inner (near-wall) and outer (channel centre) region, the data is found to support a minus one-third power law for the mean temperature in an overlap region. Using the inner and outer temperature profiles, an implicit heat transfer equation is obtained and we show that a correction term is non-negligible for the present Ra range when compared with explicit equations found in literature. In addition, we determined that the mean streamwise velocity and normal velocity fluctuations collapse in the inner region when using the outer velocity scale. We also find that the temperature fluctuations scale in inner coordinates, in contrast to the outer scaling behaviour reported in the past. Lastly, we show evidence of an incipient proportional relationship between friction velocity, u_τ, and the outer velocity scale, u_o, with increasing Ra.     

An investigation into flow mode transition and pressure fluctuations for fluidized dense-phase pneumatic conveying of fine powders

This paper presents the results of an ongoing investigation into the fluctuations of pressure signals due to solids–gas flows for dense-phase pneumatic conveying of fine powders. Pressure signals were obtained from pressure transducers installed along different locations of a pipeline for the fluidized dense-phase pneumatic conveying of fly ash (median particle diameter 30 μm; particle density 2300 kg/m³; loose-poured bulk density 700 kg/m³) and white powder (median particle diameter 55 μm; particle density 1600 kg/m³; loose-poured bulk density 620 kg/m³) from dilute to fluidized dense-phase. Standard deviation and Shannon entropy were employed to investigate the pressure signal fluctuations. It was found that there is an increase in the values of Shannon entropy and standard deviation for both of the products along the flow direction through the straight pipe sections. However, both the Shannon entropy and standard deviation values tend to decrease after the flow through bend(s). This result could be attributed to the deceleration of particles while flowing through the bends, resulting in dampened particle fluctuation and turbulence. Lower values of Shannon entropy in the early parts of the pipeline could be due to the non-suspension nature of flow (dense-phase), i.e., there is a higher probability that the particles are concentrated toward the bottom of pipe, compared with dilute-phase or suspension flow (high velocity), where the particles could be expected to be distributed homogenously throughout the pipe bore (as the flow is in suspension). Changes in straight-pipe pneumatic conveying characteristics along the flow direction also indicate a change in the flow regime along the flow.     

Measurements of emission spectra from hot,dense germanium plasma in short pulse laser experiments

Heating of thin foil targets by an high power laser at intensities of 10¹⁷–10¹⁹ W/cm² has been studied as a method for producing high temperature, high density samples to investigate X-ray opacity and equation of state. The targets were plastic (parylene-N) foils with a microdot made of a mixture of germanium and titanium buried at depth of 1.5 μm. The L-shell spectra from the germanium and the K-shell spectra from the titanium were taken using crystal spectrometers recording onto film and an ultra fast X-ray streak camera coupled to a conical focussing crystal with a time resolution of 1 ps. The conditions in the microdot were inferred by comparing the measured spectra to synthetic spectra produced by the time-dependent collisional–radiative (CR) models FLY and FLYCHK. The data were also compared to simulated spectra from a number of opacity codes assuming local thermodynamic equilibrium (LTE). Temperature and density gradients were taken into account in the comparisons. The sample conditions were inferred from the CR modelling using FLYCHK to be 800 ± 100 eV and 1.5 ± 0.5 g/cc. The best fit to the LTE models was at a temperature 20% lower than with the CR model. Though the sample departs from LTE significantly useful spectral comparisons can still be made. The results and comparisons are discussed along with improvements to the experimental technique to achieve conditions closer to LTE.     

Experimental-computational study of fibrous particle transport and deposition in a bifurcating lung model

Experiments carried out using a lung model with a single horizontal bifurcation under different steady inhalation conditions explored the orientation of depositing carbon fibers, and particle deposition fractions. The orientations of deposited fibers were obtained from micrographs. Specifically, the effects of the sedimentation parameter (γ), fiber length, and flow rate on orientations were analyzed. Our results indicate that gravitational effect on deposition cannot be neglected for 0.0228 < γ < 0.247. The absolute orientation angle of depositing fibers decreased linearly with increasing γ for values 0.0228 < γ < 0.15. Correspondence between Stokes numbers and γ suggests these characteristics can be used to estimate fiber deposition in the lower airways. Computer simulations with sphere-equivalent diameter models for the fibers explored deposition efficiency vs. Stokes number. Using the volume-equivalent diameter model, our experimental data for the horizontal bifurcation were replicated. Results for particle deposition using a lung model with a vertical bifurcation indicate that body position also affects deposition.     

Experimental study on the direct contact condensation of the steam jet in subcooled water flow in a rectangular channel: Flow patterns and flow field

Direct contact condensation (DCC) of steam jet in subcooled water flow in a channel was experimentally studied. The main inlet parameters, including steam mass flux, water mass flux and water temperature were tested in the ranges of 200–600 kg/(m² s), 7–18 × 10³ kg/(m² s), 288–333 K, respectively. Two unstable flow patterns and two stable flow patterns were observed via visualization window by a high speed camera. The flow patterns were determined by steam mass flux, water mass flux and water temperature, and the relationship between flow patterns and flow field parameters was discussed. The results indicated that whether pressure or temperature distributions on the bottom wall of channel could represent different flow patterns. And the position of pressure peak on the bottom wall could almost represent the condensation length. The upper wall pressure distributions were mainly dependent on steam and water mass flux; and the upper wall temperature distributions were affected by the three main inlet parameters. Moreover, the bottom wall pressure and temperature distributions of different unstable flow patterns had similar characteristics while those of stable flow patterns were affected by shock and expansion waves. The underlying cause of transition between different flow patterns under different inlet parameters was reflected and discussed based on pressure distributions.