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

[mesh manipulation] Problem with topoSet and FaceZone creation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 12, 2013, 08:14
Default Problem with topoSet and FaceZone creation
  #1
Member
 
Andrew Glassby
Join Date: Sep 2009
Posts: 65
Rep Power: 17
ADGlassby is on a distinguished road
Hi,

I'm hoping someone might have come across this before and found the solution!

I'm trying to use topoSet to create two faceZoneSets inside and outside of a circle, I'm working in 2D and just using the basic cavity tutorial to work out the right method. The goal is, once topoSet haas worked, to create Baffles in order to create an AMI interface. Note that the createBaffles and mergeOrSplitBaffles commands seem to do what I expect.

basically I have created a square mesh by editting the blockMeshDict file to give me a mesh big enough to play with. I have created a .obj cylindrical geometry file for use with surfaceToCell directive in my topoSetDict.

I have built up the topoSetDict file slowly checking at each step that things have gone as I had expected them to go. I have seen the cellSet created inside and outside of the cylinder using foamToVTK -cellSet <setname> and then viewing in paraFoam. I also used checkMesh to double check that nothing was going wrong there.

Things start going wrong when I try to create the final faceZoneSet. I have created an error log for this. Note that if I try to use setsToZones -noFlipMap instead of using the faceZoneSet directives in topoSetDict things SEEM to go OK until I createBaffle then use checkMesh and I find that instead of getting 2 zones (inside and outside of the cylinder) I get multiple regions reported. I have attached a log of this checkMesh result too.

I have attached my topoSetDict, blockMesh, createBaffleDict, .obj file and an error log from the topoSet run and checkMesh. All zipped up.

I hope someone can help me!

Best Regards
Andrew
Attached Files
File Type: zip cfd-online.zip (36.1 KB, 96 views)
ADGlassby is offline   Reply With Quote

Old   August 12, 2013, 11:22
Default
  #2
Senior Member
 
mohsen kh
Join Date: Jan 2013
Location: Iran
Posts: 125
Rep Power: 14
m5m5kh is on a distinguished road
hi dear foamers
I need your help
I used Gambit for modeling of my case which was is a furnace
I want to model flow in a furnace considering it as a porous medium with constant porosity of 0.5.
I do not know how to modify the long tube-shape region as a porous media
please help me if you can.
some one told me I should modify topoSetDict but I do not know how to do it
please help me,I am not so advanced user of openfoam.
If you know please tell me in detail.

thanks
Mohsen
Attached Images
File Type: jpg Gambit_ImageFile1.jpg (11.6 KB, 37 views)
File Type: jpg Gambit_ImageFile2.jpg (36.4 KB, 40 views)
File Type: jpg Gambit_ImageFile4.jpg (55.8 KB, 27 views)
m5m5kh is offline   Reply With Quote

Old   September 13, 2013, 06:59
Default
  #3
New Member
 
Join Date: Aug 2013
Posts: 20
Rep Power: 13
KYPCK444 is on a distinguished road
I don't know how to help ..

.. but I have another question xD

about topoSet, my problem is:

I have a pipe, with inlet, outlet and wall patches;
I'm using topoSet in order to be able to set (in a second time) inlet and outlet as cyclic (I'm using snappyHexMesh, which does not allow itself to set cyclics boundaries);
the point is that after running topoSet, I obtain inlet and outlet with different number of face centres (quite the same, but not properly the same):

topoSetDict:

actions
(
{
name inlet_set;
type faceSet;
action new;
source patchToFace;
sourceInfo
{
name inlet;
}
}

{
name outlet_set;
type faceSet;
action new;
source patchToFace;
sourceInfo
{
name outlet;
}
}
);


log:

Create time

Create polyMesh for time = 0

Reading topoSetDict

Time = 0
mesh not changed.
Created set inlet_set
Applying source patchToFace
Adding all faces of patch inlet ...
Found matching patch inlet with 3597 faces.
Set inlet_set now size 3597
Created set outlet_set
Applying source patchToFace
Adding all faces of patch outlet ...
Found matching patch outlet with 3690 faces.
Set outlet_set now size 3690

End


the very problem arises when I run createPatch (which sets cyclis to inlet and outlet manipulated by topoSet), in that it gives "serious error" because of this face number "mismatch".

any hint?

thanks in advance
KYPCK444 is offline   Reply With Quote

Old   September 13, 2013, 07:31
Default
  #4
Member
 
Andrew Glassby
Join Date: Sep 2009
Posts: 65
Rep Power: 17
ADGlassby is on a distinguished road
Why don't you use createPatch instead? You can take the patch as created by sHM and redefine it as a cyclic or any other patch type. Have a look in the src directories for createPatchDict, this should help I believe.

Andrew
ADGlassby is offline   Reply With Quote

Old   September 13, 2013, 07:56
Default
  #5
New Member
 
Join Date: Aug 2013
Posts: 20
Rep Power: 13
KYPCK444 is on a distinguished road
Quote:
Originally Posted by ADGlassby View Post
Why don't you use createPatch instead? You can take the patch as created by sHM and redefine it as a cyclic or any other patch type. Have a look in the src directories for createPatchDict, this should help I believe.

Andrew

thanks for your reply.

truly, it sounds better to me, in fact I already tried this way, but unfortunately it did not work. Moreover, I found somewhere kind of instructions about how to apply cyclics with snappyHexMesh, which involved topoSet.
KYPCK444 is offline   Reply With Quote

Old   September 13, 2013, 08:18
Default
  #6
New Member
 
Join Date: Aug 2013
Posts: 20
Rep Power: 13
KYPCK444 is on a distinguished road
Quote:
Originally Posted by KYPCK444 View Post
thanks for your reply.

truly, it sounds better to me, in fact I already tried this way, but unfortunately it did not work. Moreover, I found somewhere kind of instructions about how to apply cyclics with snappyHexMesh, which involved topoSet.

D:

it seems to work now, without topoSet, just setting inlet/outlet cyclics with createPatch with constructFrom patches and not from set.

(Don't know what I could have missed the first time a tried, days ago)

thanks a lot
KYPCK444 is offline   Reply With Quote

Old   September 13, 2013, 08:21
Default
  #7
Member
 
Andrew Glassby
Join Date: Sep 2009
Posts: 65
Rep Power: 17
ADGlassby is on a distinguished road
You are welcome.... not sure I did much though ;-)

Andrew
ADGlassby is offline   Reply With Quote

Old   September 16, 2013, 12:05
Default
  #8
New Member
 
Join Date: Aug 2013
Posts: 20
Rep Power: 13
KYPCK444 is on a distinguished road
I was wrong ..

it is true that topoSet is not necessary in my case, but since number of faces are different between neighbour patches, createPatch crashes (I have opened a dedicated topic: http://www.cfd-online.com/Forums/ope...eatepatch.html).
KYPCK444 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
fvOptions, faceZone, directionalPressureGradientExplicitSource Alasir OpenFOAM Running, Solving & CFD 0 November 24, 2017 06:04
How to delete faces in faceZone and update mesh during runtime? nuttita OpenFOAM Programming & Development 0 September 15, 2017 19:08
[snappyHexMesh] snappyHexMesh problem shengqiming OpenFOAM Meshing & Mesh Conversion 0 December 15, 2016 09:35
Calculate Mass Flowrate at a faceZone with OpenFOAM4 donQi OpenFOAM Post-Processing 1 August 31, 2016 22:13
simulation of radial blower - topoSetDict problem nash OpenFOAM Pre-Processing 1 August 14, 2013 09:26


All times are GMT -4. The time now is 10:15.