关于“CFD中的能量方程”中最后的方程（16）

[CFD中文网]

@hurricane007
区别不大:expressionless: 哈哈

[hurricane007]

@cfd-china 嗯啊其实基本没区别。。。不过我有点疑惑的是一个用的是volScalarField 另一个是surfaceScalarField，所以差别应该是一个是体心的值，另一个是差值到面上了？

[CFD中文网]

@hurricane007
fvm::laplacian(A, p)中的A可以使体场也可以是面场，应该在fvm::laplacian()函数内部有重载，你可以看这个函数的C文件看看:expressionless:

[hurricane007]

@cfd-china 感谢指教！我还有个问题，为啥有的要用体场，而有的要用面场呢？这个应该是牵涉到对高斯定理散度体积分变为面积分？而就这个例子来说，把以前的体心值改成了面上的插值，数值上有什么影响？

[CFD中文网]

为啥有的要用体场，而有的要用面场呢？这个应该是牵涉到对高斯定理散度体积分变为面积分？

体场一般存储值，有限体积法大体都需要把体场变成面场，通过高斯定理处理通量求解NS方程，所以说通量是有限体积法重要的概念，也是有限体积法的优点。

而就这个例子来说，把以前的体心值改成了面上的插值，数值上有什么影响？

我建议你看看fvm::laplacian()函数重载，有可能处理体场的函数调用了处理面场的函数，这样的话结果是一样的。
https://github.com/OpenFOAM/OpenFOAM-2.2.x/blob/master/src/finiteVolume/finiteVolume/fvm/fvmLaplacian.C
大略看了下，一样的。

[fireztw]

@CFD中文网 @东岳 请问∇ ⋅ （M∇φ）该如何离散，M为volScalarField，非常数。
fvc::laplacian(M,phi)和(fvc::grad(M)&fvc::grad(phi)) + M*fvc::laplacian(phi)是等价的吗？

[李东岳]

fvc::laplacian(M,phi)和(fvc::grad(M)&fvc::grad(phi)) + M*fvc::laplacian(phi)是等价的吗？

数学上是等价的。数值上后者第一项可能会出现震荡问题，尤其在网格比较糙的情况下。

[wangfei9088]

@李东岳 老师，我发现东岳流体“CFD中的能量方程”可能的三处笔误，供您参考。

1. “到处”是不是“导出”？

2. 我觉得应该是将方程（20）代入到（18），得到的方程（22）应该缺少了热源项。因为在下面一行才忽略热源得到的方程（23）。

3. 由方程（26）（27）（18）得到的方程（29）（30）（31）应该都有粘性力做功和重力做功项的。是不是（29）和（30）缺少粘性力做功和重力做功项？（31）缺少重力做功项？

[李东岳]

非常感谢！已更新

最后是都忽略了，这样就没问题了吧

[wangfei9088]

@李东岳 哦哦，没问题了。