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

problem with an upwind div scheme

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

Like Tree1Likes
  • 1 Post By wyldckat

Reply
 
LinkBack Thread Tools Display Modes
Old   May 30, 2013, 12:26
Default problem with an upwind div scheme
  #1
Member
 
Join Date: Mar 2013
Posts: 86
Rep Power: 4
giack is on a distinguished road
Hi to all,I have a problem with an upwind div scheme

I set this scheme:
divSchemes
{
div(rho*phi,U) Gauss linearUpwindV cellLimited Gauss linear 1;
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss interfaceCompression;
}

the simulation is good but in my log file there is the following warming message:
--> FOAM Warning :
From function linearUpwindV(const fvMesh&, const surfaceScalarField& faceFlux, Istream&)
in file interpolation/surfaceInterpolation/schemes/linearUpwind/linearUpwindV.H at line 153
Reading "/home/giacomo/Scrivania/20x20x500UP/processor0/../system/fvSchemes::divSchemes::div(rho*phi,U)" at line 30
unexpected additional entries in stream.
Only the name of the gradient scheme in the 'gradSchemes' dictionary should be specified.

It have influence on simulation?
Another question: I would discretize all U terms with an upwind scheme. So i would choose an upwind scheme for the term grad(U) but I obtain an error message and the simulation don't start...How I can set an upwind scheme also for grad(U)

thank to all
giack is offline   Reply With Quote

Old   June 8, 2013, 17:11
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings giack,

This isn't my field of expertise, but the following line is clearly wrong:
Code:
div(rho*phi,U)  Gauss linearUpwindV cellLimited Gauss linear 1;
If I'm not mistaken, if should be something like:
Code:
div(rho*phi,U)  Gauss linearUpwind grad(U);
at least based on what I've seen in tutorials...

You can find more examples, by running:
Code:
find $FOAM_TUTORIALS -name "fvSchemes" | xargs grep "grad(U);"
If you want to change how "grad(U)" is done, should add the entry to the "gradSchemes" block:
Code:
gradSchemes
{
    default         Gauss linear;
    grad(U)        cellLimited Gauss linear 1;
}
Best regards,
Bruno
sharonyue likes this.
wyldckat is offline   Reply With Quote

Old   June 8, 2013, 17:30
Default
  #3
Member
 
Join Date: Mar 2013
Posts: 86
Rep Power: 4
giack is on a distinguished road
thank you for your relpy.
I already tried your last advice but I have another error and simulation don't start. This is a very strange case.
Moreover I find that with the istruction
div(rho*phi,U) Gauss linearUpwindV cellLimited Gauss linear 1;
OF neglect the command and choose the default value for gradScheme required in the div term. So i use the following line:
div(rho*phi,U) Gauss linearUpwindV default;
I think that there aren't a lot of differences beetwen two gradSchene,it's true?
giack is offline   Reply With Quote

Old   June 21, 2013, 00:25
Default
  #4
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 318
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by giack View Post
thank you for your relpy.
I already tried your last advice but I have another error and simulation don't start. This is a very strange case.
Moreover I find that with the istruction
div(rho*phi,U) Gauss linearUpwindV cellLimited Gauss linear 1;
OF neglect the command and choose the default value for gradScheme required in the div term. So i use the following line:
div(rho*phi,U) Gauss linearUpwindV default;
I think that there aren't a lot of differences beetwen two gradSchene,it's true?

Are you able to solve your problem?
Try the scheme limitedLinear and do let us know about the results.

Tushar@cfd is offline   Reply With Quote

Old   June 21, 2013, 03:26
Default
  #5
Member
 
Join Date: Mar 2013
Posts: 86
Rep Power: 4
giack is on a distinguished road
Hi,
at this time I'm simulating two case: one with limetedLinearV 1 and the other with linearUpwindV Gauss linear. Simulations are not over yet but when they end I post my results.
giack 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
How to determine face velocity for upwind scheme? quarkz Main CFD Forum 13 September 6, 2012 05:28
Accuracy problem of HO schemes on unstructured mesh, HO scheme gives 1st order result gemini Main CFD Forum 12 December 27, 2011 22:01
Use of upwind scheme for interpolation of u/v quarkz Main CFD Forum 6 August 30, 2011 04:10
Second order upwind is not UPwind!!! Far CFX 9 May 31, 2011 08:21
problem about numerical scheme in LES. libin Main CFD Forum 4 July 1, 2004 04:32


All times are GMT -4. The time now is 16:32.