三元催化器冷起动排放的数值模拟：数学模型与结果分析

This paper integrated a two-dimensional axisymmetrical transient model applicable to cold-start emission applications. The model can be used to simulate and explain effects of the flow and temperature distribution on performance of a converter. The evolutions of distribution of the temperature and concentration in the monolith during the cold-start period and the effects of flow distribution in the monolith on the cold-start performance are simulated in terms of the integrated model. The investigation indicates that the axial and radial gradients of temperature of the solid become steeper as the inlet gas temperature ramp increases; this furthermore results in the movement of reaction region in the monolith, and the flow distribution in the monolith affects the radial distribution of temperature of the solid;the radial gradients of temperature of the solid become greater as the flow uniformity index decreases, whereas the light-off time doesn‘t always increase as the flow uniformity index decreases. The analyses on the distribution of temperature and concentration in the monolith show that the catalytic reaction zone concentrates in central area near the front face. The predicted curves of the velocity distribution have a good agreement with the experimental data.     

Runaway Criterion in Fixed Bed Catalytic Reactors with Radial Temperature Profile

The discrepancy between pseudo-homogeneous one-dimensional model and peeudo-homogeneous two-dimensional model is studied. It is found that there are great differences between two models. This paper compares the maximum and minimum values of the radial temperature in the hot spot in came that a single exothermic reaction is carried out, a correlation is obtlioed with peeudo-homogeneous one-dimensional model to describe the entire reactor behavier. A new runaway criterion, based on the occurrence of inflection in the hot spot locus, is developed for the case of pseudo-homogeneous two-dimensional model. This criterion predicts the maximum allowable temperature for safe operation and the regions of runaway, respectively. The calculated results show that, compared with the results based on pseudo-homogeneous one-dimensional model, runaway will easily occur when the radial temperature gradient has to be considered.     

轴向流固定床内流场的数值模拟与实验验证

The computational fluid dynamics model with porosity and drag coefficient was used to describe fluid flow in an axial flow fixed bed according to the characteristics of fluid flow in the fixed-bed of the reactor. The commercial computational fluid dynamics (CFD) code CFX was used to simulate the flow field in the axial flow fixed bed. The simulation predictions are in good agreement with experimental results of a large cold model. The influence of gas distributor on the flow field in the axial flow fixed bed was studied. A suitable gas distributor was used to attain less than 0.06 kPa radial pressure difference and less than 5.2% radial velocity difference in fixed bed.     

HYDRODYNAMIC MODEL FOR THE IMPELLER REGION IN AN AERATED AND STIRRED TANK

The radial and axial distribution of mean 1iquid velocity were measured by a.hot-filmanemometer at the impeller region in an aerated and stirred tank 0.287m in diameter.The tangentialjet model for impeller discharge flow used for single phase flow was modified to conform with thecharacteristics of gas-liquid flow.The radial and axial velocity profiles at the impeller region in thegas-liquid stirred tank were calculated by the model The results predicted by the model were in goodagreement with those obtained in experiment.     

喷嘴进料对催化裂化提升管流动行为的影响

This paper studies the influence of feed injection on the hydrodynamic behavior of fluid catalytic cracking riser reactors. Experiments were conducted in a cold model of 186mm ID with two oppositely inclined secondary air feed nozzles. The flow structure was determined by means of the axial static pressure measurements and local radial optic fiber probe measurements on different levels with emphasis on the sections downstream of the secondary injection. The measurements reveal that the secondary injection plays a crucial role on riser hydrodynamics. Just above the secondary injection, the flow and mixing are strongly affected, while below the secondary injection the effect is weak. The radial profile just downstream of secondary injection demonstrates wavy features. The effective region of secondary injection could be estimated by the axial pressure gradient profiles and/or the radial orofiles of local solids velocity and density.     

BUBBLE SIZES AND HOLDUP STRUCTURE IN BUBBLE COLUMNS DETERMINED BY DYNAMIC GAS DISENGAGEMENT

A new model for determining bubble size distributions in bubble columns by the dynamic gasdisengagement(DGD)technique is developed.It is based on an idea of non-uniform steady statedirstribution of bubble dispersion.Interpreting the axial non-uniformity,this model gives axial gasholdup distributions.If assuming an axially homogeneous dispersion,a radial gas holdup distributioncan be obtained.The Sauter mean diameters or specific interfacial areas for several systems areestimated by the technique.The results for an air-water system agree with those measured by afive-point conductivity probe technique.The obtained axial gas holdup distributions agree well withreported measurements and the radial gas holdup distributions are also reasonable.     

Phase holdup and liquid dispersion in gas-liquid-solid circulating fluidized beds

Axial and radial profiles of gas and solids holdups have been studied in agas-liquid-solid circulating fluidized bed at 140mm i.d..Experimental results indicate that the axialand radial profiles of gas and solids holdups are more uniform than those in a conventionalfluidized bed.Axial and radial liquid dispersion coefficients in the gas-liquid-solid circulating fluidizedbed are investigated for the first time.It is found that axial and radial liquid dispersioncoefficients increases with increaes in gas velocity and solids holdup.The liquid velocity has littleinfluence on the axial liquid dispersion coefficient,but would adversely affect the redial liquiddispersion coefficient.It can be concluded that the gas-liquid-solid circulating fluidized bed hasadvantages such as better interphase contact and lower liquid dispersion along the axial directionover the expanded bed.     

THE RADIAL DISTRIBUTION OF HEAT TRANSFER COEFFICIENT IN FAST FLUIDIZED BED

This paper presents an experimental study on the axial and radial distributions of heat transfer coefficients in a fastfluidized bed operated at ambient temperature.The following formula is recommended to correlate the local heat transfer coefficients in fast fluidized beds:where a,n_1 and n_2 are functions of radial positions.     

THE PARTICLE VELOCITY DISTRIBUTION IN A RISER OF DILUTE CIRCULATING FLUIDIZED BED

This paper presents an experimental study on the radial and axial profiles of local particle velocity of flu-idized catalyst(FCC)particles in a circulating fluidized bed of 8000mm in height and 186mm in inside diameterwith a Fiber-Optic Probe Laser Doppler Velocimetry.It has been found that the radial profiles of local particlevelocity are parabolic,even in the centre region and steep in the region near the bed wall,and that the length ofparticle acceleration along the bed height is longer than expected.The following formula has been recommended to correlate the 1 ocal particle velocity with dilute gas-solidflow in the bed:where a,b,c and e are functions of radial positions.     

A MODEL FOR BED-TO-SURFACE HEAT TRANSFER IN CIRCULATING FLUIDIZED BEDS

Experiments on bed-to-surface heat transfer in the axial and radial directions were carried out in a cold flowmodel of the circulating fluidized bed of 244mm i.d.and with FCC catalysts(d_p=83.4μm)as the fluidized particles.A simplified gas-solid flow model comprising of the cluster-particulate phases upflowing cocurrentlywas proposed.By taking into account the alternative contacts of cluster and particulate phases with the immersedsurface,a composite heat transfer model was developed.The model gave good predictions as compared withthe experimental results of this work and those of others.     

Determination of parameters in a two-dimensional model for fixed-bed adsorbers

A new method is proposed for the determination of the parameters in a two-dimensionalmodel which characterizes the properties of axial and radial mixing and mass transport in afixed-bed adsorber.Parameter estimation for the model is carried out with methane-air-5A molecularsieve in a bed under the condition of step injection of tracer from a point on the main axis of thebed by the curve fitting method in the time domain.     

A Sedimentation-Dispersion Model for both Non-attached and Attached Particles in Three-Phase Batchwise Fluidized Beds

The axial concentration distribution of both particles with better wetting (forming non-attached system) and poorer wetting (forming attached system) was investigated in a vertical gas-liquid-solid fluidized bed of 4.2 cm in diameter and 130 cm in height with the solids holdup less than 0.05. The one-dimensional sedimentation-dispersion model could be used satisfactorily to describe the axial distribution of solids holdup by modifying only a model parameter, i.e. by means of the terminal settling velocity minus a certain value, which is a function of gas velocity and considers the effect of an additional drag force resulted from attached rising bubbles. The axial profiles of solid concentration predicted are in good agreement with experimental results. This model also explains reasonably the different axial distributions of solid concentration, i.e. the solids holdup decreases as the axial height increases in non-attached system, but increases with the axial height in attached system at a given gas veloc     

三相间歇流化床中非粘附与粘附体系的沉降-分散模型

The axial concentration distribution of both particles with better wetting (forming non-attached systems) and poorer wetting (forming attached system) was investigated in a vertical gas-liquid-solid fluidized bed of 4.2cm in diameter and 130cm in height with the solids holdup less than 0.05. The one-dimensional sedimentation-dispersion model could be used satisfactorily to describe the axial distribution of solids holdup by modifying only a model parameter,i.e.by means of the terminal settling velocity minus a certain value,which is a function of gas velocity and considers the effect of an additional drag force resulted from attached rising bubbles.The axial profiles of solid concentration predicted are in good agreement with experimental results.This model also explains reasonably the different axial distribution of solid concentration,i.e.the solids holdup decreases as the axial height increases in non-attached system,but increases with the axial height in attached system at a given gas velocity.     

水力旋流分离器内流动和油水分离的数值模拟

The fluid flow and oil-water separation were simulated using a Reynolds stress transport equation model of turbulence in water flow and a stochastic model of oil droplet motion,Simulation results give the axial and tangential velocity components,the pressure and turbulence intensity distribution and droplet trajectories for a hydrocyclone of F type and a hydrocyclone proposed by the present authors.The flow filed predictions are in qualitative agreement with the LDV measurements.The results show that the proposed hydrocyclone has better performance than the hydrocyclone of F type due to creating stronger centrifugal force and lower axial velocity.     

A Modified Peng-Robinson Equation State for Prediction of Gas Adsorption Isotherm

This research proposes a modified two-dimensional Peng-Robinson equation model to predict adsorp-tion isotherm in adsorbate-adsorbent systems. The parameters of the proposed model are calculated by using the op-timization of experimental data for the different single gas adsorption systems at various temperatures. The experi-mental adsorption equilibrium data of adsorbate-adsorbent systems was compared with the calculated results in our proposed model and the two-dimensional Hill-deBoer equation model. The proposed model as indicated in the re-sults shows a better prediction of the experimental results compared with two others.     

合成气制二甲醚三相淤浆床反应器的数学模型研究

A mathematical model for a bubble column slurry reactor is presented for dimethyl ether synthesis from syngas. Methanol synthesis from carbon monoxide and carbon dioxide by hydrogenation and the methanol dehydration are considered as independent reactions, in which methanol, dimethyl ether and carbon dioxide are the key components. In this model, the gas phase is considered to be in plug flow and the liquid phase to be in partly back mixing with axial distribution of solid catalyst. The simulation results show that the axial dispersion of solid catalysts, the operational height of the slurry phase in the bubble column slurry reactor, and the reaction results are influenced by the reaction temperature and pressure, which are the basic data for the scale-up of reactor.     

涡流桨搅拌槽内流动非稳定性的大涡模拟

The aim of this work is to investigate the flow instabilities in a baffled, stirred tank generated by a single Rushton turbine by means of large eddy simulation （LES）. The sliding mesh method was used for the coupling between the rotating and the stationary frame of references. The calculations were carried out on the ＂Shengcao-21C＂ supercomputer using a computational fluid dynamics （CFD） code CFX5. The flow fields predicted by the LES simulation and the simulation using standard κ-ε model were compared to the results from particle image velocimetry （PIV） measurements. It is shown that the CFD simulations using the LES approach and the standard κ-ε model agree well with the PIV measurements. Fluctuations of the radial and axial velocity are predicted at different frequencies by the LES simulation. Velocity fluctuations of high frequencies are seen in the impeller region, while low frequencies velocity fluctuations are observed in the bulk flow. A low frequency velocity fluctuation with a nondimensional frequency of 0.027Hz is predicted by the LES simulation, which agrees with experimental investigations in the literature. Flow circulation patterns predicted by the LES simulation are asymmetric, stochastic and complex, spanning a large portion of the tanks and varying with time, while circulation patterns calculated by the simulation using the standard κ-ε model are symmetric. The results of the present work give better understanding to the flow instabilities in the mechanically agitated tank. However, further analysis of the LES calculated velocity series by means of fast Fourier transform （FFT） and/or spectra analysis are recommended in future work in order to gain more knowledge of the complicated flow phenomena.     

GAS-SOLIDS FLOW PATTERNS IN A CONCURRENT DOWNFLOW FAST FLUIDIZED BED(CDFFB)

Radial profiles of solid concentration and velocity for concurrent downward gas-solid suspension in a140mm inside diameter fast fluidized bed were investigated.The influence of gas velocity,solid circulating rateand axial position on radial profiles of solid concentration and particle velocity has been examined.It hasbeen found that an annular region of high solid concentration exists at r/R=0.94.At both the center and wallregion,the solid concentration and the particle velocities are relatively low.The shape of radial solid con-centration profile curves is mainly dependent on the cross-section averaged voidage,and the shape of radialparticle velocity profile is mainly affected by the gas velocity and cross-section averaged voidage.Based on the radial profiles of solid concentration and particle velocity,the solid mass flux profile and thenonuniformity of solids flow are discussed in this paper.It is shown that solids flow in CDFFB is much moreuniform than that in UFFB.     

Turbulent Liquid Flow in a Gas-Liquid Bubble Column

The radial distribution of the axial component of time-average liquid velocity and tur-bulent intensity in a gas-liquld bubble column was investigated experimentally using Laser Doppler Anemometer (LDA). The effects of operating parameters on liquid turbulent intensity are studied and an empirlcal relationship between turbulent intensity and viscosity was established, Such a relationship can he used convenieutly in the calculation of liquid velocity profilea in bubble columns.     

分级进风燃烧室内旋流反应流的湍流特性

This paper presents an experimental investigation of the turbulent reacting flow in a swirl combustor with staged air injection. The air injected into the combustor is composed of the primary swirling jet and the secondary non-swirling jet. A three dimension-laser particle dynamic analyzer （PDA） was employed to measure the instantaneous gas velocity. The probability density functions （PDF） for the instantaneous gas axial and tangential velocities at each measuring location, as well as the radial profiles of the root mean square of fluctuating gas axial and tangential velocities and the second-order moment for the fluctuating gas axial and tangential velocities are obtained. The measured results delineate the turbulence properties of the swirling reacting flow under the conditions of staged combustion.