[sydonium]

Hi all,

I am trying to run an intensive simulation with supersonic velocities and high temperature. Whenever I attempt to increase the temperature or speed, I get an error. In the case I have below, the program ran but died after a few steps.

In the past, I have tried to no avail to solve this problem by changing the value in:
~/thermophysicalModels/specie/lnInclude/thermoI.H at line 76

Are there any ways to edit the thermophysical properties to allow for more heat? Or is this a problem of disk space? Thanks!

Time = 3.5e-05

Courant Number mean: 0.000722438 max: 0.632063
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.0839927, Final residual = 4.41532e-06, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.0510721, Final residual = 3.37884e-06, No Iterations 4
DILUPBiCG: Solving for e, Initial residual = 0.038889, Final residual = 8.06531e-07, No Iterations 3
DILUPBiCG: Solving for p, Initial residual = 0.0416342, Final residual = 3.71647e-11, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.1378e-12, global = 9.65176e-14, cumulative = -1.36379e-11
PIMPLE: iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 0.045368, Final residual = 6.08536e-06, No Iterations 11
DILUPBiCG: Solving for Uy, Initial residual = 0.024744, Final residual = 5.53074e-06, No Iterations 9
DILUPBiCG: Solving for e, Initial residual = 0.0354167, Final residual = 7.05416e-06, No Iterations 11
DILUPBiCG: Solving for p, Initial residual = 0.026016, Final residual = 8.87562e-09, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.90647e-10, global = 8.57064e-11, cumulative = 7.20684e-11
DILUPBiCG: Solving for epsilon, Initial residual = 0.276884, Final residual = 6.09947e-09, No Iterations 20
DILUPBiCG: Solving for k, Initial residual = 0.232956, Final residual = 4.11722e-09, No Iterations 22
ExecutionTime = 65 s ClockTime = 65 s

Time = 4e-05

Courant Number mean: 0.000754311 max: 5.18107
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.0767201, Final residual = 6.22788e-06, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.0416175, Final residual = 3.18463e-06, No Iterations 4
DILUPBiCG: Solving for e, Initial residual = 0.0311735, Final residual = 7.26361e-07, No Iterations 3


--> FOAM FATAL ERROR:
Maximum number of iterations exceeded

From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /home/openfoam/OpenFOAM/OpenFOAM-2.3.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::TEs(double, double, double) const at ??:?
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
^CAborted (core dumped)

[ssss]

Your courant number is blowing up. Reduce timestep and modify controlDict to allow the timestep to adapt to the courant number