用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 ??:?
    
    

    1_1537536681529_residual.PNG 0_1537536681528_continuity.PNG


Log in to reply