CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Bounding k and epsilon & coldEngineFoam (http://www.cfd-online.com/Forums/openfoam/123064-bounding-k-epsilon-coldenginefoam.html)

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.24885e-05 max: 29.1764 velocity magnitude: 2.61013e+13
deltaT = 2.86066e-17
 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 CA-deg
deltaZ = -2.22045e-16
DICPCG:  Solving for cellMotionUz, Initial residual = 9.58306e-06, Final residual = 9.58306e-06, 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.06381e-06, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.999997, Final residual = 1.14924e-09, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.997449, Final residual = 1.51092e-08, No Iterations 3
DILUPBiCG:  Solving for h, Initial residual = 0.334811, Final residual = 4.57511e-08, No Iterations 2
DICPCG:  Solving for p, Initial residual = 1.10574e-05, Final residual = 2.18399e-21, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.27216e-07, global = -9.61115e-09, cumulative = -2.47977e-07
DILUPBiCG:  Solving for h, Initial residual = 0.123253, Final residual = 5.05499e-08, No Iterations 2
DICPCG:  Solving for p, Initial residual = 2.87271e-05, Final residual = 1.66698e-19, 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.4933e-20 max: 1.62696e+51 average: 4.68755e+46
DILUPBiCG:  Solving for epsilon, Initial residual = 0.000725064, Final residual = 1.36341e-20, 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.29296e-08, 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.26432e-22
 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 CA-deg
deltaZ = 0
DICPCG:  Solving for cellMotionUz, Initial residual = 9.58306e-06, Final residual = 9.58306e-06, 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.04263e-11, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.928082, Final residual = 1.83196e-06, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.760023, Final residual = 1.22382e-13, No Iterations 4
DILUPBiCG:  Solving for h, Initial residual = 0.864411, Final residual = 2.38382e-07, No Iterations 3
DICPCG:  Solving for p, Initial residual = 1.28529e-08, Final residual = 2.50118e-24, 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.75858e-06, No Iterations 3
DICPCG:  Solving for p, Initial residual = 1.20648e-09, Final residual = 1.61411e-26, 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.72405e-20 max: 1.22189e+58 average: 1.17852e+53
DILUPBiCG:  Solving for epsilon, Initial residual = 0.319562, Final residual = 4.60129e-22, 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.07289e-08, 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.401e-25
 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 CA-deg
deltaZ = 0
DICPCG:  Solving for cellMotionUz, Initial residual = 9.58306e-06, Final residual = 9.58306e-06, 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.49836e-06, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.982864, Final residual = 1.18972e-06, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.993804, Final residual = 6.60125e-08, 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.53971e-11, Final residual = 1.53971e-11, 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.68682e-05, No Iterations 3
DICPCG:  Solving for p, Initial residual = 1.69201e-11, Final residual = 1.69201e-11, 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.72817e-20 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.78121e-13, Final residual = 4.78121e-13, 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.

liwei September 13, 2013 12:12

hi, do you set all the boundary conditions properly?
we can contact with each other. My email address is liwei1@tju.edu.cn


All times are GMT -4. The time now is 02:37.