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/)
-   -   axisymetric problem with a wedge mesh (http://www.cfd-online.com/Forums/openfoam-solving/98216-axisymetric-problem-wedge-mesh.html)

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 :

Quote:

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 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)

Best,

Mat

Mat_fr March 8, 2012 12:35

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


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