CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Mesh conversion problem from Salome to openfoam

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 28, 2010, 11:59
Default Mesh conversion problem from Salome to openfoam
  #1
New Member
 
Join Date: Feb 2010
Posts: 21
Rep Power: 16
jishnusoni is on a distinguished road
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:rintStack(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:olyMesh:olyMesh(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:olyMesh
(
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
jishnusoni is offline   Reply With Quote

Old   February 28, 2010, 16:06
Default
  #2
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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.
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   February 28, 2010, 18:24
Default
  #3
New Member
 
Join Date: Feb 2010
Posts: 21
Rep Power: 16
jishnusoni is on a distinguished road
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
Attached Files
File Type: gz impinging jet1.hdf.tar.gz (33.0 KB, 22 views)
jishnusoni is offline   Reply With Quote

Old   March 1, 2010, 03:20
Default
  #4
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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,
Attached Files
File Type: gz impinging jet1.hdf.tar.gz (39.5 KB, 32 views)
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   March 1, 2010, 06:48
Default
  #5
New Member
 
Join Date: Feb 2010
Posts: 21
Rep Power: 16
jishnusoni is on a distinguished road
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
jishnusoni is offline   Reply With Quote

Old   March 1, 2010, 06:59
Default
  #6
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   March 1, 2010, 19:09
Default
  #7
New Member
 
Join Date: Feb 2010
Posts: 21
Rep Power: 16
jishnusoni is on a distinguished road
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
jishnusoni is offline   Reply With Quote

Old   March 2, 2010, 02:46
Default
  #8
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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

Attached Images
File Type: jpg SalomeVertexMiddleFace.jpg (50.7 KB, 521 views)
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   March 2, 2010, 04:28
Default
  #9
New Member
 
Join Date: Feb 2010
Posts: 21
Rep Power: 16
jishnusoni is on a distinguished road
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
jishnusoni is offline   Reply With Quote

Old   March 2, 2010, 04:44
Default
  #10
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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.
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   March 2, 2010, 05:20
Default
  #11
New Member
 
Join Date: Feb 2010
Posts: 21
Rep Power: 16
jishnusoni is on a distinguished road
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
Attached Files
File Type: gz impingingjetstudy.tar.gz (48.4 KB, 3 views)
jishnusoni is offline   Reply With Quote

Old   March 2, 2010, 05:42
Default
  #12
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   March 2, 2010, 05:49
Default
  #13
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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.
Attached Files
File Type: gz impingingStudy1.hdf.tar.gz (29.1 KB, 7 views)
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   March 2, 2010, 08:01
Default
  #14
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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.

Attached Images
File Type: jpg impingingJetHex.jpg (96.1 KB, 518 views)
Attached Files
File Type: gz impingingHexMesh.tar.gz (6.2 KB, 18 views)
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   March 2, 2010, 19:13
Default
  #15
New Member
 
Join Date: Feb 2010
Posts: 21
Rep Power: 16
jishnusoni is on a distinguished road
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
jishnusoni is offline   Reply With Quote

Old   March 3, 2010, 02:53
Default
  #16
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
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.
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Dynamic Mesh Problem. Tom Clark FLUENT 10 June 21, 2021 05:27
Moving mesh problem OpenFoam 141 kassiotis OpenFOAM Running, Solving & CFD 30 April 15, 2015 00:10
[Other] Missing BC in mesh conversion from Star-cd to OpenFOAM piccinini OpenFOAM Meshing & Mesh Conversion 2 June 19, 2009 09:12
[Commercial meshers] Problem importing mesh in openfoam from fluent alessandr0 OpenFOAM Meshing & Mesh Conversion 3 September 4, 2008 14:41
unstructured vs. structured grids Frank Muldoon Main CFD Forum 1 January 5, 1999 11:09


All times are GMT -4. The time now is 17:22.