[ayrton88]

Hi,

I'm using WCSPH to simulate Poisseuille Flow.
I'm using Morris' method for almost everything except I'm using dummy particles as the stationary plates. I have periodic boundaries and I've checked that the number of neighbours is correct throughout the simulation. The density remains constant as well.
The problem is that even though initially the solution seems to be correct and the velocity starts to get the parabolic profile, but after a point the acceleration starts to increase (instead of decreasing)
Does anyone have any idea of what may be wrong?
If it didnt work at all it would be easier to debug, now it works and then it doesnt...
I'm trying to debug it for a couple of weeks but cant find something wrong in the code and I have run out of ideas of what can be wrong.
Thanks!