Stochastic Lagrangian model for hydrodynamic acceleration of inertial particles in gas–solid suspensions

2016 ◽  
Vol 788 ◽  
pp. 695-729 ◽  
Author(s):  
Sudheer Tenneti ◽  
Mohammad Mehrabadi ◽  
Shankar Subramaniam
Keyword(s):  
Particle Acceleration ◽  
Particle Velocity ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Density Ratio ◽  
Granular Temperature ◽  
Solid Flow ◽  
Acceleration Model ◽  
The Mean ◽  
Solid Suspensions

The acceleration of an inertial particle in a gas–solid flow arises from the particle’s interaction with the gas and from interparticle interactions such as collisions. Analytical treatments to derive a particle acceleration model are difficult outside the Stokes flow regime, but for moderate Reynolds numbers (based on the mean slip velocity between gas and particles) particle-resolved direct numerical simulation (PR-DNS) is a viable tool for model development. In this study, PR-DNS of freely-evolving gas–solid suspensions are performed using the particle-resolved uncontaminated-fluid reconcilable immersed-boundary method (PUReIBM) that has been extensively validated in previous studies. Analysis of the particle velocity variance (granular temperature) equation in statistically homogeneous gas–solid flow shows that a straightforward extension of a class of mean particle acceleration models (drag laws) to their corresponding instantaneous versions, by replacing the mean particle velocity with the instantaneous particle velocity, predicts a granular temperature that decays to zero, which is at variance with the steady particle granular temperature that is obtained from PR-DNS. Fluctuations in particle velocity and particle acceleration (and their correlation) are important because the particle acceleration–velocity covariance governs the evolution of the particle velocity variance (characterized by the particle granular temperature), which plays an important role in the prediction of the core annular structure in riser flows. The acceleration–velocity covariance arising from hydrodynamic forces can be decomposed into source and dissipation terms that appear in the granular temperature evolution equation, and these have already been quantified in the Stokes flow regime using a combination of kinetic theory closure and multipole expansion simulations. From PR-DNS data we show that the fluctuations in the particle acceleration that are aligned with fluctuations in the particle velocity give rise to a source term in the granular temperature evolution equation. This approach is used to quantify the hydrodynamic source and dissipation terms of granular temperature from PR-DNS results for freely-evolving gas–solid suspensions that are performed over a wide range of solid volume fraction ($0.1\leqslant {\it\phi}\leqslant 0.4$), Reynolds number based on the slip velocity between the solid and the fluid phase ($10\leqslant \mathit{Re}_{m}\leqslant 100$) and solid-to-fluid density ratio ($100\leqslant {\it\rho}_{p}/{\it\rho}_{f}\leqslant 2000$). The straightforward extension of drag law models does not give rise to any source in the granular temperature due to hydrodynamic effects. This motivates the development of better Lagrangian particle acceleration models that can be used in Lagrangian–Eulerian formulations of gas–solid flow. It is found that a Langevin equation for the increment in the particle velocity reproduces PR-DNS results for the stationary particle velocity autocorrelation in freely-evolving suspensions. Based on the data obtained from the simulations, the functional dependence of the Langevin model coefficients on solid volume fraction, Reynolds number and solid-to-fluid density ratio is obtained. This new Lagrangian particle acceleration model reproduces the correct steady granular temperature and can also be adapted to gas–solid flow computations using Eulerian moment equations.

Analysis of Dispersion of Small Spherical Particles in a Random Velocity Field

1990 ◽  
Vol 112 (1) ◽  
pp. 114-120 ◽  
Author(s):  
H. Ounis ◽  
G. Ahmadi
Keyword(s):  
Particle Velocity ◽  
Digital Simulation ◽  
Density Ratio ◽  
Spherical Particles ◽  
Theoretical Predictions ◽  
The Mean ◽  
Analysis Of Dispersion ◽  
Gradient Effects ◽  
Analytical Expressions ◽  
Good Agreement

The equation of motion of a small spherical rigid particle in a turbulent flow field, including the Stokes drag, the Basset force, and the virtual mass effects, is considered. For an isotropic field, the lift force and the velocity gradient effects are neglected. Using the spectral method, responses of the resulting constant coefficient stochastic integrao-differential equation are studied. Analytical expressions relating the Lagrangian energy spectra of particle velocity to that of the fluid are developed and the results are used to evaluate various response statistics. Variations of the mean-square particle velocity and particle diffusivity with size, density ratio and response time are studied. The theoretical predictions are compared with the digital simulation results and the available data and good agreement is observed.

scholarly journals Interface-resolved simulations of particle suspensions in Newtonian, shear thinning and shear thickening carrier fluids

2018 ◽  
Vol 852 ◽  
pp. 329-357 ◽  
Author(s):  
Dhiya Alghalibi ◽  
Iman Lashgari ◽  
Luca Brandt ◽  
Sarah Hormozi
Keyword(s):  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Shear Thickening ◽  
Fluid Phase ◽  
Local Strain ◽  
Reynolds Numbers ◽  
Mean Value ◽  
Strain Rates ◽  
Local Shear ◽  
The Mean

We present a numerical study of non-colloidal spherical and rigid particles suspended in Newtonian, shear thinning and shear thickening fluids employing an immersed boundary method. We consider a linear Couette configuration to explore a wide range of solid volume fractions ($0.1\leqslant \unicode[STIX]{x1D6F7}\leqslant 0.4$) and particle Reynolds numbers ($0.1\leqslant Re_{p}\leqslant 10$). We report the distribution of solid and fluid phase velocity and solid volume fraction and show that close to the boundaries inertial effects result in a significant slip velocity between the solid and fluid phase. The local solid volume fraction profiles indicate particle layering close to the walls, which increases with the nominal $\unicode[STIX]{x1D6F7}$. This feature is associated with the confinement effects. We calculate the probability density function of local strain rates and compare the latter’s mean value with the values estimated from the homogenisation theory of Chateau et al. (J. Rheol., vol. 52, 2008, pp. 489–506), indicating a reasonable agreement in the Stokesian regime. Both the mean value and standard deviation of the local strain rates increase primarily with the solid volume fraction and secondarily with the $Re_{p}$. The wide spectrum of the local shear rate and its dependency on $\unicode[STIX]{x1D6F7}$ and $Re_{p}$ point to the deficiencies of the mean value of the local shear rates in estimating the rheology of these non-colloidal complex suspensions. Finally, we show that in the presence of inertia, the effective viscosity of these non-colloidal suspensions deviates from that of Stokesian suspensions. We discuss how inertia affects the microstructure and provide a scaling argument to give a closure for the suspension shear stress for both Newtonian and power-law suspending fluids. The stress closure is valid for moderate particle Reynolds numbers, $O(Re_{p})\sim 10$.

Computational Study of Slurry Flow in Pipes to Determine Profiles Sensed in Near Infrared Slurry Measurements

2011 ◽  
Author(s):  
V. Pasangulapati ◽  
N. R. Kesana ◽  
G. Sharma ◽  
F. W. Chambers ◽  
M. E. McNally ◽  
...  
Keyword(s):  
Near Infrared ◽  
Volume Fraction ◽  
Solid Phase ◽  
Computational Study ◽  
Solid Volume Fraction ◽  
Density Ratio ◽  
Optical Measurements ◽  
Concentration Profiles ◽  
Horizontal Pipe ◽  
Volume Fractions

It is desired to perform accurate Near Infrared sensor measurements of slurries flowing in pipes leaving large batch reactors. A concern with these measurements is the degree to which the slurry sensed is representative of the material in the reactor and flowing through the pipe. Computational Fluid Dynamics (CFD) has been applied to the flow in the pipe to determine the flow fields and the concentration profiles seen by the sensors. The slurry was comprised of a xylene liquid phase and an ADP (2-amino-4, 6-dimethylpyrimidine) solid phase with a density ratio of 1.7. Computations were performed for a horizontal pipe with diameter 50.8 mm, length 2.032 m, and 1.76 m/s and 3.26 m/s mixture velocities. The corresponding pipe Reynolds numbers were 1.19E+05 and 2.21E+05. The flow through a slotted cylindrical probe inserted radially in the pipe also was considered. Spherical slurry particles with diameters from 10 μm to 1000 μm were considered with solid volume fractions of 12%, 24%, and 35%. Computations were performed with ANSYS FLUENT 12 software using the Realizable k-ε turbulence model and the enhanced wall treatment function. Comparisons of computed vertical profiles of solid volume fraction to results in the literature showed good agreement. Symmetric, nearly flat solid volume fraction profiles were observed for 38 μm particles for all three initial solid volume fractions. Asymmetric solid volume fraction profiles with greater values toward the bottom were observed for the larger particles. Changes in the profiles of turbulent kinetic energy also were observed. These changes are important for optical measurements which depend upon the mean concentration profiles as well as the turbulent motion of the slurry particles.

Lateral dispersion in random cylinder arrays at high Reynolds number

2008 ◽  
Vol 600 ◽  
pp. 339-371 ◽  
Author(s):  
YUKIE TANINO ◽  
HEIDI M. NEPF
Keyword(s):  
Turbulence Intensity ◽  
Turbulent Diffusion ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Laboratory Data ◽  
Integral Length Scale ◽  
Length Scale ◽  
Linear Superposition ◽  
Cylinder Diameter ◽  
The Mean

Laser-induced fluorescence was used to measure the lateral dispersion of passive solute in random arrays of rigid, emergent cylinders of solid volume fraction φ=0.010–0.35. Such densities correspond to those observed in aquatic plant canopies and complement those in packed beds of spheres, where φ≥0.5. This paper focuses on pore Reynolds numbers greater than Res=250, for which our laboratory experiments demonstrate that the spatially averaged turbulence intensity and Kyy/(Upd), the lateral dispersion coefficient normalized by the mean velocity in the fluid volume, Up, and the cylinder diameter, d, are independent of Res. First, Kyy/(Upd) increases rapidly with φ from φ =0 to φ=0.031. Then, Kyy/(Upd) decreases from φ=0.031 to φ=0.20. Finally, Kyy/(Upd) increases again, more gradually, from φ=0.20 to φ=0.35. These observations are accurately described by the linear superposition of the proposed model of turbulent diffusion and existing models of dispersion due to the spatially heterogeneous velocity field that arises from the presence of the cylinders. The contribution from turbulent diffusion scales with the mean turbulence intensity, the characteristic length scale of turbulent mixing and the effective porosity. From a balance between the production of turbulent kinetic energy by the cylinder wakes and its viscous dissipation, the mean turbulence intensity for a given cylinder diameter and cylinder density is predicted to be a function of the form drag coefficient and the integral length scale lt. We propose and experimentally verify that lt=min{d, 〈sn〉A}, where 〈sn〉A is the average surface-to-surface distance between a cylinder in the array and its nearest neighbour. We farther propose that only turbulent eddies with mixing length scale greater than d contribute significantly to net lateral dispersion, and that neighbouring cylinder centres must be farther than r* from each other for the pore space between them to contain such eddies. If the integral length scale and the length scale for mixing are equal, then r*=2d. Our laboratory data agree well with predictions based on this definition of r*.

A Sharp Interface Direct Forcing Immersed Boundary Approach for Fully Resolved Simulations of Particulate Flows

2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Jianming Yang ◽  
Frederick Stern
Keyword(s):  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Density Ratio ◽  
Sharp Interface ◽  
Immersed Boundary ◽  
Particulate Flows ◽  
Order Accuracy ◽  
Wide Range ◽  
Direct Forcing ◽  
Coarse Grids

In recent years, the immersed boundary method has been well received as an effective approach for the fully resolved simulations of particulate flows. Most immersed boundary approaches for numerical studies of particulate flows in the literature were based on various discrete delta functions for information transfer between the Lagrangian elements of an immersed object and the underlying Eulerian grid. These approaches have some inherent limitations that restrict their wider applications. In this paper, a sharp interface direct forcing immersed boundary approach based on the method proposed by Yang and Stern (Yang and Stern, 2012, “A Simple and Efficient Direct Forcing Immersed Boundary Framework for Fluid-Structure Interactions,” J. Comput. Phys., 231(15), pp. 5029–5061) is given for the fully resolved simulations of particulate flows. This method uses a discrete forcing approach and maintains a sharp profile of the fluid-solid interface. It is not limited to low Reynolds number flows and the immersed boundary discretization can be arbitrary or totally eliminated for particles with analytical shapes. In addition, it is not required to calculate the solid volume fraction in low density ratio problems. A strong coupling scheme is employed for the fluid-solid interaction without including the fluid solver in the predictor-corrector iterative loop. The overall algorithm is highly efficient and very attractive for simulating particulate flows with a wide range of density ratios on relatively coarse grids. Several cases are examined and the results are compared with reference data to demonstrate the simplicity and robustness of our method in particulate flow simulations. These cases include settling and buoyant particles and the interaction of two settling particles showing the kissing-drafting-tumbling phenomenon. Systematic verification studies show that our method is of second-order accuracy on very coarse grids and approaches fourth-order accuracy on finer grids.

scholarly journals Measurements of solid volume fraction and granular temperature at granular avalanches in a drum

2019 ◽  
Vol 49 (6) ◽  
pp. 067001
Author(s):  
Quan CHEN ◽  
Hui YANG ◽  
Ran LI ◽  
Ren HAN ◽  
QiCheng SUN
Keyword(s):  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Granular Temperature

scholarly journals Natural Convection over Two Superellipse Shapes with a Porous Cavity Populated by Nanofluid

Energies ◽  
2021 ◽  
Vol 14 (21) ◽  
pp. 6952
Author(s):  
Noura Alsedais
Keyword(s):  
Natural Convection ◽  
Rayleigh Number ◽  
Volume Fraction ◽  
Solid Volume Fraction ◽  
Darcy Number ◽  
Governing Equations ◽  
Porous Layers ◽  
Porous Cavity ◽  
The Mean ◽  
Volume Method

The influences of superellipse shapes on natural convection in a horizontally subdivided non-Darcy porous cavity populated by Cu-water nanofluid are inspected in this paper. The impacts of the inner geometries (n = 0.5,1,1.5,4) Rayleigh number (103 ≤ Ra ≤ 106), Darcy number (10−5 ≤ Da ≤ 10−2), porosity (0.2 ≤ ϵ ≤ 0.8), and solid volume fraction (0.01 ≤ ∅ ≤ 0.05) on nanofluid heat transport and streamlines were examined. The hot superellipse shapes were placed in the cavity’s bottom and top, while the adiabatic boundaries on the flat walls of the cavity were considered. The governing equations were numerically solved using the finite volume method (FVM). It was found that the movement of the nanofluid upsurged as Ra boosted. The temperature distributions in the cavity’s core had an inverse relationship with increasing Rayleigh number. An extra porous resistance at lower Darcy numbers limited the nanofluid’s movement within the porous layers. The mean Nusselt number decreased as the porous resistance increased (Da ≤ 10−4). The flow and temperature were strongly affected as the shape of the inner superellipse grew larger.

Sign in / Sign up

Export Citation Format

Share Document