CFDEM®coupling - User Forum

Hi All,

I am running a coupled CFDEM case, with time steps of around 10e-10 , and coupling of every 10 steps. the simultion stops half way through with below error, I appreciate any help.

Thanks
ime = 6.554e-05

Courant Number mean: 2.82015e-05 max: 0.0036778
- evolve()
Starting up LIGGGHTS
Executing command: 'run 10 '
run 10 Setting up run ...
Memory usage per processor = 15.5843 Mbytes
Step Atoms KinEng 1 Volume
655391 1 0 1.3533016e-11 5.4e-09
655401 1 0 1.3533016e-11 5.4e-09

Thanks to all of the developers of CFDEM. I have been obtaining some encouraging results- http://youtu.be/mQN9WmHIW2c . One thing that was bothering me was dealing with pizza.py. It works well, but if it generates output at each coupling step, the dump file becomes enormous. I have attached two header files that output particle positions, radii, and velocities directly to Paraview readable file directories. Each folder has the same time stamp as the OpenFoam files. Thus, both pizza.py and foamToVTK become unnecessary.

Hi everyone,

I am a new user in CFDEM foram. I don't have any idea how to install it. Can anyone help me about the step by step procedure for its installation.

Regards,

Chandi

Hi Everyone,

I am running the settlingtestMPI example, where I get this error of time step continuity. I would appreciate any advise.

Time = 0.15

Courant Number mean: 5.34588e-05 max: 0.00517684

- evolve()

Starting up LIGGGHTS

Executing command: 'run 100 '

run 100 Setting up run ...

Memory usage per processor = 3.04352 Mbytes

Step Atoms KinEng 1 Volume

14901 1 0 5.1835199e-17 0.001

15000 1 0 5.1835199e-17 0.001

15001 1 0 5.1835199e-17 0.001

Hi

To get the particle sizes that I need, I am using LIGGGHTS in the "cgs" unit system, but OpenFOAM still works in the "si" unit system. I have converted all my LIGGGHTS values to work in the cgs system, however, when I combine the data in paraview, it is obvious that the two sets of data do not match up.

I presume that the coupling software uses the (dimensionless) numerical values as output by each programme, but does not check if they are of the same unit system? Therefore I want to create a conversion routine that converts the units during the coupling.

Hi,

I am trying to install cfdem coupling. I compiled OpenFOAM-1.7.1 from source. then I compiled 'liggghts-1p4p6' as library using 'make openmpi' and 'ln -s .....lmp_openmpi /usr/bin/liggghts' . The compilation of CFDEM src was successful. But when I try to compile the solver 'cfdemDolverPiso_shared', it gives the following error. Majorly "undefined reference to `Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&)'". Please see below. I was able to compile openFOAM solver 'icoFoam'. What may be the proble? any idea? or am i missing somethong obvious?

Hi everyone,

I came up to below error while compiling cfdem. I followed all instructions for coupled installations on ubantu 10.4 lts. I would appreciate any help.
thanks

Making dependency list for source file cfdemCloud/cfdemCloud.C
Making dependency list for source file derived/cfdemCloudIB/cfdemCloudIB.C
could not open file RASModel.H for source file derived/cfdemCloudIB/cfdemCloudIB.C
could not open file RASModel.H for source file cfdemCloud/cfdemCloud.C
Making dependency list for source file subModels/forceModel/forceModel/forceModel.C

Dear CFDEMers,

I would like to simulate the flow of a viscous (or even visco-elastic) fluid filled with glassfibres (about 1-4 mm in length, no spheres). Is that possible using OpenFOAM and LIGGGHTS? If not, are there any plans to include this feature?

Thank you in advance!
Thomas

Hi

What does DNUM 0 in solver output stand for? Is it okay to get this output?

Cheers
Claudio

Hi

I'm testing the adjustments in the couplingProperties-File. Some questions arose while working.
1) When I use the voidFractionModel centre the momentum coming from the particles in each cell act as one force on the fluid, right? And the momentum acting on the particles is the same for each particle?
2) What is the divided model for? When I try to use this (with dividedProps) I get an error message "type divided, constructor not in hash table". Do I have to compile something again?