Unbounded Region
2 Dim. Helmholtz Equation 16

■例題1: 波長(λ)=1, 円柱の直径=1■
では早速計算してみましょう。 手始めに、厳密解が存在するケースを計算してみます。下図を見て下さい。 波長(λ)が1の波が直径1の円柱にアプローチしている状況を表わしています。 波長が1ということは、wave number(k) が2πになることを意味します。

入力データを作成する前に、円柱の境界を要素で分割しておく必要があります。 ここでは、上図の点Aから時計方向に、境界を32この線形要素で分割しました。 データの作成にはSETLINQ10.FOR で作成しました。 すると BEM1.DAT というファイルが作成されます。 要素分割と節点番号の様子を下図に示します。 またこのプログラムは境界の要素分割の状況がエクセルで表示できるファイルBOUNDARYELEMENTS.OUTの出力します。

入力データに誤りが無いとプログラム HLM8LINQ.FOR は 次の3つのファイルを出力します。

ファイル名 内　　　容
WHCOEFFI.CNT 円柱の境界上の点Aからの円周上の距離と wave coefficientが格納されています グラフは表計算ソフトで作れます
SOLUTION.BEM Helmholtz equationの解析に関する全ての入力/出力データが入っています
BOUNDARY.BIN 円柱の境界上の{u}と{du/dn}が収められています。バイナリーファイルです

BACK NEXT
Menu Laplace equation Helmholtz equation