Special Section: Computational Fluid Mechanics and Fluid–Structure Interaction

A DEM-FEM Coupling Approach for the Direct Numerical Simulation of 3D Particulate Flows

[+] Author and Article Information
B. Avci1

Institute of Continuum Mechanics,  Leibniz University of Hannover, Appelstrasse 11, 30167 Hannover, Germanyavci@ikm.uni-hannover.de

P. Wriggers

Institute of Continuum Mechanics,  Leibniz University of Hannover, Appelstrasse 11, 30167 Hannover, Germanywriggers@ikm.uni-hannover.de


Address all correspondence related to this paper to this author.

J. Appl. Mech 79(1), 010901 (Dec 13, 2011) (7 pages) doi:10.1115/1.4005093 History: Received March 06, 2011; Revised August 02, 2011; Published December 13, 2011; Online December 13, 2011

A computational approach is presented in this paper for the direct numerical simulation of 3D particulate flows. The given approach is based on the fictitious domain method, whereby the Discrete Element Method (DEM) and the Finite Element Method (FEM) are explicitly coupled for the numerical treatment of particle-fluid interactions. The particle properties are constitutively described by an adhesive viscoelastic model. To compute the hydrodynamic forces, a direct integration method is employed, where the fluid stresses are integrated over the particles’ surfaces. For the purpose of verifying the presented approach, computational results are shown and compared with those of the literature. Finally, the method is applied for the simulation of an agglomeration example.

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



Grahic Jump Location
Figure 3

Computed time history of the velocity in the direction of gravity for a sedimenting particle and the comparison of its terminal velocity with experimental correlations

Grahic Jump Location
Figure 2

Comparison of predicted values at different Reynolds-numbers for a flow past a fixed sphere: (left) Drag coefficient Cd versus Re, (middle) Separation length (Lw /d) versus Re and (right) Separation angle Θ versus Re

Grahic Jump Location
Figure 1

(Left) Distribution of the Lebedev quadrature points for NL  = 302. (Middle) Classification of the elements: 0 =: fluid element, i =: element of Pi, −i =: boundary element of Pi. (Right) Nonlinear weighted velocity update (○:= element center point E, ×:= fluid velocity node V).

Grahic Jump Location
Figure 4

Drafting-kissing-tumbling of two particles under gravity: t = 0.0, t = 0.4, t = 0.7, t = 0.9, t = 1.1, t = 1.4

Grahic Jump Location
Figure 5

The model of the scaffold consisting of 1544 adhesive spherical particles

Grahic Jump Location
Figure 6

Numerical results of the simulation showing the velocity field of the fluid and the particles at three points in time: (left, top) t=25, (left, bottom) t=34, (right, both) t=36. The inflow velocity at the inlet boundary is applied with uf  = 5. The flow direction is from left to right.



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.

Related Journal Articles
Related eBook Content
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