Effect of particle breakage on the behavior of simulated angular particle assemblies

Many attempts have been made to find various relationships for different parameters and some kinds of constitutive models for studying the behavior of particulate media. All these models are based on concepts of continuous media. Using a numerical method such as discrete element method, one can figure out what is happening through a discontinuous media where soil particles play the main role in introducing the shear strength and deformation characteristics. The behavior of the media with breakable particles is studied in this paper and compared with that of the assembly with non-breakable particles. In this paper, the hyperbolic elastic model is investigated for the assembly of polygon shaped particles in two different test series. In addition, evolution of different macro parameters of the assembly such as volume strain, angle of friction, angle of dilatancy and elastic modulus is studied during the simulation tests both for non-breakable and breakable soil particles. At the end, a parametric study is performed on the effect of strength of particle breakage on the assembly behavior.     

Modeling fluid–particle interaction in dilute-phase turbulent liquid–particle flow simulation

The present work examines the predictive capability of a two-fluid CFD model that is based on the kinetic theory of granular flow in simulating dilute-phase turbulent liquid–particle pipe flows in which the interstitial fluid effect on the particle fluctuating motion is significant. The impacts of employing different drag correlations and turbulence closure models to describe the fluid–particle interactions (i.e. drag force and long-range interaction) are examined at both the mean and fluctuating velocity l...     

A LBM–DEM solver for fast discrete particle simulation of particle–fluid flows

The lattice Boltzmann method (LBM) for simulating fluid phases was coupled with the discrete element method (DEM) for studying solid phases to formulate a novel solver for fast discrete particle simulation (DPS) of particle–fluid flows. The fluid hydrodynamics was obtained by solving LBM equations instead of solving the Navier–Stokes equation by the finite volume method (FVM). Interparticle and particle–wall collisions were determined by DEM. The new DPS solver was validated by simulating a three-dimensional gas–solid bubbling fluidized bed. The new solver was found to yield results faster than its FVM–DEM counterpart, with the increase in the domain-averaged gas volume fraction. Additionally, the scalability of the LBM–DEM DPS solver was superior to that of the FVM–DEM DPS solver in parallel computing. Thus, the LBM–DEM DPS solver is highly suitable for use in simulating dilute and large-scale particle–fluid flows.     

Measurement on particle rotation speed in gas–solid flow based on identification of particle rotation axis

When a two-dimensional (2D) imaging system is used to visualize particle motion in a 3D gas–solid flow, the particle rotation speed was found extremely difficult to be accurately measured due to the fact that the direction of rotation axis was usually random and hard to be distinguished. The paper presents a method to calculate the particle rotation speed from particle images based on the identification of its rotation axis using two or more characteristic points on its surface. The idea was analyzed and realized in a mathematical way and based on which a calculation program was given. The measurement method was verified with an experiment using a small sphere with known rotation axis and rotation speed. The effects of several factors, including the direction of the particle rotation axis, the particle image resolution, the types and positions of characteristic points, etc., on the measurement error are discussed. The error is found to be acceptable for most cases. The measurement method was finally applied to those small glass beads in a real 3D gas–solid flow inside a cold circulating fluidized bed (CFB) riser, which indicates that the problems of 2D imaging system applying to 3D particulate system could be solved by using this mathematical method.     

XPIV–Multi-plane stereoscopic particle image velocimetry

We introduce the three-dimensional measurement technique (XPIV) based on a Particle Image Velocimetry (PIV) system. The technique provides three-dimensional and statistically significant velocity data. The main principle of the technique lies in the combination of defocus, stereoscopic and multi-plane illumination concepts. Preliminary results of the turbulent boundary layer in a flume are presented. The quality of the velocity data is evaluated by using the velocity profiles and relative turbulent intensity of the boundary layer. The analysis indicates that the XPIV is a reliable experimental tool for three-dimensional fluid velocity measurements.More information at:

---

Cellular detonation diffraction in gas–particle mixtures

Diffraction of cellular heterogeneous detonation out of a channel into open half-space in a mixture of aluminum particles and oxygen is investigated numerically. The flow is found to be very similar to gas detonation diffraction. The detonation weakening behind the step results in combustion front deceleration and decoupling from the leading shock wave. Subsequent re-initiation takes place in a transverse wave. New transverse waves are generated along the expanding front. The computations that were performed show that the critical number of cells is several times less than that for gases. This is confirmed by theoretical estimates based upon the Mitrofanov–Soloukhin approach.     

Influence of a homogeneous axial magnetic field on Taylor–Couette flow of ferrofluids with low particle–particle interaction

Experimental results concerning the stability of Couette flow of ferrofluids under magnetic field influence are presented. The fluid cell of the Taylor–Couette system is subject to a homogeneous axial magnetic field and the axial flow profiles are measured by ultrasound Doppler velocimetry. It has been found that an axial magnetic field stabilizes the Couette flow. This effect decreases with a rotating outer cylinder. Moreover, it could be observed that lower axial wave numbers are more stable at a higher axial magnetic field strength. Since the used ferrofluid shows a negligible particle–particle interaction, the observed effects are considered to be solely based on the hindrance of free particle rotation.     

Transient phenomena in one-dimensional compressible gas–particle flows

The transient behavior of compressible gas– particle flows produced in shock tubes with particle-laden driver section is studied. Particular attention is focused on the time scales with which the solution approaches the equilibrium state. Theoretical estimates indicate that the gas and particle contact surfaces equilibrate first, followed by the shock wave, and finally by the expansion fan. The estimates are in good agreement with numerical simulations. The simulations also show that the approach to equilibrium condition of the shock speed is non-monotonic (monotonic) if the mass fraction of particles initially located in the driver section is below (above) a particle-diameter dependent critical value. For the speed of the particle contact surface, the reverse trends are observed.      

Discrete particle modeling of granular Rayleigh–Taylor instability

The gravitational air–grain Rayleigh–Taylor (RT) flow instability in a Hele-Shaw cell was studied using a parallel three-dimensional discrete particle model (DPM). The onset of flow instability and the development of fingering flow structures were well captured by the model. Power spectra analysis of solid volume fraction field indicated the non-linear coarsening process of the fingering flow structures. The sensitivity of the flow patterns to the initial porosity, the Atwood number, and the ratio of particle size to the Hele-Shaw cell width was also demonstrated. The excellent agreement of DPM simulation results with the reported experimental observations proved the robustness and reliability of the numerical approach to model complex multiphase flows such as granular RT instability.     

Numerical investigation of fluid–particle interactions for embolic stroke

Roughly one-third of all strokes are caused by an embolus traveling to a cerebral artery and blocking blood flow in the brain. The objective of this study is to gain a detailed understanding of the dynamics of embolic particles within arteries. Patient computed tomography image is used to construct a three-dimensional model of the carotid bifurcation. An idealized carotid bifurcation model of same vessel diameters was also constructed for comparison. Blood flow velocities and embolic particle trajectories are resolved using a coupled Euler–Lagrange approach. Blood is modeled as a Newtonian fluid, discretized using the finite volume method, with physiologically appropriate inflow and outflow boundary conditions. The embolus trajectory is modeled using Lagrangian particle equations accounting for embolus interaction with blood as well as vessel wall. Both one- and two-way fluid–particle coupling are considered, the latter being implemented using momentum sources augmented to the discretized flow equations. It was observed that for small-to-moderate particle sizes (relative to vessel diameters), the estimated particle distribution ratio—with and without the inclusion of two-way fluid–particle momentum exchange—were found to be similar. The maximum observed differences in distribution ratio with and without the coupling were found to be higher for the idealized bifurcation model. Additionally, the distribution was found to be reasonably matching the volumetric flow distribution for the idealized model, while a notable deviation from volumetric flow was observed in the anatomical model. It was also observed from an analysis of particle path lines that particle interaction with helical flow, characteristic of anatomical vasculature models, could play a prominent role in transport of embolic particle. The results indicate therefore that flow helicity could be an important hemodynamic indicator for analysis of embolus particle transport. Additionally, in the presence of helical flow, and vessel curvature, inclusion of two-way momentum exchange was found to have a secondary effect for transporting small to moderate embolus particles—and one-way coupling could be used as a reasonable approximation, thereby causing substantial savings in computational resources.     

Modeling of the unsteady force for shock–particle interaction

The interaction between a particle and a shock wave leads to unsteady forces that can be an order of magnitude larger than the quasi-steady force in the flow field behind the shock wave. Simple models for the unsteady force have so far not been proposed because of the complicated flow field during the interaction. Here, a simple model is presented based on the work of Parmar et al. (Phil Trans R Soc A 366:2161–2175, 2008). Comparisons with experimental and computational data for both stationary spheres and spheres set in motion by shock waves show good agreement in terms of the magnitude of the peak and the duration of the unsteady force.      

A combined viscoelastic–viscoplastic behavior of particle reinforced composites

This study formulates a micromechanical model for predicting effective viscoelastic–viscoplastic responses of composites. The studied composites consist of solid spherical particle reinforcements dispersed in a homogeneous matrix. The particle constituent is assumed linear elastic, while the matrix exhibits combined viscoelastic–viscoplastic responses. The Schapery integral model is used for the 3D isotropic non-linear viscoelastic responses. Two viscoplastic models are considered: the Perzyna model, having a rate-independent yield surface and an overstress function, and the Valanis endochronic model based on an irreversible thermodynamics without a yield surface. The Valanis model is suitable for materials when viscoplastic responses occur at early loadings (small stress levels). A unit-cell model with four particle and polymer sub-cells is generated to obtain homogenized responses of the particle-reinforced composites. Available micromechanical models and experimental data in the literature are used to verify the proposed micromechanical model in predicting effective time-dependent and inelastic responses of composites. Field variables in the homogenized composites are compared to the ones in heterogeneous composites. The heterogeneous composites, having detailed particle geometries, are modeled using finite element (FE) method.     

Positron emission particle tracking——Application and labelling techniques

The positron emission particle tracking （PEPT） technique has been widely used in science and engineering to obtain detailed information on the motion and flow fields of fluids or granular materials in multiphase systems, for example, fluids in rock cracks, chemical reactors and food processors; dynamic behaviour of granular materials in chemical reactors, granulators, mixers, dryers, rotating kilns and ball mills. The information obtained by the PEPT technique can be used to optimise the design, operational conditions for a wide range of industrial process systems, and to evaluate modelling work. The technique is based on tracking radioactively labelled particles （up to three particles） by detecting the pairs of backto-back 511 keV γ-rays arising from annihilation of emitted positrons. It therefore involves a positron camera, location algorithms for calculating the tracer location and speed, and tracer labelling techniques. This paper will review the particle tracking technique from tracking algorithm, tracer labelling to their application.     

Two-phase turbulence models for simulating dense gas–particle flows

The two-fluid model is widely adopted in simulations of dense gas–particle flows in engineering facilities. Present two-phase turbulence models for two-fluid modeling are isotropic. However, turbulence in actual gas–particle flows is not isotropic. Moreover, in these models the two-phase velocity correlation is closed using dimensional analysis, leading to discrepancies between the numerical results, theoretical analysis and experiments. To rectify this problem, some two-phase turbulence models were proposed by the authors and are applied to simulate dense gas–particle flows in downers, risers, and horizontal channels; Experimental results validate the simulation results. Among these models the USM-Θ and the two-scale USM models are shown to give a better account of both anisotropic particle turbulence and particle–particle collision using the transport equation model for the two-phase velocity correlation.     

Validation of fluid–particle interaction force relationships in binary-solid suspensions

In this work several relationships governing solid–fluid dynamic interaction forces were validated against experimental data for a single particle settling in a suspension of other smaller particles. It was observed that force relationships based on Lattice-Boltzmann simulations did not perform as well as other interaction types tested. Nonetheless, it is apparent that, in the case of a suspension of different particle types, it is important that the correct choice is made as to how the contribution to the overall fluid–particle interaction force is split between buoyancy and drag. Experimental evidence clearly suggests that the “generalized” Archimedes’ principle (where the foreign particle is considered to displace the whole suspension and not just the fluid) provides the best result.     

Order and disorder in particle–liquid flows in a pipe

The spherical expanded polystyrene particle–oil two-phase flow in a vertical pipe was used to simulate the dispersed phase distribution in laminar bubbly flows. A three-dimensional particle image tracking technique was used to track the particles in the flow to study the ordered structure of dispersed phase distribution and its transition to disorder. The ordered structures behaved as particle strings aligned in the flow direction as induced by the flow shear. The structures were quite durable in high liquid velocity flows and dispersed gradually as the liquid velocity decreased. In lower velocity flows, the particles tended to form clusters in the horizontal direction, as predicted by potential theory for spherical bubbles rising in a quiescent inviscid liquid and as observed in experiments on non-shear bubbly water flows.     

Thermal energy storage： Challenges and the role of particle technology

Thermal energy is at the heart of the whole energy chain providing a main linkage between the primary and secondary energy sources. Thermal energy storage （TES） has a pivotal role to play in the energy chain and hence in future low carbon economy. However, a competitive TES technology requires a number of scientific and technological challenges to be addressed including TES materials, TES components and devices, and integration of TES devices with energy networks and associated dynamic optimization. This paper provides a perspective of TES technology with a focus on TES materials challenges using molten salts based phase change materials for medium and high temperature applications. Two key challenges for the molten salt based TES materials are chemical incompatibility and low thermal conductivity. The use of composite materials provides an avenue to meeting the challenges. Such composite materials consist of a phase change material, a structural supporting material, and a thermal conductivity enhancement material. The properties of the supporting material could determine the dispersion of the thermal con- ductivity enhancement material in the salt. A right combination of the salt, the structural supporting material, and the thermal conductivity enhancement material could give a hierarchical structure that is able to encapsulate the molten salt and give a substantial enhancement in the thermal conductivity. Understanding of the structure-property relationships for the composite is essential for the formulation design and fabrication of the composite materials. Linking materials properties to the system level performance is recommended as a key future direction of research.     

Paradox of a particle’s trajectory moving on a string

This paper deals with the paradoxical properties of the solution of string vibration under a moving mass. The solutions published to date are not simple enough and cannot be applied to investigations in the entire range of mass speeds, including the overcritical range. We propose a formulation of the problem that allows us to reduce the problem to a second-order matrix differential equation. Its solution is characteristic of all features of the critical, subcritical, and overcritical motion. Results exhibit discontinuity of the mass trajectory at the end support point, which has not been previously reported in the literature. The closed solution in the case of a massless string is analyzed and the discontinuity is proved. Numerical results obtained for an inertial string demonstrate similar features. Small vibrations are analyzed, which is why the effect discussed in the paper is of purely mathematical interest. However, the phenomenon results in complexity in discrete solutions.     

Modelling of particle layer detachment—Role of transient kinetic effects

Particle layers tend to build up on walls in many filtration and separation processes, calling for periodic removal in order to keep the apparatus running. Important factors are the adhesion of the layer on the substrate and the cohesion of the particles in the layer. Models describing such layer detachment generally assume constant and homogeneous conditions for the forces acting on the layer. But in reality detachment is extremely non-stationary concerning place and time, primarily due to changing conditi...