Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    不建议开启`体积力面插值`算法

    OpenFOAM
    2
    4
    2179
    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.
    • 李东岳
      李东岳 管理员 last edited by 李东岳

      目前OpenFOAM基金会,ESI-OpenCFD(OpenFOAMv1706)下的体积力面插值代码提出的本源在于使得较大体积力下的计算结果无震荡且稳定。然而,目前这个算法并不能够提供合理的数据,和实验甚至通常的算法偏离极大。个人觉得主要有以下缺陷:

      1. 速度处理过于激进:在UEqn.H中将时间导数项移动到压力方程,虽然从求解方程的角度是合理的。但是目前这种方法是否是数据偏离的原因暂不明确;
      2. 力的处理任意性:力是不具有守恒特性的,因此基金会仅仅将力插值到面上来处理,再一次,虽然从求解方程的角度是合理的。但是目前这种方法是否是数据偏离的原因暂不明确;

      目前这个算法需要大量的验证和提高。同时,合理的计算体积力需要更好的算法。

      CFD高性能服务器 http://dyfluid.com/servers.html

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

        @李东岳 在 不建议开启`体积力面插值`算法 中说:

        力的处理任意性:力是不具有守恒特性的,因此基金会仅仅将力插值到面上来处理,再一次,虽然从求解方程的角度是合理的。但是目前这种方法是否是数据偏离的原因暂不明确;

        保守力是可以移到压力方程中的,整到面上没啥不对。

        估计处理的时候体积力的旋度没整成0吧。可以测试一下体积力的旋度的量级。

        github: chengdi123000
        网站:chengdi123000.github.io
        本作品采用知识共享署名-非商业性使用-相同方式共享 4.0 国际许可协议进行许可。

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

          某些体积力的旋度必然不为0.

          CFD高性能服务器 http://dyfluid.com/servers.html

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

            对于单气泡的曳力,升力更倾向于是一种面力。但是由于这个面不是网格单元的面。处理起来存在一定困境。我们把这些力当做体积力?有一些模拟结果表明存在震荡。插值做面力(类似基金会这种处理)?但是结果是错的。

            需要深入研究力的处理方式。

            CFD高性能服务器 http://dyfluid.com/servers.html

            1 Reply Last reply Reply Quote
            • First post
              Last post

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