求助：关于VOF求解！

[haozz]

我使用一种基于Interfoam改良的粘弹性求解器，求解粘弹性液滴，但是计算总是终止，由于是新手log文件看的不是很明白，想知道这样的信息是代表了什么？（log文件最后几行如下）

smoothSolver:  Solving for alpha.water, Initial residual = 1.5329e-07, Final residual = 3.889e-13, No Iterations 1
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45203e-114  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45203e-114  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 1.53284e-07, Final residual = 3.86063e-13, No Iterations 1
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45175e-114  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45175e-114  Max(alpha.water) = 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.476091, Final residual = 7.54132e-11, No Iterations 27
DILUPBiCG:  Solving for Uy, Initial residual = 0.213985, Final residual = 4.74027e-11, No Iterations 26
DILUPBiCG:  Solving for Uz, Initial residual = 0.0136921, Final residual = 7.92877e-11, No Iterations 23
DICPCG:  Solving for p_rgh, Initial residual = 0.0156261, Final residual = 9.4989e-11, No Iterations 274
time step continuity errors : sum local = 3.66088e-12, global = 5.11974e-14, cumulative = -5.10157e-09
DILUPBiCG:  Solving for theta.waterxx, Initial residual = 0.0142458, Final residual = 3.43368e-11, No Iterations 6
DILUPBiCG:  Solving for theta.waterxy, Initial residual = 0.0904855, Final residual = 7.73445e-12, No Iterations 7
DILUPBiCG:  Solving for theta.waterxz, Initial residual = 0.0162719, Final residual = 2.1972e-11, No Iterations 6
DILUPBiCG:  Solving for theta.wateryy, Initial residual = 0.0303967, Final residual = 8.85733e-11, No Iterations 6
DILUPBiCG:  Solving for theta.wateryz, Initial residual = 0.00425512, Final residual = 4.27331e-12, No Iterations 6
DILUPBiCG:  Solving for theta.waterzz, Initial residual = 0.00268117, Final residual = 3.1176e-11, No Iterations 6
ExecutionTime = 810.68 s  ClockTime = 825 s

Courant Number mean: 0.0013613 max: 9.43911
Interface Courant Number mean: 5.0886e-06 max: 0.0878234
deltaT = 2.42949e-07
Time = 0.0594116

smoothSolver:  Solving for alpha.water, Initial residual = 1.62073e-08, Final residual = 5.71768e-13, No Iterations 1
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45164e-114  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45164e-114  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 1.61967e-08, Final residual = 5.48116e-13, No Iterations 1
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45153e-114  Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.111576  Min(alpha.water) = -1.45153e-114  Max(alpha.water) = 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.999514, Final residual = 2.07976e-11, No Iterations 8
DILUPBiCG:  Solving for Uy, Initial residual = 0.996533, Final residual = 2.09498e-11, No Iterations 8
DILUPBiCG:  Solving for Uz, Initial residual = 0.943074, Final residual = 1.95211e-11, No Iterations 8
DICPCG:  Solving for p_rgh, Initial residual = 0.999913, Final residual = 9.31371e-11, No Iterations 290
time step continuity errors : sum local = 1.72859e-10, global = -5.15326e-13, cumulative = -5.10208e-09
DILUPBiCG:  Solving for theta.waterxx, Initial residual = 0.743477, Final residual = 8.08487e-11, No Iterations 21
DILUPBiCG:  Solving for theta.waterxy, Initial residual = 0.910044, Final residual = 3.06345e-11, No Iterations 22
DILUPBiCG:  Solving for theta.waterxz, Initial residual = 0.845722, Final residual = 3.71998e-11, No Iterations 21
DILUPBiCG:  Solving for theta.wateryy, Initial residual = 0.968814, Final residual = 3.85885e-11, No Iterations 25
DILUPBiCG:  Solving for theta.wateryz, Initial residual = 0.899894, Final residual = 7.95995e-11, No Iterations 24
DILUPBiCG:  Solving for theta.waterzz, Initial residual = 0.595067, Final residual = 6.12874e-11, No Iterations 20

[李东岳]

求解粘弹性流体的剪切应力传输方程：$\frac{\p\tau}{\p t}+...$

[sungda]

您好，我最近也在做这方面的工作，但是遇到了一些问题想向您请教下可以吗？