Multiphase flows

The motion of dispersed objects (particles, droplets or bubbles) with a diameter larger than 10 micrometers is dominated by gravity, drag forces with the surrounding continuous medium (often air or water), and collisional contacts between the objects and with walls. Because such macroscopic objects collide inelastically, large-scale heterogeneous structures emerge. For many engineering applications it is important to be able to predict these dispersed flows through pipes, chutes, reactors and other vessels.

We are developing simulation models for macroscopic (non-Brownian) dispersed flows, where we use a multi-scale approach to cover the entire range of length scales from the particle level to the level of industrial size reactors. On the detailed level, we take into account the fundamental processes of momentum, heat and mass transfer, whereas on the larger levels we use effective (coarse-grained) interactions and correlations obtained on the detailed levels, often borrowing methods from statistical mechanics and soft matter science. Because a model is only as good as its assumptions (garbage in: garbage out), an important part of our work is the experimental validation of our models, through advanced non-invasive techniques such as Particle Image Velocimetry.

Granular flow in rotating equipment

mixing of granular particles in a rotating drum In exceptional cases, when granular particles are relatively large and the particle velocities are relatively low, the effect of the surrounding gas can be neglected. In one application, we have studied the mixing of granular materials in a rotating drum as a function of filling height and baffle size. We found that only relatively large, straight baffles perpendicular to the drum wall (67% of the drum radius) increase the mixing performance of the rotating drum.

experimental rotating table setup experimental rotating table setup More recently, we have studied the flow of granular particles down an inclined and rotating chute. The goal of this work is to study the influence of chute rotation on the flow behaviour and segregation in granular flows, and to obtain high quality experimental data with which discrete particle model simulations can be validated. The figures to the right show the experimental setup on a rotating table, and a snapshot (top view) of our simulations. We measure the surface particle velocity field through particle image velocimetry and particel tracking velocimetry. We find that, compared to a non-rotating chute, the surface flow velocity in a rotating chute is decreased initially due to compaction against the side walls (Coriolis effect) but ultimately increased due to the increasing influence of centrifugal forces.

Fluidized beds

In a fluidized bed, an uprising gas is flowing so fast that it causes the solid particles to "float". As a result, the particle emulsion has many properties of a liquid, including the formation of flow vortices and bubbles. Fluidized beds are used frequently in chemical reactors because they allow an intense contact between the gas (which may contain reactants) and the solid particles (which may be a growing product or contain catalysts).

contour plots of the solids volume fraction in a fluidized bed In one line of research, we are focussing on developing large scale computational models for gas-solid fluidized beds. One direction is to use stochastic Lagrangian coarse-grained multi-particle collision techniques borrowed from soft matter simulations, see the section on colloids. Another direction is to use and improve the so-called Two-Fluid Model, in which the gas and solid phases are treated as two interpenetrating continua. The latter approach requires accurate expressions for the rheology of the solid phase, for which we are using the kinetic theory of granular flow. The picture above shows contour plots of the solids volume fraction at three successive times in a cylindrical fluidized bed, obtained from these TFM simulations.

an intruder impacting onto a pre-fluidized bed of granular particles In another line of resarch, we are focussing on accurate simulation methods to measure the viscosity of the granular particle emulsion and to model the influence of intruders (which may also be baffles or stirrers). The viscoscity is obtained both through Green-Kubo relations and more directly through measurement of the stress in Couette flow and measurement of the drag on a dragged sphere. We have developed a hybrid method in which granular particles are coupled with the surrounding gas through effective drag correlations, whereas large intruders are coupled more directly through immersed boundaries. The influence of an impacting intruder is first studied for a pre-fluidized bed, in which the gas flow has slowly been turned off again. The picture to the right shows the effect of a large spherical intruder impacting into a pre-fluidized bed of much smaller granular particles.


Dry separation of materials through triboelectric charging

Our goal is to develop new separation technologies for dry fractionation of food products. The dry fractionation is based on electrostatic charging of powders in a gas stream when they collide with the walls of a charging tube, and subsequent separation by an external electric field (see figure above). Our aim is to accurately model the dynamics and acquired charge of the particles in the charging tube, where triboelectric charging takes place when a particle collides with the metal walls. The challenge of this project is that we need to combine triboelectric charging of the particles (about which little is known) with electrostatic interactions and gas flow. The first results are encouraging.


droplets emerging from a spray atomizer Spray drying is used to convert liquid feed materials into a dry powder form, e.g. to produce milk powder. The liquid feed is sprayed (atomized) through high pressure nozzles in a spray chamber (see picture on the left) and the resulting droplets are mixed with hot gas to evaporate the liquid content.

The main objective of this work is to develop a reliable simulation tool that can predict the droplet and particle flow, agglomeration, and size distribution for a section of a large-scale spray dryer. To be able to handle many billions of droplets and particles, we have adapted a direct simulation monte carlo method, in which the droplets and particles collide stochastically according to predictions from kinetic theory. Full two-way coupling between the droplets/particles and the gas has been implemented. Collisions between droplets and/or particles lead to different outcomes, depending on dimensionless numbers (Weber number, Ohnesorge number and impact parameter). The figure below shows the simulation results for a spray of pure droplets.

the outcome of a droplet-droplet collision depends on Weber number and impact parameter
  • S.K. Pawar, J.T. Padding, N.G. Deen, J.A.M. Kuipers, A. Jongsma and F. Innings, Eulerian-Lagrangian modelling with stochastic approach for droplets-particles collision, proceedings of the CSIRO Ninth International Conference on Computational Fluid Dynamics in the Minerals and Process Industries, Melbourne (2012).
  • S.K. Pawar, J.T. Padding, N.G. Deen, A. Jongsma, F. Innings and J.A.M. Kuipers, Lagrangian modelling of dilute granular flow-modified stochastic DSMC versus deterministic DPM, Chem. Eng. Sci. 105, 132 (2014).
  • Sandip Pawar, Ruud Abrahams, Niels Deen, Johan Padding, Gert-Jan van der Gulik, Alfred Jongsma, Fredrik Innings and Hans Kuipers, An Experimental Study of Dynamic Jet Behavior in a Scaled Cold Flow Spray Dryer Model Using PIV, Can. J. Chem. Eng. 92, 2013 (2014).
  • S. Pawar, J.T. Padding, N.G. Deen, A. Jongsma, F. Innings and J.A.M. Kuipers, Agglomeration study in the inlet section of a large scale spray dryer using stochastic euler-lagrange modelling, proceedings of the SINTEF 10th International Conference on Computational Fluid Dynamics in the Minerals and Process Industries, Trondheim (2014).
  • S.K. Pawar, J.T. Padding, N.G. Deen, A. Jongsma, F. Innings and J.A.M. Kuipers, Drying kinetics of droplets containing solids using reaction engineering approach, proceedings of the 19th International Drying Symposium, Lyon, France (2014).
  • S.K. Pawar, J.T. Padding, N.G. Deen, A. Jongsma, F. Innings and J.A.M. Kuipers, Numerical and Experimental Investigation of Induced Flow and Droplet-Droplet Interactions in a Liquid Spray, Chem. Eng. Sci. 138, 17 (2015).
  • S.K. Pawar, F. Henrikson, G. Finotello, J.T. Padding, N.G. Deen, A. Jongsma, F. Innings and J.A.M. Kuipers, An experimental study of droplet-particle collisions, Powder Technology 300, 157 (2016).
  • H. Patel, S. Das, J.A.M. Kuipers, J.T. Padding and E.A.J.F. Peters, A coupled VOF-IBM method for simulating 3D multiphase flows with contact line dynamics in complex geometries, Chem. Eng. Sci. 166, 28 (2017).
  • G. Finotello, J.T. Padding, N.G. Deen, A. Jongsma, F. Innings and J.A.M. Kuipers, Effect of viscosity on droplet coalescence and breakup, to be submitted (2016).
  • S.K. Pawar, J.T. Padding, N.G. Deen and J.A.M. Kuipers, A Novel Hybrid Computational Model Combining DSMC and a Kinetic Particle Model for Simulation of Spray Drying of Milk, to be submitted (2017).

Bubbly flows

evolution of bubbles in a bubble column Bubble column reactors are widely used as gas-liquid contactors in the chemical and energy industries. The bubble dynamics inside the column dictate the liquid phase flow. Modeling of such bubble columns is very challenging and typically empirical correlations are used even today.

The main objective of this work is to develop a stochastic Euler-Lagrange method that can predict the hydrodynamics of dense bubbly flows. For this we are adapting the direct simulation monte carlo approach to handle bubble collisions stochastically instead of deterministically, enhancing the speed of bubble collision calculations by a factor of 10 to 100 relative to more detailed discrete bubble model simulations, while still maintaining the same level of accuracy.

  • S. Kamath, J.T. Padding, K.A. Buist and J.A.M. Kuipers Stochastic DSMC model for large scale dense bubbly flows, proceedings of the SINTEF 12th International Conference on Computational Fluid Dynamics in the Minerals and Process Industries, Trondheim (2017).
  • S.S. Kamath, J.T. Padding, K.A. Buist and J.A.M. Kuipers, DSMC for large scale dense bubbly flows, to be submitted (2017).


The research on granular flow was performed by my former PhD student Sushil Shirsath. The research on granular viscosity and intruder impact is performed with my PhD student Yupeng Xu. Within these projects we collaborate intensely with the Turbulence and Vortex Dynamics group (Prof. Herman Clercx) at Eindhoven University of Technology, with Prof. Detlef Lohse, Prof. Devaraj van der Meer and Dr. Martin van der Hoef of the Physics of Fluids group and with Prof. Stefan Luding and Dr. Anthony Thornton of the Multiscale Mechanics group, both at the University of Twente, and Prof. Onno Bokhove at University of Leeds. The research on fluidized beds is performed by my current and former PhD students Luuk Seelen and Vikrant Verma, and is financially supported by an ERC-advanced grant of Prof. Hans Kuipers. The research on triboelectric charging is performed with my PhD student Martin Korevaar, in collaboration with Dr. Maarten Schutyser of Wageningen University, and is financially supported by an STW grant. The project on spray-drying was performed by our previous PhD student Sandip Pawar, continued by our current PhD student Giulia Finotello, and is financially supported by Tetrapak Heerenveen; we are indebted to intense collaborations with Dr. Alfred Jongsma and Dr. Fredrik Innings.