[alimea]

Hi all

I have combined viscoelastic solver with buoyantBuossinesqPimpleFoam solver to solve free convection flow of viscoelastic fluid in a cavity. when I run the program after some time steps the run is stopped and I get this error:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5  Foam::fvMatrix<Foam::SymmTensor<double> >::solveSegregated(Foam::dictionary const&) at ??:?
#6  Foam::fvMatrix<Foam::SymmTensor<double> >::solve(Foam::dictionary const&) at ??:?
#7  Foam::Oldroyd_B::correct() at ??:?
#8  
 at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  
 at ??:?
./Allrun: line 15: 14497 Floating point exception(core dumped) buoyantBoussinesqPimpleViscoFoam > log

meanwhile there is no problem when I run such a problem for Newtonian fluid.
when I wanted to find the reason, I saw high values of stress (in order of 10^120) in results.
Could you plz help me?

Thanks

[Bloerb]

This is a crash from the solver. Without further information this could be due to bad boundary conditions mesh, initilization, schemes, etc and of course your programming. Check if the case runs for both solvers seperatly and then try your new one with heat transfer disabled nad then enabled.

[alimea]

Originally Posted by Bloerb

This is a crash from the solver. Without further information this could be due to bad boundary conditions mesh, initilization, schemes, etc and of course your programming. Check if the case runs for both solvers seperatly and then try your new one with heat transfer disabled nad then enabled.

Hi

Thanks for your answer.
I tried many amounts of mesh, but it isn't solved yet. Also you said that I had to solve the problem in both solver seperately, but it's impossible. because free convection is possible just with solving energy eqn. so I can't solve free convection problem by viscoelastic solver without energy eqn.

When I check the obtained values of stress tensor, I find large value! What is the problem?