[Forbraer]

Hello everyone.

I'm trying to simulate the Bird's problem in section of 13.4, that must use the dsmcFoam solver.

But, in this case we must build a wall, something linke that:


|--------------|
|____|_____|
The thin wall is in the bottom middle. Image the hifens as the domain roof.

I'm new in OpenFoam and I noted that it have a problem in declaring internal faces, then I created four points twice declared as walls1 and walls2, as shown below.

vertices
(
    (-5 0 0)  // 0
    (0 0 0)   // 1
    (0 10 0)  // 2
    (-5 10 0) // 3 
    (-5 0 36) // 4
    (0 0 36)  // 5
    (0 10 36) // 6
    (-5 10 36)// 7
    (-5 0 18) // 8
    (0 0 18)  // 9
    (0 10 18) // 10
    (-5 10 18)// 11
    (20 0 36) // 12
    (20 10 36)// 13
    (20 0 18) // 14
    (20 10 18)// 15
    (20 0 0)  // 16
    (20 10 0) // 17
    
    (0 0 0)   // 18 -1
    (0 0 18)  // 19 -9
    (0 10 18) // 20 -10
    (0 10 0)  // 21 -2
    

);

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

boundary
(
    flow
    {
        type patch;
        faces
        (
            (4 5 6 7)
            (5 12 13 6)
            (0 8 11 3)
            (8 4 7 11)
            (3 11 10 2)
            (11 7 6 10)
            (6 13 15 20)
            (21 20 15 17)
            (0 8 9 1)
            (8 4 5 9)
            (18 19 14 16)
            (19 5 12 14)
            (12 13 15 14)
            (14 15 17 16)
        );
    }

    plan1
    {
        type symmetryPlane;
        faces
        (
            (0 1 2 3)
        );
    }

    plan2
    {
        type symmetryPlane;
        faces
        (
            (18 16 17 21)
        );
    }

    walls1
    {
        type wall;
        faces
        (
            (2 10 9 1)
        );
    }
    
    walls2
    {
        type wall;
        faces
        (
            (18 19 20 21)
        );
    }
)

The problem runs, blockMesh works, dsmcInitialise works, but the code stops in first "particles inserted" count.

I think is a problem about the "thin wall" inserted.
Can anybody tell me if it is ok? Or, I must change it?