想知道能量方程粘性项作用该怎么写
-
有▽·(τ)为turbulence.divDevRhoReff(U),那▽·(τ·u)该怎么写。。
rhocentral里的τ·u是这样写的。
我仿照的时候好像因为我是用的VOF的湍流模型,就没有muEff
有老师同学会写吗,或者说 我是不是出啥问题了。 -
τ·u这个是矢量,所以你的sigmaDotU应该是矢量,
-
@李东岳 但是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
-
@李东岳
李老师,报错说我sigmaDotU数据类型和方程对不上,该怎么解决
-
@lyt 解决了 我把-fvc::div(sigmaDotU)放在fvc::div(rhophi,k)后面去了,就对了 无语。
-
@lyt 难道是后面一定要加两个()(),才行?这是什么逻辑