CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   axisymetric problem with a wedge mesh (

Mat_fr March 6, 2012 07:02

axisymetric problem with a wedge mesh
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/"
#1 Foam::sigFpe::sigHandler(int) in "/***/platforms/linux64GccDPOpt/lib/"
at sigaction.c:0
#3 Foam::wedgePolyPatch::initTransforms() in "/***/platforms/linux64GccDPOpt/lib/"
#4 Foam::wedgePolyPatch::wedgePolyPatch(Foam::word const&, Foam::dictionary const&, int, Foam::PolyBoundaryMesh const&) in "/***/platforms/linux64GccDPOpt/lib/"
#5 Foam::PolyPatch::adddictionaryConstructorToTable<F oam::wedgePolyPatch>::New(Foam::word const&, Foam::dictionary const&, int, Foam::PolyBoundaryMesh const&) in "/***/platforms/linux64GccDPOpt/lib/"
#6 Foam::PolyPatch::New(Foam::word const&, Foam::dictionary const&, int, Foam::PolyBoundaryMesh const&) in "/***/platforms/linux64GccDPOpt/lib/"
#7 Foam::PolyBoundaryMesh::PolyBoundaryMesh(Foam::IOo bject const&, Foam::PolyMesh const&) in "/***/platforms/linux64GccDPOpt/lib/"
#8 Foam::PolyMesh::PolyMesh(Foam::IOobject const&) in "/***/platforms/linux64GccDPOpt/lib/"
#9 Foam::fvMesh::fvMesh(Foam::IOobject const&) in "/***/platforms/linux64GccDPOpt/lib/"
in "/software/oss/OpenFOAM/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/simpleFoam"
#11 __libc_start_main in "/lib64/"
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 !


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 March 7, 2012 05:18

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)



Mat_fr March 8, 2012 12:35

I finally found the solution of my problem, here:

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


All times are GMT -4. The time now is 01:48.