通量在求解区域应该满足平衡方程,这种平衡关系反应了输运方程的物理本质,其中包含以下几点
- 相邻的两个节点通量会互相交换
- 一个节点上增加的量必须从另一个节点上减去
- 通量可以通过边界互相交换
- 通量能够通过反应减少
- 通量是由源产生的
前三点通过输运方程第一项(泄露项)反应,第四点通过第二项(移出项)反应,第五点由方程右端项反应。同时第三点也通过边界条件进行反应。
数值通量
通过定义数值通量Φ来反应相邻节点间的通量交换关系。其中Φ是一个稀疏矩阵,其稀疏性与网格的稀疏模式一致。其元素Φij表示从节点j转移到节点i的通量,显然Φij=−Φji,即从节点i转移到节点j的等于负的节点j转移到节点i的量。即矩阵Φ是一个反对称矩阵。
通量在节点间的转移量通过输运方程第一项 (泄露项) 表示,有限元离散之后的半离散格式为 ∑jcijϕj, 由于∑jcij=0,所以转移到节点 i 的量为
∑jΦij=∑jcijϕj=∑jcijϕj+ϕi∑jcij=∑jcij(ϕi+ϕj)
可以将系数cij分解为反对称的aij和对称的bij,其中
cij=(aij+bij)/2
aij=cij−cji
bij=cij+cji
所以
∑j(aij+bij)(ϕi+ϕj)/2=∑jΦij
通量分解为
Φij=aij(ϕi+ϕj)/2,∀j≠i
Φii=∑jbij(ϕi+ϕj)/2
##