メッシュ作成について

[Y.U]

OpenFOAM初心者です。

現在DEXCS2012の中のOpenFOAMを使って、管内流の計算を使用としています。
simpleFOAMを使う予定のため、tutorialのpizeDairyを元にやっています。
そこでBlockMeshを行うとメッシングは終わるのですが、ParaFoamで確認しようとすると
ParaFoamが落ちてしまいます。
またsnappyHexMeshを実行すると、stlファイルを読めないと出てきてしまいます。
ちなみに解析に用いるstlファイルはblenderのswiftsnapを用いて作成しました。

なにが原因で動かないのかがわからず、困っています。
ご教授頂きたく思います。

[tomihara]

Y.U様

勉強会＠関西幹事の冨原と申します。

>ParaFoamが落ちてしまいます。

というのは、

paraFoamコマンドを実行したのち、paraviewが起動せずに落ちてしまうのか、

paraviewが起動した後、”Applyボタン”を押すとparaviewが落ちてしまうのか、

どのような状況でしょうか？

（ちなみにparaFoamコマンドはOpenFOAM用にparaviewを起動するコマンドです。）

ちなみに後者ならば、私もよくやってしまうのは

blockMeshで作成した境界条件（constant/polyMesh/boundary）と

0ディレクトリ内の初期条件（Uやp）の中の境界条件の設定がずれているケースです。

メッシュを確認するだけならば、

”Volume Fields”内のチェックを全てオフにしてから”Applyボタン”を押すと解決できます。

>snappyHexMeshを実行すると、stlファイルを読めないと出てきてしまいます。

こちらの件は、snappyHexMeshに詳しくないので詳しことはお伝えできないのですが、

「読めない」にもいろいろあるかと思います。

なので、エラーメッセージのログを張り付けると、ほかの方の返答に役立つかと思います。

以上、ご参考までに。

冨原

2013年1月21日月曜日 20時50分06秒 UTC+9 Y.U:

[Y.U]

冨原様

 

ご教授頂きましてありがとうございます。

 

早速その部分をチェックしてみます。

 

またsnappyHexMeshについてはおっしゃるとおり、

エラーの表示を貼り付けるようにします。

 

今後ともよろしくお願いいたします。

 

Y.U

2013年1月21日月曜日 22時32分32秒 UTC+9 tomihara:

[Y.U]

お世話になります。
Y.Uです。

先日投稿した問題のなかでsnappyHexMeshを実行した
際のエラーメッセージを以下に記します。

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.x-fa2c23b05ebe
Exec   : snappyHexMesh
Date   : Jan 22 2013
Time   : 22:54:29
Host   : "Ueda01"
PID    : 3562
Case   : /home/yoritomo/Desktop/pipe
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Read mesh in = 3.46 s

Overall mesh bounding box  : (-0.0164677 -0.0207032 -0.0249155) (0.103532 0.0992968 0.0350845)
Relative tolerance         : 1e-06
Absolute matching distance : 1.8e-07



--> FOAM FATAL ERROR:
Cannot find triSurfaceMesh starting from "/home/yoritomo/Desktop/pipe/constant/triSurface/outlet01.stl"

    From function triSurfaceMesh::checkFile(const fileName&, const fileName&)
    in file searchableSurface/triSurfaceMesh.C at line 116.

FOAM exiting

triSurfaceにはoutlet01.stlは置いているのですが、このようなメッセージと
なりました。

stlファイルがasciiでない等が問題なのでしょうか。
ちなみにstlファイルは先の投稿の通り、Blenderのswiftsnapを用いて作成して
おります。

以上よろしくお願いいたします。

Y.U



2013年1月21日月曜日 20時50分06秒 UTC+9 Y.U:

[E.Mogura]

E.Mogura です。

swiftsnap　を用いて作成したstlファイルは、

通常、snappyHexMeshDictを出力した場所（systemフォルダの下）

に作成されます。triSurfaceフォルダ以下を、

constantフォルダの下へ移動するなり、

シンボリックリンクの作成が必要になりますが、

出来ていますか？

2013年1月21日月曜日 20時50分06秒 UTC+9 Y.U:

[Y.U]

E.Mogura様

いつもお世話になります。
Y.Uです。

ご指摘頂いたところを確認したところ、
triSurfaceのフォルダがsystemフォルダ
以下に入っていました。
これをconstant以下に置いたところ、
snappyHexMeshは起動してメッシュが切れました。
ありがとうございました。

ただその後ポテンシャル流の計算をさせるべく
境界条件設定後にpotentialFoamを実行したところ
以下のエラーが発生してしまいました。

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 3

Reading field p



--> FOAM FATAL IO ERROR:
cannot find file

file: /home/yoritomo/Desktop/pipe/3/p at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting

yoritomo@Ueda01 ~/Desktop/pipe $ potentialFoam

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.x-fa2c23b05ebe

Exec   : potentialFoam
Date   : Jan 23 2013
Time   : 22:50:02
Host   : "Ueda01"
PID    : 4962

Case   : /home/yoritomo/Desktop/pipe
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 3

Reading field p



--> FOAM FATAL IO ERROR:
cannot find file

file: /home/yoritomo/Desktop/pipe/3/p at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting


またどこかでファイルの置き場所に間違いがあるため
のエラーでしょうか。

エラーの中に3フォルダ以下のpファイルを読み出している
部分がありますが、p,u,k等の境界条件のファイルは0フォルダ
に入れています。

何度もお聞きして申し訳ありませんが、
よろしくお願い致します。

Y.U


2013年1月23日水曜日 5時25分08秒 UTC+9 E.Mogura:

[tomihara]

Y.U様

冨原です。

横から失礼いたします。

>Create time
>Create mesh for time = 3

とありますので、3 sから計算を始めています。

controlDict内のstart timeが3になっているのではないでしょうか。

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

2013年1月23日水曜日 23時00分47秒 UTC+9 Y.U:

[ONO Hiroki]

小野です。
私も横から失礼いたします。

snappyHexMeshを実行すると、通常の場合、新たに1、2、3のディレクトリが作成されます。
これらにはそれぞれ、snappyHexMeshの各段階において生成されたメッシュのデータが入っています。
1→STLに合わせて指定回数refineし、重なる部分のメッシュを抜いただけの凸凹メッシュ
2→1からさらに、STLの表面形状に沿うようにメッシュ形状を変化させたもの
3→2をもとに、表面に沿うレイヤーメッシュを挿入したもの

(ここでは、3のメッシュを以降の計算に使用したいとします)

なので、ここから通常計算を行う場合は、たとえば
①3のディレクトリに初期条件(p,U)などを移す
もしくは、
②constant/polyMeshを、あらたにsnappyで生成された3/constant/polyMeshの内容に置き換える。(その後1,2,3は消去)
のような作業が必要です。
①、②はどちらかで構いません。やりやすい方で良いと思います。

ちなみに、snappyHexMeshを実行する際、-overwriteオプションをつけておくと、②が勝手に実行されます。





2013年1月23日水曜日 23時00分47秒 UTC+9 Y.U:

[Y.U]

小野様

いつもお世話になっております。
Y.Uです(これからは苗字の植田で投稿致します)。

ご指摘いただいた点を修正したところ
計算することが出来ました。
ファイルの置き場所には今後も気をつけます。

ありがとうございました。

植田



2013年1月24日木曜日 1時30分15秒 UTC+9 ONO Hiroki: