ある面に働く力を求める方法について

[ryutaitaro]

ある面に働く力を求める方法について

初めて質問させて頂きます。ryutaitaroと申します。
単独翼に働く、揚力および抗力を求めたいのですが
うまくいきません。過去の質問を見たところ、simpleFunctionObject
を用いれば計測できそうなことが分かったのでダウンロードしてみた
のですが、うまくインストールできませんでした。
私のOpenFOAMのヴァージョンは2.0.1で、ubuntu 10.4において使用しています。
関係ないのかもしれませんがbisonもインストールしています。

みなさんはどうやって面に働く力などを求めていますか？
教えて頂けたら幸いです。

[Sakuma]

OpenFOAM初心者のsakumaです。

次のディスカッションが参考になりませんでしょうか。

http://groups.google.com/group/openfoam/browse_thread/thread/ff04285304fb24e/8a993080584e413c?lnk=gst&q=%E9%9D%A2%E3%81%AB%E5%8A%A0%E3%82%8F%E3%82%8B%E8%8D%B7%E9%87%8D%E3%81%AE%E7%AE%97%E5%87%BA%E6%96%B9%E6%B3%95#8a993080584e413c

OpenFOAM 2.0.xを利用していますが，このディスカッションで教えて頂きました方法で
面に加わる力やモーメントを計算しています。

OpenFOAM 2.0のUser GuideのU-101ページにあるforcesの利用が出来るのではないでしょうか。
controlDictの中にforces関数を入れて利用しています。

forces
{
type forces;
functionObjectLibs ( "libforces.so" );
outputControl timeStep;
outputInterval 10;
patches (face_1);
pName p;
UName U;
rhoName rhoInf;
log true;
valueOutput true;
rhoInf 1.225;
CofR (0.25 0.0 0.0);
}

計算を流すと力，モーメントがこのように書き出されます。

forces output:
forces(pressure, viscous)((-0.0653452 -1.25823 -0.428035)
(-0.000366532 0.00303802 -0.00694473))
moment(pressure, viscous)((0.0147672 -0.133257 0.367372) (0.000269752
-0.00209551 -0.00127223))

もし間違っていたらすみません。

[ryutaitaro]

OpenFOAM初心者のryutaitaroです。

素早い返信ありがとうございます。
同じ内容のディスカッションがもう行われていたのですね。
見つけ出せず恥ずかしいかぎりです。

これを使って揚力、抗力を求めたいと思います。
どうもありがとうございました。

ryutaitaro
On 11月13日, 午後10:03, Sakuma <sakura.maru...@gmail.com> wrote:
> OpenFOAM初心者のsakumaです。
>
> 次のディスカッションが参考になりませんでしょうか。
>
> http://groups.google.com/group/openfoam/browse_thread/thread/ff042853...