CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    block的划分对质量分数是否有影响

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

      各位前辈好,有一个问题困扰我很久了。计算中途总是会报错 “质量分数总和大于1”,无论把来流设成空气(氧气,氮气,水)还是单组分气体(纯氧),都会出现这个报错,并且报错的都是三角形的上下两个block, 同样的情况下如果我关闭粘性和湍流模型则能完成计算。

      我已经重新划分了block开始计算,目前已经算完了一半,暂时还没有报错,不知道大家有没有遇到过类似的问题?

      然后还要谢谢东岳老师,观摩了这么久,终于可以来这里发帖了:mihu:

      模型如图所示
      model.png

      报错的block
      错误的块.png

      计算停止的地方
      1111.png

      报错信息
      2222_pro.png

      高超声速燃烧

      liujunCFD 1 条回复 最后回复 回复 引用
      • liujunCFD
        liujunCFD @liujunCFD 最后由 编辑

        @liujunCFD 对了,可能有的前辈会说三角形上下两个block的节点和相邻的点没有对上,这是因为靠近三角形的地方是边界层,另外就算对上了,也是会报同样的错,这个我已经试过了

        高超声速燃烧

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

          没见过你这个报错,你用rhoCentralFoam或者rhoSimpleFoam试试报错么

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

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

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

            @李东岳 老师好!我重新划分了block 就能完成计算了:jingya:

            话说我的计算结果还是不理想(下图均为密度)

            刚开始看着还挺正常的
            Snipaste_2020-11-12_10-41-03.png

            但是随着计算进行,激波不断向前移动
            Snipaste_2020-11-12_10-42-37_pro.png

            而且边界层好像被激波拉起来了?

            壁面条件用的都是绝热壁面(考虑粘性),现在正在尝试壁面设成固定温度或者所有壁面不考虑粘性,不知道能不能改善结果

            这个是原始文献的图
            Snipaste_2020-11-12_11-07-47.png

            话说第一次见到李老师真人,激动:chitang:

            高超声速燃烧

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

              你的激波看起来被耗散掉了。感觉是数值格式的问题。你用的什么高阶格式。也可以尝试去掉粘度看看会不会更尖锐。

              另外你这个为什么不做结构网格,为什么要用交界面?

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

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

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

                @李东岳 老师好 我用的这个求解器虽然基于OF 不过还是有些差别 我用的是该求解器自带的网格处理工具,做的是结构网格,由“点→线→面”生成不同的block,然后block组合成整个几何,您说的交界面我不太理解是什么意思。网格代码比较长,如果老师对该求解器划分网格的方式感兴趣的话我可以把网格代码发上来。

                这是我的一些参数选择 双破折号后面的为注释
                代码.png

                另外像您说的,壁面去掉粘性后,计算结果正常多了
                无粘性密度-tuya-tuya.png

                忽略左下角哈哈哈 那个block的壁面我忘记去掉粘性了:chouchou:

                高超声速燃烧

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