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. Cyclic周期性边界条件

Cyclic周期性边界条件

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

    another same issue https://www.cfd-online.com/Forums/openfoam-solving/199582-twophaseeulerfoam-dispersed-phase-not-conservative.html

    1 条回复 最后回复
  • R 离线
    R 离线
    random_ran 大神
    写于 最后由 编辑
    #22

    @cxzhmg

    我对 OF 理解也不是特别深。特别是你这里还有多相流的问题,求解器我也是从来没用过的。

    我只能从非常有限的角度去解读,或许能帮助到你,但不能保证你能解决问题。最终还是要靠你自己的。

    今天晚上,我来好好研究一下。

    Yours in CFD,

    Ran

    1 条回复 最后回复
  • R 离线
    R 离线
    random_ran 大神
    写于 最后由 编辑
    #23

    @cxzhmg

    你这个 case 我 run 的好长时间啊。

    半个小时才跑了 1/4。

    我明天继续研究。

    Yours in CFD,

    Ran

    C 1 条回复 最后回复
  • C 离线
    C 离线
    cxzhmg
    在 中回复了 random_ran 最后由 编辑
    #24

    @random_ran 谢谢您的关注。我在controlDict里面输出了全场的固含率、速度、压力,不需要计算完就可以看到质量明显的不守恒。

    如果您把fvOption加的源项去掉,仅在重力作用下计算还会出现颗粒相沿重力反方向运动的非物理现象……

    但是如果我不使用周期边界,上下边界也设成壁面,是会出现比较正常的固相堆积在底部,气相被挤压到上部运动趋势和结果,且质量是守恒的。

    所以我考虑问题应该还是出在cyclic边界上,具体需要看一下cyclic怎么实现,对于单向流和多相流的差别在什么地方,我最近也好好研究一下。如果您有什么想法和发现,欢迎讨论和指导。

    C 1 条回复 最后回复
  • C 离线
    C 离线
    Calf.Z-DNS
    在 中回复了 cxzhmg 最后由 编辑
    #25

    @cxzhmg 我之前跟同学讨论的是泄压可能是因为设置cyclic的效果就像是一道虚设的面,面上可以流进流出,而真正有阻挡作用的是wall壁面,因此channel就像是上下两个平板之前存在一个自由流动的空间一样,于是留着留着就泄压。但如果像圆管那样四周都是wall,就不会泄压。另外我采用了低马赫数假设,压力的变动不影响物性计算。因此无论压力是否泄压,我的物性计算还是正常的。

    C 1 条回复 最后回复
  • C 离线
    C 离线
    cxzhmg
    在 中回复了 Calf.Z-DNS 最后由 编辑
    #26

    @Calf-Z-DNS 圆管不会泄压吗?我计算二维的box,左右壁面,也是会出现泄压的问题。如果cyclic确实是周期边界的话,我觉得即使是3维无限大平板间流动,或者三维全周期流动,也都不应该出现泄压或者质量不守恒的问题,因为从任何边界流出去的都会从对应的neighbour流进来,所以我怀疑openFOAM里面的cyclic并不是准确的periodic,当然具体还是要看代码,我最近正在看。

    1 条回复 最后回复
  • X 离线
    X 离线
    XieXiaoyang
    在 中回复了 气固两相流 最后由 编辑
    #27

    @气固两相流 您好,我在simpleFoam里加了能量方程想测试一下fixedJump边界,是三维槽道算例,边界上会出现如下图的情况,在槽道出入口处有些网格上的结果感觉比较异常。
    e073929a-9715-4902-a7ab-0ef36e611985-5e59dbe54fe19e859d198c523d01e88.png

    该算例中采用fvoptions添加动量源项实现动量方程的周期边界,因此速度的边界就是cyclic,在出入口处云图就很平滑:
    73d9cd0a-4354-4047-936f-eac46dfc40da-006f1aa775c92625d6b2577d4f1180e.png

    除此以外,我在对一个横掠管束的算例中对压力采用fixedJump边界的时候,也出现了这种现象(压力场):
    c95b6f25-b40d-47bc-9509-e33b1ccaeb5c-3786a922e1bf9af8b00c8ff47a7484f.png
    但是在采用fvoptions添加动量源项时边界处结果就比较正常(虽然对管束这么实现周期边界是不正确的):
    1ecd55db-46bf-4756-a240-398e9e227944-7b2e764ea37dffc321d394f746dd1b2.png

    还有就是采用这种边界,温度场好像会因为数值误差的累计(存疑)使得整体计算区域温度越算越高,下图分别是采用simpleFoam刚开始计算和很久之后的温度场分布,明显发现计算很久之后的温度场温度整体偏高,按理说稳态计算应该是不变的才对。(另:这个网格原来是用来算DNS的,所以我觉得算很久之后出现的这种温度场波动应该是算出湍流脉动了):
    d3ae393b-6c7c-482e-ad35-a2ff1ada56f9-602d0132708b0bc47da524728bfd084.png
    477fd162-9ba3-40dc-968a-958ac0f2f145-19dcbff2a854ca3f48e5b2e83df3c83.png

    这是我最近调试fixedJump边界遇到的一些问题,烦请大佬有时间帮忙瞅瞅~
    如有描述不清的地方,还请多多指教!

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

    @xiexiaoyang simpleFoam怎么加入的温度场

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

    X 1 条回复 最后回复
  • X 离线
    X 离线
    XieXiaoyang
    在 中回复了 李东岳 最后由 编辑
    #29

    @李东岳 就是在simpleFoam求解器中添加了一个温度的标量输运方程27ef07d3-19fa-408d-b116-681cf2cca174-image.png

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

    @xiexiaoyang https://openfoam.org/release/2-2-0/boundary-conditions/ 我看了下这个边界条件的介绍,看起来更倾向于在计算域内部添加baffle然后实现某种jump。类似计算域内部添加的fan边界条件。所以你目前这种前后面肯定不是baffle,经过你测试,如果用起来有问题,那应该就是不能用这个fixedJump边界条件。

    simpleFoam加温度,为何不用rhoSimpleFoam。

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

    X 3 条回复 最后回复
  • X 离线
    X 离线
    XieXiaoyang
    在 中回复了 李东岳 最后由 编辑
    #31

    @李东岳 但是就计算结果来说,确实实现了温度边界的阶跃,而且楼上有老师对二维算例的压力边界计算结果,效果也都还好。
    就是不太明白为什么会出现这样部分网格计算结果不好:136:

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

    @李东岳 没用rhosimpleFoam是因为它解出来的好像是内能e,还涉及到热物性相关内容,不敢乱来hhh

    1 条回复 最后回复
  • X 离线
    X 离线
    XieXiaoyang
    在 中回复了 李东岳 最后由 编辑
    #33

    @李东岳 东岳老师,jump边界上的问题解决了,可能是因为在并行计算的时候,没有将关联的周期边界放在一起引起的,加了:
    b062b143-55b5-4f91-abcb-06ee8002b0c8-image.png
    这几句以后的计算结果种jump关联边界上的问题基本就没有了。
    ba8a96d1-9ada-4035-9ada-94abaa7307bc-image.png
    6b75a1ff-15cf-4790-98db-92929a33835d-image.png

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

    @xiexiaoyang 感谢反馈!

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

    李东岳李 X 2 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    在 中回复了 李东岳 最后由 编辑
    #35

    我测试了一下这个fixedJump,会一直一直的导致速度上升,你多算几个周期看看,

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

    X 1 条回复 最后回复
  • X 离线
    X 离线
    XieXiaoyang
    在 中回复了 李东岳 最后由 编辑
    #36

    @李东岳 是对压力边界采用fixedJump边界吗?
    我先用层流稳态算例尝试一下看看~

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

    @xiexiaoyang 是的

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

    1 条回复 最后回复
  • X 离线
    X 离线
    XieXiaoyang
    在 中回复了 李东岳 最后由 编辑
    #38

    @李东岳 李老师,对于周期边界下的压力场,piso循环中有pEqn.setReference(pRefCell,pRefValue);用来确定压力场。那在计算周期边界下的被动标量时,应该用什么来对温度场进行修正呢?
    :135:

    1 条回复 最后回复

  • 登录

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