How to add a Saffman force for particle in CFDEM@coupling

Submitted by yeting on Tue, 08/28/2018 - 05:00

hi! Everyone!
I am a beginner I want to add a Saffman force for particle in CFDEM@coupling, Who can help me ?Thanks !

yeting | Thu, 09/06/2018 - 13:26

Thanks! I noticed it ! i add this lift in /constant/couplingProperties. Then, The error was reported after calculating several steps, (Error in 'cfdemSolverIB': free(): invalid pointer: 0x000000000000004c8a3b8*****...........), the cfdemIB.C of this solver i was not modified, I don't know if it needs to be changed it . i have no idea,. what should i do next ? thanks again.

paul | Thu, 09/06/2018 - 15:57

This force model is not for use with the IB solver. Some guy on here did consider rotation in his IB calculations and published it - keep looking!

yeting | Sat, 09/08/2018 - 07:41

Thank you very much for you reply ! Now I add gradForce、viscForce etc,the error was reported after computing several steps. all these error is from cfdemSolver. However ,only ShigaokarIB and ArchimedesIB they can add. What if I want to add a new force ? I am a novice,so I'm not very familiar with how IB solver is used . How exactly is the force model added in IB? can you give me some advice? thank you!