CFDEMsolverIB FOAM fatal error

Submitted by bnaar on Mon, 02/01/2016 - 09:26

Dear all,

I am running the CFDEMsolverIB tutorial and I want to visualize my CFD and DEM results in paraview. To visualize this I obtained VTK files via adding reconstructCase="true" to the parCFDEMrun script & command and as I am running on multiple processors I had to add the "reconstructParmesh" command to the parCFDEMrun function.

Now when running the simulation I fall into the following error:

--> FOAM FATAL ERROR:
cell, tetFace and tetPt search failure at position (-0.00079004318 -0.001127873 0.0044064325)
for requested cell 0

From function void Foam::particle::initCellFacePt()
in file particle/particleI.H at line 758.

Have you experienced this error before and would you know how to solve this?

Much appreciated!

hunger's picture

hunger | Mon, 02/01/2016 - 12:27

One I had a similar error. In my case it was connected with lagrangian data dumped by the IOModel in CFDEM. If I remember right, it was a precision problem.
Since I am using the DEM data for postprocessing, I switched off the IOModel (IOModel_noIO). Then I could reconstruct my case.
You can give it a try with IOModel_noIO and rerun your simulation or simply remove the lagrangian data in your processor* directories by hand and give it another try with reconstructPar.

Best regards
Harald

Nucleophobe's picture

Nucleophobe | Fri, 02/26/2016 - 05:12

I've had this problem too. You can just run
'reconstructPar -noLagrangian'

as a quick work around (instead of manually removing the lagrangian data).

bnaar | Fri, 02/26/2016 - 14:54

Dear Nucleophobe,

Thank you for your response.

Could you maybe take a look at your previous post "cfdemSolverIB & rotating particles" as well. I would be very interested in applying your findings in my problem.

Best

Nucleophobe's picture

Nucleophobe | Mon, 03/07/2016 - 18:09

bnaar,

I shared my code in the other thread. I wanted to wait until I had some validation data to upload, but since it's taking me a while to get back to this I went ahead and shared the code with a simple example case.

-Nuc