倾斜波纹板上液膜流动的CFD研究

利用VOF法建立了液膜在倾斜波纹板上的气-液两相流CFD模型,并根据液膜流动特点提出了表面张力动量源项和气液界面作用力动量源项.模拟结果与文献实验值吻合较好,表明本文提出的CFD液膜流动模型具有一定的可靠性.通过模拟不同性质的液体在不同表面结构波纹板上的流动过程发现,波纹板表面微观结构以及液体性质尤其是液体的表面张力对连续液膜的形成有重要作用,表明通过改变波纹板面微观结构以及降低液体的表面张力可以促进连续液膜的形成,对提高气液之间的传质效率有重要意义.     

Liquid Flow on Structured Packing: CFD Simulation and Experimental Study

A CFD model of the two‐phase countercurrent flow in the geometry of the plate‐type structured packing Mellapak 250.Y was built, tested and verified. The model was applied to determine the effect of liquid and gas flow rates and physicochemical properties of the flowing liquids on the interfacial area formed on structured packing. The CFD model allowed us to determine the minimum liquid flow rate at which an unbroken liquid film was observed on the packing surface. The simulations confirmed that with an increase of the wetting rate the surface of the packing covered with a liquid film increased until the surface was totally covered up, while further slight changes of an interfacial area were the result of wave formation. The effect of gas load (F factor) on the film surface was in the range of a calculation error. Results of the CFD simulation allow us to predict the stages of film formation during liquid flow, to follow local velocity oscillations, film thickness and velocity profiles of phases.     

Computational Fluid Dynamics Modeling of Gas‐Liquid Two‐Phase Flow around a Spherical Particle

Microscale studies, which can provide basic information for meso‐ and macroscale studies, are essential for the realization of flow characteristics of a packed bed. In the present study, the effects of gas velocity, liquid velocity, liquid‐solid contact angle, and liquid viscosity on the flow behavior were parametrically investigated for gas‐liquid two‐phase flow around a spherical particle, using computational fluid dynamics (CFD) methodology in combination with the volume‐of‐fluid (VOF) model. The VOF model was first validated and proved to be in good agreement with the experimental data. The simulation results show that the film thickness decreases with increasing gas velocity. This trend is more obvious with increasing operating pressure. With increasing liquid velocity, the film thickness tends to be uniform on the particle surface. The flow regime can change from film flow to transition flow to bubble flow with increasing contact angle. In addition, only at relatively high values does the liquid viscosity affect the residence time of the liquid on the particle surface.     

多排平直翅片管换热器表面气液降膜流动特性的三维数值模拟

基于VOF模型建立了考虑重力、表面张力及界面摩擦力源项的多排平直翅片管换热器表面气液两相降膜流动三维瞬态CFD模型。不同气流速度下液膜厚度模拟结果与文献中试验值吻合较好,最大偏差小于5%,表明所建立CFD模型是可靠的。通过研究壁面接触角为30°时不同气液Reynolds数下液膜流动特性,结果表明：翅片管表面满膜流的临界Reynolds数Re_l为239,临界喷淋密度为0.06 kg/(m·s);在239 ≤ Re_l ≤ 995内,其平均液膜厚度较Nusselt理论解高16.8%～35.1%;气液逆流和顺流时气相Reynolds数Re_g应分别小于2190.7和3286.0,其主要原因在于过高的Re_g会导致气液界面摩擦力快速增大,从而引发液膜破裂和液滴脱落等现象恶化设备性能。总之,气液顺流更有利于在较高气相Reynolds数下实现翅片管表面的较薄满膜流动。     

Investigating the liquid film characteristics of gas–liquid swirling flow using ultrasound doppler velocimetry

A novel gas–liquid two‐phase flow metering method was proposed. A spiral vane mounted in the inner pipe was used to transform inlet flow patterns into gas–liquid swirling annular flow. The thickness and velocity profile of liquid film were measured by ultrasound Doppler velocimetry. The liquid flow rates were obtained by integrating of velocity profile during the liquid film zone. Experiments were carried out in an air–water two‐phase flow loop and an ultrasonic transducer was installed under the bottom of the test section with the Doppler angle of 70°. The flow patterns included stratified wavy, annular, and slug flows. Compared with non‐swirling flow, the liquid film thickness at the bottom reduces greatly. The measurement accuracy of liquid flow rate was independent of inlet flow patterns, gas and liquid velocities. © 2016 American Institute of Chemical Engineers AIChE J, 63: 2348–2357, 2017     

Particle‐scale simulation of the flow and heat transfer behaviors in fluidized bed with immersed tube

A kind of new modified computational fluid dynamics‐discrete element method (CFD‐DEM) method was founded by combining CFD based on unstructured mesh and DEM. The turbulent dense gas–solid two phase flow and the heat transfer in the equipment with complex geometry can be simulated by the programs based on the new method when the k‐ε turbulence model and the multiway coupling heat transfer model among particles, walls and gas were employed. The new CFD‐DEM coupling method that combining k‐ε turbulence model and heat transfer model, was employed to simulate the flow and the heat transfer behaviors in the fluidized bed with an immersed tube. The microscale mechanism of heat transfer in the fluidized bed was explored by the simulation results and the critical factors that influence the heat transfer between the tube and the bed were discussed. The profiles of average solids fraction and heat transfer coefficient between gas‐tube and particle‐tube around the tube were obtained and the influences of fluidization parameters such as gas velocity and particle diameter on the transfer coefficient were explored by simulations. The computational results agree well with the experiment, which shows that the new CFD‐DEM method is feasible and accurate for the simulation of complex gas–solid flow with heat transfer. And this will improve the farther simulation study of the gas–solid two phase flow with chemical reactions in the fluidized bed. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

Hydrogenation of 2‐ethylanthraquinone under Taylor flow in single square channel monolith reactors

The hydrogenation of 2‐ethylanthraquinone (EAQ) to 2‐ethylanthrahydroquinone (EAHQ) was carried out under Taylor flow in single square channel monolith reactors. The two opening ends of opaque reaction channel were connected with two circular transparent quartz‐glass capillaries, where Taylor flow hydrodynamics parameters were measured and further used to obtain practical flow state of reactants in square reaction channels. A carefully designed gas‐liquid inlet mixer was used to supply steady gas bubbles and liquid slugs with desired length. The effects of various operating parameters, involving superficial gas velocity, superficial liquid velocity, gas bubble length, liquid slug length, two‐phase velocity and temperature, on EAQ conversion were systematically researched. Based on EAQ conversion, experimental overall volumetric mass transfer coefficients were calculated, and also studied as functions of various parameters as mentioned earlier. The film model, penetration model, and existing semi‐empirical formula were used to predict gas‐solid, gas‐liquid, and liquid‐solid volumetric mass transfer coefficients in Taylor flow, respectively. The predicted overall volumetric mass transfer coefficients agreed well with the experimental ones. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

A CFD‐PBM‐PMLM integrated model for the gas–solid flow fields in fluidized bed polymerization reactors

Although the use of computational fluid dynamics (CFD) model coupled with population balance (CFD‐PBM) is becoming a common approach for simulating gas–solid flows in polydisperse fluidized bed polymerization reactors, a number of issues still remain. One major issue is the absence of modeling the growth of a single polymeric particle. In this work a polymeric multilayer model (PMLM) was applied to describe the growth of a single particle under the intraparticle transfer limitations. The PMLM was solved together with a PBM (i.e. PBM‐PMLM) to predict the dynamic evolution of particle size distribution (PSD). In addition, a CFD model based on the Eulerian‐Eulerian two‐fluid model, coupled with PBM‐PMLM (CFD‐PBM‐PMLM), has been implemented to describe the gas–solid flow field in fluidized bed polymerization reactors. The CFD‐PBM‐PMLM model has been validated by comparing simulation results with some classical experimental data. Five cases including fluid dynamics coupled purely continuous PSD, pure particle growth, pure particle aggregation, pure particle breakage, and flow dynamics coupled with all the above factors were carried out to examine the model. The results showed that the CFD‐PBM‐PMLM model describes well the behavior of the gas–solid flow fields in polydisperse fluidized bed polymerization reactors. The results also showed that the intraparticle mass transfer limitation is an important factor in affecting the reactor flow fields. © 2011 American Institute of Chemical Engineers AIChE J, 58: 1717–1732, 2012     

Regenerierbare feststoff-reaktions-schichten zur reinigung heisser schadstoffhaltiger gase

A method for the purification of a hot off-gas stream by a gas-solid reaction, e.g. neutralisation, is described. A series of metal plates are continuously coated with a solution of a suitable reactant, removed from the solution and dried in the off-gas stream.The reaction between aqueous absorbing solutions or continuously generated solid reactant films and the noxious gases HCl, SO₂ and NO_x has been studied with a small test model in flowing gas streams. In the test model the solid reactant film was generated on metal support plates having 800 cm² area, operating between 25°C and 150°C, the solid reactant films being typically regenerated 10 times per hr.With respect to purifying gas streams the solid films behave similarly, sometimes superior to aqueous scrubbers. Only the reaction of NO_x with alkaline films was not satisfactory. It is likewise that it will take place only at higher temperatures. Removal of 99% of the other noxious gases was demonstrated under the above conditions for 500 1/h of air containing 2500 ppm HCl and for 100 1/h of air containing 2000 ppm SO₄. As the figures indicate the reactive surface area of the solid film is very high. The effective area calculated for a mono-molecular layer of reactants is about 10⁴ times the area of the metal carrier plates.     

竖直多孔平板上液膜流动特性的研究

孔结构被广泛应用于传质塔填料中,对填料上的液膜流动和传质行为影响较大。对竖直光板和多孔板上的液膜流动进行了三维模拟,并通过实验验证了模拟的准确性。通过模拟研究了孔结构对液膜流动特性的影响。结果表明,干燥孔会阻碍液膜的铺展,而润湿孔促进液膜的铺展。与光板相比,多孔板上的液膜具有起伏波,这将影响液膜的厚度分布和速度分布。液膜厚度波动和水平方向的速度波动随着孔径的增加而增加,而竖直流动方向的速度随着孔径的增加而降低。当孔径增加到一定值时,毛细波将出现在孔中的液膜中,这大大增加液膜水平方向上的波动速度,而降低流动方向上的速度。当孔径继续增加到临界值时,液膜将破裂。多孔板上孔内和气侧区域存在涡旋,能够促进内部液体交换和增大气侧扰动,从而增强传质能力。     

Flow characteristics of blade unit of a tridimensional rotational flow sieve tray under concurrent gas–liquid flow

The flow characteristics of the blade unit of a tridimensional rotational flow sieve tray were investigated. First, the flow patterns are defined under different liquid arrangement methods. They are bilateral film flow, continuous perforated flow, and dispersion-mixing flow in overflow distribution, film and jet flow and jet and mixed flow in spray distribution. Second, the time and frequency domain analysis of the differential pressure pulsation signal in the blade unit is carried out, the main frequencies range is 2.0−5.5 Hz. The influence of perforation and mixing intensity on the flow pattern transition is clarified. Third, the rotational flow ratio of the gas–liquid phase is measured, the gas phase rotational flow range is 0.55–0.78, and the liquid phase range is 0.15–0.42. The influence of the operating conditions on the distribution of the rotational and perforated flow is investigated. Finally, a prediction model for the rotational flow ratio is proposed.     

Effect of annealing on morphologies and performances of hydroxypropyl methylcellulose/hydroxypropyl starch blends

The morphologies and performances of the blends of hydroxypropyl methylcellulose (HPMC), which is a thermal gel, and hydroxypropyl starch (HPS), which is a cooling gel, provides a scientific model to understand the relationship between the microstructure and performance of the unique thermal/cooling gel system. The films based on these blends have been developed for various applications, such as food packaging and medicine capsules, mainly for improving processibility. The morphologies and performance of the blends strongly depend on annealing or storage conditions and time since they are neither miscible nor compatible. This work focuses on the effects of storage conditions and time on the morphologies and performance of the blend films. It was found that the morphologies of the blend system depend on temperature, moisture, and annealing time. Correspondently, the gas permeability and mechanical properties of the film were also environmentally and time dependent. When HPS is continuous phase, the annealing results in starch recrystallization, which increases rigidity. When HPMC is the continuous phase, the materials showed less moisture sensitivity. Starch recrystallization and phase separation, which resulted in microcracks on the surface of the films, are the main reasons for the reduction in gas barrier and elongation. This work will study these issues using scanning electron microscope, X-ray diffraction, mechanical testing and establish their relationship.     

A model for collisional exchange in gas/liquid/solid fluidized beds

A new model is presented for numerical simulations of collisional transfer of mass, momentum and energy in gas/liquid/solid fluidized beds. The mathematical formulation uses a collision model similar to that of Bhatnagar, Gross, and Krook (BGK), in a particle distribution function transport equation, in order to approximate the rates at which collisions bring about local equilibration of particle velocities and the masses, compositions, and temperatures of liquid films on bed particles. The model is implemented in the framework of the computational-particle fluid dynamics (CPFD) numerical methodology, in which the particle phase is represented with computational parcels and the continuous phase is calculated on Eulerian finite-difference grid. Computational examples using the Barracuda^® code, a commercial CFD code owned by CPFD Software, LLC, show the ability of the model to calculate spray injection and subsequent liquid spreading in gas/solid flows.     

传热管外液体流动数值模拟研究

采用计算流体力学方法(CFD),对传热管外液体流动过程进行了模拟研究.建立了传热管三维网格模型,考察了不同布液密度下液体在传热管表面铺展形态、流型、液膜厚度以及绕管停留时间的变化规律.该研究结果对于加深理解降膜流动过程,指导蒸发器优化设计和工艺条件优化有一定的借鉴意义.     

Computational Approach to Characterize the Mass Transfer between the Counter‐Current Gas‐Liquid Flow

Structured packed columns are widely used in the chemical industry for distillation and absorption. However, the understanding of the transfer mechanism behind the counter‐current gas‐liquid flow in structured packed columns is still limited. In this work, a three‐dimensional CFD model that considers the local absorption and the local momentum transfer mechanism is developed for a film flow on a small plate with a counter‐current gas flow. The model, based on the Volume of Fluid (VOF) method, is built up on the basis of a pressure drop model and the penetration theory to quantitatively investigate the instantaneous hydrodynamics and mass transfer characteristics of the liquid phase. Simulations and experiments are carried out for a system consisting of propane and toluene. A comparison of the simulation results with the experimental data for the outlet concentrations shows good agreement.     

3‐D numerical simulations on flow and mixing behaviors in gas–liquid–solid microchannels

Three‐dimensional (3‐D) gas‐liquid–solid flow and mixing behaviors in microchannels were simulated by coupled volume of fluid and discrete phase method and simulations were validated against observations. The detachment time and length of gas slug are shortened in liquid–solid flow, compared with those in liquid flow due to higher superficial viscosity of liquid–solid mixture, which will move the bubble formation toward the dripping regime. Solid particles mainly distribute in liquid slug and particle flow shows obvious periodicity. With the increase of contact angle of the inner wall, gas slug (0–50°), stratified (77–120°), and liquid drop (160°) flows are observed. The residence time distributions of solid and liquid phases are similar because particles behave as tracers. The backmixing of solid and liquid phases in liquid drop flow is the weakest among the three flow patterns, and the backmixing of gas phase in slug flow is weaker than that in both stratified and liquid drop flows. The results can provide a theoretical basis for the design of microreactors. © 2013 American Institute of Chemical Engineers AIChE J, 59: 1934–1951, 2013     

Hydrodynamics of inclined packed beds under flow modulation ‐ CFD simulation and experimental validation

A three‐dimensional unsteady‐state Eulerian multi‐fluid CFD model was developed to simulate the hydrodynamic behavior of inclined gas‐liquid cocurrent downflow packed beds under ON‐OFF liquid, ON‐OFF gas, and gas/liquid alternating cyclic operations. Validation of the CFD simulation results was performed with experimental data provided by electrical capacitance tomography imaging. Incorporation in the Eulerian multifluid CFD model of capillary pressure and mechanical dispersion force was essential to accurately capture the transient spatial heterogeneities arising in tilted packed beds under different cyclic modulation strategies. The applied CFD model was able to satisfactorily predict the values of liquid holdup and pressure drop as well as the morphological characteristics of the traveling waves inside the bed for the examined flow modulations. © 2017 American Institute of Chemical Engineers AIChE J, 63: 4161–4176, 2017     

Prediction of the Induced Gas Flow Rate from a Self‐Inducing Impeller with CFD

The complex task of describing computationally two‐phase turbulent flows in aerated stirred‐tank reactors was overcome by proposing that the gas flow rate in the hollow impeller can be estimated from single‐phase flow simulations of the liquid phase in the reactor: the pressure at the impeller surface obtained from liquid phase simulations can be related to the gas induction rate. A commercial lab‐scale reactor with a radial six‐bladed hollow impeller was chosen for the study. To validate the presented methodology, the induced gas flow rate was measured experimentally from the tracking of the position of bubbles in a dynamic sequence of flow images. Notwithstanding the simplifications assumed in the presented CFD methodology, good agreement has been obtained between numerical results and experiments.