|
[Sponsors] |
Runing InterFoam for 3 D damBreak case Patch issue |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 13, 2007, 09:32 |
Hi ..
i am trying to run da
|
#1 |
Senior Member
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 166
Rep Power: 17 |
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 |
|
December 13, 2007, 12:19 |
CHECKMESH provide the followin
|
#2 |
Senior Member
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 166
Rep Power: 17 |
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 |
|
December 13, 2007, 12:43 |
Yes - reported non-orthogonali
|
#3 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,907
Rep Power: 33 |
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 |
|
December 15, 2007, 13:17 |
Yea. There was kinda basic mis
|
#4 |
Senior Member
Nishant
Join Date: Mar 2009
Location: Glasgow, UK
Posts: 166
Rep Power: 17 |
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 |
|
February 18, 2009, 01:25 |
Hi Nishant,
I am also getting
|
#5 |
New Member
Sayantan Datta Gupta
Join Date: Mar 2009
Location: Chennai, Tamilnadu, India
Posts: 16
Rep Power: 17 |
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 |
|
February 18, 2009, 01:32 |
Hi Nishant,
I am also getting
|
#6 |
New Member
Sayantan Datta Gupta
Join Date: Mar 2009
Location: Chennai, Tamilnadu, India
Posts: 16
Rep Power: 17 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Interfoam Droplet under shear test case | adona058 | OpenFOAM Running, Solving & CFD | 3 | May 3, 2010 19:46 |
Help on interFoam writeControl in damBreak tutorial | asaha | OpenFOAM Running, Solving & CFD | 1 | February 2, 2009 08:36 |
[OpenFOAM] ParaFOAM issue with cavity case | mschoenberg | ParaView | 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 11:18 |