圧力基準点の確認・変更について
東京理科大学倉渕研究室卒研生
度々の質問を失礼致します。
現在、STAR-CDとOpenFOAMで同じモデルを作成し、
その結果を比較する、ということを行っているのですが、
正確に比較するために、2種類のCFDで作成したモデルの
圧力基準点を合わせたいと思っております。
OpenFOAMにおいて、圧力基準点の確認や変更は
どのように行えばよいのでしょうか？
お忙しい中申し訳ありませんが、
よろしくお願いいたします。
Masashi Imano
流出側でpの境界条件が
_x
{
type fixedValue;
value uniform 0;
}
等と圧力のディレクレ境界条件になっている場合、この面での圧力が0基準となります。
どの境界にも圧力のディレクレ境界条件が無い場合には、
tutorials/incompressible/pisoFoam/ras/cavity
のsystem/fvSolutionで指定されるように、基準とする格子の番号と値を指定します。
pRefCell 0;
pRefValue 0;
ただ、格子番号がわからない場合には、
tutorials/multiphase/interDyMFoam/ras/damBreakWithObstacle/system/fvSolution
にあるように、 pRefPoint で基準の位置を指定することもできると思います。
pRefPoint (0.51 0.51 0.51);
pRefValue 0;
では、圧力のディレクレ境界を指定しつつ、任意の格子の圧力を基準にするには
どうするかというと、以下のようなsystem/sampleDictを作成し、計算結果から
sample ユーティリティで、任意の場所のpの値を取得して、その値を引くことを行ないます。
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object sampleDict;
}
setFormat raw;
surfaceFormat raw;
interpolationScheme cellPointFace;
fields
(
p
);
sets
(
pressureReferencePoint
{
type cloud;
axis xyz;
points
(
( 0 0 0 )
);
}
);
surfaces ();
ただし、simpleFoamといったOpenFOAMのincompressibleカテゴリのソルバーは、
流体の密度で割った圧力を解いているので、比較する時には注意してください。
もし、建築の分野で良く用いられる基準点での静圧を引いて基準点での動圧で割った
風圧係数を風洞実験値を含めて比較するのであれば、magUやkもサンプリングしてください。
fields
(
p
magU
k
);
2011/10/25 東京理科大学倉渕研究室卒研生 <miwa...@gmail.com>:
> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
> このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
> 詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。
>
>
--
IMANO Masashi, Ph.D.
Assistant Professor
Department of Architecture, Graduate School of Engineering,
The University of Tokyo
7-3-1, Hongo, Bunkyo-ku, Tokyo, Japan, 113-8656
E-mail:im...@arch.t.u-tokyo.ac.jp
Phone:+81-3-5841-6164(direct), +81-3-5841-6179(Labo.)
Facsimile:+81-3-5841-8511
東京理科大学倉渕研究室卒研生
OpenFOAMとは関係のない質問になってしまうのですが、
> もし、建築の分野で良く用いられる基準点での静圧を引いて基準点での動圧で割った
> 風圧係数を風洞実験値を含めて比較するのであれば、magUやkもサンプリングしてください。
この部分について、文献などを調べてみたのですが、理解するに至りませんでした。
申し訳ないのですが、よろしければ詳しく教えていただけないでしょうか？
参考文献などを教えていただけるだけでもありがたいです。
お手数をおかけいたしますが、よろしくお願い致します。
On 10月25日, 午後8:01, Masashi Imano <masashi.im...@gmail.com> wrote:
> 今野です。
>
> 流出側でpの境界条件が
>
> _x
> {
> type fixedValue;
> value uniform 0;
> }
>
> 等と圧力のディレクレ境界条件になっている場合、この面での圧力が0基準となります。
>
> どの境界にも圧力のディレクレ境界条件が無い場合には、
>
> tutorials/incompressible/pisoFoam/ras/cavity
>
> のsystem/fvSolutionで指定されるように、基準とする格子の番号と値を指定します。
>
> pRefCell 0;
> pRefValue 0;
>
> ただ、格子番号がわからない場合には、
>
>
>
>
>
> > お世話になっております。
> > 度々の質問を失礼致します。
>
> > 現在、STAR-CDとOpenFOAMで同じモデルを作成し、
> > その結果を比較する、ということを行っているのですが、
> > 正確に比較するために、2種類のCFDで作成したモデルの
> > 圧力基準点を合わせたいと思っております。
> > OpenFOAMにおいて、圧力基準点の確認や変更は
> > どのように行えばよいのでしょうか？
>
> > お忙しい中申し訳ありませんが、
> > よろしくお願いいたします。
>
> > --
> > このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> > このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
> > このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
> > 詳細については、http://groups.google.com/group/openfoam?hl=jaからこのグループにアクセスしてください。
>
> --
> IMANO Masashi, Ph.D.
> Assistant Professor
> Department of Architecture, Graduate School of Engineering,
> The University of Tokyo
> 7-3-1, Hongo, Bunkyo-ku, Tokyo, Japan, 113-8656
> E-mail:im...@arch.t.u-tokyo.ac.jp
> Phone:+81-3-5841-6164(direct), +81-3-5841-6179(Labo.)
>
> - 引用テキストを表示 -
Masashi Imano
magUは速度の大きさですから、基準点での動圧の計算に用います。
もし基準点での同じ高さの流入風の風速を用いるのであれば、magUは不要です。
kをサンプリングするのは、風圧係数を算出する際にpから2/3kを引くからですが、
これについては以下を参考にしてください。
http://www.ofwikija.org/index.php/WpcRAS
http://venus.iis.u-tokyo.ac.jp/personal/kato/CFD2.htm
ただしpから2/3kを引くべきかは、風洞実験での圧力の測定方法によると思います。
ここらへんは研究室で流儀が分れる所なので、研究室の方々に相談してください。
2011/10/26 東京理科大学倉渕研究室卒研生 <miwa...@gmail.com>:
東京理科大学倉渕研究室卒研生
返信が遅くなってしまい、申し訳ありません。
無事sampleを行うことができました。
参考資料まで教えていただき、本当にありがとうございました。
これからもよろしくお願いいたします。
> もし基準点での同じ高さの流入風の風速を用いるのであれば、magUは不要です。
>
> kをサンプリングするのは、風圧係数を算出する際にpから2/3kを引くからですが、
> これについては以下を参考にしてください。
>
>
> ただしpから2/3kを引くべきかは、風洞実験での圧力の測定方法によると思います。
> ここらへんは研究室で流儀が分れる所なので、研究室の方々に相談してください。
>
> 2011/10/26 東京理科大学倉渕研究室卒研生 <miwatf...@gmail.com>:
>
>
>
>
> > お世話になっております。
>
> > OpenFOAMとは関係のない質問になってしまうのですが、
>
> >> もし、建築の分野で良く用いられる基準点での静圧を引いて基準点での動圧で割った
> >> 風圧係数を風洞実験値を含めて比較するのであれば、magUやkもサンプリングしてください。
>
> > この部分について、文献などを調べてみたのですが、理解するに至りませんでした。
>
> > 申し訳ないのですが、よろしければ詳しく教えていただけないでしょうか？
> > 参考文献などを教えていただけるだけでもありがたいです。
>
> > お手数をおかけいたしますが、よろしくお願い致します。
>
> > On 10月25日, 午後8:01, Masashi Imano <masashi.im...@gmail.com> wrote:
> >> 今野です。
>
> >> 流出側でpの境界条件が
>
> >> _x
> >> {
> >> type fixedValue;
> >> value uniform 0;
> >> }
>
> >> 等と圧力のディレクレ境界条件になっている場合、この面での圧力が0基準となります。
>
> >> どの境界にも圧力のディレクレ境界条件が無い場合には、
>
> >> tutorials/incompressible/pisoFoam/ras/cavity
>
> >> のsystem/fvSolutionで指定されるように、基準とする格子の番号と値を指定します。
>
> >> pRefCell 0;
> >> pRefValue 0;
>
> >> ただ、格子番号がわからない場合には、
>