
[Sponsors] 
Maximum number of iterations exceeded chtmultiregionsimpleFoam 

LinkBack  Thread Tools  Search this Thread  Display Modes 
May 1, 2014, 15:20 
Maximum number of iterations exceeded chtmultiregionsimpleFoam

#1 
New Member
Monssif
Join Date: Mar 2014
Posts: 5
Rep Power: 12 
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.583305e05, 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.409411e05, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 7.352814e05, Final residual = 7.352814e05, No Iterations 0 time step continuity errors : sum local = 0.002046774, global = 1.761996e05, cumulative = 1.171764e05 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.80663e05, 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.788067e05, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 8.562489e05, Final residual = 8.562489e05, No Iterations 0 time step continuity errors : sum local = 0.0005046807, global = 1.07487e05, cumulative = 2.246635e05 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.347516e05, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 6.310836e05, Final residual = 6.310836e05, No Iterations 0 time step continuity errors : sum local = 0.002436507, global = 2.080632e05, cumulative = 1.660023e06 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.236421e05, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 5.190702e05, Final residual = 5.190702e05, No Iterations 0 time step continuity errors : sum local = 0.01053975, global = 5.541302e05, cumulative = 5.375299e05 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.058582e05, No Iterations 8 GAMG: Solving for p_rgh, Initial residual = 6.857858e05, Final residual = 6.857858e05, 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.927077e05, No Iterations 4 GAMG: Solving for p_rgh, Initial residual = 4.849834e05, Final residual = 4.849834e05, No Iterations 0 time step continuity errors : sum local = 0.006161369, global = 2.698653e05, 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.22613e05, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 7.160742e05, Final residual = 7.160742e05, 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.790573e05, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 5.662955e05, Final residual = 5.662955e05, 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.450961e05, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 7.896973e05, Final residual = 7.896973e05, No Iterations 0 time step continuity errors : sum local = 0.001503887, global = 1.461686e06, 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/OpenFOAM2.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" Thank you. 

May 1, 2014, 19:19 

#2  
Senior Member
Andrew Somorjai
Join Date: May 2013
Posts: 175
Rep Power: 13 
Quote:
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; } } // ************************************************************************* // 

May 19, 2014, 13:53 

#3 
Member
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 14 
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 steadystate 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; Useful discussion links can be found below: http://www.cfdonline.com/Forums/ope...exceeded.html http://www.cfdonline.com/Forums/ope...larerror.html http://www.cfdonline.com/Forums/ope...ombustion.html Cheers Lucas 

May 19, 2014, 16:18 

#4 
Senior Member
Andrew Somorjai
Join Date: May 2013
Posts: 175
Rep Power: 13 
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.


May 19, 2014, 17:14 

#5 
Member
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 14 
Hey Andrew,
So for the specific heat I am not too concerned because I am dealing with air (1004.4 J/kgK). 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 Cheers Lucas 

May 19, 2014, 23:50 

#6 
Member
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 14 
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 

April 3, 2016, 05:52 

#7 
Member
Join Date: Oct 2015
Posts: 48
Rep Power: 11 
hi
i have this problem too can anyone help me? masoud 

April 3, 2016, 06:53 

#8 
Senior Member
Derek Mitchell
Join Date: Mar 2014
Location: UK, Reading
Posts: 172
Rep Power: 13 
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 

May 17, 2016, 00:53 
Maximum number of iterations exceeded

#9 
New Member
tuananh
Join Date: Sep 2015
Posts: 4
Rep Power: 11 
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 

May 17, 2016, 02:49 

#10 
Senior Member
Derek Mitchell
Join Date: Mar 2014
Location: UK, Reading
Posts: 172
Rep Power: 13 
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 

May 31, 2017, 18:04 

#11  
New Member
Isaac
Join Date: Feb 2013
Posts: 27
Rep Power: 13 
Quote:


December 5, 2019, 06:28 

#12 
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 7 
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 ft0 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.02112e005, 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 ft11 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.18659e005, 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.43984e006, No Iterations 2 smoothSolver: Solving for k, Initial residual = 1.0733e008, Final residual = 1.0733e008, No Iterations 0 Solving for fluid region ft1 smoothSolver: Solving for Ux, Initial residual = 0.0055058, Final residual = 2.76228e011, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00287579, Final residual = 4.33547e012, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00164844, Final residual = 2.29083e012, 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.86994e006, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00266706, Final residual = 0.000110034, No Iterations 2 Solving for fluid region ft2 smoothSolver: Solving for Ux, Initial residual = 0.00391546, Final residual = 2.45767e011, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00327532, Final residual = 8.91455e012, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00320419, Final residual = 6.21128e012, 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.90312e006, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00309988, Final residual = 0.00016859, No Iterations 2 Solving for fluid region ft3 smoothSolver: Solving for Ux, Initial residual = 0.00359363, Final residual = 1.32756e011, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00403961, Final residual = 2.24922e011, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00302194, Final residual = 1.61554e011, 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.52231e005, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00407804, Final residual = 0.000334779, No Iterations 2 Solving for fluid region ft4 smoothSolver: Solving for Ux, Initial residual = 0.0038778, Final residual = 1.00521e011, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00288863, Final residual = 4.71269e012, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00155575, Final residual = 2.28635e012, 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.85172e006, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00267175, Final residual = 0.000110143, No Iterations 2 Solving for fluid region ft5 smoothSolver: Solving for Ux, Initial residual = 0.00412002, Final residual = 4.04517e012, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00327023, Final residual = 6.08837e012, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00295683, Final residual = 5.07588e012, 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.7033e006, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.0031016, Final residual = 0.000164761, No Iterations 2 Solving for fluid region ft6 smoothSolver: Solving for Ux, Initial residual = 0.00384154, Final residual = 1.90041e011, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0040261, Final residual = 1.64686e011, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00287153, Final residual = 1.43655e011, 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.49356e005, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00406731, Final residual = 0.000330538, No Iterations 2 Solving for fluid region ft6 smoothSolver: Solving for Ux, Initial residual = 0.00439894, Final residual = 1.69643e011, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00286402, Final residual = 2.20692e012, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00175316, Final residual = 1.5467e012, 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.76803e006, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00266163, Final residual = 0.000107758, No Iterations 2 Solving for fluid region ft7 smoothSolver: Solving for Ux, Initial residual = 0.00414534, Final residual = 2.97799e012, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00324775, Final residual = 1.3963e013, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00339265, Final residual = 1.05484e013, 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.45751e005, 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.4992e006, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00308845, Final residual = 0.000160447, No Iterations 2 Solving for fluid region ft8 smoothSolver: Solving for Ux, Initial residual = 0.0037457, Final residual = 1.82422e011, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00402092, Final residual = 1.64508e011, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00317998, Final residual = 1.60725e011, 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.50375e005, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00407803, Final residual = 0.000332261, No Iterations 2 Solving for fluid region ft9 smoothSolver: Solving for Ux, Initial residual = 0.00425447, Final residual = 8.96637e005, 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.96679e005, No Iterations 4 DILUPBiCG: Solving for h, Initial residual = 0.00171053, Final residual = 3.48958e005, 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 ft10 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.03338e005, 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 s1 DICPCG: Solving for h, Initial residual = 0.00350629, Final residual = 2.34866e005, 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 ft0 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.89527e006, 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 ft11 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.1421e005, 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/OpenFOAM2.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 blueCFDCore patches for OpenFOAM. job aborted: [ranks] message [02] terminated [3] application aborted aborting MPI_COMM_WORLD (comm=0x44000000), error 1, comm rank 3 [45] terminated  error analysis  [3] on pc chtMultiRegionSimpleFoam aborted the job. abort code 1  error analysis  what does it mean?? Thankyou.. Last edited by Vishsel; December 12, 2019 at 00:09. 

December 12, 2019, 16:26 

#13 
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 22 
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 

December 13, 2019, 03:32 

#14 
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 7 
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 f1 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/heattransfer/heattransfer/0/f1/T.boundaryField.wallfluidsolid1 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/OpenFOAM2.3/src/finiteVolume/lnInclude/../fields/fvPatchFields/fvPatchField/fvPatchFieldNew.C at line 143. FOAM exiting 

December 13, 2019, 04:42 

#15 
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 22 
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 

December 17, 2019, 06:10 

#16 
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 7 
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 

December 21, 2019, 08:40 

#17 
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 22 
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


December 21, 2019, 10:01 

#18 
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 7 
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 solidfluid 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 

December 21, 2019, 11:35 

#19 
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 
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! Last edited by peterhess; December 21, 2019 at 21:50. 

December 21, 2019, 11:48 

#20 
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 7 
Thank you so much Peterhess.. I will check it and come back to you soon.


Thread Tools  Search this Thread 
Display Modes  


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