■計算結果■
まず、input data を見て下さい。最後の2行に注目して下さい。数字の1は、point source の数を示します。そして、0. 0. 1000. は、それぞれ、point source の位置を示す x 座標値、y 座標値、それにpoint source のマグニチュードを示す Q₁ です。Input ファイルの名前は、BEM2.DAT です。

計算結果は、ファイル名 SOLUTION.BEM に入ります。早速、結果のFluxをチェックしてみましょう。q_n は、約16ですから、これに円周(2*π*10) を乗ずると全Flux は、1005 になります。Point souece のQ₁=1000 ですから、満足のいく値になっていますね。誤差は。0.5%です。

次に、領域内の点について、解析精度を調べてみましょう。まず、例題の場合の厳密解は、下の様になります。

ここに、R=1、T=1、Q=1000 です。

厳密解と解析値を下表にまとめてありますので、確認して下さい。表を見ると、荒っぽい要素分割でも、かなり良い結果を出しているのに気付くと思います。

全ての点で、誤差は、2% 以下になっています。

次は、有限要素法でも紹介した Free surface の計算です。