Incompressible SPH problems

ikypye · October 10, 2012, 06:08

Hello all,

Now, I'm doing research about SPH method for incompressible flow for my master degree. I use the weakly incompressible and incompressible SPH method. In incompressible SPH, which use projection method, there are 3 main steps as below
1. calculate intermediate velocity (v*) without pressure gradient term

2. calculate pressure implicitly by solve Pressure Poisson Eq. (PPE)
D(1/rho G(p)) = D(v*)/dt
D : divergence operator
G : gradient operator
dt : time step

3. correct the intermediate velocity by pressure gradient to get velocity in the next step

I have the problem in 2nd step. I have build the subroutine for this step and then I tried to test the code by simulate Poiseuille Flow. The flow is driven by body force not pressure gradient. From analytical, we know that pressure gradient is zero in flow direction and orthogonal to flow direction. Homogenous Neumann BC for pressure is applied at wall and ghost particles is used to model the wall. The result system of linear equation is singular and only have unique solution up to constant plus there is compatibility eqn for consistent system of linear equation i.e. integral volume of D(v*) must be equal to zero.

The problem is, in my code, integral volume of D(v*) not equal to zero. I think integral volume of D(v*) = 0 have physical explanation that mass of the system must be conserved.

Please give me some help. Is there something wrong with my code?

Thanks for your attention.

dbarca · January 6, 2013, 06:40

Did you get to have an answer to your problem?

FMDenaro · January 6, 2013, 17:47

Quote:

Originally Posted by dbarca

Did you get to have an answer to your problem?

That post is two-month old, however I think the error was in the BC.s that do not fulfill the compatibility contraint.

Mani Ordoubadi · May 24, 2013, 17:19

I think use of dummy particles instead of ghost particles gives better results for ISPH, for details you can refer to Lee's paper, A comparison of weakly and incompressible SPH, (I can't remember the exact name of the paper).
I did not understand your problem about Div(u*), but if you meant that you expected to have Div(u*)=0, you are wrong. u* is not a solution to the continuity eqn. It is merely a prediction for velocity field for the next time step. When you project this velocity field on the Divergence free space(that is by solving pressure poisson eqn), you will have the correct velocity field. (correction step).

October 10, 2012, 06:08	Incompressible SPH problems	#1
ikypye New Member Join Date: Oct 2012 Posts: 2 Rep Power: 0	Hello all, Now, I'm doing research about SPH method for incompressible flow for my master degree. I use the weakly incompressible and incompressible SPH method. In incompressible SPH, which use projection method, there are 3 main steps as below 1. calculate intermediate velocity (v) without pressure gradient term 2. calculate pressure implicitly by solve Pressure Poisson Eq. (PPE) D(1/rho G(p)) = D(v)/dt D : divergence operator G : gradient operator dt : time step 3. correct the intermediate velocity by pressure gradient to get velocity in the next step I have the problem in 2nd step. I have build the subroutine for this step and then I tried to test the code by simulate Poiseuille Flow. The flow is driven by body force not pressure gradient. From analytical, we know that pressure gradient is zero in flow direction and orthogonal to flow direction. Homogenous Neumann BC for pressure is applied at wall and ghost particles is used to model the wall. The result system of linear equation is singular and only have unique solution up to constant plus there is compatibility eqn for consistent system of linear equation i.e. integral volume of D(v) must be equal to zero. The problem is, in my code, integral volume of D(v) not equal to zero. I think integral volume of D(v*) = 0 have physical explanation that mass of the system must be conserved. Please give me some help. Is there something wrong with my code? Thanks for your attention.

January 6, 2013, 06:40		#2
dbarca New Member Daniel Barcarolo Join Date: Jul 2012 Posts: 11 Rep Power: 14	Did you get to have an answer to your problem? __________________ https://www.youtube.com/user/SPHflow fr.linkedin.com/in/dbarcarolo

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Incompressible boussiesq flow pressure BC piso and continuity problems	tehache	OpenFOAM Running, Solving & CFD	4	November 6, 2008 04:08
Abt: the speed of sound in incompressible flows	jinwon park	Main CFD Forum	3	August 21, 2008 21:44
Can a compressible solver solve incompressible ?	jinwon park	Main CFD Forum	17	December 13, 2007 06:00
Compressible -> incompressible.	Jinwon	Main CFD Forum	6	November 23, 2007 22:07
Some problems with Star CD	Micha	Siemens	0	August 6, 2003 14:55

May 24, 2013, 17:19		#4
Mani Ordoubadi New Member Mani Ordoubadi Join Date: May 2013 Posts: 1 Rep Power: 0	I think use of dummy particles instead of ghost particles gives better results for ISPH, for details you can refer to Lee's paper, A comparison of weakly and incompressible SPH, (I can't remember the exact name of the paper). I did not understand your problem about Div(u), but if you meant that you expected to have Div(u)=0, you are wrong. u* is not a solution to the continuity eqn. It is merely a prediction for velocity field for the next time step. When you project this velocity field on the Divergence free space(that is by solving pressure poisson eqn), you will have the correct velocity field. (correction step).