maarten-vd-sande / smoothed-particle-hydrodynamics Goto Github PK

Hi @Maarten-vd-Sande,

I'm very new to fluid simulation and I found your code very helpful. I'm wondering is it possible to make a viscous-elastic object with your code just by changing viscosity or other parameters. I'm going to implement two viscous-elastic spheres and merge them together like two bubbles (in 2D for the first step) something like the following video that I've done in python-blender:

I would appreciate if you share me your ideas on this issue.
Thanks

Hi, firstly I would like to thank you for the incredible work you've done.

I am working on implementing an sph algorithm, and your work has greatly helped.
The only issue, I will be displaying the results of this algorithm on an LED matrix, hence I am bound to a 64x64 window. I first tried reducing the simulation window's width, height and limits, however I run into an error whereby dens[i] remains zero, which causes a divide by zero in the integration script.
I traced dens[i] back to the if condition r_2 <= h_2 in compute density pressure, and have determined that the if condition was never satisfied, and so dens[i] remains zero, causing problems later on.
My next instinct was to increase the radius, h, in order to 'force' the if condition. If i set it larger than 63, I get list index out of range (obviously), setting it less I remain with the error above.

Any further suggestions would be greatly appreciated, thanks.