|
[Sponsors] |
Free surface issues with interDyMFoam for hydroturbine |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 6, 2017, 15:05 |
Free surface issues with interDyMFoam for hydroturbine
|
#1 |
New Member
oum
Join Date: Feb 2017
Posts: 14
Rep Power: 9 |
Hi, I am a beginner in OpenFOAM and my case involves a hydroturbine and a free surface. The simulation runs for a little bit of time and then crashes with high values for the velocity and pressure. Can someone help me please thank you.
Code:
Create mesh for time = 0 Selecting dynamicFvMesh dynamicRefineFvMesh Detected 5635958 cells that are protected from refinement. Writing these to cellSet protectedCells. PIMPLE: max iterations = 7 field U : relTol 0.001, tolerance 0.01 field p_rgh : relTol 0.001, tolerance 0.01 Reading field p_rgh Reading field U Using dynamicCode for patch Inlet on field U at line 36 in "/cavs/users/oumnia/test32/0/U.boundaryField.Inlet" Reading/calculating face flux field phi AMI: Creating addressing and weights between 18784 source faces and 18784 target faces AMI: Patch source sum(weights) min/max/average = 1, 1, 1 AMI: Patch target sum(weights) min/max/average = 1, 1, 1 Reading transportProperties Using dynamicCode for patch Inlet on field alpha.water at line 38 in "/cavs/users/oumnia/test32/0/alpha.water.boundaryField.Inlet" Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kOmegaSST Selecting patchDistMethod meshWave kOmegaSSTCoeffs { alphaK1 0.85; alphaK2 1; alphaOmega1 0.5; alphaOmega2 0.856; gamma1 0.555556; gamma2 0.44; beta1 0.075; beta2 0.0828; betaStar 0.09; a1 0.31; b1 1; c1 10; F3 false; } Reading g Reading hRef Calculating field g.h No MRF models present No finite volume options present GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0 GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0 GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0 GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0, global = 0, cumulative = 0 Reading/calculating face velocity Uf Courant Number mean: 0 max: 0 Starting time loop Courant Number mean: 0 max: 0 Interface Courant Number mean: 0 max: 0 Time = 3.89e-05 PIMPLE: iteration 1 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0 Phase-1 volume fraction = 0.825016 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = 0 Max(alpha.water) = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 1.37167e-08, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 1.89475e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 1.47204e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 7.23422e-05, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.000531344, Final residual = 2.38374e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.000248164, Final residual = 1.0294e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.000157466, Final residual = 7.47826e-09, No Iterations 7 time step continuity errors : sum local = 5.60832e-09, global = 1.20776e-09, cumulative = 1.20776e-09 PIMPLE: iteration 2 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 0.000109362, Final residual = 7.28451e-10, No Iterations 2 Phase-1 volume fraction = 0.825016 Min(alpha.water) = 0 Max(alpha.water) = 1.00088 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.000302918 Max(alpha.water) = 1.00087 smoothSolver: Solving for Ux, Initial residual = 0.357829, Final residual = 9.02457e-07, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 0.11343, Final residual = 1.48911e-07, No Iterations 3 smoothSolver: Solving for Uz, Initial residual = 0.0654346, Final residual = 4.70941e-07, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.0018622, Final residual = 2.05016e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00028634, Final residual = 8.39386e-09, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 0.000169115, Final residual = 6.87957e-09, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 0.000125568, Final residual = 8.18993e-09, No Iterations 5 time step continuity errors : sum local = 4.87568e-09, global = 6.93835e-10, cumulative = 1.90159e-09 PIMPLE: iteration 3 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 4.92015e-05, Final residual = 8.61516e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -2.16058e-05 Max(alpha.water) = 1.00075 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.0015319 Max(alpha.water) = 1.00075 smoothSolver: Solving for Ux, Initial residual = 0.244735, Final residual = 1.2666e-07, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 0.0267965, Final residual = 1.93774e-07, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0217409, Final residual = 1.04416e-07, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.000457317, Final residual = 7.0901e-09, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.000149006, Final residual = 8.39895e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 9.53708e-05, Final residual = 6.58676e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 7.76506e-05, Final residual = 8.30865e-09, No Iterations 5 time step continuity errors : sum local = 4.43436e-09, global = -5.25479e-10, cumulative = 1.37611e-09 PIMPLE: iteration 4 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 5.93652e-05, Final residual = 5.27156e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -2.20178e-05 Max(alpha.water) = 1.00051 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.000272849 Max(alpha.water) = 1.00051 smoothSolver: Solving for Ux, Initial residual = 0.200215, Final residual = 6.86747e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0104848, Final residual = 4.56871e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00910982, Final residual = 3.70008e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.000177334, Final residual = 8.14144e-09, No Iterations 9 GAMG: Solving for p_rgh, Initial residual = 9.83282e-05, Final residual = 6.0129e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 6.85843e-05, Final residual = 5.33475e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 5.848e-05, Final residual = 7.19041e-09, No Iterations 5 time step continuity errors : sum local = 3.7612e-09, global = -3.88068e-10, cumulative = 9.88044e-10 PIMPLE: iteration 5 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 2.35613e-05, Final residual = 1.93274e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -1.47201e-05 Max(alpha.water) = 1.00029 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.000830192 Max(alpha.water) = 1.00028 smoothSolver: Solving for Ux, Initial residual = 0.208086, Final residual = 3.73729e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00896614, Final residual = 2.3037e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00781102, Final residual = 1.64283e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.000103115, Final residual = 8.79946e-09, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 6.41372e-05, Final residual = 3.99713e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 5.20943e-05, Final residual = 6.54088e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 4.67997e-05, Final residual = 5.44332e-09, No Iterations 5 time step continuity errors : sum local = 2.85622e-09, global = -1.95179e-10, cumulative = 7.92865e-10 PIMPLE: iteration 6 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 2.35889e-05, Final residual = 1.66659e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -1.56693e-05 Max(alpha.water) = 1.00018 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.000223467 Max(alpha.water) = 1.00018 smoothSolver: Solving for Ux, Initial residual = 0.20964, Final residual = 2.92946e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00914061, Final residual = 1.65186e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0077796, Final residual = 1.26689e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 6.31289e-05, Final residual = 3.60904e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 4.73646e-05, Final residual = 6.30934e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 4.21843e-05, Final residual = 6.08149e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 3.92269e-05, Final residual = 5.87379e-09, No Iterations 5 time step continuity errors : sum local = 3.09036e-09, global = -2.87458e-11, cumulative = 7.64119e-10 PIMPLE: iteration 7 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 2.22449e-05, Final residual = 1.49512e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -7.46156e-06 Max(alpha.water) = 1.00016 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.000128277 Max(alpha.water) = 1.00016 smoothSolver: Solving for Ux, Initial residual = 0.206815, Final residual = 2.5601e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00898977, Final residual = 1.33453e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00759579, Final residual = 1.07833e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 4.51636e-05, Final residual = 3.48322e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 3.88732e-05, Final residual = 5.96812e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 3.60665e-05, Final residual = 5.66968e-09, No Iterations 5 GAMG: Solving for p_rgh, Initial residual = 3.42041e-05, Final residual = 5.51921e-09, No Iterations 5 time step continuity errors : sum local = 2.90642e-09, global = 1.05464e-10, cumulative = 8.69583e-10 smoothSolver: Solving for omega, Initial residual = 6.8874e-06, Final residual = 9.55243e-08, No Iterations 1 bounding omega, min: -4931.38 max: 135605 average: 629.038 smoothSolver: Solving for k, Initial residual = 0.999999, Final residual = 5.06221e-07, No Iterations 4 PIMPLE: not converged within 7 iterations ExecutionTime = 667.79 s ClockTime = 683 s Courant Number mean: 0.000718967 max: 4.00731 Interface Courant Number mean: 2.58857e-06 max: 1.32923 Time = 7.78e-05 PIMPLE: iteration 1 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 7.58555e-05, Final residual = 3.5615e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -3.08418e-06 Max(alpha.water) = 1.00027 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -6.17318e-05 Max(alpha.water) = 1.00027 smoothSolver: Solving for Ux, Initial residual = 0.201974, Final residual = 4.80393e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.134987, Final residual = 1.01987e-07, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.128283, Final residual = 1.12547e-07, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 8.33406e-05, Final residual = 7.48228e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00746883, Final residual = 1.15243e-07, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00217993, Final residual = 4.15288e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.0023391, Final residual = 1.27941e-08, No Iterations 10 time step continuity errors : sum local = 1.04617e-10, global = -3.71681e-11, cumulative = 8.32414e-10 PIMPLE: iteration 2 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 3.5829e-05, Final residual = 3.18921e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -9.77455e-06 Max(alpha.water) = 1.00015 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.000219096 Max(alpha.water) = 1.00015 smoothSolver: Solving for Ux, Initial residual = 0.102698, Final residual = 2.47282e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0111616, Final residual = 1.73245e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.0103395, Final residual = 1.85102e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.000225834, Final residual = 9.37608e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.0010878, Final residual = 6.83201e-09, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.000521908, Final residual = 9.29206e-09, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.000359561, Final residual = 5.87849e-09, No Iterations 7 time step continuity errors : sum local = 2.83395e-10, global = -4.90774e-11, cumulative = 7.83337e-10 PIMPLE: iteration 3 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 2.85833e-05, Final residual = 3.06568e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -5.37227e-06 Max(alpha.water) = 1.00014 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.00012904 Max(alpha.water) = 1.00014 smoothSolver: Solving for Ux, Initial residual = 0.0629293, Final residual = 1.98454e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00915738, Final residual = 1.39455e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00866795, Final residual = 1.38746e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.0010695, Final residual = 1.27654e-07, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00109508, Final residual = 6.59077e-09, No Iterations 9 GAMG: Solving for p_rgh, Initial residual = 0.000835281, Final residual = 7.00495e-09, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.000718826, Final residual = 7.62022e-09, No Iterations 8 time step continuity errors : sum local = 1.47909e-10, global = -1.88094e-11, cumulative = 7.64528e-10 PIMPLE: iteration 4 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 2.8368e-05, Final residual = 3.13901e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -3.30738e-06 Max(alpha.water) = 1.00014 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -0.000107307 Max(alpha.water) = 1.00014 smoothSolver: Solving for Ux, Initial residual = 0.0249454, Final residual = 1.67952e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00504684, Final residual = 1.04633e-08, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00472598, Final residual = 1.13174e-08, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.00186453, Final residual = 6.03328e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00272099, Final residual = 1.62492e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00231995, Final residual = 1.14266e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00217689, Final residual = 1.17336e-08, No Iterations 10 time step continuity errors : sum local = 6.592e-11, global = -2.12672e-12, cumulative = 7.62401e-10 PIMPLE: iteration 5 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 2.81914e-05, Final residual = 3.1524e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -9.71307e-06 Max(alpha.water) = 1.00016 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.825016 Min(alpha.water) = -8.93482e-05 Max(alpha.water) = 1.00016 smoothSolver: Solving for Ux, Initial residual = 0.00738903, Final residual = 1.49897e-07, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00123566, Final residual = 7.68462e-09, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00113367, Final residual = 8.72918e-09, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.00256004, Final residual = 9.47849e-08, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00143738, Final residual = 8.81855e-09, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00130557, Final residual = 8.34628e-09, No Iterations 10 GAMG: Solving for p_rgh, Initial residual = 0.00125098, Final residual = 7.26039e-09, No Iterations 10 time step continuity errors : sum local = 6.58837e-11, global = 1.09376e-12, cumulative = 7.63495e-10 PIMPLE: iteration 6 Selected 0 split points out of a possible 0. smoothSolver: Solving for alpha.water, Initial residual = 2.81377e-05, Final residual = 3.17355e-09, No Iterations 1 Phase-1 volume fraction = 0.825016 Min(alpha.water) = -1.28893e-05 Max(alpha.water) = 1.00016 Applying the previous iteration compression flux MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 127806 Min(alpha.water) = -8.85446 Max(alpha.water) = 4.21508e+14 |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
how to capture free surface cells? | MichaelVS | FLUENT | 1 | June 10, 2015 03:18 |
Sloshing Case - Probing or Sampling free surface during runtime | Esweer | OpenFOAM Running, Solving & CFD | 0 | February 17, 2015 04:34 |
Floating dok on waves: free surface divergence | Galchenko | CFX | 3 | September 11, 2014 19:04 |
Cluster ID's not contiguous in compute-nodes domain. ??? | Shogan | FLUENT | 1 | May 28, 2014 16:03 |
[snappyHexMesh] Layers don't fully surround surface | EVBUCF | OpenFOAM Meshing & Mesh Conversion | 14 | August 20, 2012 05:31 |