CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   time step continuity problem in VAWT simulation (https://www.cfd-online.com/Forums/openfoam-solving/160905-time-step-continuity-problem-vawt-simulation.html)

lpz_michele October 15, 2015 06:13
time step continuity problem in VAWT simulation
 
Hi everyone! I'm trying to simulate the flowfield on a VAWT Darrieus Turbine. I generated 2 different mesh from Pointwise, one for the rotor, another for the volume in which i'm running the simulation.
Than i merged that with mergeMeshes. Set the rotation using Pointwise cell and dynamicMeshDict. I'm using pimpleDyMFoam with GAMG, the calculation starts, but the number of iteration for the first time step is huge and time step continuity errors gets bigger and bigger.

Note that:
  • pala*, is the airfoil of the wind turbine
  • cubo*, is just part of the volume in which the case is running
  • fisso*, volume* are cyclicAMI
Here are all the attachments:

lpz_michele October 15, 2015 06:14
log file

Code:
   
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 2.49637e-14, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 2.85471e-14, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 2.78977e-14, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00808571, No Iterations 13
time step continuity errors : sum local = 4.41e-11, global = -5.01669e-12, cumulative = -5.01669e-12
GAMG:  Solving for p, Initial residual = 0.166509, Final residual = 0.000825118, No Iterations 7
time step continuity errors : sum local = 2.09824e-10, global = -4.30853e-11, cumulative = -4.8102e-11
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.616245, Final residual = 1.41262e-05, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.129979, Final residual = 4.87507e-06, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.130425, Final residual = 3.96162e-06, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.132593, Final residual = 0.000787422, No Iterations 7
time step continuity errors : sum local = 2.12629e-10, global = -4.30272e-11, cumulative = -9.11292e-11
GAMG:  Solving for p, Initial residual = 0.115701, Final residual = 0.000790563, No Iterations 5
time step continuity errors : sum local = 3.43266e-10, global = -5.17038e-11, cumulative = -1.42833e-10
PIMPLE: iteration 3
smoothSolver:  Solving for Ux, Initial residual = 0.0630906, Final residual = 0.00458348, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0110506, Final residual = 0.000677898, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.00817843, Final residual = 0.000674459, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.660562, Final residual = 0.00658876, No Iterations 64
time step continuity errors : sum local = 7.75362e-09, global = -3.95022e-12, cumulative = -1.46783e-10
GAMG:  Solving for p, Initial residual = 0.34499, Final residual = 0.00343841, No Iterations 70
time step continuity errors : sum local = 1.98197e-08, global = -5.96419e-11, cumulative = -2.06425e-10
PIMPLE: iteration 4
smoothSolver:  Solving for Ux, Initial residual = 0.982638, Final residual = 0.0722246, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.988885, Final residual = 0.0726794, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.983183, Final residual = 0.072299, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.998999, Final residual = 0.00962466, No Iterations 154
time step continuity errors : sum local = 2.82917e-05, global = -2.73187e-09, cumulative = -2.9383e-09
GAMG:  Solving for p, Initial residual = 0.000466197, Final residual = 4.64261e-06, No Iterations 207
time step continuity errors : sum local = 0.000460808, global = -1.49189e-07, cumulative = -1.52128e-07
PIMPLE: iteration 5
smoothSolver:  Solving for Ux, Initial residual = 0.998628, Final residual = 0.0978001, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 0.998121, Final residual = 0.0977083, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.998746, Final residual = 0.0978102, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.0698349, Final residual = 0.000631727, No Iterations 42
time step continuity errors : sum local = 0.279906, global = 0.000685844, cumulative = 0.000685692
GAMG:  Solving for p, Initial residual = 4.2465e-08, Final residual = 4.2465e-08, No Iterations 0
time step continuity errors : sum local = 73.2725, global = -0.000797158, cumulative = -0.000111465
PIMPLE: iteration 6
smoothSolver:  Solving for Ux, Initial residual = 0.997507, Final residual = 0.0408674, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.997316, Final residual = 0.0408327, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.997554, Final residual = 0.0408676, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00142366, Final residual = 1.39848e-05, No Iterations 364
time step continuity errors : sum local = 1086.86, global = -42.1764, cumulative = -42.1765
GAMG:  Solving for p, Initial residual = 4.33906e-11, Final residual = 4.33906e-11, No Iterations 0
time step continuity errors : sum local = 368842, global = -15.6452, cumulative = -57.8217
PIMPLE: iteration 7
smoothSolver:  Solving for Ux, Initial residual = 0.98951, Final residual = 0.00299387, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.981899, Final residual = 0.00318746, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.987728, Final residual = 0.0030354, No Iterations 2
GAMG:  Solving for p, Initial residual = 2.87655e-07, Final residual = 8.20569e-08, No Iterations 2
time step continuity errors : sum local = 2.94219e+06, global = 50046.2, cumulative = 49988.4
GAMG:  Solving for p, Initial residual = 3.15004e-13, Final residual = 3.15004e-13, No Iterations 0
time step continuity errors : sum local = 3.81806e+07, global = 16551.4, cumulative = 66539.8
PIMPLE: iteration 8
smoothSolver:  Solving for Ux, Initial residual = 0.99324, Final residual = 0.0202374, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.984123, Final residual = 0.0197882, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.991156, Final residual = 0.0201318, No Iterations 2
GAMG:  Solving for p, Initial residual = 1.223e-07, Final residual = 7.29218e-08, No Iterations 2
time step continuity errors : sum local = 2.09341e+09, global = 2.06686e+07, cumulative = 2.07351e+07
GAMG:  Solving for p, Initial residual = 3.83492e-13, Final residual = 3.83492e-13, No Iterations 0
time step continuity errors : sum local = 2.31154e+10, global = 5.76172e+06, cumulative = 2.64969e+07
PIMPLE: iteration 9
smoothSolver:  Solving for Ux, Initial residual = 0.998567, Final residual = 0.0721183, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.997347, Final residual = 0.0725961, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.998323, Final residual = 0.0721812, No Iterations 2
GAMG:  Solving for p, Initial residual = 5.87591e-08, Final residual = 5.87591e-08, No Iterations 0
time step continuity errors : sum local = 5.03518e+12, global = -851499, cumulative = 2.56454e+07
GAMG:  Solving for p, Initial residual = 6.40719e-08, Final residual = 6.40719e-08, No Iterations 0
time step continuity errors : sum local = 5.49045e+12, global = -852827, cumulative = 2.47925e+07
PIMPLE: iteration 10
smoothSolver:  Solving for Ux, Initial residual = 0.991257, Final residual = 0.00126755, No Iterations 4
smoothSolver:  Solving for Uy, Initial residual = 0.991242, Final residual = 0.0012675, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.991255, Final residual = 0.00126755, No Iterations 4
GAMG:  Solving for p, Initial residual = 2.83486e-05, Final residual = 2.78443e-07, No Iterations 476
time step continuity errors : sum local = 7.89681e+13, global = 7.90082e+09, cumulative = 7.92561e+09
GAMG:  Solving for p, Initial residual = 7.28708e-14, Final residual = 7.28708e-14, No Iterations 0
time step continuity errors : sum local = 1.69118e+17, global = 1.67098e+12, cumulative = 1.67891e+12
PIMPLE: iteration 11
smoothSolver:  Solving for Ux, Initial residual = 0.994182, Final residual = 0.0014273, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.99709, Final residual = 0.00132306, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.997727, Final residual = 0.00131261, No Iterations 2
GAMG:  Solving for p, Initial residual = 1.98488e-05, Final residual = 1.98321e-07, No Iterations 759
time step continuity errors : sum local = 2.34008e+16, global = 5.8166e+12, cumulative = 7.49551e+12
GAMG:  Solving for p, Initial residual = 2.26884e-12, Final residual = 2.26884e-12, No Iterations 0
time step continuity errors : sum local = 1.06405e+19, global = -7.31893e+13, cumulative = -6.56938e+13
PIMPLE: iteration 12
smoothSolver:  Solving for Ux, Initial residual = 0.995118, Final residual = 0.0175633, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.993056, Final residual = 0.017419, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.994499, Final residual = 0.0175166, No Iterations 2
GAMG:  Solving for p, Initial residual = 4.89944e-07, Final residual = 7.80626e-08, No Iterations 8
time step continuity errors : sum local = 1.91622e+20, global = 7.60396e+16, cumulative = 7.59739e+16
GAMG:  Solving for p, Initial residual = 9.02203e-13, Final residual = 9.02203e-13, No Iterations 0
time step continuity errors : sum local = 1.71062e+22, global = 1.58506e+18, cumulative = 1.66103e+18
PIMPLE: iteration 13
smoothSolver:  Solving for Ux, Initial residual = 0.990097, Final residual = 0.0229477, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.987019, Final residual = 0.0225001, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.989434, Final residual = 0.022828, No Iterations 2
GAMG:  Solving for p, Initial residual = 7.82812e-08, Final residual = 7.82812e-08, No Iterations 0
time step continuity errors : sum local = 6.06662e+23, global = 1.39824e+18, cumulative = 3.05927e+18
GAMG:  Solving for p, Initial residual = 8.48612e-08, Final residual = 8.48612e-08, No Iterations 0
time step continuity errors : sum local = 6.57656e+23, global = 1.36088e+18, cumulative = 4.42015e+18
PIMPLE: iteration 14
smoothSolver:  Solving for Ux, Initial residual = 0.99537, Final residual = 0.0344139, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.99524, Final residual = 0.0344185, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.995337, Final residual = 0.0344148, No Iterations 2
GAMG:  Solving for p, Initial residual = 1.17968e-05, Final residual = 1.16841e-07, No Iterations 110
time step continuity errors : sum local = 2.17112e+24, global = 9.08021e+22, cumulative = 9.08065e+22
GAMG:  Solving for p, Initial residual = 1.01278e-13, Final residual = 1.01278e-13, No Iterations 0
time step continuity errors : sum local = 7.83079e+27, global = 3.95176e+22, cumulative = 1.30324e+23
PIMPLE: iteration 15
smoothSolver:  Solving for Ux, Initial residual = 0.999463, Final residual = 0.00546958, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.999025, Final residual = 0.00546167, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.999773, Final residual = 0.00547299, No Iterations 2
GAMG:  Solving for p, Initial residual = 1.00373e-06, Final residual = 9.77242e-08, No Iterations 23
time step continuity errors : sum local = 4.63235e+28, global = 3.52165e+27, cumulative = 3.52178e+27
GAMG:  Solving for p, Initial residual = 2.59211e-13, Final residual = 2.59211e-13, No Iterations 0
time step continuity errors : sum local = 4.76462e+30, global = 1.2969e+27, cumulative = 4.81868e+27
PIMPLE: iteration 16
smoothSolver:  Solving for Ux, Initial residual = 0.994007, Final residual = 0.0166685, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.992732, Final residual = 0.0166767, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.99359, Final residual = 0.0166783, No Iterations 2
GAMG:  Solving for p, Initial residual = 2.62176e-06, Final residual = 9.86938e-08, No Iterations 21
time step continuity errors : sum local = 2.04098e+31, global = -3.87112e+30, cumulative = -3.8663e+30
GAMG:  Solving for p, Initial residual = 9.45924e-13, Final residual = 9.45924e-13, No Iterations 0
time step continuity errors : sum local = 4.57275e+33, global = -2.73623e+30, cumulative = -6.60252e+30

lpz_michele October 15, 2015 06:15
Code:
PIMPLE: iteration 17
smoothSolver:  Solving for Ux, Initial residual = 0.985753, Final residual = 0.0394114, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.984547, Final residual = 0.0391394, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.985481, Final residual = 0.0392419, No Iterations 1
GAMG:  Solving for p, Initial residual = 8.73274e-07, Final residual = 9.91197e-08, No Iterations 25
time step continuity errors : sum local = 1.92743e+34, global = 1.00819e+33, cumulative = 1.00159e+33
GAMG:  Solving for p, Initial residual = 9.71117e-12, Final residual = 9.71117e-12, No Iterations 0
time step continuity errors : sum local = 1.16766e+36, global = 3.02168e+32, cumulative = 1.30376e+33
PIMPLE: iteration 18
smoothSolver:  Solving for Ux, Initial residual = 0.994286, Final residual = 0.0377732, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.993101, Final residual = 0.0381006, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.993921, Final residual = 0.0377407, No Iterations 2
GAMG:  Solving for p, Initial residual = 4.75613e-07, Final residual = 9.60022e-08, No Iterations 7
time step continuity errors : sum local = 1.23822e+38, global = -2.28473e+36, cumulative = -2.28343e+36
GAMG:  Solving for p, Initial residual = 4.29534e-11, Final residual = 4.29534e-11, No Iterations 0
time step continuity errors : sum local = 1.50703e+39, global = -7.85489e+35, cumulative = -3.06891e+36
PIMPLE: iteration 19
smoothSolver:  Solving for Ux, Initial residual = 0.993611, Final residual = 0.00693789, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.988574, Final residual = 0.0311961, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.993369, Final residual = 0.00466954, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.000346911, Final residual = 3.40028e-06, No Iterations 35
time step continuity errors : sum local = 8.90928e+39, global = 3.20522e+38, cumulative = 3.17453e+38
GAMG:  Solving for p, Initial residual = 3.77969e-11, Final residual = 3.77969e-11, No Iterations 0
time step continuity errors : sum local = 1.64211e+42, global = 1.01557e+38, cumulative = 4.1901e+38
PIMPLE: iteration 20
smoothSolver:  Solving for Ux, Initial residual = 0.991959, Final residual = 0.0579415, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.991925, Final residual = 0.087055, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.990317, Final residual = 0.0502125, No Iterations 2
GAMG:  Solving for p, Initial residual = 1.72486e-05, Final residual = 1.67883e-07, No Iterations 23
time step continuity errors : sum local = 3.53028e+43, global = 9.35766e+41, cumulative = 9.36185e+41
GAMG:  Solving for p, Initial residual = 8.72465e-12, Final residual = 8.72465e-12, No Iterations 0
time step continuity errors : sum local = 5.19123e+45, global = 3.43068e+41, cumulative = 1.27925e+42
PIMPLE: iteration 21
smoothSolver:  Solving for Ux, Initial residual = 0.985264, Final residual = 0.0592161, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.985617, Final residual = 0.0595835, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.985759, Final residual = 0.0590704, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.000486053, Final residual = 4.69196e-06, No Iterations 22
time step continuity errors : sum local = 1.0733e+47, global = -3.43595e+45, cumulative = -3.43467e+45
GAMG:  Solving for p, Initial residual = 4.96428e-08, Final residual = 4.96428e-08, No Iterations 0
time step continuity errors : sum local = 2.213e+49, global = -1.70873e+45, cumulative = -5.1434e+45
PIMPLE: iteration 22
smoothSolver:  Solving for Ux, Initial residual = 0.995176, Final residual = 0.0486652, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.995195, Final residual = 0.0486729, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.995156, Final residual = 0.0486722, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.682619, Final residual = 0.00681379, No Iterations 51
time step continuity errors : sum local = 8.17385e+50, global = -5.25098e+47, cumulative = -5.30241e+47
GAMG:  Solving for p, Initial residual = 1.08342e-08, Final residual = 1.08342e-08, No Iterations 0
time step continuity errors : sum local = 1.32546e+53, global = 2.03494e+48, cumulative = 1.5047e+48
PIMPLE: iteration 23
smoothSolver:  Solving for Ux, Initial residual = 0.984067, Final residual = 0.0829088, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.984554, Final residual = 0.0829242, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.9846, Final residual = 0.0829219, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.398936, Final residual = 0.00369698, No Iterations 25
time step continuity errors : sum local = 2.19833e+54, global = 3.12079e+52, cumulative = 3.12094e+52
GAMG:  Solving for p, Initial residual = 9.1991e-08, Final residual = 9.1991e-08, No Iterations 0
time step continuity errors : sum local = 3.75852e+56, global = 3.88852e+52, cumulative = 7.00946e+52
PIMPLE: iteration 24
smoothSolver:  Solving for Ux, Initial residual = 0.994046, Final residual = 0.090959, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.994048, Final residual = 0.0909252, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.993961, Final residual = 0.0909041, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.848587, Final residual = 0.00769635, No Iterations 22
time step continuity errors : sum local = 1.10779e+58, global = 1.58049e+56, cumulative = 1.58119e+56
GAMG:  Solving for p, Initial residual = 2.41036e-08, Final residual = 2.41036e-08, No Iterations 0
time step continuity errors : sum local = 1.78147e+60, global = -8.98121e+54, cumulative = 1.49138e+56
PIMPLE: iteration 25
smoothSolver:  Solving for Ux, Initial residual = 0.985409, Final residual = 0.0976268, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.985872, Final residual = 0.0976392, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.985929, Final residual = 0.0976339, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.34435, Final residual = 0.00318646, No Iterations 22
time step continuity errors : sum local = 3.72033e+61, global = 5.45526e+59, cumulative = 5.45675e+59
GAMG:  Solving for p, Initial residual = 2.17059e-08, Final residual = 2.17059e-08, No Iterations 0
time step continuity errors : sum local = 6.97103e+63, global = -5.90904e+58, cumulative = 4.86585e+59
PIMPLE: iteration 26
smoothSolver:  Solving for Ux, Initial residual = 0.995057, Final residual = 0.0272908, No Iterations 3
smoothSolver:  Solving for Uy, Initial residual = 0.99507, Final residual = 0.0273262, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 0.995016, Final residual = 0.0273368, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.61391, Final residual = 0.00608321, No Iterations 68
time step continuity errors : sum local = 3.72131e+65, global = 2.00756e+62, cumulative = 2.01242e+62
GAMG:  Solving for p, Initial residual = 4.74907e-08, Final residual = 4.74907e-08, No Iterations 0
time step continuity errors : sum local = 8.03038e+67, global = 5.25702e+63, cumulative = 5.45826e+63
PIMPLE: iteration 27
smoothSolver:  Solving for Ux, Initial residual = 0.97703, Final residual = 0.0743789, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.977785, Final residual = 0.0743152, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.977813, Final residual = 0.0742994, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.756536, Final residual = 0.00725901, No Iterations 37
time step continuity errors : sum local = 8.15355e+68, global = -1.12884e+67, cumulative = -1.1283e+67
GAMG:  Solving for p, Initial residual = 1.09514e-07, Final residual = 2.08685e-08, No Iterations 1
time step continuity errors : sum local = 2.46675e+70, global = 2.2231e+67, cumulative = 1.0948e+67
PIMPLE: iteration 28
smoothSolver:  Solving for Ux, Initial residual = 0.99716, Final residual = 0.0768253, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.997224, Final residual = 0.0767676, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.99719, Final residual = 0.076746, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.557303, Final residual = 0.00536858, No Iterations 30
time step continuity errors : sum local = 7.4905e+72, global = 1.56968e+70, cumulative = 1.57078e+70
GAMG:  Solving for p, Initial residual = 6.42801e-08, Final residual = 6.42801e-08, No Iterations 0
time step continuity errors : sum local = 9.6001e+74, global = -1.31788e+71, cumulative = -1.16081e+71
PIMPLE: iteration 29
smoothSolver:  Solving for Ux, Initial residual = 0.986285, Final residual = 0.0673124, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.986847, Final residual = 0.0672444, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.986977, Final residual = 0.0672249, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.695817, Final residual = 0.00638954, No Iterations 15
time step continuity errors : sum local = 8.18519e+75, global = -4.5645e+74, cumulative = -4.56566e+74
GAMG:  Solving for p, Initial residual = 1.84291e-07, Final residual = 3.74045e-08, No Iterations 1
time step continuity errors : sum local = 2.23164e+77, global = -1.55497e+74, cumulative = -6.12064e+74
PIMPLE: iteration 30
smoothSolver:  Solving for Ux, Initial residual = 0.996761, Final residual = 0.0606023, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.99677, Final residual = 0.0612115, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.996829, Final residual = 0.0611781, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.613433, Final residual = 0.00607106, No Iterations 53
time step continuity errors : sum local = 4.48328e+79, global = -1.18692e+78, cumulative = -1.18753e+78
GAMG:  Solving for p, Initial residual = 5.84873e-08, Final residual = 5.84873e-08, No Iterations 0
time step continuity errors : sum local = 6.19905e+81, global = -9.97691e+77, cumulative = -2.18522e+78
PIMPLE: iteration 31
smoothSolver:  Solving for Ux, Initial residual = 0.952807, Final residual = 0.058891, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.954521, Final residual = 0.0592418, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.954898, Final residual = 0.0593542, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.371946, Final residual = 0.00360844, No Iterations 20
time step continuity errors : sum local = 2.7809e+82, global = -4.09669e+80, cumulative = -4.11854e+80
GAMG:  Solving for p, Initial residual = 1.74061e-07, Final residual = 4.09465e-08, No Iterations 1
time step continuity errors : sum local = 8.05945e+83, global = 6.17798e+79, cumulative = -3.50074e+80
PIMPLE: iteration 32
smoothSolver:  Solving for Ux, Initial residual = 0.9926, Final residual = 0.0569642, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.992885, Final residual = 0.0577423, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.993074, Final residual = 0.0579465, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.0399935, Final residual = 0.000381604, No Iterations 32
time step continuity errors : sum local = 5.67533e+85, global = -5.85781e+84, cumulative = -5.85816e+84
GAMG:  Solving for p, Initial residual = 1.16347e-08, Final residual = 1.16347e-08, No Iterations 0
time step continuity errors : sum local = 7.63996e+87, global = -2.45486e+84, cumulative = -8.31302e+84
PIMPLE: iteration 33
smoothSolver:  Solving for Ux, Initial residual = 0.989217, Final residual = 0.0461757, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.985905, Final residual = 0.0457226, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.984847, Final residual = 0.0456371, No Iterations 2

lpz_michele October 15, 2015 06:17
controlDict
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application    pimpleDyMFoam;

startFrom      startTime;

startTime      0;

stopAt          endTime;

endTime        0.7853976;

deltaT          0.0000109083;

writeControl    adjustableRunTime;

writeInterval  0.0000109083;

purgeWrite      0;

writeFormat    ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision  6;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo          1;

// ************************************************************************* //

functions
(
forces
{
type forces;
functionObjectLibs ("libforces.so");
outputControl timeStep;
outputInterval 1;
patches (pala1 pala1sopra pala1sotto pala2 pala2sopra pala2sotto pala3 pala3sopra pala3sotto);
// pname p;
// Uname U;
rhoName rhoInf;
log true;
rhoInf 1.205;
CofR (0 0 0);
}
);
fvSchemes
Code:
  /*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default        Euler;
}

gradSchemes
{
    default                    leastSquares;
    grad(p)                    leastSquares;
    grad(U)                    leastSquares;

}

divSchemes
{
    default                    none;

    div(phi,U)                  Gauss linearUpwind grad(U);
    div(phi,k)                  Gauss upwind;
    div(phi,omega)            Gauss upwind;
    div((nuEff*dev(T(grad(U)))))    Gauss linear;
}

laplacianSchemes
{
    default                    Gauss linear corrected;
}

interpolationSchemes
{
    default                    linear;
    interpolate(HbyA)            linear;

}

snGradSchemes
{
    default                    corrected;
}

fluxRequired
{
    default        no;
    p              ;
    pcorr          ;
}


// ************************************************************************* //
fvSolution
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    pcorr
    {
        solver GAMG;
    smoother GaussSeidel;
    nPreSweeps 0;
    nPostSweeps 2;
    cacheAgglomeration off;
    agglomerator faceAreaPair;
    nCellsInCoarsestLevel 10;
    mergeLevels 1;

    tolerance 0.02;
    relTol 0;
    }

    p
    {
    $pcorr;
    tolerance 1e-07;
    relTol 0.01;
    } 
   
    pFinal
    {
        $p;
    tolerance    1e-07;
        relTol          0;
    }


    U
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
   
    UFinal
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
   
    "(k|omega)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance      1e-7;
        relTol          0.1;
    }

    "(k|omega)Final"
    {
        $U;
    tolerance    1e-7;
        relTol          0;
    }

    cellMotionUx
    {
    solver        PCG;
    preconditioner    DIC;
    tolerance    1e-08;
    relTol        0;
    }

}

PIMPLE
{
    correctPhi        no;
    nOuterCorrectors    100;
    nCorrectors        2;
    nNonOrthogonalCorrectors 0;
   
    pRefCell        0;
    pRefValue        0;

residualControl
  {
    "(p|U|k|omega)"
    {
      tolerance 1e-2;
      relTol 0;
      }
  }
}

relaxationFactors
{
    fields
    {
      "p.*"        0.3;
    }
    equations
    {
   
    "U.*"        0.7;
    "k.*"        0.8;
    "omega.*"        0.8;
    }
}


// ************************************************************************* //
boundary
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      polyBoundaryMesh;
    location    "0.000109083/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

21
(
    pala1
    {
        type            wall;
        inGroups        1(wall);
        nFaces          9506;
        startFace      5799010;
    }
    pala1sopra
    {
        type            wall;
        inGroups        1(wall);
        nFaces          301;
        startFace      5808516;
    }
    pala1sotto
    {
        type            wall;
        inGroups        1(wall);
        nFaces          301;
        startFace      5808817;
    }
    pala2
    {
        type            wall;
        inGroups        1(wall);
        nFaces          9506;
        startFace      5809118;
    }
    pala2sopra
    {
        type            wall;
        inGroups        1(wall);
        nFaces          301;
        startFace      5818624;
    }
    pala2sotto
    {
        type            wall;
        inGroups        1(wall);
        nFaces          301;
        startFace      5818925;
    }
    pala3
    {
        type            wall;
        inGroups        1(wall);
        nFaces          9506;
        startFace      5819226;
    }
    pala3sopra
    {
        type            wall;
        inGroups        1(wall);
        nFaces          301;
        startFace      5828732;
    }
    pala3sotto
    {
        type            wall;
        inGroups        1(wall);
        nFaces          301;
        startFace      5829033;
    }
    volumelato
    {
        type            cyclicAMI;
    inGroups    1(cyclicAMI);
    nFaces          1026;
        startFace      5829334;
    matchTolerance    0.0001;
    transform    noOrdering;
    neighbourPatch    fisso;
   
    }
    volumesopra
    {
        type            cyclicAMI;
    inGroups    1(cyclicAMI);
        nFaces          744;
        startFace      5830360;
    matchTolerance    0.0001;
    transform    noOrdering;
    neighbourPatch    fissosopra;
    }
    volumesotto
    {
        type            cyclicAMI;
    inGroups    1(cyclicAMI);
        nFaces          744;
        startFace      5831104;
    matchTolerance    0.0001;
    transform    noOrdering;
    neighbourPatch    fissosotto;
    }
    cubodx
    {
        type            patch;
        nFaces          592;
        startFace      5831848;
    }
    cubosopra
    {
        type            patch;
        nFaces          588;
        startFace      5832440;
    }
    cubosotto
    {
        type            patch;
        nFaces          588;
        startFace      5833028;
    }
    cubosx
    {
        type            patch;
        nFaces          590;
        startFace      5833616;
    }
    fisso
    {
        type            cyclicAMI;
    inGroups    1(cyclicAMI);
        nFaces          576;
        startFace      5834206;
    matchTolerance    0.0001;
    transform    noOrdering;
    neighbourPatch    volumelato;
    }
    fissosopra
    {
        type            cyclicAMI;
    inGroups    1(cyclicAMI);
        nFaces          180;
        startFace      5834782;
    matchTolerance    0.0001;
    transform    noOrdering;
    neighbourPatch    volumesopra;
    }
    fissosotto
    {
        type            cyclicAMI;
    inGroups    1(cyclicAMI);
        nFaces          180;
        startFace      5834962;
    matchTolerance    0.0001;
    transform    noOrdering;
    neighbourPatch    volumesotto;
    }
    inlet
    {
        type            patch;
        nFaces          478;
        startFace      5835142;
    }
    outlet
    {
        type            patch;
        nFaces          476;
        startFace      5835620;
    }
)

// ************************************************************************* //
dynamicMeshDict
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh    solidBodyMotionFvMesh ;

motionSolverLibs ("libfvMotionSolvers.so");

solidBodyMotionFvMeshCoeffs
{
  cellZone        rotor;
 
  solidBodyMotionFunction rotatingMotion;
  rotatingMotionCoeffs
  {
    origin (0 0 0);
    axis (0 0 1);
    omega 40;// rad/s
  }
}


// ************************************************************************* //
decomposeParDict
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      decomposeParDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

numberOfSubdomains    4;

preservePatches    (fisso fissosopra fissosotto volumelato volumesopra volumesotto);

method              scotch;

distributed        no;

roots            ();


// ************************************************************************* //

lpz_michele October 15, 2015 06:19
k
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [ 0 2 -2 0 0 0 0 ];

internalField  uniform 0.375;

boundaryField
{
 
 
    inlet
    {
        type                turbulentIntensityKineticEnergyInlet;
    intensity        0.05;
    value            $internalField;
    }
   
   
    outlet
    {
        type                inletOutlet;
        inletValue          $internalField;
        value              $internalField;     
    }
 
 
    pala1
    {
        type            kqRWallFunction;
        value              $internalField;
    }
   
   
    pala1sopra
    {
        type              kqRWallFunction;
        value              $internalField;     
    }
   
   
    pala1sotto
    {
    type            kqRWallFunction;
        value                  $internalField;
    }
   
   
    pala2
    {
        type              kqRWallFunction;
        value              $internalField; 
    }
   
   
    pala2sopra
    {
       
      type            kqRWallFunction;
      value                  $internalField;
    }
   
   
    pala2sotto
    {
      type              kqRWallFunction;
        value              $internalField;
    }
   
   
    pala3
    {
        type            kqRWallFunction;
        value              $internalField;
    }
   
   
    pala3sopra
    {
        type              kqRWallFunction;
        value              $internalField;   
    }
   
   
    pala3sotto
    {
    type            kqRWallFunction;
    value                  $internalField;
    }
   
      volumelato
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
   
    volumesopra
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
   
    volumesotto
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
    cubodx
    {
        type        zeroGradient;
       
    }
   
   
    cubosopra
    {
        type        zeroGradient;
       
    }
   
   
    cubosotto
    {
        type        zeroGradient;
     
    }
   
   
    cubosx
    {
        type        zeroGradient;
       
    }
   
   
    fisso
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
   
    fissosopra
    {
        type        cyclicAMI;
        value        $internalField;
    }
   
   
    fissosotto
    {
        type        cyclicAMI;
        value        $internalField;
    }
 
}


// ************************************************************************* //
omega
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //



dimensions      [0 0 -1 0 0 0 0];

internalField  uniform 1118.034;

boundaryField
{
    cubodx
    {
        type        zeroGradient;
       
    }
   
   
    cubosopra
    {
        type        zeroGradient;
       
    }
   
   
    cubosotto
    {
        type        zeroGradient;
     
    }
   
   
    cubosx
    {
        type        zeroGradient;
       
    }
   
    inlet
    {
        type            fixedValue;
        value          $internalField;   
    }
   
   
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value          $internalField;     
    }
   
   
    fisso
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
   
    fissosopra
    {
        type        cyclicAMI;
        value        $internalField;
    }
   
   
    fissosotto
    {
        type        cyclicAMI;
        value        $internalField;
    }
   
    volumelato
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
   
    volumesopra
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
   
    volumesotto
    {
        type        cyclicAMI;
    value        $internalField;
    }
   
      pala1
    {
        type        omegaWallFunction;
    value        $internalField;
    }
   
   
    pala1sopra
    {
        type        omegaWallFunction;
    value        $internalField;
    }
   
   
    pala1sotto
    {
    type        omegaWallFunction;
    value        $internalField;
    }
   
   
    pala2
    {
        type        omegaWallFunction;
    value        $internalField;
    }
   
   
    pala2sopra
    {
       
      type        omegaWallFunction;
      value        $internalField;
    }
   
   
    pala2sotto
    {
      type          omegaWallFunction;
      value        $internalField;
    }
   
   
    pala3
    {
        type        omegaWallFunction;
    value        $internalField;
    }
   
   
    pala3sopra
    {
        type        omegaWallFunction;
    value        $internalField;
    }
   
   
    pala3sotto
    {
    type        omegaWallFunction;
    value        $internalField;
    }
   
}


// ************************************************************************* //
U
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField  uniform (0 0 0);

boundaryField
{
    cubodx
    {
        type slip;       
    }
   
   
    cubosopra
    {
        type slip;       
    }
   
   
   
    cubosotto
    {
        type slip;     
    }
   
   
   
    cubosx
    {
        type slip;       
    }
   
   
   
    fisso
    {
        type cyclicAMI;
    value $internalField;
    }
   
   
   
    fissosopra
    {
        type cyclicAMI;
        value $internalField;
    }
   
   
   
    fissosotto
    {
        type cyclicAMI;
        value $internalField;
    }
   
   
   
    inlet
    {
        type fixedValue;
    value uniform (10.0 0 0);     
    }
   
   
   
    outlet
    {
        type zeroGradient;
       
    }
 
    pala1
    {
        type fixedValue;
    value $internalField;
    }
   
   
    pala1sopra
    {
        type fixedValue;
    value $internalField;
       
    }
   
   
    pala1sotto
    {
    type fixedValue;
    value $internalField;
    }
   
   
   
    pala2
    {
        type fixedValue;
    value $internalField;
     
    }
   
   
   
    pala2sopra
    {
       
      type    fixedValue;
      value    $internalField;
    }
   
   
   
    pala2sotto
    {
      type  fixedValue;
      value $internalField;
    }
   
   
   
    pala3
    {
        type    fixedValue ;
    value $internalField;
    }
   
   
   
    pala3sopra
    {
        type fixedValue;
    value $internalField;
     
    }
   
   
   
    pala3sotto
    {
    type fixedValue;
    value $internalField;
    }
   
   
    volumelato
    {
        type cyclicAMI;
    value $internalField;
    }
   
   
    volumesopra
    {
        type cyclicAMI;
    value $internalField;
    }
   
   
    volumesotto
    {
        type cyclicAMI;
    value $internalField;
    }
   
 
}

// ************************************************************************* //
p
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField  uniform 0;

boundaryField
{
    cubodx
    {
        type zeroGradient;
       
    }
   
   
    cubosopra
    {
        type zeroGradient;
       
    }
   
   
    cubosotto
    {
        type zeroGradient;
     
    }
   
   
    cubosx
    {
        type zeroGradient;
       
    }
   
   
    fisso
    {
        type cyclicAMI;
    value $internalField;
    }
   
   
    fissosopra
    {
        type cyclicAMI;
        value $internalField;
    }
   
   
    fissosotto
    {
        type cyclicAMI;
        value $internalField;
    }
   
   
    inlet
    {
        type zeroGradient;     
    }
   
   
    outlet
    {
        type fixedValue;
    value $internalField;       
    }
 
 
    pala1
    {
        type zeroGradient;
    }
   
   
    pala1sopra
    {
        type zeroGradient;       
    }
   
   
    pala1sotto
    {
    type zeroGradient;
    }
   
   
    pala2
    {
        type zeroGradient;     
    }
   
   
    pala2sopra
    {
       
      type    zeroGradient;
    }
   
   
    pala2sotto
    {
      type  zeroGradient;
    }
   
   
    pala3
    {
        type    zeroGradient;
    }
   
   
    pala3sopra
    {
        type zeroGradient;     
    }
   
   
    pala3sotto
    {
    type zeroGradient;
    }
   
   
    volumelato
    {
        type cyclicAMI;
    value $internalField;
    }
   
   
    volumesopra
    {
        type cyclicAMI;
    value $internalField;
    }
   
   
    volumesotto
    {
        type cyclicAMI;
    value $internalField;
    }
   
}

// ************************************************************************* //

jiahui_93 February 22, 2018 19:50
Hi, i am facing the same problem. May b you had found out the solution to this, do you mind to share it? thanks a lot.


All times are GMT -4. The time now is 04:41.