Critical raw materials (CRMs) are one of the enablers of a sustainable future due to their importance in green technologies. Yet, their own circularity and end-of-life recycling rates have been lacking as their concentrations are too low in waste products to be efficiently recycled. This is not the case, however, for discarded printed circuit boards where different types of electronics components (ECs) use specific CRMs in high concentrations. Furthermore, due to worldwide manufacturing standards these ECs are consistent between different printed circuit boards (PCBs) in their physical characteristics such as size, shape and material composition. Yet at the moment no sorting methods exist that can separate ECs from modern PCBs. Therefore, we aim to evaluate multiple simple yet effective mechanical separation methods to sort said ECs with CRM recovery in mind. First of all, ECs from flat-panel displays were sieved into a small (<4 mm), medium (4 - 10 mm) and large fraction (>10 mm). This was followed by a roll sorter to separate thinner components, like IC chips, from similarly sized thicker ECs. In order to separate the components based on ferromagnetic composition an innovative overbelt ferromagnetic separator was developed where the magnetic field strength continuously decreased over the length of the belt. Lastly, the different types of ECs were analysed by laser-induced breakdown spectroscopy to identify the presence of any CRMs. Our study shows that by combining these three different sorting technologies it is possible to sort the ECs in a way that the majority of CRM containing components are concentrated in only 21.09 wt% of the total weight. This in turn results in significantly higher CRM concentrations, thus removing a major limitation to their recovery and improving CRM circularity for a more sustainable future.
@en

The agglomeration of cohesive particles can deteriorate fluidization quality and cause the defluidization of a bed, which is a common issue found in the applications of fluidized beds. This study aims to gain a better understanding of particle cohesion on agglomeration/fluidization behaviors and the effective methods for achieving a better fluidization quality, through numerical simulations based on the coupled approach of computational fluid dynamics and discrete element method (CFD-DEM). The effects of particle cohesion, gas velocities or flow conditions, and the bed geometry on the agglomeration and fluidization behaviors are analyzed. It is shown that the increase of particle cohesion can lead to deteriorated particle mixing, significant agglomeration of particles, and defluidization of the bed; the agglomeration-induced defluidization of highly cohesive particles is difficult to mitigate in a conventional flat-bottom fluidized bed. As large-sized agglomerates are more frequently found in the bottom of the bed, the spouted gas flow is then utilized and demonstrated to be effective in assisting the deagglomeration and fluidization of highly cohesive particles. Through the comparison of various spouted beds and spouted fluidized beds, the effective design of the bed bottom is identified for achieving a higher fluidization quality. Corresponding mechanisms underlying spout-assisted deagglomeration and fluidization are found to be much related to not only the enhanced particle-fluid but also particle-wall interactions in the confined space of a conical bed bottom, thus explaining the effectiveness and the importance of the bottom conical geometry of spouted beds. The obtained findings may help to understand the agglomeration-induced defluidization of fluidized beds and assist the fluidization of highly cohesive particles by the effective design of spouted beds.

@en

The high-precision scrap sorting for effective metal recycling can bring substantial environmental and economic benefits. This article presents a magnetic image sensor that can help to identify the ferrous contaminants inside nonferrous scraps of large sizes. First, the concept and the theory for detecting ferrous contaminants are described. In particular, an inversion algorithm is proposed to characterize the size and position of ferrous contaminants inside the main scrap bodies. Then, based on computed and measured results, the feasibility of sensor design using either 1-D Hall arrays or 1-D pickup coils is demonstrated. Finally, methods are suggested to minimize disturbing signals from very large nonferrous pieces passing through the slightly uneven magnetic field. The obtained findings in this study may apply not only to nonferrous scraps but many other materials of which the mass ratio of the ferrous contaminant to the main material is small.

@en

The growing demand for aluminium worldwide makes aluminium recycling critical to realising a circular economy and increasing the sustainability of our world. One effective way to improve the impact of aluminium recycling is to develop cost-efficient automated sorting technologies for obtaining pre-defined high-quality aluminium scrap products, thus reducing undesirable downcycling and increasing environmental/economic benefits. In this work, an innovative facility, which includes singulation, sensor scanning, and ejection, is optimised for the automated sorting of aluminium scraps. The sorting facility is computationally studied by a virtual experiment model based on the discrete element method. The model considers particle-scale dynamics of complex-shaped scraps and mimics the automated operation of the facility. Based on virtual experiment modelling, the flow of scrap is optimized by computation, with the feasible operation of the sorting facility being proposed. Accordingly, the sorting facility has been built and model predictions are confirmed in actual operation.@en

An in-depth exploration of the reaction kinetics and thermo-chemical behaviors of the raceway can offer practical insights for optimizing the operations of blast furnace (BF), thus achieving a more effective iron and steel production process. In this study, the dynamic characteristics and the flow, heat and mass transfer behaviors in the BF raceway were simulated by Discrete Element Method-Computational Fluid Dynamics (DEM-CFD) method at a particulate scale. The effects of coke size distribution and blast velocity on coke combustion characteristics, thermochemical behavior (particle volume fraction, raceway size, carbon loss, and coke temperature) and microscopic properties (coordination number (CN), contact normal force, pore structure and stress) were systematically investigated. The results show that as the blast velocity decreases or the size ratio λ (the largest coke particle size divided by the smallest coke particle size) increases, the raceway size becomes smaller, resulting in a smaller area of high oxygen (O₂) concentration and low carbon monoxide (CO) concentration in the raceway, and higher CO concentration in the packed bed. For the thermal-chemical behaviors, a lower blast velocity or a higher λ value decreases the number of particles experiencing mass loss, as well as increases individual particle mass loss, the average coke temperature and its variance. For microscopic properties, the CN distribution becomes wider as λ increases. The contact normal force in the coke bed with λ > 1 is significantly higher than that of λ = 1. As λ increases or blast velocity decreases, the pore distribution curve shifts to the left and the average pore volume decreases. The stress acting on the particles in the raceway increases with the blast velocity or λ. These new understandings of the complex reactive flow behaviors in the raceway will shed light on energy utilization and process optimization.

@en

Mixing structures and characteristics are crucial to the mixing quality of spherical/cylindrical binary granular systems like the biomass-coal mixtures which can affect energy release and carbon emissions. In this work, the mixing of binary spheres/cylinders in a rotating drum was numerically reproduced by using discrete element method (DEM). Systematic parametric studies were conducted to identify the role of various parameters such as rotation speed (ω) of the drum, aspect ratio (AR), mass fraction (φ_c) and volume fraction (φ_v) of cylindrical particles, and the density difference (φ_ρ) between the binary particles; meanwhile, corresponding mechanisms were also explored by analyzing the kinetic energy. Results show that when the flow regime is rolling/cascading, the mixing quality can be effectively improved; however, when the flow regime is cataracting, the mixing quality becomes worse. With AR close to 1.0, the interlocking effect between particles becomes weaker and the porosity becomes smaller, which leads to the higher contact efficiency and thus improves the mixing quality. Binary mixtures with different φ_c are synergistically affected by energy input and interlocking structure. The volume of spherical particles is more conducive to improving the mixing quality than that of cylindrical particles when the volume of the granular system is at the same level. The φ_ρ can cause segregation behavior of particles so as to make the mixing quality worse.

@en

The densification of mono-sized regular octahedral particles under vibration is simulated by a DEM model. The model is validated by physical experiments. The effects of operating parameters including vibration amplitude and frequency are investigated. The obtained packing structures are characterized, and corresponding densification mechanisms are identified. The results indicate that vibration amplitude is more critical for densification compared to frequency, where a wide range of suitable amplitudes exists to realize the dense packing structures. Ordered packing structure exists in boundary region. The central region possesses good particle position randomness and orientation randomness, and the packing density for poured packing (initial loose packing) and vibrated packing (dense packing) in the central region are 0.670 and 0.684, respectively. During densification, contacts and constraints between particles decrease, and the stress within the structure becomes larger. The typical arrangement of neighboring particles is close to but not perfect face to face (F–F) contact. Particles in the boundary region experience more rotations and can form local ordered structures. These structures prevent the relative movement of individual particles but favor their movement as a whole. Particles in the central region have relatively little movement and rotation. The ordered structures formed at boundary is difficult to extend to the central region due to the restricted change of particle contact conditions.

@en

A coupled approach of discrete element method and pore network model (DEM-PNM) is developed for simulating particle–fluid flow. By this approach, the particle/solid flow is described at a particle scale by DEM and the fluid flow is described at a pore scale by PNM. The coupling scheme between DEM and PNM and the boundary conditions for realizing particle–fluid flows under different conditions are explicitly introduced. The capability is examined in three systems: a dynamic cell, a gas-fluidized bed, and a packed bed with lateral gas injection. The simulated results by the DEM-PNM approach agree with those results by either DEM-LBM (lattice Boltzmann method) approach or experiments. In particular, the simulated reasonable flow patterns in gas fluidization and lateral gas injection demonstrate the applicability of the proposed approach to complex particle–fluid flow. However, the current model presents less accuracy in dilute flow region and does not consider some important factors, which should be further improved in the future.

@en

The packing of cohesive particles is of paramount importance in many industries because the packing structure is closely related to process performance. A general relation between packing density and interparticle force was previously proposed based on packing structures formed without dynamic fluid flows. Its universality is examined here in two different packings, formed in settling and defluidization of static and dynamic fluids, respectively. First, it is shown that the packings of the same particles formed by two different methods have different structures because of different impact-induced pressures. Nevertheless, a one-to-one relationship between packing density and structural properties still holds regardless of the different packing methods, and the force distribution in those packings obeys similar rules. Finally, the packing densities obtained by the different methods are demonstrated to be universally correlated with the ratio of the interparticle force to the effective gravity. These findings indicate that different phenomena of particulate systems at a macro- or meso-scale may share similar microscopic origins, with the interparticle force playing a crucial role.

@en

Discrete particle simulation can explicitly consider interparticle forces and obtain microscopic properties of the fluidized cohesive particles, but it is computationally expensive. It is thus pivotal to link the microscopic discrete properties to the macroscopic continuum description of the system for large scale applications. This work studies the fluidization of cohesive particles through the coupled computational fluid dynamics and discrete element method (CFD-DEM). First, discrete CFD-DEM results show the increased particle cohesion leads to the severe particle agglomeration which affects the fluidization quality significantly. Then, continuum properties are attained by a weighted time-volume averaging method, showing that tensile pressure becomes significant as particle cohesion increases. By incorporating Rumpf correlation into the solid pressure equation, the tensile pressure could be predicted consistently with the averaged CFD-DEM results for different particle cohesion. Finally, those overall steady averaged properties of the bed are obtained for understanding the general macroscopic properties of the system.

@en

To identify the dense packing of cylinder–sphere binary mixtures (spheres as filling objects), the densification process of such binary mixtures subjected to three-dimensional (3D) mechanical vibrations was experimentally studied. Various influential factors including vibration parameters (such as vibration time t, vibration amplitude A, frequency ω, vibration acceleration Γ) as well as particle size ratio r (small sphere vs. large cylinder), composition of the binary mixtures X_L (volume fraction of cylinders), and container size D (container diameter) on the packing density ρ were systematically investigated. The results show that the optimal vibration parameters for different binary cylinder–sphere mixtures are different. The smaller the size ratio, the less vibration acceleration is needed to form a stable dense packing. For each binary mixture, high packing density can be obtained when the volume fraction of large cylindrical particles is dominant. Meanwhile, increasing the container size can decrease the container wall effect and get higher packing density. The proposed analytical model has been proved to be valid in predicting the packing densification of current cylinder–sphere binary mixtures.

@en

In this paper, the transition from random to ordered packings of mono-sized cylindrical particles under 3D mechanical vibration was simulated by discrete element method (DEM). The effects of particle aspect ratio, size of the particulate system and container wall on the granular crystallization were investigated. And the mechanisms were analyzed through the characterization of the order transition process in terms of packing density, coordination number (CN), orientational ordering parameter (O), nucleation and growth of cluster and granular temperature (θ). The results show that nearly perfect crystallization of mono-sized cylindrical particles can be achieved with specific aspect ratio and proper vibration conditions in a cylindrical container. The orientational ordering parameter demonstrates that the crystallization firstly starts from the container wall and then propagates inward gradually. The lower granular temperature in a cuboid container indicates less vibration energy transferred to the granular assembly compared with that in a cylindrical container, illustrating the critical role of container shape in crystallization. The vibrated crystallization of mono-sized cylindrical particles is analogous to entropy-driven process, which can be gradually achieved by the self-assembly of particles with parallel alignments along the already formed ordered clusters (nuclei)or the container wall. These results can help design complex and ordered granular structures.

@en

Understanding the relationship between a fluid flow and the structural properties of a bed is important for various industrial applications and a proper description of particle-fluid flow. This work presents a pore-scale study of fluid flow through randomly packed beds by a network model. The model incorporates the inertial effect and thus can be applied to flows of high Reynolds numbers. The results show the structure heterogeneity associated with the connection of pores is important in determining the distribution of fluid flow. The statistical distribution of normalized drag forces on individual particles is Gaussian, related to the statistical distribution of pore properties. The mean drag force of a bed by the pore-scale model agrees well with that from a previous study by the lattice-Boltzmann method. These findings support the application of the pore network model for simulating particle-fluid flows under a wider range of Reynolds numbers.

@en

This paper presents a systematic numerical investigation of static stress profile within a confined granular packing using discrete element method. The vertical and horizontal stress profiles are carefully identified within the packing structure, and the effects of container diameter and friction coefficient on the stress distributions are systematically investigated and analyzed. The results show a quantitative agreement between the vertical stress profile gauged at the bottom and the stress profile calculated within the packing structure. I.e. beyond a hydrostatic region, the static vertical stresses saturate exponentially with the packing depth. A positive correlation between the saturation stress and the container size and a negative correlation with the friction coefficient can be observed. Besides, it is also found that the well-known Janssen coefficient (or Janssen constant) is not a constant value; it decreases asymptotically to a constant with the packing depth, which implies a non-constant friction effects existed between the granular assemblies and the boundaries. And the analysis of contact forces within the packing structure also indicates a stabilized condition beyond the hydrostatic region. It is believed that these findings could develop a comprehensive insight on the stress distribution of granular matter.

@en

The dynamic crystallization of cubic granular particles under three-dimensional mechanical vibration is numerically investigated by the discrete element method. The effects of operational conditions (vibration, container shape and system size) and particle properties (gravity and friction) on the formation of crystals and defects are discussed. The results show that the formation and growth of clusters with face-to-face aligned cubic particles can be easily realized under vibrations. Especially, a single crystal with both translational and orientational ordering can be reproduced in a rectangular container under appropriate vibrations. It is also found that the gravitational effect is beneficial for the ordering of a packing; the ordering of frictional particles can be improved significantly with an enlarged gravitational acceleration. The flat walls of a rectangular container facilitate the formation of orderly layered structures. The curved walls of a cylindrical container contribute to the formation of ring-like structures, whereas they also cause distortions and defects in the packing centers. Finally, it is shown that the crystallization of inelastic particles is basically accomplished by the pursuit of a better mechanical stability of the system, with decreasing kinetic and potential energies.

@en

Bottom pressure of confined granular assemblies saturates at a certain value even this packing bed is being continuously charged. Corresponding formulation has been established to describe this interesting phenomenon. In this work, the influences of particle size and friction on the bottom stresses of granular matter were numerically investigated by discrete element method (DEM). It is found that the Janssen model can well predict the stress profile only when the size ratio of the container versus the particle is larger than 16. Moreover, a hydrostatic linear relation between apparent mass and filling mass can be obtained when the friction coefficient becomes insignificant (μ ≤ 0.01). To further interpret the Janssen effects, the granular assemblies are characterized and evaluated from the overall interactions with sidewalls, angular distribution function, void size distribution, coordination number, contact networks, contact orientation and distributions of contact forces within the packing structure. It is believed that these results will be helpful to comprehend the granular behaviors and may offer instructive reference to industrial processes in related fields.@en

The packing densification of mono-sized equilateral cylindrical particles under mechanical vibration is numerically reproduced using discrete element method (DEM). The influences of vibration frequency, amplitude and container size on the macro property (e.g. packing density) of each packing are studied. Meanwhile, various micro-properties including coordination number (CN), radial distribution function (RDF), local structures, contact types, particle position/orientation distributions, forces/stresses of the vibrated dense packing are characterized and compared with those of the loose initial poured packing. The results show that properly controlling vibration conditions can realize the transition of equilateral cylindrical particles from random loose packing (RLP) to random close packing (RCP). The maximum packing density without wall effects can reach about 0.7166, which agrees with experimental and numerical results in literature. Micro property analyses demonstrate that the average CN increases slightly after vibration. The RDF curves indicate three obvious peaks for both poured and vibrated packings. The distributions of intersection angles imply that the perpendicular arrangements of the cylindrical particles are common in both packing structures. From RLP to RCP, the probability of side-edge, bottom-edge and edge-edge contacts between two particles decreases, while that of side-side, side-bottom and bottom-bottom contacts increases. Both particle position and orientation distributions illustrate that the structure in the center of the vibrated packing is disordered. The distribution of strong forces follows the exponential law, while that of weak forces follows the power law. The static vertical stresses in both packings can be predicted by Janssen model. After vibration, a much denser and more stable structure with a larger saturation stress is obtained due to the reducing inter-particle friction effects.

@en

Densification of cubic particles under three-dimensional (3D) mechanical vibration was studied experimentally. Effects of vibration time (t), frequency (f), vibration amplitude (A), vibration acceleration (Г), container size (D) and particle sphericity (φ) on packing density (ρ) were comprehensively analyzed. To identify the effects of particle shape, the packing densification of cuboid 1 (12 mm × 12 mm × 28 mm) and cuboid 2 (4 mm × 8 mm × 16 mm) were systematically studied under the same conditions. The results show that the structure of cubic particles can be first densified from random loose packing (RLP) to random close packing (RCP) and then to ordered packing (OP) gradually. In comparison, cuboid 1 and cuboid 2 particles can only form RCP structure, but cannot achieve OP even under further vibration. Vibration parameters (t, f, A and Г) are shown to be important to the packing densification. Based on the results of varying f and A, A - f phase diagrams are established for choosing the optimal vibration parameters to achieve the desired dense packing structures. Besides, it is shown the size of container (wall effect) has monotonic influence on packing density, i.e., the larger container size corresponds to the less wall effect and higher packing density. Cubic particles can form ordered packing because of the geometrical symmetry with aspect ratio of 1. After eliminating the wall effect by extrapolating the packing densities in different sized containers, typical packing densities of the three types of particles are obtained with ρ_RLP = 0.658, ρ_RCP = 0.830 and ρ_OP = 0.965 for cubic particles; ρ_RLP = 0.625 and ρ_RCP = 0.740 for cuboid 1 particles; ρ_RLP = 0.591 and ρ_RCP = 0.731 for cuboid 2 particles. These findings indicate cubic particles are efficient in densification which can be readily realized through proper 3D mechanical vibration.

@en

Bed structure is important to heat and mass transfer in fluidization. This study investigates solid structural transition for cohesive particles by the combined approach of computational fluid dynamics and discrete element method, supported by Voronoi tessellation. Macroscopic behaviors of fluidized beds are first discussed in connection with microscopic structural transitions. Then, the properties of Voronoi polyhedra are analyzed at both bed and particle scales. In addition, the correlation of coordination number and porosity is examined at a particle scale, and their combined probability distribution is analyzed to identify dominant structure of particles in different flow regimes. Finally, an index is proposed as the ratio of coordination number and porosity to differentiate microscopic states of a particle. These findings should be useful for the study of flow, heat, and mass transfer in fluidized beds and for the understanding of solid-like to fluid-like transition of granular materials.

@en

Packing densification of cubical particles under mechanical vibrations was dynamically simulated by using discrete element method (DEM). Effects of operating parameters such as vibration amplitude, frequency, vibration intensity, and container size (container wall) on packing densification were comprehensively investigated. The macro property such as packing density and micro properties such as coordination number (CN), radial distribution function (RDF), and particle orientation of the packings were analyzed and compared. It is found that mechanical vibration with proper vibration amplitude and frequency is effective for the densification of cubical particle packing. Packing structures of different packing densities display different properties, based on which random loose packing (RLP) and random close packing (RCP) of cubical particles are identified with the packing density of 0.591 and 0.683, respectively. Two densification mechanisms are discussed as the particle rearrangement is dominant for the transition from RLP to RCP and the crystallization along the container wall is dominant for the transition beyond RCP to ordering. The obtained results are useful for optimizing vibration conditions to generate dense packings and understanding the structural information of some fixed beds with cubical particles.

@en