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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Experimental study on steam chugging phenomenon in a vertical sparger

Unstable direct contact condensation called ``Chugging'' that occurs in certain conditions in the pressure suppression pool of Primary Containment Vessel of Boiling Water Reactors (BWRs) was studied experimentally. The mechanisms of every phase of the chugging was described, and experimental results useful for the development and validation of more accurate CFD models were provided. The experiment was conducted with a transparent pool and a transparent polycarbonate pipe or a stainless steel pipe with inner diameter of 27 mm under the conditions of the steam mass flux of 5.5–19.5 kg/m²s and the pool temperature of 19–46.5 °C. Pressure pulses were measured and synchronized with a high speed video camera for images acquisition. It was identified that the bubble implosion occurred while the pressure in the bubble quickly decreased. This condition might establish instability in the interfacial area which grew abruptly causing the implosion. Moreover the transparent apparatus allowed to interpret and relate internal condensations, generating pressure spikes of around 1.2 MPa because of the condensation-induced water hammer. Finally, the chugging condensation regime map was created from the experimental data.     

Experimental investigation of natural convection in a vertical rib-roughened channel with asymmetric heating

An experimental study in an open-ended vertical channel is carried out in order to describe the fluid dynamics and heat transfer of transient free convection inside a vertical rib-roughened channel asymmetrically heated at various uniform heat fluxes (650, 700, and 780 W/m²) corresponding to various modified Rayleigh numbers (3.65 × 10⁶, 3.93 × 10⁶ and 4.4 × 10⁶). Two ribs are symmetrically located on each wall. The investigations focused more specifically on the influence of the ribs positions inside the channel and the modified Rayleigh number (Ra^*) both in steady-state regime and during the transitional phase occurring just after the start of the heating on the flow structure and the heat transfer performance. The results showed the appearance of large-scale flow instabilities which will develop and propagate until the development of the pocket-like vortex (reversed flow). Also, the formation and breakup of recirculation eddies, vortex banishment, besides that a separation and shifting of the boundary layer from one wall to another are identified. The best position of the ribs for heat extraction depends on the magnitude of the Rayleigh number. In that case, the top position is the optimal position for the small and the moderate modified Rayleigh numbers.     

An investigation of a model of the flow pattern transition mechanism in relation to the identification of annular flow of R134a in a vertical tube using various void fraction models and flow regime maps

In the present study, new experimental data are presented for literature on the prediction of film thickness and identification of flow regime during the co-current downward condensation in a vertical smooth copper tube having an inner diameter of 8.1 mm and a length of 500 mm. R134a and water are used as working fluids in the tube side and annular side of a double tube heat exchanger, respectively. Condensation experiments are done at mass fluxes of 300 and 515 kg m^?2 s^?1. The condensing temperatures are between 40 and 50 °C; heat fluxes are between 12.65 and 66.61 kW m^?2. The average experimental heat transfer coefficient of the refrigerant HFC-134a is calculated by applying an energy balance based on the energy transferred from the test section. A mathematical model by Barnea et al. based on the momentum balance of liquid and vapor phases is used to determine the condensation film thickness of R134a. The comparative film thickness values are determined indirectly using relevant measured data together with various void fraction models and correlations reported in the open literature. The effects of heat flux, mass flux, and condensation temperature on the film thickness and condensation heat transfer coefficient are also discussed for the laminar and turbulent flow conditions. There is a good agreement between the film thickness results obtained from the theoretical model and those obtained from six of 35 void fraction models in the high mass flux region of R134a. In spite of their different valid conditions, six well-known flow regime maps from the literature are found to be predictive for the annular flow conditions in the test tube in spite of their different operating conditions.     

New contributions on laminar flow of inelastic non-Newtonian fluid in the two-dimensional symmetric expansion: Creeping and slowly moving flow conditions

The steady flow of generalized Newtonian fluid in a two-dimensional 1:3 sudden expansion was studied numerically. Finite volume method was applied to solve the momentum equations along with the continuity equation and the Power law rheological model within the laminar flow regime for a range of Reynolds number and Power law index values. The values of generalized Reynolds number, based on physical and rheological properties, upstream channel height and bulk velocity, were varied between 0.0001 ≤ Re_gen ≤ 10, while the Power law index values mapped the 0.60 ≤ n ≤ 1.40 range, allowing for the investigation of both shear-thinning and shear-thickening effects at creeping as well as slowly moving fluid flow conditions. We report accurate results of a systematic study with a focus on most important characteristics of recirculating fluid flow in the downstream section of sudden expansion geometry. It is shown that for the creeping flow regime there exist finite sized redevelopment length, extra pressure drop (Couette correction) and recirculation zones (also called as Moffatt vortices) that are influenced by the non-Newtonian viscous behaviour.