Thanks cjl for the reply.
That's correct, I do neglect gravity (and I also assume the fluid in question has no viscosity). I definitely haven't gotten to the full Navier-Stokes equations yet. But you say that incompressible inviscid flow can ignore energy and be represented by the system I already have if only I added in one more equation, and I'd be happy with that for the moment, so that's encouraging. But what is that equation exactly? Would that need to be this "equation of state" that you spoke of? All of this stuff so far is coming from a brief aside in my vector calc book, and I don't quite feel like tearing into any fluid mechanics books yet (though I have some). Is there anywhere I might look for a quick explanation, or is the explanation simply not quick?
Also, I'm confused when you say that the momentum equation I gave is for an incompressible fluid. To my understanding an incompressible fluid is defined as one for which \nabla \cdot V = 0, and this is by way of some interesting logic using the transport theorem which I don't feel the need to repeat here.
Yet, the way I understand the derivation of the momentum equation doesn't require using \nabla \cdot V = 0, or if it does then I don't see how. Nor does the resulting momentum equation imply \nabla \cdot V = 0, or if it does then (to repeat) lack of experience prevents me from seeing how. I'll give the derivation for that momentum equation here, and maybe someone can help me see why and how the result implies or requires incompressibility.
Firstly, here is the transport theorem just for reference. V is the velocity field of a fluid, \phi is the flow line function over that velocity field, W is some region of space, d^3V is a differential volume, and f(x, y, z, t) is some scalar function. The transport theorem is then
<br />
\frac{d}{dt}\int_{\phi(W)} f \; d^3V =<br />
\int_{\phi(W)}\frac{\partial f}{\partial t} + \nabla \cdot (fV) \; d^3V<br />
If we assume conservation of mass and apply the transport theorem with f = \rho, where \rho is the mass density function, then the continuity equation pops out. Specifically that is
\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho V) = 0
If the fluid is incompressible then this simplifies to
\frac{\partial \rho}{\partial t} + \nabla\rho \cdot V = 0
but we need not necessarily make that assumption for the rest of the derivations here in this post.
Now, to derive the momentum equation we first consider a volume of space which we will call W, the boundary surface of which we will call \partial W. If we approximate W as a polyhedron then \partial W will simply be a collection of nice flat surfaces which are very easy to work with. The force on a given surface caused by the fluid pressure will be approximately
\vec{F} = - \, p \cdot \vec{n} \cdot A
Where \vec{n} is the normal vector to the surface, A is the surface's area, and p is the pressure sampled at some arbitrary point on the surface and then assumed to be constant across the entire surface.
The total force on the whole volume as the polyhedron's sides increase in number to form a better approximation of the volume's surface is then found by the vector-integral equation
\vec{F} = -\int_{\partial W} p \cdot \vec{n} \cdot d^2A
If we momentarily consider only this equation's x component then we can do some interesting manipulations. We will fiddle things into a form amenable to the application of the divergence theorem, and to do this we will use a vector \vec{e} defined as \vec{e} = (1, 0, 0)
<br />
\begin{align}<br />
F_x =<br />
&-\int_{\partial W} p \cdot n_x \cdot d^2A = <br />
-\int_{\partial W} p (\vec{e} \cdot \vec{n}) \cdot d^2A = <br />
-\int_{\partial W} p \vec{e} \cdot d^2\vec{A} = <br />
-\int_{W} \nabla \cdot (p \vec{e}) \; d^3V = \\<br />
\\<br />
\\<br />
\\<br />
\\<br />
&-\int_{W} \nabla p \cdot \vec{e} + p(\nabla \cdot \vec{e}) \; d^3V = <br />
-\int_{W} \nabla p \cdot \vec{e} \; d^3V = <br />
-\int_{W} \frac{\partial p}{\partial x} \; d^3V<br />
\end{align}<br />
Similar arguments for the y and z dimensions using \vec{e} = (0, 1, 0) and \vec{e} = (0, 0, 1) then show that overall we have
\vec{F} = -\int_W \nabla p \; d^3V
Yet, the force on a moving region is also defined as the rate of change of its momentum over time, which is yet another vector-integral, which in this case is one we can apply the transport theorem to. We will write out the case for the x dimension, and other dimensions are similar.
<br />
F_x =<br />
\frac{d}{dt}\int_{\phi(W)} \rho V_x \; d^3V =<br />
\int_{\phi(W)} \frac{\partial}{\partial t}(\rho V_x) +<br />
\nabla \cdot (\rho V_xV) \; d^3V<br />
And so we have
<br />
\frac{\partial}{\partial t}(\rho V_x) + \nabla \cdot (\rho V_xV) =<br />
-\frac{\partial p}{\partial x}<br />
If we expand out the \frac{\partial}{\partial t}(\rho V_x) then we can use the resulting \frac{\partial \rho}{\partial t} term to combine this result with the continuity equation, and after the combination we have
\rho\left(\frac{\partial V_x}{\partial t} + \nabla V_x \cdot<br />
V\right) = -\frac{\partial p}{\partial x}
The more general vector form then reads as
<br />
\rho\left(\frac{\partial V}{\partial t} + \nabla V \cdot V\right) = - \nabla p<br />
Now, to my (perhaps untrained) eyes it seems that \nabla \cdot V = 0 is not a required assumption nor an implication for this momentum equation.
Over and out, and still a bit confused.
But wow, this stuff is fun.
