Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    reactingFoam计算速度的问题

    OpenFOAM
    2
    6
    663
    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.
    • W
      Wust.liao last edited by Wust.liao

      OpenFOAM版本:8
      求解器:reactingFoam
      发现一个很头疼的问题,在使用reactingFoam计算甲烷-空气预混燃烧时,计算的速度似乎跟我设置的当量比有很大的关系,当我把当量比设为0.6时,计算会越算越慢,从最开始1s能算几个步长到后面10几秒才能算一个步长,当计算工况到10s以后基本就卡住不动了。但是当量比设为1的话计算就很正常,每个步长计算的时间基本都是一样。
      另外一个就是关于constant/chemistryProperties文件里的这个设置:equilibriumRateLimiter on/off,(图1)我没有找到关于这个开关的资料,但是代码说明里的建议是保持关闭。但是当我把这个开关打开时发现计算速度就会变得很快,不会随当量比发生明显的变化。只是计算的结果看上去很难受,火焰锋面处的体积放热率、组分浓度是锯齿状的(图2)。而正常情况因是图3的情况。

      d63462a7-247b-49f4-8e66-d2b8ce43e35c-图片.png
      图1

      eca3fdd1-6575-406b-aa01-96e8cfa5447b-图片.png
      图2

      3962fe71-0428-4710-b89f-c7dd67a00117-图片.png
      图3

      香柏树 1 Reply Last reply Reply Quote
      • W
        Wust.liao last edited by

        目前还没有找到关于计算速度变慢的原因,但是经过这段时间的调试发现了一些规律并尝试着把问题解决了。
        规律:当计算变慢时火焰峰值温度都在2500k左右,而实际上合理的温度应该在1800附近。只要当计算随着迭代到达了2500k速度马上就降下来了,最后基本就算不动了。
        解决办法:调小时间步长,比如用0.1的库朗数去计算,温度就不会到这么高,随着迭代温度回到1800k后可以适当增大库朗数。

        1 Reply Last reply Reply Quote
        • 香柏树
          香柏树 @Wust.liao last edited by

          @wust-liao 你好,同学。我刚开始用reactingFoam计算甲烷预混燃烧,一直出现mpirun报错,请问你遇到过么,知道怎么解决么?万分感谢~
          32f091c0-95ea-41dc-bde2-0e0b3163fa28-图片.png
          7fdad4cb-f8e2-41d0-8ab0-8615ba8cf784-图片.png

          W 1 Reply Last reply Reply Quote
          • W
            Wust.liao @香柏树 last edited by

            @香柏树 你有试过自带的案例吗,先排除软件自身的问题

            香柏树 1 Reply Last reply Reply Quote
            • 香柏树
              香柏树 @Wust.liao last edited by

              @wust-liao 自带sandia的案例算了没问题,我是把sandia的rans改为les,换了一个复杂的燃烧室计算的。不知道哪块出问题了,导致只能断断续续地计算

              W 1 Reply Last reply Reply Quote
              • W
                Wust.liao @香柏树 last edited by

                @香柏树看上去就是发散了,换个简单的燃烧室算算,模型一个一个往上加。不容易一步到位的

                1 Reply Last reply Reply Quote
                • First post
                  Last post

                CFD中文网 | 东岳流体 | 京ICP备15017992号-2