化学反応と動粘度係数について

[kk]

数値解析並びにOpenFoam初心者のkkと申します。

周囲にOpenFoamを使っている人がいないため苦戦しております。
どなたかアドバイスをお願いいたします。

・化学反応に伴う流体の粘度変化によって流体がどのように動くのかを研究しているのですが、OpenFoamの中のどのソルバーを使うのがいいでしょう
か？？。層流で非圧縮性流体を仮定し、流体の流速や圧力などを計算しようと考えています。

・OpenFoamでは、時間や場所によって変化する動粘度係数の指定はできるのでしょうか？

以上２点について詳しい方いらっしゃいましたらお願いいたします。

[ohb...@amber.plala.or.jp]

こんばんは。とても難しい問題にチャレンジなさっているのですね。

> ・化学反応に伴う流体の粘度変化によって流体がどのように動くのかを研究しているのですが、OpenFoamの中のどのソルバーを使うのがいいでしょう
> か？？。層流で非圧縮性流体を仮定し、流体の流速や圧力などを計算しようと考えています。

私の知る範囲で、正規版、dev版ともに、kkさんの考えておられる問題をそのまま扱えるソルバはありません。
reactingFoamを元にして、カスタムソルバを作るのが良いと思います。

カスタムソルバを作ること自体はそれほど難しくないと思います。
applications/solverのreactingFoamのソースを自分のディレクトリにコピーし、名前を変えて編集します。
変更点は、粘性係数のフィールド定義と、化学種濃度から粘性係数を求める部分くらいだと思います。

＃結果が出たら是非公開してください。宜しくお願いします。

> ・OpenFoamでは、時間や場所によって変化する動粘度係数の指定はできるのでしょうか？

通常、動粘性係数は、constant/transportPropertiesに定数として定義します。普通、時空間的に変動する
ケースはないと思います。
しかし、乱流の渦粘性や、非ニュートン流体も扱えてるわけですし、フィールドデータとして定義すれば変動を
考慮すること自体は問題ありません。

[Youhei Takagi]

温度依存する場合はSutherlandの式が使えます。

2010年10月1日16:11 kk <kiyono...@hotmail.co.jp>:

--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。

[kk]

ohb...@amber.plala.or.jpさん

丁寧な返信大変ありがとうございます。早速reactingFoamでカスタムソルバを作ってみようかと思います。

＞＃結果が出たら是非公開してください。
シミュレーションは初心者なのでなかなか時間がかかるとは思いますが結果が出たら是非公開させていただきます。
また何かありましたらこの掲示板で質問することがあるかと存じますがその際はよろしくお願いいたします。
ありがとうございました。

[kk]

Youhei Takagiさん

返信ありがとうございます。
自分の方でSutherlandの式について調べてみたいと思います。