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

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

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

@李东岳响应老师号召

CFD中文网

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