A FEM-DEM model for immersed granular flows

From unresolved to semi-resolved scales

Journal Article (2025)
Author(s)

Michel Henry (Université Catholique de Louvain)

Laura Valentina Cote Martinez (TU Delft - Civil Engineering & Geosciences)

Cristina Jommi (TU Delft - Civil Engineering & Geosciences, Politecnico di Milano)

Michael Hicks (TU Delft - Civil Engineering & Geosciences)

Jonathan Lambrechts (Université Catholique de Louvain)

Vincent Legat (Université Catholique de Louvain)

Miguel Cabrera (TU Delft - Civil Engineering & Geosciences)

Research Group
Geo-engineering
DOI related publication
https://doi.org/10.1051/epjconf/202534009016 Final published version
More Info
expand_more
Publication Year
2025
Language
English
Research Group
Geo-engineering
Journal title
EPJ Web of Conferences
Volume number
340
Article number
09016
Event
10th International Conference on Micromechanics on Granular Media, Powders and Grains 2025 (2025-12-08 - 2025-12-12), Candolim, Goa, India
Page Views
59
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

This work presents a coupling between the Finite Element Method (FEM) and the Discrete Element Method (DEM) to simulate immersed granular flows, transitioning from an unresolved to a semi-resolved representation. This refers to fluid discretisations ranging from a scale much larger than the grain size to one comparable to the grain size. The fluid phase is modelled using the Volume-Averaged Navier-Stokes (VANS) equations, which are solved with the FEM, while the granular phase is represented by the DEM with Non-Smooth Contact Dynamics (NSCD). An overlap-wise spatial coupling is proposed to bridge the gap between unresolved and semi-resolved scales. The method is validated by numerically reproducing experimental work on the fluidisation of a sand layer. Accuracy and stability are assessed by varying the mesh size down to half the grain diameter.