■Seepage Face の点Aにq_n=0 を与える方法■
では、どうすれば seepage face の点 A に q_n=0 を与えられるかを説明します。まず、計算のステップ1として、下図の境界条件で h(x) と q_n を計算します。つまり、普通の境界要素法の計算ですね。

ステップ2 として、次の2つを行います。

計算結果より、Free surface の z 座標を z=h(x) とします。

次に、seepage face 側の点 A の qn=0 をチェックします。もし、qn＜0 (流れが外から領域へ) の場合は、点 A の z 座標を z=zA-α(zA-zB)/2 とします。定数α＝1.0です。すると次のページの図の様な新しい境界が出来上がります。参考までに、上の図の場合、必ず qn＜0 になります。理由は、aquifer の中に上図の水位を保つために、外から (seepage face を通過して) 水を補給しなくてはならないからです。