sasanghomi 
September 3, 2013 18:31 
Bounding k and epsilon & coldEngineFoam
Hi Dear foamers,
I am using coldEngineFoam for simulation of engine fluid flow and I have a problem for getting convergence. Several last time steps are as below :
Code:
Courant Number mean: 4.24885e05 max: 29.1764 velocity magnitude: 2.61013e+13
deltaT = 2.86066e17
MassFlows: presin = 0.00795114
> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 692
Increased the timePrecision from 12 to 13 to distinguish between timeNames at time 0.00129011
Crank angle = 23.2219 CAdeg
deltaZ = 2.22045e16
DICPCG: Solving for cellMotionUz, Initial residual = 9.58306e06, Final residual = 9.58306e06, No Iterations 0
clearance: 0.0121013
Piston speed = 7.76201 m/s
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.953778, Final residual = 2.06381e06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.999997, Final residual = 1.14924e09, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.997449, Final residual = 1.51092e08, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.334811, Final residual = 4.57511e08, No Iterations 2
DICPCG: Solving for p, Initial residual = 1.10574e05, Final residual = 2.18399e21, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.27216e07, global = 9.61115e09, cumulative = 2.47977e07
DILUPBiCG: Solving for h, Initial residual = 0.123253, Final residual = 5.05499e08, No Iterations 2
DICPCG: Solving for p, Initial residual = 2.87271e05, Final residual = 1.66698e19, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0144206, global = 0.0141697, cumulative = 0.0141699
bounding epsilon, min: 7.4933e20 max: 1.62696e+51 average: 4.68755e+46
DILUPBiCG: Solving for epsilon, Initial residual = 0.000725064, Final residual = 1.36341e20, No Iterations 1
bounding epsilon, min: 1.48494e+56 max: 1.22189e+58 average: 1.12163e+53
DILUPBiCG: Solving for k, Initial residual = 0.690174, Final residual = 4.29296e08, No Iterations 2
bounding k, min: 2.58323e+35 max: 4.41432e+35 average: 1.61653e+31
Mean pressure:8.54579e+22
Mean temperature:109312
Mean u':1.74587e+12
ExecutionTime = 1380.08 s ClockTime = 1384 s
Courant Number mean: 0.139494 max: 181009 velocity magnitude: 4.11781e+17
deltaT = 1.26432e22
MassFlows: presin = 0.00795114
> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 692
Increased the timePrecision from 13 to 14 to distinguish between timeNames at time 0.00129011
Crank angle = 23.2219 CAdeg
deltaZ = 0
DICPCG: Solving for cellMotionUz, Initial residual = 9.58306e06, Final residual = 9.58306e06, No Iterations 0
clearance: 0.0121013
Piston speed = 0 m/s
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.99933, Final residual = 4.04263e11, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.928082, Final residual = 1.83196e06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.760023, Final residual = 1.22382e13, No Iterations 4
DILUPBiCG: Solving for h, Initial residual = 0.864411, Final residual = 2.38382e07, No Iterations 3
DICPCG: Solving for p, Initial residual = 1.28529e08, Final residual = 2.50118e24, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0156733, global = 0.0142217, cumulative = 0.0283917
DILUPBiCG: Solving for h, Initial residual = 0.100925, Final residual = 2.75858e06, No Iterations 3
DICPCG: Solving for p, Initial residual = 1.20648e09, Final residual = 1.61411e26, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0152062, global = 0.0140404, cumulative = 0.0424321
bounding epsilon, min: 7.72405e20 max: 1.22189e+58 average: 1.17852e+53
DILUPBiCG: Solving for epsilon, Initial residual = 0.319562, Final residual = 4.60129e22, No Iterations 1
bounding epsilon, min: 2.40945e+69 max: 1.20735e+71 average: 4.0323e+66
DILUPBiCG: Solving for k, Initial residual = 0.959559, Final residual = 5.07289e08, No Iterations 3
bounding k, min: 1.51071e+38 max: 1.7841e+39 average: 5.51239e+34
Mean pressure:5.93685e+25
Mean temperature:1.7358e+06
Mean u':7.90401e+13
ExecutionTime = 1381.48 s ClockTime = 1385 s
Courant Number mean: 0.0002019 max: 136.664 velocity magnitude: 5.30945e+19
deltaT = 7.401e25
MassFlows: presin = 0.00795114
> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 692
Increased the timePrecision from 14 to 15 to distinguish between timeNames at time 0.00129011
Crank angle = 23.2219 CAdeg
deltaZ = 0
DICPCG: Solving for cellMotionUz, Initial residual = 9.58306e06, Final residual = 9.58306e06, No Iterations 0
clearance: 0.0121013
Piston speed = 0 m/s
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.967025, Final residual = 3.49836e06, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.982864, Final residual = 1.18972e06, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.993804, Final residual = 6.60125e08, No Iterations 4
DILUPBiCG: Solving for h, Initial residual = 0.933219, Final residual = 0.000104695, No Iterations 3
DICPCG: Solving for p, Initial residual = 1.53971e11, Final residual = 1.53971e11, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.24475e+12, global = 2.99164e+10, cumulative = 2.99164e+10
DILUPBiCG: Solving for h, Initial residual = 0.78095, Final residual = 7.68682e05, No Iterations 3
DICPCG: Solving for p, Initial residual = 1.69201e11, Final residual = 1.69201e11, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.36544e+12, global = 3.96878e+11, cumulative = 3.66962e+11
bounding epsilon, min: 8.72817e20 max: 1.20735e+71 average: 4.59937e+66
DILUPBiCG: Solving for epsilon, Initial residual = 0.996955, Final residual = 352777, No Iterations 1000
bounding epsilon, min: 1.35106e+110 max: 1.92799e+112 average: 2.02929e+107
DILUPBiCG: Solving for k, Initial residual = 4.78121e13, Final residual = 4.78121e13, No Iterations 0
Mean pressure:5.93685e+25
Mean temperature:1.96123e+08
Mean u':7.90401e+13
ExecutionTime = 1395.93 s ClockTime = 1400 s
1) I think bounding k and epsilon says that the convergence has a problem. Any idea for overcoming ?
2) continuity error gets high value.
3)piston speed gets 0 m/s Why?
Can anybody help me for getting convergence?
I appreciate any help from you.
Thanks and best regards,
Sasan.
