浮体に作用する周面摩擦の出力方法について

unread,

Jun 19, 2023, 1:26:31 AM6/19/23

to OpenFOAM

若原と申します。

OpenFOAM(interdymfoam)を用いて、静止流体中の垂直円柱を強制振動させ、流体力を評価しています。

現在、円柱をYaw方向に回転させて、周面摩擦による減衰を評価したいと考えています。

いくつかのサイトからwall shear stressを行うことで壁面せん断応力を出力できると拝見したのですが、上手く作用しません。

【やりたいこと】

①OpenFOAMでwall shear stressを直接出力する。

②Paraviewを用いて、計算を行い間接的に出力する。

非常に初歩的な質問だと思いますが、

どなたかにご教授いただければ幸いです。よろしくお願い致します。

unread,

Jun 19, 2023, 4:46:16 AM6/19/23

to OpenFOAM

若原様

kominamiです。

★１の方法について

下記の２つのサイトを見てください。

OpenFOAMのpostProcessについて

wallShearStressの計算

・・・ただし、controlDictファイルのfunctions { }内に下記のように記載する必要があります。

functions
{
wallShearStress1
{
type ...;
...
}
}

★２の方法について

ParaviewでフィルターGradientを使うとU(速度場)の勾配を計算できますが、壁面の値にするのが結構面倒です。

たしか、セル値でなくノード値を使って勾配を計算しないといかなかったと思います。昔のことだから詳しいことは忘れました。

１の方法で作成されるwallShearStressというField量のファイルは、Uやpのファイルと同様にparaviewでポスト処理ができます。

こちらの方法のほうが楽でしょう。

2023年6月19日月曜日 14:26:31 UTC+9 若原巧実:

unread,

Jun 20, 2023, 2:40:33 AM6/20/23

to OpenFOAM

kominami様

ご返信誠にありがとうございます。

また、ご丁寧に回答いただき心から感謝申し上げます。

ご教授いただいた知見を参考に、wallShearStressによる計算と可視化を行ってみます。

もし、また不明な点や行き詰ったことがあった場合は、ご相談に乗っていただけたら幸いです。