CFDEM-simulation with snappyHexMesh-Mesh

Submitted by Schrader on Mon, 03/06/2017 - 10:31

Hello,
I'm trying to run a simulation with OF 3.0.0 and LIGGGHTS 3.6.0 and CFDEM 3.6.0. The fluid mesh was created by the snappyHexMesh-tool. After a few timesteps I end up with an error, which is shown below. I'm using the snappyHex-Mesh-Tool for the first time. Is it possible that my mesh has a bad quality or is the cfdemSolverPiso not compatible with a mesh created by snappyHexMesh?


timeStepFraction() = 1
[3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ? in "/lib/x86_64-linux-gnu/libpthread.so.0"
[3] #3 Foam::interpolationCellPointFace >::findTet(Foam::Vector const&, int, Foam::Vector*, int*, int*, double*, double*, int&, double&) const at ??:?
[3] #4 Foam::interpolationCellPointFace >::interpolate(Foam::Vector const&, int, int) const at ??:?
[3] #5 Foam::viscForce::setForce() const at ??:?
[3] #6 Foam::cfdemCloud::setForces() at ??:?
[3] #7 Foam::cfdemCloud::evolve(Foam::GeometricField&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) at ??:?
[3] #8 ? at ??:?
[3] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #10 ? at ??:?
[user:16225] *** Process received signal ***
[user:16225] Signal: Floating point exception (8)
[user:16225] Signal code: (-6)
[user:16225] Failing at address: 0x3e800003f61
[user:16225] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x10330) [0x7f4c161bf330]
[user:16225] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0x2b) [0x7f4c161bf1fb]
[user:16225] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x10330) [0x7f4c161bf330]
[user:16225] [ 3] /home/user/OpenFOAM/OpenFOAM-3.0.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam26interpolationCellPointFaceINS_6VectorIdEEE7findTetERKS2_iPS2_PiS7_PdS8_RiRd+0x200) [0x7f4c1b76acc0]
[user:16225] [ 4] /home/user/OpenFOAM/OpenFOAM-3.0.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam26interpolationCellPointFaceINS_6VectorIdEEE11interpolateERKS2_ii+0x4af) [0x7f4c1b76bbcf]
[user:16225] [ 5] /home/user/OpenFOAM/user-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZNK4Foam9viscForce8setForceEv+0x42f) [0x7f4c19bc652f]
[user:16225] [ 6] /home/user/OpenFOAM/user-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam10cfdemCloud9setForcesEv+0xc6) [0x7f4c19b4b566]
[user:16225] [ 7] /home/user/OpenFOAM/user-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam10cfdemCloud6evolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERNS1_INS_6VectorIdEES2_S3_EES9_+0x15c0) [0x7f4c19b529d0]
[user:16225] [ 8] cfdemSolverPiso() [0x41fc11]
[user:16225] [ 9] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f4c15e0bf45]
[user:16225] [10] cfdemSolverPiso() [0x423b15]
[user:16225] *** End of error message ***

Best regards,

Marcel

Schrader | Mon, 03/06/2017 - 10:41

Hello,
I'm trying to run a CFDEM simulation with CFDEM 3.6.0. The fluid mesh was created by the snappyHexMesh-tool. After a few timesteps I end up with an error, which is shown below. I'm trying to use the snappyHexMesh tool for the first time. Is there any general problem to use this tool in a CFDEM simulation?


timeStepFraction() = 1
[3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ? in "/lib/x86_64-linux-gnu/libpthread.so.0"
[3] #3 Foam::interpolationCellPointFace >::findTet(Foam::Vector const&, int, Foam::Vector*, int*, int*, double*, double*, int&, double&) const at ??:?
[3] #4 Foam::interpolationCellPointFace >::interpolate(Foam::Vector const&, int, int) const at ??:?
[3] #5 Foam::viscForce::setForce() const at ??:?
[3] #6 Foam::cfdemCloud::setForces() at ??:?
[3] #7 Foam::cfdemCloud::evolve(Foam::GeometricField&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField, Foam::fvPatchField, Foam::volMesh>&) at ??:?
[3] #8 ? at ??:?
[3] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #10 ? at ??:?
[user:16225] *** Process received signal ***
[user:16225] Signal: Floating point exception (8)
[user:16225] Signal code: (-6)
[user:16225] Failing at address: 0x3e800003f61
[user:16225] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x10330) [0x7f4c161bf330]
[user:16225] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0x2b) [0x7f4c161bf1fb]
[user:16225] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x10330) [0x7f4c161bf330]
[user:16225] [ 3] /home/user/OpenFOAM/OpenFOAM-3.0.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam26interpolationCellPointFaceINS_6VectorIdEEE7findTetERKS2_iPS2_PiS7_PdS8_RiRd+0x200) [0x7f4c1b76acc0]
[user:16225] [ 4] /home/user/OpenFOAM/OpenFOAM-3.0.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam26interpolationCellPointFaceINS_6VectorIdEEE11interpolateERKS2_ii+0x4af) [0x7f4c1b76bbcf]
[user:16225] [ 5] /home/user/OpenFOAM/user-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZNK4Foam9viscForce8setForceEv+0x42f) [0x7f4c19bc652f]
[user:16225] [ 6] /home/user/OpenFOAM/user-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam10cfdemCloud9setForcesEv+0xc6) [0x7f4c19b4b566]
[user:16225] [ 7] /home/user/OpenFOAM/user-3.0.0/platforms/linux64GccDPInt32Opt/lib/liblagrangianCFDEM-PUBLIC-3.0.0.so(_ZN4Foam10cfdemCloud6evolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERNS1_INS_6VectorIdEES2_S3_EES9_+0x15c0) [0x7f4c19b529d0]
[user:16225] [ 8] cfdemSolverPiso() [0x41fc11]
[user:16225] [ 9] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f4c15e0bf45]
[user:16225] [10] cfdemSolverPiso() [0x423b15]
[user:16225] *** End of error message ***

Best regards,
Marcel

j-kerbl's picture

j-kerbl | Wed, 03/08/2017 - 16:51

Hi Marcel,

no I use snappy too and it is not an issue at all. It is tough to say what happened, the errors aren't really readable. Maybe something went wrong during the meshing, or maybe the simulation diverges? How high are the Courant Numbers? What does the output of checkMesh say?

Cheers,
Josef

Schrader | Mon, 03/20/2017 - 10:39

Hello Josef,
thank you for your response. I have changed some parameters for the meshing and now the the simulation is running.

Best regards,
Marcel