Register Blogs Members List Search Today's Posts Mark Forums Read

 August 30, 2012, 18:49 Derivation of Gradient in SPH #1 New Member   Scott Kulp Join Date: Feb 2012 Posts: 3 Rep Power: 6 I'm having a tricky time with the gradient using the SPH method. The standard SPH formulation makes sense: To find the gradient, integrate by parts: The first term is zero, since W(r,h)=0 for r>=h. So, we are left with However, the papers I am looking at [1] state That is, there is no minus sign. I feel that I'm missing some stupid mistake, but I've been looking at it for a while, and I don't see it. Thanks in advance! [1] J.J. Monaghan. "Smoothed Particle Hydrodynamics and Its Diverse Applications." Annual Review of Fluid Mechanics, Vol. 44: 323 -346. January 2012

 August 31, 2012, 07:51 #2 New Member   mehdi Join Date: Nov 2010 Location: Tehran Posts: 16 Rep Power: 7 Dear Scott Kulp It came from a little astute trick, I advise you to read page 43 and 44 of the book "Smoothed particle hydrodynamic particle a meshfree particle method" by G.R.Liu and M.B Liu in order to understand what happened.

 August 31, 2012, 09:46 #3 New Member   Scott Kulp Join Date: Feb 2012 Posts: 3 Rep Power: 6 This book is really nice and very informative. I see that equation 2.22 is basically the same as I wrote above (though 2.22 is the divergence), with the negative, and I do understand that. However, in the summary on page 44, equation 2.25 is not negative. I don't see any explanation on why the negative was removed... Any thoughts? Thanks!

 August 31, 2012, 11:20 #4 New Member   mehdi Join Date: Nov 2010 Location: Tehran Posts: 16 Rep Power: 7 I understood your question first time,however I think the answer is in page 44 where the authors wrote : "it should be noted that grad(i)Wij is taken with respect to particle i so the negative sign in equation (2.21) is removed in equation (2.25)." I agree with you there is an vagueness in the notation formulas,maybe we can say the mistake, but the concept is true, W is an even function so grad(W) is an odd one it means W in the equation 2.25 is the symmetry as W in equation 2.22 because the first is about all neighbors(which detected by j counters) and second is about the determined particle itself (i), maybe if you put Wji instead of Wij in eq 2.22 everything become correct.

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post fruitkiwi Main CFD Forum 4 June 12, 2012 01:12 diedro Main CFD Forum 0 September 7, 2011 07:45 sega OpenFOAM Running, Solving & CFD 14 August 12, 2009 06:21 Jay CFX 0 August 6, 2008 23:01 Paul Hsieh Main CFD Forum 3 November 11, 2003 05:52

All times are GMT -4. The time now is 21:39.