[Andrea_85]

Hi V.

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

Thanks

andrea

[Andrea_85]

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

[vkrastev]

Originally Posted by Andrea_85

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]

Originally Posted by Andrea_85

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.

[Andrea_85]

Thanks a lot V., you have been very helpfull!!

andrea

[Andrea_85]

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

[markusrehm]

Hello Andrea,

could you post your fvSchemes for the best case?

Thank you.

[Andrea_85]

Hi Markus,

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

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)

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