Search

• Semi-resolved CFD–DEM is extended to particulate flows with thermal convection. • Background fluid properties are corrected using kernel-based approximation. • Model A is adopted to address fluid-particle interaction with more force models. • Turbulence model is incorporated into the semi-resolved CFD–DEM. • The semi-resolved CFD–DEM provides better predictions and more details in modeling thermal particulate flows in fluidized bed. • Several phenomena observed in thermal fluidized bed are discussed. Particulate flows with heat transfer widely exist in industry including fluidized bed reactor and powder-based 3D printing, and are frequently modeled with CFD–DEM coupled approaches. Recently, we developed a semi-resolved CFD–DEM approach, which bridges the simulation gap between the resolved and unresolved CFD–DEM (Wang et al., 2019). The semi-resolved CFD–DEM is as efficient as the conventional unresolved CFD–DEM and as accurate as the resolved CFD–DEM, while it is validated for particulate flows with momentum exchange with Model B without heat exchange. In this paper, we further extend the semi-resolved CFD–DEM approach to model particulate flows with momentum exchange and thermal convection. Firstly, similar to background fluid velocity, background fluid temperature is corrected through kernel-based approximations on the neighboring fluid cells rather than simply taking values in the local cell containing the concerned particle. Secondly, Model A is adopted in the semi-resolved CFD–DEM, leading to more powerful capabilities in handling complicated flows. Force models like the Magnus force, virtual mass force are also included. Thirdly, turbulence model is incorporated into the semi-resolved CFD–DEM approach to address possible turbulence effects. Finally, the developed semi-resolved CFD–DEM is applied to modeling the complex particulate flows in fluidized bed. The obtained numerical results are compared with experimental data and results from other sources. It is demonstrated that the present semi-resolved CFD–DEM is effective in modeling particulate flows with thermal convection. Compared to the conventional unresolved CFD–DEM, this semi-resolved CFD–DEM can provide better predictions and more details on particle temperature distribution, mean temperature, solid/fluid fraction, velocity autocorrelation and vorticity field. [ABSTRACT FROM AUTHOR]

• Natural thermal convection induces fluid motion. • Particle shape-effects do play an important role in the particle and overall system dynamics. • Thermal heat exchange significantly affects the motion of the two-phase flow. • Multi-grid FEM-FBM is robust and potentially powerful tool for real particulate flow simulations. In particulate flows, particle-particle collisions play a very important role in determining the flow behavior of the fluid-particle two-phase systems. Thus, it is very important in the numerical simulations of particulate flows to treat the particle-particle interactions with a felicitous method. In this paper, a recently developed direct numerical simulation (DNS) technique, the Finite Element Fictitious Boundary Method (FEM-FBM), for thermal convective particulate flows is used for the simulations of dense particulate suspensions of elliptic shaped particles. The momentum and temperature flow fields are coupled with the aid of Boussinesq approximation. The thermal and momentum interactions between solid and fluid phases are handled by using the Fictitious boundary method (FBM). The continuity, momentum, and energy equations are solved on a fixed Eulerian mesh which is independent of flow features by using a multi-grid finite element scheme. A modified collision model is proposed that can handle not only the interactions between the circular particles but also effectively treat the collisions between the elliptic shaped particles. Firstly, we validate the newly developed collision model for isothermal circular particles together with a comparative study of "drafting, kissing and tumbling" (DKT) motion of both elliptic and circular shaped isothermal particles. Then we investigated the DKT motions of the hot and cold elliptic particles with energy exchange and studied the lateral behavior of the particles in numerous settings. Further, we studied the DKT motion of catalyst particles and investigate the particles behavior at different Grashof numbers. Numerical tests are performed to show that the present method is robust and provides an efficient approach for the simulations of particulate flows with a large number of elliptic particles. [ABSTRACT FROM AUTHOR]

This paper presents a direct numerical simulation (DNS) technique, the Finite Element Fictitious Boundary Method (FEM-FBM) for the simulation of fluid-solid two-phase flows with heat transfer. The heat transfer equation is introduced to study thermal convection in fluid-solid two-phase flows. The Boussinesq approximation is considered for the coupling of momentum and temperature flow fields. Multi-grid finite element solver is used to compute flow equations of mass, momentum, and energy on a fixed Eulerian mesh which is independent of time and the solid particles are allowed to move freely in the whole computational domain. Fictitious boundary method (FBM) is used to treat the particles inside the fluid, FBM takes account of the thermal and momentum interaction between the fluid and the particles. The accuracy and stability of presented method are validated by comparing our test cases results with the results reported in the available literature. Numerical tests are performed to show that this method is potentially powerful and provides an efficient approach to simulate complex thermal convective particulate flows with a large number of particles. [ABSTRACT FROM AUTHOR]

• The semi-resolved CFD-DEM with a multi-smoothing distance scheme is further developed for the polydispersed particle system. • Grain-scale reconstruction by μCT images is employed to quantitatively characterize the heterogeneous porous media. • Particle-scale migration process and distribution characteristics of fine particles in porous media are examined. • The influences of powder properties and anisotropy effects on migration processes and heat transfer are studied. A comprehensive understanding of the migrating phenomena of microscale particles associated with fluid flows and heat exchange in porous media is essential for developing unconventional geo-resources. In this work, we further develop the kernel-based semi-resolved CFD-DEM and combine it with grain-scale reconstruction to numerically investigate particulate flows and associated heat transfer processes in porous media. To characterize wide grain size distributions and heterogeneity of rock skeleton, a 3D spherical-packed model is quantitatively assembled based on CT images of real rock. After a series of validations, it is found that the results from the improved semi-resolved CFD-DEM, including motion, heat transfer, and pressure drop behaviors of particulate systems, match better with experimental observations and analytical solutions than that obtained from the unresolved CFD-DEM. On this basis, the particle-scale migration and distribution characteristics of fine particles are further studied. It is found that the penetration distance, temperature, and energy distributions of particle systems rely on the variation of powder properties (i.e., concentration and size) and the effects of matrix anisotropy. In a higher powder concentration system, more deposited particles at the upper of the packed bed prevent the downward migration of powder particles and lower the quantity of heat exchange between the fine and matrix particles. In addition, this can also lead to a more pronounced jamming phenomenon and reduce the seepage capacity of the packed bed. As the powder diameter increases, the heat transfer capacity of fine particles increases and gradually dominates the temperature of particulate systems. The deposition patterns and temperature variation behaviors of fine particles along the migration direction caused by the anisotropy of the porous medium are also identified. Therefore, the developed computational framework can be a powerful tool to reproduce the dynamics of fine particles in heterogeneous porous materials and reveal the underlying mechanisms. [ABSTRACT FROM AUTHOR]

• Four integrated SPH models are presented and compared for simulating natural convection. • The most suitable SPH model for natural convection is provided. • Natural convection for R a = 10 9 and Pr = 0.71 is successfully modeled for the first time by using the SPH method. • Mechanisms of the natural convection in a square cavity at different Rayleigh numbers are discussed. This paper conducts the simulation of natural convection in a differentially heated square cavity at high Rayleigh numbers by using the smoothed particle hydrodynamics (SPH) method. Due to the decrease of the accuracy and stability, it is challenging for the SPH method to simulate natural convection at high Rayleigh numbers, and there are few reported SPH literatures of natural convection at R a > 10 6 for air (Pr = 0.71). In this study, four integrated SPH models are presented to simulate the natural convection and their accuracy and stability are assessed. These four SPH models are associated with Kernel Gradient Correction (KGC) to improve approximation accuracy and Particle Shifting Technology (PST) to regularize particle distribution while they are different in treating density diffusion and calculating the pressure term. The numerical results show that SPH model_4 (KGC, PST, δ -SPH and asymmetric pressure approximation) is the most suitable for simulating the closed natural convection problems, especially at high Rayleigh numbers. Good agreements with reference solutions are obtained by SPH model_4 for the natural convection at 10 4 ≤ R a ≤ 10 8. Furthermore, the simulation of natural convection at R a = 10 9 is conducted by SPH model_4. The evolutions of thermal convection are described in detail. It is found that dynamics characteristic reveals that the dominant force is the pressure gradient, rather than the buoyancy force before the quasi-steady state. In addition, the chaotic motion at R a = 10 9 has significant influence to the heat transfer characteristic in the vertical boundary layers. [ABSTRACT FROM AUTHOR]