CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Compressible Flow Solvers

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 28, 2005, 16:33
Default Hello In compressible liqu
  #1
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
Hello

In compressible liquid/gas flows, the pEqn reads :

(
fvm::ddt(psi,p)+fvc::div(phi)+fvm::div(phid,p,".." )-fvm::laplacian(..,p)
)

where as in an incompressible code pEqn is given as,

(
fvm::laplacian(1/A, p)==fvc::div(phi)
)

I tried to look into pressure correction methodologies for compressible flows, but was not able to figure out the reason of this implementation procedure (why the need for pressre time derivative) . Can someone help me out with this ?



Thanks



Vatant
vatant is offline   Reply With Quote

Old   April 28, 2005, 16:39
Default > why the need for pressre tim
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
> why the need for pressre time derivative

Because in compressible flow there is a rate of change of density and hence pressure, or vice versa depending on how you look at it.
henry is offline   Reply With Quote

Old   April 28, 2005, 16:40
Default Is the implementation same as
  #3
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
Is the implementation same as artificial compressibility method ?
vatant is offline   Reply With Quote

Old   April 28, 2005, 16:43
Default Certainly not, it's real compr
  #4
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
Certainly not, it's real compressibility, there is nothing artificial about the compressibility in the compressible codes in OpenFOAM.
henry is offline   Reply With Quote

Old   April 28, 2005, 16:51
Default does the term 'dp/dt' vanish a
  #5
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
does the term 'dp/dt' vanish after reaching certain tolerance level ? what is this method for solving compressible flow called?
vatant is offline   Reply With Quote

Old   April 28, 2005, 16:56
Default For steady-state flows all tim
  #6
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
For steady-state flows all time derivatives will vanish but not for transient flows. If you are solving steady-state flows you might be better of using one of the steady-state solvers.
henry is offline   Reply With Quote

Old   April 28, 2005, 17:10
Default I am solving transient flows a
  #7
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
I am solving transient flows and was wondering if rho and p relations are not explicity given (unlike ideal gas relation), how does this pressure correction method works?

Is the pressure correction equation derived from combination of continuity and momentum equation ?

If i want to search in literature about this method, could you tell me whats this compressible solver schemes called?


Vatant
vatant is offline   Reply With Quote

Old   April 28, 2005, 17:17
Default > Is the pressure correction e
  #8
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
> Is the pressure correction equation derived from combination of continuity and momentum equation ?

Yes in the same manner as for incompressible PISO. I guess it would be called compressible PISO but I am not sure I didn't implement it from a paper I derived it from first principles and implemented it. However I am sure other people have implemented the same and named and published it.
henry is offline   Reply With Quote

Old   April 28, 2005, 17:25
Default I've got an old report on this
  #9
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,907
Rep Power: 33
hjasak will become famous soon enough
I've got an old report on this which contains the derivation of the pressure equation for compressible flow - send me an E-mail if you want it.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   April 28, 2005, 17:38
Default For the analysis of compressib
  #10
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
For the analysis of compressible flows, it is possible that some regions of the flow domain
such as in the boundary layers have low speeds and thus are incompressible. This low velocities result in numerical disorders..due to eigenvalue mismatch..

So, if we have an first principle formulation, would we require any kind of preconditioning for handling low speed regions? or by any means numerical stability is inherently guaranteed...



Vatant
vatant is offline   Reply With Quote

Old   April 28, 2005, 17:46
Default Eliptic pressure equation base
  #11
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
Eliptic pressure equation based compressible flow solvers do not suffer from problems handling low speed regions like density based solvers which is why we choose them. We are able to solve subsonic, transonic and supersonic flows with or without boundary layers without difficulty.
henry is offline   Reply With Quote

Old   April 28, 2005, 18:08
Default why is that the convergence pr
  #12
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
why is that the convergence problem be easily handled with elliptic systems ?

The derived pEqn from continuity and momentum eqn was of the form d/dt(p) + laplacian(..,p)+div(phi)...=0 , is it a completely elliptic version (with the time derivative terms) ?




Vatant
vatant is offline   Reply With Quote

Old   April 28, 2005, 18:13
Default > why is that the convergence
  #13
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
> why is that the convergence problem be easily handled with elliptic systems ?

This is well known and explained in books on the subject.

The pressure equation may take several forms depending on how the compressibility effects are handled. At least it is a Poisson's equation but it may also include a "convection" term as in sonicFoam.
henry is offline   Reply With Quote

Old   April 28, 2005, 18:25
Default Incompressible flows would rep
  #14
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
Incompressible flows would represent elliptic, solving laplacian (p) == div (phi)..

In the sonicFoam, the pEqn has the extra time derivative term with compressiblity effect,

ddt(psi,p)+ div (phid,p,..)-laplacian(rho/A,p)=0

could the equation still be called a poisson ?

Since, the time derivative appears (with compressiblity), im not sure how this works.



Vatant
vatant is offline   Reply With Quote

Old   April 28, 2005, 18:35
Default > could the equation still be
  #15
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
> could the equation still be called a poisson ?

No it isn't a Poisson's equation with the implicit "convection" term but it might be called one if that part is explicit as it is in our new compressible flow solver, it depends on your definition of Poisson's equation. I guess with an explicit "convection" term but implicit compressibility term it would be a linear combination of a Poisson's equation and a Helmholtz equation; I am not sure if that kind of equation has a special name.
henry is offline   Reply With Quote

Old   April 28, 2005, 19:03
Default I am still running openFoam 1.
  #16
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
I am still running openFoam 1.0 which has implicit convection solves..you think i might have any problem in convergence since convection is not done explicit taking it away from elliptic ?

Does a combination of poisson's and helmholtz equation in behavior result close to poisson ? Since that the flow needs show elliptic behavior for better convergence.


Vatant
vatant is offline   Reply With Quote

Old   April 28, 2005, 19:31
Default There are no convergence probl
  #17
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
There are no convergence problems with the compressible codes in OpenFOAM, the equations are what they are, the pressure equation is implemented implicitly and it doesn't matter what the name of this type of equation is, it is what it is.
henry is offline   Reply With Quote

Old   April 28, 2005, 19:44
Default I actually ran lot of cases wi
  #18
Member
 
Join Date: Mar 2009
Posts: 43
Rep Power: 17
vatant is on a distinguished road
I actually ran lot of cases with OpenFoam compressible codes with different flow velocities and the convergence was good . I did experience some problems in injector flows with sharp corners...
Hence I was interested in the numerical essence of the implementation.

With good convergence available with OpenFoam, I shall study some acoustics due to compressiblity effects at high speed flows with low speed regions in the domain.


Vatant
vatant is offline   Reply With Quote

Old   June 8, 2005, 15:23
Default Hi Guys I would like an opi
  #19
Member
 
Wladimyr Mattos da Costa Dourado
Join Date: Mar 2009
Location: Sao Jose dos Campos, SP, Brazil
Posts: 36
Rep Power: 17
mattos is on a distinguished road
Send a message via Skype™ to mattos
Hi Guys

I would like an opinion about the following subject. We would like to implement a pseudo-compressibility solver or some kind of preconditioned method based solver in the OpenFoam. In my knowledge, the non segregated methods is one of the best way to do that. Is possible to implement some non-segragated method in the openfoam? How hard work is necessary to do it? In other words, what is necessary to implement a non-segragated method in openfoam? Are the solvers implemented in openfoam able to solve the matrix generated in a non-segragated method?

About wiki-openfoam: Which is better to use for these type of discussions: this site or wiki site? Is wiki site as up to dated such much this site?

Many Tanks in advance!

Wladimyr
mattos is offline   Reply With Quote

Old   June 8, 2005, 15:37
Default The subject of segregated vs p
  #20
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
The subject of segregated vs point-coupled vs block-coupled compressible solvers has already been discussed on this site, you might find it useful to search for and read those threads.
henry is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Yplus in compressible solvers dinonettis OpenFOAM Post-Processing 7 October 6, 2011 06:51
Incompressible vs compressible solvers James Main CFD Forum 5 January 19, 2009 06:15
NS-incompressible and compressible flow solvers ag Main CFD Forum 2 September 27, 2005 07:18
External Flow-compressible flow solver-lift/drag Tom Brown Main CFD Forum 7 December 29, 2000 14:41
MHD flow solvers Manu Main CFD Forum 2 October 31, 2000 03:11


All times are GMT -4. The time now is 13:37.