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. Stokes流中的圆球绕流的阻力系数,怎么设置才能算准。

Stokes流中的圆球绕流的阻力系数,怎么设置才能算准。

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

    @李东岳 东岳老师,稳态之前我用simpleFoam和icoFoam算过(但算的无量纲情况,得到的Cd值也是13.5左右)
    刚才做了个带量纲的simpleFoam,结果显示如下,最终Cd值基本持平在13.496左右;
    5cfa8bf2-6f97-4351-8246-388ef92d068f-image.png
    这里是带量纲的simpleFoam案例文件
    Re.1.MB.WMesh.AllSlip.Dimension.simpleFoam.zip

    当时没有选择继续用simpleFoam的原因主要是:

    • ①看到其他文献,比如庞明军老师团队的文献,用的瞬态求解器(他们的雷诺数稍大一些,Re最大到200,这个雷诺数范围内依然是符合轴对称流场结构)
      59cee156-df87-4fc6-8fac-78883448ce87-image.png

    • ②找了个wolfdynamics的培训视频中的圆柱绕流的simpleFoam案例来改,但算气泡的收敛判据我不知道怎么给才合适,只是看到simpleFoam一直跑,Cd值大概持平到13.5左右;主要责任是我自己也没往这里多想,因此之前主要都是用pimpleFoam在算。

    再次感谢东岳老师:xinxin: :xinxin2:

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

    我看了一下,调节了一下,确实有这个问题。
    我还试了一下Re=10的情况,无滑移没问题,slip就不对。
    我尝试了调试了一下,确实没效果。
    有点进坑的感觉 :136:

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

    Prometheus10P 1 条回复 最后回复
  • Prometheus10P 离线
    Prometheus10P 离线
    Prometheus10
    在 中回复了 李东岳 最后由 编辑
    #20

    @李东岳 我在办公室差点哭出声。
    12e28590-8149-4c59-8bdd-8dafb3f31267-image.png

    依然十分感谢东岳老师,我想再转到Fluent算一下看看情况,有新的进展我再来补充上。:xinxin:

    李东岳李 1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    在 中回复了 Prometheus10 最后由 编辑
    #21

    @Prometheus10 openfoam这面是没看出啥破绽。我捉摸能不能换个附加平均速度的周期性的网格试试。但是牵涉到画网格工作量比较大。不过你先试试fluent吧。我们再讨论 :-)

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

    Prometheus10P 1 条回复 最后回复
  • Prometheus10P 离线
    Prometheus10P 离线
    Prometheus10
    在 中回复了 李东岳 最后由 编辑
    #22

    @李东岳 好的,东岳老师,再次感谢!:xinxin:
    我正在画Fluent网格,有新进展时继续来补充。

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

    foamMeshToFluent可以把openfoam的网格转化成fluent,不过我不确定fluent能不能用wedge,我估计不能用,那你画吧..

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

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

    老铁,fluent那面后来怎么样

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

    Prometheus10P 1 条回复 最后回复
  • Prometheus10P 离线
    Prometheus10P 离线
    Prometheus10
    在 中回复了 李东岳 最后由 编辑
    #25

    @李东岳 东岳老师好,fluent我用Axissymmetric的2D算起来并不太顺利,收敛不了,我还在调试。

    这段时间主要在看相关文献。后续有进展我会更新上来。

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

    我不太懂你说的H-R理论解,那能不能是openfoam算的就是对的?那你看看文献吧,到时候再讨论

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

    1 条回复 最后回复
  • P 离线
    P 离线
    pengdi
    在 中回复了 Prometheus10 最后由 李东岳 编辑
    #27

    @Prometheus10 算了十秒还在33,我是不是应该放弃了,网格质量是不是不是特别好呀 感觉不应该呀,我之前用icoFoam算过雷诺数40的圆柱绕流,和文献对的还挺准

    1 条回复 最后回复
  • P 离线
    P 离线
    pengdi
    在 中回复了 Prometheus10 最后由 编辑
    #28

    @Prometheus10 看到你后面解决啦,可能定常流算不准最可能的问题就是网格问题

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

    Stokes定律Cd=24/Re,H-R理论解Cd=16/Re都是基于球形假设

    @pengdi 我也在等更新。总结一下,目前就是OpenFOAM计算的颗粒Cd趋向于理论解Cd=24/Re,气泡不趋向于H-R理论解Cd=16/Re。

    我看OpenFOAM计算的没有问题。但是因为我不太熟悉这个H-R理论解,我也在想会不会是H-R理论解的事?

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

    P 1 条回复 最后回复
  • P 离线
    P 离线
    pengdi
    在 中回复了 李东岳 最后由 编辑
    #30

    @李东岳 东岳老师,我仔细看了一下,个人感觉理论解应该没问题。我看了一下算例设置,感觉有可能出问题的地方有两条:
    1)边界条件:气泡压力边界设置为零梯度是不是不太合理?我感觉气泡内是不是应该设置为常压?
    2)网格及时间步长。这个出问题的可能性比较小,可以试一下网格和时间步的无关性分析。

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

    2)网格及时间步长。这个出问题的可能性比较小,可以试一下网格和时间步的无关性分析。

    这个我试了一下加密网格,网格还是他那个网格,没啥影响。我看他说他尝试了3D算例也是这样

    1)边界条件:气泡压力边界设置为零梯度是不是不太合理?我感觉气泡内是不是应该设置为常压?

    这个感觉有点不太套路呢。不知道会怎么样。我感觉应该不是这个问题。但又没有啥方向,你试试?

    你也做这方面研究么

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

    P 1 条回复 最后回复
  • P 离线
    P 离线
    pengdi
    在 中回复了 李东岳 最后由 李东岳 编辑
    #32

    @李东岳 我试了一下刚才两种方法,都不太行。

    我用symmetry边界算了一下,和理论解符合的很不错。可能slip边界还有一些问题。

    83f93549-acef-4d01-839f-bccadb281ec2-image.png

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

    厉害啊 你计算算一下 看看能稳定到多少 看起来可能是bug了 symmetry是物理限定性条件 限定性更强 slip只是速度的条件

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

    P 2 条回复 最后回复
  • P 离线
    P 离线
    pengdi
    在 中回复了 李东岳 最后由 编辑
    #34

    @李东岳 好的,东岳老师。我的内心是崩溃的,我有一篇正在写的论文用的slip边界已经算了100多个工况了

    1 条回复 最后回复
  • P 离线
    P 离线
    pengdi
    在 中回复了 李东岳 最后由 李东岳 编辑
    #35

    @李东岳 安慰一下自己,有可能是symmetry边界对法向速度的影响太强了

    @Prometheus10 @李东岳 稳定在17.78了 感觉符合的挺好的

    89a714ed-20c5-4ff8-a74d-fb0fd831ef99-image.png

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

    妥了,实锤了,就是slip与symmetry的问题,估计slip没有做到完全的无剪切。我这几天详细看一下

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

    P 1 条回复 最后回复
  • P 离线
    P 离线
    pengdi
    在 中回复了 李东岳 最后由 编辑
    #37

    @李东岳 好,期待东岳老师的深度指导

    Prometheus10P 1 条回复 最后回复

  • 登录

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