clear penalization force and pressure gradient force each loop
two particles could overlap and contain the same grid center
calculate volume fraction of a grid
restrict pressure gradient grids
apply principal axis decomposition
output width 15
Modified for the Q8P8 Element
It is modifed for the Q8P8 Element
*** empty log message ***