CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    SIMPLE 算法中密度如何求解?

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

      SIMPLE 算法中连续方程被改为压力修正方程来求解,请问密度改怎么办?不可压缩流体

      1 条回复 最后回复 回复 引用
      • 赵
        赵一铭 最后由 编辑

        如果是不可压缩,为什么要求解密度?

        谷 李东岳 2 条回复 最后回复 回复 引用
        • 谷
          谷柏辰 @赵一铭 最后由 编辑

          @赵一铭 虽然是不可压缩,但其中涉及很多反应,气体组分会大大改变,所以密度也会变,当下的处理方程是用理想气体状态方程,rho=MP/RT 但结果并不理想,是否有其他更新密度手段?

          1 条回复 最后回复 回复 引用
          • C
            CFD中文网 最后由 编辑

            0_1492182483521_草图.png

            是否是这个图的意思?The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM® and Matlab®,page 666。一般这种用来求解大马赫数问题,或者那种那个所谓的全速流动。

            CFD中国标准用户测试帐号
            目前由徐笑笑登录

            1 条回复 最后回复 回复 引用
            • 谷
              谷柏辰 最后由 编辑

              非常非常感谢,我的算法流程图跟这个基本类似,不过密度的更新我是的理想其他状态方程,因为马赫数很低,假定为不可压缩流体,我的程序也是用matlab写的,让我百思不得其解的是我的程序都是按照标准方程写的,也能启动运气,但随着计算的进行,速度方程会首先不满足boundedness条件,求解的速度值超大,然后就是整个程序的不收敛。我现在是怀疑我的密度的处理有问题,所以才有此问,不知道该往哪个方向debug:upset: 0_1492187127917_QQ截图20170414182045.jpg!

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

                @赵一铭 不可压缩一般指的是密度针对压力的不可压缩,但是可能会因为温度的变化而改变,比如HVAC中的浮力驱动流。 @谷柏辰 很遗憾至今你的问题还没有解决,不过错位网格+matlab编程+数值不准的问题很难判断到底是哪里错了。 @cfd-china 那个书中的方法应该是错位网格所用的方法。

                只能祝好运:big_mouth: ...

                CFD高性能服务器 http://dyfluid.com/servers.html
                2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

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

                  @李东岳已经非常感谢了,我再继续努努力吧!:expressionless:

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