[Mhmnabil]

Hi everyone,

I am trying to solve an internal (square channel) compressible flow using bouyantPimpleFOAM. At the initial time steps, everything seems to converge fast, but the pressure is blowing up after some time steps. I have the screen printout of the solution:

Beginnig:
...
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 5.15734e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.19565e-07, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 5.07416e-09, No Iterations 4
DILUPBiCG: Solving for h, Initial residual = 0.999967, Final residual = 5.12982e-07, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.988257, Final residual = 0.00963005, No Iterations 13
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.59437e-07, global = -2.56886e-07, cumulative = -2.56886e-07
GAMGPCG: Solving for p_rgh, Initial residual = 0.0382374, Final residual = 5.06303e-08, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.7811e-13, global = -8.18409e-14, cumulative = -2.56887e-07
DILUPBiCG: Solving for nuTilda, Initial residual = 1, Final residual = 6.83767e-08, No Iterations 4
ExecutionTime = 10.36 s ClockTime = 11 s
...

When the simullation crashes:
...
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.419311, Final residual = 5.53834e-07, No Iterations 273
DILUPBiCG: Solving for Uy, Initial residual = 0.373357, Final residual = 9.16423e-07, No Iterations 265
DILUPBiCG: Solving for Uz, Initial residual = 0.241971, Final residual = 9.92968e-07, No Iterations 207
DILUPBiCG: Solving for h, Initial residual = 0.117909, Final residual = 9.75382e-07, No Iterations 170
[3] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
...

I have attached my boundary conditions files as well as fvSolution and fvSchemes. Also, I have attached a screenshot of pressure field. As you can see the pressure field is not continuous which is weird!!! Please take a look at them and let me know if anything is wrong.

Thanks a lot.