CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    关于网格局部加密问题

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

      在ICEM中局部加密使用了
      f00a8079-1bcb-419f-af18-f7e2376c299d-image.png
      这样生成的网格并不是点对点,而是一边对应多点
      323eb693-2224-4f09-a768-0cbceeb349b3-image.png
      在OpenFOAM中的paraview中显示效果
      c00320e5-5893-4a47-b0ed-87b5ec8aaf97-image.png
      请教下,这样的网格,在OpenFOAM计算的时候,是否可行呢?
      谢谢~~

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        是可行的,
        但是可能会引起很大误差,通常建议1对2,逐步过渡

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

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

          @东岳 谢谢老师~~

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

            @东岳
            老师,已经进行了对网格的1对2改进:
            1592640955(1).jpg
            在ICEM转化到OpenFOAM中时,在ployMesh文件夹中出现了set文件夹里面是一个nonOrthoFaces文件

            2f9e8f66-ef5b-49e9-8a76-ab6507f8729f-image.png
            出现这种文件,是由于1对2网格导致的,和下图这种网格非正交的是不同的
            a262e836-9e44-4399-8255-2b270336b97e-image.png

            所以在我这个案例的fvSolution文件中

            PIMPLE 
            {
                momentumPredictor   no;
                nOuterCorrectors    1;
                nCorrectors         2;
                nNonOrthogonalCorrectors 0;
                maxCo               10;
                maxAlphaCo          5;
                rDeltaTSmoothingCoeff 0.05;
                rDeltaTDampingCoeff 0.5;
                nAlphaSpreadIter    0;
                nAlphaSweepIter     0;
                maxDeltaT           1;
            }
            

            是否需要把nNonOrthogonalCorrectors 0改为nNonOrthogonalCorrectors 1,或者其它地方需要进行改动呢?

            1 条回复 最后回复 回复 引用
            • 李东岳
              李东岳 管理员 最后由 编辑

              这不好说,个人认为不需要改为1

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

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

                @东岳 谢谢老师

                1 条回复 最后回复 回复 引用
                • W
                  WangGH 最后由 编辑

                  周老师,这里我测试发现好像这种结构网格的加密方式对于三维模型可以,对于二维网格,在经过局部加密后再导入到OF里面是会报错的,错误是:--> FOAM FATAL ERROR:
                  face 11528151 does not belong to cell 210100. Face owner: -1 neighbour: 210099

                  From Foam::cellShape Foam::extrudedQuadCellShape(Foam::label, const labelList&, const faceList&, const labelList&, const labelList&, Foam::label, Foam::faceList&)
                  in file extrudedQuadCellShape.C at line 105.
                  

                  FOAM aborting
                  感觉好像是在z方向拉伸出来的面没有识别到具体的cell中

                  1 条回复 最后回复 回复 引用
                  • First post
                    Last post