
[Sponsors] 
November 1, 2013, 17:40 
incompressible limit for compressible codes

#1 
New Member
Join Date: Nov 2013
Posts: 4
Rep Power: 4 
People have always been talking about the incompressible limit for compressible codes.
It seems to me that there are some issues there, but it is not very clear. What would happen if one runs a compressible code, using pressure and velocity as variables, at very low mach number? What would they encounter in numerics? Is there anything diverging or blowing up? If possible, is there any journal paper that discuss this issue? 

November 1, 2013, 17:49 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,311
Rep Power: 28 
Quote:
Think about dp/drho, for very low Mach number, even very small errors in the density solution will be amplified in the pressure field... Generally, preconditioning can be used 

November 1, 2013, 17:55 

#3  
New Member
Join Date: Nov 2013
Posts: 4
Rep Power: 4 
Quote:
So I proposed to solve, for example isothermal NS, in terms of pressure and velocity. Hence, we may avoid solving the density. Is there any fundamental trouble there? 

November 1, 2013, 18:02 

#4  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,311
Rep Power: 28 
Quote:
To solve in terms of velocity and pressure alone you must enforce the constraint Div V = 0. 

November 1, 2013, 18:14 

#5 
New Member
Join Date: Nov 2013
Posts: 4
Rep Power: 4 
It is not necessary to enforce div(u) = 0.
Considering an isothermal flow, perfect gas law gives that pressure is a linear function of density. Hence, you may write density as a function of pressure: rho = rho(p). In doing so, all density terms could be replaced by pressure and its derivatives. For example, the mass balance equation: rho,t + \nabal(rho u) = 0 can be written as drho/dp p_,t + \nabla( rho(p) u) = 0. As incompressible limit, rho becomes a constant, hence the drho/dp becomes zero. It does not blow up, and is still welldefined. In my understanding, if there is any unstability, the source could only come from the locking phenomena, same as the Stokes problem. 

November 1, 2013, 18:28 

#6  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,311
Rep Power: 28 
Quote:
Sorry but I still see some thing I am not convinced about... if you are simoultaneously using constant temperature and density assumptions and also the gas law relation, then pressure must be also constant. The density equation can be rewritten in terms of the pressure field by using dp = a^2 drho but this is valid for isoentropic flows 

November 1, 2013, 18:37 

#7  
New Member
Join Date: Nov 2013
Posts: 4
Rep Power: 4 
Quote:
OK. This is the part not clear to me as well. In textbooks, it says the pressure is no more thermodynamic pressure, but becomes pure mechanical. Some math book says here pressure becomes a Lagrangian multiplier. Maybe someone can help me understanding what is mechanical pressure. In anyway, I think there is no problem in solving the comressible NS equation in terms of pressure. 

November 1, 2013, 18:47 

#8  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,311
Rep Power: 28 
Quote:
yes, in the incompressible (omothermal) limit, pressure is not thermodynamic exactly because p=rho*R*T would simply drive to a constant pressure. Of course that does not make sense in a pressuredriven flow. The system degenerates in the momentum quantity where a gradient of a scalar function appears, supplied by the cinematic constrain Div V=0. Therefore, any scalar function producing a gradient in the momentum such that the velocity field is divergencefree ensures a solution. You can simply see that there is no thermodynamic meaning of this function. It can be shown it is a Lagrangian multiplier, see for example the book of Peric and Ferziger. Furthermore, you can solve the compressible form of the equations by using a suitable preconditioner, for example http://www.grc.nasa.gov/WWW/5810/rvc...nditioning.pdf 

November 1, 2013, 23:58 

#9 
Super Moderator

In the limit of zero mach number, solutions of compressible equations converge to solution of incompressible equations. This is not necessarily true for numerical solutions. For an excellent discussion see
http://citeseerx.ist.psu.edu/viewdoc...=rep1&type=pdf http://hal.archivesouvertes.fr/docs...DF/RR4189.pdf
__________________
http://twitter.com/cfdlab 

Tags 
compreesible codes, low mach flow 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Incompressible turbulence models: strange implementations?  AleDR  OpenFOAM  3  November 18, 2014 12:23 
Is it worth it?  Jason Bardis  Main CFD Forum  47  July 27, 2011 04:52 
Incompressible codes  snegan  Main CFD Forum  0  January 2, 2006 13:51 
Comparison of CFD Codes  Kerem  Main CFD Forum  9  May 9, 2003 04:29 
New List of Free CFD Codes  Bert Laney  Main CFD Forum  5  September 15, 1999 15:24 