Combine SPH and DEM

Submitted by alejandro.gutierrez on Thu, 10/19/2017 - 08:15

Is possible resolve problems that have fluid and particles, for example mills, with a PUBLIC-LIGGGHTS

alejandro.gutierrez | Thu, 10/19/2017 - 08:27

Dear everyone
Is possible solve problems that have fluid and particles, for example mills, with a PUBLIC-LIGGGHTS.
I am running a SAG Mill and need incorporate water, is possible thats?. If anybody anybody have a script, I will be eternally grateful. I really want to know if a stl file can be used for SPH combined with DEM
Thanks

Alejandro Gutierrez
DIMEC-USACH

FelipeL | Fri, 10/20/2017 - 15:04

Hi

I´m Felipe Ligeti, how are you?. Oficcially, SPH + DEM is not implemented in PUBLIC-LIGGGHTS. The example of SPH that liggghts provide only works for a special atom style which is atom style sph, so you can´t add granular particles in that example because its a differente atom style (atom style granular). If you want to see a video of LIGGGHTS SPH example, i have done one:

https://www.youtube.com/watch?v=RVogLbsD4Ps

Also, is not possible to implemen stl files in that example (stl files cant act like a wall), however exist a user package called USER-SMD package for LAMMPS. It provides many material models (Johnson-Cook, plasticity with hardening, Mie-Grueneisen, Polynomial EOS) and allows new material models to be added. It implements rigid boundary conditions (walls) which can be specified as surface geometries from *.STL files. Maybe that could be a good starting point for your purpose. You can have a look of that package here:

http://lammps.sandia.gov/doc/Section_packages.html#user-smd-package

With this package, you can simulate something like this: https://www.youtube.com/watch?v=CwSqIuYRzMI

If you want more information about it, you can see this paper http://lammps.sandia.gov/doc/PDF/SMD_LAMMPS_userguide.pdf

Hope you are doing well, regards!

Nathan | Wed, 10/25/2017 - 04:44

Hi FelipeL,

This is interesting, thank for sharing. I will place this into my to-do-list.

Regards,

alejandro.gutierrez | Fri, 10/27/2017 - 07:46

I am running an example in Ubuntu-window 10, and have the following error:
ERROR on proc 1: Too many neighbor bins {../neighbor.cpp:1776)
ERROR on proc 2: Too many neighbor bins {../neighbor.cpp:1776)
ERROR on proc 3: Too many neighbor bins {../neighbor.cpp:1776)
ERROR on proc 4: Too many neighbor bins {../neighbor.cpp:1776)
MPI_ABORT was invoked on ranj 0 in communicator MPI_COMM_WORLD
with errorcode 1.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
Really I don't understand.
Thanks for you help

Alejandro Gutierrez
DIMEC-USACH