|
[Sponsors] |
October 25, 2010, 03:34 |
Compressible epsilon blows up
|
#1 |
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 15 |
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: Code:
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). -------------------------------------------------------------------------- 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 |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
SimpleFoam k and epsilon bounded | nedved | OpenFOAM Running, Solving & CFD | 16 | March 4, 2017 08:30 |
SimpleFoam k and epsilon bounded | nedved | OpenFOAM Running, Solving & CFD | 1 | November 25, 2008 20:21 |
Compressible -> incompressible. | Jinwon | Main CFD Forum | 6 | November 23, 2007 21:07 |
Problem with LaunderSharma compressible turbulence model wrong formulation | luca | OpenFOAM Running, Solving & CFD | 4 | July 3, 2006 06:24 |
Compressible vs. Incompressible formulations | Fernando Velasco Hurtado | Main CFD Forum | 3 | January 7, 2000 16:51 |