Research Papers: Multiphase Flows

Modeling Blockage of Particles in Conduit Constrictions: Dense Granular-Suspension Flow

[+] Author and Article Information
A. J. Parry

Department of Modeling and Simulation, Schlumberger Riboud Product Center, 1 Rue Henri Becquerel, Clamart 92140, Franceaparry@clamart.oilfield.slb.com

O. Millet

Laboratoire d’Etude des Phénomènes de Transfert et d’Instantanéité: Agro-Ressources et Bâtiment, Université de la Rochelle, Avenue Michel Crepeau, Cedex 1, La Rochelle 17042, France

A refinement of this transfer uses a weighted projection (10).

Although not pursued in this work, methods are available to include the effects of instantaneous turbulent velocity fluctuations on the drag (12).

J. Fluids Eng 132(1), 011302 (Jan 07, 2010) (10 pages) doi:10.1115/1.4000691 History: Received December 10, 2008; Revised November 06, 2009; Published January 07, 2010

This paper presents a numerical simulation study of dense granular-suspension flow in a conduit with constriction. An empirical function of solid concentrations and Reynolds number prescribes the force between a particle and the fluid. This simplification reduces the computing load of the fine flow-field details around each particle. In the fluid-momentum equation, a source term distributes the force over the particle volume. The study addresses particle-laden flow at constant liquid-flow rate. Two different algorithms of the interparticle contact show that the bridging phenomenon causing the blockage of the particles persists in the presence of the fluid flow. While the particles are in movement, there are frequent interparticle and particle-wall impacts and vigorous fluctuations of the net reaction force on the wall from the particle phase. There is close correlation between the component of this reaction oriented in the flow direction and the rate of change in the pressure drop across the constricted conduit.

Copyright © 2010 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 10

Evolution of the pressure drop and total particle-wall reaction for the NSCD (top) and smooth DEM (bottom)

Grahic Jump Location
Figure 11

Time rate of change in the system pressure drop with total particle-wall vertical reaction for the NSCD (top) and smooth DEM (bottom)

Grahic Jump Location
Figure 12

Fluid speed (left: legend in m/s) and particle configuration in bridged state (right), showing force chains; smooth DEM method at 0.75 s

Grahic Jump Location
Figure 13

Smooth DEM in deep geometry: graph showing correlation of the rate of change in the system pressure drop with wall vertical reaction (left bottom), evolution of energy repartition (left top), and particle configuration at 0.215 s (right)

Grahic Jump Location
Figure 1

Dense granular-suspension flowing through a conduit constriction; combined solid-solid and solid-fluid interactions

Grahic Jump Location
Figure 2

Figure showing particle i in contact with particle j; local normal-tangential coordinates at contact point; normal unit vector collinear with particle centers

Grahic Jump Location
Figure 3

Distribution of the particle-fluid force onto multiple fluid cells: schematic of the particle on multiple fluid cells with force acting from particle to fluid (left), particle discretized into small cubes and force distributed uniformly onto the cubes (middle), and projection of the cube forces onto fluid cells (right)

Grahic Jump Location
Figure 4

Parameters for the wall proximity tests. Channel of length is equal to 12.5 sphere diameters, containing a single sphere at mid-length. The channel section parameters are L1, L2, and H (top). Seven sections are studied (middle and bottom). The solid volume fraction α was calculated according to Sec. 22 (bottom).

Grahic Jump Location
Figure 5

Wall proximity test, streamwise component of drag force; resolved and unresolved methods (Res=8000)

Grahic Jump Location
Figure 6

Drag forces on particles corresponding to the configuration studied in Sec. 3; test with particles fixed in the initial positions Resolved and unresolved methods (Res=8000)

Grahic Jump Location
Figure 7

Fluid and particle-phase coupling sequences. Exchange of data between fluid (top line) and particle (bottom line) solvers.

Grahic Jump Location
Figure 8

Geometry and initial particle configuration: 135 particles of 8 mm diameter

Grahic Jump Location
Figure 9

Particle-phase-energy balance for NSCD (top) and smooth DEM (bottom): ke=kinetic energy, pe=potential energy, wd_f=work done on particles by fluid, wd_c=energy dissipated by particle interaction, and sum=ke+pe−wd_f+wd_c




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In