CFD Online Discussion Forums - multiphaseInterFoam blow up at the third PIMPLE iteration.

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - multiphaseInterFoam blow up at the third PIMPLE iteration. (https://www.cfd-online.com/Forums/openfoam-solving/178124-multiphaseinterfoam-blow-up-third-pimple-iteration.html)

multiphaseInterFoam blow up at the third PIMPLE iteration.

Hello dear Foamers,

i'm using now the solver multiphaseInterFoam to solve a problem.(in OF40) But i always met the same problem no matter what i change. So i really need your helps.

Here is the error report:

Code:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?

#1  Foam::sigFpe::sigHandler(int) at ??:?

#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"

#3  Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

#4  Foam::operator*(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

#5  Foam::fixedValueFvPatchField<double>::gradientBoundaryCoeffs() const at ??:?

#6  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?

#7  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?

#8  ? at ??:?

#9  ? at ??:?

#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

#11  ? at ??:?

Gleitkomma-Ausnahme (Speicherabzug geschrieben)

It crashes at the third PIMPLE iteration. I'd like to paste my case-file here too, but it will be too much. So i'll just up load it. the file in triSurface is too big and i also leave it. Please have a look and any help is appreciated!

And it crashes both when i use parallel and serial.

Reading the error it's coming from a boundary condition, looking at your initial files all the boundary conditions for alphas are zero. I'm not an expert in multiphase flow to know what is the correct boundary condition, but if I remember right you cannot solve an ODE or PDE with all zeroGradient boundary conditions.

Quote:

Originally Posted by anishtain4 (Post 619752)

Hey Mahdi, thank you very much for your reply. I set it all zeroGradient because in the tutorial it is always zeroGrdient(and sometimes it doesn't exist in the 0 file). And one colleague of mine had already worked with this settings and he got no problem. I will try without alphas and see if it is the problem. Thank you very much for your hint!