[wersoe]

Hello Foamers,

I am running a case with interDyMFoam with OF-2.0.x
It is simulating a shake flask, which is use in bioscience for cultivation of microorganisms and other cells... so, in simple words, its a conical flask fixed at a rotating table. If you look at the flask form one point, you see always the same point at the flask, a little like the movement of the moon... anyway, I manage to change the tankMixer tutorial so that I get the movement rite...

The problem:
With every timestep it looses a certain amount of fluid due to "mathematical diffusion" or something... it became obvious after about 5s simulation time, after 10-15 s no water is in the flask at all!
How can this happen? And can I avoid it?

The mesh is quite good, checkMesh shows no problems at all.
The residues ar quite low, see below...

I attached some of the important configuration files.

Hopefully sb has an idea what could happen here, since I am really stuck here at the moment...


solver.log, just two time step, it looks all time quite the same, time step size will be bigger later on during the simulation...

Interface Courant Number mean: 0.0066323 max: 0.498321
Courant Number mean: 0.0398762 max: 0.498321
deltaT = 7.7658e-07
Time = 0.000208016105042192333

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.000208016 transformation: ((0.000163375 2.66915e-07 0) (0.999999 (0 0 0.00163375)))
solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.000208016 transformation: ((0 0 0) (0.999999 (0 0 -0.00163375)))
solidBodyMotionFunctions::multiMotion::transformation(): Time = 0.000208016 transformation: ((0.000163375 2.66915e-07 0) (1 (0 0 0)))
Execution time for mesh.update() = 0.26 s
time step continuity errors : sum local = 1.96807e-13, global = -1.09833e-18, cumulative = -2.77411e-13
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 6.41289e-07, No Iterations 19
GAMGPCG:  Solving for pcorr, Initial residual = 0.0710413, Final residual = 5.60197e-07, No Iterations 10
GAMGPCG:  Solving for pcorr, Initial residual = 0.00429145, Final residual = 5.82585e-07, No Iterations 7
time step continuity errors : sum local = 3.62563e-18, global = -1.11019e-18, cumulative = -2.77412e-13
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 2.30226e-132  Max(alpha1) = 1
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 3.02039e-132  Max(alpha1) = 1
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 3.96591e-132  Max(alpha1) = 1
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 5.21193e-132  Max(alpha1) = 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.00237312, Final residual = 4.19164e-09, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.00328753, Final residual = 7.49817e-09, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.0020948, Final residual = 3.63414e-09, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 0.0253511, Final residual = 0.000118152, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.00189087, Final residual = 1.84691e-05, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 0.000119528, Final residual = 9.13697e-07, No Iterations 6
time step continuity errors : sum local = 2.13139e-10, global = -7.34516e-19, cumulative = -2.77413e-13
GAMG:  Solving for p_rgh, Initial residual = 0.00091004, Final residual = 2.82139e-06, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 4.74559e-05, Final residual = 2.29515e-07, No Iterations 3
GAMGPCG:  Solving for p_rgh, Initial residual = 3.12808e-06, Final residual = 1.99243e-09, No Iterations 4
time step continuity errors : sum local = 4.0557e-13, global = -7.34599e-19, cumulative = -2.77414e-13
DILUPBiCG:  Solving for omega, Initial residual = 0.00892048, Final residual = 7.19565e-09, No Iterations 4
DILUPBiCG:  Solving for k, Initial residual = 0.00479677, Final residual = 1.3686e-10, No Iterations 5
ExecutionTime = 902.75 s  ClockTime = 905 s

Interface Courant Number mean: 0.00666334 max: 0.498319
Courant Number mean: 0.0400127 max: 0.498319
deltaT = 7.79193e-07
Time = 0.000208795297929659983

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.000208795 transformation: ((0.000163987 2.68919e-07 0) (0.999999 (0 0 0.00163987)))
solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.000208795 transformation: ((0 0 0) (0.999999 (0 0 -0.00163987)))
solidBodyMotionFunctions::multiMotion::transformation(): Time = 0.000208795 transformation: ((0.000163987 2.68919e-07 0) (1 (0 0 0)))
Execution time for mesh.update() = 0.26 s
time step continuity errors : sum local = 4.06964e-13, global = 2.12725e-18, cumulative = -2.77411e-13
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 8.2677e-07, No Iterations 26
GAMGPCG:  Solving for pcorr, Initial residual = 0.0782515, Final residual = 7.38951e-07, No Iterations 21
GAMGPCG:  Solving for pcorr, Initial residual = 0.00468234, Final residual = 3.12732e-07, No Iterations 11
time step continuity errors : sum local = 4.6662e-18, global = 2.12035e-18, cumulative = -2.77409e-13
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 6.8496e-132  Max(alpha1) = 1
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 9.01104e-132  Max(alpha1) = 1
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 1.18666e-131  Max(alpha1) = 1
MULES: Solving for alpha1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.193844  Min(alpha1) = 1.5643e-131  Max(alpha1) = 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.00237884, Final residual = 3.78435e-09, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.00328912, Final residual = 7.99669e-09, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00209302, Final residual = 3.10299e-09, No Iterations 4
GAMG:  Solving for p_rgh, Initial residual = 0.0255634, Final residual = 0.000116615, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 0.00189883, Final residual = 1.77012e-05, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 0.000119078, Final residual = 8.93598e-07, No Iterations 6
time step continuity errors : sum local = 2.10798e-10, global = 2.13129e-19, cumulative = -2.77409e-13
GAMG:  Solving for p_rgh, Initial residual = 0.000931977, Final residual = 2.85761e-06, No Iterations 3
GAMG:  Solving for p_rgh, Initial residual = 4.78142e-05, Final residual = 2.35694e-07, No Iterations 3
GAMGPCG:  Solving for p_rgh, Initial residual = 3.16514e-06, Final residual = 1.28317e-09, No Iterations 4
time step continuity errors : sum local = 2.60862e-13, global = 2.1641e-19, cumulative = -2.77409e-13
DILUPBiCG:  Solving for omega, Initial residual = 0.00886702, Final residual = 5.8744e-09, No Iterations 4
DILUPBiCG:  Solving for k, Initial residual = 0.00479084, Final residual = 1.15469e-10, No Iterations 5
ExecutionTime = 910.15 s  ClockTime = 912 s