CFD-DEM Fluidized Bed Drying Study Using a Coarse-Graining Technique.

Fluidized beds are commonly applied to industrial drying applications. Modeling using the computational fluid dynamics-discrete element method (CFD-DEM) can be employed to increase the fundamental understanding of solids drying. A large drawback of CFD-DEM is the computational requirements, leading to a limitation regarding the system size. Coarse-grained CFD-DEM is an approach to reduce computational costs, allowing one to simulate larger fluidized beds. In this article, coarse-graining CFD-DEM scaling laws are used for fluidized bed solids drying. Three superficial gas velocities are investigated. The particle temperature and density are accurately described. Besides, the Sherwood number is well captured by the coarse-graining simulations.

Effect of microchannel structure and fluid properties on non-inertial particle migration.

In this work, we investigate the influence of channel structure and fluid rheology on non-inertial migration of non-Brownian polystyrene beads. Particle migration in this regime can be found in biomedical, chemical, environmental and geological applications. However, the effect of fluid rheology on particle migration in porous media remains to be clearly understood. Here, we isolate the effects of elasticity and shear thinning by comparing a Newtonian fluid, a purely elastic (Boger) fluid, and a shear-thinning elastic fluid. To mimic the complexity of geometries in real-world application, a random porous structure is created through a disordered arrangement of cylindrical pillars in the microchannel. Experiments are repeated in an empty channel and in channels with an ordered arrangement of pillars, and the similarities and differences in the observed particle focusing are analyzed. It is found that elasticity drives the particles away from the channel walls in an empty microchannel. Notably, particle focusing is unaffected by curved streamlines in an ordered porous microchannel and particles stay away from pillars in elastic fluids. Shear-thinning is found to reduce the effect of focusing and a broader region of particle concentration is observed. It is also noteworthy that the rheological characteristics of the fluid are not important for the particle distribution in a randomly arranged pillared microchannel and particles have a uniform distribution for all suspending fluids. Moreover, discussion on the current discrepancy in the literature about the equilibrium positions of the particles in a channel is extended by analyzing the results obtained in the current experiments.

Nonspherical particles in a pseudo-2D fluidized bed: Experimental study.

Fluidization is widely used in industries and has been extensively studied, both experimentally and theoretically, in the past. However, most of these studies focus on spherical particles while in practice granules are rarely spherical. Particle shape can have a significant effect on fluidization characteristics. It is therefore important to study the effect of particle shape on fluidization behavior in detail. In this study, experiments in pseudo-2D fluidized beds are used to characterize the fluidization of spherocylindrical (rod-like) Geldart D particles of aspect ratio 4. Pressure drop and optical measurement methods (Digital Image Analysis, Particle Image Velocimetry, Particle Tracking Velocimetry) are employed to measure bed height, particle orientation, particle circulation, stacking, and coordination number. The commonly used correlations to determine the pressure drop across a bed of nonspherical particles are compared to experiments. Experimental observations and measurements have shown that rod-like particles are prone to interlocking and channeling behavior. Well above the minimum fluidization velocity, vigorous bubbling fluidization is observed, with groups of interlocked particles moving upwards, breaking up, being thrown high in the freeboard region and slowly raining down as dispersed phase. At high flowrates, a circulation pattern develops with particles moving up through the center and down at the walls. Particles tend to orient themselves along the flow direction.

Viscoelastic effects on residual oil distribution in flows through pillared microchannels.

HYPOTHESIS: Multiphase flow through porous media is important in a number of industrial, natural and biological processes. One application is enhanced oil recovery (EOR), where a resident oil phase is displaced by a Newtonian or polymeric fluid. In EOR, the two-phase immiscible displacement through heterogonous porous media is usually governed by competing viscous and capillary forces, expressed through a Capillary number Ca, and viscosity ratio of the displacing and displaced fluid. However, when viscoelastic displacement fluids are used, elastic forces in the displacement fluid also become significant. It is hypothesized that elastic instabilities are responsible for enhanced oil recovery through an elastic microsweep mechanism. EXPERIMENTS: In this work, we use a simplified geometry in the form of a pillared microchannel. We analyze the trapped residual oil size distribution after displacement by a Newtonian fluid, a nearly inelastic shear thinning fluid, and viscoelastic polymers and surfactant solutions. FINDINGS: We find that viscoelastic polymers and surfactant solutions can displace more oil compared to Newtonian fluids and nearly inelastic shear thinning polymers at similar Ca numbers. Beyond a critical Ca number, the size of residual oil blobs decreases significantly for viscoelastic fluids. This critical Ca number directly corresponds to flow rates where elastic instabilities occur in single phase flow, suggesting a close link between enhancement of oil recovery and appearance of elastic instabilities.

Magnetic particle tracking for nonspherical particles in a cylindrical fluidized bed.

In granular flow operations, often particles are nonspherical. This has inspired a vast amount of research in understanding the behavior of these particles. Various models are being developed to study the hydrodynamics involving nonspherical particles. Experiments however are often limited to obtain data on the translational motion only. This paper focusses on the unique capability of Magnetic Particle Tracking to track the orientation of a marker in a full 3-D cylindrical fluidized bed. Stainless steel particles with the same volume and different aspect ratios are fluidized at a range of superficial gas velocities. Spherical and rod-like particles show distinctly different fluidization behavior. Also, the distribution of angles for rod-like particles changes with position in the fluidized bed as well as with the superficial velocity. Magnetic Particle Tracking shows its unique capability to study both spatial distribution and orientation of the particles allowing more in-depth validation of Discrete Particle Models.

Viscoelastic flow past mono- and bidisperse random arrays of cylinders: flow resistance, topology and normal stress distribution.

We investigate creeping viscoelastic fluid flow through two-dimensional porous media consisting of random arrangements of monodisperse and bidisperse cylinders, using our finite volume-immersed boundary method introduced in S. De, et al., J. Non-Newtonian Fluid Mech., 2016, 232, 67-76. The viscoelastic fluid is modeled with a FENE-P model. The simulations show an increased flow resistance with increase in flow rate, even though the bulk response of the fluid to shear flow is shear thinning. We show that if the square root of the permeability is chosen as the characteristic length scale in the determination of the dimensionless Deborah number (De), then all flow resistance curves collapse to a single master curve, irrespective of the pore geometry. Our study reveals how viscoelastic stresses and flow topologies (rotation, shear and extension) are distributed through the porous media, and how they evolve with increasing De. We correlate the local viscoelastic first normal stress differences with the local flow topology and show that the largest normal stress differences are located in shear flow dominated regions and not in extensional flow dominated regions at higher viscoelasticity. The study shows that normal stress differences in shear flow regions may play a crucial role in the increase of flow resistance for viscoelastic flow through such porous media.

Particle-based modeling of heterogeneous chemical kinetics including mass transfer.

Connecting the macroscopic world of continuous fields to the microscopic world of discrete molecular events is important for understanding several phenomena occurring at physical boundaries of systems. An important example is heterogeneous catalysis, where reactions take place at active surfaces, but the effective reaction rates are determined by transport limitations in the bulk fluid and reaction limitations on the catalyst surface. In this work we study the macro-micro connection in a model heterogeneous catalytic reactor by means of stochastic rotation dynamics. The model is able to resolve the convective and diffusive interplay between participating species, while including adsorption, desorption, and reaction processes on the catalytic surface. Here we apply the simulation methodology to a simple straight microchannel with a catalytic strip. Dimensionless Damkohler numbers are used to comment on the spatial concentration profiles of reactants and products near the catalyst strip and in the bulk. We end the discussion with an outlook on more complicated geometries and increasingly complex reactions.

Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) Study of Mass-Transfer Mechanisms in Riser Flow.

We report a computational fluid dynamics-discrete element method (CFD-DEM) simulation study on the interplay between mass transfer and a heterogeneous catalyzed chemical reaction in cocurrent gas-particle flows as encountered in risers. Slip velocity, axial gas dispersion, gas bypassing, and particle mixing phenomena have been evaluated under riser flow conditions to study the complex system behavior in detail. The most important factors are found to be directly related to particle cluster formation. Low air-to-solids flux ratios lead to more heterogeneous systems, where the cluster formation is more pronounced and mass transfer more influenced. Falling clusters can be partially circumvented by the gas phase, which therefore does not fully interact with the cluster particles, leading to poor gas-solid contact efficiencies. Cluster gas-solid contact efficiencies are quantified at several gas superficial velocities, reaction rates, and dilution factors in order to gain more insight regarding the influence of clustering phenomena on the performance of riser reactors.

Numerical investigation of the vertical plunging force of a spherical intruder into a prefluidized granular bed.

The plunging of a large intruder sphere into a prefluidized granular bed with various constant velocities and various sphere diameters is investigated using a state-of-the-art hybrid discrete particle and immersed boundary method, in which both the gas-induced drag force and the contact force exerted on the intruder can be investigated separately. We investigate low velocities, where velocity dependent effects first begin to appear. The results show a concave-to-convex dependence of the plunging force as a function of intruder depth. In the concave region the force fits to a power law with an exponent around 1.3, which is in good agreement with existing experimental observations. Our simulation results further show that the force exerted on the frontal hemisphere of the intruder is dominant. At larger intruder velocities, friction with the granular medium causes a velocity-dependent drag force. As long as the granular particles have not yet closed the gap behind the intruder, this drag force is independent of the actual intruder depth. In this regime, the drag force experienced by intruders of different diameter moving at different velocities all fall onto a single master curve if plotted against the Reynolds number, using a single value for the effective viscosity of the granular medium. This master curve corresponds well to the Schiller-Naumann correlation for the drag force between a sphere and a Newtonian fluid. After the gap behind the intruder has closed, the drag force increases not only with velocity but also with depth. We attribute this to the effect of increasing hydrostatic particle pressure in the granular medium, leading to an increase in effective viscosity.

Simulation of density segregation in vibrated beds.

We have investigated by numerical simulation the density segregation of fine equal-sized bronze and glass particles subject to vertical vibrations. The model was found to be capable of predicting the two main segregation forms ("bronze on top" and "sandwich") in roughly the same regions of the phase diagram as was found experimentally by Burtally We investigated the effects of pressure air forcing, friction and restitution of kinetic energy in collisions, and box size on the segregation behavior. We find that next to the interstitial air friction also has a large influence on the formation of the sandwich structure.

Simulation study of air-induced segregation of equal-sized bronze and glass particles.

We have investigated the effect of air on the density segregation of fine, equal-sized bronze and glass particles under vertical sinusoidal vibrations by numerical simulation, using a hybrid granular dynamics-computational fluid dynamics model. We find both the bronze-on-top and the sandwich configurations as observed in previous experiments. The simulations allow us to give a detailed explanation of the phenomena. The key observation is the precise timing of the decompaction of the bed with the oscillating gas flow, which will cause the bronze to move to the top, due to the difference in acceleration from the air drag. For higher frequencies, the different inertia in collisions will cause the bronze clusters on top to sink through the bed, leading to a sandwich formation.