[tjt]

Just wondering if anyone is aware of a forum for SPH anywhere?

I am looking at the SPH formulation for a pressure gradient. I have defined a x-direction gradient on a unit cube defined by 10000 particles. I have no boundaries. The SPH formulation approximates the gradient correction in the centre of the domain but gives an incorrect value at the x-dir edges due to truncation of the kernel radius.

Is there a way to correct for this?

Thanks in advance :-)

[beer]

Hi

The correction part is nicely described in this paper:
http://www.ibwpan.gda.pl/docs/ahem/ahem57str061.pdf
It works perfectly for gradients. I found it to be kind of unstable when used for the pressure gradient though. I'll might give it another shot later.

[tjt]

Many Thanks for this :-)