[kiran]

Hi FOAMers

i am solving solving backward facing step problem with rhoTurbFoam( transient solver for compressible turbulent flows).

Inlet conditions are

P(static) = 39KPa

T (total)= 300k

Mach no= 2

and i am using K-Epsilon turbulence model.

for the above stated conditions

solver is not running till the given time, it is giving the error as follows


Courant Number mean: 0.12239 max: 0.673613
Time = 0.00153982

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.00028022, Final residual = 7.95395e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.000296888, Final residual = 4.06123e-06, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.000420226, Final residual = 6.53906e-06, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.000579806, Final residual = 2.0317e-07, No Iterations 1
DICPCG: Solving for p, Initial residual = 2.03148e-07, Final residual = 2.03148e-07, No Iterations 0
DICPCG: Solving for p, Initial residual = 2.03148e-07, Final residual = 2.03148e-07, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.27198e-07, global = -4.81332e-08, cumulative = -7.72037e-05
DILUPBiCG: Solving for h, Initial residual = 0.00744003, Final residual = 2.16506e-07, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.00145661, Final residual = 9.72992e-09, No Iterations 2
DICPCG: Solving for p, Initial residual = 9.98167e-09, Final residual = 9.98167e-09, No Iterations 0
DICPCG: Solving for p, Initial residual = 9.98167e-09, Final residual = 9.98167e-09, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.08022e-09, global = 2.2917e-09, cumulative = -7.72014e-05
DILUPBiCG: Solving for epsilon, Initial residual = 0.167629, Final residual = 6.67002e-07, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.161715, Final residual = 3.20616e-06, No Iterations 1
bounding k, min: -125860 max: 58146.2 average: 717.796
ExecutionTime = 207.04 s ClockTime = 239 s

Courant Number mean: 0.122003 max: 0.764888
Time = 0.00154009

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.0134332, Final residual = 4.50173e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0112129, Final residual = 7.02826e-06, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.0107112, Final residual = 2.1616e-06, No Iterations 2


Maximum number of iterations exceeded#0 Foam::error:rintStack(Foam::Ostream&) in "/home/user/OpenFOAM/1.5/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/user/OpenFOAM/1.5/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::hThermo<Foam:ureMixture<Foam::constTranspo rt<Foam::specieThermo<Foam::hConstThermo<Foam:er fectGas> > > > >::calculate() in "/home/user/OpenFOAM/1.5/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#3 Foam::hThermo<Foam:ureMixture<Foam::constTranspo rt<Foam::specieThermo<Foam::hConstThermo<Foam:er fectGas> > > > >::correct() in "/home/user/OpenFOAM/1.5/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#4 main in "/home/user/OpenFOAM/1.5/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/rhoTurbFoam"
#5 __libc_start_main in "/lib/libc.so.6"
#6 Foam::regIOobject::readIfModified() in "/home/user/OpenFOAM/1.5/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/rhoTurbFoam"


From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 78.

FOAM aborting

this is the error

i have given outlet conditions as zerogradients

so please help me out!!!!!!!!!!!!!!!!!.