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

rhoCentralFoam linearUpwind

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

Like Tree2Likes
  • 1 Post By vkrastev
  • 1 Post By vkrastev

Reply
 
LinkBack Thread Tools Display Modes
Old   November 14, 2011, 02:42
Default rhoCentralFoam linearUpwind
  #1
New Member
 
Join Date: Nov 2011
Posts: 28
Rep Power: 5
Andy_bm is on a distinguished road
Dear foamers.
I want to try the schemes linearUpwind in interpolationSchemes for rhoCentralFoam solver.
In the begining my fvSchemes file looks so:
Code:
fluxScheme Tadmor;
ddtSchemes
{
default Euler;
}
gradSchemes
{
  default cellLimited leastSquares 1;
}
divSchemes
{
default none;
div(tauMC) Gauss linearUpwind;
div(phi,k) Gauss linearUpwind;
div(phi,omega) Gauss linearUpwind;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
reconstruct(rho) upwind;
reconstruct(U) upwind;
reconstruct(T) upwind;
}
snGradSchemes
{
default corrected;
}
Then i replace interpolationSchemes
Code:
interpolationSchemes
{
default linear;
reconstruct(rho) linearUpwind;
reconstruct(U) linearUpwind;
reconstruct(T) linearUpwind;
}
and OpenFoam gives out an error:

Code:
[3] --> FOAM FATAL IO ERROR: 
[3] Grad scheme not specified

Valid grad schemes are :

8
(
Gauss
cellLimited
cellMDLimited
extendedLeastSquares
faceLimited
faceMDLimited
fourth
leastSquares
)
Then i replace interpolationSchemes to

Code:
interpolationSchemes
{
default linear;
reconstruct(rho) Gauss linearUpwind;
reconstruct(U) Gauss linearUpwind;
reconstruct(T) Gauss linearUpwind;
}
OpenFoam gives out an error:

Code:
[17] --> FOAM FATAL IO ERROR: 
[17] Unknown discretisation scheme Gauss

Valid schemes are :

51
(
Gamma
Gamma01
MUSCL
MUSCL01
Minmod
OSPRE
QUICK
SFCD
SuperBee
UMIST
biLinearFit
blended
clippedLinear
cubic
cubicUpwindFit
downwind
filteredLinear
filteredLinear2
filteredLinear3
fixedBlended
harmonic
limitWith
limitedCubic
limitedCubic01
limitedGamma
limitedLimitedCubic
limitedLimitedLinear
limitedLinear
limitedLinear01
limitedMUSCL
limitedVanLeer
linear
linearFit
linearPureUpwindFit
linearUpwind
localBlended
localMax
localMin
midPoint
outletStabilised
quadraticFit
quadraticLinearFit
quadraticLinearUpwindFit
quadraticUpwindFit
reverseLinear
skewCorrected
upwind
vanAlbada
vanLeer
vanLeer01
weighted
)
How probably to use linearUpwind?

Thanks.
Andy_bm is offline   Reply With Quote

Old   November 14, 2011, 04:46
Default
  #2
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
To use the linearUpwind scheme you have to specify also a gradient scheme. The correct syntax in OF-2.0.x should be (for instance):

div(phi,k) Gauss linearUpwind grad(k);

where grad(k) is the scheme specified for k in the gradSchemes section (in your case, the default cellLimited leastSquares 1 scheme will be selected).
In OF-1.7.x this syntax will not be recognized, because you have to specify explicitly the gradient scheme, like in the following example:

div(phi,k) Gauss linearUpwind cellLimited leastSquares 1;

Finally, just a note about div(tauMC): it is an explicitly calculated divergence term, not a convective flux term, thus the only scheme that makes sense is Gauss linear.

Hope this helps

V.
hua1015 likes this.
vkrastev is offline   Reply With Quote

Old   November 14, 2011, 11:11
Default
  #3
New Member
 
Join Date: Nov 2011
Posts: 28
Rep Power: 5
Andy_bm is on a distinguished road
Thanks for your fast reply,
You tip has really helped me, but i have same questions:
If i use
Code:
interpolationSchemes
{
default linear;
reconstruct(rho) upwind cellLimited leastSquares 1;
reconstruct(U) upwind cellLimited leastSquares 1;
reconstruct(T) upwind cellLimited leastSquares 1;
}
What order of the scheme turns out, the first or the second?
Andy_bm is offline   Reply With Quote

Old   November 14, 2011, 11:30
Default
  #4
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
Quote:
Originally Posted by Andy_bm View Post
Thanks for your fast reply,
You tip has really helped me, but i have same questions:
If i use
Code:
interpolationSchemes
{
default linear;
reconstruct(rho) upwind cellLimited leastSquares 1;
reconstruct(U) upwind cellLimited leastSquares 1;
reconstruct(T) upwind cellLimited leastSquares 1;
}
What order of the scheme turns out, the first or the second?
Actually this syntax makes no sense, because the first order upwind scheme does not need any gradient scheme specification (I think that the code will simply ignore anything you write after "upwind"). Anyway, I advice you to use something more accurate for the reconcstruction schemes, like for instance:

Code:
interpolationSchemes
{
default linear;
reconstruct(rho) Gamma 1;
reconstruct(U) GammaV 1;
reconstruct(T) Gamma 1;
}
Regards

V.
hua1015 likes this.
vkrastev is offline   Reply With Quote

Old   November 14, 2011, 11:43
Default
  #5
New Member
 
Join Date: Nov 2011
Posts: 28
Rep Power: 5
Andy_bm is on a distinguished road
I already tried to use this:
Code:
interpolationSchemes
{
default linear;
reconstruct(rho) Gamma 1;
reconstruct(U) GammaV 1;
reconstruct(T) Gamma 1;
}
but my case has dispersing decision

my case has the decision only
Code:
interpolationSchemes
{
default linear;
reconstruct(rho) upwind;
reconstruct(U) upwind;
reconstruct(T) upwind;
}
And now I try find the scheme which will raise accuracy of the decision.

Thanks.
Andy_bm is offline   Reply With Quote

Old   November 14, 2011, 11:51
Default
  #6
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
Quote:
Originally Posted by Andy_bm View Post
I already tried to use this:
Code:
interpolationSchemes
{
default linear;
reconstruct(rho) Gamma 1;
reconstruct(U) GammaV 1;
reconstruct(T) Gamma 1;
}
but my case has dispersing decision

my case has the decision only
Code:
interpolationSchemes
{
default linear;
reconstruct(rho) upwind;
reconstruct(U) upwind;
reconstruct(T) upwind;
}
And now I try find the scheme which will raise accuracy of the decision.

Thanks.
It depends from what you are looking for: rhoCentralFoam is designed as a transonic shock-capturing solver and the solution scheme of Tadmor and Kurganov is also constructed for this kind of flows, so using a first order scheme in my opinion introduces too much diffusion. I've noticed that sometimes rhoCentralFoam returns an oscillating behavior, but if you set the gradientScheme to a limited option (as for instance the cellLimited leastSquares 1 option) the oscillations are largely smeared out. It's up to you.

V.
vkrastev is offline   Reply With Quote

Old   November 14, 2011, 12:02
Default
  #7
New Member
 
Join Date: Nov 2011
Posts: 28
Rep Power: 5
Andy_bm is on a distinguished road
I solve a problem of the expiration stream (М~2-5) in vacuum what solver can you recommend? I used sonicFoam, but he required too much time for the decision.
Andy_bm is offline   Reply With Quote

Old   November 14, 2011, 12:27
Default
  #8
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
Quote:
Originally Posted by Andy_bm View Post
I solve a problem of the expiration stream (М~2-5) in vacuum what solver can you recommend? I used sonicFoam, but he required too much time for the decision.
rhoCentralFoam is the best OF solver for this kind of flow: just play a bit with the schemes until the (eventual) oscillations will be sufficiently small for your requirements.

V.
vkrastev is offline   Reply With Quote

Old   November 14, 2011, 12:34
Default
  #9
New Member
 
Join Date: Nov 2011
Posts: 28
Rep Power: 5
Andy_bm is on a distinguished road
Thanks a lot, I will try select schemes.
Andy_bm 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
A few simple questions about linearUpwind and limitedLinear chegdan OpenFOAM 12 April 10, 2015 04:54
rhoCentralFoam with totalPressure/totalTemperature at inlet of subsonic channel deepblue17 OpenFOAM Running, Solving & CFD 5 February 11, 2013 03:42
localEuler in rhoCentralFoam praveen OpenFOAM 4 October 21, 2012 02:19
linearUpwind scheme in OpenFOAM 2.0.1 ??? cabul OpenFOAM 8 November 9, 2011 07:57
Always crash when solve a C-D nozzle flow field using rhoCentralFoam hawklion OpenFOAM Running, Solving & CFD 0 March 9, 2011 07:13


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