6

I have a hard time to wrap my head around pressure in the Navier-Stokes equation! It may sounds ridiculous but still I cannot understand the true meaning of pressure in the Navier-Stokes equation. Let's do some math to explain my purpose more accurately! Let's start from basics of the physics and in my opinion that would be the first equation in the classical thermodynamics as equation of state. We assume: there is a fluid, which have a equation of state as:

$$\rho = \rho(P,T)$$

Where $\rho$ is the density of the fluid, $P$ is the pressure, and $T$ is the temperature. Let's take a derivative from this equation to have:

$$d\rho = (\frac{\partial \rho}{\partial P})_{T} dP + (\frac{\partial \rho}{\partial T})_{P} dT$$

Let's assume that our fluid is in the thermal equilibrium and its temperature will not change, as a result: $d T = 0$

So, we have:

$$d \rho = (\frac{\partial \rho}{\partial P})_{T} dP$$

I know it's a lot of assumption but again let's assume that density change because of pressure change is not nonlinear and our fluid in fact behaves like a ideal gas. As a result, I call $(\frac{\partial \rho}{\partial P})_{T}$ the inverse square of the speed of sound, which is a constant number, as:

$$(\frac{\partial \rho}{\partial P})_{T} = c_{s}^{-2}$$

So, finally we have:

$$d \rho = c_{s}^{-2} d P$$

Or:

$$\Delta \rho = c_{s}^{-2} \Delta P$$

Or again:

$$(\rho - \rho_{f}) = c_{s}^{-2} (P - P_{0})$$

Where $\rho_{f}$ is the density of the fluid at the rest or reference, which is a tabulated value for each fluid, and $P_{0}$ is the reference pressure.

Now, I would assume my fluid is an incompressible fluid and it means (density is constant and it is really constant!):

$$\rho = \rho_{f}$$

As a result, because, every fluid regardless of its compressibility or incompressibility has a finite speed of sound, I would argue that:

$$P = P_{0}$$

Or in other word, strictly speaking pressure should be equal to the reference pressure.

Now, I proved that for an incompressible fluid as long as density is constant, pressure should also be a constant. So in incompressible Navier-Stokes equation we have:

$$\rho_{f} \frac{\partial \mathbf{u}}{\partial t} + \rho_{f} (\mathbf{u} \cdot \nabla)\mathbf{u} = -\nabla P + \nabla \cdot \tau$$

And I showed that for incompressible fluid, P is just constant, so: $\nabla P = 0$!

As a result, I could simplify the Navier-Stokes equation as:

$$\rho_{f} \frac{\partial \mathbf{u}}{\partial t} + \rho_{f} (\mathbf{u} \cdot \nabla)\mathbf{u} = \nabla \cdot \tau$$

Now let's back to my original question:

Based on these calculations I would say that pressure in the incompressible Navier-Stokes equation is just a dummy variable, which does not have any physical meaning! I appreciate if someone could explain this to me!

Qmechanic
  • 201,751
GGG
  • 331
  • 2
    there are relevant cases in which $dT$is not neglectible even for incompressible fluids. This is the case for e.g. describing air movement on large scales. So P as a function of T is relevant in this case and $\nabla P$ is not necessary 0 – Katermickie Sep 05 '18 at 21:40
  • @Katermickie I know what you mean but let's just limit our analysis to those cases where $d T = 0$, which are not rare in the literature! I didn't say that my analysis could cover all the cases but as long as my assumptions are hold true, I'm interested to study the consequences. I'm truely interested in these assumptions because the application that I use Navier-Stokes for, is completely match to these assumptions! – GGG Sep 05 '18 at 21:44
  • @Katermickie And also please read the question more carefully because your example when air moves on large scales probably should be studied by using compressible Navier-Stokes equation and you need to couple thermal energy balance equation (i.e. heat transfer equation) to the Navier-Stokes! As a result, your example is not relevant to my assumptions! – GGG Sep 05 '18 at 21:47
  • Related, if not dupe of, https://physics.stackexchange.com/q/319577/25301 – Kyle Kanos Sep 05 '18 at 21:49
  • @Katermickie In another word, I have no doubt about compressible Navier-Stokes equation but I believe the incompressible version is just a crap! I mean a lot of people including myself use it for research but is it possible someone explain the logic behind it to me?! – GGG Sep 05 '18 at 21:49
  • The only loop hole may exist in my derivation is that you can argue about my finite speed of sound! For a truely incompressible fluid, speed of sound should be infinitely large! (at least theoretically!!!!), but then you need to justify why that's the case for an incompressible fluid to have a infinite thermal velocity!? If someone could explain this also, I would appreciate his/her answer! – GGG Sep 05 '18 at 22:06
  • Actually on large scales for small velocities you can treat air as incompressible fluid but besides that I would assume that youre right if you leave out those cases you can just ignore the pressure term in most cases. – Katermickie Sep 05 '18 at 22:06
  • One case in which you cannot ignore it that I could Think of is a Tank filled with Water and a tilted Water Surface i.e. the height difference between the bottom of the Tank and the Surface of the Water becomes a function of Space. Since the pressure P on the bottom of the Tank depends on the height of the water above it also becomes nom constant ans $\nabla P$ is again nonzero. Regarding this I am Not Sure if your derivation using $\rho=\rho(P)$ is valid if you want to assume incompressible fluids where $d\rho/dP=0$... – Katermickie Sep 05 '18 at 22:06
  • If you are talking about hydrostatic the things even could get more complicated! Yes, let's say our fluid even does not move and as a result there should be pressure gradient because of the gravity force but how you can change my derivation to show that pressure should not be a constant! I mean otherwise you will end up arguing that because I have pressure gradient along a tank, I should have a density gradient! which is in my opinion is very unlikely! – GGG Sep 05 '18 at 22:10
  • Or I would say every conservative force field (i.e. $f = -\nabla \phi$) where $\phi$ is the potential of that force, which gravity is an example for that force field, its potential $\phi$ could be interpreted as a dummy pressure. I mean we could call it pressure but it's not really thermodynamics pressure! – GGG Sep 05 '18 at 22:16
  • Since for incompressible fluid $d\rho/dP=0$ you do not get a density gradient. So if you assume a tilted Water Surface the height changes along lets say x direction . The hidrostatic pressure at a point is $P=\rho gh$ where h is the height of the overlying Water. If $h=h(x)$ you get a pressure Gradient along x – Katermickie Sep 05 '18 at 22:19
  • An example of $f = -\nabla \phi$ but that's $\phi$ or potential of conservative force. So why I should call it pressure?! Thermodynamics pressure is something like a internal energy. You know it should depend on what you impose from the outside to the matter. The thermodynamics pressure is just a internal energy per volume and that's it. Regardless of presence of gravity or other forces, there is a pressure which does not depend on kinetics variables and only depends on equilibrium potentials of density and temperature. I admit external force will add something to thermodynamics pressure... – GGG Sep 05 '18 at 22:25
  • But I take account for all external forces in the Navier-Stokes equation separately: $\rho_{f} \frac{\partial \rho}{\partial t} + \rho_{f} (\mathbf{u} \cdot \nabla) \mathbf{u} = -\nabla P + \nabla \cdot \tau + \mathbf{f}$, where $\mathbf{f}$ is the all external force fields. So why I should include it in P itself? I could say ok at the hydrostatic case the external field is in the balance with the internal pressure. – GGG Sep 05 '18 at 22:27
  • 2
    Why are you saying that the pressure is constant for an incompressible fluid? – Chet Miller Sep 06 '18 at 00:26
  • Related (the role of pressure in incompressible hydrodynamics): https://physics.stackexchange.com/q/507526/226902 and https://physics.stackexchange.com/q/187889/226902. – Quillo Jan 16 '23 at 11:03

2 Answers2

8

There are two pressures: Thermodynamic pressure $p_\text{thermo}$, and Mechanical pressure $p_\text{mech}$. Thermodynamic pressure, a concept from equilibrium thermodynamics and therefore applicable only to a static fluid, is given by an equation of state: $p_\text{thermo}=f(\rho,T)$, where $\rho$ is fluid density and $T$ its temperature. A moving fluid is not in equilibrium and its $p_\text{thermo}$ is not defined. Mechanical pressure is the isotropic part of the stress tensor and is defined for a moving fluid too; $p_\text{mech}$ appears in the Navier-Stokes equation.

If a static fluid is isothermal and has constant density ($\rho,T$ fixed) then $p_\text{thermo}$ is also fixed. But mechanical pressure given by hydrostatic equation varies with depth in an isothermal constant-density fluid.

Your derivation mixes up the two pressures. The relations are: \begin{align} \mathrm{d}\rho&=c_s^{-2}\,\mathrm{d}p_\text{mech}\\ \mathrm{d}\rho&=\left(\frac{\partial\rho}{\partial T}\right)\,\mathrm{d}T+\left(\frac{\partial\rho}{\partial p_\text{thermo}}\right)\,\mathrm{d}p_\text{thermo}. \end{align} The latter equation from thermodynamics is applicable only to a static fluid. The former is not a thermodynamic equation. An incompressible fluid is defined to be one whose density doesn't depend on its mechanical pressure $p_\text{mech}$; it doesn't say that $p_\text{mech}$ can't vary. Therefore as you approach the limit of an incompressible fluid, $\mathrm{d}\rho\to0$, we must necessarily have $c_s\to\infty$. It is incorrect to say "...every fluid regardless of its compressibility or incompressibility has a finite speed of sound..."; incompressible fluids don't exist, so a priori you wouldn't know what sound speed should be assigned to a hypothetical fluid; to be consistent with the definition of incompressibility however a variation in $p_\text{mech}$ must be allowed which demands that sound speed in an hypothetical incompressible fluid be infinite.

P.S. Here's article1 and article2 that may interest you (NB: both are PDFs).

Kyle Kanos
  • 28,229
  • 41
  • 68
  • 131
Deep
  • 6,504
  • You said the equation $d \rho = c_{s}^{-2} d p_{mech}$ is not a thermodynamic equation. If it's not how you can derive it? I mean it's an ideal gas equation of state, so why there should be two definitions for pressure as thermodynamical and mechanical?! My approach to Navier-Stokes equation is a little bit different from conventional fluid dynamicists and as a result I see that equation ($d \rho = c_{s}^{-2} d p_{mech}$) everyday because I use lattice Boltzmann method for solving Navier-Stokes, and still I can't digest why ideal gas equation of state should give me pressure in Navier-Stokes! – GGG Sep 06 '18 at 13:25
  • @MehrdadYousefi Here's a derivation. What is needed for the derivation is the concept of bulk modulus. For an ideal gas, if you assumed that thermodynamic and mechanical pressures are the same (this assumption is often made in other contexts as well), then you can use its equation of state to explicitly compute the bulk modulus in terms of known properties of the gas. Also no matter how you derive the Navier-Stokes equation, a distinction between thermodynamic and mechanical pressure is necessary... – Deep Sep 07 '18 at 05:23
  • ... because a moving fluid is not in equilibrium and the concept of thermodynamic pressure is inapplicable to it, while its mechanical pressure is well defined. "My approach to Navier-Stokes equation is a little bit different from conventional fluid dynamicists ..." Could you post a relevant link? I would be interested in an alternate derivation. – Deep Sep 07 '18 at 05:27
  • It's a good starting point to be familiar with my approach "lattice Boltzmann method" to solve Navier-Stokes instead of using conventional FEM or FVM flow solvers: http://www.dartmouth.edu/~cushman/papers/2018-Boltzmann-to-NS.pdf – GGG Sep 07 '18 at 13:24
  • @MehrdadYousefi Thanks for the link. I ought to discuss with my colleagues who use LBM. In any case I think the pressure computed using LBM can't be thermodynamic pressure because a moving fluid is not in equilibrium. But applying equilibrium thermodynamic relations to a non-equilibrium situation (moving fluid) is an excellent approximation if the deviation from equilibrium is not too large. Usually this means that time scale of the flow is much larger than the molecular time scale (for eg. mean collision time). – Deep Sep 08 '18 at 04:27
  • Yes, in fact LBM leads to a slowly time varying solution from Chapman-Enskog expansion where diffusive and convective time-scales are decoupled. Ok, I can understand when non-equilibrium condition is not too far from equilibrium condition we could apply thermodynamics relation, but the problem is why this equilibrium should suddenly coincide with an ideal gas?! I mean it means any fluid near the equilibrium behaves like an ideal gas?! That's strange to me! I compared LBM to FEM and more strangely the pressure from LBM is equal to FEM which is derived from pressure Poisson equation!? – GGG Sep 08 '18 at 13:52
  • Also look at this discussion here which is not convincing in my opinion: http://www.palabos.org/forum/read.php?3,4809 – GGG Sep 08 '18 at 13:57
3

For an incompressible fluid, the pressure in the Navier Stokes equation represents the isotropic part of the stress tensor. It is determined up to an arbitrary constant value; that is, adding an arbitrary constant to the pressure at all location throughout the flow field still enables it to satisfy the NS equation. The arbitrariness in removed by specifying the pressure at any location on the boundary.

Chet Miller
  • 33,586
  • 3
  • 21
  • 48