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. 动网格和overset物体速度快速震荡导致崩溃的问题

动网格和overset物体速度快速震荡导致崩溃的问题

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

    @tower 老哥,不知道你找到问题原因没有?:135: 求分享!我是模拟波物相互作用问题的,对你遇到的这个问题很好奇o_o

    towerT 1 条回复 最后回复
  • towerT 离线
    towerT 离线
    tower
    在 中回复了 xiezhuoyu 最后由 编辑
    #3

    @xiezhuoyu 没有:135: 如果解决了 会在这里贴上结果的

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

    @tower 老哥,我也遇到了同样的问题:136:
    用overset,浮体z方向受力正负交替变化,然后越来越大。
    我看结果文件是因为浮体底部的压力变化剧烈,忽正忽负。所以想着加密一下浮体附近的网格,试试能不能解决问题。
    可是!我加密网格后,浮体壁面的Y+小于30了(浮体顶部(接触空气的部分)尤其小),其实我的网格也不是很密,网格尺寸大概就是浮体吃水的1/8,而且如果可能的话,我想用壁面函数啊,请问,老哥你遇到这种问题一般如何解决?谢谢
    也希望有其他前辈能够指教!谢谢

    towerT 1 条回复 最后回复
  • towerT 离线
    towerT 离线
    tower
    在 中回复了 xiezhuoyu 最后由 编辑
    #5

    @xiezhuoyu 我也询问过别人,给的建议也是加密网格,但是我试着加密了下,还是相同的错误。 其实我不是很理解的是,从我po在上面的log文件来看压力虽然迭代的次数多,但是最终残差不是已经满足条件了吗,为何压力积分后的力还会不准确,出现这种震荡问题呢

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

    最近在做一个波浪中物体二自由度的运动案例,最开始使用普通的动网格来计算的,因为计算幅度可以接受,但是结果却发现在开始计算0.2s左右,计算就发散了,发散的原因是物体Z方向上的速度急剧增大,

    这个算例要是小的话,可以发上来给你看看

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

    towerT 1 条回复 最后回复
  • towerT 离线
    towerT 离线
    tower
    在 中回复了 李东岳 最后由 编辑
    #7

    多谢@东岳老师,我把动网格和overset的案例文件都上传了,您看一下。
    动网格我是用的waves2foam+v1712的组合 overset用的1906的版本
    overset: oversetMesh.zip
    动网格: dynamicMesh.zip
    :baobao:

    1 条回复 最后回复
  • towerT 离线
    towerT 离线
    tower
    写于 最后由 李东岳 编辑
    #8

    各位,
    我尝试增大PIMPLE里面循环次数来增加p_rgh的迭代次数,但是发现每次迭代的次数还是很多 不是一个递减的趋势
    如下:

    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.72710724614e-07, Final residual = 8.72514217379e-09, No Iterations 382
    time step continuity errors : sum local = 7.03336340027e-09, global = 1.9150124018e-09, cumulative = 3.04929773866e-07
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.69104939702e-07, Final residual = 1.44876414351e-09, No Iterations 335
    time step continuity errors : sum local = 1.11124908301e-09, global = 1.95050292461e-10, cumulative = 3.05124824159e-07
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.2782727071e-07, Final residual = 1.06023228234e-09, No Iterations 607
    time step continuity errors : sum local = 7.68585702259e-10, global = -3.36686161497e-10, cumulative = 3.04788137997e-07
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.47947887396e-07, Final residual = 9.72671115058e-10, No Iterations 464
    time step continuity errors : sum local = 6.73366904995e-10, global = 7.62059709157e-11, cumulative = 3.04864343968e-07
    smoothSolver:  Solving for epsilon, Initial residual = 0.000146170038677, Final residual = 5.6866402182e-09, No Iterations 8
    smoothSolver:  Solving for k, Initial residual = 0.00238502525139, Final residual = 9.25193944196e-09, No Iterations 20
    ExecutionTime = 689.06 s  ClockTime = 695 s
    

    尝试修改为层流之后,连续迭代,次数基本成一个下降趋势,但是受力情况还是不变,错误依旧

    如下:

    Updating StokesII wave model for patch inlet
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.33810842593e-07, Final residual = 1.05803239609e-09, No Iterations 95
    time step continuity errors : sum local = 2.20723157353e-08, global = 6.73895228141e-09, cumulative = 7.30414848751e-06
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.58642796578e-09, Final residual = 9.05674029848e-10, No Iterations 23
    time step continuity errors : sum local = 1.87310827339e-08, global = 5.04129159013e-09, cumulative = 7.3091897791e-06
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.88155187954e-09, Final residual = 9.84985998753e-10, No Iterations 1
    time step continuity errors : sum local = 1.99973561049e-08, global = 4.66563789626e-09, cumulative = 7.313855417e-06
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.11812493747e-09, Final residual = 9.39334225107e-10, No Iterations 1
    time step continuity errors : sum local = 1.91664654171e-08, global = 4.75521938246e-09, cumulative = 7.31861063638e-06
    ExecutionTime = 8152.77 s  ClockTime = 8159 s
    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #9

    我这没装wave2foam,所以.. :135: 看看其他大佬有什么意见

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

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

    虽然我算例的浮体和你不一样,但是问题很类似:在波浪还没有传播到浮体附近之前,浮体就受力振荡,p_rgh的迭代次数不是递减的趋势,感觉我就是在踩你踩过的坑:135: @tower
    overset: oversetMesh.zip 这个算例用的是OF1906自带的造波消波模块吧,不需要waves2Foam造波消波,@东岳 老师能不能给我们一点建议o_0,看看应该从什么角度着手尝试解决这个问题?谢谢!
    也希望有其他前辈能够就这个问题指点一下,缩短小白踩坑的时间,感谢!!

    1 条回复 最后回复
  • C 离线
    C 离线
    CFDngu
    写于 最后由 编辑
    #11

    我用overset计算过浮体大幅度运动的问题,没遇到你们说的问题。
    第一种可能是设置问题,可以仔细检查一下算例的设置,包括overset网格的设置;
    第二中可能就是waves2Foam与overset网格结合的问题。这个我也不确定,因为我之前算的时候采用了自己写的造波程序。

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

    @CFDngu 谢谢建议!!

    1 条回复 最后回复
  • 卡 离线
    卡 离线
    卡洛
    写于 最后由 编辑
    #13

    老哥,请问最终问题怎么解决的,我最近也遇到这个问题了

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

    老哥,我来歪个楼,问你一下,wave2foam作用是啥?原版OpenFOAM不就有各种波浪么,或者自己自建个波浪边界条件也行吧,为啥还要额外一个wave2foam呢?

    1 条回复 最后回复
  • J 离线
    J 离线
    jlx
    在 中回复了 tower 最后由 编辑
    #15

    @tower 老哥,请问下问题解决了么:143:

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

    老哥们,你们是如何解决这个问题
    我在使用两套overset做浮体自衰减的时候也遇到类似的问题

    1 条回复 最后回复
  • W 离线
    W 离线
    wangwang
    写于 最后由 编辑
    #17

    你好,请问这个问题解决了吗?

    1 条回复 最后回复
  • W 离线
    W 离线
    wzlr2023
    写于 最后由 编辑
    #18

    大家好,请问这个问题解决了吗?我用overinerDyMFoam算一个二维方块浮体也遇到了这个问题,更改了很多参数都没解决。。。

    1 条回复 最后回复

  • 登录

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