CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Maximum number of iterations exceeded chtmultiregionsimpleFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree2Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 1, 2014, 16:20
Default Maximum number of iterations exceeded chtmultiregionsimpleFoam
  #1
New Member
 
Monssif
Join Date: Mar 2014
Posts: 5
Rep Power: 11
Moncef is on a distinguished road
Hi everybody,
I'm trying to make a simulation using ChtMultiregionSimpleFoam, I have 3 fluid regions and 4 solid regions.
the simulation start and run for few time steps and then it gives fatal error as below:
Code:
Time = 0.0006


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.7921945, Final residual = 0.0140039, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.8933977, Final residual = 0.02410022, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.00809518, Final residual = 2.583305e-05, No Iterations 1
Min/max T:268.031 317.383
GAMG:  Solving for p_rgh, Initial residual = 0.3162556, Final residual = 0.001765501, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0003172292, Final residual = 7.409411e-05, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 7.352814e-05, Final residual = 7.352814e-05, No Iterations 0
time step continuity errors : sum local = 0.002046774, global = 1.761996e-05, cumulative = 1.171764e-05
Min/max rho:1 1

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.7179811, Final residual = 0.0156361, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.8607746, Final residual = 0.01933143, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.008032033, Final residual = 1.80663e-05, No Iterations 1
Min/max T:271.081 350.6211
GAMG:  Solving for p_rgh, Initial residual = 0.6059028, Final residual = 0.004522184, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 0.001513172, Final residual = 8.788067e-05, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 8.562489e-05, Final residual = 8.562489e-05, No Iterations 0
time step continuity errors : sum local = 0.0005046807, global = 1.07487e-05, cumulative = 2.246635e-05
Min/max rho:1 1

Solving for fluid region midAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9779546, Final residual = 0.06721385, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.9820961, Final residual = 0.06067047, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.1571378, Final residual = 0.002911904, No Iterations 1
Min/max T:295.8254 339.5323
GAMG:  Solving for p_rgh, Initial residual = 0.2170845, Final residual = 0.001151605, No Iterations 9
GAMG:  Solving for p_rgh, Initial residual = 0.0001684687, Final residual = 6.347516e-05, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 6.310836e-05, Final residual = 6.310836e-05, No Iterations 0
time step continuity errors : sum local = 0.002436507, global = -2.080632e-05, cumulative = 1.660023e-06
Min/max rho:1 1

Solving for solid region wall1
DICPCG:  Solving for h, Initial residual = 0.04669195, Final residual = 0.003361089, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03401784, Final residual = 0.000762692, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03388482, Final residual = 0.0005884278, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 273.866 max(T) [0 0 0 1 0 0 0] 273.9941

Solving for solid region wall2
DICPCG:  Solving for h, Initial residual = 0.02323714, Final residual = 0.000526991, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02318229, Final residual = 0.00044485, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02324979, Final residual = 0.0005001818, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 316.2246 max(T) [0 0 0 1 0 0 0] 320.5685

Solving for solid region wall3
DICPCG:  Solving for h, Initial residual = 0.02690195, Final residual = 0.0005019707, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02702122, Final residual = 0.000484908, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.0271, Final residual = 0.0005217254, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 315.7717 max(T) [0 0 0 1 0 0 0] 320.3459

Solving for solid region wall4
DICPCG:  Solving for h, Initial residual = 0.6028142, Final residual = 0.01325084, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.3001546, Final residual = 0.01108584, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.09297661, Final residual = 0.007668644, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 272.9786 max(T) [0 0 0 1 0 0 0] 273.006
ExecutionTime = 20.39 s  ClockTime = 21 s

Time = 0.0007


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9800156, Final residual = 0.05335714, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.9557493, Final residual = 0.04190719, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.03768833, Final residual = 0.0004185221, No Iterations 1
Min/max T:98.43921 326.1157
GAMG:  Solving for p_rgh, Initial residual = 0.3617805, Final residual = 0.002638625, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0004697464, Final residual = 5.236421e-05, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 5.190702e-05, Final residual = 5.190702e-05, No Iterations 0
time step continuity errors : sum local = 0.01053975, global = -5.541302e-05, cumulative = -5.375299e-05
Min/max rho:1 1

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9873265, Final residual = 0.02694005, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.9716581, Final residual = 0.01859983, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.02583104, Final residual = 0.0007988098, No Iterations 1
Min/max T:228.4156 385.8373
GAMG:  Solving for p_rgh, Initial residual = 0.8851469, Final residual = 0.002629838, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 0.001207289, Final residual = 7.058582e-05, No Iterations 8
GAMG:  Solving for p_rgh, Initial residual = 6.857858e-05, Final residual = 6.857858e-05, No Iterations 0
time step continuity errors : sum local = 0.005234156, global = -0.0001167688, cumulative = -0.0001705218
Min/max rho:1 1

Solving for fluid region midAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9606028, Final residual = 0.07672253, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.9498412, Final residual = 0.07503273, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.5635709, Final residual = 0.02903341, No Iterations 1
Min/max T:129.8964 500.3314
GAMG:  Solving for p_rgh, Initial residual = 0.4279954, Final residual = 0.00187809, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0003805035, Final residual = 4.927077e-05, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 4.849834e-05, Final residual = 4.849834e-05, No Iterations 0
time step continuity errors : sum local = 0.006161369, global = 2.698653e-05, cumulative = -0.0001435352
Min/max rho:1 1

Solving for solid region wall1
DICPCG:  Solving for h, Initial residual = 0.03559617, Final residual = 0.001771068, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03220624, Final residual = 0.0007789431, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03201392, Final residual = 0.0006846124, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 273.8496 max(T) [0 0 0 1 0 0 0] 273.992

Solving for solid region wall2
DICPCG:  Solving for h, Initial residual = 0.02313799, Final residual = 0.0005387188, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02303216, Final residual = 0.0005201394, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02304905, Final residual = 0.0005488426, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 316.0678 max(T) [0 0 0 1 0 0 0] 320.2905

Solving for solid region wall3
DICPCG:  Solving for h, Initial residual = 0.02722268, Final residual = 0.0006218687, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02725845, Final residual = 0.0006519343, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02727843, Final residual = 0.0006740936, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 315.5116 max(T) [0 0 0 1 0 0 0] 320.168

Solving for solid region wall4
DICPCG:  Solving for h, Initial residual = 0.8693363, Final residual = 0.05344126, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.1648544, Final residual = 0.01580605, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.1061645, Final residual = 0.009374127, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 271.8428 max(T) [0 0 0 1 0 0 0] 273.1181
ExecutionTime = 21.76 s  ClockTime = 22 s

Time = 0.0008


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.777095, Final residual = 0.03359564, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.6043706, Final residual = 0.03027286, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.05477926, Final residual = 0.0006119246, No Iterations 1
Min/max T:-1538.403 324.6424
GAMG:  Solving for p_rgh, Initial residual = 0.1068792, Final residual = 0.0005514734, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0001685353, Final residual = 7.22613e-05, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 7.160742e-05, Final residual = 7.160742e-05, No Iterations 0
time step continuity errors : sum local = 0.02476415, global = 0.0005105658, cumulative = 0.0003670305
Min/max rho:1 1

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.2123471, Final residual = 0.004336961, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.2831057, Final residual = 0.007708903, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.01623925, Final residual = 0.0002112766, No Iterations 1
Min/max T:-290.1514 407.4514
GAMG:  Solving for p_rgh, Initial residual = 0.3359547, Final residual = 0.002359242, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 0.0003813179, Final residual = 5.790573e-05, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 5.662955e-05, Final residual = 5.662955e-05, No Iterations 0
time step continuity errors : sum local = 0.0131631, global = 0.0003268909, cumulative = 0.0006939215
Min/max rho:1 1

Solving for fluid region midAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.4976941, Final residual = 0.03899976, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.5826669, Final residual = 0.04823762, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.05494696, Final residual = 0.001551656, No Iterations 1
Min/max T:-554.9509 505.8643
GAMG:  Solving for p_rgh, Initial residual = 0.2912953, Final residual = 0.001223252, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.001785834, Final residual = 7.450961e-05, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 7.896973e-05, Final residual = 7.896973e-05, No Iterations 0
time step continuity errors : sum local = 0.001503887, global = -1.461686e-06, cumulative = 0.0006924598
Min/max rho:1 1

Solving for solid region wall1
DICPCG:  Solving for h, Initial residual = 0.1972671, Final residual = 0.001589662, No Iterations 2
DICPCG:  Solving for h, Initial residual = 0.1347692, Final residual = 0.0009945654, No Iterations 2
DICPCG:  Solving for h, Initial residual = 0.1031463, Final residual = 0.009859974, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 271.879 max(T) [0 0 0 1 0 0 0] 274.4202

Solving for solid region wall2
DICPCG:  Solving for h, Initial residual = 0.0334022, Final residual = 0.0002006326, No Iterations 2
DICPCG:  Solving for h, Initial residual = 0.03326292, Final residual = 0.002789902, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03303815, Final residual = 0.002069406, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 314.5827 max(T) [0 0 0 1 0 0 0] 321.0482

Solving for solid region wall3
DICPCG:  Solving for h, Initial residual = 0.02801664, Final residual = 0.001188034, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02788079, Final residual = 0.0008740843, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02783197, Final residual = 0.0008029646, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 315.3055 max(T) [0 0 0 1 0 0 0] 320.3934

Solving for solid region wall4
DICPCG:  Solving for h, Initial residual = 0.1125148, Final residual = 0.009700842, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.06501171, Final residual = 0.004624643, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.05446337, Final residual = 0.00333266, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 270.9044 max(T) [0 0 0 1 0 0 0] 273.5055
ExecutionTime = 23.22 s  ClockTime = 24 s

Time = 0.0009


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.356727, Final residual = 0.02233952, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.6559313, Final residual = 0.05531945, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.05865561, Final residual = 0.0009535334, No Iterations 1


--> 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/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
It seems that temperature diverge! I tried a lot of solutions but with no results ( changing the relax. factor, changing schemes, checking my BCs many times, changing time step). I have no idea! Any help please.
Thank you.
Moncef is offline   Reply With Quote

Old   May 1, 2014, 20:19
Default
  #2
Senior Member
 
Andrew Somorjai
Join Date: May 2013
Posts: 175
Rep Power: 12
massive_turbulence is on a distinguished road
Quote:
Originally Posted by Moncef View Post
Hi everybody,
I'm trying to make a simulation using ChtMultiregionSimpleFoam, I have 3 fluid regions and 4 solid regions.
the simulation start and run for few time steps and then it gives fatal error as below:
Code:
Time = 0.0006


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.7921945, Final residual = 0.0140039, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.8933977, Final residual = 0.02410022, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.00809518, Final residual = 2.583305e-05, No Iterations 1
Min/max T:268.031 317.383
GAMG:  Solving for p_rgh, Initial residual = 0.3162556, Final residual = 0.001765501, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0003172292, Final residual = 7.409411e-05, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 7.352814e-05, Final residual = 7.352814e-05, No Iterations 0
time step continuity errors : sum local = 0.002046774, global = 1.761996e-05, cumulative = 1.171764e-05
Min/max rho:1 1

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.7179811, Final residual = 0.0156361, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.8607746, Final residual = 0.01933143, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.008032033, Final residual = 1.80663e-05, No Iterations 1
Min/max T:271.081 350.6211
GAMG:  Solving for p_rgh, Initial residual = 0.6059028, Final residual = 0.004522184, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 0.001513172, Final residual = 8.788067e-05, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 8.562489e-05, Final residual = 8.562489e-05, No Iterations 0
time step continuity errors : sum local = 0.0005046807, global = 1.07487e-05, cumulative = 2.246635e-05
Min/max rho:1 1

Solving for fluid region midAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9779546, Final residual = 0.06721385, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.9820961, Final residual = 0.06067047, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.1571378, Final residual = 0.002911904, No Iterations 1
Min/max T:295.8254 339.5323
GAMG:  Solving for p_rgh, Initial residual = 0.2170845, Final residual = 0.001151605, No Iterations 9
GAMG:  Solving for p_rgh, Initial residual = 0.0001684687, Final residual = 6.347516e-05, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 6.310836e-05, Final residual = 6.310836e-05, No Iterations 0
time step continuity errors : sum local = 0.002436507, global = -2.080632e-05, cumulative = 1.660023e-06
Min/max rho:1 1

Solving for solid region wall1
DICPCG:  Solving for h, Initial residual = 0.04669195, Final residual = 0.003361089, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03401784, Final residual = 0.000762692, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03388482, Final residual = 0.0005884278, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 273.866 max(T) [0 0 0 1 0 0 0] 273.9941

Solving for solid region wall2
DICPCG:  Solving for h, Initial residual = 0.02323714, Final residual = 0.000526991, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02318229, Final residual = 0.00044485, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02324979, Final residual = 0.0005001818, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 316.2246 max(T) [0 0 0 1 0 0 0] 320.5685

Solving for solid region wall3
DICPCG:  Solving for h, Initial residual = 0.02690195, Final residual = 0.0005019707, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02702122, Final residual = 0.000484908, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.0271, Final residual = 0.0005217254, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 315.7717 max(T) [0 0 0 1 0 0 0] 320.3459

Solving for solid region wall4
DICPCG:  Solving for h, Initial residual = 0.6028142, Final residual = 0.01325084, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.3001546, Final residual = 0.01108584, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.09297661, Final residual = 0.007668644, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 272.9786 max(T) [0 0 0 1 0 0 0] 273.006
ExecutionTime = 20.39 s  ClockTime = 21 s

Time = 0.0007


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9800156, Final residual = 0.05335714, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.9557493, Final residual = 0.04190719, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.03768833, Final residual = 0.0004185221, No Iterations 1
Min/max T:98.43921 326.1157
GAMG:  Solving for p_rgh, Initial residual = 0.3617805, Final residual = 0.002638625, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0004697464, Final residual = 5.236421e-05, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 5.190702e-05, Final residual = 5.190702e-05, No Iterations 0
time step continuity errors : sum local = 0.01053975, global = -5.541302e-05, cumulative = -5.375299e-05
Min/max rho:1 1

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9873265, Final residual = 0.02694005, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.9716581, Final residual = 0.01859983, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.02583104, Final residual = 0.0007988098, No Iterations 1
Min/max T:228.4156 385.8373
GAMG:  Solving for p_rgh, Initial residual = 0.8851469, Final residual = 0.002629838, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 0.001207289, Final residual = 7.058582e-05, No Iterations 8
GAMG:  Solving for p_rgh, Initial residual = 6.857858e-05, Final residual = 6.857858e-05, No Iterations 0
time step continuity errors : sum local = 0.005234156, global = -0.0001167688, cumulative = -0.0001705218
Min/max rho:1 1

Solving for fluid region midAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.9606028, Final residual = 0.07672253, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.9498412, Final residual = 0.07503273, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.5635709, Final residual = 0.02903341, No Iterations 1
Min/max T:129.8964 500.3314
GAMG:  Solving for p_rgh, Initial residual = 0.4279954, Final residual = 0.00187809, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0003805035, Final residual = 4.927077e-05, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 4.849834e-05, Final residual = 4.849834e-05, No Iterations 0
time step continuity errors : sum local = 0.006161369, global = 2.698653e-05, cumulative = -0.0001435352
Min/max rho:1 1

Solving for solid region wall1
DICPCG:  Solving for h, Initial residual = 0.03559617, Final residual = 0.001771068, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03220624, Final residual = 0.0007789431, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03201392, Final residual = 0.0006846124, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 273.8496 max(T) [0 0 0 1 0 0 0] 273.992

Solving for solid region wall2
DICPCG:  Solving for h, Initial residual = 0.02313799, Final residual = 0.0005387188, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02303216, Final residual = 0.0005201394, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02304905, Final residual = 0.0005488426, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 316.0678 max(T) [0 0 0 1 0 0 0] 320.2905

Solving for solid region wall3
DICPCG:  Solving for h, Initial residual = 0.02722268, Final residual = 0.0006218687, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02725845, Final residual = 0.0006519343, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02727843, Final residual = 0.0006740936, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 315.5116 max(T) [0 0 0 1 0 0 0] 320.168

Solving for solid region wall4
DICPCG:  Solving for h, Initial residual = 0.8693363, Final residual = 0.05344126, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.1648544, Final residual = 0.01580605, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.1061645, Final residual = 0.009374127, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 271.8428 max(T) [0 0 0 1 0 0 0] 273.1181
ExecutionTime = 21.76 s  ClockTime = 22 s

Time = 0.0008


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.777095, Final residual = 0.03359564, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.6043706, Final residual = 0.03027286, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.05477926, Final residual = 0.0006119246, No Iterations 1
Min/max T:-1538.403 324.6424
GAMG:  Solving for p_rgh, Initial residual = 0.1068792, Final residual = 0.0005514734, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.0001685353, Final residual = 7.22613e-05, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 7.160742e-05, Final residual = 7.160742e-05, No Iterations 0
time step continuity errors : sum local = 0.02476415, global = 0.0005105658, cumulative = 0.0003670305
Min/max rho:1 1

Solving for fluid region topAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.2123471, Final residual = 0.004336961, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.2831057, Final residual = 0.007708903, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.01623925, Final residual = 0.0002112766, No Iterations 1
Min/max T:-290.1514 407.4514
GAMG:  Solving for p_rgh, Initial residual = 0.3359547, Final residual = 0.002359242, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 0.0003813179, Final residual = 5.790573e-05, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 5.662955e-05, Final residual = 5.662955e-05, No Iterations 0
time step continuity errors : sum local = 0.0131631, global = 0.0003268909, cumulative = 0.0006939215
Min/max rho:1 1

Solving for fluid region midAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.4976941, Final residual = 0.03899976, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.5826669, Final residual = 0.04823762, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.05494696, Final residual = 0.001551656, No Iterations 1
Min/max T:-554.9509 505.8643
GAMG:  Solving for p_rgh, Initial residual = 0.2912953, Final residual = 0.001223252, No Iterations 7
GAMG:  Solving for p_rgh, Initial residual = 0.001785834, Final residual = 7.450961e-05, No Iterations 6
GAMG:  Solving for p_rgh, Initial residual = 7.896973e-05, Final residual = 7.896973e-05, No Iterations 0
time step continuity errors : sum local = 0.001503887, global = -1.461686e-06, cumulative = 0.0006924598
Min/max rho:1 1

Solving for solid region wall1
DICPCG:  Solving for h, Initial residual = 0.1972671, Final residual = 0.001589662, No Iterations 2
DICPCG:  Solving for h, Initial residual = 0.1347692, Final residual = 0.0009945654, No Iterations 2
DICPCG:  Solving for h, Initial residual = 0.1031463, Final residual = 0.009859974, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 271.879 max(T) [0 0 0 1 0 0 0] 274.4202

Solving for solid region wall2
DICPCG:  Solving for h, Initial residual = 0.0334022, Final residual = 0.0002006326, No Iterations 2
DICPCG:  Solving for h, Initial residual = 0.03326292, Final residual = 0.002789902, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.03303815, Final residual = 0.002069406, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 314.5827 max(T) [0 0 0 1 0 0 0] 321.0482

Solving for solid region wall3
DICPCG:  Solving for h, Initial residual = 0.02801664, Final residual = 0.001188034, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02788079, Final residual = 0.0008740843, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.02783197, Final residual = 0.0008029646, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 315.3055 max(T) [0 0 0 1 0 0 0] 320.3934

Solving for solid region wall4
DICPCG:  Solving for h, Initial residual = 0.1125148, Final residual = 0.009700842, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.06501171, Final residual = 0.004624643, No Iterations 1
DICPCG:  Solving for h, Initial residual = 0.05446337, Final residual = 0.00333266, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 270.9044 max(T) [0 0 0 1 0 0 0] 273.5055
ExecutionTime = 23.22 s  ClockTime = 24 s

Time = 0.0009


Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.356727, Final residual = 0.02233952, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.6559313, Final residual = 0.05531945, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.05865561, Final residual = 0.0009535334, No Iterations 1


--> 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/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
It seems that temperature diverge! I tried a lot of solutions but with no results ( changing the relax. factor, changing schemes, checking my BCs many times, changing time step). I have no idea! Any help please.
Thank you.

I got errors like these running buoyantPimpleFoam recently even with super small timesteps (deltaT .0000001 ) and I still got this error.

When I played around with my pressure variables (I had huge differences of initial pressure and final pressure because I used a BC table to setup the system inlet in 0/p_rgh) I noticed that I had a problem when the internalField (declared on the top in 0/p_rgh) got close to the pressure in the inlet, usually 3/4 through the simulation.

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{    
    inlet
    {
        type            uniformTotalPressure;// this is supposed to be  buoyantPressure;
        pressure        table
        (
            (0.000010  202650.000000)        
            (0.010000  0.000000)
        );
        p0              0;
        U               U;
        phi             phi;      
        psi             none;
        gamma           1;       
    }
    outlet
    {
        type            buoyantPressure;      
        value           uniform 242650;     
    }    
    wall
    {
        type            buoyantPressure;    
        value           uniform 0;
    }
    cylinder_wall
    {
        type            buoyantPressure;        
        value           uniform 0;         
    }
    bottom
    {
        type            buoyantPressure;    
        value           uniform 0;   
    }
}

// ************************************************************************* //
Hamoon likes this.
massive_turbulence is offline   Reply With Quote

Old   May 19, 2014, 14:53
Default
  #3
Member
 
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 13
lramutti is on a distinguished road
Hey guys,

I just came across this same issue and I am not sure if you have solved it but it seems that I got it to work. Try using the following script for the controlDict instead and test it for both steady-state and transient problems by changing the fvSchemes to steadyState and Euler. Let me know how it goes for you guys.

Code:
application     chtMultiRegionSimpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         1000;

deltaT          0.011;

writeControl    runTime;

writeInterval   0.05;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep off;

maxCo 0.3;
** Please note that the name of the application is irrelevant therefore there is no need to change it between a transient and steady state solver.

Useful discussion links can be found below:
http://www.cfd-online.com/Forums/ope...-exceeded.html

http://www.cfd-online.com/Forums/ope...lar-error.html

http://www.cfd-online.com/Forums/ope...ombustion.html

Cheers

Lucas
lramutti is offline   Reply With Quote

Old   May 19, 2014, 17:18
Default
  #4
Senior Member
 
Andrew Somorjai
Join Date: May 2013
Posts: 175
Rep Power: 12
massive_turbulence is on a distinguished road
If you try to use paraview and plot your simulation up to the point of crashing you can almost always see something going out of whack, in my case the pressure variable along the inlet was negative??? You should also check your your Cp (specific heat capacity at constant pressure) because if it's too small that could be another fault.
massive_turbulence is offline   Reply With Quote

Old   May 19, 2014, 18:14
Default
  #5
Member
 
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 13
lramutti is on a distinguished road
Hey Andrew,

So for the specific heat I am not too concerned because I am dealing with air (1004.4 J/kg-K). After doing some research, I have realized that this problem has multiple reasons such as: mesh refinement and settings in the systems folder. My first suggestion for you is to make sure you have a refined mesh. Type in the command line:

Code:
checkMesh -allTopology -allGeometry
For my case, I notice that there are 3 failures that were not pointed out on Fluent. If even after solving the mesh the problem persists, then I would play around with the systems scripts.

Cheers

Lucas
lramutti is offline   Reply With Quote

Old   May 20, 2014, 00:50
Default
  #6
Member
 
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 13
lramutti is on a distinguished road
Hey Andrew,

I am still trying to solve my problem but after some testing I have noticed that the only way to control the code is by under relaxing the parameters. Still, this is not the solution of the problem but a simple "artificial delay". I have been playing with the preconditioners but I am not getting good results. Let me know if you figure this out.

Cheers
lramutti is offline   Reply With Quote

Old   April 3, 2016, 06:52
Default
  #7
Member
 
Join Date: Oct 2015
Posts: 48
Rep Power: 10
masoudsh is on a distinguished road
hi
i have this problem too can anyone help me?

masoud
masoudsh is offline   Reply With Quote

Old   April 3, 2016, 07:53
Default
  #8
Senior Member
 
Derek Mitchell
Join Date: Mar 2014
Location: UK, Reading
Posts: 172
Rep Power: 12
derekm is on a distinguished road
If you are using chtmultiregionsimplefoam for steady state why the small time steps?
in the air regions are you using a zero or close to zero initial velocity? zero velocity can be a source of instability.
have you constrained the heat flux either explicitly or implicitly by thermal diffusivity?
solutions which do not prevent near infinite heat flux can be unstable.
__________________
A CHEERING BAND OF FRIENDLY ELVES CARRY THE CONQUERING ADVENTURER OFF INTO THE SUNSET
derekm is offline   Reply With Quote

Old   May 17, 2016, 01:53
Default Maximum number of iterations exceeded
  #9
New Member
 
tuananh
Join Date: Sep 2015
Posts: 4
Rep Power: 10
tuananhtran56 is on a distinguished road
Hi all.
I'm using chtMultiRegionSimpleFoam to simulating heat transfer between solid and fluid. I got the same error of maximum no of iteration exceeded. Can anybody plz suggest me where i am wrong ??? thanks in advance .
Here is my case file and report. the coil_solid is inside the tank_water
Attached Images
File Type: jpg Screenshot from 2016-05-17 11:31:40.jpg (99.6 KB, 68 views)
File Type: jpg Screenshot from 2016-05-17 11:31:50.jpg (95.6 KB, 57 views)
Attached Files
File Type: gz 0.tar.gz (1.7 KB, 13 views)
File Type: gz Desktop.tar.gz (3.2 KB, 14 views)
tuananhtran56 is offline   Reply With Quote

Old   May 17, 2016, 03:49
Default
  #10
Senior Member
 
Derek Mitchell
Join Date: Mar 2014
Location: UK, Reading
Posts: 172
Rep Power: 12
derekm is on a distinguished road
i use this to help me fix things like this
http://www.dicat.unige.it/guerrero/o...sandtricks.pdf
__________________
A CHEERING BAND OF FRIENDLY ELVES CARRY THE CONQUERING ADVENTURER OFF INTO THE SUNSET
derekm is offline   Reply With Quote

Old   May 31, 2017, 19:04
Default
  #11
New Member
 
Isaac
Join Date: Feb 2013
Posts: 27
Rep Power: 13
isaac1326 is on a distinguished road
Quote:
Originally Posted by derekm View Post
i use this to help me fix things like this
http://www.dicat.unige.it/guerrero/o...sandtricks.pdf
The link doesn't work
isaac1326 is offline   Reply With Quote

Old   December 5, 2019, 07:28
Default
  #12
Member
 
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 6
Vishsel is on a distinguished road
Hi everyone,

My case has 13 fluid region and 1 solid region & i am using chtMultiRegionSimpleFoam solver. Openfoam v2.3. and about the case are solid fluid coupling. Heat flux applied on the solid body and need to check the outlet fluid temperature as how much heat rejected.

And my problem is after some iterations, temperature getting diverged in one fluid region.

Error is mentioned below..

Code:
Time = 77


Solving for fluid region f-t-0
smoothSolver:  Solving for Ux, Initial residual = 0.12817, Final residual = 0.00164259, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0932991, Final residual = 0.00111304, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.121568, Final residual = 0.00146671, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.127605, Final residual = 0.0119586, No Iterations 1
Min/max T:298.231 3016.98
GAMG:  Solving for p_rgh, Initial residual = 0.26606, Final residual = 0.00188142, No Iterations 5
time step continuity errors : sum local = 1044.55, global = 127.369, cumulative = -1026.26
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.107336, Final residual = 1.02112e-005, No Iterations 2
bounding epsilon, min: -728524 max: 1.9846e+010 average: 335954
smoothSolver:  Solving for k, Initial residual = 0.146899, Final residual = 0.00187524, No Iterations 2

Solving for fluid region f-t-11
smoothSolver:  Solving for Ux, Initial residual = 0.0672419, Final residual = 0.000491632, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0404616, Final residual = 7.18659e-005, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0586714, Final residual = 0.000133782, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.0440309, Final residual = 0.00214984, No Iterations 1
Min/max T:-63.5777 1052.14
GAMG:  Solving for p_rgh, Initial residual = 0.0714801, Final residual = 0.000530161, No Iterations 7
time step continuity errors : sum local = 449.977, global = -83.8794, cumulative = -1110.14
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.068333, Final residual = 3.43984e-006, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 1.0733e-008, Final residual = 1.0733e-008, No Iterations 0

Solving for fluid region f-t-1
smoothSolver:  Solving for Ux, Initial residual = 0.0055058, Final residual = 2.76228e-011, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00287579, Final residual = 4.33547e-012, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00164844, Final residual = 2.29083e-012, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.00692024, Final residual = 0.000186959, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0245849, Final residual = 0.000212541, No Iterations 2
time step continuity errors : sum local = 0.407565, global = -0.0970675, cumulative = -1110.24
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000292477, Final residual = 2.86994e-006, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00266706, Final residual = 0.000110034, No Iterations 2

Solving for fluid region f-t-2
smoothSolver:  Solving for Ux, Initial residual = 0.00391546, Final residual = 2.45767e-011, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00327532, Final residual = 8.91455e-012, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00320419, Final residual = 6.21128e-012, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.00610461, Final residual = 0.000224566, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0254049, Final residual = 0.000210484, No Iterations 2
time step continuity errors : sum local = 0.149835, global = -0.0478098, cumulative = -1110.29
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000351238, Final residual = 4.90312e-006, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00309988, Final residual = 0.00016859, No Iterations 2

Solving for fluid region f-t-3
smoothSolver:  Solving for Ux, Initial residual = 0.00359363, Final residual = 1.32756e-011, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00403961, Final residual = 2.24922e-011, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00302194, Final residual = 1.61554e-011, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.007511, Final residual = 0.000433068, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0271094, Final residual = 0.000245011, No Iterations 2
time step continuity errors : sum local = 0.273242, global = -0.0899534, cumulative = -1110.38
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000605403, Final residual = 1.52231e-005, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00407804, Final residual = 0.000334779, No Iterations 2

Solving for fluid region f-t-4
smoothSolver:  Solving for Ux, Initial residual = 0.0038778, Final residual = 1.00521e-011, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00288863, Final residual = 4.71269e-012, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00155575, Final residual = 2.28635e-012, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.0069016, Final residual = 0.000186978, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0245686, Final residual = 0.000217835, No Iterations 2
time step continuity errors : sum local = 0.419197, global = -0.0826873, cumulative = -1110.46
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000292629, Final residual = 2.85172e-006, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00267175, Final residual = 0.000110143, No Iterations 2

Solving for fluid region f-t-5
smoothSolver:  Solving for Ux, Initial residual = 0.00412002, Final residual = 4.04517e-012, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00327023, Final residual = 6.08837e-012, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00295683, Final residual = 5.07588e-012, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.00617927, Final residual = 0.000274063, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0254413, Final residual = 0.000214345, No Iterations 2
time step continuity errors : sum local = 0.153339, global = -0.0464748, cumulative = -1110.51
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000351493, Final residual = 4.7033e-006, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0031016, Final residual = 0.000164761, No Iterations 2

Solving for fluid region f-t-6
smoothSolver:  Solving for Ux, Initial residual = 0.00384154, Final residual = 1.90041e-011, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0040261, Final residual = 1.64686e-011, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00287153, Final residual = 1.43655e-011, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.00775241, Final residual = 0.000240444, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0271797, Final residual = 0.000232384, No Iterations 2
time step continuity errors : sum local = 0.257471, global = -0.108293, cumulative = -1110.61
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000599872, Final residual = 1.49356e-005, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00406731, Final residual = 0.000330538, No Iterations 2

Solving for fluid region f-t-6
smoothSolver:  Solving for Ux, Initial residual = 0.00439894, Final residual = 1.69643e-011, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00286402, Final residual = 2.20692e-012, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00175316, Final residual = 1.5467e-012, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.006988, Final residual = 0.000260046, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0246083, Final residual = 0.000204224, No Iterations 2
time step continuity errors : sum local = 0.39059, global = -0.0783273, cumulative = -1110.69
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000291615, Final residual = 2.76803e-006, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00266163, Final residual = 0.000107758, No Iterations 2

Solving for fluid region f-t-7
smoothSolver:  Solving for Ux, Initial residual = 0.00414534, Final residual = 2.97799e-012, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00324775, Final residual = 1.3963e-013, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00339265, Final residual = 1.05484e-013, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.00612513, Final residual = 0.000182331, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0253415, Final residual = 9.45751e-005, No Iterations 3
time step continuity errors : sum local = 0.0668698, global = -0.0170405, cumulative = -1110.71
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000348788, Final residual = 4.4992e-006, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00308845, Final residual = 0.000160447, No Iterations 2

Solving for fluid region f-t-8
smoothSolver:  Solving for Ux, Initial residual = 0.0037457, Final residual = 1.82422e-011, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00402092, Final residual = 1.64508e-011, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.00317998, Final residual = 1.60725e-011, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.00776677, Final residual = 0.000241656, No Iterations 1
Min/max T:300 323.15
GAMG:  Solving for p_rgh, Initial residual = 0.0271495, Final residual = 0.000194803, No Iterations 2
time step continuity errors : sum local = 0.215836, global = -0.0726273, cumulative = -1110.78
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.000601963, Final residual = 1.50375e-005, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.00407803, Final residual = 0.000332261, No Iterations 2

Solving for fluid region f-t-9
smoothSolver:  Solving for Ux, Initial residual = 0.00425447, Final residual = 8.96637e-005, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 0.00588838, Final residual = 0.00012303, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.00502057, Final residual = 9.96679e-005, No Iterations 4
DILUPBiCG:  Solving for h, Initial residual = 0.00171053, Final residual = 3.48958e-005, No Iterations 2
Min/max T:300 350.106
GAMG:  Solving for p_rgh, Initial residual = 0.0816832, Final residual = 0.000427319, No Iterations 3
time step continuity errors : sum local = 0.0076865, global = -0.00186211, cumulative = -1110.78
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.0326387, Final residual = 0.00194519, No Iterations 8
smoothSolver:  Solving for k, Initial residual = 0.0690091, Final residual = 0.00566434, No Iterations 8

Solving for fluid region f-t-10
smoothSolver:  Solving for Ux, Initial residual = 0.0114538, Final residual = 0.00101102, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0155309, Final residual = 0.000325799, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.00695301, Final residual = 0.000651431, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.00351635, Final residual = 8.03338e-005, No Iterations 2
Min/max T:300 327.784
GAMG:  Solving for p_rgh, Initial residual = 0.0816719, Final residual = 0.000632341, No Iterations 2
time step continuity errors : sum local = 0.031818, global = 0.00420144, cumulative = -1110.78
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.00403982, Final residual = 0.000315038, No Iterations 6
smoothSolver:  Solving for k, Initial residual = 0.0395869, Final residual = 0.00316137, No Iterations 8

Solving for solid region s-1
DICPCG:  Solving for h, Initial residual = 0.00350629, Final residual = 2.34866e-005, No Iterations 4
Min/max T:min(T) [0 0 0 1 0 0 0] 323.151 max(T) [0 0 0 1 0 0 0] 329.967
ExecutionTime = 706.244 s  ClockTime = 707 s

Time = 78


Solving for fluid region f-t-0
smoothSolver:  Solving for Ux, Initial residual = 0.128627, Final residual = 0.00351613, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0872831, Final residual = 0.00155013, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.108134, Final residual = 0.00238099, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.113183, Final residual = 0.000762394, No Iterations 3
Min/max T:298.361 2928.11
GAMG:  Solving for p_rgh, Initial residual = 0.291059, Final residual = 0.00282524, No Iterations 4
time step continuity errors : sum local = 1714.97, global = 399.552, cumulative = -711.228
Min/max rho:
smoothSolver:  Solving for epsilon, Initial residual = 0.115352, Final residual = 9.89527e-006, No Iterations 2
bounding epsilon, min: -1.06864e+006 max: 2.63891e+010 average: 407062
smoothSolver:  Solving for k, Initial residual = 0.140569, Final residual = 0.00155209, No Iterations 2
bounding k, min: -0.215904 max: 21616.5 average: 1.51868

Solving for fluid region f-t-11
smoothSolver:  Solving for Ux, Initial residual = 0.0594532, Final residual = 0.000278845, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.0389371, Final residual = 5.1421e-005, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0569155, Final residual = 0.000106508, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.0381555, Final residual = 0.00181037, No Iterations 1
[3]
[3]
[3] --> FOAM FATAL ERROR:
[3] Maximum number of iterations exceeded
[3]
[3]     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
[3]     in file P:/OPENFO~1.3/OpenFOAM-2.3/src/thermophysicalModels/specie/lnInclude/../thermo/thermo/thermoI.H at line 76.
[3]
FOAM parallel run aborting
[3]
[3] We're sorry, but the application crashed and safe stack tracing isn't available in this current implementation of blueCFD-Core patches for OpenFOAM.

job aborted:
[ranks] message

[0-2] terminated

[3] application aborted
aborting MPI_COMM_WORLD (comm=0x44000000), error 1, comm rank 3

[4-5] terminated

---- error analysis -----

[3] on pc
chtMultiRegionSimpleFoam aborted the job. abort code 1

---- error analysis -----
Anybody got the solution for ''Maximum number of iterations exceeded'' ??
what does it mean??

Thankyou..
Attached Files
File Type: zip case_file.zip (18.0 KB, 14 views)

Last edited by Vishsel; December 12, 2019 at 01:09.
Vishsel is offline   Reply With Quote

Old   December 12, 2019, 17:26
Default
  #13
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 21
jherb is on a distinguished road
You could try to use a fixedTemperatureConstraint to start up your simulation. Perhaps also using dpdt off might help.



https://github.com/OpenFOAM/OpenFOAM...olid/fvOptions
jherb is offline   Reply With Quote

Old   December 13, 2019, 04:32
Default
  #14
Member
 
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 6
Vishsel is on a distinguished road
Hi jherb,

Thanks for your reply.. But with fixedTemperatureConstraint is not working. I am getting unknown patchField type error..

Code:
*** Reading fluid mesh thermophysical properties for region f-1

    Adding to thermoFluid

Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}



--> FOAM FATAL IO ERROR:
Unknown patchField type fixedTemperatureConstraint for patch type wall

Valid patchField types are :

102
(
MarshakRadiation
MarshakRadiationFixedTemperature
advective
calculated
codedFixedValue
codedMixed
compressible::alphatJayatillekeWallFunction
compressible::alphatWallFunction
compressible::epsilonLowReWallFunction
compressible::epsilonWallFunction
compressible::fWallFunction
compressible::kLowReWallFunction
compressible::kqRWallFunction
compressible::omegaWallFunction
compressible::thermalBaffle1D<hConstSolidThermoPhysics>
compressible::thermalBaffle1D<hExponentialSolidThermoPhysics>
compressible::turbulentHeatFluxTemperature
compressible::turbulentMixingLengthDissipationRateInlet
compressible::turbulentMixingLengthFrequencyInlet
compressible::turbulentTemperatureCoupledBaffleMixed
compressible::turbulentTemperatureRadCoupledMixed
compressible::v2WallFunction
convectiveHeatTransfer
cyclic
cyclicACMI
cyclicAMI
cyclicSlip
directionMixed
empty
energyJump
energyJumpAMI
externalCoupled
externalCoupledTemperature
externalWallHeatFluxTemperature
fan
fanPressure
fixedEnergy
fixedFluxPressure
fixedGradient
fixedInternalValue
fixedJump
fixedJumpAMI
fixedMean
fixedPressureCompressibleDensity
fixedUnburntEnthalpy
fixedValue
freestream
freestreamPressure
gradientEnergy
gradientUnburntEnthalpy
greyDiffusiveRadiation
greyDiffusiveRadiationViewFactor
inletOutlet
inletOutletTotalTemperature
mapped
mappedField
mappedFixedInternalValue
mappedFixedPushedInternalValue
mixed
mixedEnergy
mixedUnburntEnthalpy
mutLowReWallFunction
mutURoughWallFunction
mutUSpaldingWallFunction
mutUWallFunction
mutkRoughWallFunction
mutkWallFunction
nonuniformTransformCyclic
oscillatingFixedValue
outletInlet
outletMappedUniformInlet
partialSlip
phaseHydrostaticPressure
prghPressure
processor
processorCyclic
rotatingTotalPressure
sliced
slip
symmetry
symmetryPlane
syringePressure
timeVaryingMappedFixedValue
totalFlowRateAdvectiveDiffusive
totalPressure
totalTemperature
turbulentInlet
turbulentIntensityKineticEnergyInlet
uniformDensityHydrostaticPressure
uniformFixedGradient
uniformFixedValue
uniformInletOutlet
uniformJump
uniformJumpAMI
uniformTotalPressure
variableHeightFlowRate
wallHeatTransfer
waveSurfacePressure
waveTransmissive
wedge
wideBandDiffusiveRadiation
zeroGradient
)


file: c:/case/heat-transfer/heattransfer/0/f-1/T.boundaryField.wall-fluid-solid-1 from line 30 to line 34.

    From function fvPatchField<Type>::New(const fvPatch&, const DimensionedField<Type, volMesh>&, const dictionary&)
    in file P:/OPENFO~1.3/OpenFOAM-2.3/src/finiteVolume/lnInclude/../fields/fvPatchFields/fvPatchField/fvPatchFieldNew.C at line 143.

FOAM exiting
Vishsel is offline   Reply With Quote

Old   December 13, 2019, 05:42
Default
  #15
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 21
jherb is on a distinguished road
fixedTemperatureConstraint must by used in fvOptions.


Have a look at the different tutorials, where it is used, e.g. https://github.com/OpenFOAM/OpenFOAM...olid/fvOptions


see also: https://caefn.com/openfoam/fvoptions
jherb is offline   Reply With Quote

Old   December 17, 2019, 07:10
Default
  #16
Member
 
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 6
Vishsel is on a distinguished road
Hi jherb,

By using fixedTemperatureConstraint boundary type, just the temperature values are fixed at cell zones and for every iterations the temperature value are fixed as 'Min/max T:300 300'. But my problem was temperature tends to diverge after some iterations. That showing the error ''Maximum number of iterations exceeded''.
Whether the boundary types are correct for my boundary condition in attached case file??
Any other clue to solve my error??

Thanks in advance
Vishsel
Vishsel is offline   Reply With Quote

Old   December 21, 2019, 09:40
Default
  #17
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 21
jherb is on a distinguished road
I am not sure what you are doing. fixedTemperatureConstraint is not a boundary condition/type. It has to be applied in a fvOptions file. See e.g. Building a solver with fixedTemperatureConstraint using fvOptions
jherb is offline   Reply With Quote

Old   December 21, 2019, 11:01
Default
  #18
Member
 
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 6
Vishsel is on a distinguished road
Hi Jherb,

Sorry for not seeing this problem before I have posting.. Exactly this is not a problem.. because the solid and fluid bodies are not even coupled.. there is no interface formation between both of it..
So I need a clue or solution to create an interface between solid and fluid body.
I had both solid and fluid body. I could get solid - solid interface as mappedWall by using any of this commands splitMeshRegions -cellZones -overwrite & splitMeshRegions -cellZonesOnly -useFcaeZones
But i don't get solid-fluid interface by those commands. But my mesh was also conformal. Only i am getting wall type, instead of mappedWall.
May i know why this problem occurs ?? what is the solution for it ?

Thanks in advance
Vishsel
Vishsel is offline   Reply With Quote

Old   December 21, 2019, 12:35
Default
  #19
Senior Member
 
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 16
peterhess is on a distinguished road
Hello Vishsel!

Unfortunately you did not include the mesh in your case, so I could understand, which boundary is inlet, outlet, wall, or mapped wall...

That why I will try to give a general recommendation, how to setup the boundaries.

---------------------------------------------------

FLUIDS!

---------------------------------------------------

Temperature T:

- Inlets:

fixedValue is ok as you defined


- Outlets:

fixedValue or zeroGradient is not recommended!

Instead use:

type inletOutlet;
value uniform 300;
inletValue uniform 300;

zeroGradient ist not recommended cause it is less stable than inletOutlet!

- Walls:

zeroGradient is acceptable.

fixedValue is also acceptable if the value known


- Mapped walls without radiation:

type compressible::turbulentTemperatureCoupledBaffleMix ed; // there is no space between Mix and ed! --> Mixed
Tnbr T;
kappaMethod fluidThermo;
value uniform 300;

---------------------------------------------------------

Pressure p:

All boundaries should be defined as:

type calculated;
value uniform 1e5; // don't define zero!!!

--------------------------------------------------------

Pressure p_rgh:

- Inlets:

Zero Gradient is not ok, as in the case you uploaded!

Instead use:

type fixedFluxPressure;
value uniform 1e5;

- Outlets:

type fixedValue;
value uniform 1e5;

- Walls and mapped walles:

type fixedFluxPressure;
value uniform 1e5;

--------------------------------------------------

Velocity U:

- Inlets:

flowRateInletVelocity is ok, but it is better to start with fixedValue, be sure that the simulatiuon works, and then change to flowRateInletVelocity, cause this one is more complex to be defined...

type fixedValue;
value uniform (0.5 0 0);

- Outlets:

The best is to use:

type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);

zeroGradient ist not recommended cause it is less stable than inletOutlet!

- Walls:

type fixedValue;
value uniform (0 0 0);

------------------------------------------------------------------

Turn turbulence and radiation off at the start to be sure that the simulation runs and then turn them on after being sure that the simulation works. I will not comment them here cause I want to keep my answer as short as possible!

------------------------------------------------------------------

SOLIDS:

-----------------------------------------------------------------

Temperature T:

- Walls solids to fluids - or - solids to solids:

No radiation:

type compressible::turbulentTemperatureCoupledBaffleMix ed;
Tnbr T;
kappaMethod solidThermo;
value uniform 300;

- Wall not neighbor to fluid:

fixedValue is ok!

or as you defined is also ok!

-----------------------------------------------------------------

Pressure p:

All boundaries as:

type calculated;
value uniform 1e5;


If the wall is an outer boundary (not neighbor to a fluid):

type fixedValue;
value uniform 1e5;

----------------------------------------------------------------

All the values given here are just examples and could be changed as needed!

----------------------------------------------------------------

Most huge mistake is to define a high value in p_rgh (as example 1e5) and 0 in p!

you will get in most cases a divergence!

----------------------------------------------------------------

Give please a feedback if those suggestions works and please be free to ask anything you want.

Regards

Peter

PS: send the case complete to fix it for you, if there still a problem!
PSander likes this.

Last edited by peterhess; December 21, 2019 at 22:50.
peterhess is offline   Reply With Quote

Old   December 21, 2019, 12:48
Default
  #20
Member
 
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 6
Vishsel is on a distinguished road
Thank you so much Peterhess.. I will check it and come back to you soon.
Vishsel is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Foam::error::PrintStack almir OpenFOAM Running, Solving & CFD 91 December 21, 2022 05:50
buoyantSimpleFoam and watertank Tobi OpenFOAM Running, Solving & CFD 100 December 18, 2022 09:15
Help for the small implementation in turbulence model shipman OpenFOAM Programming & Development 25 March 19, 2014 11:08
should Courant number always be kept below 1? wc34071209 OpenFOAM Running, Solving & CFD 16 March 9, 2014 20:31
Unstabil Simulation with chtMultiRegionFoam mbay101 OpenFOAM Running, Solving & CFD 13 December 28, 2013 14:12


All times are GMT -4. The time now is 06:44.