SPH numerical algorithm

Almin900

Hello everyone,
I am interested for the numerical algorithm of SPH. After the identifying particle neighbors, what is the next thing that should be calculated?
Is that density of the particle using initial density and after denisty pressure?
Regards

Alex

In a basic simulation (without motion, floatings or any special treatment):

1) You load initial state of the particles (initial position, velocity, density...)

2) You compute the neighbour list where we only store a list of particles ordering according to cells. (the domain is initially divided into cells of 2h)

3) You compute particle interactions, which means that you solve momentum equation (dv/dt) and continuity equation (drho/dt) for each particle.

4) You compute value of the time step (starting from forces, CFL conditions, viscosity...)

5) You update system: i) you compute new positions using new velocity which is computed using dv/dt, ii) new densities using drho/dt and iii) new pressure using the equation of state.

6) repeat 2-3-4-5

Following the structure of the code you can see this.

Regards

Almin900

Thank you!