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/)
-   -   [Other] topoSet Wrong number of faces. (https://www.cfd-online.com/Forums/openfoam-meshing/152992-toposet-wrong-number-faces.html)

a19910112a May 14, 2015 09:56

topoSet Wrong number of faces.
 
2 Attachment(s)
Hi Foamers. I came across the same problem as describeded in http://www.cfd-online.com/Forums/openfoam-meshing-blockmesh/128416-problems-holes-gaps-baffle-internal-wall-zero-thickness.html.

I did the same thing as described in the referred post.

To solve this problem, I located the problem at topoSet. It recognise a wrong number of faces.

My blockMeshDict is as :
Code:

convertToMeters 1;

vertices
(

        (0    0    0    )    //    0
        (2.4424    0    0.1065158)    //    1
        (2.4424    0    -0.1065158)    //    2
        (0    1    0    )    //    3
        (2.4424    1    0.1065158)    //    4
        (2.4424    1    -0.1065158)    //    5
        (4.25    0    0.1853472)    //    6
        (4.25    0    -0.1853472)    //    7
        (4.25    1    0.1853472)    //    8
        (4.25    1    -0.1853472)    //    9 

 
);

//local axes (z x y)

blocks
(

        hex    (1      2    0    0    4    5    3    3    )    innerCore        (1    26    10)    simpleGrading    (1    1    1)
        hex    (6    7    2    1    8    9    5    4    )    radialReflector        (1    18    10)    simpleGrading    (1    1    1)
);

edges
(
);
...

This generate a mesh of a wedge which is shown in the first attached pic. There are 10 meshes in the y direction.

I want to add a baffle wall in the middle of the wedge, which is shown in the second attachment.

The toposet file :
Code:

actions
(
 
    {
        name    baffle;
        type    faceSet;
        action  new;
        source  boxToFace;
        sourceInfo     
    box (2.4424    0 -0.107) (2.4424 1 0.107); //24
        }
    }

    {
        name    baffle;
        type    faceZoneSet;
        action  new;
        sourceInfo
        source  setToFaceZone;
        {
            faceSet baffle;
        }
    }

);

the log it generates:
Code:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh fluid for time = 0

Reading topoSetDict

Time = 0
    mesh not changed.
Created faceSet baffle
    Applying source boxToFace
    Adding faces with centre within boxes 1((2.4424 0 -0.107) (2.4424 1 0.107))
    faceSet baffle now size 9
Created faceZoneSet baffle
    Applying source setToFaceZone
    Adding all faces from faceSet baffle ...
    faceZoneSet baffle now size 9

End

As can be seen, topoSet recognises only 9 faces where 10 are expected, which gives the baffle shown in the second pic.

Can anyone tell me where the problem could be and how can I solve it ?

By the way, I also tried to use points 42, 15, 51, as the two points after the box item in topoSet. In those cases only 0 face could be recognised. This seems to be due to some similar or related reasons.

Thanks

Howard June 9, 2015 05:41

Quote:

Originally Posted by a19910112a (Post 546315)
Hi Foamers. I came across the same problem as describeded in http://www.cfd-online.com/Forums/openfoam-meshing-blockmesh/128416-problems-holes-gaps-baffle-internal-wall-zero-thickness.html.

I did the same thing as described in the referred post.

To solve this problem, I located the problem at topoSet. It recognise a wrong number of faces.

My blockMeshDict is as :
Code:

convertToMeters 1;
 
vertices
(
 
        (0    0    0    )    //    0
        (2.4424    0    0.1065158)    //    1
        (2.4424    0    -0.1065158)    //    2
        (0    1    0    )    //    3
        (2.4424    1    0.1065158)    //    4
        (2.4424    1    -0.1065158)    //    5
        (4.25    0    0.1853472)    //    6
        (4.25    0    -0.1853472)    //    7
        (4.25    1    0.1853472)    //    8
        (4.25    1    -0.1853472)    //    9 
 
 
);
 
//local axes (z x y)
 
blocks
(
 
        hex    (1      2    0    0    4    5    3    3    )    innerCore        (1    26    10)    simpleGrading    (1    1    1)
        hex    (6    7    2    1    8    9    5    4    )    radialReflector        (1    18    10)    simpleGrading    (1    1    1)
);
 
edges
(
);
...

This generate a mesh of a wedge which is shown in the first attached pic. There are 10 meshes in the y direction.

I want to add a baffle wall in the middle of the wedge, which is shown in the second attachment.

The toposet file :
Code:

actions
(
 
    {
        name    baffle;
        type    faceSet;
        action  new;
        source  boxToFace;
        sourceInfo     
    box (2.4424    0 -0.107) (2.4424 1 0.107); //24
        }
    }
 
    {
        name    baffle;
        type    faceZoneSet;
        action  new;
        sourceInfo
        source  setToFaceZone;
        {
            faceSet baffle;
        }
    }
 
);

the log it generates:
Code:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
 
Create polyMesh fluid for time = 0
 
Reading topoSetDict
 
Time = 0
    mesh not changed.
Created faceSet baffle
    Applying source boxToFace
    Adding faces with centre within boxes 1((2.4424 0 -0.107) (2.4424 1 0.107))
    faceSet baffle now size 9
Created faceZoneSet baffle
    Applying source setToFaceZone
    Adding all faces from faceSet baffle ...
    faceZoneSet baffle now size 9
 
End

As can be seen, topoSet recognises only 9 faces where 10 are expected, which gives the baffle shown in the second pic.

Can anyone tell me where the problem could be and how can I solve it ?

By the way, I also tried to use points 42, 15, 51, as the two points after the box item in topoSet. In those cases only 0 face could be recognised. This seems to be due to some similar or related reasons.

Thanks

Hi friend, I also have some problem that when I use the toposet it shows cellzoneset size 0. but I think my boxes is defined right. Have you solve your problem?


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