CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    rhoPimpleFoam计算飞机跑几步就发散

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

      各位大佬,我参考文献 3rd High-Lift Workshop Summary Paper - OpenFOAM, STAR-CCM+ & LAVA simulations on Unstructured Grids链接https://arc.aiaa.org/doi/10.2514/6.2018-1253。用rhoPimpleFoam计算飞机,但是算了几步就发散了,希望各位大佬帮忙看看是fvSchemes和fvSolution的问题嘛,还是nut、nuTilda设置的问题。

      目前用的是作者上传到会议网站的网格,按照文章里说的,采用的是SA模型,以及动量和湍流项采用的是二阶迎风。

      Courant Number mean: 1.4881475e-08 max: 1.8237221
      Time = 3.2e-09
      
      
      * 
      * PIMPLE: Iteration 1
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      DILUPBiCGStab:  Solving for Ux, Initial residual = 2.421018e-11, Final residual = 2.421018e-11, No Iterations 0
      DILUPBiCGStab:  Solving for Uy, Initial residual = 2.1538356e-10, Final residual = 2.1538356e-10, No Iterations 0
      DILUPBiCGStab:  Solving for Uz, Initial residual = 2.0979924e-10, Final residual = 2.0979924e-10, No Iterations 0
      DILUPBiCGStab:  Solving for e, Initial residual = 4.1760218e-08, Final residual = 2.1963671e-13, No Iterations 1
      GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
      GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
      GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
      GAMG:  Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      time step continuity errors : sum local = 5.0250634e-11, global = 5.0210593e-11, cumulative = 3.7505015e-09
      GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
      GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
      GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
      GAMG:  Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      time step continuity errors : sum local = 5.0250633e-11, global = 5.0210593e-11, cumulative = 3.8007121e-09
      PIMPLE: Iteration 2
      DILUPBiCGStab:  Solving for Ux, Initial residual = 1.3321419e-08, Final residual = 1.002751e-14, No Iterations 1
      DILUPBiCGStab:  Solving for Uy, Initial residual = 6.2117311e-09, Final residual = 6.2117311e-09, No Iterations 0
      DILUPBiCGStab:  Solving for Uz, Initial residual = 1.4225061e-10, Final residual = 1.4225061e-10, No Iterations 0
      DILUPBiCGStab:  Solving for e, Initial residual = 1.5129023e-07, Final residual = 2.4037934e-13, No Iterations 1
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 3.900971e-09
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      GAMG:  Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 4.00123e-09
      PIMPLE: Iteration 3
      DILUPBiCGStab:  Solving for Ux, Initial residual = 2.4133e-08, Final residual = 2.8587861e-14, No Iterations 1
      DILUPBiCGStab:  Solving for Uy, Initial residual = 1.1420948e-08, Final residual = 1.9949995e-14, No Iterations 1
      DILUPBiCGStab:  Solving for Uz, Initial residual = 1.4399722e-10, Final residual = 1.4399722e-10, No Iterations 0
      DILUPBiCGStab:  Solving for e, Initial residual = 1.3956424e-07, Final residual = 2.3199397e-13, No Iterations 1
      [73] #0  Foam::error::printStack(Foam::Ostream&)sh: addr2line: command not found
       addr2line failed
      [73] #1  Foam::sigFpe::sigHandler(int)sh: addr2line: command not found
       addr2line failed
      [73] #2  ?sh: addr2line: command not found
       addr2line failed
      [73] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate()sh: addr2line: command not found
       addr2line failed
      [73] #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct()sh: addr2line: command not found
       addr2line failed
      [73] #5  ?sh: addr2line: command not found
       addr2line failed
      [73] #6  __libc_start_mainsh: addr2line: command not found
       addr2line failed
      [73] #7  ?sh: addr2line: command not found
       addr2line failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x3c75cd0, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=22, req_array=0x3966230, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=52, req_array=0x304c450, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=46, req_array=0x355b680, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=36, req_array=0x3dfabb0, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x349e540, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=26, req_array=0x3bd5ce0, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=12, req_array=0x3f18790, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=18, req_array=0x3f97b30, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=14, req_array=0x345ac20, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=32, req_array=0x3df8960, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x4316400, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x2fe0f30, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=50, req_array=0x3a36c60, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x3f6c5f0, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      Fatal error in PMPI_Waitall: Unknown error class, error stack:
      PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x520edf0, status_array=0x1) failed
      MPIR_Waitall_impl(181)................: fail failed
      MPIDU_Complete_posted_with_error(1710): Process failed
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
      slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
      
      

      fvSchmes

      // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
      
      ddtSchemes
      {
          default         Euler;
      }
      
      gradSchemes
      {
          default         Gauss linear;
      
          limited         cellLimited Gauss linear 1;
          grad(U)         $limited;
          grad(k)         $limited;
          grad(omega)     $limited;
      }
      
      divSchemes
      {
          default         none;
      
          div(phi,U)      Gauss linearUpwind grad(U);//Gauss linearUpwind limited;
      
          turbulence      Gauss linearUpwind default;//Gauss linearUpwind limited;
          energy          Gauss linearUpwind default;//Gauss linearUpwind limited;
      
          div(phi,k)      $turbulence;
          div(phi,omega)  $turbulence;
      
          div(phi,e)      $energy;
          div(phi,K)      $energy;
          div(phi,Ekp)    $energy;
          div(phi,nuTilda)  $energy;
      
          div(phiv,p)     Gauss upwind;
          div((phi|interpolate(rho)),p) Gauss upwind;
      
          div(((rho*nuEff)*dev2(T(grad(U)))))    Gauss linear;
      }
      
      laplacianSchemes
      {
          default         Gauss linear limited 0.5;//Gauss linear corrected;
      }
      
      interpolationSchemes
      {
          default         linear;
      }
      
      snGradSchemes
      {
          default         limited 0.5;//corrected;
      }
      
      wallDist
      {
          method meshWave;
      }
      
      

      fvSolution

      
      // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
      
      solvers
      {
          p
          {
              solver          GAMG;
              smoother        GaussSeidel;
              tolerance       1e-5;
              relTol          0.01;
      		nPreSweeps      0;
      		nPostSweeps     2;
      		cacheAgglomeration  on;
      		agglomerator      faceAreaPair;
      		ncellsInCoarsetLevel      1000;
      		mergeLevels          1;
          }
      
          pFinal
          {
              $p;
              relTol          0;
          }
      
          "(rho|U|k|omega|e)"
          {
              solver          PBiCGStab;
              preconditioner  DILU;
              tolerance       1e-8;
              relTol          0.0;
          }
      
          "(rho|U|k|omega|e|nuTilda)Final"
          {
              $U;
              relTol          0;
          }
      }
      
      SIMPLE
      {
          residualControl
          {
              p               1e-5;
              U               1e-5;
              "(k|omega|e)"   1e-5;
          }
      
          nNonOrthogonalCorrectors 3;
          pMinFactor      0.1;
          pMaxFactor      2;
      }
      
      PIMPLE
      {
          nCorrectors              2;
          nNonOrthogonalCorrectors 3;
          nOuterCorrectors         3;
          pMinFactor      0.1;
          pMaxFactor      2;
      }
      
      relaxationFactors
      {
      	p       0.3;
      	U       0.7;
      	k       0.7;
      	omega   0.7;
       
      }
      
      // ************************************************************************* //
      
      

      nut

      
      // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
      
      dimensions      [0 2 -1 0 0 0 0];
      
      internalField   uniform 0;
      
      boundaryField
      {
          freestream
          {
              type            calculated;
              value           uniform 0;
          }
      
          wall
          {
              type            nutkWallFunction;
              type            nutkWallFunction;
              value           uniform 0;
          }
      
          #includeEtc "caseDicts/setConstraintTypes"
      }
      
      // ************************************************************************* //
      
      

      nuTilda

      
      // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
      
      dimensions      [0 2 -1 0 0 0 0];
      
      internalField   uniform 0;
      
      boundaryField
      {
          freestream
          {
              type            inletOutlet;
              inletValue      $internalField;
              value           $internalField;
          }
      
          wall
          {
              type            fixedValue;
              value           uniform 0;
          }
      
          symmetry
          {
              type            symmetry;
          }
      
          #includeEtc "caseDicts/setConstraintTypes"
      }
      
      // ************************************************************************* //
      
      
      1 条回复 最后回复 回复 引用
      • First post
        Last post