Changes in OpenFOAM8 concerning the buoyantSimpleFoam solver ?
 

Hi foamers,
i am currently trying to validate the OpenFOAM v8 with some test-cases and i am facing a problem using the buoyantSImpleFoam.
I simulate a forced convection pipe, which works without problems with OpenFOAM v6. But when i use OpenFOAM v8, the job crashes like that :




#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::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) at ??:?
#4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<Foam::DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) at ??:?
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10 Foam::fvMatrix<double>::solve() at ??:?
#11 ? at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? at ??:?
Exception en point flottant (core dumped)


Before this message, i can see that epsilon and k are very high, with for example :



bounding epsilon, min: -4.2649481e+68 max: 5.2595018e+72 average: 1.8901264e+69



Any ideas ?


Thanks a lot.


L.

Has anybody encountered the same issue with buoyantSimpleFoam ?
The problem happens again with OpenFOAM v9.

Please tell me about your experience if you have faced it during your OpenFOAM work.


Have a good day.


Laurent

Well since these versions do not have buoyantSimpleFoam but only buoyantFoam you are likely not using the SIMPLE algorithm correctly. So your fvSolution settings need to be changed to adapt to this. Post those here

Thank you for your help,


i post here my fvSolution to give you more nformations.
Please ask me if you need some additional informations to help me.
And thanks again.


Laurent