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

    最近想添加一个湍动能的源项,不知道是否可以直接在k方程添加源项呢?

    1 条回复 最后回复
  • W 离线
    W 离线
    wwzhao 超神
    写于 最后由 编辑
    #2

    可以用 fvm::susp 添加。

    M 1 条回复 最后回复
  • M 离线
    M 离线
    mohui
    在 中回复了 wwzhao 最后由 编辑
    #3

    @wwzhao 准备怼of湍流模型的调用机制,导致扫了一下K-equation 是可以添加源项的,但是具体的还不是很清楚,比如说如何引入源项(因为源项是通过某个数学模型计算得到的)。

    W 1 条回复 最后回复
  • W 离线
    W 离线
    wwzhao 超神
    在 中回复了 mohui 最后由 编辑
    #4

    @mohui 打开源码,找到kEqn,添加源项,完成。

    M 1 条回复 最后回复
  • M 离线
    M 离线
    mohui
    在 中回复了 wwzhao 最后由 编辑
    #5

    @wwzhao 谢谢,我看到这个方程了,现在我在想如何引进源项。

    W 1 条回复 最后回复
  • W 离线
    W 离线
    wwzhao 超神
    在 中回复了 mohui 最后由 编辑
    #6

    @mohui 你可以把源项 IOobject 注册到内存中,再在湍流模型里通过 lookupObject 查找得到源项。

    M 1 条回复 最后回复
  • M 离线
    M 离线
    mohui
    在 中回复了 wwzhao 最后由 编辑
    #7

    @wwzhao 好的,谢谢。现在稍微了解到3.0以后版本的湍流模型的编译较为复杂,需要仔细看看。

    Y 1 条回复 最后回复
  • zym604Z 离线
    zym604Z 离线
    zym604
    写于 最后由 编辑
    #8

    用pimpleFoam的话可以直接在fvoptions里添加momentum source项。

    W 1 条回复 最后回复
  • W 离线
    W 离线
    wwzhao 超神
    在 中回复了 zym604 最后由 编辑
    #9

    @zym604 动量方程和湍动能方程是两个不同的概念,fvOptions只能在动量方程中使用。因为只有动量方程的源项中加了fvOptions(U),而湍动能方程没有。

    zym604Z Y 2 条回复 最后回复
  • zym604Z 离线
    zym604Z 离线
    zym604
    在 中回复了 wwzhao 最后由 编辑
    #10

    @wwzhao 啊,sorry,没看清题目。

    1 条回复 最后回复
  • Y 离线
    Y 离线
    yuan_neu
    在 中回复了 mohui 最后由 编辑
    #11

    @mohui 也不复杂吧,比如对于kepsilon,用incompressible turbulence模型,
    1 在ras目录下直接在kepsilon做更改,
    2 回到incompressible目录下,编译
    3 回到求解器,在编译一下
    ok 还算挺方便的
    当然也可以添加新的文件 mykepsilon,需要在incompressible目录下修改某子文件夹下的 xxxmodels 添加两句代码,就可以了

    1 条回复 最后回复
  • Y 离线
    Y 离线
    yuan_neu
    在 中回复了 wwzhao 最后由 编辑
    #12

    @wwzhao 其实很多模型里面已经有ksource项了,当然默认下这项为零,可以直接将这个source给激活就行了。虽然我没用过

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

    of5.0里面已经有fvOption了,

        tmp<fvScalarMatrix> epsEqn
        (
            fvm::ddt(alpha, rho, epsilon_)
          + fvm::div(alphaRhoPhi, epsilon_)
          - fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_)
         ==
            C1_*alpha()*rho()*G*epsilon_()/k_()
          - fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_)
          - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_)
          + epsilonSource()
          + fvOptions(alpha, rho, epsilon_)
        );
    

    或者复制模仿realizableKE的格式,赋值于epsilonSource() 或者 kSource().

    1 条回复 最后回复

  • 登录

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