CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [mesh manipulation] Problem with topoSet and FaceZone creation (https://www.cfd-online.com/Forums/openfoam-meshing/122095-problem-toposet-facezone-creation.html)

ADGlassby August 12, 2013 07:14

Problem with topoSet and FaceZone creation
 
1 Attachment(s)
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

m5m5kh August 12, 2013 10:22

3 Attachment(s)
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

KYPCK444 September 13, 2013 05:59

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

ADGlassby September 13, 2013 06:31

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

KYPCK444 September 13, 2013 06:56

Quote:

Originally Posted by ADGlassby (Post 451497)
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 September 13, 2013 07:18

Quote:

Originally Posted by KYPCK444 (Post 451504)
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 :)

ADGlassby September 13, 2013 07:21

You are welcome.... not sure I did much though ;-)

Andrew

KYPCK444 September 16, 2013 11:05

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


All times are GMT -4. The time now is 19:57.