1. Home
  2. OpenFOAM
  3. bluecfd运算出错

bluecfd运算出错

  • L

    使用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

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

  • L

    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删掉

    L 2 Replies
  • L

    @东岳 还是一样的错误

  • L

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

    D 1 Reply
  • D
    OpenFOAM教授

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

    L 1 Reply
  • L

    @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
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2