[R_21]

Hi Foamers,

I am trying to simulate a swirling flow for compressible flow conditions. I was able to do the same using SimpleFoam. But while using rhoCentralFoam , I am getting the following error.

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam:perator/(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#5 Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
#9
at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11

I have tried changing the values in 0 folder and the BC. Initially I thought it was due to Swirl , but when I used uniform flow conditions also the same error pops up.

Any Suggestions helpful .

Regards,

Rafath

[chriss85]

Try lower timestep or other scheme/matrix solver?

[R_21]

Thanks a lot for the help .
The problem was with the DDt definition that is had given in FVschemes.

Although that issue is fixed.
later on when I run the Simulation this pops :
"
--> FOAM FATAL ERROR:
incompatible dimensions for operation
[rhoE[0 2 -3 0 0 0 0] ] - [div(sigmaDotU)[1 -1 -3 0 0 0 0] ]
"

I tried changing the rhoE solver type in FVsolution. But none seems to work.

[chriss85]

Have you checked that your units are ok in each of the initial fields?

[R_21]

Thanks chriss85 for the help.
I didn't check my unit of pressure assuming it would not change in compressible flow.

One more thing. Is there a way to understand printstack errors.

Regards,
Rafath

[chriss85]

Only somewhat. It will tell you in which function it crashes and the call stack. Most solvers are written in a linear fashion in a single function so it's not always helping. If you have a reproducible problem you're best of inserting some console output throughout the code and tracing the line that crashes.

[mkraposhin]

Originally Posted by R_21

Thanks chriss85 for the help.
I didn't check my unit of pressure assuming it would not change in compressible flow.

Regards,
Rafath

In OpenFOAM compressible and incompressible pressures have different units

Also, try to set viscousity to zero, this can help you to find place of error, because this will switch model to Euler equations