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

fvScheme and NonOrthogonal Correctors

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

Like Tree13Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 11, 2011, 11:18
Default
  #21
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17
Andrea_85 is on a distinguished road
Hi V.

where i have to put the reconCentral files? /src/finiteVolume/....?

Thanks

andrea
Andrea_85 is offline   Reply With Quote

Old   April 11, 2011, 13:10
Default
  #22
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17
Andrea_85 is on a distinguished road
Hi,
another question...

both schemes (reconCentral and pointLinear) are interpolation scheme, if i'm not wrong. So if i use one of these for the gradient Should I use the same for interpolation?
like

gradSchemes ->Gauss pointLinear;
interpolationSchemes ->pointLinear;

or, with reconCentral:

gradSchemes -> Gauss reconCentral;
interpolation -> reconCentral;

Thanks again

andrea
Andrea_85 is offline   Reply With Quote

Old   April 11, 2011, 13:10
Default
  #23
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
Quote:
Originally Posted by Andrea_85 View Post
Hi V.

where i have to put the reconCentral files? /src/finiteVolume/....?

Thanks

andrea

If my memory does not fail, you can compile it dynamically by typing wmake libso inside the source folder (wherever you choose to place it). After that you have to set properly your fvSchemes dictionary (search for some examples in the forum) and include in your controlDict something like this:

libs ( "libOpenFOAM.so" "reconCentral.so" )
Hope this helps

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

Old   April 11, 2011, 13:23
Default
  #24
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
Quote:
Originally Posted by Andrea_85 View Post
Hi,
another question...

both schemes (reconCentral and pointLinear) are interpolation scheme, if i'm not wrong. So if i use one of these for the gradient Should I use the same for interpolation?
like

gradSchemes ->Gauss pointLinear;
interpolationSchemes ->pointLinear;

or, with reconCentral:

gradSchemes -> Gauss reconCentral;
interpolation -> reconCentral;

Thanks again

andrea
You have to set something like this (of course cellLimited leastSquares 1 is only one of the possible choices, you can also try for instance cellLimited Gauss linear 1):

divSchemes
{
div(phi,U) Gauss reconCentral cellLimited leastSquares 1.0;
//...
}

interpolationSchemes
{
default linear;
interpolate(HbyA) linear;
interpolate(U) reconCentral phi cellLimited leastSquares 1.0;
}

Hope this helps too

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

Old   April 12, 2011, 04:27
Default
  #25
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17
Andrea_85 is on a distinguished road
Thanks a lot V., you have been very helpfull!!

andrea
Andrea_85 is offline   Reply With Quote

Old   April 12, 2011, 12:58
Default
  #26
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17
Andrea_85 is on a distinguished road
Hi all, just to let you know...
I finally was able to find a good solution with the unstructured mesh. The best choice for my cases was the pointLinear scheme suggested by Markus. The solution now is a (quite) perfect circle and the behavior during the simulation is very similar to the structured case.
The drop is not perfectly still in the middle of the domain but is slightly moving in a sort of circle (probably due to the parasite current) and the solution varies slightly around the analytical one, but for the moment i'm quite happy.

Hope this will help someone else!


andrea
Andrea_85 is offline   Reply With Quote

Old   April 14, 2011, 11:01
Default
  #27
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 184
Rep Power: 17
markusrehm is on a distinguished road
Hello Andrea,

could you post your fvSchemes for the best case?

Thank you.
markusrehm is offline   Reply With Quote

Old   April 14, 2011, 12:35
Default
  #28
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17
Andrea_85 is on a distinguished road
Hi Markus,

i tried different setting (always using the pointLinear schemes) and i found good results with almost all. here two of these:
1)

Code:
ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         cellMDLimited Gauss pointLinear 0.333;
}

divSchemes
{
    div(rho*phi,U)  Gauss pointLinear 0.333;
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
}

laplacianSchemes
{
    default         Gauss linear limited 0.333;
}

interpolationSchemes
{
    default         pointLinear;
}

snGradSchemes
{
    default         limited 0.333;
}

fluxRequired
{
    default         no;
    p_rgh;
    pcorr;
    alpha1;
}
2)
Code:
ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         cellMDLimited Gauss pointLinear 0.333;
}

divSchemes
{
    div(rho*phi,U)  Gauss limitedLinearV 0.333;
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
}

laplacianSchemes
{
    default         Gauss linear limited 0.333;
}

interpolationSchemes
{
    default         pointLinear;
}

snGradSchemes
{
    default         limited 0.333;
}

fluxRequired
{
    default         no;
    p_rgh;
    pcorr;
    alpha1;
}
i also found good results without cellMDLimited on the gradient and with limit 1 (corrected).


hope this help

andrea
fumiya, hua1015, SailorLiu and 3 others like this.
Andrea_85 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On



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