2017-11-24 09:12 GMT

View Issue Details Jump to Notes ]
IDProjectCategoryView StatusLast Update
0002745OpenFOAM[All Projects] Bugpublic2017-11-17 12:25
Reporterjohnk 
Assigned To 
PriorityurgentSeveritycrashReproducibilityalways
StatusnewResolutionopen 
PlatformlinuxOSubuntuOS Version16.04
Product Version5.x 
Target VersionFixed in Version 
Summary0002745: interDyMFoam does not work in parallel when using
DescriptionNow I run the case on OF 5.x

In my case I have rectangular box with BCs as follows:

front, back, left, right, up down

I initialize the velocity field (using swak4Foam) which corresponds to two
counter rotating wake vortices over water surface in the domain. Their axes are
perpendicular to front and back faces.
I need to have cyclic BC on front and back faces. The case works fine without
dynamicRefineFvMesh and interFoam. However, when I use dynamicRefineFvMesh (I
want the grid to be refined on the boarder between air and water), I set up
cyclicAMI on front and back patches. But the case works only with one processor.
When I try to run in after decomposition, it gives me an error:

AMI: Creating addressing and weights between 120990 source faces and 120990 target faces
--> FOAM Warning :
    From function void Foam::AMIMethod<SourcePatch, TargetPatch>::checkPatches() const [with SourcePatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double> >&>; TargetPatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double> >&>]
    in file lnInclude/AMIMethod.C at line 57
    Source and target patch bounding boxes are not similar
    source box span : (0 300 92)
    target box span : (0 80.91743 76.77211)
    source box : (0 0 0) (0 300 92)
    target box : (1.08 172.844 15.22789) (1.08 253.7615 92)
    inflated target box : (-4.497093 167.2669 9.650795) (6.657093 259.3386 97.57709)
--> FOAM Warning :
    From function bool Foam::AMIMethod<SourcePatch, TargetPatch>::initialise(Foam::labelListList&, Foam::scalarListList&, Foam::labelListList&, Foam::scalarListList&, Foam::label&, Foam::label&) [with SourcePatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double> >&>; TargetPatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double> >&>; Foam::labelListList = Foam::List<Foam::List<int> >; Foam::scalarListList = Foam::List<Foam::List<double> >; Foam::label = int]
    in file lnInclude/AMIMethod.C at line 97
    12078 source faces but no target faces
AMI: Patch source sum(weights) min/max/average = 0, 1, 0.1004546
AMI: Patch target sum(weights) min/max/average = 0, 1, 0.1004546
Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Selecting laminar stress model Stokes

Reading g

Reading hRef
Calculating field g.h

No MRF models present

No finite volume options present

DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 0.0978869, No Iterations 515
time step continuity errors : sum local = 8.905907e-06, global = 2.555122e-20, cumulative = 2.555122e-20
Reading/calculating face velocity Uf

Courant Number mean: 0.005911138 max: 0.7790846

Starting time loop

Courant Number mean: 0.003789191 max: 0.4994132
Interface Courant Number mean: 0 max: 0
deltaT = 0.001282051
Time = 0.001282051

PIMPLE: iteration 1
Selected 0 split points out of a possible 0.
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.1086957 Min(alpha.water) = 0 Max(alpha.water) = 1.000795
[0] #0 Foam::error::printStack(Foam::Ostream&)[3] #0 Foam::error::printStack(Foam::Ostream&)[4] #0 Foam::error::printStack(Foam::Ostream&)[7] #0 Foam::error::printStack(Foam::Ostream&)[9] #0 Foam::error::printStack(Foam::Ostream&)[2] #0 Foam::error::printStack(Foam::Ostream&)[5] #0 Foam::error::printStack(Foam::Ostream&)[6] #0 Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
 at ??:?
[3] #1 Foam::sigFpe::sigHandler(int)[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
 at ??:?
 at ??:?
[4] #1 Foam::sigFpe::sigHandler(int)[9] #1 Foam::sigFpe::sigHandler(int)[5] #1 Foam::sigFpe::sigHandler(int) at ??:?
[7] #1 Foam::sigFpe::sigHandler(int) at ??:?
 at ??:?
[6] #1 Foam::sigFpe::sigHandler(int)[2] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ? at ??:?
[0] #2 ? at ??:?
[9] #2 ? at ??:?
[5] #2 ? at ??:?
[6] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&)[3] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[7] #2 ? at ??:?
[4] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
 in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #3 at ??:?
[1] #2 ?[6] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[2] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[0] #4 in "/lib/x86_64-linux-gnu/libc.so.6"
 at ??:?
[3] #4 [1] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[9] #4 in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[5] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[6] #4 at ??:?
[7] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[4] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[3] #5 at ??:?
[0] #5 at ??:?
[1] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[9] #5 at ??:?
[2] #4 ? at ??:?
[5] #5 ?Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[6] #5 ?Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[4] #5 at ??:?
[7] #5 at ??:?
[3] #6 __libc_start_main? at ??:?
[0] #6 __libc_start_main at ??:?
[1] #5 ? at ??:?
[9] #6 __libc_start_main?? at ??:?
[2] #5 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #7 at ??:?
[5] #6 __libc_start_main? at ??:?
[6] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #7 at ??:?
[4] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #7 ? at ??:?
[7] #6 __libc_start_main??? in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #7 in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #7 at ??:?
[1] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #7 at ??:?
[2] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

I have already tried to put preservePatches (front back); in decomposeparDict,
but the issue remains

Could you please help me with that?
I attach the case.
TagsNo tags attached.
Attached Files

-Relationships
+Relationships

-Notes

~0009049

henry (manager)

It is not clear how to run this case. I unpacked the tgz file which generated the directories in the current directory rather than creating the case directory so I had to create the case directory by hand and move everything into it. Now I have a case directory inter_2D_check_2 containing

0.orig constant system

but no scripts or instructions to run the case.

Please provide a case which can be run in OpenFOAM-dev and without the need to install third-party software.
+Notes

-Issue History
Date Modified Username Field Change
2017-10-31 21:31 johnk New Issue
2017-10-31 21:41 johnk File Added: inter_2D_check_2.tar.gz
2017-11-17 12:25 henry Note Added: 0009049
+Issue History