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

    边界条件:
    进口:速度进口,压力零梯度
    出口:速度零梯度,压力固定值0
    壁面:速度固定值(0 0 0),压力零梯度

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

    运算的log文件贴一下?simpleFoam的log

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

    1 条回复 最后回复
  • C 离线
    C 离线
    cfd-cxb
    写于 最后由 编辑
    #6

    您好,算例运行不起来,只能跑一步

    Create time
    
    Create mesh for time = 0
    
    Initializing the GGI interpolator between master/shadow patches: GUIDE-OUTLET/IMPELLER-INLET
    Initializing the GGI interpolator between master/shadow patches: IMPELLER-OUTLET/VOLUTE-INLET
    
    SIMPLE: no convergence criteria found. Calculations will run for 1400 steps.
    
    Reading field p
    
    Reading field U
    
    Reading/calculating face flux field phi
    
    Selecting incompressible transport model Newtonian
    Selecting RAS turbulence model RNGkEpsilon
    RNGkEpsilonCoeffs
    {
        Cmu             0.0845;
        C1              1.42;
        C2              1.68;
        alphak          1.39;
        alphaEps        1.39;
        eta0            4.38;
        beta            0.012;
        sigmak          0.71942;
        sigmaEps        0.71942;
    }
    
    Creating MRF for cell zone IMPELLER. rpm = 2900
    Excluded patches: 2(3 4)
    
    Starting time loop
    
    Time = 1
    
    smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0108375, No Iterations 2
    smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0104409, No Iterations 2
    smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0104084, No Iterations 2
    GAMG:  Solving for p, Initial residual = 1, Final residual = 0.027634, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.0357851, Final residual = 0.000575861, No Iterations 2
    time step continuity errors : sum local = 7.76946, global = 0.352411, cumulative = 0.352411
    smoothSolver:  Solving for epsilon, Initial residual = 0.13384, Final residual = 0.00135303, No Iterations 2
    bounding epsilon, min: -2.55034 max: 52.4243 average: 2.34921
    smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0250789, No Iterations 2
    ExecutionTime = 12.85 s  ClockTime = 13 s
    
    Time = 2
    
    smoothSolver:  Solving for Ux, Initial residual = 0.20708, Final residual = 0.00402441, No Iterations 2
    smoothSolver:  Solving for Uy, Initial residual = 0.249236, Final residual = 0.0056065, No Iterations 2
    smoothSolver:  Solving for Uz, Initial residual = 0.222885, Final residual = 0.00479877, No Iterations 2
    1 条回复 最后回复
  • C 离线
    C 离线
    cfd-cxb
    写于 最后由 编辑
    #7

    我换用linearUpwind后到200多步还是发散了,谢谢您!

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

    调节一下你的k和epsilon初始的值,从log来看你的epsilon不可实现。

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

    1 条回复 最后回复
  • C 离线
    C 离线
    cfd-cxb
    写于 最后由 编辑
    #9

    首先,非常感谢您的回复
    我调整了k以及epsilon,但是还是不可以,当我将div(phi,U)调为upwind或者linearUpwind Gauss时,div(phi,k)与div(phi,epsilon)不管是那种格式都能算,但是若是我将div(phi,U)设为limitedlinearV 1格式,div(phi,k)与div(phi,epsilon)不管那种格式都不行,算例都只算一步。

    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵一铭
    写于 最后由 编辑
    #10

    我发现个问题,你那个同时调用GGI和MRF?

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

    恩恩,是的,交界面处理使用的是GGI

    1 条回复 最后回复
  • 赵 离线
    赵 离线
    赵一铭
    写于 最后由 编辑
    #12

    我记得ggi和MRF不能同时用的。。

    1 条回复 最后回复
  • C 离线
    C 离线
    cfd-cxb
    写于 最后由 编辑
    #13

    谢谢您的回复,我用的是foam-extend版本,里面有同时调用的例子的

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

    我调整了k以及epsilon,但是还是不可以,当我将div(phi,U)调为upwind或者linearUpwind Gauss时,div(phi,k)与div(phi,epsilon)不管是那种格式都能算,但是若是我将div(phi,U)设为limitedlinearV 1格式,div(phi,k)与div(phi,epsilon)不管那种格式都不行,算例都只算一步。

    如果是用稳态算例,可以先迎风算出来流场之后,再加阶数,SIMPLEFoam很稳健,一般不会出现问题。我觉得还是设置的问题。

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

    1 条回复 最后回复
  • C 离线
    C 离线
    cfd-cxb
    写于 最后由 编辑
    #15

    明白了,谢谢前辈

    1 条回复 最后回复

  • 登录

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