CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Native Meshers: blockMesh (http://www.cfd-online.com/Forums/openfoam-meshing-blockmesh/)
-   -   FOAM FATAL ERROR: face 0 in patch 0 does not have neighbour cell face (http://www.cfd-online.com/Forums/openfoam-meshing-blockmesh/112905-foam-fatal-error-face-0-patch-0-does-not-have-neighbour-cell-face.html)

am9109 February 7, 2013 08:27

FOAM FATAL ERROR: face 0 in patch 0 does not have neighbour cell face
 
Hello Everyone,

I am very new to OpeFoam as I only started to use it a few days ago, I am trying to generate a Mesh which is made of 3 rectangular blocks, one next to the other with the middle one slightly smaller than the other two.
I get the following message:

FOAM FATAL ERROR:
face 0 in patch 0 does not have neighbour cell face: 4(0 12 16 3)

What am I doing wrong?
Thanks a lot

P.S.: here is the code:
convertToMeters 1;

vertices
(
(-5 -1 -1)
(5 -1 -1)
(5 1 -1)
(-5 1 -1)
(-5 -1 1)
(5 -1 1)
(5 1 1)
(-5 1 1)
(-2 -1 -0.5)
(-2 1 -0.5)
(2 1 -0.5)
(2 -1 -0.5)
(-2 -1 -1 )
(-2 1 -1)
(2 1 -1)
(2 -1 -1)
(-2 -1 1)
(-2 1 1)
(2 1 1)
(2 -1 1)
);

blocks
(
hex (0 4 3 7 12 13 16 17) (100 4 4) simpleGrading (1 1 1)
hex (8 9 10 11 16 17 18 19) (100 4 4) simpleGrading (1 1 1)
hex (1 2 5 6 14 15 18 19) (100 4 4) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
leftWall
{
type wall;
faces
(
(0 12 16 3)
(8 11 19 16)
(1 2 19 15)
);
}
rightWall
{
type wall;
faces
(
(4 7 17 13)
(9 17 18 10)
(14 18 6 5)
);
}
topWall
{
type wall;
faces
(
(7 3 16 17)
(17 16 19 18)
(18 19 2 6)
);
}
bottomWall
{
type wall;
faces
(
(0 4 13 12)
(13 9 8 12)
(9 10 11 8)
(15 11 10 14)
(14 5 1 15)
);
}
inlet
{
type patch;
faces
(
(0 3 7 4)
);
}
outlet
{
type patch;
faces
(
(6 2 1 5)
);
}
);
patch slave_1
(
(13 17 16 12)
);

patch slave_2
(
(15 19 18 14)
) ;

mergePatchPairs
(
(16 17 9 8)
(10 18 19 11)
);

colinB February 7, 2013 10:30

Hi and welcome to the forum,

face 0 in patch 0 does not have neighbour cell face: 4(0 12 16 3)

actually means that you messed up the order of the points
in your patch definition.
In blockMesh you have to strictly follow the orders which are
well explained in the user guide.

However sometimes it is hard to see where the error is therefore I
recommend you to type

paraFoam -block

which will display your blockMeshDict file without compiling it so
it is very handy for error handling.

Another hint for the future is try to use the forum search or google
on your error message, which gives you already a lot of answers
regarding your problems.

I hope I could contribute
regards

apk1509 August 30, 2013 05:41

FOAM FATAL ERROR: face 7 in patch 0 does not have neighbour cell face: 4(23 24 20 19
 
Hello Everyone...

I am very new to OF and trying to mesh a cylinder for Hagen-Poiseuille flow. I am getting error like:

FOAM FATAL ERROR:

face 7 in patch 0 does not have neighbour cell face: 4(19 20 24 23)

From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 127.

blockMeshDict is :

convertToMeters 0.1;

vertices
(
(0 0 0)
(0.15 0.15 0)
(-0.15 0.15 0)
(-0.15 -0.15 0)
(0.15 -0.15 0) //4
(0.3 0.3 0)
(-0.3 0.3 0)
(-0.3 -0.3 0)
(0.3 -0.3 0)//8
(0.35 0.35 0)
(-0.35 0.35 0)
(-0.35 -0.35 0)
(0.35 -0.35 0)//12

(0.15 0.15 30)
(-0.15 0.15 30)
(-0.15 -0.15 30)
(0.15 -0.15 30)
(0.3 0.3 30)//17
(-0.3 0.3 30)
(-0.3 -0.3 30)
(0.3 -0.3 30)
(0.35 0.35 30) //21
(-0.35 0.35 30)
(-0.35 -0.35 30)
(0.35 -0.35 30)
(0 0 30) //25

);

blocks
(
hex (1 2 3 4 13 14 15 16) (10 10 150) simpleGrading (1 1 1)
hex (1 5 6 2 13 17 18 14) (10 10 150) simpleGrading (1 1 1)
hex (2 6 7 3 14 18 19 15) (10 10 150) simpleGrading (1 1 1)
hex (3 7 8 4 15 19 20 16) (10 10 150) simpleGrading (1 1 1)
hex (4 8 5 1 16 20 17 13) (10 10 150) simpleGrading (1 1 1)
hex (5 9 10 6 17 21 22 18) (10 10 150) simpleGrading (1 1 1)
hex (6 10 11 7 18 22 23 19) (10 10 150) simpleGrading (1 1 1)
hex (12 9 5 8 24 21 17 20) (10 10 150) simpleGrading (1 1 1)
);

edges
(
arc 5 6 (0 0.42426 0) // back face
arc 6 7 (-0.42426 0 0)
arc 7 8 (0 -0.42426 0)
arc 8 5 (0.42426 0 0)
arc 9 10 (0 0.5 0)
arc 10 11 (-0.5 0 0)
arc 11 12 (0 -0.5 0)
arc 12 9 (0.5 0 0)

arc 17 18 (0 0.42426 30) // front face
arc 18 19 (-0.42426 0 30)
arc 19 20 (0 -0.42426 30)
arc 20 17 (0.42426 0 30)
arc 21 22 (0 0.5 30)
arc 22 23 (-0.5 0 30)
arc 23 24 (0 -0.5 30)
arc 24 21 (0.5 0 30)

);

boundary
(
inlet
{
type inlet;
faces
(
(13 16 15 14)
(13 14 18 17)
(14 15 19 18)
(15 16 20 19)
(16 13 17 20)
(17 18 22 21)
(18 19 23 22)
(19 20 24 23)
(20 17 21 24)
);
}

outlet
{
type outlet;
faces
(
(1 5 6 2)
(2 6 7 3)
(3 7 8 4)
(4 8 5 1)
(1 2 3 4)
(5 9 10 6)
(6 10 11 7)
(7 11 12 8)
(8 12 9 5)
);
}

fixedWalls
{
type wall;
faces
(
(21 22 10 19)
(22 23 11 10)
(23 24 12 11)
(24 21 9 12)
);
}
);

mergePatchPairs
(
);

colinB August 30, 2013 05:47

Dear apk1509,

please read my post of Feb. 7th which is right before your
post. (post no. 2 of this thread)
The complete answer to your question is posted there as well
as hints on how to solve your problem!

regards
Colin

apk1509 August 30, 2013 05:57

Actually I tried to correct the order of points and still i got the same error..
Also, when i tried with "paraFoam -block", paraview opened and closed before showing anything.

How to rectify this..

colinB August 30, 2013 06:41

Hi,

this brings us closer to your problem.
You see why it is important to explain your problem in detail
as suggested here: http://www.cfd-online.com/Forums/ope...-get-help.html

However two remarks to your blockMeshDict:

- the way you count is wrong: blockMesh starts counting from 0
so the first point in the list with points is 0 the second is 1 and so on!

- the valid entry for patchtypes are

wall
patch
symmetry

the patchtype like inlet outlet and so on will be specified in the 0 folder.
see the examples for further details (like dambreak)

When you have fixed that let us know what the results are!

regards
Colin

sdivye February 22, 2015 05:39

FOAM FATAL ERROR: face 0 in patch 4 does not have neighbour cell face: 4(2 4 14 12)
 
Hello Everyone
I am very new to OF and trying to mesh. I am getting following error:

FOAM FATAL ERROR:
face 0 in patch 4 does not have neighbour cell face: 4(2 4 14 12)

From function polyMesh::facePatchFaceCells(const faceList& patchFaces,const labelListList& pointCells,const faceListList& cellsFaceShapes,const label patchID)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 127.

I even tried with "paraFoam -block", but that too gave an error.
Please help me rectify this error

My blockMeshDict is:

convertToMeters 1;

vertices
(
(0 0 0) // 0
(1.5 0 0) //1
(2 0.5 0) //2
(0 0.5 0) //3
(2.5 0 0) //4
(4 0 0) //5
(4 0.5 0) //6
(0 2 0) //7
(2 2 0) //8
(4 2 0) //9
(0 0 0.1) // 10
(1.5 0 0.1) //11
(2 0.5 0.1) //12
(0 0.5 0.1) //13
(2.5 0 0.1) //14
(4 0 0.1) //15
(4 0.5 0.1) //16
(0 2 0.1) //17
(2 2 0.1) //18
(4 2 0.1) //19
);

blocks
(
hex (0 1 2 3 10 11 12 13) (10 10 1) simpleGrading (1 1 1)
hex (3 2 8 7 13 12 18 17) (10 10 1) simpleGrading (1 1 1)
hex (4 5 6 7 14 15 16 17) (10 10 1) simpleGrading (1 1 1)
hex (3 6 9 8 13 16 19 18) (10 10 1) simpleGrading (1 1 1)
);

edges
(
arc 1 2 (1.64644 0.35355 0)
arc 2 4 (2.35355 0.35355 0)
arc 11 12 (1.64644 0.35355 0.1)
arc 12 14 (2.35355 0.35355 0.1)
);

boundary
(
down
{
type symmetryPlane;
faces
(
(0 1 11 10)
(4 5 15 14)
);
}
right
{
type patch;
faces
(
(5 6 16 15)
(6 9 19 16)
);
}
up
{
type symmetryPlane;
faces
(
(9 8 18 19)
(8 7 17 18)
);
}
left
{
type patch;
faces
(
(3 0 10 13)
(7 3 13 17)
);
}
cylinder
{
type symmetry;
faces
(
(2 4 14 12)
(1 2 12 11)
);
}
);

mergePatchPairs
(
);

alexB March 6, 2015 07:01

Hi Divye,

check your points, is (2 4 14 12) really a face of your boundary?
there is no hex with all the 4 points.

So my suggestion is that (2 4 14 12) is just an imaginary plane between two hexes.

greetings
Alex

curiosity July 1, 2015 08:33

Issue with multiple blocks of different sizes
 
1 Attachment(s)
Hello,

Iīve found the same error of the title of the thread:

FOAM FATAL ERROR:
face 0 in patch 9 does not have neigbour cell face


...And I know the reason thanks to reading this forum :)

However, I do not know how to solve it in a "fast way".

Attached it is a picture of my geometry: two solid blocks and two air blocks of different sizes.

And my question is: what should I do with the face Air1_to_Solid2? That is, with the portion of the Solid2 face coupled with the block called Air1?

Should I make a new block dividing the block Air1 so that the face Air1_to_solid2 belongs to that new block? This is probably the answer, but is there any faster way to debug the error? The utility mergepatchpairs maybe?

Thanks in advance

alexB July 1, 2015 17:24

Hi Paula,

First of all... what are you trying to simulate?

Secondly, to help you, at least the blockmeshdict would be useful too.

kind regards
Alex

curiosity July 1, 2015 17:39

Hello AlexB,

Sorry, I didnīt give information enough :S

Iīm doing this tutorial: http://es.slideshare.net/ARPITSINGHA...regiontutorial

Iīm having troubles with the boundary Air1_to_Solid2, as it doesnīt belong to an independent block.

I havenīt the blockMeshDict here, but I can attach it tomorrow :D

Thanks for your reply!

curiosity July 2, 2015 03:07

1 Attachment(s)
Hello Alexander,

This is the blocMeshDict file:
HTML Code:

convertToMeters 1;

vertices
(
//Solid1
    (0 0 0) // 0
    (10 0 0) // 1=8
    (10 0.3 0) // 2=17
    (0 0.3 0) // 3
    (0 0 1) // 4
    (10 0 1) // 5=12
    (10 0.3 1) // 6=21
    (0 0.3 1) // 7=20
//Solid2
    (10 0 0) // 8=1
    (20 0 0) // 9
    (20 0.5 0) //10=25
    (10 0.5 0) // 11=24
    (10 0 1) //12=5
    (20 0 1) //13
    (20 0.5 1) // 14=29
    (10 0.5 1) // 15=28
//Air1
    (0 0.3 0) // 16=3
    (10 0.3 0) // 17=2
    (10 1 0) // 18=27
    (0 1 0) // 19
    (0 0.3 1) //20=7
    (10 0.3 1) //21=6
    (10 1 1) //22=31
    (0 1 1) //23
//Air2
    (10 0.5 0) // 24=11
    (20 0.5 0) // 25=10
    (20 1 0) // 26
    (10 1 0) // 27=18
    (10 0.5 1) //28=15
    (20 0.5 1) //29
    (20 1 1) //30
    (10 1 1) //31=22
);



blocks
(
//Solid1
    hex (0 1 2 3 4 5 6 7) Solid1 (10 10 10) simpleGrading (1 1 1)
//Solid2
    hex (8 9 10 11 12 13 14 15) Solid2 (10 10 10) simpleGrading (1 1 1)
//Air1
    hex (16 17 18 19 20 21 22 23) Air1 (10 10 10) simpleGrading (1 1 1)
//Air2
    hex (24 25 26 27 28 29 30 31) Air2 (10 10 10) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    top
    {
        type wall;
        faces
        (
            (4 5 6 7) //counter clockwise from the outside of the block
            (12 13 14 15)
            (28 29 30 31)
            (20 21 22 23)
        );
    }
    bottom
    {
        type wall;
        faces
        (
            (3 2 1 0) //counter clockwise from the outside of the block
            (11 10 9 8)
            (19 18 17 16)
            (27 26 25 24)
        );
    }   
    rear
    {
        type wall;
        faces
        (
            (0 4 7 3) //counter clockwise from the outside of the block
            (16 20 23 19)
        );
    }
    left
    {
        type wall;
        faces
        (
            (18 19 23 22) //clockwise from the inside of the block
            (26 27 31 30)
        );
    } 
    right
    {
        type wall;
        faces
        (
            (0 1 5 4) //counter clockwise from the outside of the block
            (8 9 13 12)
        );
    }   
    front
    {
        type wall;
        faces
        (
            (9 10 14 13) //counter clockwise from the outside of the block
            (25 26 30 29)
        );
    }
    Air1_to_Solid1
    {
        type wall;
        faces
        (
            (16 17 21 20) //counter clockwise from the outside of the block
        );
    }
    Air2_to_Solid2
    {
        type wall;
        faces
        (
            (24 25 29 28) //counter clockwise from the outside of the block
        );
    }
    Air2_to_Air1
    {
        type wall;
        faces
        (
            (24 28 31 27) //counter clockwise from the outside of the block
        );
    }
    Air1_to_Solid2
    {
        type wall;
        faces
        (
            (17 11 15 21) //counter clockwise from the outside of the block
        );
    }
    Solid1_to_Solid2
    {
        type wall;
        faces
        (
            (1 2 6 5) //counter clockwise from the outside of the block
        );
    }

);
mergePatchPairs
(
);

// ************************************************************************* //

Attached you can also find the picture with the labels of the vertices in each block.

Thanks!

zfaraday July 2, 2015 10:32

Hi Paula,

after a quick glance to your blockMeshDict file I think the problem is that you need to use a couple of extra blocks in order to define properly your geometry. Besides that, it makes no sense to use the same point twice with a diferent label when you define all the vertices, i.e., the point 2 is tha same than the point 17... Define each point only once in the vertices subdictionary! Your vertices list will get muc shorter.

I would define your geometry by using 24 vertices and 6 blocks instead of 4, something like shown below would work.
Code:

only the vertices from the bottom of the geometry are shown

    9  ______ 10 _______ 11
    6 |_______7_|_______| 8
    3 |_______4_|_______| 5
    0 |_______1_|_______| 2

With this geometry the set of blocks would be the ones shown below

Code:

(...)
// the block divisions should be checked in order to be consistent, I only copy-pasted the blocks subdict from the previous post

blocks
(
//Solid1
    hex (0 1 4 3 12 13 16 15) Solid1 (10 10 10) simpleGrading (1 1 1)
//Solid2
    hex (1 2 5 4 12 14 17 16) Solid2 (10 10 10) simpleGrading (1 1 1)
    hex (4 5 8 7 16 17 20 19) Solid2 (10 10 10) simpleGrading (1 1 1)
//Air1
    hex (3 4 7 6 x x x x) Air1 (10 10 10) simpleGrading (1 1 1)
    hex (6 7 10 9 x x x x) Air1 (10 10 10) simpleGrading (1 1 1)
//Air2
    hex (7 78 11 10 x x x x) Air2 (10 10 10) simpleGrading (1 1 1)
);
(...)

Hope it helps!

Best regards,

Alex

curiosity July 3, 2015 01:40

Quote:

Originally Posted by zfaraday (Post 553008)

With this geometry the set of blocks would be the ones shown below

Code:

(...)
// the block divisions should be checked in order to be consistent, I only copy-pasted the blocks subdict from the previous post

blocks
(
//Solid1
    hex (0 1 4 3 12 13 16 15) Solid1 (10 10 10) simpleGrading (1 1 1)
//Solid2
    hex (1 2 5 4 12 14 17 16) Solid2 (10 10 10) simpleGrading (1 1 1)
    hex (4 5 8 7 16 17 20 19) Solid2 (10 10 10) simpleGrading (1 1 1)
//Air1
    hex (3 4 7 6 x x x x) Air1 (10 10 10) simpleGrading (1 1 1)
    hex (6 7 10 9 x x x x) Air1 (10 10 10) simpleGrading (1 1 1)
//Air2
    hex (7 78 11 10 x x x x) Air2 (10 10 10) simpleGrading (1 1 1)
);
(...)

Hope it helps!

Best regards,

Alex

Oh, sure it helps! That I what I was looking for! I didnīt know if I could define two blocks of hex for the same zone and now I see that I can. Thank you very much! :)


Iīve been reading the user guide and Iīd like to ask you another question:

If I do not use the utility mergePatchPairs would it be the "face matching" method and if I define the patches of type Air_to_solid as boundary and then use mergePatchPairs utility, would I be using the face merging method?

To clarify my question, an example of the code for only one external face (top) and one internal face (Solid_to_Air), for the face merging method:

HTML Code:

boundary
(
    top
    {
        type wall;
        faces
        (
            (X X X X)
            (X X X X)
            (X X X X)
            (X X X X)
        );
    }
...
...
...
    Solid1_to_Air1
    {
        type patch;
        faces
        (
            (X X X X)
        );
    }
    Air1_to_Solid1
    {
        type patch;
        faces
        (
            (X X X X)
        );
    }
...
...
...

);
mergePatchPairs
(
            (Solid1_to_Air1 Air1_to_Solid1)
);

// ************************************************************************* //

Maybe Iīm mixing concepts... but I fear that I donīt understand when is the utility "mergePatchPairs" necessary :confused:

Thanks again!

zfaraday July 3, 2015 08:04

Hi Paula,

Forget about the coupling patches between regions at this stage, you only need to define the external ones in blockMesh, the other ones get defined after you execute splitMeshRegions and changeDictionary utilities.

I see you are a bit lost when it comes to multi region cases, you need to take a look at the planeWall2D tutorial and try to understand the steps to reproduce a multi region case. I strongly recomend you to study this case before setting up any other multi region case on your own.

Best regards,

Alex

curiosity July 6, 2015 09:03

2 Attachment(s)
Hello Alex,

Iīve made the tutorial you suggested (planeWall2D) and some others by Iīm afraid they havenīt clarify some of my doubts.

Iīve investigated which are the files OpenFOAM ask you for, in order to set a multiRegion case step by step (without using Allrun).
Attached you can find a scheme of the process.

As you can see, my plan is not to use the ./Allrun command. What I want to do is to enter the following commands one by one:

1. blockMesh
2. topoSet
3. splitMeshRegions
4. script to remove fluid fields from solid regions (contained in the Allrun file)
5. chtMultiregionFoam

However, I do not know how to set the changeDictionaryDict and it doesnīt help me so much to look into the example of the tutorial.

Please, could you help me?

zfaraday July 6, 2015 11:22

Quote:

Originally Posted by curiosity (Post 553999)
Hello Alex,

Iīve made the tutorial you suggested (planeWall2D) and some others by Iīm afraid they havenīt clarify some of my doubts.

Iīve investigated which are the files OpenFOAM ask you for, in order to set a multiRegion case step by step (without using Allrun).
Attached you can find a scheme of the process.

As you can see, my plan is not to use the ./Allrun command. What I want to do is to enter the following commands one by one:

1. blockMesh
2. topoSet
3. splitMeshRegions
4. script to remove fluid fields from solid regions (contained in the Allrun file)
5. chtMultiregionFoam

However, I do not know how to set the changeDictionaryDict and it doesnīt help me so much to look into the example of the tutorial.

Please, could you help me?

Hi Paula,

I don't like the fact that you don't want to use the Allrun script since it is totally helpful in any case you set up. I strongly suggest you to get used to make use of them usually!

I don't understand either what you say you don't understand...:confused: I mean, the use of changeDictionary is quite obvious and it is explained in the planeWall2D tutorial

Quote:

system - In here are the files that configure or help configure the case.
  • controlDict - Here you control the number of iterations and so on...
  • topoSetDict - Here you defined the extents of each region. It's here that you can select which cells of the mesh that belong in each region.
  • bottomAir - In here are the runtime configurations for the bottom region.
    • changeDictionaryDict - This is the main file where you define the actual boundary conditions to be used for the bottom region.
  • topAir - In here are the runtime configurations for the top region.
    • changeDictionaryDict - This is the main file where you define the actual boundary conditions to be used for the top region.
  • wall - In here are the runtime configurations for the wall region.
    • changeDictionaryDict - This is the main file where you define the actual boundary conditions to be used for the wall region.

You only need to copy the changeDictionaryDict file from the tutorial case to your own case and adapt it to suit your requirements!

Hope it helps!

Alex

curiosity July 7, 2015 02:43

Hello Alex,

Of course the Allrun script is useful. I was only entering command by command in order to know what does it execute internally when I enter: ./Allrun in the console. I mean, what folders are created when I enter blockMesh and then SplitMeshRegions, for example. But Iīll use it, of course! :)

Regarding to the changeDictionary of the tutorial WallPlane2D, what I donīt understand is (comments/questions in the own code):

HTML Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      changeDictionaryDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dictionaryReplacement
{
    boundary
    {
        bottomAir_bottom
        {
            type            symmetryPlane;
        }
    }

    U
    {
        internalField  uniform (0.01 0 0);

        boundaryField
        {
            leftLet
            {
                type            fixedValue;
                value          uniform ( 0.1 0 0 );
            }

            rightLet
            {
                type            inletOutlet;
                inletValue      uniform ( 0 0 0 );
                value          uniform ( 0.1 0 0 );
            }

            bottomAir_bottom
            {
                type            symmetryPlane;
            }
           
            "bottomAir_to_.*"
            {
                type            fixedValue;
                value          uniform (0 0 0);
            }
        }
    }

    T
    {
        internalField  uniform 500;

        boundaryField
        {
            leftLet
            {
                type            fixedValue;
                value          uniform 500;
            }

            rightLet
            {
                type            zeroGradient;
            }

            bottomAir_bottom
            {
                type            symmetryPlane;
            }

            "bottomAir_to_.*"
            {
                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                Tnbr            T;
                kappa          fluidThermo;
                kappaName      none;
                value          uniform 300;
            }
        }
    }

    epsilon
    {
        internalField  uniform 0.01;

        boundaryField
        {
            leftLet
            {
                type            inletOutlet;
                inletValue      uniform 0.01;
                value          uniform 0.01;
            }

            rightLet
            {
                type            fixedValue;
                value          uniform 0.01;
            }

            bottomAir_bottom
            {
                type            symmetryPlane;
            }

            "bottomAir_to_.*"
            {
                type            compressible::epsilonWallFunction;
                value          uniform 0.01;
            }
        }
    }

    k
    {
        internalField  uniform 0.1;

        boundaryField
        {
            leftLet
            {
                type            inletOutlet;
                inletValue      uniform 0.1;
                value          uniform 0.1;
            }

            rightLet
            {
                type            fixedValue;
                value          uniform 0.1;
            }

            bottomAir_bottom
            {
                type            symmetryPlane;
            }

            "bottomAir_to_.*"
            {
                type            compressible::kqRWallFunction;
                value          uniform 0.1;
            }
        }
    }

    p_rgh
    {
        internalField  uniform 1e5;

        boundaryField
        {
            ".*"  -----------------------------------> what does this mean?
            {
                type            fixedFluxPressure;
                value          uniform 1e5;
            }

            bottomAir_bottom
            {
                type            symmetryPlane;
            }

            leftLet
            {
                type            fixedValue;
                value          uniform 1e5;
            }

            frontAndBack -------------> Why doesnīt "frontAndBack" appear for other variables: epsilon, k, T and U definition?
            {
                type            empty;
            }
         
        }
    }

    p
    {
        internalField  uniform 1e5;

        boundaryField
        {
            ".*"
            {
                type            calculated;
                value          uniform 1e5;
            }

            bottomAir_bottom
            {
                type            symmetryPlane;
            }

            leftLet
            {
                type            calculated;
                value          uniform 1e5;
            }

            frontAndBack
            {
                type            empty;
            }
         
        }
    }
}

// ************************************************************************* //

Thank you!

zfaraday July 7, 2015 05:39

Hi,

".*" is a wildcard that means "any character chain". Thus, "bottomAir_to_.*" refers to any coupling boundary condition for the region bottomAir to any solid; ".*" refers to any boundary condition found in the geometry.

Hope it helps!

curiosity July 7, 2015 05:43

Ok!

I think that little by little Iīm solving the mistery ;)

Thank you very much Alex!


All times are GMT -4. The time now is 01:20.