CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    VOF方程的疑问

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

      各位老师同学大家好,最近遇到了点小问题,真心求帮助。

      假设流体不可压缩。MULES算法中,添加了人工压缩项以保证界面尖锐:
      1.jpg

      按我的理解,其实就相当于在VOF方程中引入了一个源项S
      68236bea-8edc-485a-83ca-225a76b1c4cc-image.png

      虽然界面压缩项只在自由表面附近起作用,但是严格来讲,连续性方程应该修改成如下形式才对吧?(因为对于不可压缩流体,相体积分数的控制方程与连续性方程是一致的)
      493abdf4-2eee-494a-9f07-1a2adeca4654-image.png

      但是在interfoam里面并未发现类似的代码,请问这是怎么回事呢?是对于这个问题,这个源项可以忽略吗?

      推广一下,如果想在VOF方程中添加其他的源项,是不是也要修改质量守恒方程?

      此外,对于相体积分数为0的网格,上式中除以相体积分数这个操作是不是需要特殊处理一下?我试了一下,假设这样
      fc008287-68c8-4b11-9798-9b3081be73df-image.png
      程序运行倒是没什么问题,但是这个意义我没搞清楚,这么一除不是得到无限大的速度散度了吗?目前,我倾向于这么做:对于相体积分数为0的网格,我直接将源项S的值设为0,即源项不起作用。

      说的可能有点乱,但主要的意思就是这样。。。

      谢谢大家!

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

        额 好像发现个问题,质量守恒方程应该是这个
        7973d26e-efef-4af5-9e55-686a447c4996-image.png

        其中𝑆(𝜌)与 𝑆(𝛾)形式一样,只不过自变量由相体积分数变成了密度,不知道这样对不对?

        感觉有点弄乱了,求大佬指点,不胜感激:135:

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