CFD–DEM simulation of a pseudo-two-dimensional spouted bed comprising coarse particles

To investigate the energy dissipation mechanisms within the pump and improve the computational accuracy of the solid–liquid flow numerical simulations, in this study, an improved CFD-DEM (Computational Fluid Dynamics - Discrete Element Method) method has been presented. First, the improved method of CFD-DEM is introduced, which mainly considers the turbulent dissipation of particles in the near-wall region and velocity field reconstruction. Then, the simulation results before and after the method's enhancement are compared. Finally, the analysis of the energy characteristics of the liquid phase flow field in the solid–liquid flow is conducted. Research shows that the modified CFD-DEM method significantly improves the accuracy of the particle distribution predictions, with the numerical results for head and efficiency being much closer to experimental values. In the high-speed regions of the impeller flow field, primarily located behind the pressure side of the blades, the liquid phase flow velocity and pressure fluctuations are less affected by changes in solid phase concentration. In the fluid region of the centrifugal pump, the energy loss caused by entropy production is significantly concentrated in the volute and impeller regions. Specifically, the entropy production dissipation in the volute region accounts for the substantial portion of the total entropy production, approximately 67%–68%, while the entropy production dissipation in the impeller region accounts for about 19.7%–20.4%. As the solid phase concentration increases, the energy dissipation within the pump gradually rises, and the total vorticity at the impeller inlet also increases correspondingly, with the vorticity distribution being related to the number of blades. The findings provide a reference for further exploring solid–liquid flow within centrifugal pumps.

Computational fluid dynamics (CFD) combined with the discrete element method (DEM) are powerful tools for analyzing dense gas–solid flows. However, the computational cost of CFD–DEM will be unfeasibly great when simulating large-scale engineering applications with billions of particles. Accordingly, the coarse-grained (CG) CFD–DEM method is applied to solve this problem. This investigated method replaces several smaller particles with larger ones called parcels, aiming to reduce the number of particles and fully consider the collision of particles between composition parcels and the collision of particles within composition parcels. First, high-speed photography verifies the numerical simulation's reliability. Then, the CG CFD–DEM was used to analyze the transient spatial distribution, transient average velocity, pressure drop, bed height, and the mixing state of particles in a dense gas–solid fluidized bed. The CG CFD–DEM was also compared with the CFD–DEM results, which showed a good agreement with the calculation results and proved the accuracy and applicability of the method. Finally, the computation time of the CG CFD–DEM was evaluated, showing a significant decrease in computation time with an increasing coarse ratio (k). This investigation can provide theoretical reference for the numerical simulation of the CG CFD–DEM method in dense gas–solid flow.

Residence time distributions of different size particles in the spray zone of a Wurster fluid bed studied using DEM-CFD

Flow behavior of particles is simulated by means of discrete element method (DEM) in a screw conveyor. The gas flow is considered by the Euler model. The subgrid scale model turbulence model is applied to solve gas turbulence. A detailed description of the model equations applied has been presented. The comparison between computational fluid dynamics (CFD)‐DEM (i.e., considering gas flow) and DEM (i.e., not considering gas flow) simulation results is performed. Predicted mass flow rate with CFD‐DEM model is in better agreement with experiments by Roberts and Wills comparing with DEM results. The effects of incline angle and filling rate on the flow behavior of particles are predicted. The simulation results show that the gas flow has a major impact on particle flow behavior at a high rotation speed, whereas the velocity of particles increases with an increase of rotation speed. As the filling rate rises, both the translational granular temperature and the rotational granular particle temperature decline, which means that the particles fluctuation turns weak. While with an increase of incline angle, from the horizontal to vertical, the translational granular temperature is increased, but the rotational granular temperature reduces and then goes up. In brief, translational motion is a main mode of motion in a screw conveyor.

Bubbles as mesoscale structures play an important role in the development of the inhomogeneous drag force model, which is key to the prediction of fluidized beds. A coupled computational fluid dynamics–discrete element method is employed in this work to investigate the distribution of interphase drag force around a clouded bubble in an incipiently fluidized bed. A method for the determination of the bubble boundary is proposed. The results reveal that there is a stronger positive effect of drag force around the clouded bubble because of the presence of bubble crown and wake, which is different from the cloudless bubble. The bubble-induced local inhomogeneous coefficient is also proposed to evaluate the ensemble inhomogeneous impact around the bubble.

Particle shape effects on dynamic behaviors in a spouted bed: CFD-DEM study

To improve the operational efficiency of a hammer mill and delve into a high-efficiency, energy-saving grinding mechanism, the crucial parameters influencing the grinding of corn straw were identified as the spindle speed, hammer–sieve gap, and sieve pore diameter. According to the force analysis and kinematics analysis, the key factors affecting corn straw grinding were the spindle speed, the hammer–sieve gap, and the sieve pore diameter. The grinding process of corn straw was studied using computational fluid dynamics (CFDs) and the discrete element method (DEM) gas–solid coupling numerical simulation and experiment. The numerical simulation results showed that with the growth of time, the higher the spindle speed, the faster the bonds broke in each part, and the higher the grinding efficiency. When the energy loss of the hammer component was in the range of 985.6~1312.2 J, and the total collision force of the corn straw was greater than 47,032.5 N, the straw grinding effect was better, and the per kW·h yield was higher. The experimental results showed that the optimum combination of operating parameters was a spindle speed of 2625 r/min, a hammer-screen gap of 14 mm, and a sieve pore diameter of 8 mm. Finally, the CFD–DEM gas–solid coupling numerical simulation validation tests were performed based on the optimal combination of the operating parameters. The results showed that the energy loss of the hammer component was 1189.5 J, and the total collision force of the corn straw was 49,523.5 N, both of which were within the range of better results in terms of numerical simulation. Thus, the CFD–DEM gas–solid coupling numerical simulation could accurately predict the corn straw grinding process. This study provides a theoretical basis for improving a hammer mill’s key components and grinding performance. Meanwhile, the proposed gas–solid two-phase flow method provided theoretical references for other research in agricultural machinery.

The existence of defects in the enclosure structure is the primary cause of water and sand leakage in foundation pits, as well as being a significant source of danger in pit construction, but current research lacks an in-depth investigation of the generation mechanism and gestation process. In this paper, which comprehensively considers the microscopic particles and macroscopic level, the development mechanism of a water and sand leakage disaster in a foundation pit with a water-rich sand layer was studied using the principle of computational fluid dynamics and discrete element method coupled analysis (CFD–DEM); moreover, based on the anisotropy of the particle force and fluid energy analysis, the deformation of the stratum and ground stress field were analyzed. The results show that the stress field will produce a plugging effect at a certain distance from the defect, and the strata exhibit a dominant displacement tendency in the vertical direction, resulting in the emergence of a gradually concave stress relaxation zone and an elliptical contour in the strata displacement map near the defect. The fluid energy describes the displacement of the sand layer very well, and it is separated into the sand layer’s centralized loss region and the major loss area based on the high and low levels of the fluid energy class. The impact of fluid at the defect reaches the maximum kinetic energy, which penetrates the structural weakness and causes the loss of sand particles, and the cross-section of the water influx near the defect gradually expands with the loss of particles, indicating that there is a danger of further expansion of the defect under the impact of water flow. These results have technical implications for the management of water and sand leakage disasters in foundation pit engineering.

Acetylene hydrogenation reactions (AHR) are typically carried out in fixed bed reactors (FBRs) with high diameter‐to‐particle diameter (D/d). In this study, first, a real catalyst bed was built using the discrete element method, and the suitability of the bed was assessed using empirical equations. Second, a discrete element method‐computational fluid dynamics (DEM‐CFD) model was developed using CFD in conjunction with the reaction mechanism to simulate the AHR in the FBR and was compared with the experimental results. The results show that the DEM‐CFD can well predict the 3D interactions of velocity, temperature, and species for different D/d. The results can provide favorable theoretical guidance for the optimization of FBR and process intensification.

Recent studies highlight coarse-grained discrete element method (CGDEM) as a helpful tool for mitigating the computational cost associated with original discrete element method (DEM) simulations in large-scale industrial applications such as fluidized beds. This approach achieves cost savings by reducing the number of particles in the domain, while increasing their size. In the current work, CGDEM-LES (large-eddy simulation) numerical simulations are conducted on a 3D (three-dimensional) cylindrical pilot-scale fluidized bed in the bubbling regime, containing 9.6 M Geldart B-type particles. Macroscopic and mesoscopic analyses are performed, revealing qualitative and quantitative effects of coarse-graining on bed behavior. Among these effects, a global homogenization of the fluidized region is observed, marked by soaring bed surfaces, lower solid velocity, and solid fraction gradients. Additionally, a decrease in the bubble population is reported. These effects intensify as the coarse-graining factor increases. Despite influencing some results, the impact of mesh size is deemed negligible compared to that of particle coarse-graining. Some of the observed issues can be alleviated by incorporating additional mechanisms from the literature, aiming to dissipate the extra energy inherently present in coarse-grained systems. However, these mechanisms prove to be barely sufficient to replicate DEM results with the smallest coarse-graining factor tested. A thorough analysis allows identifying a side effect of one of these approaches, which is to slow particles down all the more as they move fast, causing a macroscopic misprediction of particle vertical velocity in turn. This model is then deemed less useful in the context of this study.

The slurry transport pump is the key equipment of deep-ocean mining systems. The motion law of coarse particles in the pump is not clear enough. In this paper, a hydraulic model of a laboratory-scale two-stage slurry transport pump is constructed, and the motion characteristics of coarse particles in the pump are numerically studied by using the computational fluid dynamics–discrete element method (CFD-DEM) method. The performance curve of the pump is obtained by experimental measurement, and the reliability of the calculated results is verified. Due to the application of the amplification flow rate design method, the optimum efficiency point of the pump is shifted to the large flow rate condition. Differences in particle swarm within two stages are compared. The position distribution, velocity variation and trajectory of particles in the impeller and bowl diffuser are studied in detail. The velocity of particles leaving the impeller depends on whether they collide with the impeller blade. The motion of particles in the bowl diffuser is divided into three periods. Collision between particles and blades in the bowl diffuser not only leads to energy loss but also gradually transforms the circumferential velocity of particles into axial velocity in the second period. This work can provide a reference for the study of wear and blockage prevention of slurry transport pumps.

Boulders in riverbeds significantly influence sediment transport and deposition by altering local flow patterns. This study employs coupled Computational Fluid Dynamics–Discrete Element Method (CFD–DEM) simulations to examine how boulder concentration affects sediment dynamics. Three boulder spacing scenarios, representing isolated, wake-interference, and skimming flow regimes, are evaluated for their impacts on flow structure, bed shear stress distribution, and sediment transport and deposition. The fluid phase is modelled using a large eddy simulation in a finite volume framework, while individual sediment grains are tracked with the discrete element method. The results reveal that increasing boulder concentration transforms the flow regime from isolated wakes behind individual boulders to a more coherent recirculation zone among neighboring boulders. This transition substantially reduces near-bed shear stresses between boulders and leads to a decline in total sediment transport rates. At high boulder concentration, sediment particles preferentially accumulate in sheltered inter-boulder corridors, forming stable depositional belts. In contrast, widely spaced boulders lead to isolated, localized deposition around individual boulders. To our knowledge, this is the first attempt to apply a coupled finite volume–DEM approach to simulate boulder–sediment interactions in open-channel flow, enabling analysis of sediment transport and deposition under varying boulder concentrations.

CFD-DEM investigation of the interaction between a particle swarm and a stationary bubble: Particle-bubble collision efficiency

This paper discusses the thermal and exergy efficiency analysis of the hydrothermal liquefaction (HTL) process, which converts sewage sludge into biocrude oil in a continuous plug–flow reactor using a linear Fresnel solar collector. The investigation focuses on the influence of key operational parameters, including slurry flow rate, temperature, pressure, residence time, and the external heat transfer coefficient, on the overall efficiency of biocrude oil production. A detailed thermodynamic evaluation was conducted using process simulation principles and a kinetic model to assess mass and energy balances within the HTL reaction, considering heat and mass momentum exchange in a multiphase system using UDF. The reactor’s receiver, a copper absorber tube, has a total length of 20 m and is designed in a coiled configuration from the base to enhance heat absorption efficiency. To optimize the thermal performance of biomass conversion in the HTL process, a Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) coupling numerical method approach was employed to investigate improved thermal performance by obtaining a heat source solely through solar energy. This numerical modeling approach allows for an in-depth assessment of heat transfer mechanisms and fluid-particle interactions, ensuring efficient energy utilization and sustainable process development. The findings contribute to advancing solar-driven HTL technologies by maximizing thermal efficiency and minimizing external energy requirements.

A novel pneumatic recovery method was proposed to curb the problem of high losses caused by side-cutting in a rapeseed combine harvester header. The influence of recovery method and material status changes on the recovery effect was studied via the CFD-DEM (computational fluid dynamic- discrete element method) coupling simulation. The effect of airflow action on the recovery effect was compared and analyzed, and the composite pneumatic recovery method was determined. In addition, the influence of material status changes and material feeding rate on the recovery effect was explored, and the critical condition of material blockage in the recovery device was configured. As such, the relationship model between air velocity and recovery rate was constructed and the air distribution ratio of the flow field in the device under this condition was optimized, had verified the rationality of this pneumatic recovery method was verified by a series of field tests. The average rapeseed recovery rate of 92.95% was achieved with the application of the recovery device, and the total loss rate of the header reduced by 52.26%, which is of great significance in reducing the total loss rate of the combine harvesters and improving the operation performance of machinery. The research results can provide a reference for the design of the header structure of a rape combine harvester. Keywords: rapeseed, combine harvester, side-cutting loss, coupled simulation, pneumatic recovery DOI: 10.25165/j.ijabe.20221502.6267 Citation: Jiang T, Zhang M, Guan Z H, Mu S L, Wu C Y, Wang G, et al. Simulation and analysis of the pneumatic recovery for side-cutting loss of combine harvesters with CFD-DEM coupling approach. Int J Agric & Biol Eng, 2022; 15(2): 117–126.