\sum {\left( {\left( {{\bf{HbyA}}_f^{n + 1} - \frac{1}{{{A_{P,f}}}}{\nabla _f}{p^{n + 1}}} \right) \cdot {{\bf{S}}_f}} \right)} = 0
\sum {\left( {{\bf{HbyA}}_f^{n + 1} \cdot {{\bf{S}}_f}} \right)} = \sum {\left( {\frac{1}{{{A_{P,f}}}}{\nabla _f}{p^{n + 1}} \cdot {{\bf{S}}_f}} \right)}
{\bf{Hby}}{{\bf{A}}^r} = {\bf{U}}_P^r - \frac{1}{{{A_P}}}\left( { - \sum {{A_N}{\bf{U}}_N^r + {{\bf{E}}_P}} } \right)
如上式，计算压力泊松方程时需要用到 HbyA 和 Ap 在面上的值。如果是边界单元（该单元包含边界面），则在求 \bf{HbyA}^r 时就要用到边界面上的 \bf{U}^r 。
如果是内部单元，用到内部面上的 \bf{HbyA}^r 值时直接用两边的单元中心 \bf{Hby}^r 插值即可。
如果是边界单元会包含边界面，则需要用到边界面上的值。 \bf{HbyA}^r 如何处理？直接用相邻单元中心的 \bf{HbyA}^r 代替？
特别的，如果给定边界面上的速度时如何处理？