CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Axisymmetric cases Test with blockmesh and icoFoam (http://www.cfd-online.com/Forums/openfoam-solving/59574-axisymmetric-cases-test-blockmesh-icofoam.html)

frackowi July 11, 2007 16:39

Hello, I try to do CFD with
 
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#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 July 11, 2007 17:47

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

frackowi July 11, 2007 19:02

Hi Mattijs Thanks a lot for y
 
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


All times are GMT -4. The time now is 12:05.