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

Difference SFCD and GammaV

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

Like Tree5Likes
  • 1 Post By bboonacker
  • 3 Post By hjasak
  • 1 Post By hjasak

Reply
 
LinkBack Thread Tools Display Modes
Old   October 23, 2008, 05:52
Default Hello All and especially Prof.
  #1
New Member
 
Bart Boonacker
Join Date: Mar 2009
Posts: 7
Rep Power: 8
bboonacker is on a distinguished road
Hello All and especially Prof. Jasak,

I am currently working on my MSc thesis and I would like to know something;
I cannot find any reference on the SFCD scheme and how it works. I can try to find out using the .C and .H, but my C++ knowledge does not extend that far I regret.
1) Does anyone have a reference or is anyone willing to share some thoughts/opinions about SFCD ?
2) Also, what is the main difference between SFCD and the Gamma scheme ? Please correct me if I'm wrong, Is the difference that the Gamma scheme contains a scaling factor (which is 2*beta_m in Jasaks paper(?)) and SFCD does not, or is it also in the blending function or other parts?

Best regards,
Bart
immortality likes this.
bboonacker is offline   Reply With Quote

Old   October 23, 2008, 07:28
Default Hello Bart, All these schem
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,763
Rep Power: 21
hjasak will become famous soon enough
Hello Bart,

All these schemes are blended upwind-central type, using a limiter. The only true difference is the form of the limiter and its basis.

The original reference for SFCD is a PhD from Dr. Harry Ziman, from mid 1990-s at Imperial College London. Basically, the limited looks like MINMOD and you can revrese-engineer it from the code:


scalar phict = LimiterFunc::phict
(
faceFlux, phiP, phiN, gradcP, gradcN, d
);

scalar limitPhict = min(max(phict, 0), 0.5);
return limitPhict/(1 - limitPhict);


Here, phict is Phi-C-tilde (have a look at the Gamma paper) and the limiter is written out. For Gamma, the limiter is:

min(max(phict/k_, 0), 1) but k is re-scaled on the constructor to be between 0 and 1 (in reality is is 0-0.5).

k_ = max(k_/2.0, SMALL);

Enjoy,

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

Old   October 23, 2008, 09:08
Default Thanks for the quick clarifyin
  #3
New Member
 
Bart Boonacker
Join Date: Mar 2009
Posts: 7
Rep Power: 8
bboonacker is on a distinguished road
Thanks for the quick clarifying awnser, I'll see if I can get a hold off Dr. Ziman's paper.

What I don't get so far, k is equal to beta_m ? Then why is Phi-C-tilde limited ?
Is the limiting done to easily select between UD, blending or CD, or something a like ?
bboonacker is offline   Reply With Quote

Old   October 23, 2008, 09:36
Default @PhdThesis{Ziman:PhD, autho
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,763
Rep Power: 21
hjasak will become famous soon enough
@PhdThesis{Ziman:PhD,
author = {Ziman, H.J.},
title = {A computer prediction for chemically reacting flows
in stirred tanks},
school = {Imperial College, University of London},
year = {1990}
}

Yes, it's all about automaticaly blending UD and CD (well, a second-order scheme of your choice) .

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

Old   October 24, 2008, 09:12
Default Sorry, I do not do support wit
  #5
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,763
Rep Power: 21
hjasak will become famous soon enough
Sorry, I do not do support without a support contract. You simulation failed way before this point.

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

Old   December 2, 2010, 05:10
Default
  #6
Senior Member
 
Dr. Alexander Vakhrushev
Join Date: Mar 2009
Posts: 213
Rep Power: 10
makaveli_lcf is on a distinguished road
Send a message via ICQ to makaveli_lcf
Is SFCD a NVD scheme? Here it is states as a "General convection scheme"
How is it possible to understand from the OpenFOAM sources, is scheme TVD or NVD? I found only the directory structure:

Quote:

Directories


directory limitedSchemes
directory multivariateSchemes
directory schemes
directory surfaceInterpolation
directory surfaceInterpolationScheme
Regards,

Alexander
__________________
Best regards,

Dr. Alexander VAKHRUSHEV

Christian Doppler Laboratory for "Advanced Process Simulation of
Solidification and Melting"

Simulation and Modelling of Metallurgical Processes
Department of Metallurgy
University of Leoben

Franz-Josef-Str. 18
A - 8700 Leoben
Österreich / Austria
Tel.: +43 3842 - 402 - 3125
http://smmp.unileoben.ac.at
makaveli_lcf is offline   Reply With Quote

Old   March 19, 2014, 18:07
Default
  #7
Senior Member
 
Join Date: Jan 2013
Posts: 204
Rep Power: 5
openfoammaofnepo is on a distinguished road
Dear Dr. Jasak,

Your reply is very helpful for me to understand the Gamma scheme. However, can I ask you other several questions?

1, what does k_ stand for in the code you mentioned?

2, about the limiter function for gamma scheme, did you have any reference for it? I did not find any information about it from the literatures?

3, this questions is not linked to Gamma scheme, but I ask here:
in openfoam, the numerical fluxes on the inter-cell face are for primitive variables (e.g. U_f) or conservative variables (e.g. (rho*U)_f)?

Thank you so much!

OFFO

Quote:
Originally Posted by hjasak View Post
Hello Bart,

All these schemes are blended upwind-central type, using a limiter. The only true difference is the form of the limiter and its basis.

The original reference for SFCD is a PhD from Dr. Harry Ziman, from mid 1990-s at Imperial College London. Basically, the limited looks like MINMOD and you can revrese-engineer it from the code:


scalar phict = LimiterFunc:hict
(
faceFlux, phiP, phiN, gradcP, gradcN, d
);

scalar limitPhict = min(max(phict, 0), 0.5);
return limitPhict/(1 - limitPhict);


Here, phict is Phi-C-tilde (have a look at the Gamma paper) and the limiter is written out. For Gamma, the limiter is:

min(max(phict/k_, 0), 1) but k is re-scaled on the constructor to be between 0 and 1 (in reality is is 0-0.5).

k_ = max(k_/2.0, SMALL);

Enjoy,

Hrv
openfoammaofnepo is offline   Reply With Quote

Old   March 20, 2014, 12:51
Default
  #8
Senior Member
 
Join Date: Jan 2013
Posts: 204
Rep Power: 5
openfoammaofnepo is on a distinguished road
Dear Prof. Jasak,

I have some understanding about the limiter for Gamma function used in Openfoam. But I found there is still some difference compared with your Gamma Scheme paper. The scalar k_ should be the constant beta_m in that paper. The difference is:

Code:
In the paper:
if 0<phict/k_<1, then blended upwind and central differencing 
if phict/k_ <= 0, then upwind
if  phict/k_ >= 1, then central difference

In OF2.1.1:
if 0<phict/k_<1, then blended upwind and central differencing 
if phict/k_ <= 0 .or. phict/k_ >=1/K-, then upwind
if  1/k_>= phict/k_ >= 1, then central difference
Is what I am saying correct? Thank you for your any comments.

OFFO
openfoammaofnepo is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SFCD Scheme sampaio OpenFOAM Running, Solving & CFD 19 November 25, 2008 13:22
SFCD in tetrahedral meshes cosimobianchini OpenFOAM Running, Solving & CFD 2 March 6, 2007 12:21
SFCD turb Main CFD Forum 0 March 16, 2005 19:25
Difference between 2D and 3D ! Christophe FLUENT 2 January 26, 2005 04:10
what is the difference cfd student Main CFD Forum 0 October 7, 2004 16:33


All times are GMT -4. The time now is 01:21.