用InterPhaseChangeDyMFoam求解器计算崩溃
-
各位老师同学好。本人最近在用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 ??:?
