チュートリアル(dnsFoam)

[hanihani]

いま，dnsFoamのチュートリアルをもとに流れ場を変えてるのですがなかなかうまくいきません．

 

チュートリアルのdnsFoam/boxTurb16をもとにblockMeshDictを編集し計算領域を平行平板間チャネル流れに拡張しました(メッシュも細かくしました)．又、初期状態を静止流体として境界条件も上壁と下壁をcyclic→wallへと変更しました．

 

メッシュの作成，初期境界条件のコンパイルは何事も無かったのですが，turbulencePropertiesでERRORを返されてしまい，原因をつかみきれていません．

 

端末に表示されたものです・・・

\*---------------------------------------------------------------------------*/
Build  : 2.2.x-45beea08cbf7
Exec   : dnsFoam
Date   : Oct 27 2013
Time   : 00:28:21
Host   : "*****"
PID    : 3888
Case   : /home/eustia/Desktop/研究/channel/dnsFoam/channel
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, no clear-out for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Reading turbulenceProperties

 

--> FOAM FATAL ERROR:
calculated number of cells is incorrect

 

    From function Kmesh::Kmesh(const fvMesh& mesh)
    in file Kmesh/Kmesh.C at line 76.

 

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/home/eustia/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/home/eustia/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::Kmesh::Kmesh(Foam::fvMesh const&) in "/home/eustia/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/librandomProcesses.so"
#3 
 in "/home/eustia/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/dnsFoam"
#4  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#5 
 in "/home/eustia/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/dnsFoam"

中止(コアダンプ)

 

ERROR内容を見ると計算されたセルの個数が違うとありますが，turbulencePropertiesにはセル数のかかわるようなものはないのではないでしょうか？

その後のFOAM abortingについてもいまいちよくわかっていません・・・

すみませんが原因のほうを教えていただけませんでしょうか？

[大渕真志]

こんにちは。

チャネル流のDNSを行うならpisoFoamかpimpleFoamを用いることをお勧めします。

dnsFoamは等方性乱流計算用の特殊処理が組み込まれています。

エラーはdnsFoam内でFFTの計算するときにセル数が２のべき乗でないために生じている

と思います。

2013年10月25日 21:49 hanihani <kanade...@kke.biglobe.ne.jp>:

--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
このグループから退会し、メールの受信を停止するには、openfoam+u...@googlegroups.com にメールを送信します。
このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
http://groups.google.com/group/openfoam からこのグループにアクセスしてください。
その他のオプションについては、https://groups.google.com/groups/opt_out にアクセスしてください。

[hanihani]

ありがとうございます．

pimpleoFoamを用いるということですと，pimpleFoam/channel395がLESでの記述となっているのでこれの諸々をDNSに変えてやればいいわけですね？

がんばります・・・

 

 

[Tatsuya SHIMIZU]

横から失礼します。
dnsはスペクトル法以外に有りませんので、チャンネルをスペクトル法で解くすべは有りません。pisoFoamやpimpleFoam（les）で解いてもマルチスケールの全エネルギーはほぞんされませんのでdnsとは呼ばないと思います。当然、openfoamのdnsfoamはスペクトル法です

2013/10/25 23:29 "大渕真志" <ohbu...@gmail.com>:

[hanihani]

LongGe 様

 

そうなんですか・・・

では自分でSMAC法や圧力方程式を解くプログラムを組まないといけないんですね・・・

[Tatsuya SHIMIZU]

お疲れ様です
smac法はスペクトル法とは全く関係無いと思います。それともスペクトル法の時間進行法にsmac法を使った文献あるのてすか？だとしたら、私の不勉強ですね。スペクトル法に時間進行があったかな？周波数領域の話ですから

2013/10/26 0:22 "hanihani" <kanade...@kke.biglobe.ne.jp>:

LongGe 様

 

そうなんですか・・・

では自分でSMAC法や圧力方程式を解くプログラムを組まないといけないんですね・・・

[ohbuchi]

おはようございます。dnsFoamは単なるPISO法のソルバです。

boxTurb16チュートリアルでは対流項離散化にGausss cubic（4次中心差分）

を使っています。

スペクトル法は確かに最も精度の高い計算手法ですが、その特徴である

スペクトル精度（分解能を上げた時、誤差が指数的に減少すること）を

発揮するには変数が分解能に応じた回数だけ微分可能であることを要求

します。実際の流れではこれは満たせないので、現実的には高次精度の

中心差分スキームでも十分な場合が多いと思います。4次で十分かは別

問題ですが。

ということで、チャネル流の計算では、対流項の離散化をなるべく高次

精度のものを使い、メッシュ分解能に注意して計算すれば良いと思います。

2013年10月26日土曜日 0時25分55秒 UTC+9 LongGe:

お疲れ様です
smac法はスペクトル法とは全く関係無いと思います。それともスペクトル法の時間進行法にsmac法を使った文献あるのてすか？だとしたら、私の不勉強ですね。スペクトル法に時間進行があったかな？周波数領域の話ですから

2013/10/26 0:22 "hanihani" <kanade...@kke.biglobe.ne.jp>:

LongGe 様

 

そうなんですか・・・

では自分でSMAC法や圧力方程式を解くプログラムを組まないといけないんですね・・・

--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。

このグループから退会し、メールの受信を停止するには、openfoam+unsubscribe@googlegroups.com にメールを送信します。

[Tatsuya SHIMIZU]

すいません、私の確認不足でした。dnsFoamは乱流モデルを使わないというだけのソルバーでした。たぶん、これだとレイノルズ数が大きく、せん断乱流になると破綻します

2013/10/26 5:47 "ohbuchi" <ohbu...@gmail.com>:

このグループから退会し、メールの受信を停止するには、openfoam+u...@googlegroups.com にメールを送信します。

[hanihani]

Ohbuchi 様

LongGe 様

様々なご指摘ありがとうございます．

まだまだ勉強することがいっぱいありますね・・・

また質問させていただくと思いますのでその時はよろしくお願いします．

[Tatsuya SHIMIZU]

こちらこそ不勉強で申し訳ありません。ただ、DNSをやる理由で方法論が変わると思います。乱流現象を追うのであればスペクトル法以外に方法は無いと現状では思います。ご指摘があった問題をスペクトル法は持っているものの、それ以外の方法はもっとレベルの低いところでつまずきます。よって、目的をはっきりと持って方法を選択してください。

2013/10/26 19:52 "hanihani" <kanade...@kke.biglobe.ne.jp>:

--