• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Process Intensification in Pneumatically Agitated Slurry Reactors

    2021-07-08 03:51:46ShujunGengZaiShaMaoQingshanHuangChaoYang
    Engineering 2021年3期

    Shujun Geng,Zai-Sha Mao,Qingshan Huang,*,Chao Yang,,*

    a Key Laboratory of Biofuels,Qingdao Institute of Bioenergy and Bioprocess Technology,Chinese Academy of Sciences,Qingdao 266101,China

    b Dalian National Laboratory for Clean Energy,Dalian Institute of Chemical Physics,Chinese Academy of Sciences,Dalian 116023,China

    c Key Laboratory of Green Process and Engineering,Institute of Process Engineering,Chinese Academy of Sciences,Beijing 100190,China

    Keywords:

    ABSTRACT Pneumatically agitated slurry reactors,including bubble column reactors and airlift loop reactors(ALRs),are important gas–liquid–solid multiphase reactors.These reactors have been widely applied in many processes,especially in the biological fermentation and energy chemical industry,due to their low shear stress,good mixing,perfect mass-/heat-transfer properties,and relatively low costs.To further improve the performance of slurry reactors(i.e.,mixing and mass/heat transfer)and to satisfy industrial requirements(e.g.,temperature control,reduction of back-mixing,and product separation),the process intensification of slurry reactors is essential.This article starts by reviewing the latest advancements in the intensification of mixing and mass/heat transfer in these two types of reactors.It then summarizes process-intensification methods for mixing and separation that allow continuous production in these slurry reactors.Process-intensification technology that integrates directional flow in an ALR with simple solid–liquid separation in a hydrocyclone is recommended for its high efficiency and low costs.This article also systematically addresses vital considerations and challenges,including flow regime discrimination,gas spargers,solid particle effects,and other concerns in slurry reactors.It introduces the progress of numerical simulation using computational fluid dynamics(CFD)for the rational design of slurry reactors and discusses difficulties in modeling.Finally,it presents conclusions and perspectives on the design of industrial slurry reactors.

    1.Introduction

    Pneumatically agitated slurry reactors are characterized by buoyancy-driven flow with an acceptable energy input for mixing but without moving parts.The most common examples of this type of reactors are the bubble column and the airlift loop reactor(ALR).These reactors have been widely used in processes for gas–liquid and gas–liquid–solid multiphase reactions,such as biological fermentation[1],wastewater treatment[2],the mass cultivation of photosynthetic organisms[3,4],and the energy chemical industry[5]due to their simple construction,excellent mass-and heattransfer properties,and relatively low energy consumption[6].In addition to these advantages,ALRs possess low and uniform shear stress[7,8],which is necessary for biological reactions[7,9].There are two types of ALR:the internal airlift loop reactor(IALR)and the external airlift loop reactor(EALR).The hydrodynamics in ALRs(e.g.,gas holdup,solid holdup,liquid velocity,bubble-size distribution,flow regime,mixing time,residence time distribution),masstransfer properties,heat-transfer coefficient,and reaction rate,which represent the operating performance,have been thoroughly investigated with varying structural parameters(e.g.,reactor height[10],downcomer-to-riser cross-sectional area ratio[11],and gas distributor[12])and operating conditions(e.g.,temperature,pressure[5],superficial gas velocity[13],and solid concentration[14]).Based on the experimental data and energy-balance equations,various theoretical and empirical correlations have been deduced to predict the hydrodynamics and transport properties in these slurry reactors[15,16].It should be noted that the proposed relations have some restrictions on operating conditions and structural parameters due to limited mathematical statistics or the simplification of balance equations.

    Along with numerous investigations in traditional pneumatically agitated reactors,various studies have focused on the design and optimization of novel slurry reactors for process intensification for different industrial purposes.Specifically,internal tubes can be installed for heat exchange in order to maintain reactors at the desired temperature;baffles,static mixers,and mechanical internals can be equipped to enhance the overall mixing[17,18];perforated plates can be used to reduce liquid back-mixing[19,20].With the application of coated structured catalysts,no extra separation of liquid and solid catalysts is necessary[21].Hydrodynamics,mass-transport properties,and mixing characteristics have been investigated in these novel pneumatically agitated reactors,and various empirical correlations have been proposed to predict and optimize their structures.

    Great endeavors have been undertaken to improve the performance of slurry reactors and to meet the desired requirements for different industrial processes.For the Fischer–Tropsch synthesis,fixed-bed reactors,multitubular reactors,slurry reactors,and circulated fluidized-bed reactors have been commercialized[22].Industrial gas–liquid–solid slurry reactors with internal cooling tubes,which are suitable for producing diesel and paraffin wax[23],are prospective slurry reactors due to their excellent mixing and mass-/heat-transfer properties.Suspended catalytic particle separation is one of the biggest challenges in slurry reactors,igniting the development of integral catalyst design(i.e.,honeycombs,structured packings,or foams)and the integration of mixing and separation[24,25].The slurry reactor design and concerns are similar to those of other gas-to-liquid processes,such as methanol synthesis and dimethyl ether synthesis[26].Compared with conventional stirred tanks and bubble columns,the ALR has been verified to be an effective device for wastewater treatment due to its high adsorption efficiency and low shear stress to cells[27].Furthermore,various techniques,such as biofilm,membrane,electrocoagulation,ultrasound,and photochemical methods,can be conveniently combined with ALRs,providing numerous alternatives for wastewater treatment and biological fermentation[28].

    This article reviews intensification concepts and techniques for both gas–liquid and gas–liquid–solid flows in aerated reactors and presents related industrial processes.It also summarizes the relevant empirical and theoretical relationships of the hydrodynamics and transport properties in intensified slurry reactors.The article is organized as follows:Section 2 provides an overview of the intensification methods of mixing and mass/heat transfer in bubble columns and ALRs,including internals(i.e.,internal tubes,perforated plates,baffles,static mixers,structuring or packed beds,and mechanical internals),vibrating excitement,and combination methods.It also describes and compares intensified methods and their influence on hydrodynamics and transport properties.In light of the challenges of separating solid catalyst particles from continuous liquid products,Section 3 outlines separation methods applied in slurry reactors,with a focus on settling(both gravitational and centrifugal),filtration,and composite methods.Section 4 then elaborates important design considerations and challenges for the industrial design of pneumatically agitated slurry reactors.With the rapid development of computational fluid dynamics(CFD),mathematical modeling has become an efficient method for the design,optimization,and scale-up of slurry reactors.It can also be considered as an intensification method for substantial improvement of the performance in pneumatically agitated slurry reactors.Thus,Section 5 focuses on CFD modeling.Finally,Section 6 provides conclusions and perspectives on pneumatically agitated slurry reactors.

    2.Process intensification of mixing and mass/heat transfer

    It is widely accepted that the vital feature in a bubble column reactor is the non-ideal flow pattern for each phase,which can significantly influence reactant conversion and selectivity[29].However,directional flow for each phase is an outstanding advantage in the ALR[7].Therefore,the technologies of process intensification should be diversified due to the nature of their flow.

    2.1.Bubble columns

    The bubble column has been widely applied in chemical,biochemical,and petrochemical processes,such as hydrogenation,fermentation,and Fischer–Tropsch synthesis,due to its simple construction,low energy cost,and good mass-and heat-transfer properties[19,30,31].It has been reported that the bubble column reactor is now competitive with the traditional fixed-bed reactor for converting syngas into liquid fuels[32].Most reactions conducted in a bubble column must be maintained within a suitable temperature range by the addition of internal heat-exchanging tubes.Moreover,the flow structure should be reasonably controlled to enlarge the gas–liquid interfacial area and reduce backmixing in order to improve the gas–liquid mass transfer and reduce byproducts.Installing different types of internals and enforcing constant-frequency flow pulsation are two main approaches to intensify hydrodynamic interactions,mass-/heattransfer properties,and reactant conversion and selectivity in bubble columns.

    2.1.1.Intensification with internals

    2.1.1.1.Internal tubes.Many of the chemical reactions occurring in a bubble column are highly exothermic,such as Fischer–Tropsch synthesis,methanol synthesis,and acetic acid production;thus,a bundle of heat-exchange tubes is necessary to maintain the desired reaction temperature and avoid local overheating[22,33].In the 1990s,Sasol developed a slurry reactor with internally cooled tubes for Fischer–Tropsch synthesis,with a capacity of 2500 bbl?d-1(1 bbl=158.9873 L)[23].For the methanol synthesis process,a bubble column with a diameter of 457.2 mm equipped with internal tubes was developed by Alternate Fuels Development Unit in La Porte,Texas,USA[19].Fig.1[34]provides a schematic diagram of a bubble column equipped with internal tubes.It has been shown that the gas holdup and bubble passage frequency increase with the installation of internal tubes inside the column[34].At the same time,the bubble chord length and bubble-rise velocity decrease due to the enhancement of bubble breakup,especially for dense internals with high occupation of the cross-sectional area[35,36].The bundle of internal tubes is commonly arranged in a hexagonal or circular pattern(Fig.1).Different kinds and sizes of internals can influence the hydrodynamics and transfer properties in different ways and subsequently affect the reactant conversion and selectivity in bubble column reactors.Investigations on bubble columns equipped with different arrangements of internal tubes are summarized in Table 1[20,34,35,37–44].

    Fig.1.Schematic diagram of a bubble column equipped with inner tubes.Reproduced from Ref.[34]with permission of American Chemical Society,?2009.

    Kagumba and Al-Dahhan[35]investigated the hydrodynamics through four-point fiber optical probes in a bubble column equipped with a bundle of two different-diameter vertical tubes(i.e.,12.7 and 25.4 mm,respectively),with the same occupation of the cross-sectional area(25%).It was found that the local gas holdup at the column center was enhanced by up to 40% in the bubbly flow regime when the superficial gas velocity was 0.03 m?s-1,and that the promotion of gas holdup was greater with 25.4 mm internals than with 12.7 mm internals in the bubble column.However,the gas holdup only attained an average increase of 2%for both the 25.4 and 12.7 mm tubes in the churn turbulent flow regime at a superficial gas velocity of 0.45 m?s-1.Regarding the specific interfacial area and bubble passage frequency,it was observed that the values obtained when using 12.7 mm internals attained a higher increase than the values obtained when using 25.4 mm internals with a superficial gas velocity ranging from 0.03 to 0.45 m?s-1.However,Sultan et al.[45]obtained different conclusions indicating that two different-sized internals had negligible influence on the local gas holdup in the bubbly flow regime,while a higher gas holdup was obtained with 25.4 mm internals in the region close to the wall in the turbulent regime when the superficial gas velocity was 0.45 m?s-1.It was also demonstrated that the gas holdup showed a wavy radial distribution in a bubble column with internals,compared with a parabolic distribution in acolumn without internals.In addition,the shape and steepness of the radial profile of the gas holdup had a relationship with the number of internal bundles and sizes,which was also reported by Al Mesfer et al.[46].

    Table 1 Summary of investigations on bubble columns equipped with internals.

    Pradhan et al.[37]investigated the influences of the volume fraction occupied by internal tubes on the overall gas holdup by varying the superficial gas velocity.When the volume fraction of the internals was increased from 1.4% to 19.3%,the gas holdup increased monotonously,regardless of whether straight tubes or helical coils were installed.Moreover,a higher gas holdup was obtained with helical coils than with straight tubes.To be specific,the gas holdup was increased by about 19.5% when the straight tubes occupied a volume fraction of 10.8%,compared with 48.4%enhancement when the helical coils occupied a volume fraction of 3.9%,for a superficial gas velocity of 0.41 m?s-1.Youssef and Al-Dahhan[34]also reported that the overall and local gas holdup and the interfacial area increased when the occupation of the cross-sectional area covered by internals was increased,especially at the center of the bubble column.

    Saxena et al.[38]investigated the gas holdup of both air–water and air–water–glass systems in a bubble column equipped with five sparsely arranged tubes,and with seven tubes and 37 tubes with tightly hexagonal arrangements,occupying 1.9%,2.7%,and 14.3%of the column cross-sectional area,respectively.Their results demonstrated that the overall gas holdup in the bubble column equipped with 37 tubes was notably higher than that in the columns with five and seven tubes,respectively,when the superficial gas velocity was above 0.1 m?s-1.However,when the superficial gas velocity was less than 0.1 m?s-1,the differences in gas holdup between columns with differently arranged tubes were obscure,and the column equipped with seven central tubes achieved the highest gas holdup at ambient temperature.It was also observed that the centrally arranged internal tubes were apt to increase the gas holdup due to the enhancement of bubble breakup in the center region.For a three-phase system,the gas holdup decreased with increasing particle size from 125 to 212μm when the solid concentration was greater than 5%,with temperature varying from 298 to 353 K.In contrast,the variation tendency of the gas holdup versus the particle diameter was sensitive to temperature for dilute slurries(about 5 wt%).The gas holdup decreased with increasing particle diameter at ambient temperature,with no variation being observed above 323 K.

    Based on the discussion above,the enhancement of hydrodynamics by internal tubes is related to their size,occupation area of the cross-section,arrangement,and operational conditions.In general,greater promotion of gas holdup is obtained with greater occupation of the cross-sectional area,and regular centrally arranged internal tubes are more beneficial for bubble breakup.

    2.1.1.2.Perforated plates.Perforated plates are one of the most traditional internals used to intensify the hydrodynamics and mass/heat transfer in bubble columns[47].As early as 1962,Fair et al.[39]investigated gas holdup and heat transfer in pilot-scale bubble columns with an air–water system.The bubble columns were 457.2 or 1016 mm in diameter and 3048 mm in height and were equipped with 1–20 static or moving perforated plates.The pore diameter of the perforated baffles varied from 3.175 to 7.925 mm;the opening area ranged from 9.5% to 33.0% of the cross-sectional area;the numbers of equipped baffles were 5,10,15,and 20,respectively.With the assembly of static plates,the gas holdup was increased by 40%–50% compared with that in an empty bubble column and further increased by 25%–30% caused by the moving perforated plates when the superficial gas velocity was changed from 1.524 to 10.06 cm?s-1.Analogously,heat transfer was promoted by 10%–15% as a result of both stationary and moving perforated plates.The enhancements of perforated plates in bubble columns are summarized in Table 1.

    Intense liquid back-mixing is a distinct feature of the slurry bubble column reactor and may reduce the selectivity and conversion of reactant.It has been reported that the liquid axial dispersion can be regulated by sectioning a bubble column with perforated plates[48].The extent of liquid back-mixing decreases with a decrease in pore diameter and in the free area percentage of the orifice plates[40].This decrease has been ascribed to the uniform redistribution of the dispersed phase at each plate with the increase of gas flow resistance and minimization of the density gradient.Nevertheless,the new problem of choking has been found to occur much more easily with a low opening fraction and small pore diameter of sieve plates.Thus,the appropriate design of the number,pore diameter,and free area percentage of perforated plates is important to achieve the desired hydrodynamics and liquid dispersion characteristics[48].The axial dispersion model is commonly employed to depict the residence time distribution of the liquid phase in bubble columns[49,50].It was found that the axial dispersion coefficient was related to an exponential function of the superficial liquid velocity,with an exponent ranging from 0.234 to 0.85[40]when the net liquid velocity relative to the rising bubbles was less than 10%of the bubble-rising velocity[51].However,the liquid axial dispersion coefficient was independent of the superficial liquid velocity if the circulation velocity was ten times higher than the superficial liquid velocity[52].Furthermore,the liquid axial dispersion coefficient was found to be dependent on an exponential function of the superficial gas velocity,with an exponent ranging from 0.12 to 0.33[40,52].Dreher and Krishna[20]studied the liquid residence time distribution using the axial dispersion model with the liquid exchange velocity for quantitative calculation.Their results showed that the liquid exchange velocity was independent of the column diameter,which was important for the scale-up of multi-stage bubble column reactors.

    2.1.1.3.Static mixers.Static mixers have been widely applied to disperse the gas phase in liquid bulk in chemical and industrial processes due to their minimal space requirement and low equipment cost,and the fact that they do not require extra power supply[53].The commercial Inliner mixer,produced by Lightnin Inc.,has been applied in petrochemical industrial processes such as hydrocarbon refining[53].With a static mixer inserted,the gas–liquid interfacial area in the column is increased,and the mass/heat transfer and chemical reactions are intensified.Fig.2(a)[41]shows a schematic diagram of a bubble column equipped with static mixers,and Fig.2(b)[53]depicts common types of static mixer.The rate of gas–liquid mass transfer was nearly double in bubble columns packed with Koch static mixers,compared with that in an unpacked bubble column,as investigated by Fan et al.[54]and Wang and Fan[55]with a liquid velocity ranging from 6.7 to 39.9 cm?s-1.These researchers also deduced a correlation of the volumetric mass-transfer coefficient(kLa),with Vl(where Vlis the liquid velocity)as an independent variable.

    Gaspillo and Goto[42]studied the influence of a static mixer in a draft tube on the minimum gas velocity to suspend the solid particles and mass-transfer coefficient in both gas–liquid and gas–liquid–solid slurry bubble columns.It was demonstrated that the minimum gas velocity for slurry suspension decreased with the presence of a draft tube but slightly increased with the addition of a static mixer,which might be caused by the turbulent motion of solid particles in the static mixer.When a single-nozzle gas distributor was used,the gas–liquid mass-transfer coefficient increased by about 34% with the static mixer installed.Nevertheless,when an orifice ball distributor was applied to produce extremely fine bubbles,the gas–liquid mass-transfer coefficient decreased due to bubble coalescence caused by the static mixer.The rate of gas–liquid mass transfer increased by about 20% with the static mixer equipped in the draft tube,regardless of which kind of gas distributor was adopted.Table 1 depicts the improvement in hydrodynamics and mass transfer.

    Fig.2.(a)Schematic diagram of a bubble column equipped with static mixers;(b)common types of static mixer.LPD:low pressure drop;HEV:high efficiency vortex.(a)Reproduced from Ref.[41]with permission of Elsevier B.V.,?2014;(b)reproduced from Ref.[53]with permission of Institution of Chemical Engineers,?2003.

    Recently,Rabha et al.[41]evaluated the length of the helical mixer on the gas holdup,gas-phase flow structure,interfacial area,and bubble-size distribution using ultrafast electron beam X-ray tomography.Three,six,and nine helical static mixers that were 80 mm in diameter,80 mm in length,and had a 180-degree rotation of the blade,were inserted in the column to investigate their impact.The liquid velocity and helical structure induced extra shear,forcing bubble breakup.At a low flow ratio(Ug/Ul),the bubble size distribution was notably narrowed,and large bubbles suddenly decreased in size after flowing through these static mixers.Compared with the bubble diameter at the same axial position before the mixers,the bubble Sauter diameter decreased by about 42% regardless of the number of static mixers,as displayed in Fig.3(a).With an increase of superficial gas velocity,the quantitative reduction of the bubble Sauter diameter decreased and the number of large bubbles increased due to bubble coalescence during turbulent swirling motions,as shown in Fig.3(b)[41].Nonetheless,when the flow ratio was increased with a low liquid velocity,the large bubbles could not be clearly dispersed by the static mixers,resulting in a slight increase in the bubble Sauter diameter.Similarly,it was observed that the specific gas–liquid interfacial area varied with the mixer elements because the interfacial area was directly related to the bubble diameter and liquid circulation.

    Thus,to realize process intensification in a slurry bubble column,the structure and scale of the static mixer elements should be appropriately designed while considering the gas distributor,other internals,and the operating conditions.

    2.1.1.4.Structuring.Liquid–solid separation and reducing backmixing are typical challenges in a pneumatically agitated slurry reactor because the separation process can increase the operating and capital costs,while back-mixing decreases the selectivity and conversion of reactant.Structuring the reactor by equipping regular internals has been identified as an efficient method to intensify the slurry bubble column reactor to avoid solid–liquid separation and to reduce back-mixing.On the one hand,catalysts are structured into monoliths or coated on regular internals in the reactor so that it is not necessary to separate solid catalyst particles from the liquid products.On the other hand,the exact shape and size of column internals are designed to guide the flow to reduce liquid back-mixing.Furthermore,the internal structure can be designed precisely according to catalytic,reaction kinetics,and heattransfer characteristics,making scale-up much easier[21,56].Structured packing has been successfully applied in industrial absorption and distillation processes and then applied in the chemical and petrochemical industries[57].Meanwhile,a packed bed is suitable for acetic acid production due to its high surfacespecific area[58].

    Structured packing can be classified into four kinds:monoliths with parallel straight channels,corrugated sheets with open or closed flow structure,knitted wire packing,and open-celled foams[59].High porosity,high surface area,and a low pressure drop are common and essential features of structured packing,in order to provide a sufficient specific surface area for heat and mass transfer and to restrict energy consumption.Structured packing has been proven to promote radial heat[21]and mass-transfer rates[60],thus increasing the reactant conversion and productivity[56].Although the structured bubble column has been widely investigated,most relevant research has focused on the co-current tricked bed[59,61,62]rather than on the countercurrent bubble column.

    Using open cross-flow structure(KATAPAK-S elements,supplied by Sulzer Chemtech)packing in two bubble column reactors with diameters of 0.1 and 0.24 m,respectively,Urseanu et al.[43]investigated the hydrodynamics and mixing characteristics in comparison with those in an empty bubble column.Fig.4[43,59]shows a schematic diagram of the experimental setup.The experimental results showed that the total gas holdup in the structured bubble column was almost the same as that in the empty bubble column at the same superficial gas velocity based on the free cross-sectional area.It could then be inferred that the existence of structured packing did not promote additional bubble coalescence and breakup.In regard to axial mixing characteristics,the axial dispersion factor—which was developed by Krishna et al.[49]—decreased by about one order of magnitude at a superficial gas velocity of up to 1.25 m?s-1due to inhibition of the overall axial recirculation of the liquid phase.Moreover,axial dispersion could be fitted into a linear correlation with the liquid velocity for both the empty and structured bubble columns.

    The gas holdup and liquid-to-wall mass transfer in structured and empty bubble columns were measured in a square sectional reactor,as described in Ref.[44].The results showed that the gas holdup increased with an increase of the superficial gas velocity and liquid viscosity in the packed bubble column.Under the same operating conditions,the gas holdup in the structured column was higher than that in the empty column.Similar correlations of mass transfer were derived for the structured packing and empty column,regardless of whether Newtonian or non-Newtonian fluids were used.This result indicated that structured packing had a negligible influence on the liquid–wall mass transfer.

    2.1.2.Intensification with vibrating excitement

    Vibration exciters,such as pulsators,bellows,pistons,electric cams,and electromagnetic vibrators,can be installed at the bottom of the column to generate the necessary frequencies of vibration directly in the slurry within the reactor.Fig.5 provides a schematic diagram of a common pulsing bubble column.The introduction of overall pulsations to the liquid in a bubble column has been substantially confirmed as an effective method for improving the gas–liquid mass-and heat-transfer properties[63,64].Numerous studies have been conducted to explore the influence of oscillations on the hydrodynamics and mass-transfer properties,with the aim of revealing the relationships between the performance and operating conditions in a bubble column.Table 2[64–73]summarizes relevant investigations and critical findings on intensification by vibration.

    Fig.3.Influence of static mixer length on bubble diameter(a)quantitatively and(b)qualitatively under different liquid and gas superficial velocities.Reproduced from Ref.[41]with permission of Elsevier B.V.,?2014.

    Fig.4.(a)Schematic diagram of a structured bubble column including a packed element;(b)typical examples of cordierite monoliths.(a)Reproduced from Ref.[43]with permission of Elsevier Science B.V.,?2001;(b)reproduced from Ref.[59]with permission of American Chemical Society,?2008.

    Fig.5.Schematic diagram of a bubble column intensified by vibration.

    Ellenberger et al.[64–67]systematically investigated the bubble diameter,gas holdup,and mass-transfer characteristics in air–water and air–water–silica bubble columns subjected to lowfrequency vibration(40–200 Hz).For a single-orifice gas distributor under a constant superficial gas velocity,the bubble diameter first decreased sharply with an increase in vibration frequency and amplitude,and then ultimately reached a relative plateau at a vibration frequency of around 100 Hz.Similar trends of bubble diameter versus vibration amplitude were observed;however,the inflection point varied with the superficial gas velocity.The bubble diameter was reduced by 40%–50% when the orifice gas velocity ranged from 0.2 to 3 m?s-1with a pulsing frequency and amplitude of around 100 Hz and 1 mm,respectively.To analyze the mechanisms of bubble breakup in a pulsing bubble column,Knopf et al.[68,69]applied a high-speed camera to observe the bubble-formation process in a glass gas T-injector.The liquid phase was forced by a flexible piston at low frequencies(0–30 Hz)with a low amplitude.When the gas flow rate was lower than 5 m?s-1(the corresponding superficial gas velocity was lower than 0.0064 m?s-1),mixing of the gas and liquid phases in both the ingesting and expulsion parts caused the gas phase to break into a wide range of bubbles.Meanwhile,external sinusoidal pulsations resulted in water suck-back into the injector at a high gas velocity,causing the gas slugs to fragment intensely into tiny bubbles.With bubble breakup and the formation of smaller bubbles,the gas holdup and gas–liquid specific interfacial surface area increased,and the mass-transfer process between phases was subsequently enhanced.

    Table 2 Summary of investigations on a bubble column with vibrating excitement.

    Several researchers[74,75]have focused on seeking a general theory to explain the experimental observations and provide guidance for the design,optimization,and scale-up of the pulsing bubble column.Waghmare et al.[70]developed a one-dimensional theory of the mass-transfer coefficient by taking the superficial gas velocity,liquid viscosity,and vibration frequency and amplitude as independent variables.The balance of the buoyancy and drag force exerted on gas bubbles in liquid pulsations could be achieved by considering three additional effects.The first effect is the Bjerknes force[70],which is caused by interactions between the forcing liquid and the pulsing bubbles.This is known as a kinetic buoyancy force acting in the opposite direction to the upward movements of the bubbles,resulting in a decrease of bubble rise velocity.The second effect is bubble breakage enhanced by pulsation based on the Hinze theory[74],which considers both the gas injection and the oscillation power.The third effect is the decrease in mass diffusivity resulting from the increase in liquid viscosity,which can be obtained according to the Stokes–Einstein equation.When these three effects were considered,the final prediction model of the mass-transfer coefficient provided a fairly reasonable alignment with the experimental data for the pulsation bubble column.Table 2 displays the correlation for predicting the volumetric mass-transfer coefficient.A defect of this model was its assumption that the amplitude of the liquid oscillation was equal to that of the exciter vibrations.In fact,it has been shown that the liquid oscillation amplitude is related to the types of vibration exciter[69].With a correction based on the ratio of the gas injector power to the total power(including the injector and pulsation power),a prediction model of bubble size was deduced based on Davidson and Schuler’s relation[71].This correlation was successfully validated by experimental data with a deviation of 10%,as shown in Table 2.

    According to the analysis mentioned above,the intensification of mass transfer between phases has two reasons.First,the presence of oscillations provides the bubbles with an additional downward Bjerknes force,leading to a decrease in bubble-rise velocity and an increase in gas–liquid contact time[70,71].Second,vibrations introduced to the continuous phase enhance the breakage rates of bubbles in the bubble column.This promotes the formation of many smaller bubbles and magnifies the specific interfacial area for gas–liquid interactions.Finally,the mass transfer can be significantly improved by pulsations.This enhancement of mass transfer is related to the superficial gas velocity and the amplitude and frequency of the pulsations.There are usually optimum operating conditions under which great enhancement of the gas holdup and mass transfer can be achieved while requiring a low additional pulsation power.It should be noted that the type of pulsing piston can influence the amplitude of the liquid oscillation,and subsequently influence the bubble breakage dynamics.

    2.1.3.Combinational process intensification

    2.1.3.1.Multiple types of internals.A combination of perforated plates and straight tubes can reduce the overall back-mixing in a bubble column and provide relatively isothermal conditions,which is suitable for strongly exothermic reactions such as Fischer–Tropsch synthesis and methanol synthesis[76].Maretto and Krishna[77]simulated the syngas conversion and productivity of the Fischer–Tropsch synthesis process in a multistage bubble column with a diameter of 7 m and a height of 30 m,and heat-exchange tubes.Fig.6[77]shows the schematic diagram.The flow regime of the gas phase was assumed as the plug flow in each stage,with the dense phase(i.e.,slurry phase)remaining completely mixed under the isothermal condition.The modeling results demonstrated that the selectivity and productivity of the Fischer–Tropsch reaction notably increased with an increase in the number of perforated plates from 1 to 4,and the flow regime changed to plug flow when the number of sieves was greater than 4.Moreover,additional heat-exchange tubes needed to be installed to ensure the desired isothermal conditions in each compartment of the column.

    2.1.3.2.Combination of different intensification methods.Installing orifice baffles and producing oscillations at a constant frequency are conventional methods for intensifying the mixing in bubble columns,as stated in Sections 2.1.1.2 and 2.1.2.A combination of these two methods has also been widely applied to intensify the hydrodynamics and mass transfer in column reactors.Table 3[78–85]summarizes the research in this area.

    Fig.6.Schematic diagram of a bubble column equipped with vertical tubes and sieve plates.Reproduced from Ref.[77]with permission of Elsevier Science B.V.,?2001.

    Lucas et al.[78]applied an oscillatory baffled column reactor to the ozonation of wastewater containing p-hydroxybenzoic acid(p-HBA).In comparison with the results from a traditional bubble column,a 20% increase in the p-HBA degradation rate and a 4.5-to 5.0-fold increase in the p-HBA mineralization rate per mole of ozone were obtained.The intensification of the ozonation process was attributed to the enhanced gas holdup and mass transfer due to the presence of orifice baffles and pulsations.Dissolution of carbon dioxide in the water was also enhanced by the use of a multi-orifice oscillatory baffled column[79].It was shown that the baffle design parameters—including the orifice diameter,orifice number,and open area—had a significant effect on the bubble-size distribution and mass-transfer rate.Table 3[78–85]illustrates the relevant investigations and the main conclusions or relationships that were obtained.

    Five different types of oscillatory bubble columns,including a helical baffle reactor,smooth periodic constriction reactor,single-orifice baffled reactor,multi-orifice baffled reactor,and oscillatory bubble column reactor without baffles,were designed and investigated by Ahmed et al.[80].The kLa in each oscillatory column reactor was diligently measured and compared for different oscillatory conditions and superficial gas velocities.It was demonstrated that a significant enhancement of the masstransfer rate was obtained with the addition of orifice baffles and periodic constriction.Among these distinctive designs,the highest kLa was achieved in the oscillatory multi-orifice baffled reactor and was up to seven times higher than that in a reactor without baffles and oscillation under the same operating conditions.

    Ni and Gao[81]first investigated the scale-up parameters in two different oscillatory baffled columns and obtained a correlation of kLa that adopted the power density as an independent variable,as shown in Table 3.It was found that there was an optimal baffle space where the vortices could be delivered without dispersion and suppression between the baffles.Their research indicated that the best baffle space was 1.8 times the tube diameter.Similarly,the optimal gas rate was obtained when the interfacial area reached a maximum without the occurrence of gas channeling.With an oscillatory frequency ranging from 2 to 6 Hz and an amplitude varying from 4 to 12 mm,the kLa monotonously increased with the pulsing frequency and amplitude.Oliveira et al.[82–84]investigated a pulsed baffle column with single-orifice baffles using a high-speed camera and inferred empirical relations between the gas holdup and bubble size,with power density as avariable.Furthermore,the correlation of kLa was substantially improved,as shown in Table 3.Smith and Mackley[85]investigated the axial dispersion and liquid back-mixing in three oscillatory column reactors with different scales(i.e.,24,54,and 150 mm)and deduced a correlation in which the axial dispersion was independent of the tube diameter under the conditions of a constant net flow Reynolds number(Re),oscillatory Reynolds number(Reo),and Strouhal number(Str).The minimum axial dispersion occurred in the column when the value of Reoe-0.4Strwas equal to 1.37-fold Re0.8.

    Table 3 Summary of investigations on oscillatory baffled bubble columns.

    2.1.4.General remarks

    Numerous methods have been elaborated thus far,including equipping different internals,vibrating excitement,and combinational techniques to intensify bubble column performance.All of these intensification methods can be applied to promote gas holdup and mixing properties.Mass-transfer processes are intensified by these methods due to the generation of small bubbles,increase of the interfacial area,and generation of vortices.

    For internal tubes,the enhancement of gas holdup increases with an increase in the occupation of the cross-sectional area.However,in regard to the influence of tube size on hydrodynamics,controversy still exists for different flow regimes.In a highly exothermic industrial process,internal tubes with 22%–25% occupation of the cross-sectional area are essential for heat exchange[35].With the application of perforated plates,liquid backmixing can be significantly reduced,and the liquid superficial velocity decreases with the decrease in open area.Static mixers installed in the reactor are applied to disperse large bubbles,resulting in an increase in the gas holdup and volumetric masstransfer rate.Structured packing with a specialized design can also be used as static mixers and can simultaneously reduce undesirable back-mixing.Moreover,catalytic particles can be bound or loaded on the packing,and then equipped as a monolith in the reactor,eliminating the problems of liquid–solid separation.It has been estimated that structured packing could be applied in 25% of catalytic processes worldwide[86].With the presence of oscillations at a constant frequency,the gas–liquid mass transfer is enhanced by the increase of the gas–liquid interfacial area resulting from the pulsation,and the contact time between gas and liquid can be extended due to the additional downward Bjerknes force.With the use of this vibrating excitement method,no additional intense liquid back-mixing or high shear—which should be avoided in some biotechnological processes—is induced in the reactor.

    In comparison with a single intensification technique,a combination of different methods holds more potential for efficient process intensification from multiple aspects.It should be noted that no matter which method is adopted,the optimal process intensification can only be achieved using a suitable design and appropriate matching with the operating conditions in the bubble column.

    2.2.Airlift loop reactor

    Based on bubble column technology,the ALR was developed by installing concentric tubes,splitting a vessel into two vertical halves,or connecting two separated columns at the top and bottom.Compared with a bubble column,an ALR possesses better mixing ability,better energy-saving,and lower shear stress[87].Thanks to these dominant advantages,the ALR has been widely employed in wastewater treatment,microbial fermentation,and Fischer–Tropsch synthesis[88].To further improve ALR performance,numerous methods and techniques for process intensification have been continuously developed.The gas holdup,bubblesize distribution,mass transfer,mixing time,and extent of backmixing are the principal parameters that determine the hydrodynamics,mass/heat transfer,and finally reactant conversion and selectivity in the reactor.Thus,the process intensification method focuses on these main parameters.Similar to the bubble column,internals insertion is one of the most common methods for process intensification in an ALR.The following sections concisely elaborate on the methods used for process intensification in the ALR.

    2.2.1.Baffles

    Baffles can be applied in ALRs to improve the mixing properties,enhance the mass-transfer rate,and transform the flow direction.In order to enhance the solids mixing,especially at the bottom of the column,and improve the heat exchange at the external wall in a photosynthetic process,helical flow promoters with an angle of 35°with respect to the column axis were installed in the downcomer of an ALR[89,90].With these inserted helical baffles,radial mixing was enhanced by the helical movement of the gas–liquid–solid slurry,resulting in a more homogenous distribution of energy and suspended solid particles.Due to the resistance effects of the baffles,the liquid circulation time increased by less than 15% for both the water and dilute carboxy-methyl cellulose(CMC)solutions.Compared with a bubble column without internals,half of the gas flow rate was sufficient to arrive at a complete fluidization state with a helical flow promoter.The effects of solid concentration on the gas holdup,liquid velocity,and mass-transfer rate were thoroughly investigated.It was observed that the gas holdup and liquid velocity decreased slightly with an increase in the solid concentration.However,the gas–liquid volumetric mass-transfer rate decreased significantly.Quantitatively,with the addition of 206.5 g?L-1of solid particles,the volumetric mass-transfer rates in water and in 0.2 wt% CMC solutions decreased by 37% and 47%,respectively.Based on substantial experimental data under different operating conditions,the empirical relations of the gas holdup,liquid velocity,and kLa with various influencing factors were summarized,as shown in Table 4[17,90–97].

    To observe the detailed flow structure and identify the flow trajectories in an IALR,Wu and Merchuk[98]applied an optical trajectory tracking system.Fluorescent particles were employed as tracers,and successive pictures were taken by two cameras connected to computers.This allowed the three-dimensional coordinates,velocity vector,and local shear strain rates of the tracer particles to be recorded and analyzed by digital image analysis;the flow paths were then drawn.Through this optical working system,the existence of secondary flow in an IALR with a helical promoter was first verified,and the mechanism of helix baffles improving radial mixing and heat exchange with the external walls was identified.R?s?nen et al.[91]combined a helical flow promoter with a gas sparger to intensify the mixing and masstransfer properties.They attached helical tubes with 0.6 mm needle to rotate on the outer side of the downcomer and/or on the inside of the riser column of the IALR.Different helical flowpromoter/gas sparger combinations were designed to investigate their influence on the hydrodynamics and mass-transfer rate.It was found that the gas holdup and kLa could be significantly promoted without additional energy input.The value of kLa for oxygen could be increased threefold by the use of both a helical flowpromoter and gas sparger,compared with that in an IALR without internals.

    In addition to enhancing the hydrodynamics using helical internals,baffles of different configurations can be applied to guide the flow in the column.Pi et al.[99]equipped a trumpet-shaped hood above the internal tube in an ALR.They found that the liquid circulated in two-layer cyclic routes both below and above the hood,and that an extra circulating route was formed in a suitable installation position.With the trumpet-shaped riser and an imaged riser,the flow structure tended to be more uniform and the dead space in the reactor decreased.With the installation of this novel riser,All the symbols in this table are defined in the Nomenclatures list at the end of this paper.the overall gas holdup and kLa were enhanced by 16.2%and 10.2%,respectively.Li et al.[18]set an inverted bell-shaped internal between two stages in a two-stage IALR.The slurry from the riser of the bottom stage was compressed into the second stage with the constriction of the internals,resulting in greater gas holdup in the upper stage.Different bubble circulation regimes existed in the two stages of the reactor under the low superficial gas velocity;that is,a complete bubble circulation regime occurred in the first stage while a transition regime occurred in the second stage.With an increase in the superficial gas velocity and solid concentration,the bubble circulation gradually transformed into a complete circulation regime.It was reported that the axial distribution of the solid concentration tended to be homogeneous when the superficial gas velocity and solid loading were increased in the whole reactor.

    Table 4 Summary of theoretical and empirical correlations in intensified ALRs.

    2.2.2.Perforated plates

    Horizontal perforated plates with multiple orifices have been demonstrated to be an efficient method to induce bubble breakup and enhance the mixing extent in aerated reactors.Krichnavaruk and Pavasant[100]investigated the influence of sieve plates with different pore diameters and pore numbers on the kLa,gas holdup,and riser liquid velocity.It was found that the value of the kLa increased to as much as double,compared with the value without perforated plates.This substantial improvement in the kLa was attributed to the increase of the gas–liquid interfacial area with frequent bubble breakup by the sieve plates.Nevertheless,the gas–liquid mass-transfer coefficient(kL)decreased when perforated plates were installed because of the decreased liquidcirculation velocity resulted from the hindrance caused by the sieve plates in the flow pathway.These two contrary effects made it difficult to create an optimal design of the sieve plates in order to achieve the highest improvement of kLa.In that study,it was observed that three perforated plates with 21 holes that were 4 mm in diameter were the best choice to increase the kLa at a superficial gas velocity ranging from 1.889 to 7.515 cm?s-1.

    Luo et al.[92]reported that a larger orifice diameter on the sieve plates,varying from 2.5 to 4.5 mm,provided greater enhancement of the kLa for the same free area ratio.This phenomenon was attributed to two aspects.First,with a smaller pore diameter,large bubbles tended to congregate to form a bubble layer below the sieve plates,which led to a decrease of the gas–liquid interfacial area.Second,the flow resistance was strengthened with the decrease in sieve diameter,resulting in a greater decrease in liquid velocity.Furthermore,correlations of the overall gas holdup and kLa related to the sieve plate structural parameters and operating conditions were derived based on the experimental data,as shown in Table 4.

    Aside from the significant intensification techniques discussed above,Vorapongsathorn et al.[101]reported that perforated plates(smaller than the riser cross-section)with eight 3 mm holes placed along the riser of the IALR slightly improved the gas holdup and kLa.This was ascribed to the fact that the orifice plate occupied only half of the riser sectional area,with a fairly small open area.When these perforated plates were equipped,the fluids flowed around the sieve plates following a path of lower resistance.Under these circumstances,the internals had a limited positive influence on the hydrodynamics,except for acting as a hindrance to the fluid.Thus,these findings indicated that a suitable perforated plate design was crucial to positively enhance the hydrodynamics in a pneumatically agitated reactor.

    Zhang et al.[102]set up a novel perforated internal structure containing baffles maintained at 45° with the vertical axis in the riser tube of an EALR.The influence of these internals on the bubble Sauter diameter and rise velocity,as well as on the local gas holdup and its radial distribution,was investigated by comparing the hydrodynamics before and after the installation of sieve plates.It was verified that the perforated plates enhanced the gas holdup,narrowed the bubble-size distribution,and reduced the radial maldistribution.Nevertheless,the influence of the internals on the hydrodynamics in that study was effective only for a limited distance of about 1.1 m above the sieve plate.

    Zheng et al.[93]assembled helical sieve plates in the riser of an IALR.With the helical plates,the gas holdup was enhanced by 38%–63% in comparison with that of a setup without plates.The researchers proposed that the enhancement was independent of the free area ratio,which ranged from 35% to 63%,and the helix angle,which ranged from 10°to 31°.Similar to the effect of planar sieve plates,the bubble diameter was decreased,and there was a more uniform size distribution due to frequent bubble breakup through the plates,which occurred every time.As a result,the gas–liquid interfacial area increased,and the volumetric mass transfer was improved.Quantitatively,the kLa increased by about 20% with the use of helical sieve plates in the IALR,in comparison with an empty IALR.Empirical correlations of the gas holdup and kLa,with the influences of the superficial gas velocity,free area ratio,and helical angle,were deduced by fitting the experimental data by the least square method.Table 4 summarizes these empirical equations and valid ranges.It was shown that the homogenous flow regime of the IALR was broadened with the use of helix sieve plates.The mixing time was extended due to the space partition and resistance to the gas–liquid flow provided by the helical plates,while the tangential mixing performance was improved.This phenomenon was also observed in the ALR assembly with a helical promoter in the downcomer[89,90].

    Perforated plates can also be applied between stages in multistage ALRs.Yu et al.[103,104]designed a perforated plate equipped not only with multiple orifices but also with three long tubes to improve the distribution of solid particles and to avoid liquid flooding in the column.It was found that the novel inter-stage internals separated the gas–liquid–solid slurry flow into two parts:The gas flowed upward through the orifices,as in an ALR with traditional sieve plates,while the liquid and solids flowed upward through the tubes.The flow regime and gas layer height below the sieve plates were significantly influenced by the superficial gas velocity and the opening area of the gas orifice for both the co-current and countercurrent flow.The gas layer height should be strictly controlled within an appropriate range when the slurry flow circulates normally within each stage.

    The effects of various configurations of perforated plates(i.e.,planar or helical)with different open areas and orifice diameters on the hydrodynamics and mass-transfer properties in ALRs have been investigated under a wide range of operational conditions.Due to contrary influences such as bubble breakup and flow resistance,a suitable design of the equipped sieve plates is a key parameter for positive process intensification.Thus,a suitable design and mounting position should be carefully considered in order to achieve the desired purpose.

    The old man was fond of all creatures, and every morning he used to open the cage door, and the sparrow flew happily about until it caught sight of a cat or a rat or some other fierce beast, when it would instantly return to the cage, knowing that there no harm could come to it

    2.2.3.Static mixers or packed beds

    As discussed in Section 2.1.1.3,static mixers have been widely employed to intensify the hydrodynamics and interphase mass transfer,due to their enhancement of the bubble breakup and flow disturbance.Likewise,static mixers can be applied in ALRs to achieve similar effects[96].Chisti et al.[94]conducted experiments in an EALR equipped with static mixers in the riser and demonstrated that the kLa could be significantly enhanced by static mixers.Furthermore,the kLa enhancement increased with an increase in fluid viscosity because larger gas bubbles coalesced more easily in viscous liquid without motionless internals.Quantitatively,the kLa was enhanced up to sixfold compared with that in the same setup without static mixers in a NaCl solution with 0.6%w/v(gram per 100 mL)CMC.Prominent enhancement of the kLa by static mixers in the riser was also observed by Goto and Gaspillo[95]in a gas–liquid–solid system.Their results showed that the minimum gas velocity required to completely suspend the solid was decreased by about 30% when static mixers were used,which is beneficial for extending the application range of an EALR.Table 4 shows the related empirical correlations based on the experimental data.

    Packed beds have been widely applied as a support for catalysts or immobilized microorganisms and enzymes in multiphase chemical reactors[58].Porous materials,such as fibrous cotton,glass,nylon,and polymer foam,have been identified as suitable for packing in aerated reactors,due to their high specific surface area,low mass-transfer resistance,and small pressure drops[105,106].Meng et al.[97]and Nikakhtari and Hill[107]thoroughly investigated the hydrodynamics and mass-transfer properties in EALRs equipped with woven nylon packing.With the placement of mesh packing internals,the bubble Sauter diameter decreased significantly due to frequent bubble breakup.It should be noted that only an appropriate arrangement of the porosity and height of the packed bed resulted in a positive enhancement of the hydrodynamic properties.When the packing porosity was 0.9,the gas holdup in the column was decreased by installing the packing internals.In that study,the maximum gas holdup was obtained with a porosity of 0.99.When the packing porosity was increased,the gas holdup and liquid velocity in the riser increased monotonously.The liquid velocity in the riser decreased with increasing packing height,while the gas holdup was nearly invariant.With a packing porosity of 0.99 and a packing height of 1.2 m,the gas holdup and kLa between oxygen and water at a superficial gas velocity of 0.005 m?s-1were improved by about 74%and 174%,respectively.For the volatile organic chemicals toluene and benzene,the overall kLa was enhanced on average by 65.1% and 33.4%,respectively[108].This novel packed bed has also been reported to be efficient for the bioremediation of phenol-polluted air steam,with all the phenol being totally removed at one-third of the height of the bubble column[109].

    Mixing properties versus superficial liquid velocity were explored using non-intrusive electrical resistance tomography in a circulated EALR with equipped two-stage packing internals and a gas redistributor[110].It was reported that the liquid velocity decreased with the installation of packing and a gas redistributor,due to hindrance from the internals[111].Correspondingly,the mixing time increased by up to 44% with one bed of packing and further increased by 270% with the combination of two beds of packing and a gas redistributor,because the resistance of the liquid flow was significantly increased.Moreover,the overall gas holdup increased by about 36%due to the bubble breakup and fluid resistance[112].Ultimately,the kLa was clearly enhanced because of the improved gas dispersion.

    2.2.4.Mechanical internals

    Mechanical internals,such as blades or impellers,have been extensively applied in stirred tanks and bubble columns to intensify the mixing properties[113,114].In addition,mechanical impellers can be applied in an IALR in combination with static mixers[96].The influences of agitation speed and static mixer length on the gas holdup,liquid velocity,and kLa were thoroughly investigated by Lu et al.[96].It was found that the kLa enhancement derived from the mechanical impeller was greater than that derived from the static mixers in a highly viscous fluid system.It should be noted that most of the studied mechanical impellers were motor-driven,requiring significant additional energy input.

    With an appropriate design of the shafts and impellers,a novel IALR inserted with constructed self-agitated impellers was established,and its hydrodynamics and mass transfer coefficient were explored[115].Fig.7[115]provides a schematic diagram of the IALR.The number and arrangement of impellers were carefully chosen in order to have them be agitated only by the gas throughput and liquid circulation.It was observed that the impellers in the riser of the IALR began rotating under an extremely low gas superficial velocity of 0.00629 m?s-1.With an increase of the superficial gas velocity,the rotating speed first increased and then tended to remain at a relatively constant value when the superficial gas velocity was above 0.45 m?s-1.With the installation of impellers,the gas holdup in the riser increased by about 45%.The enhancement of the gas holdup in the riser was 25%–48%greater than that induced by certain other internals,such as baffles[101]and perforated plates[100].In contrast,the enhancement of the overall gas holdup was insignificant,due to the offset effects in the downcomer.Based on this research,Lukic′et al.[17,116]built a novel EALR with self-agitated impellers inserted in the riser.They observed that the kLa was enhanced by up to 82% under a low superficial gas velocity(about 0.01 m?s-1),while the value of the kLa increased by 20%–30% under higher aerated conditions.The reduced improvement in the kLa when the superficial gas velocity increased was attributed to a decrease in impeller efficiency.

    2.2.5.General remarks

    Fig.7.Schematic diagram of an IALR equipped with self-agitated impellers[115].

    3.Process intensification of mixing and separation

    One of the major challenges in improving ALR performance is the efficient separation of solids from liquids in three-phase working systems.A considerable number of investigations have been conducted on separating solid catalyst particles from the liquid product to achieve high efficiency and low energy consumption.

    3.1.Settling

    Compared with the solid catalysts in fixed beds,which usually have a particle diameter greater than 1 mm,the solid particles in slurry reactors typically range from 10 to 200μm in diameter[25,26].The particles in the slurry can be easily separated by settling.Benham et al.[117]reported that a passive dynamic settler could be applied to divide the catalysts from a wax product in the downcomer of an EALR without a pump.The catalysts were delivered back to the column below the settler,with the light wax products being collected at the overflow.Unfortunately,the separation mechanism of the dynamic settler was not clearly presented.

    The solid particles in most slurry reactors have diameters ranging from 20 to 200μm[118,119]and can be cheaply separated using a hydrocyclone[120].Based on this fact,Yang et al.[25]proposed a continuous separation technique in an IALR,in which a compact,efficient,and energy-efficient hydrocyclone[121]was equipped in the bottom of the downcomer.Through this built-in separating device,clear liquid products were collected from the overflow of the hydrocyclone,with dense slurry catalysts concentrated at the underflow and continually circulated in the reactor.The division of clear liquid from the concentrated slurry in the hydrocyclone was accomplished by two driving forces:first,the hydraulic potential energy derived from the altitude difference between the liquid level and the hydrocyclone inlet;second,the kinetic energy of the slurry at the hydrocyclone entrance.It was shown that the distribution of solid particles inside the reactor varied slightly for four hours of operation in this novel slurry reactor,indicating that quite clean liquid products were obtained at the overflow,with few solid catalysts entrained.During this process,clear products were obtained,and a continuous catalytic reaction was ensured with no extra energy consumption.Therefore,a considerable amount of capital costs(including an expensive slurry pump,solids feeder,and pipes)and operating costs(including electric charges for the liquid–solid separation and solids feeding,and maintenance costs of moving devices)can be saved.Moreover,catalyst attrition can be reduced due to suspension by the continuous liquid phase,and catalyst deactivation due to a broken catalyst can be effectively prevented.It has been demonstrated that this new slurry reactor with an elaborate design functions very well when employing solid particles with an appropriate particle-size distribution and density.Process intensification of separation and mixing has also been applied in an EALR;Fig.8 illustrates the schematic diagram of this reactor.

    3.2.Filtration

    Filtration is one of the simplest and most convenient methods for separating solids from liquid flow.Rytter et al.[122]reported on a bubble column with a filtration element installed,which was free of catalysts,for separating the liquid products from the slurry.A constant-level controller was placed inside the filter unit to maintain a constant filtrate level below the slurry.The stable pressure differential between the slurry and the filtration device drove the slurry into the element for separation.Jager et al.[123]also set filters inside a bubble column with a back-flushed device to prevent blocking of the filtration medium.Anderson[124]described a separation method for a liquid product and solid catalyst using a novel internal microfilter with three parts:a porous reactor-side metal-cylinder outer surface,a porous reactor-side metal-cylinder inner surface,and a filter medium.Moreover,the outer and inner surfaces had a greater porosity than the filter medium.With this filter,catalysts with diameters from 0.5 to 100μm could be separated from wax products inside the slurry reactor.

    3.3.Composite methods

    In order to obtain a clear liquid product by separating solid particles from the slurry,Clerici and Belmonte[125]developed a hydrocyclone combined with single or multiple filtration units.After coarse separation by the hydrocyclone,the overflow containing fine particles was further separated by a micro/ultrafiltration element.The concentrated slurry from the underflow of the hydrocyclone and the filtration units flowed back to the reactor through the outboard pipeline.White et al.[126]described a method for separating a wax product from the catalysts by combining extraction with settling.The slurry was first delivered to a stationary mixer to extract the organic products and was then pumped to a centrifugal separator for further separation.An additional separator was highly desired for the substantial segregation of the catalyst particles.Hu et al.[127,128]combined settling and filtration for continuous solid–liquid separation in a bubble column and achieved high efficiency and long-term stable operation.

    Fig.8.Schematic diagram of an external loop slurry reactor integrating mixing and separation.

    Settling,which is induced by gravity or by centrifugal force,is one of the most common primary liquid–solid separation methods.Although filtration can be used for relatively complete separation,it often requires considerable energy input and introduces the problem of blocking.For the complete separation of a wax product and solid catalysts,a combination of different separation schemes is usually essential.Moreover,the separation method should be selected appropriately based on the chemical and physical properties of the slurry.

    4.Additional design considerations and challenges

    4.1.Flow regime

    Similar to a gas–liquid system,the hydrodynamics in industrial pneumatically agitated slurry reactors are characterized by three typical flow regimes:homogeneous(bubbly flow),transitional,and heterogeneous(churn-turbulent flow)regimes,mainly depending on the superficial gas velocity employed[8,26,87].The flow regimes in the reactor feature different bubble-size distributions,from which the flow patterns can be discriminated.In a homogeneous regime,the superficial gas velocity is relatively low and the bubble size is small(1–7 mm in diameter[76,129]);a narrow bubble-size distribution is manifested,and the coalescence and breakup phenomena are correspondingly negligible due to minor bubble–bubble interaction.At a higher gas velocity,some large bubbles(20–70 mm in diameter[76,129])are present due to bubble coalescence,and then a heterogeneous regime is developed;a wide bubble-size distribution is demonstrated,and coalescence and breakup phenomena are notable.These two regimes are bordered by a transition regime.

    It is vital to understand the hydrodynamic characteristics in slurry reactors under different flow regimes for the purposes of proper design,operation,control,and scale-up.In the last decades,many studies have been dedicated to the regime transition.Several methods for identifying the transition in the flow regime have been put forward,and details can be found in the related literature[8,26,87].It is notable that although there are publications on the regime transition for the gas–liquid two-phase flow,charts of the flow regime for the gas–liquid–solid three-phase flow are very scarce,especially for high solid concentrations.In a homogeneous regime,the gas holdup in the reactor increases linearly with an increase of the superficial gas velocity,and the line of the gas holdup versus the superficial gas velocity passes through the origin with a slope greater than or close to 1[8,130].However,for a heterogeneous flow regime,the slope is markedly suppressed due to large-bubble formation.Both the experimental and numerical results indicate that large bubbles mainly ascend in the central region of the column with large rise velocities,while small bubbles were much more uniform in the whole crosssection[26].Therefore,the radial profile of the gas holdup in the column is relatively uniform in a homogeneous regime,while it becomes parabolic in a heterogeneous regime.

    Gas-agitated slurry reactors can be operated in both the homogeneous flow regime and heterogeneous flow regime[131].In the heterogeneous regime,the small bubbles coalesce into large bubbles,which rise up at high velocities in a plug-flow manner.Under these circumstances,there is no axial mixing for the large bubbles,while the small bubbles have the same back-mixing characteristics as that in the slurry phase.Therefore,the heterogeneous flow regime has been thought to be the most optimal condition for the Fischer–Tropsch synthesis[132].

    It should be pointed out that the flow regime depends not only on the reactor design parameters(i.e.,distributor design and column diameter),but also on the operating parameters(i.e.,superficial gas and liquid velocity,temperature,and pressure)and the physical properties of the working system(i.e.,liquid viscosity,solid holdup,surface tension,density,and coalescing nature of the liquid phase)[133,134].It is noteworthy that when a poor gas distributor was used,heterogeneous regimes prevailed at all the superficial gas velocities[8].Van Baten and Krishna[135]found that when the solid particle concentrations exceeded 30 vol%,fast-rising large bubbles belonging to the spherical cap family occurred almost exclusively.

    4.2.Gas sparger

    The initial gas bubble size in a slurry reactor column is strongly affected by the gas sparger,which in turn influences the bubblesize distribution,gas holdup,and mass-transfer properties in the slurry reactor.Commonly applied types of gas sparger include single-orifice,sinter plate,perforated plate,porous plate,membrane,ring-type,and arm-type distributors.The spatial arrangement and orifice diameter of the gas sparger have been shown to influence the gas holdup in the reactor[136].It has been indicated that higher orifice numbers with smaller orifice diameters show better performance in mass transfer[137].It is noteworthy that bubble formation at multiple submerged orifices in a gas–liquid apparatus under industrial jetting conditions has been systematically investigated.A simple correlation of the influence of the orifice diameter,orifice superficial gas velocity,and liquid properties on the initial bubble diameter was proposed and was validated against extensive experimental data[138].

    For the design of a gas sparger,Lin et al.[139]investigated the radial distribution of the gas holdup and bubble size in EALRs equipped with gas distributors of porous plates and perforated plates.It was shown that different radial profiles of the gas holdup were obtained—that is,a wall-peaking profile was obtained for a column with a porous sinter plate distributor,and a core-peaking profile was obtained for a perforated plate distributor.Wei et al.[140]stated that a novel membrane-tube sparger with a bunch of porous cylindrical tubes enhanced the overall gas holdup by as much as 48.8% under a superficial gas velocity ranging from 0.004 to 0.04 m?s-1in an air–water system,with an improvement of the kLa of up to 84.4%.

    Aside from hydrodynamics and mass transfer,the flow structure in the column can be affected by the gas sparger design.Hooshyar et al.[141,142]adopted a structured gas distributor with a uniform needle sparger in order to structure the flow pattern and reduce back-mixing in the slurry bubble column reactors.It was shown that when applying this sparger,vortical structures were significantly reduced and the homogeneous regime was broadened.Vial et al.[12,143]investigated the flow regime in EALRs equipped with single-orifice,multiple-orifice,and porous-plate gas spargers through pressure fluctuations analysis.It was shown that the column with a single orifice always operated in the heterogeneous regime with a superficial gas velocity(Ug)ranging from 0.01 to 0.24 m?s-1.For the multiple-orifice sparger,the heterogeneous regime was achieved at a Ugof about 0.11 m?s-1.Nevertheless,the transition velocity was dependent on the start-up method for the column.A heterogeneous flow pattern in a column using porous plate with a‘‘dry”start-up was reached when the superficial gas velocity was about 0.11 m?s-1(i.e.,the same value by using a multiple-orifice distributor),whereas it was reached after about 0.07 m?s-1with a‘‘wet”start-up.

    The extent of the influence of the gas sparger on the hydrodynamics and transfer properties is related to the structural and operational conditions.Cao et al.[144]argued that sparger design had a noticeable effect on the gas holdup and axial dispersion at a low superficial gas velocity(Ug<0.025 m?s-1)with a low solid volume fraction(φs<2%),while slight effects were found at a high superficial gas velocity and high solid loading.It is the consensus that sparger design remarkably affects the hydrodynamics in a bubbly flow regime but only affects the region around the distributor in a churn-turbulent flow regime[145].Michele and Hempel[146]indicated that the sparger influence was more pronounced at low superficial gas velocities than at higher ones because large bubbles commonly formed at high superficial gas velocities regardless of the sparger type.

    Generally speaking,for porous-type gas distributors,the pores in the plate or tubes are usually designed to be large enough for the passage of gas yet small enough to keep solid particles from leaking through the pores.However,severe pore clogging may occur due to solid particle attrition or failure of the gas supply[26].For an industrial slurry reactor operating in the heterogeneous regime with a high superficial gas velocity and a high solid concentration,the ring-type distributor and arm-type distributor have been recommended,and some holes facing downward are highly desired to ensure a full solid suspension and clean discharge[147,148].When operating at a low superficial gas velocity,a gas sparger with small orifices,such as a porous-plate,perforatedplate,or membrane sparger,is preferred for high efficiency.It should be noted that the initial bubble size increases with an increase in solid volume fraction[149],and severe problems of mixing and mass/heat transfer may arise at high solid volume fractions[147].

    4.3.Influence of solid particles

    For industrial pneumatically agitated slurry reactors,the solid volume fraction in the column usually ranges from 15% to 30%,and can even go up to 50%[150].The presence of solid particles has a significant influence on the bubble size and distribution,bubble-rising velocity,gas holdup,and flow regime transition,which are pivotal parameters for mass-and heat-transfer properties[151–153].Although extensive studies have been conducted to investigate the effects of solid particles on the hydrodynamics and transfer efficiency in slurry reactors,the underlying mechanisms are still not clarified due to their complexity.

    For particles in the slurry that have a diameter smaller than 100μm,the solids can be easily suspended in the liquid phase without the problem of particle settling[26].Therefore,the slurry phase in the reactor can be simplified by assuming the liquids and solids to be in a pseudo-slurry phase,and the main influence of the solid concentration can be taken as an increase in the apparent viscosity of the liquid phase.In order to predict the hydrodynamics in slurry reactors using the apparent physical properties,numerous empirical and theoretical correlations of the effective viscosity have been proposed,as shown in Table 5[154–169].For a dilute suspension of solids,Einstein first proposed a correlation of the effective slurry viscosity based on the liquid viscosity and solid concentration;that is,μslurry=μl(1+2.5φs).However,this is only suitable for extremely low loading of solid particles(i.e.,φsis up to 0.01)because it neglects the hydrodynamic interactions.For a concentrated solids suspension in a slurry reactor,the relative viscosity(μslurry/μl)is commonly expressed as a polynomial function with solid concentration as an independent variable in order to extend the application range.The value of each coefficient,which reflects the rotation,collisions,and hydrodynamic interactions of the solids,is fitted using experimental data.Thus,different empirical correlations are achieved due to various operating conditions and physical properties[154].Furthermore,several phenomenal and theoretical relations have been proposed,with validation against experimental results[155,156].However,as pointed out byRabha et al.[170],the assumption of a pseudo-homogeneous slurry phase is an oversimplification,although it may be reasonable for smaller particles(ds<50μm)when Res<0.3 and the Stokes law assumption is valid.Krishna and Sie[76]also proposed that the assumption of a pseudo-slurry was valid only for the operation of large columns(i.e.,a diameter larger than 0.5 m)at high superficial gas velocities(Ug>0.2 m?s-1).

    Table 5 Theoretical and empirical correlations of the slurry’s apparent viscosity.

    Variations in the slurry’s apparent physical properties due to the solid presence influence the bubble formation,bubble size,and distribution in the reactor.Krishna et al.[50]observed that the gas voidage of small bubbles significantly decreased with an increase of the solid fraction from 0 to 0.36 due to the enhancement of bubble coalescence.Moreover,empirical correlations of gas holdup and bubble size in the churn-turbulent flow regime were proposed,based on statistics of experimental results[50].Chilekar et al.[171]proved that the relation of bubble size with slurry concentration was suitable when the solid volume fraction was up to 0.78.The bubble-size distribution has also been shown to be dependent on solid concentration[151].

    With an enhanced bubble size due to the solid existence,the gas holdup generally decreases with increasing solid concentration[25,172].Vandu et al.[173]also indicated that the flow transition velocity decreased and the homogenous regime narrowed with the presence of solids,due to the early formation of large bubbles.Nevertheless,the contrary effects of solids have been found under different operating conditions and solid physical properties.Kelkar et al.[174]showed that solid physical properties,size,and concentration had no significant influence on the gas holdup.Sada et al.[175]argued that the gas holdup decreased with an increased concentration of solids larger than 50μm,while it increased with the addition of solids smaller than 10μm.Jamialahmadi and Müller-Steinhagen[176]concluded that the gas holdup increased with wettable solids added but decreased with the addition of nonwettable solids.Mena et al.[177]clarified that the gas holdup first increased with increasing solid concentration and then decreased when the solid volume fraction was higher than 0.03.Correspondingly,the homogeneous regime stabilizes with a solid volume fraction of less than 0.03 and destabilizes when the solid volume fraction is larger than 0.03.Similar dual influences on the liquid velocity have also been demonstrated by Milivojevic et al.[178].The dual influence of solids on the aforementioned gas holdup and flow transition might result from the steric effect,physical properties,bubble formation and coalescence,and spatial distribution with the addition of solid particles in slurry reactors.Bubble diameter variations have been reported by Rabha et al.[151],who observed large bubbles with the addition of solid particles due to bubble coalescence but found that bubble slugs broke up at high solid volume fractions.In addition,the presence of suspended solids decreases the driving forces for liquid circulation and affects the flow pattern.It is noteworthy that the liquid circulation velocity can be influenced not only by the reactor geometry and superficial gas velocity but also by the flow regime and solids loading[6].

    Another noteworthy problem in the three-phase slurry reactor is the obvious axial solid concentration distribution with the highest concentration at the column base,which might result in local overheating.Li et al.[18]found that the axial difference of solid volume fraction can be as high as 76% only 1.4 m away from the column base with 5% solid loading at a superficial gas velocity of 0.03 m?s-1in a two-stage IALR.The non-uniformity of the solid holdup was thought to be caused by advection and circulations of buoyant clusters formed randomly.Murray and Fan[179]argued that the solid axial distribution tended to be more uniform with increasing superficial gas velocity and decreasing particle size in batch and continuous slurry columns.Zhang[180]investigated the axial distribution of the solid holdup in both tapped and cylindrical slurry columns.It was found that the uniformity of the solid axial distribution increased with an increasing superficial gas velocity ranging from 0.02 to 0.28 m?s-1and solid loading varying from 53 to 159 kg?m-3.Meanwhile,decreasing the particle diameter(70–180μm)and static height(0.6–2.4 m)of the slurry also enhanced the homogeneity of the solid axial distribution.Based on experimental results,an empirical correlation of the dimensionless Pelect number(Pe)was obtained through the least-squares method in order to predict the axial distribution of both monodispersed particles and binary mixtures.

    4.4.Other concerns

    It has been widely accepted that the hydrodynamics in a slurry reactor is independent of the column size and sparger layout,if D,H,and H/D are above certain critical values.Although different critical values have been proposed[26,181,182],the following criteria were recommended for the sake of safe design:①a diameter greater than 0.2 m;②an aspect ratio(H/D)greater than 5;③a sparger hole size greater than 1–2 mm[182].In slurry reactors,the separation problem is aggravated when fine particles are produced by either mechanical or chemical attrition;foam formation is another annoying problem that should be prevented[76].

    Xing et al.[134]reported that when the liquid viscosity was less than 10 mPa?s,the influence of liquid viscosity on the bubble breakup rate was negligible,while for a liquid viscosity exceeding 10 mPa?s,the bubble breakup rate decreased with an increase of liquid viscosity.The superficial liquid velocity had little influence on the gas holdup and the kLa[26].Yang et al.[149]argued that the influence of pressure on the initial bubble diameter was significant under relatively low pressure(below 2.0 MPa),while the influence of pressure on the initial bubble diameter became negligible under high pressure.In addition,high pressure decreased the bubble size and delayed the flow regime transition from the dispersed bubbly flow to the coalesced regime.Wang et al.[26]found that the influence of temperature on the gas–liquid mass transfer was much more prominent than that on the gas voidage,due to a higher liquid diffusivity at high temperatures.Although qualitative conclusions have been drawn,the quantitative results and underlying mechanism involved should be investigated systematically in the future.

    5.Computational fluid dynamics modeling

    Thanks to rapid developments in mathematics and highperformance computing techniques,CFD has been successfully employed as an efficient tool to investigate the hydrodynamics and transfer properties in gas-agitated slurry reactors[183–185].The Eulerian–Eulerian model[186],in which all phases are treated as interpenetrating continua,has been widely applied to predict the hydrodynamics in slurry reactors due to its low computational cost.Besides the dispute over turbulence[8],another controversy is the closure of the interphase momentum transfer,which is usually closed by drag,turbulent dispersion,lift,virtual mass,and wall lubrication forces.However,various theoretical and empirical models for these interphase forces have been proposed[187–189].Thus far,although many researchers have contributed to ensuring the most suitable closure correlations in their simulation cases[190],there is no definite single set of closure models for reliable simulation of the hydrodynamics and mass/heat transfer in multiphase slurry reactors due to their diversified configurations and wide range of operating conditions.

    In order to improve the computational accuracy and reduce its expense,different numerical schemes and algorithms have been proposed,as summarized in Table 6[112,191–201].Huang et al.[191]simulated the overall gas holdup and mass transfer in an IALR by applying a steady-two-dimensional axisymmetric scheme with an improved decoupling algorithm and specially treated outlet boundary conditions[192].Oey et al.[193]treated the threephase system as a pseudo-two-phase one,in which the dispersed solids were included in the liquid slurry phase,by taking into account solid physical properties such as density and viscosity.The solid volume fraction was computed by discretizing the mass-balance equation,in which the solid velocity was the sum of the liquid velocity,slip velocity due to settling,and turbulent dispersion term.It was confirmed that the predicted gas fraction at different gas flow rates aligned with that in the literature.By adopting this pseudo-two-phase model,Wang et al.[194]successfully simulated the local gas holdup and bubble-rise velocity in an EALR.

    Constant bubble size with no breakup or coalescence has been commonly employed when calculating the interfacial forces in the balance equations of CFD modeling to reduce the computational cost[202].However,bubble breakup and the coalescence phenomenon frequently occur in a column,resulting in a wide variation of bubble size and velocity,especially in transition andturbulent regimes.To remedy this situation,a population balance model(PBM)with a description of bubble coalescence and breakup is incorporated in the CFD scheme to predict the bubble size distribution in the column.The bubble breakup model proposed by Luo and Svendsen[203]and the bubble coalescence model proposed by Prince and Blanch[204]have been verified as being able to reliably describe the bubble distribution in aerated reactors[195,196].It was first reported by Lo[205]that the hydrodynamics in a bubble column can be well predicted with the CFD-PBM model combined with the developing multi-size group model(MUSIG),in which the bubbles are assumed to have the same motion velocity,for a moderate calculation cost.Frank et al.[197]further developed the inhomogeneous MUSIG,which divided the dispersed phase size into several groups regarding the mass and momentum balance equations.They found that it was sufficient to capture the hydrodynamics in the bubble columns with typically 3–5 groups.Considerable investigations have continually been conducted to obtain a more accurate prediction of the bubble-size distribution.To validate simulation results with experimental data,Chen et al.[196]had to increase the value of the bubble breakage rate by ten times so that the calculated breakage rate from the closure functions could match the experimental data.Bhole et al.[206]suggested that the coalescence kernel should be corrected by a coefficient related to the Stokes number,in which the slip velocity between the liquid and the bubble was considered.Yang and Xiao[198]coupled the energy-minimization multi-scale model with PBM equations to correct the bubble coalescence rate and obtained good agreement with the experimental data.

    Table 6 Numerical investigations of pneumatically agitated slurry reactors.

    With numerous studies having developed the mathematical model and discretization scheme,CFD modeling has been widely applied to predict the hydrodynamics and transfer features of gas–liquid and gas–liquid–solid flows,as well as the structural and operational optimization of pneumatically agitated reactors.Lehr et al.[199]successfully predicted the bubble-size distribution and flow regimes in a bubble column using a CFD-PBM model,while treating both large and small bubbles as pseudocontinuous phases in the liquid phase.Lestinsky et al.[201]explored the effects of the geometrical parameters of IALR—that is,the inner diameter,outer diameter,and equipping location of the draft tube—on the liquid circulation velocity and mixing time.Ni et al.[200]performed both two-and three-dimensional unsteady modeling on the flow patterns in an oscillatory baffled column and obtained reasonable agreement with the experimental data using digital particle image velocimetry measurements.Moraveji et al.[112]modeled the intensification of turbulence,hydrodynamics,and mass-transfer properties in a packed bed using a two-fluid turbulent model.They also provided a clear description of the flow patterns,gas holdup distribution,and flow regime transition in the column.Thus,to intensify the hydrodynamics and transfer in a pneumatically agitated reactor,CFD modeling has become an efficient and powerful tool for design,optimization,and scale-up.

    6.Conclusions and perspectives

    Pneumatically agitated slurry reactors(including bubble columns and ALRs)have a promising future due to their outstanding advantages,which include easy heat removal,excellent temperature control/removal for highly exothermic reactions,low production costs,high interphase mass-transfer rates,and a high online factor for catalyst addition and withdrawal[148,207,208].Mixing is better in ALRs than in bubble column reactors,due to the intensified well-defined fluid flow.In terms of mass-transfer capability,bubble column reactors are marginally superior to ALRs,due to a relatively higher gas holdup under the same conditions[209,210].

    For slurry reactors,challenges such as considerable back-mixing in both the continuous phase and the dispersed gas phase,low volumetric catalyst loading,bubble coalescence,and difficulties in scaling-up remain and should be considered[148].Meanwhile,the distribution of solids in slurry reactor has been experimentally verified to be non-uniform,especially in the axial direction.This may result in different catalytic reaction rates in different regions,and hot spots and runaways may even occur.It is noteworthy that a margin of temperature between the top and bottom is still observed in slurry reactors[211].Moreover,the catalyst deactivation rate and attrition may increase with an increase in slurry concentration.Furthermore,foaming,coking,and solid sediment can become severe problems in some applications.Therefore,despite the enormous advantages of pneumatically agitated slurry reactors over multitubular fixed-bed reactors,the slurry reactor technology should be systematically evaluated using an economically and industrially feasible method.

    Process intensification is one of the most commonly used methods to improve performance or solve problems in traditional slurry reactors.There are two main categories of process intensification.Methods in the first category involve intensifying the mixing and mass/heat transfer in the slurry reactor and include the use of internals(i.e.,tubes,perforated plates,static mixers,and structuring),vibrating excitement,and combination methods for bubble columns,and the use of baffles,perforated plates,static mixers or a packed bed,and mechanical internals for ALRs.It is noteworthy that enhancing the hydrodynamics and transfer in a slurry reactor is highly dependent on the design and operating conditions.Operating a churn-turbulent flow regime at high solid loading and high superficial gas velocity is usually preferred in order to obtain a high space–time yield.To remove the high exothermic heat of reaction efficiently,large amounts of equipment for heat removal are very desirable,and the hydrodynamics in the slurry may be correspondingly affected to a significant degree.The second category involves combining the mixing and the solid–liquid separation within the same slurry reactor in order to realize continuous production and reduce cost.The new technology of integrating directional flow in IALRs with low-cost solid–liquid separation through a hydrocyclone is one of the most competitive methods.For this novel technology,a compact and efficient hydrocyclone with a low pressure drop is one of the decisive factors[212].It should be noted that the separation of very fine particles from liquid products usually requires multiple separation methods.

    Although great endeavors have been made in this field,it remains difficult to design and scale-up industrial slurry reactors due to limited experimental information on hydrodynamics and mass transfer over a wide range of operating conditions.CFD has been verified as a powerful tool for investigating the hydrodynamics and transfer properties in pneumatically agitated slurry reactors.The CFD-PBM model is regarded as the most promising model for numerical simulations.However,it has been widely accepted that current simulations cannot predict the bimodal bubble-size distribution in a heterogeneous regime,and different bubble coalescence and breakup models may even give different bubble-size distributions[8,134].Therefore,CFD models should be developed and validated against experimental data.Otherwise,they will be dangerous,and the predicted results will be more art than science[213].

    Nomenclatures

    Acknowledgements

    This work was supported by the National Key Research and Development Program of China(2016YFB0301701),the National Natural Science Foundation of China(21808234,21878318,and 21938009),the DNL Cooperation Fund,Chinese Academy of Sciences(CAS)(DNL201902),the Strategic Priority Research Program of the CAS(XDA21060400),the QIBEBT and Dalian National Laboratory for Clean Energy of the CAS(QIBEBT ZZBS201803 and QIBEBT I201907),and the CAS Key Technology Talent Program.

    Compliance with ethics guidelines

    Shujun Geng,Zai-Sha Mao,Qingshan Huang,and Chao Yang declare that they have no conflict of interest or financial conflicts to disclose.

    欧美日韩国产mv在线观看视频| 如日韩欧美国产精品一区二区三区 | 午夜激情久久久久久久| 中国三级夫妇交换| 久久久久视频综合| 我的老师免费观看完整版| 国精品久久久久久国模美| 美女国产视频在线观看| 视频在线观看一区二区三区| 三级国产精品欧美在线观看| 久久热精品热| 22中文网久久字幕| 中文欧美无线码| 国产精品一国产av| 99热国产这里只有精品6| 成人国产av品久久久| 日日摸夜夜添夜夜添av毛片| 日本91视频免费播放| av线在线观看网站| av线在线观看网站| 建设人人有责人人尽责人人享有的| 久久精品夜色国产| 国产精品99久久99久久久不卡 | 极品少妇高潮喷水抽搐| 久久久a久久爽久久v久久| 国产爽快片一区二区三区| 飞空精品影院首页| 免费黄网站久久成人精品| 丝袜喷水一区| 日韩一区二区视频免费看| 国产色爽女视频免费观看| 亚洲欧美日韩另类电影网站| 青青草视频在线视频观看| 久久99热这里只频精品6学生| 久久久精品94久久精品| 久久人人爽av亚洲精品天堂| 丝袜美足系列| 欧美日韩av久久| 久热久热在线精品观看| 99热这里只有精品一区| 亚洲欧美色中文字幕在线| 在线观看一区二区三区激情| 在线观看人妻少妇| 久久国产精品男人的天堂亚洲 | 国产乱来视频区| 免费少妇av软件| 九九在线视频观看精品| 免费久久久久久久精品成人欧美视频 | 精品久久国产蜜桃| 欧美日韩av久久| 精品一品国产午夜福利视频| 亚洲精品日韩在线中文字幕| 91久久精品国产一区二区三区| 国产精品一区www在线观看| av在线app专区| 99热这里只有是精品在线观看| 亚洲精品一二三| 亚洲精品中文字幕在线视频| 亚洲国产精品成人久久小说| 久久青草综合色| 国产高清有码在线观看视频| av不卡在线播放| 久久99热6这里只有精品| 亚洲人成网站在线播| 国产探花极品一区二区| 欧美bdsm另类| 亚洲精品国产av成人精品| 亚洲少妇的诱惑av| 国产精品久久久久久久电影| 午夜激情久久久久久久| 天堂俺去俺来也www色官网| 91精品一卡2卡3卡4卡| 日韩av在线免费看完整版不卡| 国产极品天堂在线| 高清黄色对白视频在线免费看| 久久久久久久久久久免费av| 久久毛片免费看一区二区三区| 97超碰精品成人国产| 看免费成人av毛片| 亚洲欧美一区二区三区黑人 | 免费黄网站久久成人精品| 九草在线视频观看| 色婷婷av一区二区三区视频| 午夜免费男女啪啪视频观看| 欧美+日韩+精品| 亚洲色图综合在线观看| 亚洲伊人久久精品综合| 国产精品三级大全| av国产久精品久网站免费入址| 毛片一级片免费看久久久久| 99九九线精品视频在线观看视频| 18禁在线无遮挡免费观看视频| 国产视频内射| 乱码一卡2卡4卡精品| 国产精品99久久久久久久久| 少妇高潮的动态图| 视频在线观看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 国产成人免费观看mmmm| 免费不卡的大黄色大毛片视频在线观看| 青青草视频在线视频观看| 在现免费观看毛片| 久久精品国产鲁丝片午夜精品| 亚洲精品一二三| 一级毛片黄色毛片免费观看视频| 亚洲欧洲国产日韩| 一区二区三区四区激情视频| 韩国高清视频一区二区三区| 91成人精品电影| 午夜视频国产福利| 免费观看av网站的网址| av.在线天堂| 中文天堂在线官网| 美女主播在线视频| 一区二区三区精品91| 日韩中文字幕视频在线看片| 亚洲国产av新网站| av福利片在线| 美女国产视频在线观看| 久久精品国产自在天天线| 亚洲美女视频黄频| 日韩大片免费观看网站| 免费av中文字幕在线| 3wmmmm亚洲av在线观看| 亚洲国产av影院在线观看| 国产男女内射视频| 精品少妇内射三级| 亚洲欧美日韩卡通动漫| 如何舔出高潮| 国产又色又爽无遮挡免| 成年美女黄网站色视频大全免费 | 成人免费观看视频高清| 91午夜精品亚洲一区二区三区| 五月伊人婷婷丁香| 十分钟在线观看高清视频www| 人人澡人人妻人| av国产精品久久久久影院| 午夜老司机福利剧场| 狂野欧美激情性bbbbbb| 日韩 亚洲 欧美在线| 亚洲精品乱码久久久v下载方式| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线看a的网站| 免费观看a级毛片全部| 日韩三级伦理在线观看| 久久精品久久精品一区二区三区| 国产免费视频播放在线视频| 在线免费观看不下载黄p国产| 国产一级毛片在线| 97精品久久久久久久久久精品| 特大巨黑吊av在线直播| 99久国产av精品国产电影| 蜜桃久久精品国产亚洲av| 国产精品久久久久久精品古装| 狂野欧美激情性bbbbbb| 成年人午夜在线观看视频| 人人妻人人添人人爽欧美一区卜| 男女边吃奶边做爰视频| av一本久久久久| 国产精品久久久久久av不卡| 国产精品熟女久久久久浪| 人妻少妇偷人精品九色| 亚洲人与动物交配视频| 久久久久网色| 亚洲国产欧美日韩在线播放| 国产日韩欧美在线精品| 日韩成人伦理影院| 亚洲av欧美aⅴ国产| 亚洲精品,欧美精品| 日本猛色少妇xxxxx猛交久久| 2021少妇久久久久久久久久久| 下体分泌物呈黄色| 亚洲精品一区蜜桃| 亚洲精品一二三| 纵有疾风起免费观看全集完整版| 精品久久久久久久久亚洲| 亚洲精品国产色婷婷电影| 国产成人aa在线观看| 精品人妻在线不人妻| 大陆偷拍与自拍| 欧美精品一区二区免费开放| 亚洲av福利一区| 久久国产亚洲av麻豆专区| 18禁动态无遮挡网站| 午夜福利,免费看| 男女免费视频国产| 日本午夜av视频| 777米奇影视久久| 男女国产视频网站| 插阴视频在线观看视频| 一级黄片播放器| 国产成人精品一,二区| 亚洲在久久综合| 国产永久视频网站| 国产伦精品一区二区三区视频9| av天堂久久9| 中文字幕精品免费在线观看视频 | 午夜激情久久久久久久| 卡戴珊不雅视频在线播放| 欧美少妇被猛烈插入视频| 中国国产av一级| 男人爽女人下面视频在线观看| 亚洲国产欧美日韩在线播放| 伊人久久精品亚洲午夜| 国产免费福利视频在线观看| 五月开心婷婷网| 欧美xxⅹ黑人| 夫妻午夜视频| 亚洲欧美中文字幕日韩二区| 如何舔出高潮| 免费黄色在线免费观看| 国产老妇伦熟女老妇高清| 国产老妇伦熟女老妇高清| 美女国产高潮福利片在线看| 热re99久久精品国产66热6| 精品视频人人做人人爽| 国产不卡av网站在线观看| 国产欧美亚洲国产| 国产黄频视频在线观看| 我的女老师完整版在线观看| 少妇猛男粗大的猛烈进出视频| 自拍欧美九色日韩亚洲蝌蚪91| 爱豆传媒免费全集在线观看| 免费av中文字幕在线| 国产av国产精品国产| 国产视频内射| 中文字幕久久专区| 国产高清不卡午夜福利| 制服诱惑二区| 777米奇影视久久| 精品少妇久久久久久888优播| 美女cb高潮喷水在线观看| 超碰97精品在线观看| 欧美变态另类bdsm刘玥| 精品一品国产午夜福利视频| 女性生殖器流出的白浆| 天天操日日干夜夜撸| 99热国产这里只有精品6| 亚洲av福利一区| 亚洲,一卡二卡三卡| 免费黄色在线免费观看| 老熟女久久久| 99国产精品免费福利视频| 全区人妻精品视频| 日本av免费视频播放| 午夜久久久在线观看| 中文字幕av电影在线播放| 久久精品夜色国产| 午夜福利视频在线观看免费| 亚洲久久久国产精品| 日韩人妻高清精品专区| 亚洲av男天堂| 日日爽夜夜爽网站| 狂野欧美白嫩少妇大欣赏| 久久久久久久国产电影| av视频免费观看在线观看| 伦理电影大哥的女人| 久久精品久久久久久噜噜老黄| 亚洲国产最新在线播放| 在线观看一区二区三区激情| 高清黄色对白视频在线免费看| 少妇丰满av| 91久久精品国产一区二区三区| 高清欧美精品videossex| 国产 精品1| 18禁动态无遮挡网站| 亚洲精品一二三| 成年女人在线观看亚洲视频| 成人免费观看视频高清| 超色免费av| 一级a做视频免费观看| 中文字幕精品免费在线观看视频 | 久久青草综合色| 麻豆精品久久久久久蜜桃| 国产精品麻豆人妻色哟哟久久| 中文字幕人妻丝袜制服| 久久久久精品久久久久真实原创| 黄片播放在线免费| 欧美日韩国产mv在线观看视频| 欧美亚洲日本最大视频资源| 赤兔流量卡办理| 这个男人来自地球电影免费观看 | 久久久久精品久久久久真实原创| 丰满迷人的少妇在线观看| 人妻系列 视频| 欧美激情国产日韩精品一区| 中文天堂在线官网| 精品国产一区二区久久| av天堂久久9| av女优亚洲男人天堂| 国产一区有黄有色的免费视频| 亚洲欧美成人精品一区二区| 亚洲欧美精品自产自拍| 精品久久久久久电影网| 人人澡人人妻人| 精品一区二区三区视频在线| 中文字幕免费在线视频6| 熟女av电影| 亚洲av国产av综合av卡| 午夜激情久久久久久久| 91精品伊人久久大香线蕉| 人妻少妇偷人精品九色| 黑人巨大精品欧美一区二区蜜桃 | 日本欧美视频一区| 99精国产麻豆久久婷婷| 一区二区av电影网| 国产免费视频播放在线视频| 一个人免费看片子| 国产成人精品在线电影| 少妇人妻 视频| 中文字幕av电影在线播放| 婷婷色av中文字幕| 伊人久久精品亚洲午夜| 91aial.com中文字幕在线观看| 激情五月婷婷亚洲| 九色亚洲精品在线播放| 成年女人在线观看亚洲视频| 三级国产精品欧美在线观看| 97在线视频观看| a 毛片基地| 伦理电影大哥的女人| 一级,二级,三级黄色视频| 多毛熟女@视频| 中文乱码字字幕精品一区二区三区| 男女无遮挡免费网站观看| 伊人亚洲综合成人网| 2018国产大陆天天弄谢| 老司机亚洲免费影院| 另类亚洲欧美激情| av播播在线观看一区| 午夜91福利影院| 国产精品麻豆人妻色哟哟久久| 永久网站在线| 亚洲av电影在线观看一区二区三区| 特大巨黑吊av在线直播| 免费看光身美女| 欧美亚洲 丝袜 人妻 在线| 91久久精品国产一区二区三区| 91在线精品国自产拍蜜月| av卡一久久| 亚洲av不卡在线观看| 伦精品一区二区三区| 美女国产视频在线观看| 亚洲欧美清纯卡通| 欧美老熟妇乱子伦牲交| 婷婷色麻豆天堂久久| 中文字幕av电影在线播放| 美女中出高潮动态图| 美女脱内裤让男人舔精品视频| 视频区图区小说| 免费观看无遮挡的男女| 国产一级毛片在线| 国产有黄有色有爽视频| 亚洲av电影在线观看一区二区三区| 亚洲精品中文字幕在线视频| 亚洲,欧美,日韩| 欧美人与性动交α欧美精品济南到 | 亚洲av欧美aⅴ国产| 九色成人免费人妻av| 女性被躁到高潮视频| 亚洲三级黄色毛片| 高清黄色对白视频在线免费看| 亚洲综合色网址| 国产亚洲一区二区精品| 女人久久www免费人成看片| 日韩中文字幕视频在线看片| a级毛色黄片| 国产精品一区二区三区四区免费观看| 国产精品人妻久久久影院| 丰满饥渴人妻一区二区三| 免费高清在线观看日韩| 日韩强制内射视频| av线在线观看网站| 高清黄色对白视频在线免费看| 日韩强制内射视频| 蜜桃久久精品国产亚洲av| 亚洲av欧美aⅴ国产| 国产精品一区二区在线不卡| 亚洲精品456在线播放app| 欧美日韩综合久久久久久| 又粗又硬又长又爽又黄的视频| 久久精品久久久久久久性| 老司机影院毛片| 国产精品人妻久久久影院| 五月天丁香电影| 国产在线视频一区二区| 精品一区二区免费观看| 成年人午夜在线观看视频| 亚洲性久久影院| av免费在线看不卡| 美女cb高潮喷水在线观看| 久久精品国产鲁丝片午夜精品| 满18在线观看网站| 色视频在线一区二区三区| 日日啪夜夜爽| 亚洲国产精品一区二区三区在线| 国产男女内射视频| 国产成人aa在线观看| 中国美白少妇内射xxxbb| 少妇被粗大猛烈的视频| 一区二区三区免费毛片| 中文字幕人妻熟人妻熟丝袜美| 在线天堂最新版资源| 九色成人免费人妻av| 亚洲精品成人av观看孕妇| 女性生殖器流出的白浆| 最近最新中文字幕免费大全7| 97精品久久久久久久久久精品| 久久女婷五月综合色啪小说| 黑人猛操日本美女一级片| 三级国产精品片| 久久久久久久久久成人| 亚洲精品日韩av片在线观看| 人妻人人澡人人爽人人| 精品亚洲成国产av| a级毛片在线看网站| 少妇猛男粗大的猛烈进出视频| 又黄又爽又刺激的免费视频.| 男女啪啪激烈高潮av片| 日日爽夜夜爽网站| 免费人成在线观看视频色| 在线观看免费日韩欧美大片 | 亚洲精品久久成人aⅴ小说 | 天堂8中文在线网| 亚洲国产毛片av蜜桃av| 欧美精品人与动牲交sv欧美| av国产精品久久久久影院| 一区二区日韩欧美中文字幕 | 丁香六月天网| 国产精品成人在线| 熟女人妻精品中文字幕| 国产在线视频一区二区| 一级片'在线观看视频| 亚洲国产欧美日韩在线播放| 国产免费一级a男人的天堂| 久久韩国三级中文字幕| 国产淫语在线视频| 少妇丰满av| 麻豆乱淫一区二区| 在线看a的网站| 亚洲精品av麻豆狂野| 国产亚洲av片在线观看秒播厂| 在线观看一区二区三区激情| 成人毛片60女人毛片免费| av天堂久久9| 日日啪夜夜爽| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看| av线在线观看网站| 亚洲国产精品999| 七月丁香在线播放| 久久精品国产亚洲网站| 夜夜看夜夜爽夜夜摸| 日韩成人av中文字幕在线观看| 一个人免费看片子| 国产 一区精品| 99久久精品国产国产毛片| 国产成人精品福利久久| 高清av免费在线| 女的被弄到高潮叫床怎么办| 久久99一区二区三区| 中国三级夫妇交换| 国产片特级美女逼逼视频| 蜜桃在线观看..| 99久久中文字幕三级久久日本| 国产欧美日韩一区二区三区在线 | 日本av手机在线免费观看| 欧美变态另类bdsm刘玥| 免费不卡的大黄色大毛片视频在线观看| 国产精品一二三区在线看| 成人亚洲欧美一区二区av| 伦理电影免费视频| 国产高清三级在线| 亚洲精品日韩av片在线观看| 97精品久久久久久久久久精品| 久久人妻熟女aⅴ| 国产精品久久久久久av不卡| 大又大粗又爽又黄少妇毛片口| 欧美人与性动交α欧美精品济南到 | 国产精品欧美亚洲77777| 国产精品久久久久久久久免| 91久久精品电影网| 国产成人免费无遮挡视频| h视频一区二区三区| 一区二区av电影网| 韩国av在线不卡| 亚洲中文av在线| 丰满乱子伦码专区| 亚洲少妇的诱惑av| 婷婷色综合大香蕉| 简卡轻食公司| 在线天堂最新版资源| 黄色怎么调成土黄色| 久久久久久久亚洲中文字幕| 五月伊人婷婷丁香| 丝袜脚勾引网站| 欧美最新免费一区二区三区| 多毛熟女@视频| 久久国产精品男人的天堂亚洲 | 久久99一区二区三区| 精品少妇内射三级| 99国产精品免费福利视频| 国产精品免费大片| 久久久久久久久久久丰满| 亚洲精品,欧美精品| a级毛片黄视频| 久久国产精品大桥未久av| 天堂8中文在线网| 人妻一区二区av| 伦理电影免费视频| 天天影视国产精品| 精品一区二区免费观看| 成年av动漫网址| 国产精品国产三级专区第一集| 狂野欧美白嫩少妇大欣赏| 久久久国产精品麻豆| 亚洲情色 制服丝袜| 国产av一区二区精品久久| 亚洲国产精品专区欧美| 蜜桃国产av成人99| 国产精品偷伦视频观看了| 美女国产高潮福利片在线看| 插逼视频在线观看| 中文字幕av电影在线播放| 久久久久久久久久成人| av一本久久久久| 精品一区二区三卡| 99热网站在线观看| 国产一级毛片在线| 国产精品成人在线| 秋霞在线观看毛片| 亚洲精品aⅴ在线观看| 一本久久精品| 国产男女超爽视频在线观看| 91在线精品国自产拍蜜月| 国产伦理片在线播放av一区| 多毛熟女@视频| 一级毛片aaaaaa免费看小| 国产成人91sexporn| 婷婷色av中文字幕| 亚洲精品,欧美精品| 天堂俺去俺来也www色官网| 中文天堂在线官网| 国产精品熟女久久久久浪| 成人亚洲欧美一区二区av| 亚洲成人手机| 91精品三级在线观看| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 日日啪夜夜爽| 97超视频在线观看视频| 精品一区二区三卡| 亚洲综合色网址| 久热这里只有精品99| 国产片内射在线| 人人妻人人爽人人添夜夜欢视频| 中国国产av一级| 国产精品不卡视频一区二区| 成人漫画全彩无遮挡| 久热这里只有精品99| 五月玫瑰六月丁香| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 国产午夜精品一二区理论片| 国国产精品蜜臀av免费| 蜜桃在线观看..| 观看美女的网站| 人妻制服诱惑在线中文字幕| 午夜福利视频精品| 亚洲人成77777在线视频| 在线观看免费高清a一片| 人妻制服诱惑在线中文字幕| 丰满迷人的少妇在线观看| h视频一区二区三区| 在线观看三级黄色| 国产高清三级在线| 国产伦精品一区二区三区视频9| 久久99热6这里只有精品| 日韩人妻高清精品专区| 下体分泌物呈黄色| 在线看a的网站| 一本一本综合久久| 久久99一区二区三区| 人妻人人澡人人爽人人| 99国产综合亚洲精品| 毛片一级片免费看久久久久| 日韩熟女老妇一区二区性免费视频| 国产高清有码在线观看视频| 国产成人freesex在线| 观看美女的网站| 99re6热这里在线精品视频| 精品少妇久久久久久888优播| 在线精品无人区一区二区三| 亚洲精品av麻豆狂野| 最近中文字幕2019免费版| 欧美 亚洲 国产 日韩一| 99国产综合亚洲精品| 国产在视频线精品| 在线天堂最新版资源| 国产亚洲精品久久久com| 亚洲成人av在线免费| 国语对白做爰xxxⅹ性视频网站| 久久国产精品大桥未久av| 18禁在线无遮挡免费观看视频| 久久久久久久久久人人人人人人| 一区二区三区免费毛片| 久久女婷五月综合色啪小说| 一区二区三区免费毛片| 欧美日韩综合久久久久久| 一本久久精品| 永久网站在线| 亚洲欧美中文字幕日韩二区| 精品人妻偷拍中文字幕| 亚洲精品成人av观看孕妇|