View Issue Details

IDProjectCategoryView StatusLast Update
0001157OpenFOAMBugpublic2015-06-16 15:34
Reporteruser612Assigned Tohenry  
PrioritynormalSeverityminorReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSUbuntuOS Version12.04
Summary0001157: faceAgglomerate crashes with hanging pointer error
DescriptionA geometry is meshed with snappyHexMesh and then the patches are agglomerated using the faceAgglomeration utility.

faceAgglomeration crashes with the following error:

Agglomerating patch : geometry_ascii


--> FOAM FATAL ERROR:
hanging pointer at index 4 (size 50), cannot dereference

    From function PtrList::operator[]
    in file /home/user/OpenFOAM/OpenFOAM-2.2.x/src/OpenFOAM/lnInclude/PtrListI.H at line 150.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/home/user/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/user/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::pairPatchAgglomeration::setEdgeWeights(int) in "/home/user/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libpairPatchAgglomeration.so"
#3 Foam::pairPatchAgglomeration::agglomerate() in "/home/user/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libpairPatchAgglomeration.so"
#4
 in "/home/user/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/faceAgglomerate"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6
 in "/home/user/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/faceAgglomerate"
Aborted (core dumped)
Steps To ReproduceRun the attached case with the latest git version of OpenFOAM 2.2.x
TagsNo tags attached.

Activities

user612

2014-02-12 11:06

 

user21

2014-02-17 09:30

  ~0002820

The case you sent is extremely large. We would need a small case which can
reproduce the problem for the faceAgglomeration.

Regards

user612

2014-02-17 10:33

  ~0002821

I am sorry, that I can can not give you a smaller case, but this problem only occurs with rather complex geometries. Simple geometries like boxes run smoothly.

user854

2014-02-19 22:03

  ~0002861

In viewFactorsDict, the geometry patch is called "geometry_ascii", whereas it is called "geometry" in snappyHexMeshDict. Changing back to "geometry" in viewFactorsDict does not make the faceAgglomeration run, however it changes the error that is generated.

user612

2014-02-19 22:09

  ~0002862

It depends on, which version of OpenFOAM you are using. As this bug is reported for version 2.2.x the resulting patch name will be "geometry_ascii". If you are using version 2.3.x then the name will be "geometry".

user854

2014-02-20 15:08

  ~0002870

You're right indeed, Martin.
Well, maybe the following could help better, then. The patch geometry_ascii has 10067 faces. I tried to use faceAgglomerate with different values for nFacesInCoarsestLevel, keeping all other settings unchanged.

10000, pass
5000, pass
2500, pass
2000, pass
1500, pass
1400, pass
1350, pass
1325, pass
1320, pass
1319, pass
1318, pass
1317, fail
1316, fail
1315, fail
1310, fail
1300, fail
1250, fail
1000, fail
 500, fail

derekm

2015-04-28 12:19

reporter   ~0004669

This bug still exists in 2.3.1 with meshes created by blockmesh.

wyldckat

2015-05-28 13:03

updater   ~0004825

@derekm: Can you provide a test case with such a mesh?

derekm

2015-05-29 12:26

reporter  

faceglomtest.zip (66,107 bytes)

derekm

2015-05-29 12:31

reporter   ~0004837

Last edited: 2015-05-29 12:35

see attached file faceglomtest.zip two cases one fails the other passes only difference is the file faceglomtest/fail/constant/topAir/viewFactorsDict line 32
if set to 210 it fails set to 250 it passes in single processor
if you run it in parallel the threshold changes in value upwards.

The bug is in 2.4.0 as well O/S opensuse 12.3

henry

2015-06-13 23:15

manager   ~0004923

I have spent a bit of time testing faceAgglomerate on your case with valgrind and the issue relates to the logic in

Foam::pairPatchAgglomeration:: agglomerate()

in which the results of agglomeratePatch were being used even if the return flag indicates that they are not valid. I have corrected the logic and your case now runs fine but I haven't tested any other cases. I have pushed the correction into OpenFOAM-2.4.x and OpenFOAM-dev; could you test either and report back?

user612

2015-06-15 15:17

 

faceAgglomerationTest.tar.gz (1,501,355 bytes)

user612

2015-06-15 15:25

  ~0004928

I have added a new test case (faceAgglomerationTest). I tried faceAgglomeration on region FLUID using the new logic. Now it doesn't crash anymore, but it doesn't stop either. I have the impression, that faceAgglomeration is trapped in an endless loop for patch FLUID_to_plate.
Maybe it is not possible for the algorithm to connect more faces. In this case the algorithm should quit instead of endless trying.

henry

2015-06-15 16:14

manager   ~0004929

I have just pushed a further fix into OpenFOAM-2.4.x to stop the agglomeration loop if a reduction in the number of faces is not achieved. It is not clear why the agglomeration stalls but at least this fix avoids the infinite loop. Let me know if it works.

If you need further agglomeration than this provides you will need to analyse the algorithm and look at the intermediate results to find out why the agglomerations stalls.

user612

2015-06-16 09:41

  ~0004948

Last edited: 2015-06-16 10:03

Thank you henry, this fix works. Unfortunately for complex meshes the algorithm stops the agglomeration very early. This results in a very large number of agglomerations. For me viewFactor radiation is still not usable for real world problems.

Maybe I try to analyse the pairPatchAgglomeration code sometime. I already started to develop a new agglomeration logic, but I haven't finished it yet.

henry

2015-06-16 09:57

manager   ~0004949

We have not received any funding for the development or maintenance of the viewFactor radiation functionality since the initial release in 2011. If you or others can secure funding we will do more work on it alternatively if you can provide bug-fixes or improvements we can include them.

Issue History

Date Modified Username Field Change
2014-02-12 11:06 user612 New Issue
2014-02-12 11:06 user612 File Added: faceAgglomerationBug.tar.gz
2014-02-17 09:30 user21 Note Added: 0002820
2014-02-17 10:33 user612 Note Added: 0002821
2014-02-19 22:03 user854 Note Added: 0002861
2014-02-19 22:09 user612 Note Added: 0002862
2014-02-20 15:08 user854 Note Added: 0002870
2015-04-28 12:19 derekm Note Added: 0004669
2015-05-28 13:03 wyldckat Note Added: 0004825
2015-05-29 12:26 derekm File Added: faceglomtest.zip
2015-05-29 12:31 derekm Note Added: 0004837
2015-05-29 12:32 derekm Note Edited: 0004837
2015-05-29 12:35 derekm Note Edited: 0004837
2015-06-13 23:15 henry Note Added: 0004923
2015-06-15 15:17 user612 File Added: faceAgglomerationTest.tar.gz
2015-06-15 15:25 user612 Note Added: 0004928
2015-06-15 16:14 henry Note Added: 0004929
2015-06-16 09:41 user612 Note Added: 0004948
2015-06-16 09:57 henry Note Added: 0004949
2015-06-16 10:03 user612 Note Edited: 0004948
2015-06-16 15:33 henry Status new => resolved
2015-06-16 15:33 henry Resolution open => fixed
2015-06-16 15:33 henry Assigned To => henry