[kelvinaaa]

Hello all,

I am changing the damBreak 2D tutorial to 3D case in order to perform analyses of water impacting a barrier.

I have modified the blockMeshDict file with the script shown as follows:

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

convertToMeters 1;

vertices
(
(0 0 0) //0
(1.25 0 0) //1
(1.25 0.8 0) //2
(0 0.8 0) //3
(0 0 0.2) //4
(1.25 0 0.2) //5
(1.25 0.8 0.2) //6
(0 0.8 0.2) //7
(1.43 0.0405 0) //8
(1.43 0.0405 0.2) //9
(1.637 0.144 0) //10
(1.65 0.15 0) //11
(1.551 0.352 0) //12
(1.541 0.347 0) //13
(1.637 0.144 0.2) //14
(1.65 0.15 0.2) //15
(1.551 0.352 0.2) //16
(1.541 0.347 0.2) //17
(1.532 0.365 0) //18
(1.56 0.378 0) //19
(1.56 0.8 0) //20
(1.322 0.8 0) //21
(1.532 0.365 0.2) //22
(1.56 0.378 0.2) //23
(1.56 0.8 0.2) //24
(1.322 0.8 0.2) //25
(1.745 0 0) //26
(1.745 0.8 0) //27
(1.745 0 0.2) //28
(1.745 0.8 0.2) //29
(2.5 0 0) //30
(2.5 0.8 0) //31
(2.5 0 0.2) //32
(2.5 0.8 0.2) //33
(1.373 0.265 0) //34
(1.369 0.285 0) //35
(1.25 0.285 0) //36
(1.25 0.265 0) //37
(0 0.285 0) //38
(0 0.265 0) //39
(1.373 0.265 0.2) //40
(1.369 0.285 0.2) //41
(1.25 0.285 0.2) //42
(1.25 0.265 0.2) //43
(0 0.285 0.2) //44
(0 0.265 0.2) //45
);

blocks
(
hex (0 1 37 39 4 5 43 45) (50 50 1) simpleGrading (1 1 1)
hex (39 37 36 38 45 43 42 44) (50 50 1) simpleGrading (1 1 1)
hex (38 36 2 3 44 42 6 7) (50 50 1) simpleGrading (1 1 1)
hex (1 8 34 37 5 9 40 43) (50 50 1) simpleGrading (1 1 1)
hex (37 34 35 36 43 40 41 42) (50 50 1) simpleGrading (1 1 1)
hex (36 35 2 2 42 41 6 6) (50 50 1) simpleGrading (1 1 1)
hex (8 10 13 34 9 14 17 40) (50 50 1) simpleGrading (1 1 1)
hex (34 13 18 35 40 17 22 41) (50 50 1) simpleGrading (1 1 1)
hex (35 18 21 2 41 22 25 6) (50 50 1) simpleGrading (1 1 1)
hex (10 11 12 13 14 15 16 17) (50 50 1) simpleGrading (1 1 1)
hex (18 19 20 21 22 23 24 25) (50 50 1) simpleGrading (1 1 1)
hex (19 26 27 20 23 28 29 24) (50 50 1) simpleGrading (1 1 1)
hex (26 30 31 27 28 32 33 29) (50 50 1) simpleGrading (1 1 1)

);

edges
(
arc 1 8 (1.34 0.0102 0)
arc 5 9 (1.34 0.0102 0.2)
);

boundary
(
leftWall
{
type wall;
faces
(
(0 1 37 39)
(39 37 36 38)
(38 36 2 3)
(1 8 34 37)
(37 34 35 36)
(36 35 2 2)
(8 10 13 34)
(34 13 18 35)
(35 18 21 2)
(10 11 12 13)
(18 19 20 21)
(19 26 27 20)
(26 30 31 27)
);
}
rightWall
{
type wall;
faces
(
(4 5 43 45)
(45 43 42 44)
(44 42 6 7)
(5 9 40 43)
(43 40 41 42)
(42 41 6 6)
(9 14 17 40)
(40 17 22 41)
(41 22 25 6)
(14 15 16 17)
(22 23 24 25)
(23 28 29 24)
(28 32 33 29)
);
}
lowerWall
{
type wall;
faces
(
(0 1 5 4)
(1 8 9 5)
(8 10 14 9)
(10 11 15 14)
(18 13 17 22)
(18 19 23 22)
(19 26 28 23)
(26 30 32 28)
);
}
endWall
{
type wall;
faces
(
(0 4 45 39)
(39 45 44 38)
(38 44 7 3)
);
}
frontWall
{
type wall;
faces
(
(30 32 33 31)
);
}
barrier
{
type wall;
faces
(
(11 12 16 15)
);
}
deflector
{
type wall;
faces
(
(12 13 17 16)
);
}
atmosphere
{
type patch;
faces
(
(3 2 6 7)
(2 21 25 6)
(21 20 24 25)
(20 27 29 24)
(27 31 33 29)
);
}
);

mergePatchPairs
(
);

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

And I have encountered the following problem after typing blockMesh

/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 7-ca808c8420bf
Exec : blockMesh
Date : Dec 03 2019
Time : 00:25:52
Host : "kelvin-VirtualBox"
PID : 15387
I/O : uncollated
Case : /home/kelvin/OpenFOAM/kelvin-7/run/damBreak
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Deleting polyMesh directory
"/home/kelvin/OpenFOAM/kelvin-7/run/damBreak/constant/polyMesh"
Creating block mesh from
"/home/kelvin/OpenFOAM/kelvin-7/run/damBreak/system/blockMeshDict"
Creating block edges
No non-planar block faces defined
Creating topology blocks
Creating topology patches

Creating block mesh topology
--> FOAM Warning :
From function Foam:olyMesh:olyMesh(const Foam::IOobject&, Foam:ointField&&, const cellShapeList&, const faceListList&, const wordList&, const Foam::PtrList<Foam::dictionary>&, const Foam::word&, const Foam::word&, bool)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 873
Found 1 undefined faces in mesh; adding to default patch.

Can all help to debug the problem and how to solve it?
Many Thanks for all kindly help.

Regards,
Kelvin

[simrego]

Hi!


It is not a big problem. You have missed a boundary face in the boundary ( ... ).
Open your mesh and look for the patch called "defaultPatches" (I think this is the default name). This is the face what you have missed.
Or you can open your mesh with "paraFoam -block" and you will see the vertices with their numbers in the block mesh.


But you can use this feature for example with 2D geometry. Just skip the front and the back faces, and they will be added to the "defaultFaces" patch automatically.