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/)
-   -   Free surface issues with interDyMFoam for hydroturbine (https://www.cfd-online.com/Forums/openfoam-solving/193928-free-surface-issues-interdymfoam-hydroturbine.html)

oumnion October 6, 2017 15:05

Free surface issues with interDyMFoam for hydroturbine
 
1 Attachment(s)
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



All times are GMT -4. The time now is 12:28.