Understanding the Navier-Stokes Equations for Smooth Particle Hydrodynamics

Bucky
Messages
79
Reaction score
0
Hi, hope this is the right area. Also please excuse me completley ignoring the template, I don't think it's applicable for the problem.

I'm an honours year student in a comp sci course and I've decided to do an implimentation of Smooth Particle Hydrodynamics in a 3d application as my topic. Currently I'm trying to understand the Navier-Stokes equations and figure out how to break them down into variables and functions.

My two comparisons of the equation at the moment are a paper on the subject:
http://graphics.ethz.ch/Downloads/Publications/Papers/2003/mue03b/p_Mue03b.pdf

and the wikipedia article on navier stokes:
http://en.wikipedia.org/wiki/Navier-Stokes_equations


that muller 03 paper lists the navier-stokes equation as:

ρ (∂v / ∂t + v·∇v) = −∇p+ρg+μ∇2v,

(am I right in saying that ∂v / ∂t is just acceleration?)

whereas wikipedia gives it as:

ρ (∂v / ∂t + v·∇v) = −∇p + f +μ∇2v,

wheeere...
ρ = density
v = velocity
t = time
p = pressure
μ = viscosity

the notable change is f (representing external forces) to ρg. where, rho represents density, and g represents a density field(?).

firstly, I don't get these 'field' things. are they just the weighted average of certain variables of particles surrounding the particle we're currently looking at? is this what smoothing kernels are used for?

if i "solve" this equation, what do i actually end up with?

also the upside down triangle (grad?), represents the gradient of the element? how do you even find that?
 
Last edited by a moderator:
Physics news on Phys.org
Bucky said:
also the upside down triangle (grad?), represents the gradient of the element? how do you even find that?

I believe that is the del operator...http://en.wikipedia.org/wiki/Del
 
ok i looked at that wiki article, and del is used to mean a lot of things. what's it supposed to be representing here? I noticed it has some relation to fields too :/


if i was planning to use this in an application, say to find the acceleration each frame. could i just rearrange the equation to get the acceleration term on it's own (or with velocity, since that'd be a known variable)?
 
if you didn't found anything
I recommend you the fallowing papers:
http://liu.diva-portal.org/smash/get/diva2:324983/FULLTEXT01
Here they are showing some kernel functions

also a paper with a different (hack) approach that is easier to implement:
http://www.iro.umontreal.ca/labs/infographie/papers/Clavet-2005-PVFS/pvfs.pdf
 
Last edited by a moderator:
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top