I have just implemented a particle system in XNA based on this paper and I am overall impressed with the result. Using as little as 100 particles beautiful liquid formations can still emerge.
There is a bug somewhere in my code and I could need some input from those who have implemented the paper.
Referring to Algorithm 4 'Spring Adjustment' I do not see the difference between 'L' and 'L_ij'. They are both the same, the rest length of the spring between the two particles 'i' and 'j', right?
Line 11 in Algorithm 4 is causing instability in my implementation when drops of ~10+ particles are forming. This is no matter which parameters are used in the simulation. Maybe it has to do with 'L' not being equal to 'L_ij'. Its a pain to debug as the drops start to jitter the more particles they consist of. Has someone had stability issues with line 11? (plasticity spring compression)
Last question. Are people adding a drag force to the particles to increase stability?
Best regards and thanks for the link