Skip to content
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]
Skins
  • 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

  • Default (No Skin)
  • No Skin
Collapse
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界

reactingTwoPhaseEulerFoam求解器的组分输运方程中添加质量源相导致质量分数无界

Scheduled Pinned Locked Moved OpenFOAM
15 Posts 3 Posters 12.3k Views
  • Oldest to Newest
  • Newest to Oldest
  • Most Votes
Reply
  • Reply as topic
Log in to reply
This topic has been deleted. Only users with topic management privileges can see it.
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    wrote on last edited by
    #6

    Y方程没有详细看过。粗略看一眼并没有用MULES求解。另外,带源项要有界可以参考可压缩的相方程

    9月CFD算法编程课: http://dyfluid.com/class.html

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

    Z 1 Reply Last reply
  • Y Offline
    Y Offline
    yfclark 神
    wrote on last edited by
    #7

    @zwl 会不会是源项处理有问题,求解不稳定,Interface Composition搞了个fvm::Sp:
    c86d36d2-987e-433f-b3ad-2d3f527a8772-image.png
    此外可以考虑更换对流项格式,Laneer,Gamma01等

    Z 1 Reply Last reply
  • Z Offline
    Z Offline
    zwl
    replied to yfclark on last edited by
    #8

    @yfclark 非常感谢你的建议。源相我确实是直接添加的,没有再做处理,这个Sp会有什么作用(虽然有时我也会模仿着这样处理)?越界的主要原因是什么?

    1 Reply Last reply
  • Z Offline
    Z Offline
    zwl
    replied to 李东岳 on last edited by
    #9

    @东岳 那么下面带拉普拉斯项的组分输运方程如何采用MULES求解呢?谢谢。

            fvm::ddt(alpha, rho, Yi)
          + fvm::div(alphaRhoPhi, Yi, "div(" + alphaRhoPhi.name() + ",Yi)")
    
          - fvm::laplacian
            (
                fvc::interpolate(alpha)
               *fvc::interpolate(this->muEff()/Sct_),
                Yi
            )
         ==
            alpha*this->R(Yi)
           + Source
          + fvc::ddt(residualAlpha_*rho, Yi)
          - fvm::ddt(residualAlpha_*rho, Yi)
    
    
    1 Reply Last reply
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    wrote on last edited by
    #10

    扩散项+对流项用高阶/低阶格式处理,通过反扩散算法做出通量,源项进Sp或者Su就行

    9月CFD算法编程课: http://dyfluid.com/class.html

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

    Z 1 Reply Last reply
  • Z Offline
    Z Offline
    zwl
    replied to 李东岳 on last edited by
    #11

    @东岳 我看MULES代码,感觉他就是为求解相方程定制的,没有扩散相的处理,所以我不知道这里的扩散相放在哪个位置,怎么处理。

    李东岳李 2 Replies Last reply
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    replied to zwl on last edited by
    #12

    @zwl 扩散项、梯度项、对流项都是通过高斯定律化为$\sum$加和的形式,都是一样的,以梯度项举例:

    捕获.JPG

    9月CFD算法编程课: http://dyfluid.com/class.html

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

    1 Reply Last reply
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    wrote on last edited by
    #13

    对了,扩散项不会引起越界问题,或许你都不用处理了,扩散项是好朋友

    9月CFD算法编程课: http://dyfluid.com/class.html

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

    1 Reply Last reply
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    replied to zwl on last edited by
    #14

    @zwl 你在研究E-E模型下面的组分传输么

    9月CFD算法编程课: http://dyfluid.com/class.html

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

    1 Reply Last reply
  • Z Offline
    Z Offline
    zwl
    wrote on last edited by
    #15

    是的,求解器自带的是蒸发,我想计算空化,二者的计算原理大不相同,我更改了质量源相,导致质量分数越界。

    1 Reply Last reply

  • Login

  • Login or register to search.
  • First post
    Last post
0
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]