bluecfd运算出错



  • 使用pimpleDyMFoam求解器计算风机所受转矩,使用kωSST模型,运算第一步就出现如下错误:

    Create mesh for time = 0
    
    Selecting dynamicFvMesh dynamicMotionSolverFvMesh
    Selecting motion solver: solidBody
    Selecting solid-body motion function rotatingMotion
    Applying solid body motion to cellZone inter
    
    PIMPLE: no residual control data found. Calculations will employ 2 corrector loops
    
    Reading field p
    
    AMI: Creating addressing and weights between 14390 source faces and 12742 target faces
    AMI: Patch source sum(weights) min/max/average = 0.570448, 1.45816, 1.00126
    AMI: Patch target sum(weights) min/max/average = 0.000958814, 1.18605, 0.993037
    Reading field U
    
    Reading/calculating face flux field phi
    
    Selecting incompressible transport model Newtonian
    Selecting turbulence model type RAS
    Selecting RAS turbulence model kOmegaSST
    Selecting patchDistMethod meshWave
    RAS
    {
        RASModel        kOmegaSST;
        turbulence      on;
        printCoeffs     on;
        alphaK1         0.85;
        alphaK2         1;
        alphaOmega1     0.5;
        alphaOmega2     0.856;
        gamma1          0.555556;
        gamma2          0.44;
        beta1           0.075;
        beta2           0.0828;
        betaStar        0.09;
        a1              0.31;
        b1              1;
        c1              10;
        F3              false;
    }
    
    No MRF models present
    
    Reading/calculating face velocity Uf
    
    No finite volume options present
    
    Courant Number mean: 0.0006176 max: 0.225383
    
    Starting time loop
    
    forces Forces:
        Not including porosity effects
    forces forceCoeffs:
        Not including porosity effects
    forceCoeffs forceCoeffs:
        Not including porosity effects
    Courant Number mean: 0.0006176 max: 0.225383
    Time = 0.001
    
    AMI: Creating addressing and weights between 14390 source faces and 12742 target faces
    AMI: Patch source sum(weights) min/max/average = 0.570366, 1.45766, 1.00147
    AMI: Patch target sum(weights) min/max/average = 0, 1.18298, 0.993152
    GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.0220577, No Iterations 4
    GAMG:  Solving for pcorr, Initial residual = 0.230843, Final residual = 0.0159595, No Iterations 1
    GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
    GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
    GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
    GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
    time step continuity errors : sum local = 2.90421e-009, global = -5.82674e-012, cumulative = -5.82674e-012
    PIMPLE: iteration 1
    smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0228862, No Iterations 1
    smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0203193, No Iterations 1
    smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0226345, No Iterations 1
    Generating stack trace...
    
    
    Backtrace:
            ZN10StackTraceC1Ev [0x705c1465+0x25]
                     module: D:\blueCFD-Core-2017\ThirdParty-5.x\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll
            ZN4Foam5error10printStackERNS_7OstreamE [0x18d1c88+0x218]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
            ZN4Foam6sigFpe13sigFpeHandlerEi [0x18d2af3+0x33]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
            (No symbol) [0x4047cd]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
            _C_specific_handler [0x7ffe9fcc7c58+0x98]
                     module: C:\WINDOWS\System32\msvcrt.dll
            0_chkstk [0x7ffe9ff8f7dd+0x11d]
                     module: C:\WINDOWS\SYSTEM32\ntdll.dll
            RtlWalkFrameChain [0x7ffe9fefd856+0x13f6]
                     module: C:\WINDOWS\SYSTEM32\ntdll.dll
            KiUserExceptionDispatcher [0x7ffe9ff8e70e+0x2e]
                     module: C:\WINDOWS\SYSTEM32\ntdll.dll
            ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE [0x1877fa4+0x24]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
            (No symbol) [0x42c1e7]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
            (No symbol) [0x4361db]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
            (No symbol) [0x4363dd]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
            (No symbol) [0x4546ca]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
            (No symbol) [0x4013f7]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
            (No symbol) [0x40152b]
                     module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
            BaseThreadInitThunk [0x7ffe9d393dc4+0x14]
                     module: C:\WINDOWS\System32\KERNEL32.DLL
            RtlUserThreadStart [0x7ffe9ff63691+0x21]
                     module: C:\WINDOWS\SYSTEM32\ntdll.dll
    

    求教各位大神,这是哪里的设置出现问题?



  • openfoam报错如下:

    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::sigFpe::sigHandler(int) at ??:?
    #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
    #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
    #5  ? at ??:?
    #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #7  ? at ??:?
    


  • 试试把你controlDict里面附加的function删掉



  • @东岳 还是一样的错误



  • @东岳 请问还有可能是哪里的设置有问题呢?


  • OpenFOAM教授

    @lv1995 能把完整的输出信息贴出来吗?



  • @dzw05
    k:

    dimensions      [0 2 -2 0 0 0 0];
    internalField   uniform 0.081;
    boundaryField
    { 
        inlet
        {
            type            zeroGradient;
        }
        outlet
        {
            type            zeroGradient;
        }
        frontAndBack
        {
            type            symmetry;
         }
        inter2
        {
            type            cyclicAMI;
         }
        inter1
        {
            type            cyclicAMI;
         }
    
        blades
        {
            type            kqRWallFunction;
            value           $internalField;
        }
    }
    

    nut

    dimensions      [0 2 -1 0 0 0 0];
    
    internalField   uniform 0;
    
    boundaryField
    {
            
        inlet
        {
            type            fixedValue;
            value           $internalField;
        }
    
        outlet
        {
            type             calculated;
            value           uniform 0;
    
        }
        frontAndBack
        {
            type            symmetry;
         }
        inter2
        {
            type            cyclicAMI;
         }
        inter1
        {
            type            cyclicAMI;
         }
         blades
        {
            type            nutkWallFunction;
            value           uniform 0;
        }
    }
    

    omega

    dimensions      [0 0 -1 0 0 0 0];
    
    internalField   uniform 0.095;
    
    boundaryField
    {
        
    
        outlet
        {
            type            zeroGradient;
        }
        blades
        {
            type            omegaWallFunction;
            value           $internalField;
        }
        inlet
        {
            type            zeroGradient;
        }
            frontAndBack
        {
            type            symmetry;
         }
        inter2
        {
            type            cyclicAMI;
         }
        inter1
        {
            type            cyclicAMI;
         }
    }
    

    p

    dimensions      [0 2 -2 0 0 0 0];
    
    internalField   uniform 0;
    
    boundaryField
    {
        inlet
        {
            type            zeroGradient;
        }
    
        outlet
        {
            type            fixedValue;
            value           uniform 0;
       }
            frontAndBack
        {
            type            symmetry;
         }
        inter2
        {
            type            cyclicAMI;
         }
        inter1
        {
            type            cyclicAMI;
         }
    
    
        blades
        {
            type            zeroGradient;
    
        }
    }
    

    u

    dimensions      [0 1 -1 0 0 0 0];
    
    internalField   uniform (14 0 0);
    
    boundaryField
    {
        
        inlet
        {
            type            fixedValue;
    
            value      uniform (14 0 0);
        }
    
        outlet
        {
            type            zeroGradient;
        }
            frontAndBack
        {
            type            symmetry;
         }
        inter2
        {
            type            cyclicAMI;
         }
        inter1
        {
            type            cyclicAMI;
         }
    
    
         blades
        {
            type            rotatingWallVelocity;
    
           origin           (0 0 0);
    
           axis             (0 0 1);
    
           omega           0;
        }
    }
    

    boundary

    inlet
        {
            type            patch;
            nFaces          400;
            startFace       6377766;
        }
        outlet
        {
            type            patch;
            nFaces          400;
            startFace       6378166;
        }
        frontAndBack
        {
            type            symmetry;
            inGroups        1(symmetry);
            nFaces          4800;
            startFace       6378566;
        }
        blades
        {
            type            wall;
            inGroups        1(wall);
            nFaces          80334;
            startFace       6383366;
        }
        inter1
        {
            type            cyclicAMI;
            inGroups        1(cyclicAMI);
            nFaces          13994;
            startFace       6463700;
            matchTolerance  0.0001;
            transform       noOrdering;
            neighbourPatch  inter2;
        }
        inter2
        {
            type            cyclicAMI;
            inGroups        1(cyclicAMI);
            nFaces          11572;
            startFace       6477694;
            matchTolerance  0.0001;
            transform       noOrdering;
            neighbourPatch  inter1;
        }
    

    fvSchemes

    ddtSchemes
    {
        default Euler;
    }
    
    gradSchemes
    {
        default         Gauss linear;
        grad(p)         Gauss linear;
        grad(U)         Gauss linear;
    }
    
    divSchemes
    {
        default         none;
        div(phi,U)      Gauss linearUpwind grad(U);
        div(phi,k)      Gauss limitedLinear 1;
        div(phi,omega)  Gauss limitedLinear 1;
        div((nuEff*dev2(T(grad(U))))) Gauss linear;
    }
    
    laplacianSchemes
    {
        default         Gauss linear limited corrected 0.5;
    }
    
    interpolationSchemes
    {
        default         linear;
    }
    
    snGradSchemes
    {
        default         corrected;
    }
    
    wallDist
    {
        method meshWave;
    }
    

    fvSolution

    solvers
    {
        "pcorr.*"
        {
            solver           GAMG;
            tolerance        0.02;
            relTol           0;
            smoother         GaussSeidel;
        }
    
        p
        {
            $pcorr;
            tolerance        1e-7;
            relTol           0.01;
        }
    
        pFinal
        {
            $p;
            tolerance        1e-7;
            relTol           0;
        }
    
        "(U|k|omega)"
        {
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1e-06;
            relTol          0.1;
        }
    
        "(U|k|omega)Final"
        {
            $U;
            tolerance       1e-06;
            relTol          0;
        }
    
        cellDisplacement
        {
            solver          GAMG;
            tolerance       1e-5;
            relTol          0;
            smoother        GaussSeidel;
        }
    }
    
    PIMPLE
    {
        correctPhi          yes;
        nOuterCorrectors    2;
        nCorrectors         2;
        nNonOrthogonalCorrectors 5;
        pRefCell  0;
        pRefPoint 0;
        pRefValue 0;
    }
    
    relaxationFactors
    {
        fields
        {
            p               0.2;
        }
        equations
        {
            "(U|k|omega)"   0.7;
            "(U|k|omega)Final" 0.7;
        }
    }
    
    cache
    {
        grad(U);
    }
    

Log in to reply