Evaporative liquid jets in gas–liquid–solid flow system

Some aspects of the fundamental characteristics of evaporative liquid jets in gas–liquid–solid flows are studied and some pertinent literature is reviewed. Specifically, two conditions for the solids concentration in the flow are considered, including the dilute phase condition as in pneumatic convey and the dense phase condition as in bubbling or turbulent fluidized beds. Comparisons of the fundamental behavior are made of the gas–solid flow with dispersed non-evaporative as well as with evaporative liquids.For dilute phase conditions, experiments and analyses are conducted to examine the individual phase motion and boundaries of the evaporative region and the jet. Effects of the solids loading and heat capacity, system temperature, gas flow velocity and liquid injection angle on the jet behavior in gas and gas–solid flows are discussed. For dense phase conditions, experiments are conducted to examine the minimum fluidization velocity and solids distribution across the bed under various gases and liquid flow velocities. The electric capacitance tomography is developed for the first time for three-phase real time imaging of the dense gas–solid flow with evaporative liquid jets. The images reflect significantly varied bubbling phenomenon compared to those in gas–solid fluidized beds without evaporative liquid jets.     

Measurement of real-time flow structures in gas–liquid and gas–liquid–solid flow systems using electrical capacitance tomography (ECT)

The real-time cross-sectional distributions of the gas holdups in gas–liquid and gas–liquid–solid systems are measured using electrical capacitance tomography. For the gas–liquid system, air as the gas phase and both Norpar 15 (paraffin) and Paratherm as the liquid phases are used. Polystyrene beads whose permittivity is similar to that of Paratherm are used as the solid phase in the gas–liquid–solid system. The three-phase system is essentially a dielectrically two-phase system enabling measurement of the gas holdup in the gas–liquid–solid system independent of the other two phases. A new reconstruction algorithm based on a modified Hopfield dynamic neural network optimization technique developed by the authors is used to reconstruct the tomographic data to obtain the cross-sectional distribution of the gas holdup. The real-time flow structure and bubbles flow behavior in the two- and three-phase systems are discussed along with the effects of the gas velocity and the solid particles.     

Discrete particle simulation of solids motion in a gas–solid fluidized bed

The solids motion in a gas–solid fluidized bed was investigated via discrete particle simulation. The motion of individual particles in a uniform particle system and a binary particle system was monitored by the solution of the Newton's second law of motion. The force acting on each particle consists of the contact force between particles and the force exerted by the surrounding fluid. The contact force is modeled by using the analogy of spring, dash-pot and friction slider. The flow field of gas was predicted by the Navier–Stokes equation. The solids distribution is non-uniform in the bed, which is very diluted near the center but high near the wall. It was also found that there is a single solids circulation cell in the fluidized bed with ascending at the center and descending near the wall. This finding agrees with the experimental results obtained by Moslemian. The effects of the operating conditions, such as superficial gas velocity, particle size, and column size on the solids movement, were investigated. In the fluidized bed containing uniform particles better solids mixing was found in the larger bed containing smaller size particles and operated at higher superficial gas velocity. In the system containing binary particles, it was shown that under suitable conditions the particles in a fluidized bed could be made mixable or non-mixable depending on the ratios of particle sizes and densities. Better mixing of binary particles was found in the system containing particles with less different densities and closer sizes. These results were found to follow the mixing and segregation criteria obtained experimentally by Tanaka et al.     

Development and verification of coarse-grain CFD-DEM for nonspherical particles in a gas–solid fluidized bed

Computational fluid dynamics coupled with discrete element method (CFD-DEM) has been widely used to understand the complicated fundamentals inside gas–solid fluidized beds. To realize large-scale simulations, CFD-DEM integrated with coarse-grain model (CG CFD-DEM) provides a feasible solution, and has led to a recent upsurge of interest. However, when dealing with large-scale simulations involving irregular-shaped particles such as biomass particles featuring elongated shapes, current CG models cannot function as normal because they are all developed for spherical particles. To address this issue, a CG CFD-DEM for nonspherical particles is proposed in this study, and the morphology of particles is characterized by the super-ellipsoid model. The effectiveness and accuracy of CG CFD-DEM for nonspherical particles are comprehensively evaluated by comparing the hydrodynamic behaviors with the results predicted by traditional CFD-DEM in a gas–solid fluidized bed. It is demonstrated that the proposed model can accurately model gas–solid flow containing nonspherical particles, merely the particle dynamics are somewhat lost due to the scaleup of particle size. Finally, the calculation efficiency of CG CFD-DEM is assessed, and the results show that CG CFD-DEM can largely reduce computational costs mainly by improving the calculation efficiency of DEM. In general, the proposed CG CFD-DEM for nonspherical particles strikes a good balance between efficiency and accuracy, and has shown its prospect as a high-efficiency alternative to traditional CFD-DEM for engineering applications involving nonspherical particles.     

Methanol synthesis in a countercurrent gas—solid—solid trickle flow reactor. An experimental study

The synthesis of methanol from CO and H₂ was executed in a gas—solid—solid trickle flow reactor. The reactor consisted of three tubular reactor sections with cooling sections in between. The catalyst was Cu on alumina, the adsorbent was a silica—alumina powder and the experimental range 498–523 K, 5.0–6.3 MPa and 0.2–0.33 molar fraction of CO. Complete conversion in one pass was achieved for stoichiometric feed rates, so that the gas outlet could be closed. The experimental results are compared with the model presented in the previous paper [Westerterp, K.R. and Kuczynski, M. (1987) Chem. Engng Sci.42,]; agreement is close over the entire conversion range from 15% to 100%.     

Total and gas convective heat transfer from a vertical tube to a mixed particle gas—solid fluidized bed

Experimental data were obtained for the average gas convective and total heat transfer coefficients for a vertical tube immersed in an air-fluidized bed of narrowly as well as widely distributed particle size mixtures. The gas convective heat transfer coefficient was determined by measuring the rate of mass loss from a vertical naphthalene tube 0.0262 m in diameter and 0.1012 m in length and using a heat and mass transfer analogy. These data were obtained at a bed temperature of about 330 K and superficial velocity of 0.1 to 1.1 m/s. The total heat transfer coefficients were measured under identical conditions using an electrically heated vertical tube. The total heat transfer coefficient decreased with an increase in particle diameter from 0.237 to 1.35 mm. The addition of fines was found to increase the total heat transfer coefficient. The gas convective heat transfer coefficient increased with increase in particle size and fluidizing velocity. The dependence of the gas convective heat transfer coefficient on gas velocity was more pronounced for large particles. The addition of fines resulted in decrease in gas convective coefficient. The relative contribution of the gas convective component of heat transfer coefficient was found to increase with increase in particle diameter. Its dependency on fluidizing velocity was found to be more complex. The experimental data were compared with the existing heat transfer models and correlations.     

Horizontal gas and gas/solid jet penetration in a gas–solid fluidized bed

In this paper, the real time, dynamic phenomena of the three-dimensional horizontal gas and gas/solid mixture jetting in a 0.3 m (12 in) bubbling gas–solid fluidized bed are reported. The instantaneous properties of the shape of the jets and volumetric solids holdup are qualified and quantified using the three-dimensional electrical capacitance volume tomography (ECVT) recently developed in the authors’ group. It is found that the horizontal gas jet is almost symmetric along the horizontal axis during its penetration. As the jet width expands, the total volume of the gas jet increases. A mechanistic model is also developed to account for the experimental results obtained in this study. Comparison of jet penetration length and width between the model prediction and ECVT experiment shows that both the maximum penetration length and the maximum width of the horizontal gas jet increase with the superficial gas velocity. When the horizontal gas jet coalesces with a bubble rising from the bottom distributor, it loses its symmetric shape and can easily penetrate into the bed. For the horizontal gas/solid mixture jet penetration in the bed, the tail of the jet at the nozzle shrinks and the jet loses its jet shape immediately when the jet reaches its maximum penetration length, which are different from the characteristics exhibited by the gas jet. The solids holdup in the core region of the gas/solid mixture jet is higher than that in the gas jet. The penetration length of the horizontal gas/solid mixture jet is also larger than that of the gas jet.     

Simulation analysis of gas–solid flow characteristics and water evaporation in flue gas semi-dry desulphurization process based on CPFD method

The gas–solids flow in an industrial-scale semi-dry method desulphurization tower is simulated by the computational particle fluid dynamics (CPFD) approach. Compared with previous studies on desulphurization towers, this study focuses on analyzing particle distribution characteristics such as particle volume fraction, temperature distribution, and residence time. The simulation fully considered the particle–fluid, particle–particle, and particle–wall interactions in the desulphurization tower. Based on these considerations, the effects of flue gas inlet velocity and temperature on the gas–solid distribution characteristics of the desulphurization tower are simulated. An optimization scheme for adjusting the gas–solid flow in the desulphurization tower is proposed. The research results show that the error between the CPFD simulation data and experimental data is small and the changing trend is consistent. The particles in the bed of the desulphurization tower show a typical core–annulus flow. The distribution of gas and particles in the bed has a serious deviation with the increase of the flue gas inlet velocity and temperature. As the axial height of the desulphurization tower increases, the flue gas velocity, temperature, particle concentration, and water vapour distribution in the bed become more uniform. The relatively stable operating conditions for the gas–solid flow in the desulphurization tower is that the flue gas inlet velocity and temperature are 15 m/s and 393 K, respectively. Under these operating conditions, the pressure loss caused by the venturi accounted for 73.6% of the total pressure loss of the desulphurization tower. When the particle radius is between 0–150 μm, the particle size and the flue gas inlet velocity have the greatest influence on the particle residence time. Finally, the distribution of gas and particles before and after the adjustment of the desulphurization tower is compared, which showed that adjusting the bottom structure of the desulphurization tower could optimize the gas–solid flow.     

Simulation and experiment of gas–solid flow field in short-contact cyclone reactors

A short-contact cyclone reactor has been designed for the particular case of fluid catalytic cracking. The new type reactor mainly includes two parts: a reaction chamber and a separation chamber. So the cracking reactions and the separations between the products and catalysts could occur respectively and simultaneously. A three dimensional model was used to representing key parts of a laboratory cyclone reactor. The Eulerian–Eulerian computational fluid dynamics model with the kinetic theory of granular flow was adopted to simulate the gas–solid two-phase flow. The particle concentration distribution and pressure drop were measured by a PV-6A particles velocity measure instrument and a U-manometer, respectively. Simulated results show that in the reaction chamber solids can be transformed into a homogeneous dispersed flow, particles’ concentration becomes uniform gradually while catalysts flowing down, the concentration is a little higher near the wall because of boundary effect. After the gas–solid flowing into the separation chamber, the gas phase is separated with solids completely. The new reactor has a good contact and separation effect. Simulated results make a reasonable agreement with the experimental findings.     

Flow regime transition for cocurrent gas–liquid flow in micro-channels

Experimental investigations of gas–liquid two-phase flow regimes in micro-channels were carried out in this work. Four distinctive flow patterns of slug flow, slug-annular flow, annular flow and parallel stratified flow were captured by a digital video recording (DVR) system and the transitions among different flow regimes were studied. The effect of fluid properties and wetting properties of the channel wall on the flow regime and flow pattern transition were studied as well. Novel empirical correlations for predicting flow pattern transitions during the steady gas–liquid flow in micro-channels have been developed with the fluid properties and the wetting properties of micro-channels incorporated, whereas the traditional flow criteria were based on gas and liquid superficial velocities only. The predictions of the present empirical models agree well with the experimental data.     

Phenomenological study of saltating motion of individual particles in horizontal particle–gas systems

A comprehensive use of particle–fluid conveying systems for a wide range of industries requires a deep understanding in all interactions of the particular conveying process. One of the most common particle motions occurring in conveying systems is the saltating motion of particles. Although the literature abounds with theoretical, empirical and numerical studies that investigate the saltation phenomenon, there remain many questions and misunderstandings. Some of the recently solved issues are: which non-dimensional groups are introducing the particle saltating motion, how the saltation length might be predicted, how the pipe diameter and the coefficient of restitution influence the saltation velocity and length.The present work investigates the motion of individual saltating particles and presents a wide range of experimental measurements of the conveying length for a variety of particulate solids, sizes and shapes. The total conveying length was divided into three lengths: the first flight, the rebound and the rolling/sliding and each one of them is theoretically and empirically analyzed and compared. This phenomenological study presents the theoretical evidence to previously empirical findings. The theoretical analysis is further used to define the border conditions between various mechanisms. The results show that for coarse particles the rebound and rolling/sliding motions might be presented by a simple relationship between the Reynolds and Archimedes numbers. Additionally we find that the preferred saltation mode of fine powders depends on the conveying system length and diameter. For example for large pipe diameters and short length the first flight mode is the dominant; however, for small pipe diameters and long systems length the rebound mode is the dominant.     

Experimental research on steady-state operation characteristics of gas–solid flow in a 15.5 m dual circulating fluidized bed system

Gas-solid hydrodynamic steady-state operation is the operating basis in a chemical looping dual-reactor system. This study reported the experimental results on the steady-state operation characteristics of gas–solid flow in a 15.5 m high dual circulating fluidized bed(CFB) cold test system. The effects of superficial gas velocity, static bed material height and solid returning modes on the steady-state operation characteristics between the two CFBs were investigated. Results suggest that the solid distributions in the dual CFB test system was mainly determined by the superficial gas velocity and larger solid inventory may help to improve the solid distributions. Besides, cross-returning mode coupled with self-returning is good for steady-state running in the dual-reactor test system.     

Simulation of gas–solid two-phase flow in the annulus of drilling well

This paper investigates the hydrodynamic behavior of gas–solid two-phase flow in the annular space of an air drilling well under different arrangements by using three-dimensional approach. Two-fluid model is used to solve the governing equations in the Eulerian–Eulerian framework. Effect of eccentricity and drill pipe rotation on the pressure drop, volume fraction and velocity profile are examined. The results are compared with available data in the literature and good agreement is found. The results show that the presence of solid particles in the annulus change the air velocity profile significantly and create two off-center peaks velocity close to the walls instead of one peak velocity in the middle. Eccentricity of drill pipe makes more accumulation of the cuttings in the smaller space of the annulus. Increasing the eccentricity increases pressure drop due to impact of particles with annulus wall and also particles collision with each other. Rotation of the drill pipe shifts maximum air velocity location toward smaller space of the annulus which results more uniform cutting distributions in the annulus and improvement in their transportations. Pressure drop in the annulus increases as eccentricity and rotation of drill pipe increase.     

Effect of the injection-nozzle geometry on the interaction between a gas–liquid jet and a gas–solid fluidized bed

In industrial fluid cokers, bitumen is first mixed with steam in a premixer, and then fed to the atomization nozzle. The objective of this work was to evaluate the impact of both the premixer and the nozzle geometrical configuration on the quality of the liquid–solid contact resulting from injections of liquid into a gas–solid fluidized bed. To assess the quality of the liquid–solid contact a method based on electric conductance measurements of the bed material previously developed by the authors [9] was used. Liquid atomization efficiency in open air, spray geometry, and spray stability were also characterized to evaluate their effects on the nozzle spraying performance within the fluidized bed. This study indicated that spray stability is highly beneficial to the liquid–solid contact efficiency. In particular, fluid constrictions such as the series of converging and diverging sections within the nozzle have a stabilizing effect on the spray. Future optimization of the existing liquid-injection systems should consider alternative gas–liquid premixers and nozzle geometries to enhance the jet stability.     

Investigation on gas–solid flow regimes in a novel multistage fluidized bed

Gas–solid flow regime in a novel multistage circulating fluidized bed is investigated in this study. Pressure fluctuations are first sampled from gas–solid flow systems and then are analyzed through frequency and time–frequency domain methods including power spectrum and Hilbert–Huang transform. According to the flow characteristics obtained from pressure fluctuations, it is found that the gas–solid motions in the multistage circulating fluidized bed exhibit two dominant motion peaks in low and high frequencies. Moreover, gas-cluster motions become intensive for the multistage circulating fluidized bed in comparison with the fast bed. Unlike the traditional methods, the fuzzy C-means clustering method is introduced to objectively identify flow regime in the multistage circulating fluidized bed on the basis of the flow characteristics extracted from bubbling, turbulent, fast, and multistage fluidized beds. The identification accuracy of fuzzy C-means clustering method is first verified. The identification results show that the flow regime in the multistage circulating fluidized bed is in the scope of fast flow regime under examined conditions. Moreover, the results indicate that the consistency of flow regime between two enlarged sections exists. In addition, the transition onset of fast flow regime in the multistage circulating fluidized bed is higher than that in the fast bed.     

Comparative study of gas–oil and gas–water two-phase flow in a vertical pipe

A wire-mesh sensor has been employed to study air/water and air/silicone oil two-phase flow in a vertical pipe of 67 mm diameter and 6 m length. The sensor was operated with a conductivity-measuring electronics for air/water flow and a permittivity-measuring one for air/silicone oil flow. The experimental setup enabled a direct comparison of both two-phase flow types for the given pipe geometry and volumetric flow rates of the flow constituents. The data have been interrogated at a number of levels. The time series of cross-sectionally averaged void fraction was used to determine characteristics in amplitude and frequency space. In a more three-dimensional examination, radial gas volume fraction profiles and bubble size distributions were processed from the wire-mesh sensor data and compared for both flow types. Information from time series and bubble size distribution data was used to identify flow patterns for each of the flow rates studied.     

Profiles of solid fraction and heterogeneous phase structure in a gas–solid airlift loop reactor

The time-averaged and transient local solid fractions in a gas–solid airlift loop reactor (ALR) were investigated systematically by experiments and CFD simulations. To demonstrate the macro-flow pattern, the time-averaged local solid fractions in four regions of the ALR were measured by optical fiber probe under the conditions of different superficial gas velocities and particle circulation fluxes. The experimental results show that the lateral distribution of time-averaged local solid fraction is a core-annulus or heterogeneous structure in the three regions (draft tube, bottom region, particle diffluence region), but a uniform lateral distribution in the annulus. The operating conditions have different effects on the lateral distribution of time-averaged local solid fraction in each region. In the CFD simulation, a modified Gidaspow drag model considering the formation of particle clusters was incorporated into the Eulerian–Eulerian CFD model with particulate phase kinetic theory to simulate and analyze the transient local solid fraction and the two-phase micro-structures in the gas–solid ALR. The predicted values of solid fraction were compared with the experimental results, validating the drag model. The contours of transient flow field indicate that the flow field of the ALR should be divided into five flow regions, i.e., draft tube, annulus, bottom region, particle diffluence region and constrained back-mixing region, which further improves the understanding of the airlift reactor where only four divisions were determined from the experiments. The transient local solid fraction and its probability density function profoundly reveal the two-phase micro-structures (dilute phase and emulsion phase or cluster phase in the constrained back-mixing region) and explain the heterogeneous phenomenon of solid fraction in the ALR. The dilute phase tends to exist in the center of bed, while the emulsion phase mainly appears in the wall region. The results also indicate that the gas–solid ALR has the common characteristic of aggregative fluidization similar to that in normal fluidized beds. The simulated two-phase transient micro-structures provide the appropriate explanations for the experimental core-annulus macro-structures of time-averaged local solid fraction.     

Measurement of the impurity gas pressure during combustion of a “gasless” system in a long cylindrical shell

A method for estimating the maximum impurity gas pressure averaged over the specimen cross section in the combustion wave for different cross sections of a long specimen is proposed. The impurity gas pressure in a given specimen cross section is understood as the gas pressure that induces cross-sectional fracture of the specimen enclosed into a shell and certain displacement along the shell of the pressed matter portion located ahead of the combustion wave front at the instant of fracture initiation. Special composite specimens are used to ensure constant conditions of impurity gas removal from the volume and to model the SHS process and associated impurity degassing in the case of long (L/d γ 1) cylindrical specimens under the condition L = const. The model mixtures are chosen to be typical compositions for SHS: Ti+C and Ti + 2B.