
Initialization. Set the initial condition of attributes of all SPH particles

Do the followings for all fluid particles until the run-time limit has been reached:

Calculate the current pressure of all particles by using state equation eq. (REF )

For each particle, insert it into the linked-list grid[1]} and register its nearest neighbour particles

For each particle, calculate its momentum change rate (eq. ()) by using the new calculated pressure, then calculate its new velocity and position

For each particle, calculate its density change rate (eq. (REF )) by using new calculated velocity, then calculate its new density


