
[Sponsors] 
March 13, 2015, 16:01 
Physics of an incompressible fluid

#1 
Senior Member
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 128
Rep Power: 18 
It is useful to consider the necessary properties of a hypothetical incompressible fluid in the same sense that one considers an inviscid fluid as the limiting case of vanishing viscosity. Consider a small volume of fluid. It is usually assumed that the effect of the surrounding fluid can be represented in terms of surface (contact) forces or stresses. For compressible fluids, the assumption of a normal surface force (due to pressure) correctly leads to a gradient term in the momentum balance equation, which is the basis of the NavierStokes equation. It is conventionally assumed that the same holds true in the incompressible case.
But by the HelmholtzHodge decomposition, the incompressible NavierStokes equation is composite, separable into a pressurefree equation for the divergencefree velocity and an equation for P as a function of U. So we see a contradiction. The assumption of a normal force leads to a grad P term, while the Helmholtz decomposition excludes such a term. We must conclude that an incompressible fluid does not support the assumption of a normal surface force. If not a surface or contact force, what force keeps the velocity divergencefree? We invoke the classical concept of actionatadistance, resulting in a body force. In general, the advective/convective term (U dot grad U) is not divergence free. When the divergence of (U dot grad U) is positive, this 'virtual' fluid source acts like a repulsive force, pushing adjacent fluid away. When negative, there is a 'virtual' sink, resulting in an attractive force as adjacent fluid is pulled in. We consider the form of the force on fluid particle at r due to a nonzero divergence at r’. Assume the strength of this force is proportional to the divergence of (U dot grad U), the force acting along the line joining r and r’, and decreasing with distance inversely as the distance squared in 3D and inversely as the distance in 2D. Consequently the force on the particle at r due to the divergence of the advective term at r’ would be f(r,r’)=rho Cd (div (U dot grad U)) (rr’)/rr’^d , where Cd is a constant to be determined. Summing over all the distant positions r’ gives the force on the fluid particle at r. With a(r)=f(r)/rho, a(r)=integral(Cd (div (U dot grad U)) (rr’)/rr’^d dr’). We require the divergence of ( U dot grad U+a(r) ) to be equal to zero. In free space, the integral can be evaluated to find the coefficient Cd. So the governing equation for incompressible velocity U is an integrodifferential equation. There is an elegant way to solve the variational/weak form of the pressureless incompressible equation using finite elements. Question: Can integrodifferential equations like this be solved using the methods of finite differences? 

March 13, 2015, 17:18 

#2 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73 
Hello Jonas,
this statement is not clear to me: The assumption of a normal force leads to a grad P term, while the Helmholtz decomposition excludes such a term. We must conclude that an incompressible fluid does not support the assumption of a normal surface force. However, FD can be used as way to discretize any differential operator (provided regularity of the function is ensured), of course when an integral needs to be discretized you have to complete the FD formulation with some integration rules such as trapezoidal, Simpson.. 

March 13, 2015, 18:27 

#3  
Senior Member
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 128
Rep Power: 18 
Quote:
(1) The assumption (existence) of a normal surface force (pressure) leads to a grad P term (force). True? (2) The Helmholtz decomposition excludes such a (grad P) term. True? These statements are incompatible: ((1) there is ... or (2) there isn't ... a grad p term. They can't both be true). True? How can we make them compatible except excluding some assumption about incompressible fluids? What other assumption is there than the existence of a normal surface force? It took me a long time to accept this conclusion and its consequences. I am still working on this. Quote:


March 13, 2015, 18:44 

#4 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73 
I try to clarify my doubt...Consider the momentum equation written in terms of accelerations:
a + Grad p = a* where a is divergencefree Grad p = Div (I p) is curlfree a* is convection + diffusion therefore the decomposition of a* is already present in the equation. You see also that that action of the pressure is by the isotropic part. I dont see contradiction... PS: yes, FEM, FV and FD could be seen in a unified framework of projective methods along suitable shape functions. 

March 14, 2015, 12:16 

#5 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
I'm not sure I understand the question/topic, but for steady incompressible flow, it is conservation of mass which keeps the flow divergence free.


March 14, 2015, 12:28 

#6 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
Oops, let me change my statement. Conservation of mass keeps incompressible flow divergence free regardless of whether the flow is steady or not since drho/dt is neglected.


March 14, 2015, 12:33 

#7 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
BTW, what happens in real life is that an acoustic wave goes out and sets up this divergence free flow. And for incompressible flow it moves infinitely fast. It is the acoustic wave which provides the "action at a distance." But the acoustic wave is a pressure wave.


March 14, 2015, 12:45 

#8 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73 
yes, no matter steady/unsteady, the mass equation is alway Dv v = 0, for steady flow the decomposition writes as
Grad p = a* the incompressible model tell us that a small variation of density (d rho) causes istantaneusly an infinite variation of pressure in isoentropic condition, such that dp/d rho > Infinity Of course, that is a model, not the real physics in which the compressibility is always present 

March 14, 2015, 12:55 

#9 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
And, I gather that the original posters "P as a function of U" is Bernoulli's equation.


March 14, 2015, 13:11 

#10 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
Oh, I guess the other assumption that the original poster is missing is that the pressure wave moves very fast relative to the motions of interest. Another way to view the incompressibility assumption is that the speed of sound is much much larger than the velocity of the object and/or flow.


March 14, 2015, 13:48 

#11 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73 

March 18, 2015, 14:03 
Physics of an incompressible fluid

#12 
Senior Member
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 128
Rep Power: 18 
I posted this thread hoping to prompt more discussion about the physical properties one might find for a hypothetical incompressible fluid (in contrast to incompressible flow), which might be used as a model for incompressible fluid flow. One thinks of a hypothetical inviscid fluid without viscosity. That concept is very useful, but has limitations near boundaries, where viscosity, however small, dominates. My tact was to model the force between fluid particles in terms which would lead to the incompressible NavierStokes equation.
In the classical mechanics of the 19th century, forces were viewed as due to contact, or by actionatadistance. Examples of the latter were gravitational (Newton), electric (Coulomb), and magnetic (BiotSavart). Maxwell united the latter two with the concept of an electromagnetic field as an intermediary between charged particles. The interactions of Coulomb and BiotSavart were instantaneous (effects propagated at an infinite velocity, if you wish). Einstein’s theory of relativity, with a limiting velocity of interactions, changed that. Although there were theories of direct interaction using ‘retarded’ potentials, field theories gained acceptance over actionatadistance. For fluids this means fluid particles interact via pressure fields, interactions limited by the speed of sound in the fluid. For incompressible fluids, the speed of sound would be infinite. It would seem that to accommodate this fact, we need to resurrect the classical concept of direct interactionatadistance. The form I described for the force in the post is just a ‘physical’ interpretation of the terms in a Green’s function for the projection. In response to Martin Hegedus: (1) A basic assumption is that the velocity of pressure waves (if they exist in this fluid) would be infinite. (2) Incompressibility (or div U = 0) is treated as a conservation law, not the agency of a thermodynamic pressure. (3) Bernoulli’s equation is just a statement that the incompressible part of ( U dot grad U) is invariant along the stream line. (4) The solenoidal projection of the NavierStokes equation does not contain a pressure, and divergence free flow can be computed without reference to a pressure, so what evidence exists that a pressure exists in an incompressible fluid? Isn’t it used as just another name for a subtractive term in the solenoidal projection? In response to Filippo Denaro: As you note, P^S (U dot grad U  grad P) = (U dot grad U  grad P’) . The form is conserved but grad P is not necessarily the same as grad P’. By P one might mean pressure, but P’ would be something else. This argument worries me too because it seems a little circular. The result comes from writing P^S=(1P^I), where P^S and P^I are the solenoidal and irrotational projection operators. But we could formally write P^S=(curl((curl curl)^(1))curl) and P^I=1P^S instead of P^I=grad(Delta^(1))div, in which case a grad P term would never appear. Of course we would have to find an expression for the Green’s function of the (curl curl) operator. This is much closer to the modified Hermite finite element method I use. Still unresolved in this discussion is how to handle the influence of boundaries, and the force an incompressible fluid might exert on a boundary. Last edited by Jonas Holdeman; March 18, 2015 at 14:06. Reason: easier to read 

March 18, 2015, 14:33 

#13 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
Sorry, a little lost here but this sounds like a discussion about potential flow (Panel methods, etc.)


March 18, 2015, 14:46 

#14 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
BTW, neither Bernoulli's equation or potential take into account viscosity, but, I assume, that if this discussion applies to the N.S. equations, it also applies to the Euler equations.


March 18, 2015, 14:50 

#15 
Senior Member
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 128
Rep Power: 18 
It has nothing to do with potential flow. Further short description with publication references can be found on the CFDonline\wiki\source code archive\educational, along with MATLAB code for computing pressureless flow using primitive variables.


March 18, 2015, 15:06 

#16 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
An incompressible, irrotational, inviscid flow = potential flow.
Whatever you are talking about should boil down to potential flow. I guess. And potential flow is "actionatadistance" 

March 18, 2015, 15:35 

#17  
Member
Alex
Join Date: Jan 2014
Posts: 54
Rep Power: 12 
Quote:
Consider the equation for the speed of sounds c^2 = dP/drho which describes how a wave (ie a brief density gradient) propagates with a change in pressure. So, provided there is no change in density, it is possible to have a pressure gradient, but in this case c^2 should go to infinity, and the pressure felt in response to a change in density (c^2 drho = dP), should also approach infinity. Of course, if the speed is sound is infinite then you have to have felt this force for t = 0 seconds (or your favorite unit). Im at work so I cant dig into this deeper. 

March 18, 2015, 15:56 

#18 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 
A truly incompressible fluid means that the fluid is actually at rest. So, no pressure or velocity gradients.


March 18, 2015, 15:59 

#19 
Senior Member
Martin Hegedus
Join Date: Feb 2011
Posts: 500
Rep Power: 19 

March 18, 2015, 17:03 

#20 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,850
Rep Power: 73 
Well, "incompressible fluid" is a model and as a model has approximations and implications...
What we call pressure has no thermodinamic meaning (p=rho*R*T in the compressible formulation but if we would use rho=constant and T= constant, we would get dp=0). That sounds strange but the momentum equation is mathematically a well posed problem when the divergencefree is enforced. What we call "pressure" comes from the continuity equation and is simply a scalar function having gradients such to ensure continuity. We should simply denote as "phi" function and it contains any scalar function acting as gradient. It can be shown that it acts as a lagrangian multiplier. As the BC.s are concerned, the pressure problem is well posed with one boundary condition, that is the projection of the Hodge decomposition along the normal direction to the boundary. That ensures to fulfill the compatibility condition and have a solution for the problem 

Tags 
actionatadistance, body force, incompressible fluid, physics 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
how to set periodic boundary conditions  Ganesh  FLUENT  15  November 18, 2020 07:09 
Difficulty in calculating angular velocity of Savonius turbine simulation  alfaruk  CFX  14  March 17, 2017 07:08 
derivation for special case for incompressible fluid; condition for v?  Simonee  Main CFD Forum  0  March 18, 2014 07:56 
How to apply negtive pressure to outlet  bioman66  CFX  5  June 3, 2006 02:40 
Modeling of free surface in CFD ?  Kim TaeMin  Main CFD Forum  18  July 16, 2001 12:38 