CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   interDyMFoam - change in volume fraction (http://www.cfd-online.com/Forums/openfoam-solving/64020-interdymfoam-change-volume-fraction.html)

gopala April 27, 2009 10:46

interDyMFoam - change in volume fraction
 
Hello,

I tried running the the case "sloshingTank3D6DoF" in the tutorial under interDyMFoam. I understand that this is a closed domain problem and when I run the case, I see that the volume fraction in the domain is changing every time step!! I have pasted the output for this case for few time steps. Can someone please explain what is happening here??

************************************************** ***
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
deltaT = 0.01
Time = 0.01
solidBodyMotionFunctions::SKA::transformation(): Time = 0.01 transformation: ((0.00993214 0.0235843 0.00796523) (0.999999 (0.00104291 0.000602303 0.000433999)))
Execution time for mesh.update() = 0.07 s
MULES: Solving for gamma
Liquid phase volume fraction = 0.337111 Min(gamma) = 0 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.337012 Min(gamma) = 0 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.336913 Min(gamma) = 0 Max(gamma) = 1
GAMG: Solving for pd, Initial residual = 1, Final residual = 0.00591554, No Iterations 5
time step continuity errors : sum local = 5.98113e-05, global = 2.86642e-17, cumulative = 2.86642e-17
GAMGPCG: Solving for pd, Initial residual = 0.0591541, Final residual = 1.29214e-09, No Iterations 8
time step continuity errors : sum local = 3.39181e-11, global = 2.8702e-17, cumulative = 5.73662e-17
ExecutionTime = 1.34 s ClockTime = 2 s
Courant Number mean: 0.0178644 max: 0.159383
deltaT = 0.01
Time = 0.02
solidBodyMotionFunctions::SKA::transformation(): Time = 0.02 transformation: ((0.0198643 0.0471685 0.0159305) (0.999997 (0.00208633 0.0012037 0.000869254)))
Execution time for mesh.update() = 0.07 s
MULES: Solving for gamma
Liquid phase volume fraction = 0.336914 Min(gamma) = 0 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.336914 Min(gamma) = -2.54116e-39 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.336915 Min(gamma) = -2.54116e-39 Max(gamma) = 1
GAMG: Solving for pd, Initial residual = 0.188129, Final residual = 0.00119762, No Iterations 6
time step continuity errors : sum local = 2.73703e-05, global = 2.91035e-16, cumulative = 3.48401e-16
GAMGPCG: Solving for pd, Initial residual = 0.0687751, Final residual = 1.39387e-09, No Iterations 8
time step continuity errors : sum local = 2.24666e-11, global = 2.91099e-16, cumulative = 6.39499e-16
ExecutionTime = 2.36 s ClockTime = 3 s
Courant Number mean: 0.0178365 max: 0.148186
deltaT = 0.01
Time = 0.03
solidBodyMotionFunctions::SKA::transformation(): Time = 0.03 transformation: ((0.0297964 0.0707528 0.0238957) (0.999993 (0.00313028 0.00180419 0.00130576)))
Execution time for mesh.update() = 0.07 s
MULES: Solving for gamma
Liquid phase volume fraction = 0.336915 Min(gamma) = -2.68208e-31 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.336915 Min(gamma) = -8.86794e-34 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.336916 Min(gamma) = -3.81862e-30 Max(gamma) = 1
GAMG: Solving for pd, Initial residual = 0.139871, Final residual = 0.00123457, No Iterations 3
time step continuity errors : sum local = 1.68448e-05, global = -4.18923e-16, cumulative = 2.20576e-16
GAMGPCG: Solving for pd, Initial residual = 0.0201855, Final residual = 2.87405e-10, No Iterations 8
time step continuity errors : sum local = 3.50033e-12, global = -4.18868e-16, cumulative = -1.98292e-16
ExecutionTime = 3.34 s ClockTime = 4 s
Courant Number mean: 0.0178408 max: 0.140226
deltaT = 0.01
Time = 0.04
solidBodyMotionFunctions::SKA::transformation(): Time = 0.04 transformation: ((0.0397285 0.0943371 0.0318609) (0.999987 (0.00417474 0.00240376 0.00174352)))
Execution time for mesh.update() = 0.07 s
MULES: Solving for gamma
Liquid phase volume fraction = 0.336916 Min(gamma) = -8.91948e-30 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.336917 Min(gamma) = -1.30713e-31 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.336917 Min(gamma) = -3.99721e-34 Max(gamma) = 1
GAMG: Solving for pd, Initial residual = 0.107624, Final residual = 0.00105195, No Iterations 2
time step continuity errors : sum local = 1.15421e-05, global = 9.4616e-16, cumulative = 7.47868e-16
GAMGPCG: Solving for pd, Initial residual = 0.00677147, Final residual = 4.96459e-10, No Iterations 7
time step continuity errors : sum local = 5.08561e-12, global = 9.46159e-16, cumulative = 1.69403e-15
ExecutionTime = 4.24 s ClockTime = 5 s
************************************************** *********


All times are GMT -4. The time now is 03:31.