|
[Sponsors] |
December 3, 2012, 08:19 |
Problem in solving
|
#1 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
this error tells the maximum number of iterations the solver is reached probably in e equation solving.why this is occuring?
Im engaged in this problem but it don't want to be resolved! fluxScheme: Kurganov Starting time loop Mean and max Courant Numbers = 0.10019 0.211345 Time = 1e-09 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 swak4Foam: Allocating new repository for sampledGlobalVariables smoothSolver: Solving for Ux, Initial residual = 0.878819, Final residual = 9.93747e-10, No Iterations 438 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 9.87651e-10, No Iterations 693 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for e, Initial residual = 0.94115, Final residual = 5.58403e-06, No Iterations 1002 --> FOAM FATAL ERROR: Maximum number of iterations exceeded 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/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 Foam::specieThermo<Foam::hConstThermo<Foam:erfec tGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam:erfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam:erfe ctGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam:erfe ctGas> >::*)(double) const) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #3 Foam::ePsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #4 Foam::ePsiThermo<Foam:ureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::correct() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #5 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoCentralFoam" #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/rhoCentralFoam" Aborted |
|
January 29, 2013, 01:10 |
|
#2 |
Member
Anant Diwakar
Join Date: Jan 2013
Posts: 68
Rep Power: 13 |
Hi
I am also having the same problem. Did you find the cause of its occurrence ? Thanks Anant |
|
January 29, 2013, 09:51 |
|
#3 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
hi dear Anant.I did some changes in fvSchemes and fvSolutions.using CrankNicholson .5;and div(phi,U) Gauss upwind; is better than Euler and linear respectively.
a very important thing is boundary conditions and initial conditions.I changed initial condition and increased it because it had made a supersonic flow at inlet but it must be subsonic in my problem. and also try another solver and change maxIter for e in specieThermo.c (I think it was there) can lead to more stability. whats the problem you want to model?comp or incomp? |
|
January 30, 2013, 00:45 |
|
#4 |
Member
Anant Diwakar
Join Date: Jan 2013
Posts: 68
Rep Power: 13 |
Dear Ehsan
Thanks a lot for your help. I tweaked the initial conditions and it's working now. I am solving compressible flow between rotating cylinders using rhoCentralFoam. Regards Anant |
|
January 30, 2013, 04:57 |
|
#5 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
Dear Anant
Im working with rhoCentralFoam and rhoPimpleFoam. Im interested in the problem you are working on.could you send me your case please? Thanks. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
AMI speed performance | danny123 | OpenFOAM | 21 | October 24, 2020 04:13 |
alphaEqn.H in twoPhaseEulerFoam | cheng1988sjtu | OpenFOAM Bugs | 15 | May 1, 2016 16:12 |
pimpleFoam: turbulence->correct(); is not executed when using residualControl | hfs | OpenFOAM Running, Solving & CFD | 3 | October 29, 2013 08:35 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 15:33 |
Problem with SimpleFoam for a solution around an OneraM6 wing | gara1988 | OpenFOAM Running, Solving & CFD | 3 | November 13, 2012 05:47 |