基于EMMS模型的搅拌釜内气液两相流数值模拟

采用欧拉-欧拉模型对搅拌釜内气液两相流进行了三维CFD模拟,重点研究了采用不同曳力模型时CFD模拟对搅拌桨附近排出流区两相流动的预测能力。模拟结果表明CFD能准确地预测排出流区的液相速度分布,但采用传统的Schiller-Naumann曳力一定程度上低估了排出流区的气液相间曳力,导致在完全扩散区CFD预测的分布器和桨叶下方区域气含率偏小,而基于气液非均匀结构和能量最小多尺度（EMMS）方法得到的DBS-Global曳力模型能更准确地描述完全扩散区气液搅拌釜内流动情况。与传统曳力模型相比,采用DBS-Global曳力模型能显著提高对气含率的预测。     

Experimental determination and numerical simulation of mixing time in a gas-liquid stirred tank

In this work, mixing experiments and numerical simulations of flow and macro-mixing were carried out in a 0.24 m i.d. gas-liquid stirred tank agitated by a Rushton turbine. The conductivity technique was used to measure the mixing time. A two-phase CFD (computational fluid dynamics) model was developed to calculate the flow field, k and ε distributions and holdup. Comparison between the predictions and the reported experimental data [Lu, W.M., Ju, S.J., 1987. Local gas holdup, mean liquid velocity and turbulence in an aerated stirred tank using hot-film anemometry. Chemical Engineering Journal 35 (1), 9-17] of flow field and holdup at same conditions were investigated and good agreements have been got. As the complexity of gas-liquid systems, there was still no report on the prediction of mixing time through CFD models in a gas-liquid stirred tank. In this paper, the two-phase CFD model was extended for the prediction of the mixing time in the gas-liquid stirred tank for the first time. The effects of operating parameters such as impeller speed, gas flow rate and feed position on the mixing time were compared. Good agreements between the simulations and experimental values of the mixing time have also been achieved.     

Phase distributions in a gas–liquid–solid circulating fluidized bed riser

The distributions of the three phases in gas–liquid–solid circulating fluidized beds (GLSCFB) were studied using a novel measurement technique that combines electrical resistance tomography (ERT) and optical fibre probe. The introduction of gas into a liquid–solid circulating fluidized bed (LSCFB), thus forming a GLSCFB, caused the increase of solids holdup due to the significantly decreased available buoyancy with the lower density of the gas, even with a somewhat increased liquid velocity due to the decreased liquid holdup giving space for the gas holdup. The gas passed through the riser in the form of bubbles, which tended to flow more through the central region of the riser, leading to more radial non‐uniformity in radial holdup of the phases. The gas velocity has the most significant effect on the gas phase holdup. While the gas velocity also has an obvious effect to the solids holdups, the liquid flow rate had a much more considerable effect on the phase holdups. The solids circulation rate also had a significant effect on the phase holdups, with increasing solids circulation rate causing much more increased solids holdup in the central region than close to the wall. A correlation was developed for the relative radial distributions of solids holdup in GLSCFB, as such radial profiles were found similar over a wide range of operating conditions, like those in a typical gas–solid circulating fluidized beds (GSCFB). Finally, the axial solids profiles in a GLSCFB was found to be much closer to those in an LSCFB which are very uniform, than those found in a GSCFB which are less uniform and sometime having a S shape. Water was used as the continuous and conductive phase, air was the gas phase and glass bead and lava rock particles were used as the solid and non‐conductive phase.     

γ-CT measurement and CFD simulation of cross section gas holdup distribution in a gas–liquid stirred standard Rushton tank

Cross section gas holdup distributions at 3/4 dimensionless static liquid height in a gas–liquid stirred standard Rushton tank were measured using ¹³⁷Cs γ-CT scan measuring technology at larger gas flow rates and higher impeller rotating speeds. The obtained CT scan images and digital distribution curves of gas holdup with dimensionless radius based on the CT images could explain the fluctuation changes of gas holdup distribution. The dense area of gas holdup distribution appeared in the upper space of impeller blades. Gas holdup increased both with gas flow rate and impeller rotating speed, but gas flow rate had more influence on gas holdup than impeller rotating speed. The Eulerian–Eulerian two-fluid model coupling with the bubbles' coalescence and break-up models, and the drag coefficient model were established to make CFD simulation of gas holdup distributions for the gas–liquid stirred Rushton tank under different gas flow rates and impeller rotating speeds.     

Local Turbulent Energy Dissipation Rate in an Agitated Vessel: Experimental and Turbulence Scaling

The hydrodynamics and the flow field in an agitated vessel were measured using 2-D time resolved particle image velocimetry (2-D TR PIV). The experiments were carried out in fully baffled cylindrical flat bottom vessels 300 and 400 mm in inner diameter. The 300 mm inner diameter tank was agitated by a Rushton turbine 100 mm in diameter, and the 400 mm inner diameter tank was agitated by a Rushton turbine 133 mm in diameter. Three liquids of different viscosities were used as the agitated liquid: (i) distilled water (ν = 9.35 × 10^–7 m²/s), (ii) a 28 vol % aqueous solution of glycol (ν = 2 × 10^–6 m²/s), and (iii) a 43 vol % aqueous solution of glycol (ν = 3 × 10^–6 m²/s). The velocity fields were measured at an impeller rotation speed in the range from 300 to 850 rpm, which covers the Reynolds number range from 50000 to 189000. This means that fullydeveloped turbulent flow was reached. The experiments were performed to investigate the applicability of the following relations: ε* = ε/(u⁴/ν) = const, vK/u = const, Λ/ηK = const, τ_Λ/τ_K = const, ε* = ε/((Nd)4/ν) = const, Λ/d ∝ Re^–1, ηK/d ∝ Re^–1, vK/(Nd) = const, Nτ_Λ ∝ R^–1, Nτ_K ∝ Re^–1, and ε/(Nq) ∝ Re. These formulas were theoretically derived in our previous work, using turbulence theory, in particular, using turbulence spectrum analysis. The correctness of the proposed relations is investigated by statistical hypothesis testing.     

Hydrodynamics of a slurry airlift reactor at high solid concentrations

The hydrodynamics of a slurry airlift reactor at high solid concentrations were experimentally studied. The influences of the average solid concentration, superficial gas velocity and particle size on the radial and axial profiles of the solid holdup, average gas holdup and liquid circulation velocity were investigated. The local solid holdup was measured with an electrical conductivity probe. At low solid concentrations or high superficial gas velocities, the radial profile of the solid holdup was uniform. At high solid concentrations, the radial profile of the solid holdup was nonuniform, with higher values near the wall. This radial nonuniformity increased with decreased superficial gas velocity or increased average solid concentration. The axial profile of the cross-sectional average solid holdup was uniform at all conditions in this work, even at high solid concentrations. The average gas holdup and liquid circulation velocity increased with the superficial gas velocity but decreased with the average solid concentration. A mathematical model based on the balance of the transverse lift force and turbulent dispersion force was proposed to predict the radial profile of the solid holdup. Reasonable predictions were obtained from this model with an adjustable model parameter.     

鼓泡床反应器内流动与传质行为的研究进展

总结了有关浆态鼓泡床反应器内流动、混合用气液传质特性的研究成果,详细地介绍了鼓泡床反应器内气含率、液速、液体轴向扩散系数、传质系数的测量方法,阐述了鼓泡床反应器性能的主要影响因素,如系统压力、温度、气体表观气速、液体性质及固含率等对流动、液相混合和传质特性的影响,并对鼓泡床反应器的应用前景进行了详述.     

Experimental and numerical investigations of scale-up effects on the hydrodynamics of slurry bubble columns

Experiments and simulations were conducted for bubble columns with diameter of 0.2 m(180 mm i.d.), 0.5 m(476 mm i.d.) and 0.8 m(760 mm i.d.) at high superficial gas velocities(0.12–0.62 m·s-1) and high solid concentrations(0–30 vol%). Radial profiles of time-averaged gas holdup, axial liquid velocity, and turbulent kinetic energy were measured by using in-house developed conductivity probes and Pavlov tubes. Effects of column diameter, superficial gas velocity, and solid concentration were investigated in a wide range of operating conditions. Experimental results indicated that the average gas holdup remarkably increases with superficial gas velocity, and the radial profiles of investigated flow properties become steeper at high superficial gas velocities. The axial liquid velocities significantly increase with the growth of the column size, whereas the gas holdup was slightly affected. The presence of solid in bubble columns would inhibit the breakage of bubbles, which results in an increase in bubble rise velocity and a decrease in gas holdup, but time-averaged axial liquid velocities remain almost the same as that of the hollow column. Furthermore, a 2-D axisymmetric k–ε model was used to simulate heterogeneous bubbly flow using commercial code FLUENT 6.2. The lateral lift force and the turbulent diffusion force were introduced for the determination of gas holdup profiles and the effects of solid concentration were considered as the variation of average bubble diameter in the model. Results predicted by the CFD simulation showed good agreement with experimental data.     

Influence of impeller diameter on overall gas dispersion properties in a sparged multi-impeller stirred tank☆

The impeller configuration with a six parabolic blade disk turbine below two down-pumping hydrofoil propellers, identified as PDT + 2CBY, was used in this study. The effect of the impeller diameter D, ranging from 0.30T to 0.40T (T as the tank diameter), on gas dispersion in a stirred tank of 0.48 m diameter was investigated by experimental and CFD simulation methods. Power consumption and total gas holdup were measured for the same impeller configuration PDT + 2CBY with four different D/T. Results show that with D/T increases from 0.30 to 0.40, the relative power demand (RPD) in a gas–liquid system decreases slightly. At low superficial gas velocity V_S of 0.0078 m·s^-1, the gas holdup increases evidently with the increase of D/T. However, at high superficial gas velocity, the systemwith D/T=0.33 gets a good balance between the gas recirculation and liquid shearing rate, which resulted in the highest gas holdup among four different D/T. CFD simulation based on the two-fluid model along with the Population Balance Model (PBM) was used to investigate the effect of impeller diameter on the gas dispersion. The power consumption and total gas holdup predicted by CFD simulation were in reasonable agreement with the experimental data.     

连续内环流三相反应器局部流动特性

在φ200 mm×2500 mm连续内环流三相反应器内,考察了空气 水 玻璃珠体系反应器内局部流动参数随操作条件的变化规律。结果表明,导流筒内截面平均气含率随表观气速的增大而增大,较之气液两相流,在低固含率时,加入固体对气含率影响不明显,而在较高固含率下,气含率有明显降低,但固体再增加时对气含率变化影响不大。在较低表观气速下,进料浆速对导流筒内气含率轴向分布趋势有一定的影响,但在较高表观气速下影响不大,导流筒内的气含率大于环隙内的气含率且随气速增大差别更加明显,浆相连续有利于气相分散并增大环隙内的气含率。导流筒内循环浆速径向分布呈抛物状,中心高、近壁处低,受进料浆速和入口固含率影响都不大。浆相循环强度最低为20,高可达180。固含率轴、径向分布受表观气速和进料浆速的影响,固含率轴、径向分布基本均匀,随进料浆速增加,反应器内固含率降低。     

Flow behavior of wakes in a three-phase slurry bubble column with viscous liquid medium

Flow behavior of wakes has been investigated in a three-phase slurry bubble column of 0.102 m ID and 1.5 m in height. The dependence of wake characteristics such as rising velocity, frequency, holdup and equivalent size on the operating variables was examined by employing an electric resistivity probe method. The gas velocity, liquid viscosity and solid content in the slurry phase were chosen as independent parameters. The rising velocity of wake region increased with an increase in the gas velocity (4.0–12.0 cm/s), liquid viscosity (1.0–50.0 mPa·s) or solid content (0–25 wt%) in the slurry phase. The frequency and holdup of wake phase increased with increasing gas velocity, but decreased with increasing liquid viscosity or solid content in the slurry phase. The equivalent size of wake phase increased with increasing gas velocity, liquid viscosity or solid content in the slurry phase. The wake properties and holdup were well correlated with operation variables within these experimental conditions.     

竖直上升气液泡状流数学模型封闭的研究

竖直上升管气液两相流广泛应用于相变传热、核反应堆等工业过程。本文以竖直上升气液两相流为研究对象,运用欧拉双流体模型,针对表观液速为0.45m/s、表观气速分别为0.015m/s和0.1m/s的泡状流数值模拟过程中的升力、壁面润滑力、湍流扩散力、气泡诱导湍流（BIT）等封闭模型,开展数值模拟比较研究。模拟发现：①低气速泡状流中,升力和壁面润滑力的同时加入能够改善壁面附近的气含率,气泡在这两个力作用下在径向上达到一个相对平衡,得到与实验气含率类似的壁面峰,模拟的液相速度较合理;低气速时,BIT的影响可以忽略。②高气速泡状流中,BIT对气-液两相流的模拟结果影响比较明显,湍动耗散源项的加入能使液速分布的模拟结果得到改善,Troshko模型相对Sato模型更能反映气泡诱导湍流对液相湍流的作用。③高气速时升力的引入使气含率产生壁面峰,加入湍流扩散力能使峰值略微降低,但仍没有解决高气速时引入升力出现的气含率壁面峰问题,说明在径向上湍流扩散力还不足以抵抗升力。     

Experimental analysis and development of correlations for gas holdup in high pressure slurry co-current bubble columns

The effect of liquid and gas velocities, solid concentrations, and operating pressure has been studied experimentally in a 15 cm diameter air-water-glass beads bubble column. The superficial gas and liquid velocities varied from 1.0 to 40.00 cm/s and 0 to 16.04 cm/s, respectively, while the solid loading varied from 1 to 9%. The gas holdup in the column was reduced sharply as we switched from batch to co-current mode of operation. At low gas velocity, the effect of liquid velocity was insignificant; while at high gas velocity, increasing liquid velocity decreased the gas holdup. Drift flux approach was applied to quantify the combined effect of liquid and gas velocities over gas holdup. For co-current three phase flows, the gas holdup decreased with increase in solid loading for all pressures. But for batch operations, when solid loading was 5% or more, settling started leading to higher gas holdup. Increasing pressure from atmospheric conditions increased the gas holdup significantly, flattening asymptotically.     

气体分布方式对带列管内构件的鼓泡塔流动规律的影响

在带列管内构件的鼓泡塔内测量了4种不同布气方式下的气含率和液速径向分布,并与无列管内构件的空塔中的分布进行了比较。结果表明：中心布气条件下气含率与液速的径向分布比空塔更为陡峭;环隙及近壁布气时呈现出环隙高、两边低的马鞍形分布;均匀布气时径向分布较空塔更为平坦。空塔内气体分布器的影响是局部性的,充分发展段在塔内占主要部分;而在列管塔中气体分布器的影响是全局性的,气含率与液速的初始分布决定着其全塔分布。在带列管的大型鼓泡塔中难以观察到充分发展段的存在,因此,气体分布器的设计具有比空塔更为重要的意义     

Experimental investigation on multiscale hydrodynamics in a novel gas–Liquid–Solid three phase jet-Loop reactor

Multiphase flow hydrodynamics in a novel gas–liquid–solid jet-loop reactor (JLR) were experimentally investigated at the macroscales and mesoscales. The chord length distribution was measured by an optical fiber probe and transformed for bubble size distribution through the maximum entropy method. The impacts of key operating conditions (superficial gas and liquid velocity, solid loading) on hydrodynamics at different axial and radial locations were comprehensively investigated. JLR was found to have good solid suspension ability owing to the internal circulation of bubbles and liquid flow. The gas holdup, axial liquid velocity, and bubble velocity increase with gas velocity, while liquid velocity has little influence on them. Compared with the gas–liquid JLRs, solids decrease the gas holdup and liquid circulation, reduces the bubble velocity and delays the flow development due to the enhanced interaction between bubbles and particles (Stokes number >1). This work also provides a benchmark data for computational fluid dynamics (CFD) model validation. © 2019 American Institute of Chemical Engineers AIChE J, 65: e16537, 2019     

新型气固环流反应器内颗粒流动的CFD模拟

采用基于结构的EMMS曳力模型,对一种新型气固环流反应器中的颗粒流动特性进行数值模拟。模拟的固含率与颗粒速率预测值与实验数据具有一致性,验证了模型的适用性。模拟结果表明：导流筒表观气速增加,导流筒中的床层固含率减小,向上的颗粒速率增加;反应器中存在多个颗粒逆流和错流混合区,促进了颗粒沿径向的混合;槽孔处,导流筒中的固含率以及颗粒速率分布更加均匀,而环隙中存在颗粒浓集区;进料区在0≤L≤0.058 m,0< r/R< 0.3的范围内固含率增加并且颗粒存在明显的径向流动。     

装配水平筛板的气升式反应器中气液传质特性的数值模拟

利用Turbulent–Lehr组合模型对装配水平筛板的气升式反应器进行了计算流体力学(CFD)模拟,研究水平筛板对气含率、气泡直径、体积传质系数(kLa)和气液流速的影响。结果表明,筛板对气相的囤积作用和对液相的阻碍作用增加了反应器的整体气含率。筛板对气相的二次均布作用减弱了筛板和液面之间区域的气泡聚并过程,筛板筛孔对气泡的破碎作用产生了大量小于初始直径的气泡,增加了气泡比表面积(a);筛板对液相的阻碍作用提高了筛板附近的气–液相流动速度差,从而提高了该区域的液膜传质系数(kL),强化了反应器内的气液传质效果。     

CFD analysis of two‐phase turbulent flow in internal airlift reactors

The hydrodynamic performance of three internal airlift reactor configurations was studied by the Eulerian–Eulerian k–ε model for a two‐phase turbulent flow. Comparative evaluation of different drag and lift force coefficient models in terms of liquid velocity in the riser and downcomer and gas holdup in the riser was highlighted. Drag correlations as a function of Eötvös number performed better results in comparison to the drag expressions related to Reynolds number. However, the drag correlation as a function of both Reynolds and Eötvös numbers fitted well with experimental results for the riser gas holdup and downcomer liquid velocity in configurations I and II. Positive lift coefficients increase the liquid velocity and decrease the riser gas holdup, while opposite results were obtained for negative values. By studying the effects of bubble size and their shape, the smaller bubbles provide a lower liquid velocity and a gas holdup. The effects of bubble‐induced turbulence and other non‐drag closure models such as turbulent dispersion and added mass forces were analysed. The gas velocity and gas holdup distributions, liquid velocity in the riser and downcomer, vectors of velocity magnitude and streamlines for liquid phase, the dynamics of gas holdup distribution and turbulent viscosity at different superficial gas velocities for different reactor configurations were computed. The effects of various geometrical parameters such as the draft tube clearance and the ratio of the riser to the downcomer cross‐sectional area on liquid velocities in the riser and the downcomer, the gas velocity and the gas holdup were explored. © 2011 Canadian Society for Chemical Engineering     

Longitudinal dispersion of solid particles in fluid beds with horizontal baffles

The axial pressure profiles, allowable gas velocities and temperature distributions are measured for the fluidization of air—FCC cracking catalyst systems in 12- and 19-cm-diam. eight-stage fluid beds equipped with seven horizontal baffles. From these measurements, gas bubble holdup, apparent longitudinal dispersion and intermixing velocity of solid particles through the baffles are studied as functions of baffle design. It is shown that the gas bubble holdup increases, the operational range of gas flow decreases and the flow pattern of solid particles approaches plug flow with decreasing free area of baffles.     

Experimental and modelling study of gas dispersion in a double turbine stirred tank

Gas dispersion in a double turbine stirred tank is experimentally characterised by measuring local gas holdups and local bubble size distributions throughout the tank, for three liquid media: tap water, aqueous sulphate solution and aqueous sulphate solution with PEG. For all these media, bubble coalescence generally prevails over breakage. Where average bubble size decreases, this can be attributed to the difference in slip velocity between different sized bubbles. Most of the coalescence takes place in the turbine discharge stream.A compartment model that takes into account the combined effect of bubble coalescence and breakage is used to simulate gas dispersion. The model predicts spatial distribution of gas holdup and of average bubble size, with average bubble size at the turbines as an input. Reasonable agreement between experiment and simulation is achieved with optimisation of two parameters, one affecting mainly the slip velocity, the other related mainly to the bubble coalescence/breakage balance. Different sets of parameters are required for each of the three liquid systems under study, but are independent of stirring/aeration conditions. The model only fails to simulate the smaller average bubble diameters at the bottom of the tank.