View Issue Details

Jump to NotesJump to History
IDProjectCategoryView StatusDate SubmittedLast Update
0003978OpenFOAMBugpublic2023-05-05 14:562023-05-18 09:04
Reporterhandrake0724 Assigned Towill  
PrioritynormalSeverityminorReproducibilityhave not tried
Status closedResolutionunable to reproduce 
Product Version10 
Summary0003978: createNonConformalCouples fails in propeller tutorial
DescriptionI am testing NCC features in OpenFOAM 10 with propeller tutorial in incompressible/pimpleFoam/RAS/propeller.

for the rotating zone to include fully propeller, I just replace the propeller.obj with a simple cylinder just enough to be inside of the inner cylinder.
After running Allmesh, createNonConformalCouples produced the following errors. I did it without propeller geometry but got the same results.
Since the surfaces which will be NCC interfaces is just simple geometry, I think it would not make any problems.
I am familiar with the NCC feature and not sure what the cause of the errors is from.
FYI, I uploaded the tutorial case I tested.


Pstream initialised with:
    floatTransfer : 0
    nProcsSimpleSum : 0
    commsType : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Adding nonConformalCyclic interfaces between patches:
    (nonCouple1 nonCouple2)
Named:
    (nonConformalCyclic_on_nonCouple1 nonConformalCyclic_on_nonCouple2)
With transform:
    transformType none;

patchToPatch: Calculating couplings between 22632 source faces and 22632 target faces
[3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ?[0] #0 Foam::error::printStack(Foam::Ostream&) in "/usr/lib/libc.so.6"
[3] #3 Foam::Vector<double> Foam::face::centre<Foam::FixedList<Foam::Vector<double>, 4u> >(Foam::FixedList<Foam::Vector<double>, 4u> const&)[1] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #4 Foam::patchToPatches::intersection::intersectFaces(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, int, int) at ??:?
[1] #1 at ??:?
Foam::sigFpe::sigHandler(int)[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #5 Foam::patchToPatch::intersectPatchQueue(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, bool, Foam::DynamicList<Foam::Pair<int>, 0u, 2u, 1u> const&, Foam::List<int>&, Foam::DynamicList<Foam::Pair<int>, 0u, 2u, 1u>&, Foam::List<int> const&, Foam::List<bool>&, Foam::List<bool>&) at ??:?
[0] #2 ? at ??:?
[1] #2 ? at ??:?
[3] #6 Foam::patchToPatch::intersectPatches(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&) in "/usr/lib/libc.so.6"
[0] #3 Foam::Vector<double> Foam::face::centre<Foam::FixedList<Foam::Vector<double>, 4u> >(Foam::FixedList<Foam::Vector<double>, 4u> const&) in "/usr/lib/libc.so.6"
[1] #3 Foam::Vector<double> Foam::face::centre<Foam::FixedList<Foam::Vector<double>, 4u> >(Foam::FixedList<Foam::Vector<double>, 4u> const&) at ??:?
[3] #7 Foam::patchToPatch::update(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::transformer const&) at ??:?
[0] #4 Foam::patchToPatches::intersection::intersectFaces(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, int, int) at ??:?
[1] #4 Foam::patchToPatches::intersection::intersectFaces(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, int, int) at ??:?
[3] #8 Foam::patchToPatch::update(Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::transformer const&) at ??:?
[0] #5 Foam::patchToPatch::intersectPatchQueue(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, bool, Foam::DynamicList<Foam::Pair<int>, 0u, 2u, 1u> const&, Foam::List<int>&, Foam::DynamicList<Foam::Pair<int>, 0u, 2u, 1u>&, Foam::List<int> const&, Foam::List<bool>&, Foam::List<bool>&) at ??:?
[1] #5 Foam::patchToPatch::intersectPatchQueue(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, bool, Foam::DynamicList<Foam::Pair<int>, 0u, 2u, 1u> const&, Foam::List<int>&, Foam::DynamicList<Foam::Pair<int>, 0u, 2u, 1u>&, Foam::List<int> const&, Foam::List<bool>&, Foam::List<bool>&) at ??:?
[3] #9 Foam::nonConformalCyclicPolyPatch::intersection() const at ??:?
[0] #6 Foam::patchToPatch::intersectPatches(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&) at ??:?
[1] #6 Foam::patchToPatch::intersectPatches(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&) at ??:?
[3] #10 Foam::fvMeshStitcher::connect(bool, bool, bool) at ??:?
[0] #7 Foam::patchToPatch::update(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::transformer const&) at ??:?
[1] #7 Foam::patchToPatch::update(Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitiveOldTimePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::transformer const&) at ??:?
[0] #8 Foam::patchToPatch::update(Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::transformer const&) at ??:?
[1] #8 Foam::patchToPatch::update(Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::Field<Foam::Vector<double> > const&, Foam::PrimitivePatch<Foam::SubList<Foam::face>, Foam::Field<Foam::Vector<double> > const&> const&, Foam::transformer const&) at ??:?
[3] #11 at ??:?
[0] #9 Foam::nonConformalCyclicPolyPatch::intersection() const at ??:?
[1] #9 Foam::nonConformalCyclicPolyPatch::intersection() const? at ??:?
[3] # at ??:?
12 ? at ??:?
[0] #10 [1] #10 Foam::fvMeshStitcher::connect(bool, bool, bool)Foam::fvMeshStitcher::connect(bool, bool, bool) in "/usr/lib/libc.so.6"
[3] #13 __libc_start_main at ??:?
[0] #11 at ??:?
[1] #11 in "/usr/lib/libc.so.6"
[3] #14 ??? at ??:?
[0] #12 ? at ??:?
[1] #12 ? at ??:?
TagsNo tags attached.

Activities

handrake0724

2023-05-05 14:59

viewer   ~0013009

I forgot to attach the files.
propeller.tar.gz (39,919 bytes)

peksa

2023-05-06 18:58

reporter   ~0013010

I'm not able to reproduce the error. Case runs with by executing Allrun with OF-10.

Which OF-10 commit you're having?

handrake0724

2023-05-07 23:10

viewer   ~0013011

OF-10 commit is 73bdbb2a3c7b82bbe554ca304180b20b6a931591 dated on Apr 13. I guess it is the latest commit.
For just in case, I recompiled and tried it again but still have the same result.
I am not sure whether it is related or not, I tested it with gcc 11.1 and gcc 13.1 plus openmpi 4.1.1.

peksa

2023-05-08 06:46

reporter   ~0013012

Have you tried with OF-dev?

handrake0724

2023-05-09 01:14

viewer   ~0013013

yesterday, I tried with OF-dev and it worked without any problem.
I tried OF-10 in other machine. it worked without any problem.
So, I guess the problem is specific to my machine but don't have any clue.

will

2023-05-09 08:52

manager   ~0013014

I cannot reproduce this either. The stack trace does not narrow it down sufficiently for me to make any guesses.

What I'd do next: Try a full re-build. Try the mesh from a version/build that fails in a version/build that works. If that makes the previously working version/build fail, then upload that mesh somewhere so we can try it. Try building meshTools in debug to get a better trace (*). Start printing out values just before the failure.

I will close this in a few days if no progress has been made.

(*) I.e., add `-DFULLDEBUG -O0 -g` to EXE_INC in src/meshTools/Make/options

Issue History

Date Modified Username Field Change
2023-05-05 14:56 handrake0724 New Issue
2023-05-05 14:59 handrake0724 Note Added: 0013009
2023-05-05 14:59 handrake0724 File Added: propeller.tar.gz
2023-05-06 18:58 peksa Note Added: 0013010
2023-05-07 23:10 handrake0724 Note Added: 0013011
2023-05-08 06:46 peksa Note Added: 0013012
2023-05-09 01:14 handrake0724 Note Added: 0013013
2023-05-09 08:52 will Note Added: 0013014
2023-05-18 09:04 will Assigned To => will
2023-05-18 09:04 will Status new => closed
2023-05-18 09:04 will Resolution open => unable to reproduce