[kishore208]

I am using chtmultiRegionFoam solver in simple mode with temperature dependent properties. My solver is blowing up.


Solving for fluid region fluid2
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.00686722, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.00658337, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.00899513, No Iterations 1
DILUPBiCGStab: Solving for h, Initial residual = 7.9322e-05, Final residual = 4.88093e-07, No Iterations 1
Min/max T:383 383
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00279123, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 0.012005, Final residual = 0.000101146, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00169198, Final residual = 1.25692e-05, No Iterations 4
time step continuity errors : sum local = 5.39468, global = -1.24968, cumulative = -1.24968
DILUPBiCGStab: Solving for omega, Initial residual = 0.022636, Final residual = 0.000250283, No Iterations 1
DILUPBiCGStab: Solving for k, Initial residual = 1, Final residual = 0.00675252, No Iterations 1

Solving for solid region solid
DICPCG: Solving for e, Initial residual = 1, Final residual = 0.0966584, No Iterations 12
DICPCG: Solving for e, Initial residual = 0.311321, Final residual = 0.0248797, No Iterations 11
DICPCG: Solving for e, Initial residual = 0.198824, Final residual = 0.0182485, No Iterations 10
DICPCG: Solving for e, Initial residual = 0.161084, Final residual = 0.0150139, No Iterations 10
Min/max T:382.894 412.073
ExecutionTime = 7.55 s ClockTime = 9 s

Region: fluid1 Courant Number mean: 131294 max: 237647
Region: fluid2 Courant Number mean: 146.836 max: 291.079
Region: solid Diffusion Number mean: 6789.46 max: 7.1017e+06
Time = 2


Solving for fluid region fluid1
DILUPBiCGStab: Solving for Ux, Initial residual = 0.999053, Final residual = 0.000692058, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.999802, Final residual = 0.000242225, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.999056, Final residual = 0.000794289, No Iterations 1
DILUPBiCGStab: Solving for h, Initial residual = 0.996271, Final residual = 0.000992663, No Iterations 1

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const, bool) const [with Thermo = Foam::hPolynomialThermo<Foam::icoPolynomial<Foam:: specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hPolynomialThermo<Foam ::icoPolynomial<Foam::specie> >, Foam::sensibleEnthalpy>]
in file /home/ubuntu/OpenFOAM/OpenFOAM-8/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 70
Energy -> temperature conversion failed to converge:
iter Test e/h Cv/p Tnew
0 795.66 7.62999e+08 1147.88 -343887
1 -343887 -1.85787e+15 2.15655e+10 -257737
2 -257737 -5.8784e+14 9.09795e+09 -193124
3 -193124 -1.85996e+14 3.8382e+09 -144665
4 -144665 -5.88499e+13 1.61924e+09 -108321
5 -108321 -1.86203e+13 6.83117e+08 -81062.4
6 -81062.4 -5.89137e+12 2.88186e+08 -60618.2
7 -60618.2 -1.86389e+12 1.21572e+08 -45283.7
8 -45283.7 -5.8958e+11 5.12791e+07 -33779
9 -33779 -1.86388e+11 2.16203e+07 -25141.1
10 -25141.1 -5.8818e+10 9.10326e+06 -18639.5
11 -18639.5 -1.84559e+10 3.81643e+06 -13707.4
12 -13707.4 -5.68673e+09 1.57799e+06 -9870.78
13 -9870.78 -1.65005e+09 623284 -6634.06
14 -6634.06 -3.84267e+08 208236 -3024.64
15 -3024.64 -2.53638e+07 27642.9 11181.9
16 11181.9 -1.55416e+09 -600329 7981.13
17 7981.13 -3.51707e+08 -199050 4368.72
18 4368.72 -2.00471e+07 -23570.2 -12067
19 -12067 -3.50767e+09 1.09805e+06 -8537.97
20 -8537.97 -9.63585e+08 416018 -5338.75
21 -5338.75 -1.77647e+08 116348 -654.604
22 -654.604 3.25124e+07 2098.64 158893
23 158893 -8.24753e+13 -2.08567e+09 119350
24 119350 -2.60954e+13 -8.79893e+08 89691.6
25 89691.6 -8.2565e+12 -3.71203e+08 67448.1
26 67448.1 -2.6122e+12 -1.56598e+08 50764.7
27 50764.7 -8.26332e+11 -6.60578e+07 38250
28 38250 -2.61283e+11 -2.78579e+07 28857.6
29 28857.6 -8.25047e+10 -1.17378e+07 21797.4
30 21797.4 -2.59426e+10 -4.9314e+06 16462.2
31 16462.2 -8.04938e+09 -2.05247e+06 12361.4
32 12361.4 -2.39183e+09 -828199 9029.86
33 9029.86 -6.10665e+08 -299548 5764.91
34 5764.91 -7.7122e+07 -65350 -1036.44
35 -1036.44 1.57049e+07 3469.57 100314
36 100314 -1.29634e+13 -5.20656e+08 75414.7
37 75414.7 -4.10149e+12 -2.19649e+08 56740.1
38 56740.1 -1.29755e+12 -9.26588e+07 42732.7
39 42732.7 -4.10374e+11 -3.90816e+07 32222.8
40 32222.8 -1.29676e+11 -1.64747e+07 24329.3
41 24329.3 -4.08662e+10 -6.93221e+06 18381.2
42 18381.2 -1.277e+10 -2.89974e+06 13850.7
43 13850.7 -3.88366e+09 -1.18975e+06 10277.7
44 10277.7 -1.07807e+09 -457092 7115.45
45 7115.45 -2.08328e+08 -135293 2860.45
46 2860.45 -1.49956e+07 -3484.8 -106857
47 -106857 -1.764e+13 6.55965e+08 -79964.2
48 -79964.2 -5.58122e+12 2.76731e+08 -59794.6
49 -59794.6 -1.76576e+12 1.1674e+08 -44665.8
50 -44665.8 -5.58531e+11 4.92401e+07 -33315.4
51 -33315.4 -1.76564e+11 2.07599e+07 -24792.7
52 -24792.7 -5.57098e+10 8.74003e+06 -18376.5
53 -18376.5 -1.74725e+10 3.66286e+06 -13506.1
54 -13506.1 -5.37567e+09 1.51277e+06 -9709.7
55 -9709.7 -1.55191e+09 595223 -6485.26
56 -6485.26 -3.54191e+08 195828 -2800.72
57 -2800.72 -1.92722e+07 23045.3 13975.7
58 13975.7 -4.03454e+09 -1.22415e+06 10379.8
59 10379.8 -1.12552e+09 -471992 7216.91
60 7216.91 -2.22365e+08 -141942 3062.3
61 3062.3 -1.1972e+07 -5040.62 -72189.9
62 -72189.9 -3.72131e+12 2.04191e+08 -53963.5
63 -53963.5 -1.17727e+12 8.61357e+07 -40291.5
64 -40291.5 -3.72334e+11 3.63278e+07 -30032.1
65 -30032.1 -1.17651e+11 1.5311e+07 -22324
66 -22324 -3.70701e+10 6.43923e+06 -16510.1
67 -16510.1 -1.15753e+10 2.68958e+06 -12069.7
68 -12069.7 -3.51068e+09 1.09875e+06 -8540.24
69 -8540.24 -9.64531e+08 416325 -5341.11
70 -5341.11 -1.77923e+08 116485 -660.074
71 -660.074 3.21689e+07 2113.4 157936
72 157936 -8.04964e+13 -2.04802e+09 118631
73 118631 -2.54693e+13 -8.64011e+08 89153.1
74 89153.1 -8.0584e+12 -3.64503e+08 67044.2
75 67044.2 -2.54952e+12 -1.53771e+08 50461.8
76 50461.8 -8.065e+11 -6.48653e+07 38022.7
77 38022.7 -2.55008e+11 -2.73547e+07 28686.9
78 28686.9 -8.05193e+10 -1.15254e+07 21668.8
79 21668.8 -2.53145e+10 -4.84166e+06 16364.5
80 16364.5 -7.85071e+09 -2.01443e+06 12284.9
81 12284.9 -2.32908e+09 -811919 8963.82
82 8963.82 -5.91126e+08 -292404 5685.93
83 5685.93 -7.21388e+07 -62236.9 -1375.56
84 -1375.56 7.43915e+06 5358.52 65789.8
85 65789.8 -2.36204e+12 -1.4521e+08 49520.9
86 49520.9 -7.47179e+11 -6.12534e+07 37316.7
87 37316.7 -2.36239e+11 -2.58307e+07 28156.8
88 28156.8 -7.45811e+10 -1.08821e+07 21269.5
89 21269.5 -2.34359e+10 -4.56984e+06 16060.8
90 16060.8 -7.25652e+09 -1.89919e+06 12046.5
91 12046.5 -2.14144e+09 -762559 8756.51
92 8756.51 -5.32793e+08 -270714 5431.46
93 5431.46 -5.77168e+07 -52865.5 -2609
94 -2609 -1.4738e+07 19553.2 16931.7
95 16931.7 -9.05713e+09 -2.24197e+06 12728
96 12728 -2.71018e+09 -909228 9343.28
97 9343.28 -7.10014e+08 -335029 6127.55
98 6127.55 -1.03373e+08 -80947.2 312.417
99 312.417 3.08561e+08 1048.3 56388.7
100 56388.7 -1.26529e+12 -9.0926e+07 42469.1
101 42469.1 -4.00168e+11 -3.83505e+07 32025


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

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const, bool) const [with Thermo = Foam::hPolynomialThermo<Foam::icoPolynomial<Foam:: specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hPolynomialThermo<Foam ::icoPolynomial<Foam::specie> >, Foam::sensibleEnthalpy>]
in file /home/ubuntu/OpenFOAM/OpenFOAM-8/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 106.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::species::thermo<Foam::hPolynomialThermo<Foam ::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hPolynomialThermo<Foa m::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foa m::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foa m::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double) const, bool) const at ??:?
#3 Foam::species::thermo<Foam::hPolynomialThermo<Foam ::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hPolynomialThermo<Foa m::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foa m::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foa m::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double) const, bool) const at ??:?
#4 Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::singleComponentMixture<F oam:olynomialTransport<Foam::species::thermo<Foa m::hPolynomialThermo<Foam::icoPolynomial<Foam::spe cie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::calculate() at ??:?
#5 Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::singleComponentMixture<F oam:olynomialTransport<Foam::species::thermo<Foa m::hPolynomialThermo<Foam::icoPolynomial<Foam::spe cie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::correct() at ??:?
#6 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
Aborted

[fuad]

Did you solve this problem? I am also facing the same problem for rhoCentralFoam. The error is somewhat different. I tried to change the maxIter value in fvSolution but the error remained the same. I am also not sure why its saying that iteration exceeded 100 while I changed them to 1000.

The tail of the simulation looks like this.

--> FOAM FATAL ERROR: (openfoam-2106 patch=211215)
Maximum number of iterations exceeded: 100 when starting from T0:903.501 old T:588.831 new T:627.101 f:4.51284e+06 p:1.41616e+07 tol:0.0903501

    From Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hPolynomialThermo<Foam::PengRobinsonGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hPolynomialThermo<Foam::PengRobinsonGas<Foam::specie> >, Foam::sensibleEnthalpy>]
    in file /usr/src/packages/BUILD/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::exitOrAbort(int, bool) at ??:?
#2  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::PengRobinsonGas<Foam::specie>, 8>, Foam::sensibleEnthalpy>, 8> > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool) at ??:?