用扩散流动模型分析悬浮床内的气固两相向上流动

A mathematical model of two-dimensional turbulent gas-particle two-phase flow based on the modified diffusion flux model (DFM) and a numerical simulation method to analyze the gas-particle flow structures are developed. The modified diffusion flux model, in which the acceleration due to various forces is taken into account for the calculation of the diffusion velocity of particles, is applicable to the analysis of multi-dimensional gas-particle two-phase turbulent flow. In order to verify its accuracy and efficiency, the numerical simulation by DFM is compared with experimental studies and the prediction by κ-ε-κp two-fluid model, which shows a reasonable agreement. It is confirmed that the modified diffusion flux model is suitable for simulating the multi-dimensional gas-particle two-phase flow.     

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

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.     

An integrated model for predicting the flame propagation in crimped ribbon flame arresters

Crimped ribbon flame arresters are important safety devices in the chemical industry,especially for the dangerous situations.Although proper design of arresters by the numerical simulation method is promising,its reliability and accuracy are dependent upon the mathematical model.In this work,an integrated mathematical model for the microchannel in the crimped ribbon flame arresters was set up;the fluid flow behavior and the sensitivities of four chemical kinetics mechanisms of propane-air on the accuracy were analysed.It is shown that turbulence is predominant in the microchannel of the crimped ribbon flame arresters under the deflagration and detonation conditions,and a new quenching criterion for the numerical simulation is proposed.The kinetics mechanism of Mansouri et al.among the four ones is the most accurate due to the best agreement of the predicted outlet temperature at the experimental flameproof velocity with the autoignition temperature of propane-air.The species mass fraction profiles and the temperature distribution,which are too difficult to measure due to the tiny dimension of the microchannel in experiments,are captured.The fundamental insights into chemical reactions and heat loss are well portrayed.It can be concluded that the integrated mathematical model established in this work can be used as a reliable tool for modeling,selecting and designing such type of crimped ribbon flame arresters with the propane-air medium in the future.     

数值模拟黏弹流体通过椭圆环模具的挤出胀大

The numerical simulation of extrudate swell is significant in extrusion processing.Precise prediction of extrudate swell is propitious to the control of melt flow and the quality of final products.A mathematical model of three-dimensional(3D)viscoelastic flow through elliptical ring die for polymer extrusion was investigated.The penalty function formulation of viscoelastic incompressible fluid was introduced to the finite element model to analyze 3D extrusion problem.The discrete elastic viscous split stress(DEVSS)and streamline-upwind PetrovGalerkin(SUPG)technology were used to obtain stable simulation results.Free surface was updated by updating the streamlines which needs less memory space.According to numerical simulation results,the effect of zero-shear viscosity and elongation parameter on extrudate swell was slight,but with the increase of volumetric flow rate and relax time the extrudate swell ratio increased markedly.Finally,the numerical simulation of extrudate swell flow for low-density polyethylene(LDPE)melts was investigated and the results agreed well with others’work.These conclusions provided quantitative basis for the forecasting extrudate swell ratio and the controlling of extrusion productivity shape.     

Eulerian–Lagrangian simulation of bubble coalescence in bubbly flow using the spring-dashpot model

The Eulerian–Lagrangian simulation of bubbly flow has the advantage of tracking the motion of bubbles in continuous fluid, and hence the position and velocity of each bubble could be accurately acquired. Previous simulation usually used the hard-sphere model for bubble–bubble interactions, assuming that bubbles are rigid spheres and the collisions between bubbles are instantaneous. The bubble contact time during collision processes is not directly taken into account in the collision model. However, the contact time is physically a prerequisite for bubbles to coalesce, and should be long enough for liquid film drainage. In this work we applied the spring-dashpot model to model the bubble collisions and the bubble contact time, and then integrated the spring-dashpot model with the film drainage model for coalescence and a bubble breakage model. The bubble contact time is therefore accurately recorded during the collisions. We investigated the performance of the spring-dashpot model and the effect of the normal stiffness coefficient on bubble coalescence in the simulation.The results indicate that the spring-dashpot model together with the bubble coalescence and breakage model could reasonably reproduce the two-phase flow field, bubble coalescence and bubble size distribution. The influence of normal stiffness coefficient on simulation is also discussed.     

搅拌槽内微观混合研究的回顾与展望(英文)

Mixing problems are most likely encountered and sometimes can be severe in scaling-up projects.Micro-mixing is an important aspect for fast or quasi-instantaneous reactions.Poor micro-mixing might produce more undesired by-products,leading to higher purification costs.This paper gives an extensive review and analysis of micro-mixing studies in single-and multiphase stirred tanks.The relevant experiment techniques,micro-mixing models and numerical approaches are critically reviewed and analyzed with remarks and perspectives.The reported studies on two-phase micro-mixing experiments and on the impact of the presence of the dispersed phases on turbu-lence have been limited to a narrow range of conditions.More importantly,disparities widely exist among different reports.Both Lagrangian and Eulerian models are based on oversimplified assumptions,which may lead to uncer-tainties or even unrealistic results.A heuristic model,which is from the perspective of CFD(computational fluid dynamics) and can cover the whole spectrum of scales and also focus on every sub-process,is desired in the future.     

基于精细自由基反应机理的丙烷裂解炉炉管CFD模拟(英文)

In the radiant section of cracking furnace, the thermal cracking process is highly coupled with turbulent flow, heat transfer and mass transfer. In this paper, a three-dimensional simulation of propane pyrolysis reactor tube is performed based on a detailed kinetic radical cracking scheme, combined with a comprehensive rigorous computational fluid dynamics（CFD） model. The eddy-dissipation-concept（EDC） model is introduced to deal with turbulence-chemistry interaction of cracking gas, especially for the multi-step radical kinetics. Considering the high aspect ratio and severe gradient phenomenon, numerical strategies such as grid resolution and refinement, stepping method and relaxation technique at different levels are employed to accelerate convergence. Large scale of radial nonuniformity in the vicinity of the tube wall is investigated. Spatial distributions of each radical reaction rate are first studied, and made it possible to identify the dominant elementary reactions. Additionally, a series of operating conditions including the feedstock feed rate, wall temperature profile and heat flux profile towards the reactor tubes are investigated. The obtained results can be used as scientific guide for further technical retrofit and operation optimization aiming at high conversion and selectivity of pyrolysis process.     

旋流板内两相流场的CFD模拟与分析

The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The simulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, mad a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall;the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.     

Heat transfer and friction factor of Therminol liquid phase heat transfer fluid in a ribbed tube

Experiments and simulations on flow and heat transfer behavior of Therminol-55 liquid phase heat transfer fluid have been conducted in a ribbed tube with the outer diameter and inner diameter 25.0 and 20.0 mm,pitch and rib height of 4.5 and 1.0 mm.respectively.Experimental results show that the heat transfer and thermal performance of Therminol-55 liquid phase heat transfer fluid in the ribbed tube are considerably improved compared to those of the smooth tube.The Nusselt number increase with the increase of Reynolds number.The increase in heat transfer rate of the ribbed tube has a mean value of 2.24 times.Also,the pressure drop results reveal that the average friction factor of the ribbed tube is in a range of 2.4 and 2.8 times over the smooth tube.Numerical simulations of three-dimensional flow behavior of Therminol-55 liquid phase heat transfer fluid are carried out using three different turbulence models in the ribbed tube.The numerical results show that the heat transfer of ribbed tube is improved because vortices are generated behind ribs,which produce some disruptions to fluid flow and enhance heat transfer compared with smooth tube.The numerical results prove that the ribbed tube can improve heat transfer and fluid flow performances of Therminol liquid phase heat transfer fluid.     

Numerical flow simulation of viscoplastic fluids in annuli

This paper describes numerical solutions for the laminar flow of non-Newtonian fluids in vertical annuli using the Herschel-Bulkley model to describe the rheological behaviour of such materials. Numerical solutions have been obtained when there is both axial and tangential flows in either a concentric or eccentric annulus. The tangential flow arises from the rotation of the inner cylinder of the annulus and the axial flow from a constant axial pressure gradient. The flow is analysed by solving the momentum and continuity equation numerically using the finite element method. The dimensionless velocity, deformation and stress profiles with other quantities such as the apparent viscosity and pressure distribution have been calculated for various eccentricities, radius ratios, fluid properties and flow parameters; the results give insights into the flow behaviour in the annuli. It is shown that the inclusion of rotational effects, for a fixed pressure gradient, is likely to increase the axial volumetric flowrate over non-rotating situations in concentric geometries. New results reveal that, in eccentric annuli, the situation is reversed and the flowrate gradually decreases as the rotation rate is increased.     

Viscoplastic fluid flow in irregular eccentric annuli due to axial motion of the inner pipe

Fully developed axial laminar flow of viscoplastic Herschel-Bulkley fluids in eccentric annuli between two pipes has been investigated numerically. The pipes are closed at one end and flow is due to the axial motion of the inner pipe. The annuli may be filly open or partially blocked. General non-orthogonal, boundary-fitted curvilinear coordinates have been used to accurately model the irregular annular geometry due to the presence of a flow blockage. A computer code has been developed using a second-order finite-difference scheme. An exponential model for the shear stress, valid for both yielded and unyielded regions of the flow, is used in the computation. The effects of generalized Bingham number, flow behavior index, eccentricity, and blockage height on the pressure gradient or the surge pressure have been studied and the results are presented in dimensionless form. The pressure gradient is found to decrease with increasing eccentricity. For a partially blocked eccentric annulus the pressure gradient is found to decrease with an increase in the blockage height.     

LOSS COEFFICIENTS IN ENTRANCE REGION FLOWS OF NEWTONIAN AND NON-NEWTONIAN FLUIDS IN ECCENTRIC ANNULI

Loss coefficients for laminar flow of an incompressible Newtonian and non-Newtonian fluids flowing in an eccentric annulus are calculated by modelling an eccentric annulus as a slit of variable height. While the loss coefficients for a power-law fluid are smaller than those for a Newtonian one, the reverse is true for a Bingham fluid. In the limiting case of a concentric annulus, the results are found to be in good agreement with those reported in the literature.     

Flow simulation of herschel-bulkley fluids through extrusion dies

The flow of viscoplastic materials through extrusion dies has been studied numerically using the finite element method. Rheological data for viscoplastic doughs have been fitted by the Herschel-Bulkley model, which incorporates a yield stress into the power-law model. Non-isothermal simulations show the extent and shape of yielded/unyielded regions and the development of temperature field assuming different modes of heat transfer at the boundaries. The results reveal that viscous dissipation causes appreciable temperature rises in the extrudate in agreement with measured values at the extruded material surface. The extrudate swell results show a maximum for a certain range of apparent shear rates also observed experimentally. However, the inelastic simulations based on the Herschel-Bulkley model always under-predict the experimental swelling values. A heuristic approach is also used to determine the level of elasticity required to produce the experimental values.     

Toward Understanding of Two‐Phase Eccentric Helical Reactor Performance

Mass transfer investigations in a two‐phase gas‐liquid Couette‐Taylor flow (CTF) reactor and a numerical flow simulation are reported. The CTF reactor is characterized by high values of the mass transfer parameters. Previous mass transfer investigations have yielded high values of the volumetric mass transfer coefficients (of the order of 10^–1 s^–1) and the specific interfacial area, compared to those obtained in a stirred tank (10³ m² m^–3). In order to intensify mass transfer in the CTF reactor, an eccentric rotor (rotating inner cylinder) was used. In the eccentric annulus with rotating inner cylinder, due to frequent variation of the hydrodynamic flow field parameters, nonlinear hydrodynamic conditions occurred. These conditions can influence the rate of mass transfer. The experimental results of benzaldehyde oxidation in an eccentric CTF reactor confirmed an increase in mass transfer, as against a concentric CTF reactor. Numerical simulation of the Couette‐Taylor (helical) flow was performed in a concentric and in an eccentric annulus. Calculation of parameters such as velocity, static pressure, kinetic energy and energy dissipation rate revealed a significant effect of gap eccentricity on the flow behavior.     

A COMPUTATIONAL ALGORITHM FOR CALCULATING DISPLACEMENT EFFICIENCY IN HORIZONTAL ANNULI - ECCENTRIC CASE

This paper expands the current model, for laminar displacement of one Power-law fluid by another in concentric horizontal annuli (Nguyen and Rahman, 2000), to take into account the effect of eccentricity. A calculation procedure is also presented which, can be used for computer modelling of laminar displacement process in eccentric horizontal annuli. Simulated results clearly demonstrate the detrimental effect of eccentricity on the displacement process, and suggest that the flow velocity in the narrow part of an eccentric annulus should be one of the major design considerations for an effective horizontal annular displacement.     

A computational algorithm for calculating displacement efficiency in horizontal annuli - eccentric case

This paper expands the current model, for laminar displacement of one Power-law fluid by another in concentric horizontal annuli (Nguyen and Rahman, 2000), to take into account the effect of eccentricity. A calculation procedure is also presented which, can be used for computer modelling of laminar displacement process in eccentric horizontal annuli. Simulated results clearly demonstrate the detrimental effect of eccentricity on the displacement process, and suggest that the flow velocity in the narrow part of an eccentric annulus should be one of the major design considerations for an effective horizontal annular displacement.     

A simplified couette flow solution of non-newtonian power-law fluids in eccentric annuli

Flow of non-Newtonian fluids in both the concentric and eccentric annuli is of great importance in extruders for molten plastics and wellbore fluid circulation for the removal of drilling cuttings. The steady laminar couette flow of non-Newtonian power-law fluids in eccentric annulus is employed in this study to analyze the problems of surge or swab pressures encountered when running or pulling tubular goods (pipes) in a liquid filled borehole. This is similar to the annular space created by two long co-axial cylinders with the inner cylinder in motion at a steady velocity, and a stationary outer cylinder. The solutions of the equations of motion are presented in both dimensionless form and as a family of curves for different pipe/borehole eccentricity ratios and power-law fluid index values for a more general application. The expected error in surge computation for concentric annulus as a result of eccentricity is evaluated.     

Drag of a cylinder moving near a wall in a yield stress fluid

The drag of a cylindrical obstacle moving at a constant velocity in a yield stress fluid close to a wall is studied experimentally and numerically. The wall influence has been explored for gap values between the cylinder of diameter D and the wall ranging from 0.01D to 100D, which corresponds, respectively, to hydrodynamic lubrication and to unconfined domain conditions. A model yield stress fluid (Carbopol gel) is used in the experiments. The viscous and plastic drag coefficients have been calculated and measured as depending on the Oldroyd number, in conditions where the yield stress effects are more important than those of viscosity and the inertia negligible. We have performed experimental and numerical validations in the Newtonian case and provided more specifically comparisons of our measured data on yield stress materials with those resulting from viscoplastic flow simulations. © 2018 American Institute of Chemical Engineers AIChE J, 64: 4118–4130, 2018     

幂律流体偏心环空流动的数值计算

在幂律流体偏心环空层流流动的理论分析基础上,运用控制体积法对幂律流体偏心环空轴向层流的基本方程进行了离散化处理,并利用ADI方法进行了求解,得到了流场中的速度分布、窄宽边速度比等流场参数,结果表明,当偏心度较大时,忽略周向应力的影响会使计算结果产生较大的偏差;偏心度对流场的影响最为明显,偏心度增大时,流量增大,窄宽边的速度比减小。