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. rhoCentralFoam求解shockTube算例出现数值震荡

rhoCentralFoam求解shockTube算例出现数值震荡

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

    大家好。我发现用rhoCentralFoam求解shockTube算例时x方向速度会有一些数值震荡。我调了网格数和maxCo数还是消不掉它。计算结果如下面两个图所示:
    f2.jpg f3.jpg
    根据这个网站[Shocktube – rhoCentralFoam TVD Schemes Test – curiosityFluids](Shocktube – rhoCentralFoam TVD Schemes Test – curiosityFluids)我将限制器改为vanAlbada可是仍有数值震荡。
    我将shockTube算例也上传了shockTube_test.rar 。我觉的应该哪里的参数配置需要调整一下,可以消除这个数值震荡。希望大佬们指点迷津。

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

    这个求解器有一个自己的算例,自带的算例没有震荡,你尝试一下

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

    1 条回复 最后回复
  • 李 离线
    李 离线
    李玉棠
    写于 最后由 编辑
    #3

    感谢李老师。这是我原始算例计算的结果。绿线是100网格./Allrun计算结果,蓝线是将网格改为800然后./Allrun的计算结果。其它参数都没有动。f5.jpg 这是算例shockTube_org.rar 我用的OF-4.1

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

    你时间步长减小没

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

    1 条回复 最后回复
  • 尚 在线
    尚 在线
    尚善若水
    写于 最后由 尚善若水 编辑
    #5

    时间步长太大了,0.02试试?还不行网格数改为1000

    1 条回复 最后回复
  • 李 离线
    李 离线
    李玉棠
    写于 最后由 编辑
    #6

    两位老师这是设置不同maxCo的计算结果。800网格
    f5.jpg
    调小时间步长后,我觉得振荡被抑制的很好,但没有完全消去。似乎rhoCentralFoam的鲁棒性不是很好?我印象中有些求解器的CFL数设置在0.98都没有数值振荡的。不知道我的理解对不对

    李东岳李 1 条回复 最后回复
  • 尚 在线
    尚 在线
    尚善若水
    写于 最后由 尚善若水 编辑
    #7

    rCF是这样的,库朗数一定要足够小,单纯加密网格没改善。这是我的对比计算结果
    image.png
    一些可能的分析(一家之言):
    image.png

    注:bF为blastFoam,rCF为rhoCentralFoam,sMF是我自己的求解器

    李 1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    在 中回复了 李玉棠 最后由 编辑
    #8

    这个求解器有一个自己的算例,自带的算例没有震荡,你尝试一下

    @李玉棠 我看原始算例,默认是1000个网格,也没有震荡。你尝试原始的算例没有

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

    李 1 条回复 最后回复
  • 李 离线
    李 离线
    李玉棠
    在 中回复了 李东岳 最后由 编辑
    #9

    @李东岳 李老师好。对于shockTube算例,OF-4.1和OF-6是100mesh的;OF-8是1000mesh。我用OF-4.1算了OF-8的算例,只添加了turbulenceProperties字典,计算结果如图:f9.jpg
    原始算例maxCo=0.2时Ux参数仍有震荡,调小maxCo后震荡被抑制。那位老师能给一下OF-8的计算结果吗?可能是版本的问题吗?

    1 条回复 最后回复
  • 李 离线
    李 离线
    李玉棠
    在 中回复了 尚善若水 最后由 编辑
    #10

    @尚善若水 谢谢老师。我觉得rCF这点不太好,,有没有更新的算法解决这个问题?现在有些既是高阶精度稳定性又好的格式,有能移植到OF的相关文献吗?

    尚 1 条回复 最后回复
  • 李东岳李 李东岳 被引用 于这个主题
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #11

    我用OF-4.1算了OF-8的算例,只添加了turbulenceProperties字典

    用of8添加湍流后,也有震荡么

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

    李 1 条回复 最后回复
  • 尚 在线
    尚 在线
    尚善若水
    在 中回复了 李玉棠 最后由 编辑
    #12

    @李玉棠 blastFoam好像就可以

    1 条回复 最后回复
  • 李 离线
    李 离线
    李玉棠
    在 中回复了 李东岳 最后由 编辑
    #13

    @李东岳 老师我OF-4.1版本算例constant文件夹下得有turbulenceProperties字典,里面写的是laminar,不涉及湍流。好像OF-8以后都不需要turbulenceProperties字典了。

    1 条回复 最后回复

  • 登录

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