interDyMFoam fails in parallel
 

Hello all,

first of all the outputs from the debug-version:

Create time

Create mesh for time = 0

...

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model laminar
time step continuity errors : sum local = 0, global = 0, cumulative = 0
2 additional processes aborted (not shown)


[0]
[0]
[0] index -1 out of range 0 ... 32767#0 Foam::error::printStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-1.5.x/src/OSspecific/Unix/printStack.C:203
#1 Foam::error::abort() at ~/OpenFOAM/OpenFOAM-1.5.x/src/OpenFOAM/lnInclude/error.C:229
#2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) at ~/OpenFOAM/OpenFOAM-1.5.x/src/OpenFOAM/lnInclude/errorManip.H:87
#3 Foam::UList<double>::checkIndex(int) const at ~/OpenFOAM/OpenFOAM-1.5.x/src/OpenFOAM/lnInclude/UListI.H:106
#4 Foam::UList<double>::operator[](int) at ~/OpenFOAM/OpenFOAM-1.5.x/src/OpenFOAM/lnInclude/UListI.H:118
#5 Foam::fvMatrix<double>::setReference(int, double const&) at ~/OpenFOAM/OpenFOAM-1.5.x/src/finiteVolume/lnInclude/fvMatrix.C:487
#6 main at ~/OpenFOAM/OpenFOAM-1.5.x/applications/solvers/multiphase/interDyMFoam/correctPhi.H:31
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/spyrou2/OpenFOAM/OpenFOAM-1.5.x/applications/bin/linux64GccDPDebug/interDyMFoam"
[0]
[0]
[0] From function UList<T>::checkIndex(const label)
[0] in file /home/spyrou2/OpenFOAM/OpenFOAM-1.5.x/src/OpenFOAM/lnInclude/UListI.H at line 102.
[0]
FOAM parallel run aborting


Now some details:

- the calculation fails at the point when the pressure-equation is to be solved

- boundary conditions for pressure:
all boundaries have a zeroGradient condition for the pressure
so a reference cell and reference value for pd are defined

-the calculation with interDyMFaom fails if:
- the above boundary condition for the pressure is applied and the computation is parallel

- in the following cases the calculation runs without problems:
-serial calculation with interDyMFoam
-parallel calculation (interDyMFoam) but with a fixedValue b.c. for the pressure (splitting the domain in such a way that every processor domain contains a part of that patch)
-when using interFoam (interFoam runs in parallel even with zeroGradient on all boundaries)


From the output of the debug version i see that there is an index out of bounds.

Does anybody know or have an idea why this problem occurs?

Thank You

Looks like a bug. Can you set up a testcase and post it on the OpenFOAM-bugs forum?

Thanks, Mattijs

Hi Mattijs,

on monday there will be an uploaded testcase.

Thank you