lagrangian particle tracking fail in parallel
Hi guys,
I'm trying to implement a new hybrid lagrangian particle/FV solver for pdf method. But I'm experiencing a problem with parallization. It seems that the some particle can not get through the processor patch and the the solver traped in a dead loop. (the return value of trackToFace() is zero.) Did anyone have the same experience? Cheers, Edison |
Hi, guys,
I finally figured out what's wrong. When creating your new particle class be careful about how you write the particleIO.H. I find some weired particle private variable values after particle across the processor boundary. If you write the IO wrong then the code will read wrong memory location and the code will fail at parallel |
Hi Edison,
I'm working on the coupling of interFoam and solidParticle. The solver works well in serial. It runs in parallel as well, and the fluid phase behaves, however, the results of the solid phase are not reasonable, let alone resemble those produced in serial run. So there must be something wrong with the parallelisation of solidParticle, it runs without error message though. I checked the output, origId of the particles are all -1 at every time step. I suspect this might be the reason, however I didn't manage to get origId right, let alone to know whether this is the reason to the problem. Your posts are quite interesting to me. Could you shed some light on this please? Any thoughts would be much appreciated! Thanks in advance. Sophie |
All times are GMT -4. The time now is 04:33. |