[swahono]

Hi All,

I'm simulating a conjugate heat transfer across a flat plate in a rectangular channel.
I modified the multiRegionHeater tutorial in OF 1.7.0, and used k-epsilon turbulence model.

The Case geometry setup is shown in the attached PNG file.
Basically, the top and bottom air is turbulent (Re ~ 10^6).
The top air is hot at 746K, while the bottom air is cold at 288K.
The two flow regions are separated by a solid plate.

I ran the case with realizableKE model, and the simulation blew up after 38 iterations.
Here is the last few iterations before it blowing up:

Time = 36

Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.3204164, Final residual = 0.0049385518, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.43623786, Final residual = 0.01306048, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.46554744, Final residual = 0.015560302, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.14392941, Final residual = 0.010392872, No Iterations 1
Min/max T:119.30898 433.59925
GAMG:  Solving for p_rgh, Initial residual = 0.049386189, Final residual = 0.00042936621, No Iterations 4
time step continuity errors : sum local = 68.178516, global = -12.111985, cumulative = -12.049972
Min/max rho:0.2 2
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0060067685, Final residual = 0.000382735, No Iterations 1
bounding epsilon, min: -810751.63 max: 1.6587192e+11 average: 2266471.6
DILUPBiCG:  Solving for k, Initial residual = 0.60163212, Final residual = 0.020304535, No Iterations 1
bounding k, min: -196.25661 max: 1.3913795e+12 average: 28036037
maxResidual: 0.46554744  convergence criterion: 0

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.010704973, Final residual = 0.000326993, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0189083, Final residual = 0.00060774355, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.018775159, Final residual = 0.00059849536, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.0093179606, Final residual = 9.2234893e-05, No Iterations 1
Min/max T:290.60892 746.00013
GAMG:  Solving for p_rgh, Initial residual = 0.0037615237, Final residual = 1.1872936e-05, No Iterations 3
time step continuity errors : sum local = 4.7313027e-07, global = -2.6202936e-08, cumulative = -12.049972
Min/max rho:0.47318397 1.2147625
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00086839005, Final residual = 1.9941133e-05, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.010316428, Final residual = 0.00084437086, No Iterations 1
maxResidual: 0.0189083  convergence criterion: 0

Solving for solid region solidWall
DICPCG:  Solving for T, Initial residual = 0.025633671, Final residual = 0.0021160135, No Iterations 1
Min/max T:287.64755 296.48869
maxResidual: 0.025633671  convergence criterion: 0
ExecutionTime = 200.4 s  ClockTime = 201 s

Time = 37


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.50298705, Final residual = 0.012268166, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.46978779, Final residual = 0.013510551, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.45493707, Final residual = 0.011565384, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.025038369, Final residual = 0.0019166388, No Iterations 1
Min/max T:-1557.9293 422.26051
GAMG:  Solving for p_rgh, Initial residual = 0.19809226, Final residual = 0.0019424248, No Iterations 4
time step continuity errors : sum local = 2.9035403e+10, global = 5.090356e+09, cumulative = 5.0903559e+09
Min/max rho:0.2 2
DILUPBiCG:  Solving for epsilon, Initial residual = 0.021184386, Final residual = 5.8090526e-05, No Iterations 3
bounding epsilon, min: -3578671.8 max: 5.937534e+11 average: 10175646
DILUPBiCG:  Solving for k, Initial residual = 0.99999971, Final residual = 0.050065216, No Iterations 2
maxResidual: 0.50298705  convergence criterion: 0

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.010419636, Final residual = 0.00031881238, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.018670427, Final residual = 0.00059330634, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.018306335, Final residual = 0.0005741435, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.0090171776, Final residual = 8.9433187e-05, No Iterations 1
Min/max T:290.65192 746.00013
GAMG:  Solving for p_rgh, Initial residual = 0.003663336, Final residual = 1.8532424e-05, No Iterations 3
time step continuity errors : sum local = 7.4935637e-07, global = -2.1118708e-08, cumulative = 5.0903559e+09
Min/max rho:0.47318398 1.214584
DILUPBiCG:  Solving for epsilon, Initial residual = 0.00083188985, Final residual = 1.8973471e-05, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.0098047206, Final residual = 0.0008637932, No Iterations 1
maxResidual: 0.018670427  convergence criterion: 0

Solving for solid region solidWall
DICPCG:  Solving for T, Initial residual = 0.025238773, Final residual = 0.0020939892, No Iterations 1
Min/max T:287.7893 295.63824
maxResidual: 0.025238773  convergence criterion: 0
ExecutionTime = 204.98 s  ClockTime = 206 s

Time = 38


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.3018213, Final residual = 0.004181555, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.39316677, Final residual = 0.0066244412, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.4739378, Final residual = 0.014966079, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.01415212, Final residual = 0.0010154766, No Iterations 1
[7] 
[7] 
[7] --> FOAM FATAL ERROR: 
[7] Maximum number of iterations exceeded
[7] 
[7]     From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
[7]     in file /home/stefano/OpenFOAM/OpenFOAM-1.7.0/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67.
[7] 
FOAM parallel run aborting
[7] 
[7] #0  Foam::error::printStack(Foam::Ostream&) in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libOpenFOAM.so"
[7] #1  Foam::error::abort() in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libOpenFOAM.so"
[7] #2  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
[7] #3  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
[7] #4  main in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/applications/bin/linux64GccDPOpt/chtMultiRegionSimpleFoam"
[7] #5  __libc_start_main in "/lib64/libc.so.6"
[7] #6  Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/stefano/OpenFOAM/OpenFOAM-1.7.0/applications/bin/linux64GccDPOpt/chtMultiRegionSimpleFoam"
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 7 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun has exited due to process rank 7 with PID 5379 on
node barracuda exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------

Please if anyone can help me....It is urgent.
I have attached the /system and /constant files.

The case ran well with KwSST model.

Thank you in advance.

Best Regards,
Stefano

run6-topAir-bottomAir-Streamlines2Ts.jpg

parallelPlates.tar.gz