[RodriguezFatz]

For the mapped case:
You set the velocity to some fixed value at the outlet. Of course you will not get the same as in cyclic. Also the pressure boundary condition is wrong.
So you should really try to use the boundary conditions for mapped case that I used in my old thread. Why did you change them?

[nusivares]

Hi,

It has been quite a mad week, so I struggled to reply.

I tried implementing conditions your suggested and I did not achieve much success for Re=200 case, which should be low.

It struggled with the covergence as shown below:

RESIDUALS.png


The case files are again in dropbox:

https://www.dropbox.com/s/sjd6e552fb...mapped.7z?dl=0

[RodriguezFatz]

The files you uploaded converge on my computer.

[nusivares]

Hi,

That did slightly confuse me. What version are you using of OpenFOAM?

I have just downloaded the files myself and got the same behaviour as in the post on Saturday.

[RodriguezFatz]

It says "2.3.x.".

[nusivares]

That is fairly weird - I did load 2.3.1 (normally using 2.4.0) and I still got the same result


Do you have any of the solvers modified?

Could send back the file? I would be interested to see it.

[RodriguezFatz]

I deleted the 1000, 2000, ... directories and typed "simplefoam > log". What did you do?

res.png

[nusivares]

I have just now done exactly the same and got a following:

Screenshot-1.png

This is really getting weird and hardly explanable

[RodriguezFatz]

These (1st. order) settings converge much smoother than yours:

gradSchemes
{
    default         Gauss linear;
    //grad(U)         cellLimited Gauss linear 1;
}

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

laplacianSchemes
{
    default         Gauss linear uncorrected;
}

res.png

[nusivares]

Hi,

Thanks for giving a bit of advice on the which schemes to use - I am still slightly confused how to set/modify these.

It did converge. The result seems a bit different, but again, I need to evaluate some physical quantities.

The documentation on them is slightly interesting sometimes.

[Qihao]

[QUOTE=nusivares;569006]Hello all,

As of OpenFOAM I imported the fluent mesh, created cyclic boundary conditions and used the fvOptions file to set the momentum source terms using mean velocity to make the things flow.

I was simulating a channel flow which inlet and out are cyclic boundary conditions, while I do not know how to give a flow velocity boundary condition at the inlet to let it flow, can you add a flow velocity or discharge at the base of cyclic boundary condition?