气固两相流介尺度LBM-DEM模型

LBM-DEM耦合方法通常是指一种颗粒流体系统直接数值模拟算法,即是一种不引入经验曳力模型的计算方法,颗粒尺寸通常比计算网格的长度大一个量级,颗粒的受力通过表面的粘性力与压力积分获得,其优点是能描述每个颗粒周围的详细流场,产生详细的颗粒-流体相互作用的动力学信息,可以探索颗粒流体界面的流动、传递和反应的详细信息及两相相互作用的本构关系,但其缺点是计算量巨大,无法应用于真实流化床过程模拟。本文针对气固流化床中的流体以及固体颗粒间的多相流体力学行为,建立了一种稠密气固两相流的介尺度LBMDEM模型,即LBM-DEM耦合的离散颗粒模型,实现在颗粒尺度上流化床的快速离散模拟。该耦合模型采用格子玻尔兹曼方法(LBM)描述气相的流动和传递行为,离散单元法(DEM)用于描述颗粒相的运动,并利用能量最小多尺度(EMMS)曳力解决气固耦合不成熟问题,以提高其模拟精度。通过经典快速流态化的模拟,验证了介尺度LBM-DEM耦合模型的有效性。模拟结果表明介尺度LBM-DEM模型是一种探索实验室规模气固系统的有力手段。     

EFFECT OF NOZZLE FAN ANGLE ON SPRAYS IN GAS-SOLID RISER FLOW

A three-dimensional simulation study is performed for investigating the hydrodynamic behaviors of a cross-flow liquid nitrogen spray injected into an air-fluidized catalytic cracking (FCC) riser of rectangular cross-section. Rectangular nozzles with a fixed aspect ratio but different fan angles are used for the spray feeding. While our numerical simulation reveals a generic three-phase flow structure with strong three-phase interactions under rapid vaporization of sprays, this paper tends to focus on the study of the effect of nozzle fan angle on the spray coverage as well as vapor flux distribution by spray vaporization inside the riser flow. The gas-solid (air-FCC) flow is simulated using the multi-fluid method while the evaporating sprays (liquid nitrogen) are calculated using the Lagrangian trajectory method, with a strong two-way coupling between the Eulerian gas-solid flow and the Lagrangian trajectories of spray. Our simulation shows that the spray coverage is basically dominated by the spray fan angle. The spray fan angle has a very minor effect on spray penetration. The spray vaporization flux per unit area of spray coverage is highly non-linearly distributed along the spray penetration. The convection of gas-solid flow in a riser leads to a significant downward deviation of vapor generated by droplet vaporization, causing a strong recirculating wake region in the immediate downstream area of the spray.     

Accelerating fully resolved simulation of particle-laden flows on heterogeneous computer architectures

An efficient computing framework, namely PFlows, for fully resolved-direct numerical simulations of particle-laden flows was accelerated on NVIDIA General Processing Units (GPUs) and GPU-like accelerator (DCU) cards. The framework is featured as coupling the lattice Boltzmann method for fluid flow with the immersed boundary method for fluid-particle interaction, and the discrete element method for particle collision, using two fixed Eulerian meshes and one moved Lagrangian point mesh, respectively. All the parts are accelerated by a fine-grained parallelism technique using CUDA on GPUs, and further using HIP on DCU cards, i.e., the calculation on each fluid grid, each immersed boundary point, each particle motion, and each pair-particle collision is responsible by one computer thread, respectively. Coalesced memory accesses to LBM distribution functions with the data layout of Structure of Arrays are used to maximize utilization of hardware bandwidth. Parallel reduction with shared memory for data of immersed boundary points is adopted for the sake of reducing access to global memory when integrate particle hydrodynamic force. MPI computing is further used for computing on heterogeneous architectures with multiple CPUs-GPUs/DCUs. The communications between adjacent processors are hidden by overlapping with calculations. Two benchmark cases were conducted for code validation, including a pure fluid flow and a particle-laden flow. The performances on a single accelerator show that a GPU V100 can achieve 7.1–11.1 times speed up, while a single DCU can achieve 5.6–8.8 times speed up compared to a single Xeon CPU chip (32 cores). The performances on multi-accelerators show that parallel efficiency is 0.5–0.8 for weak scaling and 0.68–0.9 for strong scaling on up to 64 DCU cards even for the dense flow (φ = 20%). The peak performance reaches 179 giga lattice updates per second (GLUPS) on 256 DCU cards by using 1 billion grids and 1 million particles. At last, a large-scale simulation of a gas-solid flow with 1.6 billion grids and 1.6 million particles was conducted using only 32 DCU cards. This simulation shows that the present framework is prospective for simulations of large-scale particle-laden flows in the upcoming exascale computing era.     

THERMAL BOUNDARY LAYER IN CFB BOILER RISER

Measurement of temperature profiles of gas-solid two-phase flow at different heights in commercial-scale circulating fluidized bed (CFB) boilers was carried out. Experimental results showed that the thickness of thermal boundary layer was generally independent of the distance from the air distributor, except when close to the riser outlet. Through analysis of flow and combustion characteristics in the riser, it was found that the main reasons for the phenomena were: 1) the hydrodynamic boundary layer was thinner than the thermal layer and hardly changed along the CFB boiler height, and 2) both radial and axial mass and heat exchanges were strong in the CFB boiler. Numerical simulation of gas flow in the outlet zone confirmed that the distribution of the thermal boundary layer was dominated by the flow field characteristics.     

On the macroscopic dynamics induced by a model wave-particle collision operator

The macroscopic dynamics of a kinetic equation involving a model wave-particle collision operator of plasma physics is investigated. The Chapman-Enskog asymptotics is first considered in the framework of a hydrodynamic scaling. The obtained macroscopic model still involves a kinetic variable, the particle energy in the rest frame of the fluid, but shares similarities with the compressible Navier-Stokes equation of gas dynamics. Then a diffusive scaling is examined under the hypothesis of small perturbations of a global equilibrium. In this case, the macroscopic model couples the usual incompressible Navier-Stokes with a diffusion equation for the energy distribution function of the particles, constrained by an extended version of the Boussinesq relation. In both cases, the effect of a Lorentz force term is developed, in the perspective of plasma physical modelling. Received June 16, 1997     

二维喷动床CFD-DEM模拟中曳力模型的影响

CFD-DEM已经广泛应用到喷动床的研究中,其模拟的准确性与用于处理颗粒-流体相互作用的曳力模型密切相关。为了探究不同曳力模型对喷动床CFD-DEM模拟结果的影响,基于非结构化网格的喷动床仿真,使用7个曳力模型分别对锥底喷动床内气固两相运动进行了数值模拟。综合床层压降、喷动高度和颗粒速度特性三个方面,Wen-Yu模型和Gibilaro模型预测的气固两相运动最剧烈,其次是Di Felice模型、Syamlal-O’Brien模型、Gidaspow模型和Huilin-Gidaspow模型,BVK模型预测的气固两相运动最平缓。由于模拟的气固两相体系属于密相体系,Huilin-Gidaspow模型的光滑过渡函数没有产生效果,所以Gidaspow模型和Huilin-Gidaspow模型在各个方面的预测结果基本一致。     

鼓泡流化床中流动特性的多尺度数值模拟

鼓泡流化床因其较高的传热特性以及较好的相间接触已经被广泛应用于工业生产中,而对鼓泡流态化气固流动特性的充分认知是鼓泡流化床设计的关键.在鼓泡流化床中,气泡相和乳化相的同时存在使得床中呈现非均匀流动结构,而这种非均匀结构给鼓泡流化床的数值模拟造成了很大的误差.基于此,以气泡作为介尺度结构,建立了多尺度曳力消耗能量最小的稳定性条件,构建了适用于鼓泡流化床的多尺度气固相间曳力模型.结合双流体模型,对A类和B类颗粒的鼓泡流化床中气固流动特性进行了模拟研究,分析了气泡速度、气泡直径等参数的变化规律.研究表明,与传统的曳力模型相比,考虑气泡影响的多尺度气固相间曳力模型给出的曳力系数与颗粒浓度的关系是一条分布带,建立了控制体内曳力系数与局部结构参数之间的关系.通过模拟得到的颗粒浓度和速度与实验的比较可以发现,考虑气泡影响的多尺度曳力模型可以更好地再现实验结果.通过A类和B类颗粒的鼓泡床模拟研究发现,A类颗粒的鼓泡床模拟受多尺度曳力模型的影响更为显著.      

气固两相流中颗粒弥散的拉格朗日模拟

本文提出了一种对于均匀,稳定及各向同性气固两相紊流场中圆形固体颗粒弥散的拉格朗日模拟计算方法,应用该方法对带有网栅的垂直与水平管道中均匀,稳定的气固两相流模拟计算结果与Snyder及Wells等人所做的相同情况下的试验结果进行了比较,以证明该模拟计算方法的有效性,。     

气相爆轰波冲击气固界面的透反射特性

为了研究气相爆轰波冲击气固界面过程中透射波和反射波的相关特性,建立爆轰波冲击气固界面的一维理论模型,对不同初始压力条件下爆轰波到达气固界面后的界面两侧的压力和界面速度变化进行分析。利用时空守恒元求解元方法对气相爆轰波冲击气固界面过程进行数值模拟,分析气体部分反射波的压力分布和速度变化规律及透射入固体中应力波的波形和波速特征,并搭建气相爆轰波冲击活塞实验装置进行进一步验证。结果表明:气体爆轰波到达气固界面后,在固体中透射指数形式的弹性波,并在界面处向气体区反射一道激波。爆轰波后的稀疏波与反射激波相交,削弱反射激波,最终形成稳定激波回传。气固界面在稀疏波和反射稀疏波的作用下,压力和速度逐渐下降,最终也形成稳定状态。在不同混气初始压力情况下,爆轰波冲击过程中产生的最高压力和爆压的比值基本保持不变。理论模型对特征点相关物理量的计算值和实验数据符合的较好。     

FLOW DYNAMICS OF GAS-SOLID FLUIDIZED BEDS WITH EVAPORATIVE LIQUID INJECTION

1. Introduction Injections of evaporative liquids into fluidized solid parti- cles are routinely practiced in industrial processes involv- ing gas-solid fluidization systems such as fluid catalytic cracking, polymerization, and plastic coating (Fan et al., 2001). In the FCC riser system, heavy oil is injected into the system to evaporate rapidly by contact with the hot catalyst particles. Simultaneously, thermal and catalytic cracking reactions take place. During a polymerization process, a …     

On the onset of instability in the wake of super-hydrophobic spheres

We report an experimental investigation of free falling super-hydrophobic (SH) spheres in glycerine-water mixtures over a wide range of Reynolds number. SH coatings have the ability to reduce the contact area between the surrounding liquid and the solid surface by entrapping an air layer in the surface roughness. We investigate the effect of this air plastron on the hydrodynamic performance of spheres, focusing our attention on the onset of wake instabilities. Our results emphasise the key role of the surface roughness properties on the triggering of wake instabilities. It is shown that, unlike what was reported in previous numerical studies on SH bluff-bodies, local deformation of the interface may act as a by-pass mechanism promoting earlier transition, yielding a decrease of the critical Reynolds number at which the wake becomes unstable. A scenario coupling the hydrodynamic instabilities (scaling with the Reynolds number) and the interface deformation (scaling with the roughness-based Weber number) is proposed to describe the different transition mechanisms in presence of SH surfaces. It is found that the promotion of wake instabilities over SH surfaces occurs when the roughness-based Weber number is larger than a critical threshold. These findings are of primary importance for guiding the design of resilient and efficient SH surfaces.     

Bubble dynamics in a two-dimensional gas-solid fluidized bed

Related referential studies on gas-solid two-phase flows were briefly reviewed. Bubble ascending in a two-dimensional （2D） gas-solid fluidized bed was studied both experimentally and numerically. A modified continuum model expressed in the conservation form was used in numerical simulation. Solid-phase pressure was modeled via local sound speed; gas-phase turbulence was described by the K-ε two-equation model. The modified implicit multiphase formulation （IMF） scheme was used to solve the model equations in 2D Cartesian/cylindrical coordinates. The bubble ascending velocity and particle motion in the 2D fluidized bed were measured using the photochromic dye activation （PDA） technique, which was based on UV light activation of particles impregnated with the dye. Effects of bed height and superficial gas velocity on bubble formation and ascent were investigated numerically. The numerically obtained bubble ascending velocities were compared with experimental measurements. Gas bubble in jetting gas-solids fluidized bed was also simulated numerically.     

Impuls- und Wärmeaustausch in Zweiphasenströmungen

Zusammenfassung Nach Schilderung der Probleme, die bei der Vorausberechnung von Zweiphasenströmungen auftreten, wird ein Zweiphasenkontinuumsmodell diskutiert. Die Wechselwirkungen zwischen beiden Phasen werden dabei durch eine hydrodynamische (in den Impulsgleichungen) und eine thermische Deborah-Zahl (in den Energiegleichungen) erfaßt. Als Anwendungsbeispiele werden Druckabfall und Wärmeübergang einer turbulenten Gas-Feststoffströmung im senkrechten Rohr behandelt.

---

Momentum and heat transfer in two-phase flows

After describing the problems concerning the calculation of two-phase flows a two-phase continuum model is discussed. The interactions between both phases are considered by a hydrodynamic (in the momentum equations) and a thermal Deborah-number (in the energy equations). As an example the pressure loss and the heat transfer of a turbulent gas-solid flow in a vertical pipe are treated.

---

---

Herrn Prof. Dr.-Ing. K. Stephan zum 60. Geburtstag gewidmet     

A kinetic–hydrodynamic simulation of microstructure of liquid crystal polymers in plane shear flow

We study the microstructure formation and defects dynamics arising in liquid crystalline polymers (LCPs) in plane shear flow by a kinetic–hydrodynamic coupled model. The kinetic model is an extension of the Doi theory with a non-local intermolecular potential, including translational diffusion and density variation. LCP molecules are ensured anchoring at the boundary by an additional boundary potential, meanwhile mass conservation of LCPs holds in the whole flow region. Plane Couette flow and Poiseuille flow are studied using the kinetic–hydrodynamic model and the molecular director is restricted in the shear plane. In plane Couette flow, the numerical results predict seven in-plane flow modes, including four in-plane modes reported by Rey and Tsuji [Macromol. Theory Simul. 7 (1998) 623–639] and three new complicated in-plane modes with inner defects. Furthermore, some significant scaling properties were verified, such as the thickness of the boundary layer is proportional to molecular length, the tumbling period is proportional to the inverse of shear rate. In plane Poiseuille flow, the micro-morph is quasi-periodic in time when flow viscosity and molecular elasticity are comparable. Different local states, such as flow-aligning, tumbling or wagging, arise in different flow region. The difference of the local states, or difference of the tumbling rates in near-by regions causes defects and form branch pattern in director spatial–temporal configuration figure.     

Numerical simulation of gas-solid flow with two fluid model in a spouted-fluid bed

The flow characteristics in a spouted-fluid bed differ from those in spouted or fluidized beds because of the injection of the spouting gas and the introduction of a fluidizing gas. The flow behavior of gas-solid phases was predicted using the Eulerian-Eulerian two-fluid model （TFM） approach with kinetic theory for granular flow to obtain the flow patterns in spouted-fluid beds. The gas flux and gas incident angle have a significant influence on the porosity and particle concentration in gas-solid spouted-fluid beds. The fluidizing gas flux affects the flow behavior of particles in the fountain. In the spouted-fluid bed, the solids volume fraction is low in the spout and high in the annulus. However, the solids volume fraction is reduced near the wall.     

On the stability of Taylor bubbles inside a confined highly porous medium

This work focuses on the local hydrodynamics of a multiphase gas–liquid flow forced into an innovative medium of high porosity (96%): an open cell solid foam. The gas (nitrogen) and liquid (ethanol) phases are injected at constant flow-rates in a millichannel to form a well-controlled Taylor flow which enters the porous medium. Based on a fluorescence technique, the apparent liquid holdup in the porous medium is quantified, and its evolution in time and along the porous medium extracted from spatiotemporal diagrams. The analysis of the main frequency, when varying the gas–liquid flow-rate ratio, leads to the identification of two hydrodynamic regimes. A model based on a scaling analysis is proposed to quantify the dimensionless numbers describing the transition between both regimes. It points out that the bubble length fixed by the Taylor flow is the control parameter. The model prediction of the critical bubble length at which the transition occurs is in good agreement with the experimental observations.     

Explosive dispersal of solid particles

Abstract. The rapid dispersal of inert solid particles due to the detonation of a heterogeneous explosive, consisting of a packed bed of steel beads saturated with a liquid explosive, has been investigated experimentally and numerically. Detonation of the spherical charge generates a blast wave followed by a complex supersonic gas-solid flow in which, in some cases, the beads catch up to and penetrate the leading shock front. The interplay between the particle dynamics and the blast wave propagation was investigated experimentally as a function of the particle size (100–925 m) and charge diameter (8.9–21.2 cm) with flash X-ray radiography and blast wave instrumentation. The flow topology during the dispersal process ranges from a dense granular flow to a dilute gas-solid flow. Difficulties in the modeling of the high-speed gas-solid flow are discussed, and a heuristic model for the equation of state for the solid flow is developed. This model is incorporated into the Eulerian two-phase fluid model of Baer and Nunziato (1986) and simulations are carried out. The results of this investigation indicate that the crossing of the particles through the shock front strongly depends on the charge geometry, the charge size and the material density of the particles. Moreover, there exists a particle size limit below which the particles cannot penetrate the shock for the range of charge sizes considered. Above this limit, the distance required for the particles to overtake the shock is not very sensitive to the particle size but remains sensitive to the particle material density. Overall, excellent agreement was observed between the experimental and computational results. Received 16 August 1999 / Accepted 26 June 2000     

Numerical investigations on gas-solid flow in circulating fluidized bed risers using a new cluster-based drag model

A cluster-based drag model is proposed for the gas-solid circulating fluidized bed (CFB) riser by including the cluster information collected from image process...