CFDEM®coupling - User Forum

I'm using OF-4.x to couple with CFDEM. While I compile CFDEM there's an error showing that the cfdemPostproc.C contains an error:

cfdemPostproc.C: In function ‘int main(int, char**)’:
cfdemPostproc.C:118:89: error: passing ‘const Internal {aka const Foam::DimensionedField}’ as ‘this’ argument of ‘void Foam::DimensionedField::operator=(const Foam::DimensionedField&) [with Type = double, GeoMesh = Foam::volMesh]’ discards qualifiers [-fpermissive]
make: *** [Make/linux64GccDPInt32Opt/cfdemPostproc.o] Error 1
waiting...

Dear All
i try to compile openfoam 4 with CFDEMcoupling-4.0. but i faced with below error:

using CFDEM_ADD_LIBS_NAME=additionalLibs_4.0 .
!!! ERROR !!!: CFDEM_ADD_LIBS_DIR/CFDEM_ADD_LIBS_NAME=/home/cfdemlab/CFDEM/CFDEMcoupling-PUBLIC-4.0/src/lagrangian/cfdemParticle/etc/additionalLibs_4.0 does not exist.
************************************
WARNING from your CFDEM code!
WARNING: Coupling with your OpenFOAM(R) version is currently under development!
In doubt switch to another version.
************************************

Hello! here is the installation tutorial which works fine for me

1)lets install OF231 following the procedure

http://openfoam.org/download/ubuntu.php

(instead of sudo apt-get install openfoam240 you should use sudo apt-get install openfoam231 or try this "cfdemSolverPiso is based on pisoFoam
and if you compare the options file of pisoFoam (version 2.4.x) with the same file from 2.3.x version you will notice that they differ by:

-I$(LIB_SRC)/meshTools/lnInclude

I have successfully installed Liggghts, cfdem and openfoam. Liggghts is perfectly running, went on to test cfdem coupling and I can execute the 'Allrun.sh' files in tutorial folders but no matter which I pick, CFDEM seems to be deleting the results files before I can have a look at them!

I check the case folder and I can see the cfdem dumps being created '0, 0.1, 0.2 etc' but when the run is complete, these are deleted....

Here is snippet of the terminal window at the end of the run:

Hi there,

i try to run a coupled simulation with servo wall pushing on some kind of particle column. An overview of the geometry is give in the attached screen shots. Actually, the particles are within a pipe, but for the sake of visibility this is not shown in the overview. The red circular face is the servo wall. The fluid region is meshed with snappy.
Running the simulation in DEM only is done without problems, however the coupled run dies at the first DEM run. From the log-file it seems that the mesh/surface/stress/servo makes some problems:

Please help me. I'm trying to compile CFDEM, I compiled with sucess openFoam 3.0.1 and LIGGGHTS 3.5.1, and the CFDEM compilation give-me this error:

I am trying to build CFDEM. I have successfully run
cfedmCompLIG
cfdemCompCFDEMsrc

cfdemCompCFDEMsol failed with error:

ld: cannot find -llagrangianCFDEM--v3.0+

Where is lagrangianCFDEM and how to build it?

Thanks,

Ping

Hi all,

I have probably got a bug when running LIGGGHTS coupled with CFD. I found that the write_start command in LIGGGHTS does not work properly when they are working in parallel with CFD. Write_start works fine in single running of DEM, but it does not save any updated state when running through CFD. I suppose specialists of CFDEM project can shed a light on this issue.

Thank you,

Regards,

Nathan,

Hallo,

I am trying to understand the implementation of the different force models in CFDEMcoupling and to extract the applied definition of the "drag coefficient" Cd0.

Looking at SchillerNaumannDrag, for me everything is clear. The drag force is given by the following expression (compare SchillerNaumannDrag.C):

drag=0.125*Cd*rho*M_PI*ds*ds*magUr*Ur

and I get the "drag coefficient" Cd0=Cd

Looking at DiFeliceDrag, I can also follow the implementation (compare DiFeliceDrag.C)

Hi there,

I have a simple test case where cfdemSolverIB gives very different results for the steady state compared to simpleFoam. The velocity field is in pretty good agreement, but the pressure field is completely different. Am I missing something fundamental about pressure in the IB solver or did I make a mistake somewhere in the case setup? Any help would be greatly appreciated.