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

linearUpwind scheme in OpenFOAM 2.0.1 ???

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

Like Tree3Likes
  • 1 Post By cabul
  • 2 Post By idrama

Reply
 
LinkBack Thread Tools Display Modes
Old   September 12, 2011, 04:24
Post linearUpwind scheme in OpenFOAM 2.0.1 ???
  #1
New Member
 
André Laß
Join Date: Mar 2010
Location: Rostock, Germany
Posts: 9
Rep Power: 7
cabul is on a distinguished road
Hello FOAMers,
i just wanted to set up a new simulation an user the following schemes:
Quote:
divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV cellLimited Gauss linear 1;
div(phi,k) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,epsilon) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,R) Gauss linearUpwind cellLimited Gauss linear 1;
div(nonlinearStress) Gauss linearUpwind cellLimited Gauss linear 1;
div(R) Gauss linearUpwind cellLimited Gauss linear 1;
div(phi,nuTilda) Gauss linearUpwind cellLimited Gauss linear 1;
div((nuEff*dev(T(grad(U))))) Gauss linearUpwind cellLimited Gauss linear 1;
}
but i ended up with the following message:

Quote:
--> FOAM FATAL IO ERROR:
Unknown discretisation scheme linearUpwind

Valid schemes are :

41
(
Gamma
MUSCL
Minmod
OSPRE
QUICK
SFCD
SuperBee
UMIST
biLinearFit
blended
clippedLinear
cubic
cubicUpwindFit
downwind
filteredLinear
filteredLinear2
filteredLinear3
fixedBlended
limitWith
limitedCubic
limitedLinear
limiterBlended
linear
linearFit
linearPureUpwindFit
localBlended
localMax
localMin
midPoint
outletStabilised
pointLinear
quadraticFit
quadraticLinearFit
quadraticLinearUpwindFit
quadraticUpwindFit
reverseLinear
skewCorrected
upwind
vanAlbada
vanLeer
weighted
)


file: ./S-0.00/system/fvSchemes::divSchemes::div((nuEff*dev(T(grad(U)))) ) at line 40.

From function surfaceInterpolationScheme<Type>::New(const fvMesh&, Istream&)
in file lnInclude/surfaceInterpolationScheme.C at line 82.

FOAM exiting
Is there a new name for the linearUpwind scheme or has it been removed? I am quite new to OpenFOAM so maybe i made a mistake myself. Can Anyone help me?

Greeting,s
André
mm.abdollahzadeh likes this.
cabul is offline   Reply With Quote

Old   September 13, 2011, 03:59
Default
  #2
New Member
 
André Laß
Join Date: Mar 2010
Location: Rostock, Germany
Posts: 9
Rep Power: 7
cabul is on a distinguished road
Does nobody have the same Problem?
cabul is offline   Reply With Quote

Old   September 13, 2011, 07:00
Exclamation
  #3
New Member
 
André Laß
Join Date: Mar 2010
Location: Rostock, Germany
Posts: 9
Rep Power: 7
cabul is on a distinguished road
Hello there,

i played around with the fvSchemes for a while and ended up with the following settings, that are working for me.

Quote:
gradSchemes
{
default cellLimited Gauss linear 1;
}

divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV default;
div(phi,k) Gauss linearUpwind default;
div(phi,epsilon) Gauss linearUpwind default;
div(phi,omega) Gauss linearUpwind default;
div(R) Gauss linear;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
Can somebody have a look at the divschemes? As far as i understood OpenFOAM 2.0.1 checks for an gradShemes right behind the definition "Gauss linearUpwind" . Is it correct to set the "default" value here, so that it uses the "cellLimited Gauss linear 1" -scheme, defined in gradSchemes?

as mentioned before, the calculation is running fine and the results look good too.

Greetings,
André
cabul is offline   Reply With Quote

Old   September 13, 2011, 10:37
Default
  #4
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 8
idrama is on a distinguished road
In OF-2 I use linearUpwind(V) in this way:

div(phi,U) Gauss linearUpwindV grad(U);
dib(phi,k) Gauss linearUpwind grad(k);

and so on.

Cheers
mm.abdollahzadeh and s.m like this.
idrama is offline   Reply With Quote

Old   November 7, 2011, 14:05
Default
  #5
Member
 
fisch
Join Date: Feb 2010
Posts: 97
Rep Power: 7
fisch is on a distinguished road
Hello,

could someome tell me what the meaning of this grad(U) in the div(phi,U) discretization is?
Thanks

rupert
fisch is offline   Reply With Quote

Old   November 8, 2011, 06:19
Default
  #6
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
Quote:
Originally Posted by fisch View Post
Hello,

could someome tell me what the meaning of this grad(U) in the div(phi,U) discretization is?
Thanks

rupert
It tells you which gradient scheme to use for the gradient part of the LUD discretisation. The specific gradient scheme is defined in the gradSchemes section. So you do not have to use "grad(U)", you can use anything you like as long as it is defined in the gradSchemes section. E.g.

gradSchemes
{
blahblahblah Gauss linear;
}

divSchemes
{
div(phi,U) Gauss linearUpwindV blahblahblah;
}

should work just fine. In all cases things like "div(phi,U)" and "grad(U)" are just identifiers, thy have no special structure or significance outside this.
eugene is offline   Reply With Quote

Old   November 8, 2011, 06:46
Default
  #7
Member
 
fisch
Join Date: Feb 2010
Posts: 97
Rep Power: 7
fisch is on a distinguished road
Thanks for the answer eugene.

now arise two questions:
1.so it should be the same as the gradient in the implicit field in the div-term?
2. what happens then if i dont specify this: just like div(phi,U) Gauss linear?

thanks a lot
fisch is offline   Reply With Quote

Old   November 9, 2011, 06:39
Default
  #8
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
1. As far as I know there is no fixed rule about the nature of the gradient portion of a LUD scheme. Its not like pressure and velocity flux interpolation that has to be consistent.
2. If you don't specify anything at all, the code should crash (or maybe not, I haven't checked). If you specify something with no match in the gradScemes section, it should use the "default" gradScheme entry without complaint. If no default entry is present, it should exit with an error.
eugene is offline   Reply With Quote

Old   November 9, 2011, 07:57
Default
  #9
Member
 
fisch
Join Date: Feb 2010
Posts: 97
Rep Power: 7
fisch is on a distinguished road
thank you eugene
fisch 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
Second order upwinding scheme in OpenFoam subash OpenFOAM Running, Solving & CFD 4 June 20, 2012 18:28
OpenFOAM 1.6.x, 1.7.0 and 1.7.x are not fully prepared to work with gcc-4.5.x wyldckat OpenFOAM Bugs 18 October 21, 2010 05:51
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 3 September 8, 2010 06:25
Parameter for the cubic upwind scheme in OpenFOAM JinBiao OpenFOAM 1 May 23, 2010 03:51
Different flow pattern between OpenFOAM and CFX AirS OpenFOAM 0 January 12, 2010 08:08


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