粗度長z0の設定について

[k-natsuaki]

お世話になっています。
現在、非線形 Shih k-εモデルを用いて計算を行っています。
境界条件で粗度長を設定して計算を行いたいと思っています。
nut内をnutkAtmRoughWallFunctionを利用しました。

        type            nutkAtmRoughWallFunction;
        z0               0.1;
        value           uniform 0;

しかし、nut内の粗度長z0を変化させて計算したところ、計算結果に変化が見られませんでした。
非線形 Shih k-εモデルで粗度長を設定することはできないのでしょうか。

お忙しいかと思いますが、よろしくお願いいたします。

[Masashi Imano]

今野です。

> 非線形 Shih k-εモデルで粗度長を設定することはできないのでしょうか。

以下のように、nonLinearKEShihモデルではnutの壁関数は非線型モデル
用がインクルードされているので、他のkEpsiklon等と異なりnutのtypeを変更しても
結果は変わらないと思います。

nonLinearKEShihモデル
src/turbulenceModels/incompressible/RAS/NonlinearKEShih/NonlinearKEShih.C

nut_ = Cmu_*sqr(k_)/epsilon_;

#include "wallNonlinearViscosityI.H"

kEpsilonモデル
src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C

// Re-calculate viscosity
nut_ = Cmu_*sqr(k_)/epsilon_;
nut_.correctBoundaryConditions();

よって、nonLinearKEShihモデルでz0型壁関数を用いる場合には、
include以下にある非線型モデル用の壁関数ソースをカスタマイズする
必要があると思います。

下記スレッドで既にS.kitaさんがこの件について詳細にコメントしているので、
こちらも参考になさってください。

https://groups.google.com/forum/#!topic/openfoam/PDbfH2DLqhU/discussion

以上です。

2012年11月9日 12:15 k-natsuaki <k-nat...@wind.ce.tokushima-u.ac.jp>:

> --
> このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> このディスカッションをウェブ上で閲覧するには、https://groups.google.com/d/msg/openfoam/-/mRu6Ln4EOngJ
> にアクセスしてください。
> このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
> このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
> 詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。



--
IMANO Masashi, Ph.D.