CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Free surface issues with interDyMFoam for hydroturbine

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 6, 2017, 14:05
Default Free surface issues with interDyMFoam for hydroturbine
  #1
New Member
 
oumnia
Join Date: Feb 2017
Posts: 12
Rep Power: 9
oumnion is on a distinguished road
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
Attached Images
File Type: png problem1.png (27.4 KB, 21 views)
oumnion is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
how to capture free surface cells? MichaelVS FLUENT 1 June 10, 2015 02:18
Sloshing Case - Probing or Sampling free surface during runtime Esweer OpenFOAM Running, Solving & CFD 0 February 17, 2015 03:34
Floating dok on waves: free surface divergence Galchenko CFX 3 September 11, 2014 18:04
Cluster ID's not contiguous in compute-nodes domain. ??? Shogan FLUENT 1 May 28, 2014 15:03
[snappyHexMesh] Layers don't fully surround surface EVBUCF OpenFOAM Meshing & Mesh Conversion 14 August 20, 2012 04:31


All times are GMT -4. The time now is 11:21.