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

Fundamental problem with Gamma scheme

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

Like Tree5Likes
  • 5 Post By hjasak

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 28, 2008, 05:08
Default Hello, do I understand the
  #1
New Member
 
Dominik Christ
Join Date: Mar 2009
Posts: 28
Rep Power: 17
dominik_christ is on a distinguished road
Hello,

do I understand the Gamma discretisation scheme right that it is a blend of central differencing (CD) and upwind differencing (UD)? But then I don't understand the following problem:
My test case converges easily with both UD and CD (also LUD runs without problems) giving me residuals in the order of 1e-07.
If I switch to Gamma scheme, I don't get convergence, the residuals remain around 0.4 and every step takes a lot of internal iterations.

My very simple test case (I attached it as a tgz-file) uses scalarTransportFoam on a 10x10 checkerboard mesh. I want to use it to evaluate convection discretisation schemes like it is done in the book by Versteeg and Malalasekera.
The uniform velocity field is diagonal from the lower left to the upper right corner. The T field is set to 1 on the left side and 0 on the bottom.
The boundary conditions are:

For T:
internalField uniform 0;

boundaryField
{
left
{
type fixedValue;
value uniform 1;
}
bottom
{
type fixedValue;
value uniform 0;
}
top
{
type zeroGradient;
}
right
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
For U:
internalField uniform (2 2 0);

boundaryField
{
left
{
type fixedValue;
value uniform (2 2 0);
}
bottom
{
type fixedValue;
value uniform (2 2 0);
}
top
{
type zeroGradient;
}
right
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

The value for the diffusivity DT is set to zero in the transportProperties dictionary.

So, do I run into a fundamental problem of the Gamma scheme, or did I just overlook a fault in my setup? Or is it a problem to use steadyState with scalarTransportFoam?

Thanks for your answers in advance!!

Regards
Dominik

checkerBoard.tgz
__________________
__________________
Dominik Christ
Providing commercial foam-extend/OpenFOAM and CFD Consulting: http://wikki.co.uk
dominik_christ is offline   Reply With Quote

Old   August 28, 2008, 05:42
Default In fact, you have not done any
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,906
Rep Power: 33
hjasak will become famous soon enough
In fact, you have not done anything wrong - you have discovered the problem at the basis of CFD convection discretisation research for the last 25 years. Basically, a scheme must be non-linear to be bounded and more than first-order accurate at the same time. With a non-linear scheme, there is no guarantee of convergence. In practice and for most schemes, you will have to play with the adjustable factor in the scheme to get what you want - it is a balance between accuracy and convergence.

Have a look at some convection differencing scheme papers and all will be clear. To make things easier, you may wish to tighten the solver tolerances, decrease the time-step size, use deferred correction schemes or all of the above. BTW, steadyState is definitely a bad idea, because you just lost the diagonal dominance.

Enjoy,

Hrv
Amir, hua1015, 1/153 and 2 others like this.
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   August 28, 2008, 11:28
Default Hi Hrv, thank you very much
  #3
New Member
 
Dominik Christ
Join Date: Mar 2009
Posts: 28
Rep Power: 17
dominik_christ is on a distinguished road
Hi Hrv,

thank you very much for your quick answer!!
Now that I know that the case setup was correct, I can dig into theory...
__________________
__________________
Dominik Christ
Providing commercial foam-extend/OpenFOAM and CFD Consulting: http://wikki.co.uk
dominik_christ is offline   Reply With Quote

Old   August 29, 2008, 08:24
Default Hi Dominik, Sometimes the a
  #4
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 21
eugene is on a distinguished road
Hi Dominik,

Sometimes the act of switching between CD to UD blended can itself introduce instability.
eugene 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
Display regions central linear and upwind for gamma scheme braennstroem OpenFOAM Post-Processing 3 May 28, 2008 03:49
Fundamental Problem Naghman Khan FLUENT 8 August 1, 2007 11:01
A fundamental problem about the UcorrectBoundaryConditions in icoFoam dbxmcf OpenFOAM Running, Solving & CFD 0 February 26, 2007 22:16
?to implement gamma differencing scheme Gonski Main CFD Forum 15 January 28, 2007 17:00
I cannot find the route from fvmdivphi gamma to Gamma201 scheme luckyluke OpenFOAM Running, Solving & CFD 1 June 17, 2005 15:58


All times are GMT -4. The time now is 12:41.