icoFoam中压力泊松方程的求解
-
压力泊松方程的离散问题
of中构建的压力泊松方程似乎是这样的形式:
$\nabla\cdot(\frac{1}{a_P}\nabla p)=\nabla\cdot(HbyA)$
(icoFoam解析中的式27、33)按我的理解,对上面的方程似乎应该写成下面这种离散形式?
fvm::laplacian(rAU, p) == fvc::div(HbyA)
为什么实际上icoFoam使用了phiHbyA作为中间变量呢?
phiHybA的问题:
phiHybA是一个surfaceScalarField 对象,即面上标量场;随后求解器执行了对这个面上标量的散度:
fvc::div(phiHbyA)
为什么可以对一个面上标量计算散度呢?是否div函数在对surfaceScalarField作用的时候采用了不同的方式?我查了div的代码,确实看到div对surfaceScalarField似乎有不同的处理,但是我没看懂。请问这是否有关系呢?
期待各位大佬指点