
[Sponsors] 
April 28, 2010, 21:39 
Pressure Gradient and Velocity Divergence

#1 
New Member
AlvinXu
Join Date: Nov 2009
Posts: 15
Rep Power: 9 
the approximation of the pressure gradient and the divergence of the velocity must be compatible if the kinetic energy conservation is to hold.
Once either approximation is chosen, the freedom to choose the other is lost. I want to know the meaning of this word " compatible" 

May 18, 2010, 10:10 

#2 
New Member
Join Date: May 2010
Posts: 8
Rep Power: 8 
I suppose you are referring to the NavierStokes equations for an uncompressible fluid.
In that case, at least for FEM, there is a well established theory that explains when two discretizations are compatible. You should look for LBB or infsup condition. It should be explained in any recent text books on uncompressible CFD. 

May 25, 2010, 17:11 

#3 
New Member
Join Date: May 2009
Posts: 28
Rep Power: 9 
You can read these papers...
http://www.imati.cnr.it/brezzi/paper...tokesCIME.pdf http://ima.umn.edu/~arnold/papers/mixed.pdf There is also a lecture of Gilbert Strang about the LBB condition at the MIT section of iTunes U. 

May 29, 2010, 00:07 

#4 
Senior Member
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 108
Rep Power: 10 
The frequentlyheard statement which Alvin Xu refers to is incorrect as it stands. However, if the word "pressure" is replaced, a very similar statement is true. The resolution is, as one author states, "not everything denoted by the symbol p is a pressure". Further, the statement is restricted to a particular class of projection methods.
The governing equation (NavierStokes) for incompressible flow is a restatement of Newton's law for conservation of momentum, perhaps after some simplifying assumptions about the fluid density. With those incompressibility assumptions, the velocity field is divergencefree. Then the governing equation begins dv/dt=... . Since the velocity is divergencefree, so the right hand side must be. Now, according to Helmholtz and his successors, a vector field can be orthogonally decomposed into divergencefree (solenoidal) and irrotational parts, subject to appropriate boundary conditions. One can see immediately that solutions of Laplace's equation are both divergencefree and irrotational. One can embrace this ambiguity or eliminate it through boundary conditions. A term like grad p is patently irrotational, and cannot appear on the RHS. That is  the fluid flow is independent of the pressure. Such a statement will not surprise those who use stream functionvorticity formulation (which does not contain a pressure) to compute fluid flow. But to be correct, one must project out the irrotational part of the RHS. There are several ways of doing this. Doing this projection by the usual mathematical analysis leads to an integrodifferential equation. I don't know anyone who has used this approach for practical computation. The most common projection method is to subtract off an irrotational function in the form of the gradient of a potential to cancel the errant terms on the RHS. When one does this, the resulting equation looks just like the original form except that the subtracted term is not the pressure. It looks like a pressure, but has different boundary conditions, and the functional form must satisfy the LBB/infsup conditions. Once one has computed an approximate (weak) divergencefree velocity, one can compute the pressure with the "pressure Poison" equation (as the stream functionvorticity practitioners do) using whatever functional form you wish, with appropriate b.c. for the pressure (and no LBB). A third approach (which I advocate) is to expand the velocity in terms of an appropriately conforming divergencefree basis, i.e. using divergencefree finite elements. These velocity elements are the curl of a Hermitetype stream function element in 2D. Such elements can be found by modifying one of the many Hermite platebending finite elements found in the literature. Projection is accomplished by the orthogonality of solenoidal and irrotational functions a la Helmholtz when computing the element matricies. Much of this carries over to 3D, but there is no shortcut to finding satisfactory elements as in 2D. Note that there is always a scalar function that appears: the projection gradient, the 2D vorticity, or the stream function. 

May 29, 2010, 11:11 

#5  
New Member
vahid velayati
Join Date: May 2010
Posts: 1
Rep Power: 0 
Quote:
hi i am looking for velicity profile over naca0012 airfoil at incompressible mach number. please send me any document that may help me. thanks 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
velocity gradient issue in cfx  rystokes  CFX  3  November 29, 2010 01:43 
Physical meaning of pressure in pressure eqn  dontknow  Main CFD Forum  2  August 30, 2010 10:12 
Pressure Gradient in channelOodles  sega  OpenFOAM Running, Solving & CFD  14  August 12, 2009 06:21 
How to calculate gradient of pressure and velocity on an outlet face  skabilan  OpenFOAM Running, Solving & CFD  2  December 8, 2008 18:26 
order of magnitude analysis  atit koonsrisuk  Main CFD Forum  3  July 27, 2000 11:59 