[Rapha]

Hi there,

I am using the rhoPorousMRFPimpleFoam solver for my base case however I am having problems with the output, getting increases in velocity and pressure when I should be getting the opposite. I have included the temperature equation into the solver.
The details are:
- Mesh imported from Salome with 5 meshes joined together in a compound mesh.
- Each boundary is specified, eg inlet, wall, regionCouple, etc
- Two porous regions either side of the inlet channel
- Inlet velocity of 30 m/s. Inlet temp of 700K
- the model should run by air flowing through the inlet at 700K, 30m/s, passing through and heating up the porous zones on either side of the inlet duct leaving the outlets on either side at a lower temperature.

The problems are several:
- The time scale I can reach is of the order 1e-4 s when really it should be of the order 1e3 s.
- The velocity increases rapidly to 267m/s at certain regions - why?
- The pressure increases from 1 bar to 1.5 bar
- The flow doesn't go from the inlet to the outlet through the porous zones.

If anyone could give any hints as to whether it's the fvschemes, solutions or solver I need to look at, I would be really grateful.

Regards,

Rapha

[Rapha]

This is the error code which I get for the last iteration:

Courant Number mean: 0.0166694 max: 0.605291
deltaT = 1.8885e-07
Time = 0.000333925

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.0115056, Final residual = 4.09179e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0112349, Final residual = 3.55671e-06, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.278506, Final residual = 1.42709e-07, No Iterations 2
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::calculate() in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#4 Foam::hPsiThermo<Foam:ureMixture<Foam::sutherlan dTransport<Foam::specieThermo<Foam::hConstThermo<F oam:erfectGas> > > > >::correct() in "/opt/openfoam200/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#5
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/HTTESS_Solver_P"
#6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7
in "/opt/openfoam200/platforms/linuxGccDPOpt/bin/HTTESS_Solver_P"
Floating point exception
ra@ubuntu:~/OpenFOAM/ra-2.0.0/run/3_08/HTTESS$