fix mesh/surface/stress not working in lastest CFDEM coupling?

Submitted by keepfit on Thu, 07/17/2014 - 19:17

Hi Chris,

Recently I upgraded the cfdem scheme from OpenFoam 2.1.1 to OpenFoam 2.2.2 + CFDEMcoupling 2.7.1+Liggghts 3.0.2, However,I am facing some issue regarding the mesh wall for CFD-DEM coupling:


fix mesh_wall all mesh/surface/stress file ../DEM/walls.stl type 1 wear finnie
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 1 meshes mesh_wall

In old version there is no problem importing STL mesh walls for coupling, but in the latest Version of cfdem scheme it seems not working.


...
Reading STL file 'walls.stl'
Resetting global state of Fix n_neighs_mesh_cad Style property/atom from restart file info
Resetting per-atom state of Fix n_neighs_mesh_cad Style property/atom from restart file info
nevery as specified in LIGGGHTS is overriden by calling external program (../cfd_datacoupling_mpi.cpp:57)
Resetting global state of Fix contacthistory Style contacthistory from restart file info
Resetting per-atom state of Fix contacthistory Style contacthistory from restart file info
Setting up run ...
[0] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigSegv::sigHandler(int) at ??:?
[0] #2 in "/lib64/libc.so.6"
[0] #3 LAMMPS_NS::FixContactHistoryMesh::sort_contacts() at ??:?
[0] #4 LAMMPS_NS::FixContactHistoryMesh::pre_exchange() at ??:?
[0] #5 LAMMPS_NS::Modify::setup_pre_exchange() at ??:?
[0] #6 LAMMPS_NS::Verlet::setup() at ??:?
...

However, if switching to primitive type wall it does work. In the restart file all of the "fix mesh/surface/stress" are "unfixed". If not unfixed, these meshes still exist in the restart file and "fix mesh/surface/stress" is not necessary any more?

Thanks for your help.

David