0
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

1

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

## Abstract

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.

<>

## Figures

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).

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

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

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

Figure 5

The model of the scaffold consisting of 1544 adhesive spherical particles

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.

## Errata

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 Proceedings Articles
Related eBook Content
Topic Collections