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

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

Influence of cycle time distribution on coating uniformity of particles in a spray fluidized bed by using CFD-DEM simulations

Cycle Time Distribution (CTD) plays a critical role for determining uniformity of particle coating in spray fluidized beds. However, the CTD is influenced by both geometrical structure and operating conditions of fluidized bed. In this study, a spray fluidized bed of coating process is simulated by a comprehensive Computational Fluid Dynamics-Discrete Element Model (CFD-DEM). To achieve different behaviors of CTD, some modifications are designed on a pseudo-2D internally circulating fluidized bed, which traditionally composes of a high-velocity upward bed and low-velocity downward bed. These modifications include making the air distributor slope and/or laying a baffle in the downward bed. First, the CTD and evolution of particle size distribution under different bed structures are compared. The CTD directly influences the coating uniformity. By making the particles flowing along a parallel direction in the downward bed through the geometrical modifications, the CTD becomes narrower and the coating uniformity is significantly improved. Second, under the optimized bed structure, the influence of operating conditions on the coating uniformity is studied. Properly increasing the fluidization gas velocity and the fluidization gas temperature and reducing the liquid spray rate can improve the coating uniformity.     

Fluid-dynamics evaluation in a conical spouted bed and characterization of foxtail millet seeds

Foxtail millet (Setaria italica) is one of the most valuable species in economic terms in the genus Setaria and plays an important role in human nutrition, animal feed, and agriculture. The present study described chemical, physical, and quality aspects of seeds of foxtail millet. Furthermore, the fluid-dynamic behavior of the seeds was evaluated in a conical spouted bed, which has advantages in terms of promoting the cyclic and regular movement of the seed particles. Dynamic parameters of spouting (minimum spouting velocity, stable and peak pressure drop) were determined and compared with those obtained from empirical correlations available in the literature. The results obtained from physical characterization showed that the seeds can be classified as belonging to Group D of Geldart, having a non-rough surface, mean diameter of 1.75 mm, and sphericity of 0.74. Fluid-dynamics analysis showed that the seeds are suitable for processing in a spouted bed, which is in agreement with the results of particle physical characterization.     

Mixing and segregation of solid particles in a conical spouted bed: Effect of particle size and density

In this work, the mixing and segregation of binary mixtures of particles with different sizes and densities in a pseudo-2D spouted bed were studied experimentally. A binary mixture of solid particles including sand, gypsum, and polyurethane was used. To determine the particles mass fraction, and their mixing and segregation in the bed, an image-processing technique was developed and used. Important hydrodynamic parameters, such as the axial and radial segregation profiles of the solid particles, were measured. The effects of air velocity, particle size, and particle mass fraction were also evaluated. The flow regime in the spouted bed and the time required for reaching the equilibrium state of the solid particles were discussed. The results showed that the segregation of solid particles and the time to equilibrium both decreased when the air velocity increased to much larger than the minimum spouting velocity. The axial segregation increased with the diameter ratio of the particles. Upon completion of the test, coarse particles were concentrated mainly in the spout region, while fine particles were aggregated in the annulus region. Examination of the flow pattern in the spouted bed showed that the particles near the wall had longer flow paths, while those near the spout region had shorter flow paths.     

Particle convective heat transfer near the wall in a supercritical water fluidized bed by single particle model coupled with CFD-DEM

Supercritical water fluidized bed (SCWFB) is a promising reactor to gasify biomass or coal. Its optimization design is closely related to wall-to-bed heat transfer, where particle convective heat transfer plays an important role. This paper evaluates the particle convective heat transfer coefficient (h_pc) at the wall in SCWFB using the single particle model. The critical parameters in the single particle model which is difficult to get experimentally are obtained by the computational fluid dynamics-discrete element method (CFD-DEM). The contact statistics related to particle-to-wall heat transfer, such as contact number and contact distance, are also presented. The results show that particle residence time (τ), as the key parameter to evaluate h_pc, is found to decrease with rising velocity, while increase with larger thermal boundary layer thickness. τ follows a gamma function initially adopted in the gas–solid fluidized bed, making it possible to evaluate h_pc in SCWFB by a simplified single particle model. The theoretical predicted h_pc tends to increase with rising thermal gradient thickness at a lower velocity (1.5 U_mf), while first decreases and then increases at higher velocity (1.75 and 2 U_mf). h_pc occupies 30%–57% of the overall wall-to-bed heat transfer coefficient for a particle diameter of 0.25 mm. The results are helpful to predict the overall wall-to-bed heat transfer coefficient in SCWFB combined with a reasonable fluid convective heat transfer model from a theoretical perspective.     

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.     

Penetration into a granular bed in the presence of upward gas flows

We present a numerical study on the penetration of spherical projectiles into a granular bed in the presence of upward gas flows. Due to the presence of interstitial fluid, the force chains between particles in the granular bed are weakened significantly, and this distinguishes the penetration behavior from that in the absence of fluid. An interesting phenomenon, namely granular jet, is observed during the penetration, and the mechanism for its formation and growth is attributed to the merging of granular vortices generated by the interaction between the intruder and primary particles. Moreover, both the final penetration depth and the maximum diameter of the crater are found to follow a power-law dependence with the impact velocity, and the maximum height reached by the granular jet tends to increase linearly as the impact velocity increases, agreeing well with the experimental results reported in the literature.     

DEM simulation of polydisperse systems of particles in a fluidized bed

Numerical simulations based on three-dimensional discrete element model (DEM) are conducted for mono-disperse, binary and ternary systems of particles in a fluidized bed. Fluid drag force acting on each particle depending on its size and relative velocity is assigned. The drag coefficient corresponding to Ergun’s correlation is applied to the system of fluidized bed with particle size ratios of 1:1 for the mono-disperse system, 1:1.2, 1:1.4 and 1:2 for the binary system and 1:1.33:2 for the ternary system b...     

Numerical simulation of fluid–particle hydrodynamics in a rectangular spouted vessel

A three-dimensional, Eulerian simulation was developed to describe isothermal, two-phase flow of the continuous (water) and dispersed (solid particles) phases in a rectangular spouted vessel. The mass and momentum conservation equations for each phase were solved using the finite volume technique, which treats each phase separately, while coupling them through drag, turbulence, and energy dissipation due to particle fluctuations. Particle–particle interactions via friction were also included.     

Numerical simulation of flow and heat transfer in a random packed bed

Random packed beds have more complex interior structure than structured beds and are widely used in industry and engineering.CFD simulation was carried out to investigate and analyze the local flow and heat transfer in a 120-sphere random packed bed.3D Navier-Stokes equation was solved with a finite volume formulation based on the Chimera meshing technique.Investigation was focused on low Reynolds number flow(Re=4.6-56.2),which typically occurs in packed bed reactors in bio-chemical fields.Detailed temperature field information was obtained.Inhomogeneity of flow and heat transfer due to the non-uniform distribution of void fraction was discussed and analyzed.     

Numerical simulation of methane-air turbulent jet flame using a new second-order moment model

A new second-order moment model for turbulent combustion is applied in the simulation of methane-air turbulent jet flame. The predicted results are compared with the experimental results and with those predicted using the wellknown EBU-Arrhenius model and the original second-order moment model. The comparison shows the advantage of the new model that it requires almost the same computational storage and time as that of the original second-order moment model, but its modeling results are in better agreement with experiments than those using other models. Hence, the new second-order moment model is promising in modeling turbulent combustion with NO_x formation with finite reaction rate for engineering application. The project sponsored by the Foundation for Doctorate Thesis of Tsinghua University, and the National Key Project in 1999–2004 sponsored by the Ministry of Science and Technology of China     

直管内胞格爆轰的基元反应数值研究

基于基元反应和二维欧拉方程,对直管内胞格爆轰进行了数值模拟。采用5阶WENO(weighted essentially nonoscillatory scheme)求解对流项,采用2阶附加半隐的龙格-库塔法处理化学反应源相引起的刚性。获得了密度、压力、温度和典型组元质量分数流场及数值胞格结构等。结果表明：网格精度的差异明显影响胞格的规则性和爆轰的平衡模数,随着网格尺度的减小,胞格由不规则变为规则。预混气组成、初压、初温及管道宽度给定,三波点数收敛为确定值。足够强度的初始扰动可再现胞格爆轰,最终形成的自持胞格爆轰模数与初始扰动的形状、大小、位置均无关。沿胞格中心线,爆轰波速度变化范围为0.88DCJ~1.5DCJ,爆轰波平均速度与CJ爆轰速度仅偏差0.88％。峰值压力与初压之比为14~50。计算爆轰波平均速度、胞格宽长比与实验值基本一致,但计算胞格宽度比实验值略小。数值模拟加深了对横波的产生和发展、未反应气囊、爆轰胞格的二次起爆等胞格爆轰特性的认识。     

Effect of a pulsating flow on hydrodynamics and fluid catalytic cracking chemical reaction in a circulating fluidized bed riser

The three-dimensional computational fluid dynamics (3D-CFD) of a pulsating flow applied to the fluid catalytic cracking (FCC) reaction was investigated in the riser of a circulating fluidized bed reactor. The kinetic parameters of the FCC and coke burning reactions for predicting the reactant conversion and product yield percentages were applied. To increase the reactant conversion level and product yield, the effect of the pulsating flow operating parameters was considered using a 2^k statistical experimental design with four factors (amplitude, frequency, types of the waveform, and amplitude ratio). The 3D-CFD simulation was successfully validated from the experimental literature data. The frequency and type of the waveform were found to be the significant operating parameters. The expression of the fitted regression model and response surface contour were derived and revealed that the pulsating flow provides a higher reactant conversion level and product yield percentages compared to a non-pulsating or steady flow.     

Numerical simulation of large amplitude oscillatory shear of a high-density polyethylene melt using the MSF model

We study the flow response in large amplitude oscillatory shear of the molecular stress function (MSF) model that has recently been proposed by Wagner et al. [M.H. Wagner, P. Rubio, H. Bastian, The molecular stress function model for polydisperse polymer melts with dissipative convective constraint release, J. Rheol. 45 (2001) 1387–1412]. The MSF model is derived from molecular theory and has only two parameters to describe the non-linear material response. The model predictions are analysed in both the frequency and time domain. It shows good agreement with experimental data for a linear high-density polyethylene melt. At low and medium strains, MSF model predictions are in excellent agreement with experimental data and predictions of a six-mode Giesekus model which has six parameters to describe the non-linear material response. At medium strains, the basic Doi–Edwards model, which has no non-linear parameters, already underpredicts the data. At high strains, the MSF model predictions agree slightly better with the experimental data than the Giesekus model. Surprisingly, however, it is the Doi–Edwards model that shows excellent agreement with experimental data at high strains. For the linear melt we consider, it outperforms the models that have non-linear parameters, both in the time and frequency domain.     

CFD-DEM simulation of the hole cleaning process in a deviated well drilling: The effects of particle shape

We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics (CFD) with the discrete element method (DEM). This numerical method allows us to incorporate the fluid–particle interactions (drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force) using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle−particle, particle−wall, and particle−drill pipe are taken into account with the Hertz–Mindlin model. We compare the transport of spheres with non-spherical particles (non-smooth sphere, disc, and cubic) constructed via the multi-sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid–particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid–solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.     

Numerical simulation of storm surges in Bangladesh using a multi-level model

A comparative study with a vertically integrated model and the multi-level model has shown that the former can remain a good substitute for the latter in the prediction of sea-surface elevation as long as the bottom friction coefficients in the vertically integrated models are properly specified.     

Numerical simulation and experimental validation of gas-solid flow in the riser of a dense fluidized bed reactor

-Ladd models. As a result, the proposed drag force model can be used as an efficient approach for the dense gas-solid two-phase flow.     

Three-dimensional multi-phase simulation of the mixing and segregation of binary particle mixtures in a two-jet spout fluidized bed

This study presents a three-dimensional numerical study of the mixing and segregation of binary particle mixtures in a two-jet spout fluidized bed based on an Eulerian–Eulerian three-fluid model. Initially, the particle mixtures were premixed and packed in a rectangular fluidized bed. As the calculation began, the gas stream was injected into the bed from the distributor and jet nozzles. The model was validated by comparing the simulated jet penetration depths with corresponding experimental data. The main features of the complex gas–solid flow behaviors and the mechanism of mixing and segregation of the binary mixtures were analyzed. Moreover, further simulations were carried out to evaluate the effects of operating conditions on the mixing and segregation of binary particle mixtures. The results illustrate that mixing can be enhanced by increasing the jet velocity or enlarging the difference of initial proportions of binary particle mixtures.     

Modelling the flow of power law fluids in a packed bed using a volume-averaged equation of motion

The development of a theoretical model for the prediction of velocity and pressure drop for the flow of a viscous power law fluid through a bed packed with uniform spherical particles is presented. The model is developed by volume averaging the equation of motion. A porous microstructure model based on a cell model is used. Numerical solution of the resulting equation is effected using a penalty Galerkin finite element method. Experimental pressure drop values for dilute solutions of carboxymethylcellulose flowing in narrow tubes packed with uniformly sized spherical particles are compared to theoretical predictions over a range of operating conditions. Overall agreement between experimental and theoretical values is within 15%. The extra pressure drop due to the presence of the wall is incorporated directly into the model through the application of the no-slip boundary condition at the container wall. The extra pressure drop reaches a maximum of about 10% of the bed pressure drop without wall effect. The wall effect increases as the ratio of tube diameter to particle diameter decreases, as the Reynolds number decreases and as the power law index increases.     

混凝土爆破损伤的SPH-FEM耦合法数值模拟

为了提高计算效率以及更好展现爆炸荷载下混凝土破坏过程,采用SPH-FEM耦合法对混凝土爆破成坑进行模拟。首先结合前人给出的C30混凝土Holmquist-Johnson-Cook（HJC）部分本构参数,通过理论推导等方法确定出剩余的参数;然后代入模型中计算,将数值解与实测数据进行对比;最后以峰值压力和峰值加速度作为考察对象,对HJC模型中21个参数敏感性进行分析。结果表明：SPH-FEM耦合法能直观地模拟爆炸荷载作用下爆坑的发展全过程,且能够较好地处理SPH边界问题;基于所给出的C30混凝土HJC本构参数,采用SPH-FEM耦合法对混凝土爆破破坏进行模拟,计算结果与实测数据吻合度高,表明HJC本构参数的确定具有合理性。此外,还发现HJC本构参数对爆破问题结果的敏感度各不相同,指出对峰值压力和峰值加速度均有较大影响的参数在确定的时候需引起足够的重视。