Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 不建议开启`体积力面插值`算法

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

已定时 已固定 已锁定 已移动 OpenFOAM
4 帖子 2 发布者 4.6k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 李东岳 编辑
    #1

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

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

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

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    程 1 条回复 最后回复
  • 程 离线
    程 离线
    程迪
    在 中回复了 李东岳 最后由 编辑
    #2

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

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

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

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

    已婚,勿扰。
    本作品采用知识共享署名-非商业性使用-相同方式共享 4.0 国际许可协议进行许可。

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #3

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

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #4

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

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

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]