Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    VOF方程的疑问

    OpenFOAM
    1
    2
    242
    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.
    • C
      CFDngu last edited by

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

      假设流体不可压缩。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 Reply Last reply Reply Quote
      • C
        CFDngu @CFDngu last edited by

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

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

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

        1 Reply Last reply Reply Quote
        • First post
          Last post

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