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
16 帖子 2 发布者 3.1k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #1

    c51cbb80-1941-44d6-a346-6c591165f0f1-b12210c4b7a572099c5d0aa1728cf5f.png 75372991-556d-4a9f-bcd8-fca775180ae5-37ff47f32b1473c7d85be81d8dfdcba.pngcc090fa1-4fde-4573-a677-0070cf3a874d-35ac3b4bded1e5fe45bb3dc84a17a9d.png 各位前辈我植入这个hosen曳力模型之后编译没有报错, 但是运行算例时候马上就会报错,这种原因一般是怎么造成的呢,对于RE数等为0的情况已经做了限制

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #2

    image.png看了一下这个表达式的空间,雷诺数限制在5以上的话基本上CD小于10,是合理的范围,但是设置residualRe为5之后还是发散

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

    你这个不是算法的问题。是植入的问题。感觉是没有植入对。bad size一般是数组超过界限了。

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

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #4

    谢谢李老师回复,我在调整之后能够跑算例了,但是跑一段时间后还是发散,在求解域上存在一些动能比较大的点,这种情况一般怎么去解决比较好呢,是直接限制Cd的大小还是有更好的办法呢

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

    一般曳力模型不应该植入不正确。你看那个Cd-Re那个图,天下曳力模型基本都没有太大差异。在这种情况下,可能你写了bug进去。

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

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #6

    李老师我又整理了一下我的问题,首先现在openfoam目前在内的曳力模型都是针对牛顿流体的,然后一般在雷诺数较大的情况下,曳力系数都是随着雷诺数增加,比如说常用的:
    image.png
    image.png等等。我目前计算这个是非牛顿的体系,然后非牛顿体系内的气泡曳力系数在一定范围内是随着雷诺数降低的,现存的关联式也基本上都是符合这个规律:
    image.png比如:
    image.png等等。然后这些实验一般针对的雷诺数范围都在几百以下,这样的话使用这些公式一方面在雷诺数较小的时候cd很大,我直接把最小雷诺数设置的比较大,一般牛顿的设置为1e-3,我设置成个位数,在雷诺数较大的入口等地方就会出现cd小的可怜的情况,就会导致计算发散。然后我看很多文献计算非牛顿系统并没有使用牛顿的曳力系数模型,或者使用一个下界来限制高雷诺数下的cd,dewsbury这个文章也提到了他认为雷诺数比较大的时候cd趋近0.95,
    image.png尝试类似设置效果并不好还是会发散。目前想问一下李老师,一是这种限制上下界是否合理,二是如果不存在高雷诺数下靠谱的公式,高雷诺数下再使用牛顿的公式进行融合是否合理呢。
    image.png

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

    你的算例,非牛顿流体,粘度场怎么样,最大值能有多少

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

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #8

    粘度最大280mpas,用幂律定律描述的话流变指数n0.8。从我实验测试的气泡速度来看剪切变稀的粘度最低能到160多,仿真的粘度场也差不多这个范围

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

    你用SN曳力模型,针对这个非牛顿流体,能计算正确的流场结果么

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

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #10

    SN模型和ishiizuber模型算出来大差不差,感觉误差算是可接受
    untitled.png

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

    看起来还是曳力模型的问题。你可以把你的模型Cd跟Re的线,以及SN的Cd跟Re的线,放到一个图里么?

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

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #12

    Fig1.png
    Fig2.png,我画的ishiizuber是看HZDR他们baselinemodel那个模型,表达式如下
    image.pngopenfoam自带的ishiizuber趋势应该是和sn模型差不多

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

    在Re比较小的时候,Cd差的很大。在算例初始的时候,速度的滑移非常小,一般Re非常小。所以你的模型可能会有稳定性问题。我建议你把你的模型,在Re比较小的时候做个clip,而不是原生的模型,来提高数值稳定性。

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

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #14

    好的李老师我去试试

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

    我在2018年左右的时候研究过曳力对稳定性的影响,当时做了一个数值稳定性的改进,但具体怎么搞的我已经记不清了。clip是个非常好处理的方法。但是我记得Re比较小的时候存在一个问题。当时用拉格朗日粒子验证的。

    这个是需要写sci的么,还是你做着玩的

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

    1 条回复 最后回复
  • WeakFormW 离线
    WeakFormW 离线
    WeakForm
    写于 最后由 编辑
    #16

    是的老师,打算评估一下各种非牛顿曳力模型作为小论文一部分内容

    1 条回复 最后回复

  • 登录

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