compressibleInterFoam能量方程疑惑



  • 今天看了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月份之后我自己推一下也可以。
    另外这个是多相的,需要适配。



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



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



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



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



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


Log in to reply