CFD Online Discussion Forums - Runing InterFoam for 3 D damBreak case Patch issue

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Runing InterFoam for 3 D damBreak case Patch issue (https://www.cfd-online.com/Forums/openfoam-solving/57946-runing-interfoam-3-d-dambreak-case-patch-issue.html)

Hi .. i am trying to run da

Hi ..

i am trying to run damBreak case for 3 D case and My blockmesh looks ok with the geometry but there are some issue with the patch. I m getting some warning related to patch. I have taken the boundary conditions of side walls (which was empty-type in 2 D case and run well on my machine) similar to the rightwall boundary condition.i.e. Zerogradient/u=0.

The error reported as:

[343880@w191-210 interFoam]$ interFoam . damBreak-expt
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : interFoam . damBreak-expt
Date : Dec 13 2007
Time : 13:25:31
Host : w191-210
PID : 6147
Root : /home/343880/OpenFOAM/343880-1.4.1/run/tutorials/interFoam
Case : damBreak-expt
Nprocs : 1
Create time

Create mesh for time = 0

Reading environmentalProperties
Reading field pd

Reading field gamma

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Calculating field g.h

time step continuity errors : sum local = 0, global = 0, cumulative = 0
DICPCG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: -0

Starting time loop

Courant Number mean: 0 max: -0
deltaT = 0.000595238
Time = 0.000595238

MULES: Solving for gamma
Liquid phase volume fraction = 0.232558 Min(gamma) = 0 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.232558 Min(gamma) = 0 Max(gamma) = 1
#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0x110420]
#3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::multiply<foam::fvspatchfield,>(Foam::Geometr icField<double,>&, Foam::GeometricField<double,> const&, Foam::GeometricField<double,> const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/interFoam"
#5 Foam::tmp<foam::geometricfield<double,> > Foam::operator*<foam::fvspatchfield,>(Foam::tmp<fo am::geometricfield<double,> > const&, Foam::tmp<foam::geometricfield<double,> > const&) in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/interFoam"
#6 main in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/interFoam"
#7 __libc_start_main in "/lib/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/343880/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/interFoam"
Floating point exception

Also the WARNINGS with blockmesh utility is like :
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 134
zero or negative pyramid volume: -0.013065 for face 4
--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 134
zero or negative pyramid volume: -0.013065 for face 5
--> FOAM Warning :
From function blockMesh::createTopology(IOdictionary& meshDescription)
in file createTopology.C at line 391
negative volume block : 22, probably defined inside-out

--> FOAM Warning :
From function cellModel::mag(const labelList&, const pointField&)
in file meshes/meshShapes/cellModel/cellModel.C at line 134
zero or negative pyramid volume: -0.015075 for face 5
--> FOAM Warning :
From function blockMesh::createTopology(IOdictionary& meshDescription)
in file createTopology.C at line 391
negative volume block : 34, probably defined inside-out

Default patch type set to empty

Check block mesh topology

Basic statistics
Number of internal faces : 70
Number of boundary faces : 70
Number of defined boundary faces : 70
Number of undefined boundary faces : 0

Checking patch -> block consistency

Creating block offsets

Creating merge list .

Creating points

Creating cells

Creating patches

Creating mesh from block mesh

Default patch type set to empty

Creating merge patch pairs

Writing polyMesh

end

BlockMeshDict file is attached.
http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif blockMeshDict
Can anyone please look into it and let me know of the problem. Looking forward for your suggestions in this regard..

Nishant

CHECKMESH provide the followin

CHECKMESH provide the following results:-

[343880@w191-210 interFoam]$ checkMesh . damBreak-expt
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : checkMesh . damBreak-expt
Date : Dec 13 2007
Time : 16:10:42
Host : w191-210
PID : 6566
Root : /home/343880/OpenFOAM/343880-1.4.1/run/tutorials/interFoam
Case : damBreak-expt
Nprocs : 1
Create time

Create polyMesh for time = constant

Time = constant

Mesh stats
points: 34481
edges: 99820
faces: 96300
internal faces: 89460
cells: 30960
boundary patches: 4
point zones: 0
face zones: 0
cell zones: 0

Number of cells of each type:
hexahedra: 30960
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface
leftWall 540 589 ok (not multiply connected)
rightWall 2628 2793 ok (not multiply connected)
lowerWall 1932 2021 ok (not multiply connected)
atmosphere 1740 1829 ok (not multiply connected)

Checking geometry...
Domain bounding box: (0 0 0) (0.805 0.25 0.25)
Boundary openness (-9.05278e-17 -1.4379e-15 1.01299e-15) OK.
***High aspect ratio cells found, Max aspect ratio: 1.31439e+196, number of cells 30960
<<Writing 30960 cells with high aspect ratio to set highAspectRatioCells
Minumum face area = 4.44444e-05. Maximum face area = 0.000285938. Face area magnitudes OK.
Min volume = 2e-300. Max volume = 2e-300. Total volume = 6.192e-296. Cell volumes OK.
Mesh non-orthogonality Max: 180 average: 180
***Number of non-orthogonality errors: 89460.
<<Writing 89460 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 185760 faces are incorrectly oriented.
<<Writing 96300 faces with incorrect orientation to set wrongOrientedFaces
Max skewness = 6.66134e-14 OK.
Min/max edge length = 0.00666667 0.01875 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1
All face flatness OK.

Failed 3 mesh checks.

End

Can anybody tell me, where am i going wrong??

Nishant

Yes - reported non-orthogonali

Yes - reported non-orthogonality is 180 deg, and it should be (well) below 90 deg. This mesh is wrong - did you actually look at it?

Hrv

Yea. There was kinda basic mis

Yea. There was kinda basic mistake of not taking x-y-z system in proper order. now my case is running fine.
Thank you hrv!
cheers!!

Hi Nishant, I am also getting

Hi Nishant,
I am also getting the same error as you reported in cavity problem with icoFoam. Your previous post tells about "xyzsystem" can you explain the solution.
Bye
Rishi
P.S: I am attaching my blockmeshdict file content

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

convertToMeters 0.1;

vertices
(
(0 0 0)
(0.5 0 0)
(1 0 0)
(0 -0.5 0)
(0.5 -0.5 0)
(1 -0.5 0)
(0.5 -1 0)
(1 -1 0)
(0 0 0.1)
(0.5 0 0.1)
(1 0 0.1)
(0 -0.5 0.1)
(0.5 -0.5 0.1)
(1 -0.5 0.1)
(0.5 -1 0.1)
(1 -1 0.1)
);

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

edges
(
);

patches
(
wall movingWall
(
(0 8 9 1)
(1 9 10 2)
)
wall fixedWalls
(
(0 3 11 8)
(3 4 12 11)
(2 10 13 5)
(4 12 14 6)
(6 7 15 14)
(5 13 15 7)
)
empty frontAndBack
(
(0 1 4 3)
(1 2 5 4)
(4 5 7 6)
(8 9 12 11)
(9 10 13 12)
(12 13 15 14)
)
);

mergePatchPairs
(
);

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

Hi Nishant, I am also getting