[fidu]

Hi foamers

I hope my question is not too trivial but I am quite stuck. I think I did not define the blocks in the right manner. I want to adapt my blockMeshDict in a manner that I split my computational domain in too, in near ground block and in a near ceiling block. In order to that I added the coordinate y_tree which is where I want to split my block in too. My original block was:

    hex(0 1 2 3 4 5 6 7)
    (
        $:backgroundMesh.xCells
        $:backgroundMesh.yCells
        $:backgroundMesh.zCells
    ) 
    simpleGrading 
    (  
        ( 
            (0.109375 0.085106383 0.3)    // l% x-dir, % cells, expansion =  
            (0.53125  0.780141844 1)    // l% x-dir, % cells, expansion =  
            (0.359375 0.134751773 10) // l% x-dir, % cells, expansion =  
        ) 
        1                  // y-direction expansion ratio
        (                  
            (0.1572376 0.375 1)    // z% x-dir, % cells, expansion =      
            (0.8427624 0.625 8)    // z% x-dir, % cells, expansion =               
        )  
    )

I did adapt it now in the following way:

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

backgroundMesh
{
    xMin   1.9; // L = 
    xMax   8.3;
    yMin   0.0054; // L = 
    yMax   1.9;
    zMin     0;   //-0.01875
    zMax   1.61675;   //1.6
    xCells  141;    //141   800
    yCells  62;     //62  200
    zCells  32;    //32   200

    zTree   15;

}

//convertToMeters 1;

vertices
(
    ($:backgroundMesh.xMin $:backgroundMesh.yMin $:backgroundMesh.zMin)  //0
    ($:backgroundMesh.xMax $:backgroundMesh.yMin $:backgroundMesh.zMin) //1
    ($:backgroundMesh.xMax $:backgroundMesh.yMax $:backgroundMesh.zMin) //2
    ($:backgroundMesh.xMin $:backgroundMesh.yMax $:backgroundMesh.zMin) //3

    ($:backgroundMesh.xMin $:backgroundMesh.yMin $:backgroundMesh.zMax) //4
    ($:backgroundMesh.xMax $:backgroundMesh.yMin $:backgroundMesh.zMax) //5
    ($:backgroundMesh.xMax $:backgroundMesh.yMax $:backgroundMesh.zMax) //6
    ($:backgroundMesh.xMin $:backgroundMesh.yMax $:backgroundMesh.zMax) //7

    ($:backgroundMesh.xMin $:backgroundMesh.yMin $:backgroundMesh.zTree)    //8    new 4 
    ($:backgroundMesh.xMax $:backgroundMesh.yMin $:backgroundMesh.zTree)    //9     new 5
    ($:backgroundMesh.xMax $:backgroundMesh.yMax $:backgroundMesh.zTree)    //10    new 6
    ($:backgroundMesh.xMin $:backgroundMesh.yMax $:backgroundMesh.zTree)    //11    new 7

);

blocks
(
    hex (0 1 2 3 8 9 10 11)               //0 1 2 3 4 5 6 7   domain of ground to zTree
    (
        $:backgroundMesh.xCells
        $:backgroundMesh.yCells
        $:backgroundMesh.zCells
    ) 
    simpleGrading   (1 1 1)

    hex (8 9 10 11 4 5 6 7)               //0 1 2 3 4 5 6 7 domain of zTree to ceilling   hex (8 9 10 11 4 5 6 7)
    (
        $:backgroundMesh.xCells
        $:backgroundMesh.yCells
        $:backgroundMesh.zCells
    )

    simpleGrading (1 1 1) 

);

edges
(
);

boundary
(
    inlet
    {
        type patch;
        faces
        (
            (0 3 11 8)
            (8 11 7 4)       //(0 3 7 4)
        );
    }

    outlet
    {
        type patch;
        faces
        (
            (1 9 10 2)
            (9 5 6 10)       //(1 5 6 2)
        );
    }

    ground
    {
        type wall;
        faces
        (
            (0 1 2 3)       //(0 1 2 3)
        );
    }

    frontAndBack
    {
        type wall;
        faces
        (
            (0 8 9 1) (8 4 5 9)                 //front (0 4 5 1) 
            (3 2 10 11) (11 10 6 7)             //back (3 2 6 7)
            (4 7 6 5)                           //ceilling (4 7 6 5)
        );
    }

);

mergePatchPairs
(
);

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

When I try to run blockMesh I get the following error:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _b45f8f6f58-20200629 OPENFOAM=2006
Arch   : "LSB;label=32;scalar=64"
Exec   : blockMesh
Date   : Nov 02 2020
Time   : 19:37:38
Host   : LAPTOP-9VDLVQ8A
PID    : 7987
I/O    : uncollated
Case   : /mnt/c/Users/david/Desktop/working_cases/work/porou_media/two_blocks
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Creating block mesh from "system/blockMeshDict"
Creating block edges
No non-planar block faces defined
Creating topology blocks
Creating topology patches

Creating block mesh topology


--> FOAM FATAL ERROR:
face 0 in patch 0 does not have neighbour cell face: 4(0 3 7 4)

    From Foam::labelList Foam::polyMesh::facePatchFaceCells(const faceList&, const labelListList&, const faceListList&, Foam::label) const
    in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 119.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::exitOrAbort(int, bool) at ??:?
#2  Foam::polyMesh::facePatchFaceCells(Foam::List<Foam::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const at ??:?
#3  Foam::polyMesh::setTopology(Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<int>&, Foam::List<int>&, int&, int&, Foam::List<Foam::cell>&) at ??:?
#4  Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Field<Foam::Vector<double> >&&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::PtrList<Foam::dictionary> const&, Foam::word const&, Foam::word const&, bool) at ??:?
#5  Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&) at ??:?
#6  Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&, Foam::blockMesh::mergeStrategy) at ??:?
#7  ? at ??:?
#8  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#9  ? at ??:?
Aborted (core dumped)

Which I intepreted that the order of the second block was not defined in the right direction/order. However I dont seem to achieve the right order. Could you help?

Thanks a lot for your help and thanks a lot for your time.

Best
fiud

[fidu]

I found the solution, I did not define z_tree in right unit....