势流算法问题
-
哪个第一个图？
另外，按照 http://dyfluid.com/potentialFoam.html 上面写的，目前我们要证明的是，链接里面的方程10（进口固定梯度，出口固定值），与方程11（进口零法向梯度，出口固定值，附加源项）的等价。考虑4网格系统，方程10与11离散出来的矩阵都是：
\begin{equation}
\begin{bmatrix}
1 & -1 & & \\
-1 & 2 & -1 & \\
& -1 & 2 & -1 \\
& & -1 & 3 \\
\end{bmatrix}
\begin{bmatrix}
\Phi_1 \\
\Phi_2 \\
\Phi_3 \\
\Phi_4 \\
\end{bmatrix} \equiv
\begin{bmatrix}
u_{inlet} \\
0 \\
0 \\
v \\
\end{bmatrix}
\end{equation}老师，方程(11)和进口固定梯度，出口固定值得到的矩阵我又手算了一遍。对于4网格系统:
1. 进口固定梯度，出口固定值
边界条件:
cell 1:
cell 2:
同理cell 3:
cell 4:
矩阵为:
2. 方程(11)
边界条件:
cell 1:
cell 2:
同理cell 3:
cell 4:
矩阵为:
和您算的系数矩阵第四行不一致，和等号右边的第一行和第四行也不一致，但计算过程我检查了好几遍，也没发现问题，问题出在哪儿了？
