Hello,
I try to do CFD with axisymmetric meshes. My test case is very simple, The domain is pentahedric (inlet, upperwall, outlet, the axis and the 2 wedges.) Blockmesh works and generate the mesh (See the script at the end)
However, when I run icoFoam, I get the following errormessage, whereas I don't get it if the same mesh is 2D (inlet instead of wedges, and no axis)

Starting time loop

Time = 2e-08

Courant Number mean: 0 max: 0
#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xffffe420]
#3 Foam::polyMesh::calcDirections() const
#4 Foam::polyMesh::directions() const
#5 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&)
#6 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&)
#7 main
#8 __libc_start_main
#9 __gxx_personality_v0 at ../sysdeps/i386/elf/start.S:122
Exception en point flottant

I try to change the axis position, but it changes nothing. What can I do?

Thank you in advance
Bruno



convertToMeters 0.001;

vertices
(
(0 10 -0.436)
(0 0 -1)
(50 0 -1)
(50 10 -0.436)
(0 10 0.436)
(50 10 0.436)
);

blocks
(
hex (0 1 2 3 4 1 2 5) (10 50 1) simpleGrading (1 1 1)
);

edges
(
);

patches
(
empty axis
(
(1 2 2 1)
)

wedge back
(
(0 1 2 3)
)

wedge front
(
(4 1 2 5)
)

patch inlet
(
(0 1 1 4)
)

wall top_flange
(
(0 3 5 4)
)

patch outlet
(
(3 2 2 5)
)
);
mergePatchPairs
(
);
// ************************************************** *********************** //

[mattijs]

See http://www.cfd-online.com/cgi-bin/Op...4944#POST14944

Hi Mattijs
Thanks a lot for your advice, the test seems to work now and provide good results.

However, is the mesh still considered as axisymmetric from the solver, without an axis definition in the boundary list?

How does the solver know, where is the axis and where is the non varying direction?

Thank you in advance
Bruno