CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[blockMesh] generating additional block

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 2, 2020, 13:42
Default generating additional block
  #1
New Member
 
David
Join Date: Oct 2020
Posts: 21
Rep Power: 5
fidu is on a distinguished road
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:
Code:
    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:
Code:
/*--------------------------------*- 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:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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

Last edited by fidu; November 3, 2020 at 10:11.
fidu is offline   Reply With Quote

Old   November 3, 2020, 11:14
Default Found Solution
  #2
New Member
 
David
Join Date: Oct 2020
Posts: 21
Rep Power: 5
fidu is on a distinguished road
I found the solution, I did not define z_tree in right unit....
fidu is offline   Reply With Quote

Reply

Tags
blockmeshdict, blocks, refine mesh

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[blockMesh] FOAM FATAL ERROR: Inconsistent number of faces blockMesh::createMergeList() line 193 Hengel OpenFOAM Meshing & Mesh Conversion 7 November 15, 2021 22:56
[ICEM] Replay not working properly - Creating the mesh for Free Block. Wingman ANSYS Meshing & Geometry 4 January 17, 2017 03:57
[ICEM] Quest for seemingly simple Blocking (degenerate block) la7low ANSYS Meshing & Geometry 0 July 12, 2011 19:48
[blockMesh] FOAM FATAL ERROR: Inconsistent number of faces blockMesh::createMergeList() line 193 Hengel OpenFOAM Meshing & Mesh Conversion 0 September 15, 2010 08:34
[ICEM] Hexa to Sweep block conversion pertupd ANSYS Meshing & Geometry 1 June 19, 2010 21:37


All times are GMT -4. The time now is 15:18.