CFD Online Discussion Forums

CFD Online Discussion Forums (
-   Open Source Meshers: Gmsh, Netgen, CGNS, ... (
-   -   Axisymmetric cylindrical meshes (

jgaricano March 3, 2008 11:33

Greetings to you all! I'm n
Greetings to you all!

I'm new with OpenFOAM. I need to model an axisymmetric cylindrical problem (very similarj to a Burke-Schumann configuration) but in2D.

Could any of you send me a Block Mesh dictionary in which appears the definition of an axisymmetric cylindrical geometry but treated as an 2D problem?
Thanks in advance.

mike_jaworski March 3, 2008 11:55

Hi Jesus, I'm not sure I un
Hi Jesus,
I'm not sure I understand your problem. Do you want to build a 2D axisymmetric mesh and need an example, or are you trying to approximate axisymmetry with a cartesian coordinate mesh?

If you want to use an axisymmetric mesh for your case, the icoDyMFoam tutorial "moving cone" has such a mesh generated with the blockMesh utility.

Good Luck,
Mike J.

jgaricano March 4, 2008 09:12

Hi Michael: I've started mo
Hi Michael:

I've started modifying the geometry related to the lesInterFoam tutorial in order to get a longitudinal section of an arrangement of two coaxial tubes (the inner one being shorter than the outer one). This intends to be a 2D model of a combustion chamber.

I've developed the blockMeshDictionary but when running blockMesh I get the following:

Reading block mesh description dictionary

Creating block mesh

Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology

Default patch type set to empty

--> FOAM FATAL ERROR : Trying to specify a boundary face 4(2 8 30 25) on the face on cell 3 which is either an internal face or already belongs to some other patch. This is face 0 of patch 4 named separatingWall.#0 Foam::error::printStack(Foam: in "/home/jgaricano/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#1 Foam::error::abort() in "/home/jgaricano/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#2 Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Field<foam::vector<double> > const&, Foam::List<foam::cellshape> const&, Foam::List<foam::list<foam::face> > const&, Foam::List<foam::word> const&, Foam::List<foam::word> const&, Foam::word const&, Foam::List<foam::word> const&, bool) in "/home/jgaricano/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/"
#3 Foam::blockMesh::createTopology(Foam::IOdictionary &) in "/home/jgaricano/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/blockM esh"
#4 Foam::blockMesh::blockMesh(Foam::IOdictionary&) in "/home/jgaricano/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/blockM esh"
#5 main in "/home/jgaricano/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/blockM esh"
#6 __libc_start_main in "/lib/"
#7 Foam::regIOobject::readIfModified() in "/home/jgaricano/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/blockM esh"

The problem is in the nonzero width wall of the inner tube.

Thanks in advanced for your interest.

mike_jaworski March 4, 2008 14:12

I don't think I've seen that p
I don't think I've seen that particular error before. If you've never used blockMesh before, I think the biggest hurdle is making sure your vertices are defined in the correct order for both the patches and the blocks. There's information in the user manual which has all the "rules" for doing so. You may need to check that the patches all have their face-normals pointing away from the computational domain.

Good luck,
Mike J.

mattijs March 4, 2008 15:07

Sometimes it helps to visualis
Sometimes it helps to visualise the points. blockMesh with -blockTopology option will dump an .obj file which can be viewed in e.g. javaview or paraview (using objToVTK).

jgaricano March 4, 2008 15:11

Thanks a lot for your interest
Thanks a lot for your interest. Fortunately I've been able to solve the problem. It was a mistaken space. I'm very sorry of having troubled you.

See you.

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