求助：关于VOF求解！
-
我使用一种基于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}+...$
-
您好，我最近也在做这方面的工作，但是遇到了一些问题想向您请教下可以吗？
