CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    bluecfd运算出错

    OpenFOAM
    3
    7
    1823
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • L
      lv1995 最后由 编辑

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

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

      1 条回复 最后回复 回复 引用
      • L
        lv1995 最后由 编辑

        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 ??:?
        
        1 条回复 最后回复 回复 引用
        • 李东岳
          李东岳 管理员 最后由 编辑

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

          线上CFD课程开始报名:http://www.dyfluid.com/class.html

          CFD高性能服务器 http://dyfluid.com/servers.html

          L 2 条回复 最后回复 回复 引用
          • L
            lv1995 @李东岳 最后由 编辑

            @东岳 还是一样的错误

            1 条回复 最后回复 回复 引用
            • L
              lv1995 @李东岳 最后由 编辑

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

              D 1 条回复 最后回复 回复 引用
              • D
                dzw05 教授 @lv1995 最后由 编辑

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

                自主匠心,普惠仿真。

                L 1 条回复 最后回复 回复 引用
                • L
                  lv1995 @dzw05 最后由 编辑

                  @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);
                  }
                  
                  1 条回复 最后回复 回复 引用
                  • First post
                    Last post