Marco Ramaioli

Fluid–particle flow simulations using two-way-coupled mesoscale SPH–DEM and validation

First, a meshless simulation method is presented for multiphase fluid–particle flows with a two-way coupled Smoothed Particle Hydrodynamics (SPH) for the fluid and the Discrete Element Method (DEM) for the solid phase. The unresolved fluid model, based on the locally averaged Navier Stokes equations, is expected to be considerably faster than fully resolved models. Furthermore, in contrast to similar mesh-based Discrete Particle Models (DPMs), our purely particle-based method enjoys the flexibility that comes from the lack of a prescribed mesh. It is suitable for problems such as free surface flow or flow around complex, moving and/or intermeshed geometries and is applicable to both dilute and dense particle flows. Second, a comprehensive validation procedure for fluid–particle simulations is presented and applied here to the SPH–DEM method, using simulations of single and multiple particle sedimentation in a 3D fluid column and comparison with analytical models. Millimetre-sized particles are used along with three different test fluids: air, water and a water–glycerol solution. The velocity evolution for a single particle compares well (less than 1% error) with the analytical solution as long as the fluid resolution is coarser than two times the particle diameter. Two more complex multiple particle sedimentation problems (sedimentation of a homogeneous porous block and an inhomogeneous Rayleigh Taylor Instability) are also reproduced well for porosities 0.6⩽∊⩽1.0, although care should be taken in the presence of high porosity gradients. Overall the SPH–DEM method successfully reproduces quantitatively the expected behaviour in these test cases, and promises to be a flexible and accurate tool for other, realistic fluid–particle system simulations (for which other problem-relevant test cases have to be added for validation).

Experiments and discrete element simulation of the dosing of cohesive powders in a simplified geometry

We perform experiments and discrete element simulations on the dosing of cohesive granular materials in a simplified geometry. The setup is a simplified canister box where the powder is dosed out of the box through the action of a constant-pitch screw feeder connected to a motor. A dose consists of a rotation step followed by a period of rest before the next dosage. From the experiments, we report on the operational performance of the dosing process through a variation of dosage time, coil pitch and initial powder mass. We find that the dosed mass shows an increasing linear dependence on the dosage time and rotation speed. In contrast, the mass output from the canister is not directly proportional to an increase/decrease in the number coils. By calibrating the interparticle friction and cohesion, we show that DEM simulation can quantitatively reproduce the experimental findings for smaller masses but also overestimate arching and blockage. With appropriate homogenization tools, further insights into microstructure and macroscopic fields can be obtained. This work shows that particle scaling and the adaptation of particle properties is a viable approach to overcome the untreatable number of particles inherent in experiments with fine, cohesive powders and opens the gateway to simulating their flow in more complex geometries.

Wicking in a Powder

We investigate the wicking in granular media by considering layers of grains at the surface of a liquid and discuss the critical contact angle below which spontaneous impregnation takes place. This angle is found to be on the order of 55° for monodisperse layers, significantly smaller than 90°, the threshold value for penetrating assemblies of tubes. Owing to geometry, impregnating grains is more demanding than impregnating tubes. We also consider the additional effects of polydispersity and pressure on this wetting transition and discuss the corresponding shift observed for the critical contact angle.

Dynamic wetting on a thin film of soluble polymer: effects of nonlinearities in the sorption isotherm.

The wetting dynamics of a solvent on a soluble substrate interestingly results from the rates of the solvent transfers into the substrate. When a supported film of a hydrosoluble polymer with thickness e is wet by a spreading droplet of water with instantaneous velocity U, the contact angle is measured to be inversely proportionate to the product of thickness and velocity, eU, over two decades. As for many hydrosoluble polymers, the polymer we used (a polysaccharide) has a strongly nonlinear sorption isotherm φ(a(w)), where φ is the volume fraction of water in the polymer and aw is the activity of water. For the first time, this nonlinearity is accounted for in the dynamics of water uptake by the substrate. Indeed, by measuring the water content in the polymer around the droplet φ at distances as small as 5 μm, we find that the hydration profile exhibits (i) a strongly distorted shape that results directly from the nonlinearities of the sorption isotherm and (ii) a cutoff length ξ below which the water content in the substrate varies very slowly. The nonlinearities in the sorption isotherm and the hydration at small distances from the line were not accounted for by Tay et al., Soft Matter 2011, 7, 6953. Here, we develop a comprehensive description of the hydration of the substrate ahead of the contact line that encompasses the two water transfers at stake: (i) the evaporation-condensation process by which water transfers into the substrate through the atmosphere by the condensation of the vapor phase, which is fed by the evaporation from the droplet itself, and (ii) the diffusion of liquid water along the polymer film. We find that the eU rescaling of the contact angle arises from the evaporation-condensation process at small distances. We demonstrate why it is not modified by the second process.

Tuning the bulk properties of bidisperse granular mixtures by small amount of fines

We study the bulk properties of isotropic bidisperse granular mixtures using discrete element simulations. The focus is on the influence of the size (radius) ratio of the two constituents and volume fraction on the mixture properties. We show that the effective bulk modulus of a dense granular (base) assembly can be enhanced by up to 20% by substituting as little as 5% of its volume with smaller sized particles. Particles of similar sizes barely affect the macroscopic properties of the mixture. On the other extreme, when a huge number of fine particles are included, most of them lie in the voids of the base material, acting as rattlers, leading to an overall weakening effect. In between the limits, an optimum size ratio that maximizes the bulk modulus of the mixture is found. For loose systems, the bulk modulus decreases monotonically with addition of fines regardless of the size ratio. Finally, we relate the mixture properties to the ‘typical’ pore size in a disordered structure as induced by the combined effect of operating volume fraction (consolidation) and size ratio.

Powder wettability at a static air–water interface

The reconstitution of a beverage from a dehydrated powder involves several physical mechanisms that determine the practical difficulty to obtain a homogeneous drink in a convenient way and within an acceptable time for the preparation of a beverage. When pouring powder onto static water, the first hurdle to overcome is the air-water interface. We propose a model to predict the percentage of powder crossing the interface in 45 s, namely the duration relevant for this application. We highlight theoretically the determinant role of the contact angle and of the particle size distribution. We validate experimentally the model for single spheres and use it to predict the wettability performance of commercial food powders for different contact angles and particles sizes. A good agreement is obtained when comparing the predictions and the wettability of the tested powders.

Recent advances in the simulation of particle-laden flows

A substantial number of algorithms exists for the simulation of moving particles suspended in fluids. However, finding the best method to address a particular physical problem is often highly non-trivial and depends on the properties of the particles and the involved fluid(s) together. In this report, we provide a short overview on a number of existing simulation methods and provide two state of the art examples in more detail. In both cases, the particles are described using a Discrete Element Method (DEM). The DEM solver is usually coupled to a fluid-solver, which can be classified as grid-based or mesh-free (one example for each is given). Fluid solvers feature different resolutions relative to the particle size and separation. First, a multicomponent lattice Boltzmann algorithm (mesh-based and with rather fine resolution) is presented to study the behavior of particle stabilized fluid interfaces and second, a Smoothed Particle Hydrodynamics implementation (mesh-free, meso-scale resolution, similar to the particle size) is introduced to highlight a new player in the field, which is expected to be particularly suited for flows including free surfaces.

In vivo observations and in vitro experiments on the oral phase of swallowing of Newtonian and shear-thinning liquids

In this study, an in vitro device that mimics the oral phase of swallowing is calibrated using in vivo measurements. The oral flow behavior of different Newtonian and non-Newtonian solutions is then investigated in vitro, revealing that shear-thinning thickeners used in the treatment of dysphagia behave very similar to low-viscosity Newtonian liquids during active swallowing, but provide better control of the bolus before the swallow is initiated. A theoretical model is used to interpret the experimental results and enables the identification of two dynamical regimes for the flow of the bolus: first, an inertial regime of constant acceleration dependent on the applied force and system inertia, possibly followed by a viscous regime in which the viscosity governs the constant velocity of the bolus. This mechanistic understanding provides a plausible explanation for similarities and differences in swallowing performance of shear-thinning and Newtonian liquids. Finally, the physiological implications of the model and experimental results are discussed. In vitro and theoretical results suggest that individuals with poor tongue strength are more sensitive to overly thickened boluses. The model also suggests that while the effects of system inertia are significant, the density of the bolus itself plays a negligible role in its dynamics. This is confirmed by experiments on a high density contrast agent used for videofluoroscopy, revealing that rheologically matched contrast agents and thickener solutions flow very similarly. In vitro experiments and theoretical insights can help designing novel thickener formulations before clinical evaluations.

SPH-DEM simulations of grain dispersion by liquid injection

We study the dispersion of an initially packed, static granular bed by the injection of a liquid jet. This is a relevant system for many industrial applications, including paint dispersion or food powder dissolution. Both decompaction and dispersion of the powder are not fully understood, leading to inefficiencies in these processes. Here we consider a model problem where the liquid jet is injected below a granular bed contained in a cylindrical cell. Two different initial conditions are considered: a two-phase case where the bed is initially fully immersed in the liquid and a three-phase case where the bed and cell are completely dry preceding the injection of the liquid. The focus of this contribution is the simulation of these model problems using a two-way coupled SPH-DEM granularliquid method [M. Robinson, M. Ramaioli, and S. Luding, submitted (2013) and http://arxiv.org/abs/1301.0752 (2013)]. This is a purely particle-based method without any prescribed mesh, well suited for this and other problems involving a free (liquidgas) surface and a partly immersed particle phase. Our simulations show the effect of process parameters such as injection flow rate and injection diameter on the dispersion pattern, namely whether the granular bed is impregnated bottom-up or a jet is formed and compare well with experiments.

Grain Sedimentation with SPH-DEM and its Validation

Our mesoscale simulation method [M. Robinson, S. Luding, and M. Ramaioli, submitted (2013)] for multiphase fluid-particle flows couples Smoothed Particle Hydrodynamics (SPH) and the Discrete Element Method (DEM) and enjoys the flexibility of meshless methods, such as being capable to handling free surface flows or flow around complex and/or moving geometries. We use this method to simulate three different sedimentation test cases and compare the results to existing analytical solutions. The grain velocity in Single Particle Sedimentation compares well (< 2% error) with the analytical solution as long as the fluid resolution is coarser than two times the particle diameter. The multiple particle sedimentation problem and Rayleigh Taylor Instability (RTI) also perform well against the theory, but it was found that the method is susceptible to fluid velocity fluctuations in the presence of high porosity gradients. These fluctuations can be damped by the addition of a dissipation term, which has no effect on the terminal velocity but can lead to slower growth rates for the RTI.