View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0002745 | OpenFOAM | Bug | public | 2017-10-31 21:31 | 2017-11-28 10:12 |
| Reporter | johnk | Assigned To | henry | ||
| Priority | urgent | Severity | crash | Reproducibility | always |
| Status | closed | Resolution | unable to reproduce | ||
| Platform | linux | OS | ubuntu | OS Version | 16.04 |
| Summary | 0002745: interDyMFoam does not work in parallel when using | ||||
| Description | Now 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. | ||||
| Tags | No tags attached. | ||||
|
|
|
|
|
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. |
|
|
Not enough information provided to reproduce the problem. |
| 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 | |
| 2017-11-28 10:12 | henry | Assigned To | => henry |
| 2017-11-28 10:12 | henry | Status | new => closed |
| 2017-11-28 10:12 | henry | Resolution | open => unable to reproduce |
| 2017-11-28 10:12 | henry | Note Added: 0009075 |
