Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    用QMOM求解出现负的Abscissas

    Algorithm
    2
    12
    4522
    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.
    • Astremon
      Astremon last edited by

      用QMOM求解Pure aggregation中出现负的Abscissas,然后在公式(37)中就出现negative number cannot be raised to a fractional power 。有经验的大神怎么处理的呢。。
      0_1513428469642_f0215486-09e5-4319-abff-c830e8e43697-image.png

      Astremon 1 Reply Last reply Reply Quote
      • Astremon
        Astremon @Astremon last edited by

        @astremon 文献传送门

        1 Reply Last reply Reply Quote
        • 李东岳
          李东岳 管理员 last edited by 李东岳

          你是怎么发现的负值?无对流问题么?减小时间步长试试。

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

          Astremon 2 Replies Last reply Reply Quote
          • Astremon
            Astremon @李东岳 last edited by Astremon

            @李东岳 总是在一个时间点出现负的abscissas,提示错误“negative number cannot be raised to a fractional power ”。减小时间步长也不管用哟,,,。李老师这个是代码可以请您看看吗?
            可能是算法的局限性?

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

              @李东岳 没有对流。就是homogeneous 常数聚并核函数。

              1 Reply Last reply Reply Quote
              • 李东岳
                李东岳 管理员 last edited by

                73, 73行是你的源项么?

                temp_sum1[i2,i3] = w_i[i2]*w_i[i3]*(r_i[i2]**3 + r_i[i3]**3)**(m/3.)  ##
                temp_sum2[i2,i3] = w_i[i2]*w_i[i3]*r_i[i2]**m
                

                对应的方程是什么?

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

                Astremon 3 Replies Last reply Reply Quote
                • Astremon
                  Astremon @李东岳 last edited by

                  @李东岳 恩,是源项。是图中的公式(37)和(38)
                  0_1513603166441_7b8c545f-aa39-4bc2-aceb-db110a0d064d-image.png

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

                    @李东岳 其中$\beta (L,\lambda) = C_0 = 1 $

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

                      @李东岳 现在只能算到1.4s左右,超过1.4s出现负特征尺寸了。0_1513603616384_62a4f560-dffb-4bc2-843d-92e525eeac2e-image.png

                      1 Reply Last reply Reply Quote
                      • 李东岳
                        李东岳 管理员 last edited by

                        出现负值对应的m具体值是多少?

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

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

                          @李东岳

                          moments  --->[0.574993, 0.670135, 0.805031, 1, 1.29021, 1.73242]
                          
                          1 Reply Last reply Reply Quote
                          • 李东岳
                            李东岳 管理员 last edited by 李东岳

                            我试试,不过感觉你的算法没毛病,M3守恒。
                            你考虑过用Wheeler求节点么?比PD要强健。

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

                            1 Reply Last reply Reply Quote
                            • First post
                              Last post

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