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. 同一个case更换网格后计算发散

同一个case更换网格后计算发散

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

    请教下各位老师,换了个其他模型的网格,网格数更多,条件全都一样,计算发散是为什么?有什么需要修改的吗?
    原网格

    Create time
    
    Create mesh for time = 0
    
    Time = 0
    
    Mesh stats
        points:           16482
        internal points:  0
        faces:            32240
        internal faces:   15760
        cells:            8000
        faces per cell:   6
        boundary patches: 7
        point zones:      0
        face zones:       0
        cell zones:       0
    
    Overall number of cells of each type:
        hexahedra:     8000
        prisms:        0
        wedges:        0
        pyramids:      0
        tet wedges:    0
        tetrahedra:    0
        polyhedra:     0
    
    Checking topology...
        Boundary definition OK.
        Cell to face addressing OK.
        Point usage OK.
        Upper triangular ordering OK.
        Face vertices OK.
        Number of regions: 1 (OK).
    
    Checking patch topology for multiply connected surfaces...
        Patch               Faces    Points   Surface topology                  
        inletTop            4        10       ok (non-closed singly connected)  
        inletBottom         4        10       ok (non-closed singly connected)  
        outletTop           4        10       ok (non-closed singly connected)  
        walls               464      936      ok (non-closed singly connected)  
        front               8000     8241     ok (non-closed singly connected)  
        back                8000     8241     ok (non-closed singly connected)  
        outletBottom        4        10       ok (non-closed singly connected)  
    
    Checking faceZone topology for multiply connected surfaces...
        No faceZones found.
    
    Checking basic cellZone addressing...
        No cellZones found.
    
    Checking geometry...
        Overall domain bounding box (-3.88578e-16 2.27374e-15 0) (10 2 0.2)
        Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
        Mesh has 2 solution (non-empty) directions (1 1 0)
        All edges aligned with or perpendicular to non-empty directions.
        Boundary openness (-2.71051e-19 2.43945e-18 6.33e-17) OK.
        Max cell openness = 1.73472e-16 OK.
        Max aspect ratio = 1 OK.
        Minimum face area = 0.0025. Maximum face area = 0.01.  Face area magnitudes OK.
        Min volume = 0.0005. Max volume = 0.0005.  Total volume = 4.  Cell volumes OK.
        Mesh non-orthogonality Max: 0 average: 0
        Non-orthogonality check OK.
        Face pyramids OK.
        Max skewness = 1.77636e-13 OK.
        Coupled point location match (average 0) OK.
    
    Mesh OK.
    
    End
    

    新网格

    Create time
    
    Create mesh for time = 0
    
    Time = 0
    
    Mesh stats
        points:           486356
        internal points:  0
        faces:            968451
        internal faces:   482097
        cells:            241758
        faces per cell:   6
        boundary patches: 5
        point zones:      0
        face zones:       0
        cell zones:       0
    
    Overall number of cells of each type:
        hexahedra:     241758
        prisms:        0
        wedges:        0
        pyramids:      0
        tet wedges:    0
        tetrahedra:    0
        polyhedra:     0
    
    Checking topology...
        Boundary definition OK.
        Cell to face addressing OK.
        Point usage OK.
        Upper triangular ordering OK.
        Face vertices OK.
        Number of regions: 1 (OK).
    
    Checking patch topology for multiply connected surfaces...
        Patch               Faces    Points   Surface topology                  
        INLET               198      398      ok (non-closed singly connected)  
        OUTLET              198      398      ok (non-closed singly connected)  
        TOP                 1221     2444     ok (non-closed singly connected)  
        WALL                1221     2444     ok (non-closed singly connected)  
        frontAndBack        483516   486356   ok (non-closed singly connected)  
    
    Checking faceZone topology for multiply connected surfaces...
        No faceZones found.
    
    Checking basic cellZone addressing...
        No cellZones found.
    
    Checking geometry...
        Overall domain bounding box (-3500 0 0) (10500 4000 15.7505)
        Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
        Mesh has 2 solution (non-empty) directions (1 1 0)
        All edges aligned with or perpendicular to non-empty directions.
        Boundary openness (9.7422e-20 2.77781e-17 4.16351e-14) OK.
        Max cell openness = 3.24713e-16 OK.
        Max aspect ratio = 629.147 OK.
        Minimum face area = 1. Maximum face area = 61981.3.  Face area magnitudes OK.
        Min volume = 15.7505. Max volume = 976235.  Total volume = 8.80452e+08.  Cell volumes OK.
        Mesh non-orthogonality Max: 32.1987 average: 12.372
        Non-orthogonality check OK.
        Face pyramids OK.
        Max skewness = 0.231851 OK.
        Coupled point location match (average 0) OK.
    
    Mesh OK.
    
    End
    
    
    1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    写于 最后由 编辑
    #2

    那就是网格导致的呗 很正常 网格会导致发散

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

    A 1 条回复 最后回复
  • A 离线
    A 离线
    Afr1yne
    在 中回复了 李东岳 最后由 编辑
    #3

    @李东岳 李老师,我看这个新网格质量也OK。是需要改离散格式边界条件这些,还是说对于这求解器我的新网格不过关?

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

    顶一下 更换更大计算域 更多网格 需要该边界条件的数值和离散格式吗?

    1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    写于 最后由 编辑
    #5

    基本不用修改,发散大概率是网格问题。看看fvSchemes和fvSolution怎么设置的,有的格式在好点的网格能算,复杂的网格就发散

    A 1 条回复 最后回复
  • A 离线
    A 离线
    Afr1yne
    在 中回复了 coolhhh 最后由 编辑
    #6

    @coolhhh 老师 我的网格就是二维矩形 大小均匀 只是数量更多 计算域更大了 会是质量问题吗

    C 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    在 中回复了 Afr1yne 最后由 编辑
    #7

    @Afr1yne 如果普通的均匀网格都发散,有点难搞,得再仔细查一查什么原因。fvSchemes和fvSolution改为更稳定的格式试试

    A 1 条回复 最后回复
  • A 离线
    A 离线
    Afr1yne
    在 中回复了 coolhhh 最后由 编辑
    #8

    @coolhhh 老师 边界条件的具体数值跟网格尺寸互相有什么影响吗 能不能是我换网格之后 u k e这种值给错了

    C 1 条回复 最后回复
  • C 离线
    C 离线
    coolhhh 神
    在 中回复了 Afr1yne 最后由 编辑
    #9

    @Afr1yne 如果模型什么都没变,只是网格加密, u k e正常也是一样的。把fvSchemes换成最稳定的格式计算看看,再逐个改二阶格式试算,看哪个设置导致的发散

    A 1 条回复 最后回复
  • A 离线
    A 离线
    Afr1yne
    在 中回复了 coolhhh 最后由 编辑
    #10

    @coolhhh 我去试下 谢谢老师指点

    1 条回复 最后回复

  • 登录

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