CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Mesh conversion problem from Salome to openfoam (http://www.cfd-online.com/Forums/openfoam/73138-mesh-conversion-problem-salome-openfoam.html)

jishnusoni February 28, 2010 11:59

Mesh conversion problem from Salome to openfoam
 
Hello,

I am trying to create an impinging jet. I have created an geometry and mesh in Salome (Automatic tetrahedralisation). When I try to import mesh from salome to openfoam using ideasUnvToFoam, I get a constant error which is below,

Can anyone please help me.

thanks in advance

regards
jish






Constructing mesh with non-default patches of size:
walls 13708
inlet 75
pipe-outlet 75
outlet 1668



Trying to specify a boundary face 3(353 5484 8) on the face on cell 37553 which is either an internal face or already belongs to some other patch. This is face 4743 of patch 0 named walls.#0 Foam::error::printStack(Foam::Ostream&) in "/home/caelinux/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/caelinux/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#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::word const&, Foam::List<Foam::word> const&, bool) in "/home/caelinux/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#3 main in "/home/caelinux/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/ideasUnvToFoam"
#4 __libc_start_main in "/lib/libc.so.6"
#5 __gxx_personality_v0 in "/home/caelinux/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/ideasUnvToFoam"


From function polyMesh::polyMesh
(
const IOobject& io,
const pointField& points,
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchTypes,
const wordList& boundaryPatchNames,
const word& defaultBoundaryPatchType
)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 483.

FOAM aborting

Aborted

linnemann February 28, 2010 16:06

Hi

The ideasUnvToFoam converter have been updated with better support in the 1.6.x release.

Although it looks like you have an internal face in you mesh from the output.

You must make sure you have no internal faces in Salome before meshing as OpenFOAM do not handle internal faces like ex. Fluent does.

In Salome you can go to the "repair" menu where you can heal/stitch/remove faces etc. See if it wont be possible to clean up the geometry even further.

Else export the Salome script and attach it here and I will give it a whirl.

jishnusoni February 28, 2010 18:24

1 Attachment(s)
Hello linnemann,

Thanks very much for you reply.

I am trying to use the repair tool, but not too sure what needs to be done. I have attached my Salome file, if you could look at it and give some feedback, that would be a great help for me.

thanks

jish

linnemann March 1, 2010 03:20

1 Attachment(s)
Hi I had a look and it is because you have an internal face where the pipe and the box is fused.

I've created a small wall with a diameter slightly larger than the inner diameter of the pipe.

See if the attached wont work for you,

jishnusoni March 1, 2010 06:48

I tried your mesh, I think its more better then mine and getting some results with it. But I wanted to specify the cylinder outlet face. When I tried to do that I am getting again the same error. But if I dont specify the cylinder outlet face, then its working fine.

Could you please suggest me how can I specify the cylinder outlet face?

regards

jish

linnemann March 1, 2010 06:59

Hi

Unfortunately you cant. It would mean an internal face which is not possible with OpenFOAM.

You have one possibility if you really want a BC at the jet outlet.

Use Ggi from OpenFOAM-1.5-dev. You create two separate meshes, one for the box and one for the cylinder. Put the mesh in separate cases and use mergeMeshes to union them into one case.

Then setup you BC's connecting the two meshes like in this article

jishnusoni March 1, 2010 19:09

Hello,

Thanks for your response.

I am just trying to work on your geometry. I am trying to make a similar geometry as you send me. Can you please explain me what did you do 'Cylinder_;vertex_12, cylinder_1;edge_7 and with vertex_2'

I am little confused about the above things.

regards
jish

linnemann March 2, 2010 02:46

1 Attachment(s)
Hi

Its done creating a new vertex from the middle of a face.

I just did it to make a point in the middle of the jet outlet face so I could create a slightly larger cylinder to cut the box with.

See attached screen-shot

http://www.cfd-online.com/Forums/att...1&d=1267512353

jishnusoni March 2, 2010 04:28

hi
Thanks for your response.
can you also explain what is disk_1, as I thought its a circle with a face given, but when I checked it, it wasn't the same.

thanks

regards
jish

linnemann March 2, 2010 04:44

Hi

disk_1 is just a disk created from the new vertex and a vector. in this instance I used the line stradling cylinder_1 as my direction vector.

Then I just extruded the disk in the z-direction. which became extrusion_1.

I could have just created a cylinder with the new vertex as base point, but this is just how I did it :).

There are numerous ways in which you could arrive at the same output.

All roads lead to Rome, or in this instance to an impinging jet.

jishnusoni March 2, 2010 05:20

1 Attachment(s)
hi,

The thing is that I have only started recently to use OF and Salome, so I don't know all the commands. I think the 'Disk' command is not available in Salome 4.1 (as I tried searching for it but couldnt find it). Moreover, I have tried in different ways to do the impinging jet but I can not Mesh it. could you please check my geom and give me feedback. thanks very much...

jish

linnemann March 2, 2010 05:42

Hi

I will look into the case, but I would suggest that you upgrade your Salome install.

Are you using CAELinux?

I would suggest you try and install CentOS and get the Salome version for Scientific Linux 64bit from Salomes website. www.salome-platform.org. (SL and CentOS are both based on the same RHEL sources, so they are compatible)

OpenFOAM is easy to install on CentOS since I've created a scripted install for it.

linnemann March 2, 2010 05:49

1 Attachment(s)
Hi

I have no problems creating a mesh using simple Netgen 3D-2D-1D parameters.

I although think you should make a hex mesh on this simple geometry. Ill post an example when I have it.

linnemann March 2, 2010 08:01

2 Attachment(s)
Hello

Here is the same geometry made so it can be hex meshed.

Hope you can figure the case structure out. I don't have anymore time to use on this as you just need experience with Salome from this point on.

The tar file is a script dump which you need to run in Salome using "file" -> "Load script". Open the file without _GEOM and _SMESH in the name.

http://www.cfd-online.com/Forums/att...1&d=1267531242

jishnusoni March 2, 2010 19:13

Hello,

Thanks very much for you time. I really appreciate it.

Now, I will play around Salome and learn it.

Just a last thing, the Hex files which you send me are not opening. Is the file corrupt as its not an .hdf files or its something else.

Yes, I am using CAElinux. Whats the difference between CentOS and Caelinux. Cant I just upgrade Salome directly here in CAElinux?

thanks

regards
jish

linnemann March 3, 2010 02:53

Hi

I answered how to open the file in my previous post. If you cant figure that out give me a private message with your email and I will send the HDF to you. It is just 50kb larger than the forum allows.

The difference between CAELinux and CentOS is that CAELinux is a complete CAE package with all the applications already installed. CentOS on the other hand is like windows. It comes only with the basic stuff and you have to install Salome and OpenFOAM yourself.

Else use the internet search function (Google, yahoo, etc) and you will find answers to your questions. I'm saying this politely, but some stuff you really need to figure out yourself otherwise people wont bother helping.


All times are GMT -4. The time now is 06:46.