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. rhoPimpleFoam求解器中EEqn.H能量方程具体形式是什么?

rhoPimpleFoam求解器中EEqn.H能量方程具体形式是什么?

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

    各位大佬好,对于能量方程,其以内能的形式为
    $$
    \frac{\partial}{\partial}(\rho(e+\frac{1}{2}\pmb{u}\cdot\pmb{u})) + \nabla\cdot [\pmb{u}(\rho e + \frac{1}{2}\rho \pmb{u}\cdot\pmb{u} + p)] = \nabla\cdot(\pmb{\sigma}\cdot\pmb{u}) + Q
    $$
    其中
    $$
    \pmb{\sigma} = \mu(\nabla\pmb{u} + \nabla\pmb{u}^T) - \frac{2}{3}\mu\nabla\cdot\pmb{u}\mathrm{I}
    $$

    rhoPimpleFoam求解器的EEqn.H中的能量方程代码如下

        volScalarField& he = thermo.he();
    
        fvScalarMatrix EEqn
        (
            fvm::ddt(rho, he) + fvm::div(phi, he)
          + fvc::ddt(rho, K) + fvc::div(phi, K)
          + (
                he.name() == "e"
              ? fvc::div
                (
                    fvc::absolute(phi/fvc::interpolate(rho), U),
                    p,
                    "div(phiv,p)"
                )
              : -dpdt
            )
          - fvm::laplacian(turbulence->alphaEff(), he)
         ==
            fvOptions(rho, he)
        );
    

    与能量公式对比发现代码中似乎没有$\nabla\cdot(\pmb{\sigma}\cdot\pmb{u})$,其被 fvm::laplacian(turbulence->alphaEff(), he)这一项替代了,这部分是怎么推导的呢?或者说是用了什么特殊的处理手法吗?

    W 1 条回复 最后回复
  • W 离线
    W 离线
    wangfei9088 大神
    在 中回复了 AppleKiller 最后由 wangfei9088 编辑
    #2

    @AppleKiller $\nabla \cdot (\sigma \cdot u)$是压力和粘性应力做功项,亚音速流动中忽略了。fvm::laplacian(turbulence->alphaEff(), he),能量he是标量,模化的湍流能量输运。

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

    http://dyfluid.com/energy.html 刚要发,发现二楼回复了,楼上大神说的对,被忽略了

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

    W A 2 条回复 最后回复
  • W 离线
    W 离线
    wangfei9088 大神
    在 中回复了 李东岳 最后由 编辑
    #4

    哈哈,东岳老师才是大神:146:

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

    @wangfei9088 各有所长各有所长 互相学习! :duang:

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

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

    @李东岳 在 rhoPimpleFoam求解器中EEqn.H能量方程具体形式是什么? 中说:

    http://dyfluid.com/energy.html 刚要发,发现二楼回复了,楼上大神说的对,被忽略了

    两位都是大神,那么rhoPimpleFoam一般不常用于求解跨音速或者超音速的流动吗?

    fvm::laplacian(turbulence->alphaEff(), he)这项对应热的部分是吧?

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

    rhoPimpleFoam能处理跨音速或者超音速,但是没有rhoCentralFoam准。rhoPimpleFoam更倾向于HVAC:

    Transient solver for turbulent flow of compressible fluids for HVAC and
         similar applications, with optional mesh motion and mesh topology changes.
    

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

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

    好的,谢谢李老师

    1 条回复 最后回复

  • 登录

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