Processing math: 100%

通量分解

通量在求解区域应该满足平衡方程,这种平衡关系反应了输运方程的物理本质,其中包含以下几点

  1. 相邻的两个节点通量会互相交换
  2. 一个节点上增加的量必须从另一个节点上减去
  3. 通量可以通过边界互相交换
  4. 通量能够通过反应减少
  5. 通量是由源产生的

前三点通过输运方程第一项(泄露项)反应,第四点通过第二项(移出项)反应,第五点由方程右端项反应。同时第三点也通过边界条件进行反应。

数值通量

通过定义数值通量Φ来反应相邻节点间的通量交换关系。其中Φ是一个稀疏矩阵,其稀疏性与网格的稀疏模式一致。其元素Φij表示从节点j转移到节点i的通量,显然Φij=Φji,即从节点i转移到节点j的等于负的节点j转移到节点i的量。即矩阵Φ是一个反对称矩阵。

通量在节点间的转移量通过输运方程第一项 (泄露项) 表示,有限元离散之后的半离散格式为 jcijϕj, 由于jcij=0,所以转移到节点 i 的量为
jΦij=jcijϕj=jcijϕj+ϕijcij=jcij(ϕi+ϕj)

可以将系数cij分解为反对称的aij和对称的bij,其中
cij=(aij+bij)/2
aij=cijcji
bij=cij+cji

所以
j(aij+bij)(ϕi+ϕj)/2=jΦij

通量分解为
Φij=aij(ϕi+ϕj)/2,ji
Φii=jbij(ϕi+ϕj)/2

##