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.     

Positive frictional pressure gradient in vertical gas-high viscosity oil slug flow

The present study describes the wall shear stress and the falling liquid film behavior in upward vertical slug flow of air and high viscosity oil. The frictional pressure gradient is directly related to the wall shear stress, and it is usually negative (opposite to the overall flow direction). However, in vertical slug flow, the average total wall shear stress of a slug unit may be negative (in the same direction of the overall flow), resulting in a positive frictional pressure gradient. However, this does not mean, by any way, generation of additional energy or violation of the second law of thermodynamics.The positive frictional pressure gradient phenomenon, reasons and required conditions were explained in this paper. A simplified model was developed and validated against recent experimental data of air-high viscosity oil slug flow in a 50.8 mm ID vertical pipe. The oil viscosity was in the range of 127 mPa s to 580 mPa s. Positive frictional pressure gradient appears when the liquid film wall shear stress supersede the wall shear stress in the slug body. The rate of increase of both wall shear stresses (with respect to the mixture Reynolds number) depend, not only, on the mixture Reynolds number but also, highly, on the liquid viscosity.     

Experimental investigation on thermal properties of silver nanofluids

This paper reports on an experimental investigation of the thermal properties behavior of 0.5 wt% silver nanoparticle-based nanofluids (NF) containing oleic acid (OA) and potassium oleate surfactant (OAK⁺) with concentrations of 0.5, 1, and 1.5 wt% respectively. The experiments were conducted from 20 °C to 80 °C. It was shown that the NF with 1 wt% OAK⁺ yielded the highest thermal behavior enhancement of about 28% at 80 °C compared to deionized water. The thermal performance had higher than the base fluid/nanofluids at approximately 80%. Moreover, the NF containing OAK⁺ showed higher thermal conductivity and dynamics of specific heat capacity than deionized water in all of the experimental conditions in this study. The rheological experiment showed that viscosity of NF was significantly dependant on temperature. As shear rate increased, the shear stress of the NF increased; however, the viscosity of the nanofluids decreased first and then stabilized. It was further found that NF containing OAK⁺ at a range of operating temperatures produced Newtonian behavior.     

An experimental and numerical study of the localization behavior of tantalum and stainless steel

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic–plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material’s local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length is 870–890 μm with deformation imposed through a split-Hopkinson pressure bar system at maximum top-to-bottom velocity in the range of 8–25 m/s. We present experimental results of the deformation response of tantalum and 316L stainless steel samples. The tantalum samples did not form shear bands but the stainless steel sample formed a late stage shear band. We have also modeled these experiments using both conductive and adiabatic continuum models. An anisotropic elasto-viscoplastic constitutive model with damage evolution was used within the finite element code EPIC. A Mie-Gruneisen equation of state and the rate and temperature sensitive MTS flow stress model together with a Gurson flow surface were employed. The models performed well in predicting the experimental data. The numerical results for tantalum suggested a maximum equivalent strain rate on the order of 7 × 10⁴ s⁻¹ in the gage section for an imposed top surface displacement rate of 17.5 m/s. The models also suggested that for an initial temperature of 298 K a temperature in the neighborhood of 900 K was reached within the shear section. The numerical results for stainless steel suggest that melting temperature was reached throughout the shear band shortly after peak load. Due to sample geometry, the stress state in the shear zone was not pure shear; a significant normal stress relative to the shear zone basis line was developed.     

Liquid jet breakup for non-circular orifices under low pressures

This study uses a high-speed visualization technique to investigate the breakup process and flow behavior of low pressure water jets issued from non-circular orifices including square, triangular, and rectangular shapes. These orifices have approximately the same sectional areas. Stability curve and Ohnesorge chart are employed to make a comparison with circular jets discharged from a circular orifice of the same sectional area. The analysis is carried out for gauge pressures varying from 0.1 psi to 70 psi with small pressure steps corresponding to a range from 0.7 kPa to 482.6 kPa in metric units. Axis-switching phenomenon is observed and analyzed through calculating the axis-switching wavelength and oscillation frequency for rectangular jets. It is found that results for circular jets agreed well with classic theory. Non-circular jets demonstrate enhanced instabilities as a whole compared to circular jets. The different behaviors of non-circular jets are reasonably explained by Rayleigh’s oscillation theory. Axis-switching and aspect-ratio effect in rectangular jets is found to slow down the increase of breakup-length in the Rayleigh breakup regime. Square and triangular jets are more susceptible to wind effects and they are more unstable especially at higher pressure conditions. This can be concluded from the shorter breakup-length and narrower transitional region from the Rayleigh regime to the wind-induced regime as compared to the circular and rectangular jets. Axis-switching wavelength of the rectangular jets is found to increase linearly with increasing jet velocity and oscillation frequency decreases correspondingly.     

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).     

Numerical evaluation of turbulence models for dense to dilute gas–solid flows in vertical conveyor

A two-fluid model (TFM) of multiphase flows based on the kinetic theory and small frictional limit boundary condition of granular flow was used to study the behavior of dense to dilute gas–solid flows in vertical pneumatic conveyor. An axisymmetric 2-dimensional, vertical pipe with 5.6 m length and 0.01 m internal diameter was chosen as the computation domain, same to that used for experimentation in the literature. The chosen particles are spherical, of diameter 1.91 mm and density 2500 kg/m³. Turbulence interaction between the gas and particle phases was investigated by Simonin's and Ahmadi's models and their numerical results were validated for dilute to dense conveying of particles. Flow regimes transition and pressure drop were predicted. Voidage and velocity profiles of each phase were calculated in radial direction at different lengths of the conveying pipe. It was found that the voidage has a minimum, and gas and solid velocities have maximum values along the center line of the conveying pipe and pressure drop has a minimum value in transition from dense slugging to dilute stable flow regime. Slug length and pressure fluctuation reduction were predicted with increasing gas velocity, too. It is shown that solid phase turbulence plays a significant role in numerical prediction of hydrodynamics of conveyor and the capability of particles turbulence models depends on tuning parameters of slip-wall boundary condition.     

Smoothed particle hydrodynamics simulation of non-Newtonian moulding flow

Smoothed particle hydrodynamics (SPH) has been widely applied in simulating fluid flow because of its attractive properties, for example, it is fully Lagrangian and mesh free. However, this method usually uses the explicit method to solve the conservation equations and in this form it is only suitable to momentum dominated flows with low viscosity. In polymer processing, the fluid is non-Newtonian with high viscosity, O(10³) to O(10⁴) Pa-s say, and the pressure is high as O(10⁶) to O(10¹⁰) Pa. The algorithm of the standard SPH is infeasible in this case, because only very small time steps can be used for a stable simulation. We have developed an implicit SPH for non-Newtonian flow, which is completely matrix free, to solve the equation system iteratively and robustly. The artificial pressure is introduced between particles to stabilize the SPH system avoiding the tensile instability. The fluid is compressible under high pressure. Realistic state equations for polymers, such as the Tait and SSY [16] equations are adopted to describe the density/pressure relations. The method is finally applied to the simulation of moulding flow of a modified power law fluid with the zero shear rate viscosity of 1.22 × 10⁴ Pa-s, Reynolds number of 3 × 10^?4 to 6 × 10^?5 and the highest pressure of O(10⁸) to O(10¹⁰) Pa.     

Buckling behavior of Kelvin open-cell foams under [0 0 1], [0 1 1] and [1 1 1] compressive loads

This paper describes buckling modes and stresses of elastic Kelvin open-cell foams subjected to [0 0 1], [0 1 1] and [1 1 1] uniaxial compressions. Cubic unit cells and cell aggregates in model foams are analyzed using a homogenization theory of the updated Lagrangian type. The analysis is performed on the assumption that the struts in foams have a non-uniform distribution of cross-sectional areas as observed experimentally. The relative density is changed to range from 0.005 to 0.05. It is thus found that long wavelength buckling and macroscopic instability primarily occur under [0 0 1] and [0 1 1] compressions, with only short wavelength buckling under [1 1 1] compression. The primary buckling stresses under the three compressions are fairly close to one another and almost satisfy the Gibson–Ashby relation established to fit experiments. By also performing the analysis based on the uniformity of strut cross-sectional areas, it is shown that the non-uniformity of cross-sectional areas is an important factor for the buckling behavior of open-cell foams.     

Extreme softness of brain matter in simple shear

We show that porcine brain matter can be modelled accurately as a very soft rubber-like material using the Mooney–Rivlin strain energy function, up to strains as high as 60%. This result followed from simple shear experiments performed on small rectangular fresh samples (2.5 cm³ and 1.1 cm³) at quasi-static strain rates. They revealed a linear shear stress–shear strain relationship $(R^2>0.97)$, characteristic of Mooney–Rivlin materials at large strains. We found that porcine brain matter is about 30 times less resistant to shear forces than a silicone gel. We also verified experimentally that brain matter exhibits the positive Poynting effect of non-linear elasticity, and numerically that the stress and strain fields remain mostly homogeneous throughout the thickness of the samples in simple shear.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

Effect of gravity on the mechanical properties of lunar regolith tested using a low gravity simulation device

Simulations of the bearing capacity and shear strength of regolith under Earth’s gravity produce different results from those under low gravity. A low-gravity simulation device was developed in this study, and an internal stress model of regolith simulant was established to correct the errors. The model revealed additional force on both shear plane in the shear test and the press plate area in the pressure–sinkage test. The sinkage and shear test results showed that low gravity decreased the deformable index n, frictional modulus k_φ and cohesion c, whereas there were no obvious changes to the cohesive modulus k_c and internal friction angle φ. The sinkage generally increased as the gravity decreased under a consistent normal load larger than 50 N, but when the wheel load was lower than 50 N, the sinkage of the TYII-1 simulant was larger under 1 G than 1/6 G. Gravity had little effect on the shear strength of the regolith. However, the tractive thrust of the TYII-1 simulant was lower under 1/6 G than 1 G. The smaller difference was due to differences in the way the soils responded to changes in the gravity level for the TYII-2 simulant.     

Measurement method of complex viscoelastic material properties

A measurement technique of viscoelastic properties of polymers is proposed to investigate complex Poisson’s ratio as a function of frequency. The forced vibration responses for the samples under normal and shear deformation are measured with varying load masses. To obtain modulus of elasticity and shear modulus, the present method requires only knowledge of the load mass, geometrical characteristics of a sample, as well as both the amplitude ratio and phase lag of the forcing and response oscillations. The measured data were used to obtain the viscoelastic properties of the material based on a 2D numerical deformation model of the sample. The 2D model enabled us to exclude data correction by the empirical form factor used in 1D model. Standard composition (90% PDMS polymer + 10% catalyst) of silicone RTV rubber (Silastic^® S2) were used for preparing three samples for axial stress deformation and three samples for shear deformation. Comprehensive measurements of modulus of elasticity, shear modulus, loss factor, and both real and imaginary parts of Poisson’s ratio were determined for frequencies from 50 to 320 Hz in the linear deformation regime (at relative deformations 10^?6 to 10^?4) at temperature 25 °C. In order to improve measurement accuracy, an extrapolation of the obtained results to zero load mass was suggested. For this purpose measurements with several masses need to be done. An empirical requirement for the sample height-to-radius ratio to be more than 4 was found for stress measurements. Different combinations of the samples with different sizes for the shear and stress measurements exhibited similar results. The proposed method allows one to measure imaginary part of the Poisson’s ratio, which appeared to be about 0.04–0.06 for the material of the present study.     

Irreversible thermodynamics of triple junctions during the intergranular void motion under the electromigration forces

A rigorous reformulation of internal entropy production and the rate of entropy flow is developed for multi-component systems consisting of heterophases, interfaces and/or surfaces. The result is a well-posed moving boundary value problem describing the dynamics of curved interfaces and surfaces associated with voids and/or cracks that are intersected by grain boundaries. Extensive computer simulations are performed for void configuration evolution during intergranular motion. In particular we simulate evolution resulting from the action of capillary and electromigration forces in thin film metallic interconnects having a “bamboo” structure, characterized by grain boundaries aligned perpendicular to the free surface of the metallic film interconnects. Analysis of experimental data utilizing previously derived mean time to failure formulas gives consistent values for interface diffusion coefficients and enthalpies of voids. 3.0 × 10⁻⁶ exp(−0.62 eV/kT) m² s⁻¹ is the value obtained for voids that form in the interior of the aluminum interconnects without surface contamination. 6.5 × 10⁻⁶ exp(−0.84 eV/kT) m² s⁻¹ is obtained for those voids that nucleate either at triple junctions or at the grain boundary-technical surface intersections, where the chemical impurities may act as trap centers for hopping vacancies.     

Set up of an experimental apparatus for the study of fragmentation of solid fuels upon severe heating

An experimental apparatus has been developed in order to perform tests of primary fragmentation of solid fuels under severe heating conditions. The device is a modified heated strip reactor, capable to reach 2000 °C in less than 0.2 s. Particles are laid on the strip and pyrolysed under inert or moderately oxidizing conditions. The char particles and their fragments, generated upon pyrolysis, can be recovered and analysed to assess the fragmentation propensity of the fuel.Some preliminary experiments have been carried out on two biomass samples in order to assess the time-temperature history of particles in the experimental apparatus. In particular biomass particles of approximately 2–3 mm have been used. The temperature of the heated strip reactor in such preliminary tests was varied between 1000 and 1600 °C, while the strip nominal heating rate was kept at 10⁴ °C/s and the holding time was set at the value of 10 s. A near infrared fast camera (38,000 frames/s) has been used to measure the temperature of the heated strip and of the particles during the tests. A heat up model was developed and validated against experimental results. The model was then used to estimate the temperature gradients across particles of biomass and of coal as well.Results show that the strip of the reactor reaches the set temperature in less than 0.2 s. When particles are laid on the strip, their bottom surface, which is in physical contact with the strip, immediately reaches the set temperature value. For 1 mm coal particles the upper surface can be considered at the same temperature as well. Under the most severe conditions tested (strip temperature of 1600 °C , biomass particles of 2 mm thickness) the temperature difference between the bottom and the upper face is 200 °C after 3 s and drops to 100 °C after 10 s. On the whole the experimental apparatus simulates uniform heating of the particles with reasonable approximation. In the next future the apparatus will be further upgraded to operate at pressures up to 20 bars.     

Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal 1.030 mm circular channel

This research focuses on acquiring accurate flow boiling heat transfer data and flow pattern visualization for three refrigerants, R134a, R236fa and R245fa in a 1.030 mm channel. We investigate trends in the data, and their possible mechanisms, for mass fluxes from 200 to 1600 kg/m²s, heat fluxes from 2.3 kW/m² to 250 kW/m² at T_sat = 31 °C and ΔT_sub from 2 to 9 K. The local saturated flow boiling heat transfer coefficients display a heat flux and a mass flux dependency but no residual subcooling influence. The changes in heat transfer trends correspond well with flow regime transitions. These were segregated into the isolated bubble (IB) regime, the coalescing bubble (CB) regime, and the annular (A) regime for the three fluids. The importance of nucleate boiling and forced convection in these small channels is still relatively unclear and requires further research.     

Hydrothermal synthesis of ultrathin hexagonal nickel hydroxide nanosheets

Nickel hydroxide, Ni(OH)₂ is widely used in electrodes of nickel-based alkaline secondary batteries. Ultrathin hexagonal Ni(OH)₂ nanosheets of space group P-3m1 were hydrothermally prepared at 200 °C for 10 h. Their diameter and thickness were 200–300 and 3–5 nm, respectively. Their formation was attributed to the oriented assembly of growing particles, which was assisted by surfactant molecules. The specific surface area of the Ni(OH)₂ nanosheets was 8.66 m²/g. Their magnetization curve exhibited linear paramagnetic behavior across the entire measurement region.     

Influences of ions and temperature on performance of carbon nano-particulates in supercapacitors with neutral aqueous electrolytes

A commercial product of carbon nano-particles, Cabot MONACH 1300 pigment black (CMPB), was studied for basic structural information and electrochemical performance in neutral aqueous electrolytes, aiming at applications in supercapacitors. As confirmed by SEM and HRTEM, the CMPB had a hierarchical structure, containing basic 10 nm nano-spheres which combined into ca. 50 nm agglomerates which further aggregated into larger particles of micrometres. The capacitance of this commercial material was found to increase with decreasing the size of hydrous cation (Li⁺ → Na⁺ → K⁺), instead of the cation crystal radius (K⁺ → Na⁺ → Li⁺) when coupled with the same anion (Cl⁻). In electrolytes with the same cation concentration (K⁺), changing the anion from the larger dianion (SO₄²⁻) to the smaller monoanion (Cl⁻) also increased the capacitance at high potential scan rates (>50 mV/s). Increasing electrolyte concentration produced expected effect, including raising the electrode capacitance, but lowering the equivalent series resistance (ESR), charge transfer resistance (CTR), and the diffusion resistance. At higher temperatures, the CMPB exhibited slightly higher capacitance, which does not agree with the Gouy–Chapman theory on electric double layer (EDL). A hypothesis is proposed to account for the capacitance increase with temperature as a result of the CMPB opening up some micro-pores for more ions to access in response to the temperature increase.     

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.     

Large Eddy Simulations of turbulent particle-laden channel flow

This paper scrutinises the Large Eddy Simulation (LES) approach to simulate the behaviour of inter-acting particles in a turbulent channel flow. A series of simulations that are fully (four-way), two-way and one-way coupled are performed in order to investigate the importance of the individual physical phenomena occurring in particle-laden flows. Moreover, the soft sphere and hard sphere models, which describe the interaction between colliding particles, are compared with each other and the drawbacks and advantages of each algorithm are discussed. Different models to describe the sub-grid scale stresses with LES are compared. Finally, simulations accounting for the rough walls of the channel are compared to simulations with smooth walls. The results of the simulations are discussed with the aid of the experimental data of Kussin J. and Sommerfeld M., 2002, Experimental studies on particle behaviour and turbulence modification in horizontal channel flow with different wall roughness, Exp. in Fluids, 33, pp. 143–159 of Reynolds number 42,000 based on the full channel height. The simulations are carried out in a three-dimensional domain of 0.175 m × 0.035 m  × 0.035 m where the direction of gravity is perpendicular to the flow. The simulation results demonstrate that rough walls and inter-particle collisions have an important effect in redistributing the particles across the channel, even for very dilute flows. A new roughness model is proposed which takes into account the fact that a collision in the soft sphere model is fully resolved and it is shown that the new model is in very good agreement with the available experimental data.