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

Differents results between OPF2.1.2 and OPF2.30

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 16, 2015, 05:29
Default Differents results between OPF2.1.2 and OPF2.30
  #1
Member
 
ynos
Join Date: Jan 2011
Posts: 39
Rep Power: 15
ynos is on a distinguished road
Hi Foamers
I have got differents results when I run the same model with OP211 and OP230, the only one file that I change is fvSchemes.

OPF211
divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV grad(U);
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

OPF230
divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwindV grad(U);
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

Can anyone help me?
Thanks in advanced
ynos is offline   Reply With Quote

Old   March 16, 2015, 05:58
Default
  #2
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

If you really like to get an answer it is better not to increase the number of posts with the same content but

1. Provide more detail on your case (from "bounded" keyword addition I can conclude that you are using steady state solver, but that is all I know about your case and it is not enough).

2. Explain real meaning of "I have got differents results".
alexeym is offline   Reply With Quote

Old   March 16, 2015, 06:35
Default
  #3
Member
 
ynos
Join Date: Jan 2011
Posts: 39
Rep Power: 15
ynos is on a distinguished road
Hi there
I meshed a geometry with OPF2.3.0. And then run this mesh with OPF2.1.1 and 2.3.0. Everyfiles are the same except the fvSchemes, you can see the diferences down.
I get different results in the integration of forces ,on the order of 10%


OPF211
divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV grad(U);
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

OPF230
divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwindV grad(U);
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

Can anyone help me?
Thanks in advanced
ynos is offline   Reply With Quote

Old   March 16, 2015, 15:06
Default
  #4
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 17
pablodecastillo is on a distinguished road
I did the same test with streambody , with 0 angle of atack and i got 7.3 % diference in pressure drag between 2.1.1 and 2.3.1 versions. Same Mesh, same fvsolution , fvschemes .... just adding bounded in 2.3.1 for divSchemes.

Coments are welcome.
pablodecastillo is offline   Reply With Quote

Old   March 17, 2015, 06:06
Default
  #5
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
I guess Henry Weller pointed out this "not-an-issue" thing in my bug report that you can find here:

http://www.cfd-online.com/Forums/ope...-versions.html
Phicau is offline   Reply With Quote

Old   March 17, 2015, 08:00
Default
  #6
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 17
pablodecastillo is on a distinguished road
Following the code linearUpwindV.C, has changes between 2.1.x and 2.3.x, i will back to the old version and recompile, i will inform if it is improving.
pablodecastillo is offline   Reply With Quote

Old   March 18, 2015, 14:47
Default
  #7
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 17
pablodecastillo is on a distinguished road
Still similar diferences, lastest times i meshed with 2.3.x and run with 2.1.1. I decided to move run to the 2.3.x also, but results are too diferents.

Really, i can not understand if other user and developers never got same problems.

At the end editing linearupwindV.C and coment lines from 137 to 216 (both included), i got really similar results. Of course you must recompile the obj.

Last edited by pablodecastillo; March 19, 2015 at 13:07.
pablodecastillo is offline   Reply With Quote

Reply


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



All times are GMT -4. The time now is 08:36.