境界条件の与え方の質問

[tmtm19...@gmail.com]

OpenFOAMを大学の研究で使っている草野と申します．

現在，超音速流れの研究をしています．

超音速ノズルと噴流を同時に解析すると解が発散してしまうので，

ノズルと噴流別々に解析を行なっています．

その手順として，ノズルをまず初めに解析し，ノズル出口中心のの物理量を

噴流入口の境界値に一様に与えて行なっています．

しかし，ノズル出口の物理量は断面で分布をもっているので，その分布を

そのまま噴流入口の境界値に与えることはできるのでしょうか．

どなたかわかる方がいらしゃいましたらご教授お願いいたします．

[大渕真志]

別ケースとして計算しているのですね。ノズル計算結果の出口パッチ値を

（座標、値）のテーブルファイルとして保存し、timeVariyingMappedFixedValue

を使って噴流計算ケースの流入パッチに取り込めば良いと思います。

timeVaryingMappedFixedValueの使い方は、過去ログを検索して下さい。

ノズル出口パッチの計算値取得はsampleユーティリティ等でできます。

ご参考まで。

2013年7月4日 18:33 <tmtm19...@gmail.com>:

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

[tmtm19...@gmail.com]

ohbuchi様 ご教授ありがとうございます．

教えていただいたことを試してみたいと思います．

詳細な設定方法で質問させていただくかも知れませんが，

その際は何卒よろしくお願い致します．

[tmtm19...@gmail.com]

ohbuchi様　ご教授ありがとうございました.

過去ログを調べ，tutorials/incompressible/simpleFoam/pitzDailyExptInlet
を参考にして．以下のように記述したのですがうまくいきません.

Points（座標）

(
(0.0000 0.0000  0.0000)
(0.0000 1.0000e-04      0.0000)
(0.0000 0.00020000      0.0000)
(0.0000 0.00030000      0.0000)
(0.0000 0.00040000      0.0000)
(0.0000 0.00050000      0.0000)
(0.0000 0.00060000      0.0000)
(0.0000 0.00070000      0.0000)
(0.0000 0.00080000      0.0000)
(0.0000 0.00090000      0.0000)
(0.0000 0.0010000       0.0000)
(0.0000 0.0011000       0.0000)
(0.0000 0.0012000       0.0000)
(0.0000 0.0013000       0.0000)
(0.0000 0.0014000       0.0000)
(0.0000 0.0015000       0.0000)
(0.0000 0.0016000       0.0000)
(0.0000 0.0017000       0.0000)
(0.0000 0.0018000       0.0000)
(0.0000 0.0019000       0.0000)
(0.0000 0.0020000       0.0000)
(0.0000 0.0021000       0.0000)
(0.0000 0.0022000       0.0000)
(0.0000 0.0023000       0.0000)
(0.0000 0.0024000       0.0000)
(0.0000 0.0025000       0.0000)
(0.0000 0.0026000       0.0000)
(0.0000 0.0027000       0.0000)
(0.0000 0.0028000       0.0000)
(0.0000 0.0029000       0.0000)
(0.0000 0.0030000       0.0000)
(0.0000 0.0031000       0.0000)
(0.0000 0.0032000       0.0000)
(0.0000 0.0033000       0.0000)
(0.0000 0.0034000       0.0000)
(0.0000 0.0035000       0.0000)
(0.0000 0.0036000       0.0000)
(0.0000 0.0037000       0.0000)
(0.0000 0.0038000       0.0000)
(0.0000 0.0039000       0.0000)
(0.0000 0.0040000       0.0000)
(0.0000 0.0041000       0.0000)
(0.0000 0.0042000       0.0000)
(0.0000 0.0043000       0.0000)
(0.0000 0.0044000       0.0000)
(0.0000 0.0045000       0.0000)
(0.0000 0.0046000       0.0000)
(0.0000 0.0047000       0.0000)
(0.0000 0.0048000       0.0000)
(0.0000 0.0049000       0.0000)
(0.0000 0.0050000       0.0000)
(0.0000 0.0051000       0.0000)
(0.0000 0.0052000       0.0000)
(0.0000 0.0053000       0.0000)
(0.0000 0.0054000       0.0000)
)
// ************************************************************************* //

U（速度）

(0 0 0)

55
(
(0.0000 0.0000  1194.4)
(0.0000 0.0000  1198.3)
(0.0000 0.0000  1197.7)
(0.0000 0.0000  1197.0)
(0.0000 0.0000  1196.4)
(0.0000 0.0000  1195.7)
(0.0000 0.0000  1195.1)
(0.0000 0.0000  1194.5)
(0.0000 0.0000  1192.2)
(0.0000 0.0000  1189.7)
(0.0000 0.0000  1187.3)
(0.0000 0.0000  1184.9)
(0.0000 0.0000  1182.5)
(0.0000 0.0000  1180.1)
(0.0000 0.0000  1176.5)
(0.0000 0.0000  1172.5)
(0.0000 0.0000  1168.6)
(0.0000 0.0000  1164.7)
(0.0000 0.0000  1160.7)
(0.0000 0.0000  1156.3)
(0.0000 0.0000  1151.1)
(0.0000 0.0000  1145.8)
(0.0000 0.0000  1140.6)
(0.0000 0.0000  1135.4)
(0.0000 0.0000  1129.3)
(0.0000 0.0000  1122.9)
(0.0000 0.0000  1116.6)
(0.0000 0.0000  1110.2)
(0.0000 0.0000  1103.1)
(0.0000 0.0000  1095.6)
(0.0000 0.0000  1088.1)
(0.0000 0.0000  1080.6)
(0.0000 0.0000  1072.1)
(0.0000 0.0000  1063.4)
(0.0000 0.0000  1054.8)
(0.0000 0.0000  1045.3)
(0.0000 0.0000  1035.5)
(0.0000 0.0000  1025.7)
(0.0000 0.0000  1014.7)
(0.0000 0.0000  1003.7)
(0.0000 0.0000  991.83)
(0.0000 0.0000  979.51)
(0.0000 0.0000  966.30)
(0.0000 0.0000  952.47)
(0.0000 0.0000  937.35)
(0.0000 0.0000  921.51)
(0.0000 0.0000  904.02)
(0.0000 0.0000  884.82)
(0.0000 0.0000  863.66)
(0.0000 0.0000  839.82)
(0.0000 0.0000  812.39)
(0.0000 0.0000  780.07)
(0.0000 0.0000  738.83)
(0.0000 0.0000  683.24)
(0.0000 0.0000  590.70)
)
// ************************************************************************* //

他の物理量はスカラーで速度と同じように記述しました．
０ファイルのinletの境界条件はすべて以下のように記述しました.

inlet
    {
       　type            timeVaryingMappedFixedValue;
        　setAverage      off;
    }

私の考えでsetAverageが怪しいと思っていますが何分初めて使う境界条件なので,
困っております．

お忙しいとは存じますが何卒よろしくお願い致します．





 

[ohbuchi]

pitzDailyExptInletチュートリアルは二次元計算の例で、メッシュはｚ方向に１層になっており
ポイントデータはセル厚みの両端(minz,maxz)で定義されている筈です。
草野さんが提示されたポイントデータは１本のライン上でのみ定義されている様に見えます。

ポイントから補間するには、パッチ中心座標がポイントで囲まれないとうまく補間できないと
思われます。

ご参考まで。


2013年7月9日火曜日 15時38分31秒 UTC+9 tmtm19...@gmail.com:

[tmtm19...@gmail.com]

ohbuchi様　ご教授ありがとうございます．

教えて頂いたとおりにしたところ、うまくいきました．

今後わからないことがあるかもしれないので、

その際は何卒よろしくお願い致します．