想知道能量方程粘性项作用该怎么写

[LYT]

@李东岳 但是rhocentralFoam里就是这么写的，这一点我也很不明白为什么。后来我理解成面标量，他后面散度可能是面积分的意思，这样去说服自己。

[李东岳]

他那个sigmaDotU是已经离散之后的了，你看他成了面矢量，类似$\nabla\cdot\bfU$，里面$\bfU$是体矢量，但是离散后的形式是$\sum \bfU_f\cdot\bfS_f$，里面$\bfU_f\cdot\bfS_f$就是标量，但是是面标量。

$\nabla\cdot(\tau\cdot\bfU)$中$\tau\cdot\bfU$是体矢量，离散后是$\sum (\tau\cdot\bfU)_f\cdot\bfS_f$，$(\tau\cdot\bfU)_f\cdot\bfS_f$为面标量，这个就是sigmaDotU

[LYT]

@李东岳
李老师，报错说我sigmaDotU数据类型和方程对不上，该怎么解决

[LYT]

@lyt 解决了 我把-fvc：：div（sigmaDotU）放在fvc：：div（rhophi，k）后面去了，就对了 无语。

[LYT]

@lyt 难道是后面一定要加两个（）（），才行？这是什么逻辑

[biaowai]

应该写成fvc::div(turbulence->devRhoReff() & U)吧，最近准备研究对流换热问题，不太清楚OF rhoSimpleFOAM、rhoPimpleFoam中不考虑粘性热项影响大吗，打算加上这一项试一试

[李东岳]

@biaowai 期待测试结果！

[biaowai]

@李东岳 在单孔平板气膜冷却里，rhoSimpleFoam+SST+展向周期性边界，这一项还是比较重要的，没有这一项温度会扩散很慢，尤其是周期面上，fvOptions中可以直接设置：

dissiapation1
{
    type            viscousDissipation;
    selectionMode   cellSet;
    cellSet         sourceCells1;
}
参考
$FOAM_TUTORIALS/compressible/rhoPimpleFoam/RAS/TJunction
$FOAM_SRC/fvOptions/sources/derived/viscousDissipation 

[李东岳]

感谢反馈！非常好！我关注一下

[biaowai]

@李东岳 响应老师号召