Drag reduction for liquid flow through micro-apertures

Pressure drops in the flow through micro-orifices and capillaries were measured for silicone oils, aqueous solutions of polyethylene glycol (PEG), and surfactant aqueous solutions. The diameter of micro-orifices ranged from 5 μm to 400 μm. The corresponding length/diameter ratio was from 4 to 0.05 and capillary diameters were 105 μm and 450 μm. The following results were obtained: silicone oils of 10^?6 m²/s and 10^?5 m²/s in kinematic viscosity generated a reduction of pressure drop (RPD), that is, drag reduction, similar to the RPD of water and a glycerol/water mixture reported in the previous paper by the present authors. When RPD occurred, the pressure drop (PD) of silicone oils of 10^?6 m²/s and 10^?5 m²/s had nearly the same magnitude. Namely, the difference in viscosity did not influence RPD. A 10³ ppm aqueous solution of PEG20000 provided almost the same PD as that of PEG8000 for the 400 μm to 15 μm orifices, but a greater PD than that of PEG8000 for the 10 μm to 5 μm orifices. A non-ionic surfactant and a cationic surfactant were highly effective in RPD compared with anionic surfactants: the non-ionic and cationic surfactant solutions had PD one order of magnitude lower than that of water under some flow conditions in the concentration range from 1 ppm to 10⁴ ppm, but the anionic surfactant solutions did not generate RPD except in the case of the smallest orifice of 5 μm in diameter. The PD of the non-ionic surfactant solution showed a steep rise at a Reynolds number (Re_t) for 400 μm to 15 μm orifices. The Re_t provides the relationship Re_t = K/D, where D is the orifice diameter, and K is a constant of 2 × 10^?2 m for the 100–20 μm orifices irrespective of liquid concentration. Capillary flow experiment revealed that the PEG, non-ionic and cationic surfactant solutions generated RPD also in a laminar flow through the capillary of 105 μm in diameter, but not in the flow through the capillary of 450 μm in diameter. In order to clarify the cause of RPD, an additional experiment was carried out by changing the orifice material from metal to acrylic resin. The result gave a different appearance of RPD, suggesting that RPD is related to an interfacial phenomenon between the liquid and wall. The large RPDs found in the present experiment are very interesting from both academic and practical viewpoints.     

Graphite dust resuspension in an HTR-10 steam generator

Graphite dust has an important effect on the safety of high-temperature gas-cooled reactors (HTR). The flow field in the steam generator was studied by the computational fluid dynamics (CFD) method, with the results indicating that the friction velocity in the windward and the leeward of the heat transfer tubes is relatively low and is higher at the sides. Further analysis of the resuspension of graphite dust indicates that the resuspension fraction reaches nearly zero for particles with a diameter less than 1 μm, whereas it will increases as the helium velocity in the steam generator increases for particle size larger than 1 μm. Moreover, the resuspension fraction increases as the particle size increases. The results also indicate that resuspension of the particles with sizes larger than 1 μm exhibited obvious differences in different parts of the steam generator.     

Removal of inhalable particles from coal and refuse combustion by agglomeration with solid nuclei

Airborne inhalable particles are a potent environmental pollutant. Formed via industrial processes, separation of these particles is difficult using conventional clean up techniques. In this work, solid nuclei particles of different chemical compositions were introduced into an agglomeration chamber with simulated flue gases to investigate their ability to remove these particles. Organic nuclei were able to capture more inhalable particles from coal-derived fly ash than inorganic nuclei, though these proved more effective for the agglomeration of inhalable particles in refuse-derived fly ash. Increasing the diameter of the solid nuclei benefitted the agglomeration process for both types of ash. Varying the local humidity changed adhesion between the particles and encouraged them to aggregate. Increasing the relative humidity consistently increased particle agglomeration for the refuse-derived ash. For coal-derived fly ash, the removal efficiency increased initially with relative humidity but then further increases in humidity had no impact on the relatively high efficiencies. After agglomeration in an atmosphere of 62% relative humidity, the mean mass diameter of inhalable particles in the coal-derived fly ash increased from 3.3 to 9.2 μm. For refuse-derived fly ash, agglomeration caused the percentage of particles that were less than 2 μm to decrease from 40% to 15%. After treatment at a relative humidity of 61%, the mean size of inhalable particles exceeded 10 μm.     

Experimental investigation of turbulence modulation in particle-laden coaxial jets by Phase Doppler Anemometry

The effect of solid particles on the flow characteristics of axisymmetric turbulent coaxial jets for two flow conditions was studied. Simultaneous measurements of size and velocity distributions of continuous and dispersed phases in a two-phase flow are presented using a Phase Doppler Anemometry (PDA) technique. Spherical glass particles with a particle diameter range from 102 to 212 μm were used in this two-phase flow, the experimental results indicate a significant influence of the solid particles and the Re on the flow characteristics. The data show that the gas phase has lower mean velocity in the near-injector region and a higher mean velocity at the developed region. Near the injector at low Reynolds number (Re = 2839) the presence of the particles dampens the gas-phase turbulence, while at higher Reynolds number (Re = 11 893) the gas-phase turbulence and the velocity fluctuation of particle-laden jets are increased. The particle velocity at higher Reynolds number (Re = 11 893) and is lower at lower Reynolds number (Re = 2839). The slip velocity between particles and gas phase existed over the flow domain was examined. More importantly, the present experiment results suggest that, consideration of the gas characteristic length scales is insufficient to predict gas-phase turbulence modulation in gas-particle flows.     

Reynolds number dependence of gas-phase turbulence in gas–particle flows

A downward flow of glass bead particles in a vertical pipe is investigated using a two-component LDV/PDPA for a range of Re (6400 < Re < 24,000) and a constant particle loading (m = 0.7). Two particle sizes of 70 and 200 μm are considered in the present work. For the 70 μm particles, the presence of the particles dampens the gas-phase turbulence intensity at the lowest value of Re investigated (8300) compared with the single-phase flow at the same Re. As Re increases, the gas turbulence increases, and for Re > 13,800 the gas turbulence is enhanced compared with the single-phase flow at the same Re. For the 200 μm particles, the intensity also increases with Re and is enhanced for all values of Re investigated, except at the lowest value of Re investigated (6400). At this value, the gas turbulence is equal to that of single-phase flow at the same Re. The observed trend in the gas-phase turbulence modulation with Re is proposed to be due to the change in the segregation patterns and in the average volume fractions of the particles with increasing Re. More importantly, the present experimental results suggest that, consideration of either the gas and particle characteristic length scales or the particle Reynolds number solely is insufficient to predict gas-phase turbulence modulation in gas–particle flows.     

Characteristics of turbulent particle transport in human airways under steady and cyclic flows

Motion of monodispersed aerosol particles suspended in air flow has been studied on realistic transparent model of human airways using Phase Doppler Particle Analyser (P/DPA). Time-resolved velocity data for particles in size range 1–8 μm were processed using Fuzzy Slotting Technique to estimate the power spectral density (PSD) of velocity fluctuations. The optimum processing setup for our data was found and recommendations for future experiments to improve PSD quality were suggested. Typical PSD plots at mainstream positions of the trachea and the upper bronchi are documented and differences among (1) steady-flow regimes and equivalent cyclic breathing regimes, (2) inspiration and expiration breathing phase and (3) behaviour of particles of different sizes are described in several positions of the airway model. Systematically higher level of velocity fluctuations in the upper part of the frequency range (30–500 Hz) was found for cyclic flows in comparison with corresponding steady flows. Expiratory flows in both the steady and cyclic cases produce more high-frequency fluctuations compared to inspiratory flows. Negligible differences were found for flow of particles in the inspected size range 1–8 μm at frequencies below 500 Hz. This finding was explained by Stokes number analysis. Implied match of the air and particle flows thereby indicates turbulent diffusion as important deposition mechanism and confirms the capability to use the P/DPA data as the air flow velocity estimate.     

Effect of particle size distribution on pressure drop and concentration profile in pipeline flow of highly concentrated slurry

The experiments were conducted in 54.9 mm diameter horizontal pipe on two sizes of glass beads of which mean diameter and geometric standard deviation are 440 μm & 1.2 and 125 μm & 1.15, respectively, and a mixture of the two sizes in equal fraction by mass. Flow velocity was up to 5 m/s and overall concentration up to 50% by volume for each velocity. Pressure drop and concentration profiles were measured. The profiles were obtained traversing isokinetic sampling probes in the horizontal, 45° inclined and vertical planes including the pipe axis. Slurry samples of the mixture collected in the vertical plane were analyzed for concentration profiles of each particle batch constituting the mixture. It was found that the pressure drop is decreased for the mixture at high concentrations except 5 m/s and a distinct change of concentration profiles was observed for 440 μm particles indicating a sliding bed regime, while the profiles in the horizontal plane remains almost constant irrespective of flow velocity, overall concentration and slurry type.     

A solids mixing rate correlation for small scale fluidized beds

A new first degree solids mixing rate is proposed to evaluate the mixing of solids in small scale fluidized beds. Particle mixing experiments were carried out in a 2D fluidized bed with a cross-section of 0.02 m × 0.2 m and a height of 1 m. White and black particles with average diameters of 850 and 450 μm were used in our experiments. Image processing was used to measure the concentration of the tracers at different times. The effects of four representative operating parameters (superficial gas velocity, ratio of tracer particles to bed particles, tracer particle position, and particle size) on mixing are discussed with reference to the mixing index. We found that the Lacey index depends on the concentration of the tracers. The position of the tracers affects the initial mixing rate but not the final degree of mixing. However, the new mixing rate equation does not depend on the initial configuration of the particles because this situation is considered to be the initial condition. Using the data obtained in this work and that found in literature, an empirical correlation is proposed to evaluate the mixing rate constant as a function of dimensionless numbers (Archimedes, Reynolds, and Froude) in small scale fluidized beds. This correlation allows for an estimation of the mixing rate under different operating conditions and for the detection of the end point and/or the time of mixing.     

Investigation of slip boundary conditions in the T-shaped microchannel

Flow regimes and mixing performance in a T-type micromixer at high Reynolds numbers were studied by numerical solution of the Navier–Stokes equations. The Reynolds number was varied from 1 to 1000. The cross section of the mixing channel was 100 μm × 200 μm, and its length was 1400 μm. The transverse inlet channels were symmetric about the mixing channel, and their cross-section was 100 μm × 100 μm, and the total length was 800 μm. Five different flow regimes were identified: (i) stationary vortex-free flow (Re < 5); (ii) stationary symmetric vortex flow with two horseshoe vortices (5 < Re < 150); (iii) stationary asymmetric vortex flow (150 < Re < 240); (iv) non-stationary periodic flow (240 < Re < 400); and (v) stochastic flow (Re > 400). Maximum mixing efficiency is obtained for nonstationary asymmetric vortex flow. In this case, an S-shaped vortex structure is formed in the flow field. The effect of the slip conditions on the flow pattern and mixing efficiency is studied. The slip length varied from 1 to 70 μm in the calculations. It was shown that the mixing can be controlled by hydrophobic coating.     

Preparation of SiO2/TiO2 Janus particles by electrostatic assembly,hydrolysis and calcination

This paper describes a novel chemical method for preparing SiO₂/TiO₂ Janus particles. First, polystyrene (PSt)/SiO₂ particles with a raspberry-like structure are prepared by electrostatic assembly. The influences of the reaction time of sulfonation and the treatment times of polyelectrolyte solutions (PDADMAC, PSS) on PSt and SiO₂ are investigated with respect to the surface charge density of the particles. SiO₂/TiO₂ Janus particles are then obtained by hydrolysis of butyl titanate on the surface of PSt/SiO₂ particles followed by a calcination process. Particle size analyzer, Zeta potential instrument, FTIR, TEM and SEM are used to characterize the particle size, the amount of charge on the surface of PSt and SiO₂ particles and the compositions and morphologies of PSt/SiO₂, SiO₂/TiO₂ and PSt/SiO₂/TiO₂. The diameters of the PSt, SiO₂, PSt/SiO₂ and SiO₂/TiO₂ particles are 2.0 μm, 303 nm, 2.7 μm and 330 nm, respectively.     

Technique for quantitative mapping of three-dimensional liquid–gas phase boundaries in microchannel flows

A diagnostic technique capable of characterizing interfaces between transparent, immiscible fluids is developed and demonstrated by investigating the morphology of liquid–gas interfaces in an adiabatic two-phase flow through a microchannel of 500 μm × 500 μm square cross section. Water seeded with 0.5 μm-diameter fluorescent polystyrene particles is pumped through the channel, and the desired adiabatic two-phase flow regime is achieved through controlled air injection. The diagnostic technique relies on obtaining particle position data through epifluorescent imaging of the flow at excitation and emission wavelengths of 532 nm and 620 nm, respectively. The particle position data are then used to resolve interface locations to within ±1 μm in the focal plane. By mapping the interface within individual focal planes at various depths within the channel, it is possible to determine the complete liquid–gas interface geometry across the channel cross section in a dynamic flow environment. Utilizing this approach, the liquid–gas phase boundaries of annular flows within a microchannel have been successfully characterized.     

Optical characterization of nano-sized organic carbon particles emitted from a small gasoline engine

The nano-sized organic carbon (NOC) particles emitted from a small gasoline engine were characterized using various ex situ optical techniques to assess their hazardous impact. The exhaust gas was sampled iso-kinetically by a quartz probe and passed through de-ionized water to gather the hydrophilic carbonaceous particulates as hydrosol. The hydrodynamic diameter of the particles ranged between 1.7 and 3.6 nm at no load, with a mean diameter of 2.4 nm. The particle size in the engine exhaust was found to increase at higher loads, which is attributed to coagulation of the particles. The chemical structure of the particles was analyzed using UV–vis and infra-red spectroscopy. Both the band gap energy and oscillator strength data evaluated from the UV–vis absorbance showed that the NOC particles contained polyaromatic hydrocarbon structures with three to five aromatic rings. Infra-red spectroscopy analysis further confirmed the presence of aliphatic and carbonyl functionalities in the aromatic structures of the particles. The fine size of the particles, their high number concentration for the type of the engine under study and their structural features, make the particles extremely hazardous for environment and health.     

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.     

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.     

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.     

Electrostatic beneficiation of fly ash in a free-falling system

A systematic study of fly ash electrostatic beneficiation in a free-falling separation system was carried out to provide fundamental understanding of the separation efficiency for the design of a suitable process for industrial applications. The parameters investigated included feeding position, electric field strength, particle size and moisture content. Particles larger than 105 μm presented the best separation efficiency among four different size fractions, whereas particles smaller than 44 μm showed minimal separation. However, sonication treatments helped separation by liberating more carbon from ash particles, although particle sizes were reduced as well. Experiments also showed that exposure to moisture significantly altered charging behavior of fly ash and its subsequent separation due to more free mobile ion-induced charge exchanges. The optimal feeding position was found to be slightly on the side of the negative electrode, leading to a 30% reduction in loss-on-ignition (LOI) and a 45% recovery in a single pass. A simplified mechanical model based on trajectory analysis for charged particles in an electrical field was in reasonable agreement with experimental results.     

Flow of mono-dispersed particles through horizontal bend

Pipeline slurry flow of mono-dispersed particles through horizontal bend is numerically simulated by implementing Eulerian two-phase model in FLUENT software. A hexagonal shape and Cooper type non-uniform three-dimensional grid is chosen to discretize the entire computational domain, and a control volume finite difference method is used to solve the governing equations. The modeling results are compared with the experimental data collected in 53.0 mm diameter horizontal bend with radius of 148.4 mm for concentration profiles and pressure drops. Experiments are performed on narrow-sized silica sand with mean diameter of 450 μm and for flow velocity up to 3.56 m/s (namely, 1.78, 2.67 and 3.56 m/s) and four efflux concentrations up to 16.28% (namely, 0%, 3.94%, 8.82% and 16.28%) by volume for each velocity. Eulerian model gives fairly accurate predictions for both the pressure drop and concentration profiles at all efflux concentrations and flow velocities.     

Saturated critical heat flux in a multi-microchannel heat sink fed by a split flow system

An extensive experimental campaign has been carried out for the measurement of saturated critical heat flux in a multi-microchannel copper heat sink. The heat sink was formed by 29 parallel channels that were 199 μm wide and 756 μm deep. In order to increase the critical heat flux and reduce the two-phase pressure drop, a split flow system was implemented with one central inlet at the middle of the channels and two outlets at either end. The base critical heat flux was measured using three HFC Refrigerants (R134a, R236fa and R245fa) for mass fluxes ranging from 250 to 1500 kg/m² s, inlet subcoolings from ?25 to ?5 K and saturation temperatures from 20 to 50 °C. The parametric effects of mass velocity, saturation temperature and inlet subcooling were investigated. The analysis showed that significantly higher CHF was obtainable with the split flow system (one inlet–two outlets) compared to the single inlet–single outlet system, providing also a much lower pressure drop. Notably several existing predictive methods matched the experimental data quite well and quantitatively predicted the benefit of higher CHF of the split flow.     

Experimental investigation into transient pressure pulses during pneumatic conveying of fine powders using Shannon entropy

This paper presents the results of an ongoing investigation into transient pressure pulses using Shannon entropy. Pressure fluctuations (produced by gas–solid two-phase flow during fluidized dense-phase conveying) are recorded by pressure transducers installed at strategic locations along a pipeline. This work validates previous work on identifying the flow mode from pressure signals (Mittal, Mallick, & Wypych, 2014). Two different powders, namely fly ash (median particle diameter 45 μm, particle density 1950 kg/m³, loosely poured bulk density 950 kg/m³) and cement (median particle diameter 15 μm, particle density 3060 kg/m³, loosely poured bulk density 1070 kg/m³), are conveyed through different pipelines (51 mm I.D. × 70 m length and 63 mm I.D. × 24 m length). The transient nature of pressure fluctuations (instead of steady-state behavior) is considered in investigating flow characteristics. Shannon entropy is found to increase along straight pipe sections for both solids and both pipelines. However, Shannon entropy decreases after a bend. A comparison of Shannon entropy among different ranges of superficial air velocity reveals that high Shannon entropy corresponds to very low velocities (i.e. 3–5 m/s) and very high velocities (i.e. 11–14 m/s) while low Shannon entropy corresponds to mid-range velocities (i.e. 6–8 m/s).     

Discrete element method simulations of Mars Exploration Rover wheel performance

Mars Exploration Rovers (MERs) experienced mobility problems during traverses. Three-dimensional discrete element method (DEM) simulations of MER wheel mobility tests for wheel slips of i = 0, 0.1, 0.3, 0.5, 0.7, 0.9, and 0.99 were done to examine high wheel slip mobility to improve the ARTEMIS MER traverse planning tool. Simulations of wheel drawbar pull and sinkage MIT data for i ⩽ 0.5 were used to determine DEM particle packing density (0.62) and contact friction (0.8) to represent the simulant used in mobility tests. The DEM simulations are in good agreement with MIT data for i = 0.5 and 0.7, with reasonable but less agreement at lower wheel slip. Three mobility stages include low slip (i < 0.3) controlled by soil strength, intermediate slip (i ∼ 0.3–0.6) controlled by residual soil strength, and high slip (i > 0.6) controlled by residual soil strength and wheel sinkage depth. Equilibrium sinkage occurred for i < 0.9, but continuously increased for i = 0.99. Improved DEM simulation accuracy of low-slip mobility can be achieved using polyhedral particles, rather than tri-sphere particles, to represent soil. The DEM simulations of MER wheel mobility can improve ARTEMIS accuracy.