[Mat_fr]

Dear all,

I wish to simulate an axisymmetric incompressible flow in a pipe (using simpleFoam and the standard k-eps model).
So, I created a wedge mesh in gambit and defined the different BC, as walls, inlet, outlet, and the two faces "front" and "back" of the wedge.
In the polyMesh/boundary file, I specified the wedge type for front and back. And I did the same for the physical variables at time 0.

However, when I'm running the simulation, I get the following error :

Create time

Create mesh for time = 0

#0 Foam::error:: printStack(Foam::Ostream&) in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2
at sigaction.c:0
#3 Foam::wedgePolyPatch::initTransforms() in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::wedgePolyPatch::wedgePolyPatch(Foam::word const&, Foam::dictionary const&, int, Foam::PolyBoundaryMesh const&) in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::PolyPatch::adddictionaryConstructorToTable<F oam::wedgePolyPatch>::New(Foam::word const&, Foam::dictionary const&, int, Foam::PolyBoundaryMesh const&) in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::PolyPatch::New(Foam::word const&, Foam::dictionary const&, int, Foam::PolyBoundaryMesh const&) in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::PolyBoundaryMesh::PolyBoundaryMesh(Foam::IOo bject const&, Foam::PolyMesh const&) in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#8 Foam::PolyMesh::PolyMesh(Foam::IOobject const&) in "/***/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#9 Foam::fvMesh::fvMesh(Foam::IOobject const&) in "/***/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10
in "/software/oss/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/simpleFoam"
#11 __libc_start_main in "/lib64/libc.so.6"
#12
in "/software/oss/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)

I check the tutorial OpenFOAM-2.0.1/tutorials/incompressible/pimpleDyMFoam/movingCone, and one can see that the axis of symmetry is also specified in the BC. Is it necessary ? I don't think so. Anyway I've defined it in Gambit, but when I'm using fluent3DMeshToFoam the axis is not created in the boundary file.
Does somebody has an idea ?

Thank you for your help !

Mat

ps :
* I'm not sure if my thread is in the good forum, but as I don't know the reason of my error...

[Mat_fr]

For information, when I'm using symmetryPlane boundary conditions, instead of wedge, with the same mesh, the computation starts. But in this case the calculation diverges (which seems logical)

Best,

Mat

[Mat_fr]

I finally found the solution of my problem, here:
http://www.ltengsoft.com/node/10

It was something quite tricky. One must take care of the mesh position in the space.


Mat