- #1

- 82

- 0

As far as I can see, they're both separate ways of resolving the problem of keeping the flow divergence free. Am I wrong in thinking this?

You are using an out of date browser. It may not display this or other websites correctly.

You should upgrade or use an alternative browser.

You should upgrade or use an alternative browser.

- Thread starter Bucky
- Start date

- #1

- 82

- 0

As far as I can see, they're both separate ways of resolving the problem of keeping the flow divergence free. Am I wrong in thinking this?

- #2

arildno

Science Advisor

Homework Helper

Gold Member

Dearly Missed

- 9,970

- 134

What "pressure equation" are you talking about?

Euler's pressure equation? That only holds for certain inviscid flows, most notably for irrotational flow, but an analogue can be made when a Clebsch decomposition* can be made (i.e, when the helicity in closed vortex tubes is zero, rather than just constant (the latter being the general case for inviscid flows).

Anyhow, none of this is applicable for general NS problems; for those, the continuity equation (rather than the condition of a solenoidal field) is an indispensable fourth equation along with the three momentum equations.

Possibly, I've misunderstood what you are referring to; it's been awhile since I was slightly acquainted with all of this.

* The Clebsch decomposition is not a general decomposition like the Helmholtz decomposition, and looks like this:

[tex]\vec{v}=\nabla\phi+\alpha\nabla\beta[/tex]

where alpha and beta are concerned with, and in their form represents a limitation on, the vorticity in the fluid field.

Euler's pressure equation? That only holds for certain inviscid flows, most notably for irrotational flow, but an analogue can be made when a Clebsch decomposition* can be made (i.e, when the helicity in closed vortex tubes is zero, rather than just constant (the latter being the general case for inviscid flows).

Anyhow, none of this is applicable for general NS problems; for those, the continuity equation (rather than the condition of a solenoidal field) is an indispensable fourth equation along with the three momentum equations.

Possibly, I've misunderstood what you are referring to; it's been awhile since I was slightly acquainted with all of this.

* The Clebsch decomposition is not a general decomposition like the Helmholtz decomposition, and looks like this:

[tex]\vec{v}=\nabla\phi+\alpha\nabla\beta[/tex]

where alpha and beta are concerned with, and in their form represents a limitation on, the vorticity in the fluid field.

Last edited:

- #3

- 82

- 0

The Navier Stokes equations for incompressible fluid flow are the following:

[tex] \frac{\partial u} {\partial t} = -(u. \nabla )u - \frac {1} {\rho} \nabla p + v \nabla ^2 u + f [/tex]

[tex] \nabla .u = 0[/tex]

The second equation (the incompressibility equation) is the one I'm curious about.

I've read that this can be resolved through the pressure equation, or through the helmholtz-hodge decomposition.

One paper suggested this was done through substituting a pressure update formula into the divergence formula.

the pressure update formula is

[tex] u^{n+1} _{i+1/2,j} = u_{i+1/2,j} - \delta t \frac{1}{\rho} \frac{p_{i+1,j} - pi,j}{\delta x}[/tex]

[tex] u^{n+1} _{i,j+1/2} = u_{i,j+1/2} - \delta t \frac{1}{\rho} \frac{p_{i,j+1} - pi,j}{\delta x}[/tex]

substituting into the divergence formula gives

[tex] \frac{\delta t}{\rho} (\frac{4p_{i,j}-p_{i,j+1}-p_{i+1,j}-p_{i-1,j}-p_{i,j-1}}{\Delta x^2}) = -\frac{u_{i+1/2,j}-u_{i-1/2,j}}{\Delta x} + \frac{v_{i,j+1/2} - v_{i,j-1/2}}{\Delta x} [/tex]

The helmholtz-hodge decomposition is

[tex]

\xi = \nabla u + \nabla .v +h

[/tex]

where u is a scalar potential field (note that [tex]\nabla * (\nabla u) = 0[/tex]

where v is a vector potential field (note that [tex]\nabla .(\nabla * v) = 0[/tex]

where h is the harmonic vector field (note that [tex]\nabla .h = 0[/tex]

One of the papers I've read states that it uses a Poisson equation to derive a height field, which is subtracted from [tex]\xi[/tex] to yield a divergence free flow (it should be noted that this paper alltogether ignores the harmonic vector field). The paper can be found here:

http://www.dgp.toronto.edu/people/stam/reality/Research/pdf/GDC03.pdf

in a few papers I've read the authors have implied a relationship between the pressure and helmholtz-hodge solutions. However I don't see how they're related; in fact in that paper I just mentioned the author doesn't even acknowledge the pressure equation in his statement of the NS equations as pertaining to his solver, so I'm finding difficulty in tieing these things together.

I hope this has made my question more clear, and look forward to any answers.

- #4

- 1

- 0

P.S. Are you from CG?

- #5

- 14

- 0

Share:

- Replies
- 14

- Views
- 2K