首页 理论教育叶片泵设计数值通用对流扩散方程的有限体积离散化

叶片泵设计数值通用对流扩散方程的有限体积离散化

【摘要】:回代入原通量密度式得其中Fe=,令aE=∣DeA∣,有Je=p-aEE 类似可得Jw=aWW-p 对南北方向类似地有Jn=p-aNN Js=aSS-p 将以上四项回代入原式得aPP=aEE+aWW+aNN+aSS+b 对界面上函数及其导数采取特定的构造格式,即确定了系数A和B的具体形式,就可以获得最终的离散表达式。

直角坐标下通用形式的二维对流扩散方程表示为

利用有限体积法对上式进行离散

将源项线性化为:S=SC+SPϕP

重新整理上式得

令式(1.2-4)中Fe=(ρueyFw=(ρuwyFn=(ρvnxFs=(ρvsx978-7-111-52131-0-Chapter01-17.jpg978-7-111-52131-0-Chapter01-18.jpg978-7-111-52131-0-Chapter01-19.jpg978-7-111-52131-0-Chapter01-20.jpg;令P=F/D,表征对流与扩散作用的相对大小。

为获得通用的离散格式,从原通用式中提取出东西方向的通量密度

再除以D

对式(1.2-6)取统一格式

这样,FD与离散格式联系起来,只要选取恰当的格式,即确定了系数BA的具体内容,就可以得到统一的离散方程。

回代入原通量密度式得

其中Fe=(DePe),令aE=∣DeAPe)∣,有

Je=(aE+Feϕp-aEϕE (1.2-9)

类似可得

Jw=aWϕW-(aw-Fwϕp (1.2-10)

对南北方向类似地有

Jn=(aN+Fnϕp-aNϕN (1.2-11)

Js=aSϕS-(aS-Fsϕp (1.2-12)将以上四项回代入原式得

aPϕP=aEϕE+aWϕW+aNϕN+aSϕS+b (1.2-13)

对界面上函数及其导数采取特定的构造格式,即确定了系数AB的具体形式,就可以获得最终的离散表达式。例如,由aE=∣DeAPe)∣可以求得aE,同样可以获得aWaNaS

上述方法同样适用于泵流动基本方程。对于流动方程,需将压力从源项中分离出,得到

AiJ为东侧或西侧的面积,AIj为南侧或北侧的面积。下标nb表示边界面,(ij)和(IJ)为单元节点编号。