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

[LYT]

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

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

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

[李东岳]

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

[LYT]

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

[李东岳]

他那个sigmaDotU是已经离散之后的了，你看他成了面矢量，类似$\mathrm{\nabla }\cdot \mathbf{U}$，里面$\mathbf{U}$是体矢量，但是离散后的形式是$\sum {\mathbf{U}}_f\cdot {\mathbf{S}}_f$，里面${\mathbf{U}}_f\cdot {\mathbf{S}}_f$就是标量，但是是面标量。

$\mathrm{\nabla }\cdot (\tau \cdot \mathbf{U})$中$\tau \cdot \mathbf{U}$是体矢量，离散后是$\sum (\tau \cdot \mathbf{U}{)}_f\cdot {\mathbf{S}}_f$，$(\tau \cdot \mathbf{U}{)}_f\cdot {\mathbf{S}}_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 期待测试结果！