[adrieno]

Hi all,

I've tried to run icofoam on a simple case (similar as a "flow through a circular pipe") but my courant number mean and max are explosing as soon as the first iteration ...
How could I get through this ?

CheckMesh told me that my mesh was ok.
But checkMesh also told something about sever non orthogonality :

Mesh non-orthogonality Max: 87.7597 average: 32.39
*Number of severely non-orthogonal (> 70 degrees) faces: 13240.
Non-orthogonality check OK.
<<Writing 13240 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 1.43914 OK.

Could it be a problem that explain divergence ?

If you have any idea of where it could comes from...

Thanks
Adrien

[adrieno]

Here is my geometry with my mesh (let's have a sight on photo at the bottom) .
It correspond to a pipe in two parts, first a square which is joining a cercle, then a standard pipe.

(I've juste design a quarter of it, in order to define symmetryplan later when this simple case will work)


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.000001;

vertices
(
(0 0 0)
(4500 0 0)
(4500 4500 0)
(0 4500 0)

(0 0 3970)
(3145 0 3970)
(2223.851 2223.851 3970)
(0 3145 3970)

(0 0 14760)
(3145 0 14760)
(2223.851 2223.851 14760)
(0 3145 14760)

);

blocks
(
hex (0 1 2 3 4 5 6 7) (20 20 70) simpleGrading (0.1 0.1 1)
hex (4 5 6 7 8 9 10 11) (20 20 140) simpleGrading (0.1 0.1 1)

);

edges
(
arc 5 6 (2723.650 1572.5 3970)
arc 6 7 (1572.5 2723.650 3970)

arc 9 10 (2723.650 1572.5 14760)
arc 10 11 (1572.5 2723.650 14760)

);

boundary
(
WALL
{
type wall;
faces
(
(3 7 6 2)
(2 6 5 1)

(6 10 9 5)
(7 11 10 6)
);
}

paroigauche
{
type wall;
faces
(
(0 4 7 3)
(4 8 11 7)
);
}

paroidroite
{
type wall;
faces
(
(1 5 4 0)
(4 5 9 8)
);
}

inlet
{
type patch;
faces
(
(0 3 2 1)
);
}
outlet
{
type patch;
faces
(
(8 9 10 11)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //

[rdbisme]

Well, did you do a checkMesh? Any problems with that?

What about the boundary conditions you're setting? Try to post them here!

PS: Put your files under the code bbCode.

[adrieno]

I've found a documentation on how to prepare your problem.
That will probably help me to understand and solve my problem by my own.
Here is the like to that document, it could help someone else than me.
http://www.dicat.unige.it/guerrero/o...sandtricks.pdf
I will come back later with my files if I can't make it.
Thank you,
Adrien