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

fvScheme and NonOrthogonal Correctors

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

Like Tree7Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   April 11, 2011, 10:18
Default
  #21
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
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, 12:10
Default
  #22
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
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, 12:10
Default
  #23
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 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.
vkrastev is offline   Reply With Quote

Old   April 11, 2011, 12:23
Default
  #24
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 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, 03:27
Default
  #25
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
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, 11:58
Default
  #26
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
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, 10:01
Default
  #27
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 176
Rep Power: 8
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, 11:35
Default
  #28
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 275
Rep Power: 7
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 and SailorLiu like this.
Andrea_85 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



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