SPH fluid simulation: 1st test review

Typical falling fluid simulation and review of the first test that I published on YouTube (by popular request, thanks at all for the comments and messages). I add the surface tension to eliminate the particle dispersion in the collision and the jet.

The color indicate the speed module (scaled to the max) using "jet color map" (extracted from matlab), the surface tension is a "spring-like" model using the surface particles only.

The code is a modification of SPH (named XSPH) using 2400 particles, and was programmed in Python-Pylab-Pygame.The animation was built in slow camera mode (1/15).

Category:

Science & Technology

Tags:

• SPH
• test
• simulation
• fluid
• particles
• dam
• breaking
• surface
• tension

[ricardomartins3]

Nice work :) I'd love o know how do you apply surface tension on the fluid i'd like to see if i coud aply them in the code i'm using

[jahkr]

Well, the model I used for surface tension is very simple (and probably wrong) and can be summarized in two steps:

1 - calculate the number of individuals that exist around each of the fluid particles in a radius equal to two or three times the length of smoothing

2 - If the number of particles is less than a certain amount (which occurs on the surface) add the calculated force (as if they were connected by springs) of the particle on its neighbors.

Thanks for your comment

[Cyber12321]

i wonder that how was your complexity?

[jahkr]

sorry but I don't understand the question.

[dibbuck]

Interesting and enjoyable. Why did you use python over matlab?

[jahkr]

I used matlab for many years and is faster than python, but I prefer to use free tools.

[ricardomartins3]

@jahkr Thank you for your reply. Apparently it works just fine. Does it affect much pressure results? Good luck for the future work

[ThaFacka]

Hi I'm into python and opengl, would it be possible to release (parts) of the source?

[Terabytekit]

I am starting the work on a SPH sim for blender based on kotsoft's code, but your code is already python so it could be a bit simpler to make it work in Blender.

Send me a message if you can

[Weatherlawyer]

Interesting stuff but some explanation to us outsiders might be considered.

I liked the view of the wave front contaminating the free flow out the beaker presumably a wave effect in the original resevoir caused the bifurcation that explains why there is always some residue?

Repeat the pour at enough height to eliminate backlash.