• 最新
  • 版块
  • 热门
  • 东岳流体
皮肤
  • 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中文网

初学CFD最容易犯的错 公众号 2023年11月23日
CFD领域俄美友谊长存 公众号 2023年11月23日
这是CFD加特林 网站 2023年11月9日

湍流模型方程和代码问题

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

    编程小白请教各位老师两个不成熟的问题。。。

    图一中湍流模型方程左右都乘了一个相分数,那如果在两相流模型中,都是空气时,相分数为0,并且相分数的时间导数和空间导数也是0,这方程怎么解?
    1.png

    图二中左边的代码和右边的公式感觉不匹配,求老师们指点一下,十分感谢!
    123.png

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

    @是嘿嘿牛啊 在 湍流模型方程和代码问题 中说:

    图二中左边的代码和右边的公式感觉不匹配,求老师们指点一下,十分感谢!

    我大略看了下,好像公式写错了,我明天详细看一下

    图一中湍流模型方程左右都乘了一个相分数,那如果在两相流模型中,都是空气时,相分数为0,并且相分数的时间导数和空间导数也是0,这方程怎么解?

    这个是著名的数值问题。每个版本的openfoam处理不一样。这个问题非常有意思,我这几天更新在ns方程笔记里面。

    这就是phase intensive方法要处理的问题。有些文章叫做奇异。

    24年3月CFD课将在杭州举办 http://dyfluid.com/class.html
    论坛随机不定时开放注册
    如果需要我给调试算例 1)算例可公开 2)网格数少 3)确保可以跑通

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

    你第二个问题是哪个版本的openfoam,我这面openfoam-9跟你的代码不一样

    24年3月CFD课将在杭州举办 http://dyfluid.com/class.html
    论坛随机不定时开放注册
    如果需要我给调试算例 1)算例可公开 2)网格数少 3)确保可以跑通

    是 1 条回复 最后回复
  • 是 离线
    是 离线
    是嘿嘿牛啊
    replied to 李东岳 on 最后由 编辑
    #4

    @李东岳 李老师您好,我目前用的版本是v1812

    李东岳李 1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    replied to 是嘿嘿牛啊 on 最后由 编辑
    #5

    @是嘿嘿牛啊

    1. openfoam9跟openfoam1812代码不一样

    2. 我又细看了一下,虽然写的不一样,但是公式没问题,公式是哪里有问题呢?

    24年3月CFD课将在杭州举办 http://dyfluid.com/class.html
    论坛随机不定时开放注册
    如果需要我给调试算例 1)算例可公开 2)网格数少 3)确保可以跑通

    是 1 条回复 最后回复
  • 是 离线
    是 离线
    是嘿嘿牛啊
    replied to 李东岳 on 最后由 编辑
    #6

    @李东岳 公式没问题,只是左边红框内代码fvm::Sp(alpha()*rho()*epsilonByk(F1, tgradU()), k_)不知道是什么意思,是公式里面的哪一项?刚开始学of编程,求李老师指点一下,十分感谢!

    李东岳李 1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    replied to 是嘿嘿牛啊 on 最后由 编辑
    #7

    @是嘿嘿牛啊

    fvm::Sp(alpha()*rho()*epsilonByk(F1, tgradU()), k_)

    可以理解为alpha()*rho()*epsilonByk(F1, tgradU())*k
    epsilonByk(F1, tgradU())的定义,要参考源代码的定义:

    kOmegaSST<MomentumTransportModel, BasicMomentumTransportModel>::epsilonByk
    (
        const volScalarField::Internal& F1,
        const volScalarField::Internal& F2
    ) const
    {
        return betaStar_*omega_();
    }
    

    alpha()*rho()*epsilonByk(F1, tgradU())*k 就变成了alpha()*rho()*betaStar_*omega_()*k 。这个跟公式就可以对上了

    24年3月CFD课将在杭州举办 http://dyfluid.com/class.html
    论坛随机不定时开放注册
    如果需要我给调试算例 1)算例可公开 2)网格数少 3)确保可以跑通

    是 1 条回复 最后回复
  • 是 离线
    是 离线
    是嘿嘿牛啊
    replied to 李东岳 on 最后由 李东岳 编辑
    #8

    @李东岳 李老师,下面那一行代码alpha()*rho()*betaStar_*omegaInf_*kInf_和alpha()*rho()*betaStar_*omega_()*k有什么区别吗?
    我看上面代码中有对omegaInf_和kInf_定义,

        kInf_
        (
            dimensioned<scalar>::lookupOrAddToDict
            (
                "kInf",
                this->coeffDict_,
                k_.dimensions(),
                0
            )
        ),
        omegaInf_
        (
            dimensioned<scalar>::lookupOrAddToDict
            (
                "omegaInf",
                this->coeffDict_,
                omega_.dimensions(),
                0
            )
        )
    

    那是否可以理解成omegaInf_和kInf_就是0?不明白alpha()*rho()*betaStar_*omegaInf_*kInf_这串代码的意义,请教一下老师,谢谢!

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

    是0 ,但是提供了一个用户自定义的接口

    24年3月CFD课将在杭州举办 http://dyfluid.com/class.html
    论坛随机不定时开放注册
    如果需要我给调试算例 1)算例可公开 2)网格数少 3)确保可以跑通

    是 1 条回复 最后回复
  • 是 离线
    是 离线
    是嘿嘿牛啊
    replied to 李东岳 on 最后由 编辑
    #10

    @李东岳 明白了,谢谢李老师!

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

    图一中湍流模型方程左右都乘了一个相分数,那如果在两相流模型中,都是空气时,相分数为0,并且相分数的时间导数和空间导数也是0,这方程怎么解?

    http://dyfluid.com/reactingTwoPhaseEulerFoam.html 更新到这里了

    24年3月CFD课将在杭州举办 http://dyfluid.com/class.html
    论坛随机不定时开放注册
    如果需要我给调试算例 1)算例可公开 2)网格数少 3)确保可以跑通

    是 1 条回复 最后回复
  • 是 离线
    是 离线
    是嘿嘿牛啊
    replied to 李东岳 on 最后由 编辑
    #12

    @李东岳 正在学习,谢谢李老师!

    1 条回复 最后回复

CFD中文网2016-2023 | 京ICP备15017992号-2

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 热门
  • 东岳流体
  • 登录

  • 登录或注册以进行搜索。