I have noticed an issue for my runs using the resolved CFDEM solver (cfdemSolverIB) when attempting to use cyclic patches on the OpenFOAM side. For some reason, the velocity values at the cyclic patches do not match.
To check, I've run a toy case based on the "twoSpheresGlowinskiMPI" tutorial. The toy case is modified with gravity turned off and just including an initial uniform z-velocity in the cells specified in 0/U. Inlet and outlet (ie the upper and lower z faces of the domain) are changed to be paired cyclic patches.
When I compare the velocity values at the paired cyclic patches, I noticed that there is a fairly significant difference in some of the paired patch point values. In fact, none of the values match perfectly. By comparison, when running the case with only fluid phase (icoFoam), there is full agreement between the paired velocity values.
I've included below a link to the toy case with the runs and including the comparison of values.
Has anyone else had this issue before? If so any advice on how to remedy it would be greatly appreciated.
Thanks and all the best,