Heroic |
December 13, 2012 02:44 |
1 Attachment(s)
It seems I got things up and running now, including the k-epsilon model. Due to some high velocities through small cells I seem to be getting very high (bounded) values for epsilon locally. These are of the size 1e5 and for some reason they do not appear for every timestep, but only every now and then. Even more weird is that these bounded epsilon also show a decreasing behavior from time to time and k is rarely bounded.
Upon checking the results in paraview, the high values are always there, even if there is no reported bounding epsilon in the log file. I have terminated my simulation for now and I am going to discuss this behavior with a colleague of mine.
It would be nice to hear some opinions on this from others with more knowledge on this subject as well. Please see the attached picture and fragment of my log file. Keeping in mind that I am actually reaching good convergence, would it be wise to continue on this simulation or is it better to start over and find a way to totally eliminate the bounded variables? Is there a way to keep going with the current simulation and eliminate this behavior?
Code:
Courant Number mean: 9.55372e-05 max: 0.219991
Time = 0.13315
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 6.73373e-07, Final residual = 3.20221e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.22739e-06, Final residual = 1.4288e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1.30619e-05, Final residual = 6.47912e-10, No Iterations 1
DILUPBiCG: Solving for alpha1, Initial residual = 8.6486e-08, Final residual = 8.6486e-08, No Iterations 0
Solid phase fraction = 4.71282e-06 Min(alpha1) = 0 Max(alpha1) = 0.4
GAMG: Solving for p_rgh, Initial residual = 0.000202456, Final residual = 1.80145e-07, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.70531e-05, Final residual = 2.34654e-08, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.42526e-14, global = -3.46809e-15, cumulative = -4.2928e-11
GAMG: Solving for p_rgh, Initial residual = 4.09015e-06, Final residual = 2.69331e-09, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 6.81026e-07, Final residual = 6.59748e-10, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.7721e-14, global = -4.55999e-15, cumulative = -4.29325e-11
GAMG: Solving for p_rgh, Initial residual = 1.44307e-07, Final residual = 1.03114e-10, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 4.03171e-08, Final residual = 3.46799e-11, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.59999e-14, global = -4.04968e-15, cumulative = -4.29366e-11
GAMG: Solving for p_rgh, Initial residual = 1.17086e-08, Final residual = 6.23799e-12, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 9.4489e-09, Final residual = 9.74704e-13, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.43245e-14, global = -4.13848e-15, cumulative = -4.29407e-11
DILUPBiCG: Solving for epsilon, Initial residual = 2.65759e-07, Final residual = 1.47193e-11, No Iterations 1
bounding epsilon, min: -0.00553507 max: 103118 average: 39.5082
DILUPBiCG: Solving for k, Initial residual = 1.35119e-06, Final residual = 5.30976e-10, No Iterations 1
ExecutionTime = 63036.2 s ClockTime = 64522 s
... Some timesteps without bounding epsilon follow. I had to cut them off to stay within the max character limit ...
Courant Number mean: 9.55379e-05 max: 0.220011
Time = 0.13325
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 6.73169e-07, Final residual = 3.24129e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.22689e-06, Final residual = 1.43334e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1.30599e-05, Final residual = 6.52675e-10, No Iterations 1
DILUPBiCG: Solving for alpha1, Initial residual = 8.64956e-08, Final residual = 8.64956e-08, No Iterations 0
Solid phase fraction = 4.71282e-06 Min(alpha1) = 0 Max(alpha1) = 0.4
GAMG: Solving for p_rgh, Initial residual = 0.000201444, Final residual = 1.60844e-07, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.61183e-05, Final residual = 9.24127e-09, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 5.61681e-14, global = -2.98193e-15, cumulative = -4.30777e-11
GAMG: Solving for p_rgh, Initial residual = 3.86619e-06, Final residual = 2.34723e-09, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 6.35014e-07, Final residual = 4.72341e-10, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.67156e-14, global = -3.68224e-15, cumulative = -4.30814e-11
GAMG: Solving for p_rgh, Initial residual = 1.34363e-07, Final residual = 8.48447e-11, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 3.58766e-08, Final residual = 2.88988e-11, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.59304e-14, global = -3.55608e-15, cumulative = -4.30849e-11
GAMG: Solving for p_rgh, Initial residual = 1.06668e-08, Final residual = 5.16055e-12, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 8.27392e-09, Final residual = 8.58417e-13, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.31626e-14, global = -3.55087e-15, cumulative = -4.30885e-11
DILUPBiCG: Solving for epsilon, Initial residual = 2.65268e-07, Final residual = 1.45611e-11, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1.3485e-06, Final residual = 5.63431e-10, No Iterations 1
ExecutionTime = 63230.5 s ClockTime = 64720 s
Courant Number mean: 9.5538e-05 max: 0.220013
Time = 0.13326
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 6.73155e-07, Final residual = 3.24443e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.22686e-06, Final residual = 1.43377e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1.30598e-05, Final residual = 6.53214e-10, No Iterations 1
DILUPBiCG: Solving for alpha1, Initial residual = 8.64965e-08, Final residual = 8.64965e-08, No Iterations 0
Solid phase fraction = 4.71282e-06 Min(alpha1) = 0 Max(alpha1) = 0.4
GAMG: Solving for p_rgh, Initial residual = 0.000206822, Final residual = 1.82456e-07, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.9564e-05, Final residual = 1.68898e-08, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.30398e-14, global = -3.45633e-15, cumulative = -4.30919e-11
GAMG: Solving for p_rgh, Initial residual = 4.34841e-06, Final residual = 3.49822e-09, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 7.17663e-07, Final residual = 4.03597e-10, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.85699e-14, global = -4.27903e-15, cumulative = -4.30962e-11
GAMG: Solving for p_rgh, Initial residual = 1.53179e-07, Final residual = 1.33318e-10, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 4.03591e-08, Final residual = 3.8995e-11, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.77545e-14, global = -3.73605e-15, cumulative = -4.31e-11
GAMG: Solving for p_rgh, Initial residual = 1.16859e-08, Final residual = 5.82896e-12, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 9.22654e-09, Final residual = 5.36028e-13, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.40814e-14, global = -3.61683e-15, cumulative = -4.31036e-11
DILUPBiCG: Solving for epsilon, Initial residual = 2.6522e-07, Final residual = 1.45454e-11, No Iterations 1
bounding epsilon, min: -0.0146594 max: 103133 average: 39.5529
DILUPBiCG: Solving for k, Initial residual = 1.34824e-06, Final residual = 5.6666e-10, No Iterations 1
ExecutionTime = 63252.2 s ClockTime = 64743 s
Courant Number mean: 9.55381e-05 max: 0.220015
Time = 0.13327
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 6.73124e-07, Final residual = 3.24828e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.22681e-06, Final residual = 1.43556e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1.30595e-05, Final residual = 6.53373e-10, No Iterations 1
DILUPBiCG: Solving for alpha1, Initial residual = 8.64974e-08, Final residual = 8.64974e-08, No Iterations 0
Solid phase fraction = 4.71282e-06 Min(alpha1) = 0 Max(alpha1) = 0.4
GAMG: Solving for p_rgh, Initial residual = 0.00020463, Final residual = 1.80884e-07, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.84468e-05, Final residual = 1.89107e-08, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.12578e-14, global = -3.17955e-15, cumulative = -4.31067e-11
GAMG: Solving for p_rgh, Initial residual = 4.26351e-06, Final residual = 3.66951e-09, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 7.33175e-07, Final residual = 5.91502e-10, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.61081e-14, global = -4.4349e-15, cumulative = -4.31112e-11
GAMG: Solving for p_rgh, Initial residual = 1.55993e-07, Final residual = 6.94e-11, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 4.29368e-08, Final residual = 3.61959e-11, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.54608e-14, global = -3.98236e-15, cumulative = -4.31152e-11
GAMG: Solving for p_rgh, Initial residual = 1.19962e-08, Final residual = 6.20707e-12, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 9.62159e-09, Final residual = 5.64534e-13, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.27451e-14, global = -3.89695e-15, cumulative = -4.31191e-11
DILUPBiCG: Solving for epsilon, Initial residual = 2.65171e-07, Final residual = 1.45298e-11, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1.34797e-06, Final residual = 5.70079e-10, No Iterations 1
ExecutionTime = 63272 s ClockTime = 64764 s
Courant Number mean: 9.55381e-05 max: 0.220017
Time = 0.13328
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 6.73101e-07, Final residual = 3.25277e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.22672e-06, Final residual = 1.43649e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1.30592e-05, Final residual = 6.53624e-10, No Iterations 1
DILUPBiCG: Solving for alpha1, Initial residual = 8.64984e-08, Final residual = 8.64984e-08, No Iterations 0
Solid phase fraction = 4.71282e-06 Min(alpha1) = 0 Max(alpha1) = 0.4
GAMG: Solving for p_rgh, Initial residual = 0.000199448, Final residual = 1.54254e-07, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.50483e-05, Final residual = 1.89519e-08, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.20065e-14, global = -3.52836e-15, cumulative = -4.31226e-11
GAMG: Solving for p_rgh, Initial residual = 3.64174e-06, Final residual = 2.16334e-09, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 6.05503e-07, Final residual = 2.69169e-10, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.72243e-14, global = -4.49821e-15, cumulative = -4.31271e-11
GAMG: Solving for p_rgh, Initial residual = 1.28108e-07, Final residual = 8.60518e-11, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 3.53406e-08, Final residual = 2.89977e-11, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.67016e-14, global = -4.02559e-15, cumulative = -4.31311e-11
GAMG: Solving for p_rgh, Initial residual = 1.03696e-08, Final residual = 6.77484e-12, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 9.30403e-09, Final residual = 5.6405e-13, No Iterations 9
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.33453e-14, global = -3.93433e-15, cumulative = -4.3135e-11
DILUPBiCG: Solving for epsilon, Initial residual = 2.65122e-07, Final residual = 1.45141e-11, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1.3477e-06, Final residual = 5.73409e-10, No Iterations 1
ExecutionTime = 63293.2 s ClockTime = 64786 s
Courant Number mean: 9.55382e-05 max: 0.220019
Time = 0.13329
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 6.73085e-07, Final residual = 3.25652e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.22669e-06, Final residual = 1.43745e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1.3059e-05, Final residual = 6.53576e-10, No Iterations 1
DILUPBiCG: Solving for alpha1, Initial residual = 8.64994e-08, Final residual = 8.64994e-08, No Iterations 0
Solid phase fraction = 4.71282e-06 Min(alpha1) = 0 Max(alpha1) = 0.4
GAMG: Solving for p_rgh, Initial residual = 0.000203005, Final residual = 1.79916e-07, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.68167e-05, Final residual = 2.67205e-08, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.5267e-14, global = -3.79928e-15, cumulative = -4.31388e-11
GAMG: Solving for p_rgh, Initial residual = 3.97084e-06, Final residual = 2.6292e-09, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 6.54607e-07, Final residual = 6.21318e-10, No Iterations 6
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.91615e-14, global = -4.95388e-15, cumulative = -4.31438e-11
GAMG: Solving for p_rgh, Initial residual = 1.37601e-07, Final residual = 9.80408e-11, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 3.82828e-08, Final residual = 2.98687e-11, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.58331e-14, global = -4.6095e-15, cumulative = -4.31484e-11
GAMG: Solving for p_rgh, Initial residual = 1.08881e-08, Final residual = 7.40095e-12, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 1.00421e-08, Final residual = 8.25128e-13, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.31451e-14, global = -4.21148e-15, cumulative = -4.31526e-11
DILUPBiCG: Solving for epsilon, Initial residual = 2.65074e-07, Final residual = 1.44984e-11, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1.34743e-06, Final residual = 5.76752e-10, No Iterations 1
ExecutionTime = 63313.1 s ClockTime = 64807 s
Courant Number mean: 9.55383e-05 max: 0.220021
Time = 0.1333
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 6.73067e-07, Final residual = 3.25873e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 2.22662e-06, Final residual = 1.4377e-10, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1.30587e-05, Final residual = 6.53616e-10, No Iterations 1
DILUPBiCG: Solving for alpha1, Initial residual = 8.65003e-08, Final residual = 8.65003e-08, No Iterations 0
Solid phase fraction = 4.71282e-06 Min(alpha1) = 0 Max(alpha1) = 0.4
GAMG: Solving for p_rgh, Initial residual = 0.000203332, Final residual = 1.55745e-07, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 2.68662e-05, Final residual = 2.01764e-08, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.53667e-14, global = -4.25321e-15, cumulative = -4.31569e-11
GAMG: Solving for p_rgh, Initial residual = 3.94813e-06, Final residual = 2.79355e-09, No Iterations 5
GAMG: Solving for p_rgh, Initial residual = 6.43015e-07, Final residual = 2.92496e-10, No Iterations 7
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.86857e-14, global = -4.21542e-15, cumulative = -4.31611e-11
GAMG: Solving for p_rgh, Initial residual = 1.36678e-07, Final residual = 9.05711e-11, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 3.68273e-08, Final residual = 2.84291e-11, No Iterations 5
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.78365e-14, global = -4.25555e-15, cumulative = -4.31653e-11
GAMG: Solving for p_rgh, Initial residual = 1.08069e-08, Final residual = 6.99844e-12, No Iterations 8
GAMG: Solving for p_rgh, Initial residual = 9.87503e-09, Final residual = 7.91717e-13, No Iterations 8
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.48573e-14, global = -4.24741e-15, cumulative = -4.31696e-11
DILUPBiCG: Solving for epsilon, Initial residual = 2.65025e-07, Final residual = 1.44828e-11, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1.34717e-06, Final residual = 5.80109e-10, No Iterations 1
ExecutionTime = 63335.1 s ClockTime = 64831 s
|