compressibleInterFoam能量方程疑惑

[mohui]

今天看了compressibleInterFoam的Teqn方程，如下：

fvScalarMatrix TEqn
    (
        fvm::ddt(rho, T)
      + fvm::div(rhoPhi, T)
      - fvm::laplacian(mixture.alphaEff(turbulence->mut()), T)
      + (
            fvc::div(fvc::absolute(phi, U), p)
          + fvc::ddt(rho, K) + fvc::div(rhoPhi, K)
        )
       *(
           alpha1/mixture.thermo1().Cv()
         + alpha2/mixture.thermo2().Cv()
        )
    );

看起来是总能方程，但是还是有点对不上，特别是后面的相乘的部分完全对不上啊。
能量方程推导过程，希望大神们给点指导意见。

[李东岳]

这个求解的是多相的温度方程，推导温度方程的时候，需要在e方程或者h方程的基础上，继续推，你尝试推一下？然后我更新到网站 :cheeky: 或者9月份之后我自己推一下也可以。
另外这个是多相的，需要适配。

[chengan.wang]

@东岳 请教群主，- fvm::laplacian(mixture.alphaEff(turbulence->mut()), T)对应的是什么方程呢？

[李东岳]

从这一行代码可以看出来是$\nabla\cdot(\alpha_{...}\nabla T)$，其中$\alpha_{...}$对应的是alphaEff

[chengan.wang]

@东岳 导热系数除以热容=alpha,表示什么意思呢？还有turbulence->mut()，表示什么呢

[李东岳]

turbulence->mut()表示调用turbulence对象的mut()函数，->是指针，主要是C++的东西。

[chengan.wang]

@东岳 谢谢群主，我大概明白这段话是什么作用了，调用湍流动力粘度，结合粘度模型，可以获得导热系数与热容的比值

[Rachel0096]

@chengan-wang
您好，想问一下，你知道如何根据openfoam中的内能守恒方程，推导得到compressibleInterFoam中的能量方程嘛？知乎上面说是乘以alpha/Cv，但我自己推导有点问题
div(rhoPhi, T)和div(Phi，T)有什么不一样呢？ fvc::div(fvc::absolute(phi, U), p)这一项中absolute是什么意思呢

[李东岳]

我插一句：openfoam这个compressibleInterFoam的温度方程，应该是多相流里面最后一个尚未完善的地方了。温度并不是守恒变量。理想情况应该求解e方程。