View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0001544 | OpenFOAM | Bug | public | 2015-02-20 15:46 | 2015-06-12 08:54 |
Reporter | michael2015 | Assigned To | henry | ||
Priority | normal | Severity | minor | Reproducibility | always |
Status | resolved | Resolution | fixed | ||
Summary | 0001544: findElements::findElements cannot find point inside domain | ||||
Description | Hello. I have a listo of 45 points which are inside of the domain. the function gives me an error form some of them: --> FOAM Warning : From function findElements::findElements(const fvMesh&) in file probes/probes.C at line 111 Did not find location (75746 23540 119.025) in any cell. Skipping location. | ||||
Tags | No tags attached. | ||||
|
Which version of OpenFOAM are you using? Are some of the point on or very near the boundary? How do we reproduce the error? Is the case available somewhere? |
|
the points are closed to a wall but there are alwayes some cell layers between the points and the wall. The version is 2.3.0. Unfortunatly I cannot send you the case because it is industrial relevant. But I saw resent similar bug reports: http://www.openfoam.org/mantisbt/view.php?id=1349 |
|
Do you get the same behavior with OpenFOAM-2.3.x? Are the points which cannot be found colocated with mesh vertices, edges or faces? |
|
Hello! I installed OpenFOAM-2.3.x and the same bug occured. Some points which are not found are close to faces bot some are well in the middle of the cell. Furthermore I relized that the cellRepresentation is of type FACEDIAGTETS. So I guess something goes wrong in the pointInCell function in polyMesh.C for the case FACEDIAGTETS |
|
Have you checked the mesh quality? Are some tets inside-out? |
|
i think everything is ok with the mesh. checkMesh does not give any wornings |
|
checkMesh.out (6,165 bytes)
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.3.0-f5222ca19ce6 Exec : checkMesh Date : Feb 23 2015 Time : 13:38:50 Host : "aiolos" PID : 22835 Case : /gpfs/home/ENERCON1/00059128/OpenFOAM/OpenFOAM-2.3.0/Projects/Askervein_Hill/Askervein_Hill_kepsilon/case210 nProcs : 1 sigFpe : Floating point exception trapping - not supported on this platform fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 2177778 faces: 5868703 internal faces: 5588894 cells: 1850745 faces per cell: 6.1908 boundary patches: 39 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 1664473 prisms: 6921 wedges: 0 pyramids: 0 tet wedges: 354 tetrahedra: 4 polyhedra: 178993 Breakdown of polyhedra by number of faces: faces number of cells 4 10622 5 6697 6 11931 7 40012 8 17040 9 83125 10 919 11 216 12 5301 13 3 14 1 15 2955 18 171 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology p1 276 332 ok (non-closed singly connected) p2 276 332 ok (non-closed singly connected) p3 276 332 ok (non-closed singly connected) p4 276 332 ok (non-closed singly connected) p5 276 332 ok (non-closed singly connected) p6 276 332 ok (non-closed singly connected) p7 276 332 ok (non-closed singly connected) p8 276 332 ok (non-closed singly connected) p9 276 332 ok (non-closed singly connected) p10 276 332 ok (non-closed singly connected) p11 276 332 ok (non-closed singly connected) p12 276 332 ok (non-closed singly connected) p13 276 332 ok (non-closed singly connected) p14 276 332 ok (non-closed singly connected) p15 276 332 ok (non-closed singly connected) p16 276 332 ok (non-closed singly connected) p17 276 332 ok (non-closed singly connected) p18 276 332 ok (non-closed singly connected) p19 276 332 ok (non-closed singly connected) p20 276 332 ok (non-closed singly connected) p21 276 332 ok (non-closed singly connected) p22 276 332 ok (non-closed singly connected) p23 276 332 ok (non-closed singly connected) p24 276 332 ok (non-closed singly connected) p25 276 332 ok (non-closed singly connected) p26 276 332 ok (non-closed singly connected) p27 276 332 ok (non-closed singly connected) p28 276 332 ok (non-closed singly connected) p29 276 332 ok (non-closed singly connected) p30 276 332 ok (non-closed singly connected) p31 276 332 ok (non-closed singly connected) p32 276 332 ok (non-closed singly connected) p33 276 332 ok (non-closed singly connected) p34 276 332 ok (non-closed singly connected) p35 276 332 ok (non-closed singly connected) p36 276 332 ok (non-closed singly connected) BOTTOM 0 0 ok (empty) TOP 261 280 ok (non-closed singly connected) Terrain 269612 281125 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (69522.7 17522.7 -1.39002e-07) (81477.2 29477.3 2000) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (-5.10488e-17 -9.79937e-17 8.78762e-15) OK. Max cell openness = 4.04411e-16 OK. Max aspect ratio = 81.1674 OK. Minimum face area = 1.31317. Maximum face area = 600167. Face area magnitudes OK. Min volume = 7.18999. Max volume = 1.1964e+08. Total volume = 2.23475e+11. Cell volumes OK. Mesh non-orthogonality Max: 64.9998 average: 17.1163 Non-orthogonality check OK. Face pyramids OK. ***Max skewness = 4.51365, 2 highly skew faces detected which may impair the quality of the results <<Writing 2 skew faces to set skewFaces Coupled point location match (average 0) OK. Failed 1 mesh checks. End |
|
I uploaded the output of the checkMesh utility |
|
Can you check with checkMesh -allGeometry -allTopology |
|
ok now the output is different. here the final output of checkmesh: Checking geometry... Overall domain bounding box (69522.7 17522.7 -1.39002e-07) (81477.2 29477.3 2000) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (-5.10488e-17 -9.79937e-17 8.78762e-15) OK. Max cell openness = 4.04411e-16 OK. Max aspect ratio = 81.1674 OK. Minimum face area = 1.31317. Maximum face area = 600167. Face area magnitudes OK. Min volume = 7.18999. Max volume = 1.1964e+08. Total volume = 2.23475e+11. Cell volumes OK. Mesh non-orthogonality Max: 64.9998 average: 17.1163 Non-orthogonality check OK. Face pyramids OK. ***Max skewness = 4.51365, 2 highly skew faces detected which may impair the quality of the results <<Writing 2 skew faces to set skewFaces Coupled point location match (average 0) OK. Face tets OK. Min/max edge length = 0.235511 1045.87 OK. *There are 776 faces with concave angles between consecutive edges. Max concave angle = 79.9294 degrees. <<Writing 776 faces with concave angles to set concaveFaces Face flatness (1 = flat, 0 = butterfly) : min = 0.10157 average = 0.999385 *There are 1394 faces with ratio between projected and actual area < 0.8 Minimum ratio (minimum flatness, maximum warpage) = 0.10157 <<Writing 1394 warped faces to set warpedFaces Cell determinant (wellposedness) : minimum: 2.45496e-05 average: 2.55151 ***Cells with small determinant (< 0.001) found, number of cells: 148 <<Writing 148 under-determined cells to set underdeterminedCells ***Concave cells (using face planes) found, number of cells: 67108 <<Writing 67108 concave cells to set concaveCells Face interpolation weight : minimum: 0.032057 average: 0.469562 ***Faces with small interpolation weight (< 0.05) found, number of faces: 31 <<Writing 31 faces with low interpolation weights to set lowWeightFaces Face volume ratio : minimum: 0.0111648 average: 0.869301 Face volume ratio check OK. are concave cells a problem? |
|
> are concave cells a problem? Yes if they cause negative face-pyramids, but this looks OK for your mesh. Could you sub-set a section of the mesh containing cells in which points are not found and those point locations and up-load it for testing? |
|
I subset the mesh with all not found points. unfortunately the security adds of my browser does not allow me to upload the file. can I sent you an email? |
|
Could you send the case to some distribution repository so that others may download it and test possible patches to the search algorithm? |
|
boundary (5,445 bytes)
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format binary; class polyBoundaryMesh; location "1/polyMesh"; object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 40 ( p1 { type patch; nFaces 0; startFace 350829; } p2 { type patch; nFaces 0; startFace 350829; } p3 { type patch; nFaces 0; startFace 350829; } p4 { type patch; nFaces 0; startFace 350829; } p5 { type patch; nFaces 0; startFace 350829; } p6 { type patch; nFaces 0; startFace 350829; } p7 { type patch; nFaces 0; startFace 350829; } p8 { type patch; nFaces 0; startFace 350829; } p9 { type patch; nFaces 0; startFace 350829; } p10 { type patch; nFaces 0; startFace 350829; } p11 { type patch; nFaces 0; startFace 350829; } p12 { type patch; nFaces 0; startFace 350829; } p13 { type patch; nFaces 0; startFace 350829; } p14 { type patch; nFaces 0; startFace 350829; } p15 { type patch; nFaces 0; startFace 350829; } p16 { type patch; nFaces 0; startFace 350829; } p17 { type patch; nFaces 0; startFace 350829; } p18 { type patch; nFaces 0; startFace 350829; } p19 { type patch; nFaces 0; startFace 350829; } p20 { type patch; nFaces 0; startFace 350829; } p21 { type patch; nFaces 0; startFace 350829; } p22 { type patch; nFaces 0; startFace 350829; } p23 { type patch; nFaces 0; startFace 350829; } p24 { type patch; nFaces 0; startFace 350829; } p25 { type patch; nFaces 0; startFace 350829; } p26 { type patch; nFaces 0; startFace 350829; } p27 { type patch; nFaces 0; startFace 350829; } p28 { type patch; nFaces 0; startFace 350829; } p29 { type patch; nFaces 0; startFace 350829; } p30 { type patch; nFaces 0; startFace 350829; } p31 { type patch; nFaces 0; startFace 350829; } p32 { type patch; nFaces 0; startFace 350829; } p33 { type patch; nFaces 0; startFace 350829; } p34 { type patch; nFaces 0; startFace 350829; } p35 { type patch; nFaces 0; startFace 350829; } p36 { type patch; nFaces 0; startFace 350829; } BOTTOM { type wall; inGroups 1(wall); nFaces 0; startFace 350829; } TOP { type wall; inGroups 1(wall); nFaces 0; startFace 350829; } Terrain { type wall; inGroups 2 ( motorBikeGroup wall ) ; nFaces 18933; startFace 350829; } oldInternalFaces { type empty; inGroups 1(empty); nFaces 4841; startFace 369762; } ) // ************************************************************************* // |
|
|
|
boundary1 (5,404 bytes)
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format binary; class polyBoundaryMesh; location "1/polyMesh"; object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 40 ( p1 { type patch; nFaces 0; startFace 73037; } p2 { type patch; nFaces 0; startFace 73037; } p3 { type patch; nFaces 0; startFace 73037; } p4 { type patch; nFaces 0; startFace 73037; } p5 { type patch; nFaces 0; startFace 73037; } p6 { type patch; nFaces 0; startFace 73037; } p7 { type patch; nFaces 0; startFace 73037; } p8 { type patch; nFaces 0; startFace 73037; } p9 { type patch; nFaces 0; startFace 73037; } p10 { type patch; nFaces 0; startFace 73037; } p11 { type patch; nFaces 0; startFace 73037; } p12 { type patch; nFaces 0; startFace 73037; } p13 { type patch; nFaces 0; startFace 73037; } p14 { type patch; nFaces 0; startFace 73037; } p15 { type patch; nFaces 0; startFace 73037; } p16 { type patch; nFaces 0; startFace 73037; } p17 { type patch; nFaces 0; startFace 73037; } p18 { type patch; nFaces 0; startFace 73037; } p19 { type patch; nFaces 0; startFace 73037; } p20 { type patch; nFaces 0; startFace 73037; } p21 { type patch; nFaces 0; startFace 73037; } p22 { type patch; nFaces 0; startFace 73037; } p23 { type patch; nFaces 0; startFace 73037; } p24 { type patch; nFaces 0; startFace 73037; } p25 { type patch; nFaces 0; startFace 73037; } p26 { type patch; nFaces 0; startFace 73037; } p27 { type patch; nFaces 0; startFace 73037; } p28 { type patch; nFaces 0; startFace 73037; } p29 { type patch; nFaces 0; startFace 73037; } p30 { type patch; nFaces 0; startFace 73037; } p31 { type patch; nFaces 0; startFace 73037; } p32 { type patch; nFaces 0; startFace 73037; } p33 { type patch; nFaces 0; startFace 73037; } p34 { type patch; nFaces 0; startFace 73037; } p35 { type patch; nFaces 0; startFace 73037; } p36 { type patch; nFaces 0; startFace 73037; } BOTTOM { type wall; inGroups 1(wall); nFaces 0; startFace 73037; } TOP { type wall; inGroups 1(wall); nFaces 0; startFace 73037; } Terrain { type wall; inGroups 2 ( motorBikeGroup wall ) ; nFaces 4738; startFace 73037; } oldInternalFaces { type empty; inGroups 1(empty); nFaces 2408; startFace 77775; } ) // ************************************************************************* // |
|
cellZones1 (96,987 bytes)
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format binary; class regIOobject; location "1/polyMesh"; object cellZones; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 1 ( Box1 { type cellZone; cellLabels List<label> 24012 ( ! " # $ % & |