|
[Sponsors] | |||||
Maximum number of iterations exceeded chtmultiregionsimpleFoam |
![]() |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
|
|
|
#1 |
|
New Member
Monssif
Join Date: Mar 2014
Posts: 5
Rep Power: 13 ![]() |
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"
Thank you. |
|
|
|
|
|
|
|
|
#2 | |
|
Senior Member
Andrew Somorjai
Join Date: May 2013
Posts: 175
Rep Power: 14 ![]() |
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;
}
}
// ************************************************************************* //
|
||
|
|
|
||
|
|
|
#3 |
|
Member
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 15 ![]() |
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; 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 |
|
|
|
|
|
|
|
|
#4 |
|
Senior Member
Andrew Somorjai
Join Date: May 2013
Posts: 175
Rep Power: 14 ![]() |
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.
|
|
|
|
|
|
|
|
|
#5 |
|
Member
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 15 ![]() |
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 Cheers Lucas |
|
|
|
|
|
|
|
|
#6 |
|
Member
Lucas Mutti
Join Date: Aug 2013
Posts: 47
Rep Power: 15 ![]() |
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 |
|
|
|
|
|
|
|
|
#7 |
|
Member
Join Date: Oct 2015
Posts: 48
Rep Power: 12 ![]() |
hi
i have this problem too can anyone help me? masoud |
|
|
|
|
|
|
|
|
#8 |
|
Senior Member
Derek Mitchell
Join Date: Mar 2014
Location: UK, Reading
Posts: 173
Rep Power: 14 ![]() |
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 |
|
|
|
|
|
|
|
|
#9 |
|
New Member
tuananh
Join Date: Sep 2015
Posts: 4
Rep Power: 12 ![]() |
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 ![]()
|
|
|
|
|
|
|
|
|
#10 |
|
Senior Member
Derek Mitchell
Join Date: Mar 2014
Location: UK, Reading
Posts: 173
Rep Power: 14 ![]() |
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 |
|
|
|
|
|
|
|
|
#11 | |
|
New Member
Isaac
Join Date: Feb 2013
Posts: 27
Rep Power: 14 ![]() |
Quote:
|
||
|
|
|
||
|
|
|
#12 |
|
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 8 ![]() |
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 ----- what does it mean?? Thankyou.. Last edited by Vishsel; December 12, 2019 at 01:09. |
|
|
|
|
|
|
|
|
#13 |
|
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 23 ![]() |
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 |
|
|
|
|
|
|
|
|
#14 |
|
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 8 ![]() |
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
|
|
|
|
|
|
|
|
|
#15 |
|
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 23 ![]() |
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 |
|
|
|
|
|
|
|
|
#16 |
|
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 8 ![]() |
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 |
|
|
|
|
|
|
|
|
#17 |
|
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 23 ![]() |
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
|
|
|
|
|
|
|
|
|
#18 |
|
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 8 ![]() |
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 |
|
|
|
|
|
|
|
|
#19 |
|
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 18 ![]() |
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 22:50. |
|
|
|
|
|
|
|
|
#20 |
|
Member
Vishnu
Join Date: May 2019
Location: Tamilnadu, India
Posts: 55
Rep Power: 8 ![]() |
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 08:56 |
| 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 |