Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    有遇见过simpleFoam达不到稳态的情况吗?

    OpenFOAM
    2
    3
    1491
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • chpjz0391
      chpjz0391 last edited by 李东岳

      最近再用snappyHexMesh画网格,网格画完了checkMesh也通过了,但是simpleFoam达不到稳态。用了几个probe来检查,发现速度一直在震荡。 同样的system下面的设置,之前跑其他的算例都没问题的。
      我贴一下我的fvSchemes

      
      /*--------------------------------*- C++ -*----------------------------------*\
      | =========                 |                                                 |
      | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
      |  \\    /   O peration     | Version:  2.3.1                                 |
      |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
      |    \\/     M anipulation  |                                                 |
      \*---------------------------------------------------------------------------*/
      FoamFile
      {
          version     2.0;
          format      ascii;
          class       dictionary;
          location    "system";
          object      fvSchemes;
      }
      // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
      
      ddtSchemes
      {
          default         steadyState;
      }
      
      gradSchemes
      {
          default         cellLimited leastSquares 1; //Gauss linear;
      }
      
      divSchemes
      {
          default         none;
          div(phi,U)      bounded Gauss linearUpwindV cellLimited leastSquares 1;
          div(phi,k)      bounded Gauss linearUpwind celllimited leastSquares 1;
          div(phi,epsilon) bounded Gauss linearUpwind celllimited leastSquares 1;
          div(phi,R)      bounded Gauss linear corrected;
          div(R)          Gauss linear;
          div(phi,nuTilda) bounded Gauss upwind;
          div((nuEff*dev2(T(grad(U))))) Gauss linear;
      }
      
      laplacianSchemes
      {
          default         Gauss linear corrected;
      }
      
      interpolationSchemes
      {
          default         linear;
      }
      
      snGradSchemes
      {
          default         corrected;
      }
      
      fluxRequired
      {
          default         no;
          p               ;
      }
      
      
      // ************************************************************************* //
      
      1 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

        simpleFoam达不到稳态的情况吗?

        发现速度一直在震荡

        这是可能的。有些物理过程,用稳态模拟就是会一直震荡。不过,用snappyHexMesh做的网格,也有可能是网格的问题,snappyHexMesh做出来的网格某些部位质量很差。但如果几何不是特别复杂,一般不会影响到全局。

        另外,如果你要降低震荡,应该增加粘性,用低阶格式,比如迎风试试。

        线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        chpjz0391 1 Reply Last reply Reply Quote
        • chpjz0391
          chpjz0391 @李东岳 last edited by

          @东岳 好的谢谢东岳,我先用一阶迎风格式试试

          1 Reply Last reply Reply Quote
          • First post
            Last post

          CFD中文网 | 东岳流体 | 京ICP备15017992号-2
          论坛登录问题反馈可联系 li.dy@dyfluid.com