各位老师同学好。本人最近在用InterPhaseChangeDyMFoam计算旋转机械（介质是水）的空化问题，边界条件是入口速度fixedValue ，压力totalpressure；出口速度inletOutlet，压力fixedFluxPressure，湍流采用标准K-e模型，额定的转速是303rad/s，入口段直径0.125m，速度给定5m/s，给定的k和epsilon初值分别是0.0284和0.09，时间步长设定1e-5，网格40W左右（非结构），用AMI连接域。在计算几步后就回出现崩溃，但是通过大幅度缩小时间步长和减少转速（大概到50rad）可以正常收敛，但是这两者都不太能接受。所以想寻求各位大神帮忙看看有没有更好的方法能稳定计算，或者如果是我的设置有问题也欢迎指正~~~~附上log和pyfoam的图
Create time
Create mesh for time = 0
Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: solidBody
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone impeller
PIMPLE: No convergence criteria found
PIMPLE: No corrector convergence criteria found
Calclations will do 3 corrections
Reading field p_rgh
Reading field U
Reading/calculating face flux field phi
AMI: Creating addressing and weights between 589 source faces and 754 target faces
AMI: Patch source sum(weights) min/max/average = 0.996846, 1, 0.999444
AMI: Patch target sum(weights) min/max/average = 0.99046, 1, 0.999443
AMI: Creating addressing and weights between 2752 source faces and 2836 target faces
AMI: Patch source sum(weights) min/max/average = 0.999976, 1.00012, 1.00003
AMI: Patch target sum(weights) min/max/average = 0.99994, 1.00012, 1.00003
Creating phaseChangeTwoPhaseMixture
Selecting phaseChange model SchnerrSauer
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
RASModel kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}
Reading g
Reading hRef
Calculating field g.h
No finite volume options present
Reading/calculating face velocity Uf
Courant Number mean: 3.46663e-05 max: 0.00375
Starting time loop
Courant Number mean: 3.46663e-05 max: 0.00375
Time = 1e-05
PIMPLE: Iteration 1
AMI: Creating addressing and weights between 589 source faces and 754 target faces
AMI: Patch source sum(weights) min/max/average = 0.99691, 1, 0.999455
AMI: Patch target sum(weights) min/max/average = 0.990625, 1, 0.999455
AMI: Creating addressing and weights between 2752 source faces and 2836 target faces
AMI: Patch source sum(weights) min/max/average = 0.999971, 1.00012, 1.00004
AMI: Patch target sum(weights) min/max/average = 0.999955, 1.00011, 1.00004
Execution time for mesh.update() = 0.28 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0046122, No Iterations 7
time step continuity errors : sum local = 7.2582e-05, global = -6.94718e-05, cumulative = -6.94718e-05
smoothSolver: Solving for alpha.water, Initial residual = 0.999973, Final residual = 1.1325, No Iterations 10
Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00343267, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.00357127, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 4.53115e-06, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 9.27694e-10, No Iterations 48
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 1, Final residual = 4.3514e-09, No Iterations 3
Phase-1 volume fraction = 0.979984 Min(alpha.water) = 0.821031 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.979982 Min(alpha.water) = 0.820777 Max(alpha.water) = 1
smoothSolver: Solving for Ux, Initial residual = 0.152882, Final residual = 0.000456751, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.158655, Final residual = 0.000456774, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.173638, Final residual = 0.000499279, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.00387877, Final residual = 9.67937e-10, No Iterations 8
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 0.000217683, Final residual = 5.3264e-11, No Iterations 2
Phase-1 volume fraction = 0.995947 Min(alpha.water) = 0.351358 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.995942 Min(alpha.water) = 0.345137 Max(alpha.water) = 1
smoothSolver: Solving for Ux, Initial residual = 0.0528343, Final residual = 6.27726e-07, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.0567381, Final residual = 7.23559e-07, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.120368, Final residual = 8.2518e-08, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 3.49149e-06, Final residual = 5.80041e-10, No Iterations 5
smoothSolver: Solving for epsilon, Initial residual = 0.0986499, Final residual = 5.16839e-07, No Iterations 4
smoothSolver: Solving for k, Initial residual = 1, Final residual = 1.29045e-07, No Iterations 6
ExecutionTime = 8.23 s ClockTime = 9 s
Courant Number mean: 0.0197059 max: 5.41516
Time = 2e-05
PIMPLE: Iteration 1
AMI: Creating addressing and weights between 589 source faces and 754 target faces
AMI: Patch source sum(weights) min/max/average = 0.996998, 1, 0.999472
AMI: Patch target sum(weights) min/max/average = 0.990859, 1, 0.999472
AMI: Creating addressing and weights between 2752 source faces and 2836 target faces
AMI: Patch source sum(weights) min/max/average = 0.999977, 1.00014, 1.00004
AMI: Patch target sum(weights) min/max/average = 0.999968, 1.00012, 1.00004
Execution time for mesh.update() = 0.26 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00504073, No Iterations 3
time step continuity errors : sum local = 0.000188015, global = -1.15051e-06, cumulative = -7.06224e-05
smoothSolver: Solving for alpha.water, Initial residual = 9.67477e-06, Final residual = 6.18883e-11, No Iterations 2
Phase-1 volume fraction = 0.995772 Min(alpha.water) = 0.211932 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.995769 Min(alpha.water) = 0.174408 Max(alpha.water) = 1
smoothSolver: Solving for Ux, Initial residual = 0.325971, Final residual = 0.0183609, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.309661, Final residual = 0.0165613, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.351942, Final residual = 0.0299024, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 1.43572e-07, Final residual = 7.23286e-10, No Iterations 5
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 0.0155352, Final residual = 2.2274e-09, No Iterations 3
Phase-1 volume fraction = 0.952704 Min(alpha.water) = 0.277133 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.952698 Min(alpha.water) = 0.228159 Max(alpha.water) = 1.00002
smoothSolver: Solving for Ux, Initial residual = 0.104737, Final residual = 0.00347544, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0999211, Final residual = 0.00465959, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.18458, Final residual = 0.0128277, No Iterations 2
GAMG: Solving for p_rgh, Initial residual = 0.0234743, Final residual = 7.77366e-10, No Iterations 5
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 0.00178812, Final residual = 5.32355e-09, No Iterations 2
Phase-1 volume fraction = 0.985629 Min(alpha.water) = 0.0661184 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.985616 Min(alpha.water) = 0.0661918 Max(alpha.water) = 1.00001
smoothSolver: Solving for Ux, Initial residual = 0.0821993, Final residual = 9.95813e-07, No Iterations 383
smoothSolver: Solving for Uy, Initial residual = 0.092795, Final residual = 9.84964e-07, No Iterations 401
smoothSolver: Solving for Uz, Initial residual = 0.174958, Final residual = 9.94153e-07, No Iterations 455
GAMG: Solving for p_rgh, Initial residual = 0.00206972, Final residual = 3.9387e-10, No Iterations 6
smoothSolver: Solving for epsilon, Initial residual = 0.393518, Final residual = 8.58378e-07, No Iterations 27
smoothSolver: Solving for k, Initial residual = 0.951336, Final residual = 9.76279e-07, No Iterations 374
ExecutionTime = 24.93 s ClockTime = 25 s
Courant Number mean: 0.0203512 max: 31.9041
Time = 3e-05
PIMPLE: Iteration 1
AMI: Creating addressing and weights between 589 source faces and 754 target faces
AMI: Patch source sum(weights) min/max/average = 0.997111, 1, 0.999493
AMI: Patch target sum(weights) min/max/average = 0.991164, 1, 0.999493
AMI: Creating addressing and weights between 2752 source faces and 2836 target faces
AMI: Patch source sum(weights) min/max/average = 0.999987, 1.00015, 1.00004
AMI: Patch target sum(weights) min/max/average = 0.999955, 1.00014, 1.00004
Execution time for mesh.update() = 0.22 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00576691, No Iterations 3
time step continuity errors : sum local = 0.000687359, global = -2.57908e-06, cumulative = -7.32014e-05
smoothSolver: Solving for alpha.water, Initial residual = 0.000155352, Final residual = 4.58115e-09, No Iterations 3
Phase-1 volume fraction = 0.985431 Min(alpha.water) = 0.00232962 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.985415 Min(alpha.water) = -0.0498938 Max(alpha.water) = 1.00005
smoothSolver: Solving for Ux, Initial residual = 0.556211, Final residual = 0.0456756, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.483701, Final residual = 0.0456852, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.54432, Final residual = 0.0380588, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 0.000953864, Final residual = 3.46933e-10, No Iterations 7
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 7.33141e-05, Final residual = 3.27551e-09, No Iterations 3
Phase-1 volume fraction = 0.985373 Min(alpha.water) = 0.00519769 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.985357 Min(alpha.water) = -0.0226666 Max(alpha.water) = 1.00002
smoothSolver: Solving for Ux, Initial residual = 0.251887, Final residual = 0.0202833, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.190243, Final residual = 0.0150066, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.216529, Final residual = 0.0209459, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 2.22683e-06, Final residual = 4.46301e-10, No Iterations 1
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 5.92958e-05, Final residual = 4.78988e-09, No Iterations 3
Phase-1 volume fraction = 0.985322 Min(alpha.water) = 0.00340138 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.985304 Min(alpha.water) = -0.0173934 Max(alpha.water) = 1.00005
smoothSolver: Solving for Ux, Initial residual = 0.300509, Final residual = 0.0025018, No Iterations 1000
smoothSolver: Solving for Uy, Initial residual = 0.398457, Final residual = 0.00194781, No Iterations 1000
smoothSolver: Solving for Uz, Initial residual = 0.353341, Final residual = 0.00647722, No Iterations 1000
GAMG: Solving for p_rgh, Initial residual = 2.82287e-06, Final residual = 4.84489e-10, No Iterations 3
smoothSolver: Solving for epsilon, Initial residual = 0.244345, Final residual = 8.33659e-07, No Iterations 29
smoothSolver: Solving for k, Initial residual = 0.933459, Final residual = 0.00774907, No Iterations 1000
ExecutionTime = 56.97 s ClockTime = 58 s
Courant Number mean: 0.039008 max: 94.2267
Time = 4e-05
PIMPLE: Iteration 1
AMI: Creating addressing and weights between 589 source faces and 754 target faces
AMI: Patch source sum(weights) min/max/average = 0.997235, 1, 0.999519
AMI: Patch target sum(weights) min/max/average = 0.99154, 1, 0.999518
AMI: Creating addressing and weights between 2752 source faces and 2836 target faces
AMI: Patch source sum(weights) min/max/average = 0.999987, 1.00015, 1.00005
AMI: Patch target sum(weights) min/max/average = 0.999955, 1.00015, 1.00005
Execution time for mesh.update() = 0.28 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00975609, No Iterations 5
time step continuity errors : sum local = 0.00591848, global = 3.92436e-06, cumulative = -6.92771e-05
smoothSolver: Solving for alpha.water, Initial residual = 0.000372595, Final residual = 7.54046e-09, No Iterations 8
Phase-1 volume fraction = 0.984937 Min(alpha.water) = 0.00110965 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.984913 Min(alpha.water) = -2.29641 Max(alpha.water) = 2.59561
smoothSolver: Solving for Ux, Initial residual = 0.541575, Final residual = 0.0121456, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.62719, Final residual = 0.0557338, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.462147, Final residual = 0.0348564, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 3.06236e-05, Final residual = 2.15551e-10, No Iterations 2
PIMPLE: Iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 0.00637181, Final residual = 5.1554e-08, No Iterations 10
Phase-1 volume fraction = 0.984205 Min(alpha.water) = 1.05723e-20 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.982971 Min(alpha.water) = -1336.22 Max(alpha.water) = 97.6268
smoothSolver: Solving for Ux, Initial residual = 0.351934, Final residual = 0.0239104, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.318012, Final residual = 0.0115367, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.337556, Final residual = 0.03052, No Iterations 1
GAMG: Solving for p_rgh, Initial residual = 4.60749e-08, Final residual = 3.66729e-10, No Iterations 1
PIMPLE: Iteration 3
smoothSolver: Solving for alpha.water, Initial residual = 0.998656, Final residual = 9.85984e-09, No Iterations 7
Phase-1 volume fraction = 0.983918 Min(alpha.water) = -0.000242279 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Liquid phase volume fraction = 0.984043 Min(alpha.water) = -16.1095 Max(alpha.water) = 105.668
[1] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[1] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[1] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[1] #6 ? at ??:?
[1] #7 ? at ??:?
[1] #8 ? at ??:?
[1] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #10 ? at ??:?
1_1537536681529_residual.PNG 0_1537536681528_continuity.PNG
