CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Runing InterFoam for 3 D damBreak case Patch issue

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 13, 2007, 09:32
Default Hi .. i am trying to run da
  #1
Senior Member
 
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 165
Rep Power: 8
nishant_hull is on a distinguished road
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:stream&) 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.
blockMeshDict
Can anyone please look into it and let me know of the problem. Looking forward for your suggestions in this regard..

Nishant
__________________
Thanks and regards,

Nishant
nishant_hull is offline   Reply With Quote

Old   December 13, 2007, 12:19
Default CHECKMESH provide the followin
  #2
Senior Member
 
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 165
Rep Power: 8
nishant_hull is on a distinguished road
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
__________________
Thanks and regards,

Nishant
nishant_hull is offline   Reply With Quote

Old   December 13, 2007, 12:43
Default Yes - reported non-orthogonali
  #3
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
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
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   December 15, 2007, 13:17
Default Yea. There was kinda basic mis
  #4
Senior Member
 
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 165
Rep Power: 8
nishant_hull is on a distinguished road
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!!
__________________
Thanks and regards,

Nishant
nishant_hull is offline   Reply With Quote

Old   February 18, 2009, 01:25
Default Hi Nishant, I am also getting
  #5
New Member
 
Sayantan Datta Gupta
Join Date: Mar 2009
Location: Chennai, Tamilnadu, India
Posts: 16
Rep Power: 8
rishi123 is on a distinguished road
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
(
);

// ************************************************** *********************** //
__________________
Rishi
rishi123 is offline   Reply With Quote

Old   February 18, 2009, 01:32
Default Hi Nishant, I am also getting
  #6
New Member
 
Sayantan Datta Gupta
Join Date: Mar 2009
Location: Chennai, Tamilnadu, India
Posts: 16
Rep Power: 8
rishi123 is on a distinguished road
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
(
);

// ************************************************** *********************** //
__________________
Rishi
rishi123 is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Interfoam Droplet under shear test case adona058 OpenFOAM Running, Solving & CFD 3 May 3, 2010 18:46
Help on interFoam writeControl in damBreak tutorial asaha OpenFOAM Running, Solving & CFD 1 February 2, 2009 08:36
ParaFOAM issue with cavity case mschoenberg OpenFOAM Paraview & paraFoam 3 November 21, 2008 04:57
Cylindrical Patch for interFoam vatant OpenFOAM Running, Solving & CFD 3 December 12, 2006 05:08
Problems starting new case interFoam billy OpenFOAM Running, Solving & CFD 3 June 21, 2006 10:18


All times are GMT -4. The time now is 11:33.