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

[LYT]

有▽·（τ）为turbulence.divDevRhoReff(U)，那▽·（τ·u）该怎么写。。
rhocentral里的τ·u是这样写的。

我仿照的时候好像因为我是用的VOF的湍流模型，就没有muEff

有老师同学会写吗，或者说 我是不是出啥问题了。

[李东岳]

τ·u这个是矢量，所以你的sigmaDotU应该是矢量，

[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 难道是后面一定要加两个（）（），才行？这是什么逻辑